Large-scale analysis of orthologs and paralogs under covarion-like and constant-but-different models of amino acid evolution.

Functional divergence between homologous proteins is expected to affect amino acid sequences in two main ways, which can be considered as proxies of biochemical divergence: a "covarion-like" pattern of correlated changes in evolutionary rates, and switches in conserved residues ("conserved but different"). Although these patterns have been used in case studies, a large-scale analysis is needed to estimate their frequency and distribution. We use a phylogenomic framework of animal genes to answer three questions: 1) What is the prevalence of such patterns? 2) Can we link such patterns at the amino acid level with selection inferred at the codon level? 3) Are patterns different between paralogs and orthologs? We find that covarion-like patterns are more frequently detected than "constant but different," but that only the latter are correlated with signal for positive selection. Finally, there is no obvious difference in patterns between orthologs and paralogs.

Comparative analysis of human and mouse expression data illuminates tissue-specific evolutionary patterns of miRNAs.

MicroRNAs (miRNAs) constitute an important class of gene regulators. While models have been proposed to explain their appearance and expansion, the validation of these models has been difficult due to the lack of comparative studies. Here, we analyze miRNA evolutionary patterns in two mammals, human and mouse, in relation to the age of miRNA families. In this comparative framework, we confirm some predictions of previously advanced models of miRNA evolution, e.g. that miRNAs arise more frequently de novo than by duplication, or that the number of protein-coding gene targeted by miRNAs decreases with evolutionary time. We also corroborate that miRNAs display an increase in expression level with evolutionary time, however we show that this relation is largely tissue-dependent, and especially low in embryonic or nervous tissues. We identify a bias of tag-sequencing techniques regarding the assessment of breadth of expression, leading us, contrary to predictions, to find more tissue-specific expression of older miRNAs. Together, our results refine the models used so far to depict the evolution of miRNA genes. They underline the role of tissue-specific selective forces on the evolution of miRNAs, as well as the potential co-evolution patterns between miRNAs and the protein-coding genes they target.

Efficacy of prescribed injectable diacetylmorphine in the Andalusian trial: Bayesian analysis of responders and non-responders according to a multi domain outcome index.

BACKGROUND The objective of this research was to evaluate data from a randomized clinical trial that tested injectable diacetylmorphine (DAM) and oral methadone (MMT) for substitution treatment, using a multi-domain dichotomous index, with a Bayesian approach. METHODS Sixty two long-term, socially-excluded heroin injectors, not benefiting from available treatments were randomized to receive either DAM or MMT for 9 months in Granada, Spain. Completers were 44 and data at the end of the study period was obtained for 50. Participants were determined to be responders or non responders using a multi-domain outcome index accounting for their physical and mental health and psychosocial integration, used in a previous trial. Data was analyzed with Bayesian methods, using information from a similar study conducted in The Netherlands to select a priori distributions. On adding the data from the present study to update the a priori information, the distribution of the difference in response rates were obtained and used to build credibility intervals and relevant probability computations. RESULTS In the experimental group (n = 27), the rate of responders to treatment was 70.4% (95% CI 53.287.6), and in the control group (n = 23), it was 34.8% (95% CI 15.354.3). The probability of success in the experimental group using the a posteriori distributions was higher after a proper sensitivity analysis. Almost the whole distribution of the rates difference (the one for diacetylmorphine minus methadone) was located to the right of the zero, indicating the superiority of the experimental treatment. CONCLUSION The present analysis suggests a clinical superiority of injectable diacetylmorphine compared to oral methadone in the treatment of severely affected heroin injectors not benefiting sufficiently from the available treatments. TRIAL REGISTRATION Current Controlled Trials ISRCTN52023186.

Analysis of Functional and Anatomic Success Following Endonasal Dacryocystorhinostomy

Objectives: The aim of this study was to assess the concomitant perioperative procedures, the causes of nasolacrimal duct obstruction, the success rate, and the complications associated with endonasal dacryocystorhinostomy (ENDCR).Methods: In this single-center retrospective study, 98 patients underwent 104 ENDCRs between January 1994 and February 2006. There were 78 patients with 84 nasolacrimal duct obstructions who were included in this study.Results: The overall functional success rate with improvement in symptoms was 94.9% for primary surgery (59 of 84 obstructions) and 63.6% for salvage surgery after failure of primary surgery performed in another hospital (25 of 84 obstructions). The mean follow-up time was 36.8 +/- 17.11 months. Primary surgery showed better results, with a complete success rate of 93.2%, than did salvage surgery, with a success rate of only 68%. Persistent symptoms, despite an open rhinostomy, were found in 1.7% of patients with primary surgery and in 12% of those with salvage surgery. Failure of ENDCR was observed in 3.4% of patients after primary surgery and in 20% after salvage surgery. We encountered only minimal perioperative complications, and these were essentially related to lacrimal intubation.Conclusions: Because of the possibility of treating concomitant sinonasal disorders, the cosmetic advantages, and the excellent results, ENDCR represents the procedure of choice for treating nasolacrimal duct obstructions. The main challenge lies in the exact preoperative assessment, as well as postoperative evaluation in case of failure.

L'analyse des renvois régionaux et interrégionaux en pratique médicale ambulatoire. [Analysis of regional and interregional referrals in ambulatory medical practice].

This article examines, in two Swiss cantons, the interdependence from a medical care point of view of various regions (health planning zones in one canton, political districts in the other). The volume and the destination of patient referrals prescribed by physicians in ambulatory practice are analyzed. The available data (on 1609 referrals) were gathered by the practitioners themselves, during a National Ambulatory Medical Care Survey type study in February-March 1981, in which 203 physicians participated. Several indicators are proposed (including an integration coefficient and an attraction coefficient for each zone); they show marked differences among the regions. This dynamic approach, based on the effective behavior of physicians, appears to be of major interest for health planning purposes (as compared with the frequent practice to use mainly parameters in relation with the availability of care services--the "supply"--numbers of professionals and/or health facilities).

Pseudogenes in the ENCODE regions: consensus annotation, analysis of transcription, and evolution

Arising from either retrotransposition or genomic duplication of functional genes, pseudogenes are “genomic fossils” valuable for exploring the dynamics and evolution of genes and genomes. Pseudogene identification is an important problem in computational genomics, and is also critical for obtaining an accurate picture of a genome’s structure and function. However, no consensus computational scheme for defining and detecting pseudogenes has been developed thus far. As part of the ENCyclopedia Of DNA Elements (ENCODE) project, we have compared several distinct pseudogene annotation strategies and found that different approaches and parameters often resulted in rather distinct sets of pseudogenes. We subsequently developed a consensus approach for annotating pseudogenes (derived from protein coding genes) in the ENCODE regions, resulting in 201 pseudogenes, two-thirds of which originated from retrotransposition. A survey of orthologs for these pseudogenes in 28 vertebrate genomes showed that a significant fraction (∼80%) of the processed pseudogenes are primate-specific sequences, highlighting the increasing retrotransposition activity in primates. Analysis of sequence conservation and variation also demonstrated that most pseudogenes evolve neutrally, and processed pseudogenes appear to have lost their coding potential immediately or soon after their emergence. In order to explore the functional implication of pseudogene prevalence, we have extensively examined the transcriptional activity of the ENCODE pseudogenes. We performed systematic series of pseudogene-specific RACE analyses. These, together with complementary evidence derived from tiling microarrays and high throughput sequencing, demonstrated that at least a fifth of the 201 pseudogenes are transcribed in one or more cell lines or tissues.

Analysis of epidermis- and mesophyll-specific transcript accumulation in powdery mildew-inoculated wheat leaves.

Powdery mildew is an important disease of wheat caused by the obligate biotrophic fungus Blumeria graminis f. sp. tritici. This pathogen invades exclusively epidermal cells after penetrating directly through the cell wall. Because powdery mildew colonizes exclusively epidermal cells, it is of importance not only to identify genes which are activated, but also to monitor tissue specificity of gene activation. Acquired resistance of wheat to powdery mildew can be induced by a previous inoculation with the non-host pathogen B. graminis f. sp. hordei, the causal agent of barley powdery mildew. The establishment of the resistant state is accompanied by the activation of genes. Here we report the tissue-specific cDNA-AFLP analysis and cloning of transcripts accumulating 6 and 24 h after the resistance-inducing inoculation with B. graminis f. sp. hordei. A total of 25,000 fragments estimated to represent about 17,000 transcripts were displayed. Out of these, 141 transcripts, were found to accumulate after Bgh inoculation using microarray hybridization analysis. Forty-four accumulated predominantly in the epidermis whereas 76 transcripts accumulated mostly in mesophyll tissue.

Measures of fit in multiple correspondence analysis of crisp and fuzzy coded data

When continuous data are coded to categorical variables, two types of coding are possible: crisp coding in the form of indicator, or dummy, variables with values either 0 or 1; or fuzzy coding where each observation is transformed to a set of "degrees of membership" between 0 and 1, using co-called membership functions. It is well known that the correspondence analysis of crisp coded data, namely multiple correspondence analysis, yields principal inertias (eigenvalues) that considerably underestimate the quality of the solution in a low-dimensional space. Since the crisp data only code the categories to which each individual case belongs, an alternative measure of fit is simply to count how well these categories are predicted by the solution. Another approach is to consider multiple correspondence analysis equivalently as the analysis of the Burt matrix (i.e., the matrix of all two-way cross-tabulations of the categorical variables), and then perform a joint correspondence analysis to fit just the off-diagonal tables of the Burt matrix - the measure of fit is then computed as the quality of explaining these tables only. The correspondence analysis of fuzzy coded data, called "fuzzy multiple correspondence analysis", suffers from the same problem, albeit attenuated. Again, one can count how many correct predictions are made of the categories which have highest degree of membership. But here one can also defuzzify the results of the analysis to obtain estimated values of the original data, and then calculate a measure of fit in the familiar percentage form, thanks to the resultant orthogonal decomposition of variance. Furthermore, if one thinks of fuzzy multiple correspondence analysis as explaining the two-way associations between variables, a fuzzy Burt matrix can be computed and the same strategy as in the crisp case can be applied to analyse the off-diagonal part of this matrix. In this paper these alternative measures of fit are defined and applied to a data set of continuous meteorological variables, which are coded crisply and fuzzily into three categories. Measuring the fit is further discussed when the data set consists of a mixture of discrete and continuous variables.

The investigation of deaths in custody: a qualitative analysis of problems and prospects.

The right to be treated humanely when detained is universally recognized. Deficiencies in detention conditions and violence, however, subvert this right. When this occurs, proper medico-legal investigations are critical irrespective of the nature of death. Unfortunately, the very context of custody raises serious concerns over the effectiveness and fairness of medico-legal examinations. The aim of this manuscript is to identify and discuss the practical and ethical difficulties encountered in the medico-legal investigation following deaths in custody. Data for this manuscript come from a larger project on Death in Custody that examined the causes of deaths in custody and the conditions under which these deaths should be investigated and prevented. A total of 33 stakeholders from forensic medicine, law, prison administration or national human rights administration were interviewed. Data obtained were analyzed qualitatively. Forensic experts are an essential part of the criminal justice process as they offer evidence for subsequent indictment and eventual punishment of perpetrators. Their independence when investigating a death in custody was deemed critical and lack thereof, problematic. When experts were not independent, concerns arose in relation to conflicts of interest, biased perspectives, and low-quality forensic reports. The solutions to ensure independent forensic investigations of deaths in custody must be structural and simple: setting binding standards of practice rather than detailed procedures and relying on preexisting national practices as opposed to encouraging new practices that are unattainable for countries with limited resources.

Integrated Analysis of Molecular and Clinical Prognostic Factors in Stage II/III Colon Cancer.

Background The prognostic potential of individual clinical and molecular parameters in stage II/III colon cancer has been investigated, but a thorough multivariable assessment of their relative impact is missing. Methods Tumors from patients (N = 1404) in the PETACC3 adjuvant chemotherapy trial were examined for BRAF and KRAS mutations, microsatellite instability (MSI), chromosome 18q loss of heterozygosity (18qLOH), and SMAD4 expression. Their importance in predicting relapse-free survival (RFS) and overall survival (OS) was assessed by Kaplan-Meier analyses, Cox regression models, and recursive partitioning trees. All statistical tests were two-sided. Results MSI-high status and SMAD4 focal loss of expression were identified as independent prognostic factors with better RFS (hazard ratio [HR] of recurrence = 0.54, 95% CI = 0.37 to 0.81, P = .003) and OS (HR of death = 0.43, 95% CI = 0.27 to 0.70, P = .001) for MSI-high status and worse RFS (HR = 1.47, 95% CI = 1.19 to 1.81, P < .001) and OS (HR = 1.58, 95% CI = 1.23 to 2.01, P < .001) for SMAD4 loss. 18qLOH did not have any prognostic value in RFS or OS. Recursive partitioning identified refinements of TNM into new clinically interesting prognostic subgroups. Notably, T3N1 tumors with MSI-high status and retained SMAD4 expression had outcomes similar to stage II disease. Conclusions Concomitant assessment of molecular and clinical markers in multivariable analysis is essential to confirm or refute their independent prognostic value. Including molecular markers with independent prognostic value might allow more accurate prediction of prognosis than TNM staging alone.

Analysis of CD4+ and CD8+ T cells by defined MHC-peptide complexes

MHC class II-peptide multimers are important tools for the detection, enumeration and isolation of antigen-specific CD4+ Τ cells. However, their erratic and often poor performance impeded their broad application and thus in-depth analysis of key aspects of antigen-specific CD4+ Τ cell responses. In the first part of this thesis we demonstrate that a major cause for poor MHC class II tetramer staining performance is incomplete peptide loading on MHC molecules. We observed that peptide binding affinity for "empty" MHC class II molecules poorly correlates with peptide loading efficacy. Addition of a His-tag or desthiobiotin (DTB) at the peptide N-terminus allowed us to isolate "immunopure" MHC class II-peptide monomers by affinity chromatography; this significantly, often dramatically, improved tetramer staining of antigen-specific CD4+ Τ cells. Insertion of a photosensitive amino acid between the tag and the peptide, permitted removal of the tag from "immunopure" MHC class II-peptide complex by UV irradiation, and hence elimination of its potential interference with TCR and/or MHC binding. Moreover, to improve loading of self and tumor antigen- derived peptides onto "empty" MHC II molecules, we first loaded these with a photocleavable variant of the influenza A hemagglutinin peptide HA306-318 and subsequently exchanged it with a poorly loading peptide (e.g. NY-ESO-1119-143) upon photolysis of the conditional ligand. Finally, we established a novel type of MHC class II multimers built on reversible chelate formation between 2xHis-tagged MHC molecules and a fluorescent nitrilotriacetic acid (NTA)-containing scaffold. Staining of antigen-specific CD4+ Τ cells with "NTAmers" is fully reversible and allows gentle cell sorting. In the second part of the thesis we investigated the role of the CD8α transmembrane domain (TMD) for CD8 coreceptor function. The sequence of the CD8α TMD, but not the CD8β TMD, is highly conserved and homodimerizes efficiently. We replaced the CD8α TMD with the one of the interleukin-2 receptor a chain (CD8αTac) and thus ablated CD8α TMD interactions. We observed that ΤΙ Τ cell hybridomas expressing CD8αTacβ exhibited severely impaired intracellular calcium flux, IL-2 responses and Kd/PbCS(ABA) P255A tetramer binding. By means of fluorescence resonance energy transfer experiments (FRET) we established that CD8αTacβ associated with TCR:CD3 considerably less efficiently than CD8αβ, both in the presence and the absence of Kd/PbCS(ABA) complexes. Moreover, we observed that CD8αTacβ partitioned substantially less in lipid rafts, and related to this, associated less efficiently with p56Lck (Lck), a Src kinase that plays key roles in TCR proximal signaling. Our results support the view that the CD8α TMD promotes the formation of CD8αβP-CD8αβ dimers on cell surfaces. Because these contain two CD8β chains and that CD8β, unlike CD8α, mediates association of CD8 with TCR:CD3 as well as with lipid rafts and hence with Lck, we propose that the CD8αTMD plays an important and hitherto unrecognized role for CD8 coreceptor function, namely by promoting CD8αβ dimer formation. We discuss what implications this might have on TCR oligomerization and TCR signaling. - Les multimères de complexes MHC classe II-peptide sont des outils importants pour la détection, le dénombrement et l'isolation des cellules Τ CD4+ spécifiques pour un antigène d'intérêt. Cependant, leur performance erratique et souvent inadéquate a empêché leur utilisation généralisée, limitant ainsi l'analyse des aspects clés des réponses des lymphocytes Τ CD4+. Dans la première partie de cette thèse, nous montrons que la cause principale de la faible efficacité des multimères de complexes MHC classe II-peptide est le chargement incomplet des molécules MHC par des peptides. Nous montrons également que l'affinité du peptide pour la molécule MHC classe II "vide" n'est pas nécessairement liée au degré du chargement. Grâce à l'introduction d'une étiquette d'histidines (His-tag) ou d'une molécule de desthiobiotine à l'extrémité N-terminale du peptide, des monomères MHC classe II- peptide dits "immunopures" ont pu être isolés par chromatographic d'affinité. Ceci a permis d'améliorer significativement et souvent de façon spectaculaire, le marquage des cellules Τ CD4+ spécifiques pour un antigène d'intérêt. L'insertion d'un acide aminé photosensible entre l'étiquette et le peptide a permis la suppression de l'étiquette du complexe MHC classe- Il peptide "immunopure" par irradiation aux UV, éliminant ainsi de potentielles interférences de liaison au TCR et/ou au MHC. De plus, afin d'améliorer le chargement des molécules MHC classe II "vides" avec des peptides dérivés d'auto-antigènes ou d'antigènes tumoraux, nous avons tout d'abord chargé les molécules MHC "vides" avec un analogue peptidique photoclivable issu du peptide HA306-318 de l'hémagglutinine de la grippe de type A, puis, sous condition de photolyse, nous l'avons échangé avec de peptides à chargement faible (p.ex. NY-ESO-1119-143). Finalement, nous avons construit un nouveau type de multimère réversible, appelé "NTAmère", basé sur la formation chélatante reversible entre les molécules MHC-peptide étiquettés par 2xHis et un support fluorescent contenant des acides nitrilotriacetiques (NTA). Le marquage des cellules Τ CD4+ spécifiques pour un antigène d'intérêt avec les "NTAmères" est pleinement réversible et permet également un tri cellulaire plus doux. Dans la deuxième partie de cette thèse nous avons étudié le rôle du domaine transmembranaire (TMD) du CD8α pour la fonction coréceptrice du CD8. La séquence du TMD du CD8α, mais pas celle du TMD du CD8β, est hautement conservée et permet une homodimérisation efficace. Nous avons remplacé le TMD du CD8α avec celui de la chaîne α du récepteur à l'IL-2 (CD8αTac), éliminant ainsi les interactions du TMD du CD8α. Nous avons montré que les cellules des hybridomes Τ T1 exprimant le CD8αTacβ présentaient une atteinte sévère du flux du calcium intracellulaire, des réponses d'IL-2 et de la liaison des tétramères Kd/PbCS(ABA) P255A. Grâce aux expériences de transfert d'énergie entre molécules fluorescentes (FRET), nous avons montré que l'association du CD8αTacβ avec le TCR:CD3 est considérablement moins efficace qu'avec le CD8αβ, et ceci aussi bien en présence qu'en absence de complexes Kd/PbCS(ABA). De plus, nous avons observé que le CD8αTacβ se distribuait beaucoup moins bien dans les radeaux lipidiques, engendrant ainsi, une association moins efficace avec p56Lck (Lck), une kinase de la famille Src qui joue un rôle clé dans la signalisation proximale du TCR. Nos résultats soutiennent l'hypothèse que le TMD du CD8αβ favorise la formation des dimères de CD8αβ à la surface des cellules. Parce que ces derniers contiennent deux chaînes CD8β et que CD8β, contrairement à CD8α, favorise l'association du CD8 au TCR:CD3 aussi bien qu'aux radeaux lipidiques et par conséquent à Lck, nous proposons que le TMD du CD8α joue un rôle important, jusqu'alors inconnu, pour la fonction coreceptrice du CD8, en encourageant la formation des dimères CD8αβ. Nous discutons des implications possibles sur l'oligomerisation du TCR et la signalisation du TCR.

Asymptotic robust inferences in the analysis of mean and covariance structures

Structural equation models are widely used in economic, socialand behavioral studies to analyze linear interrelationships amongvariables, some of which may be unobservable or subject to measurementerror. Alternative estimation methods that exploit different distributionalassumptions are now available. The present paper deals with issues ofasymptotic statistical inferences, such as the evaluation of standarderrors of estimates and chi--square goodness--of--fit statistics,in the general context of mean and covariance structures. The emphasisis on drawing correct statistical inferences regardless of thedistribution of the data and the method of estimation employed. A(distribution--free) consistent estimate of $\Gamma$, the matrix ofasymptotic variances of the vector of sample second--order moments,will be used to compute robust standard errors and a robust chi--squaregoodness--of--fit squares. Simple modifications of the usual estimateof $\Gamma$ will also permit correct inferences in the case of multi--stage complex samples. We will also discuss the conditions under which,regardless of the distribution of the data, one can rely on the usual(non--robust) inferential statistics. Finally, a multivariate regressionmodel with errors--in--variables will be used to illustrate, by meansof simulated data, various theoretical aspects of the paper.

Can we afford to add chemotherapy to radiotherapy for glioblastoma multiforme? Cost-identification analysis of concomitant and adjuvant treatment with temozolomide until patient death.

BACKGROUND: Adding temozolomide (TMZ) to standard radiotherapy as a first-line therapy for glioma may increase costs to a disproportionate degree compared with the resulting survival benefits. METHODS: Forty-six consecutive patients (28 males and 18 females; median age, 52 years; age range, 24-70 years) received concomitant TMZ with radiotherapy for 6 weeks followed by adjuvant TMZ for 6 cycles, and they were followed until disease recurrence and then until death. The authors assessed the costs associated with the four phases of treatment from a hospital-centered perspective. RESULTS: Treatment was discontinued early in 3 patients, 9 patients, and 15 patients during concomitant TMZ, before adjuvant TMZ, and during adjuvant TMZ, respectively. Karnofsky index values varied between 85% (at the beginning of treatment) and 76% (at the end of treatment). The nature of care after disease recurrence was diverse. Overall survival ranged from 1.4 months to 64.3 months (median, 15.8 months) and was better if surgical debulking could be carried out before treatment. Global costs amounted to Euros 39,092 +/- Euros 21,948 (concomitant TMZ, Euros 14,539 +/- Euros 4998; adjuvant TMZ, Euros 13,651 +/- Euros 4320; follow-up, Euros 6363 +/- Euros 6917; and recurrence, Euros 12,344 +/- Euros 18,327), with 53% of these costs being related to the acquisition of TMZ; this represented an eightfold increase in cost compared with radiotherapy alone. CONCLUSIONS: TMZ may be an effective but costly adjuvant outpatient therapy for patients with glioblastoma multiforme. Definite cost-effectiveness/utility must be assessed in a randomized Phase III trial.

Analysis of past and future dam formation and failure in the Santa Cruz River (San Juan province, Argentina)

Around 11.5 * 106 m3 of rock detached from the eastern slope of the Santa Cruz valley (San Juan province, Argentina) in the first fortnight of January 2005. The rockslide?debris avalanche blocked the course, resulting in the development of a lake with maximum length of around 3.5 km. The increase in the inflow rate from 47,000?74,000 m3/d between April and October to 304,000 m3/d between late October and the first fortnight of November, accelerated the growing rate of the lake. On 12 November 2005 the dam failed, releasing 24.6 * 106 m3 of water. The resulting outburst flood caused damages mainly on infrastructure, and affected the facilities of a hydropower dam which was under construction 250 km downstream from the source area. In this work we describe causes and consequences of the natural dam formation and failure, and we dynamically model the 2005 rockslide?debris avalanche with DAN3D. Additionally, as a volume ~ 24 * 106 m3of rocks still remain unstable in the slope, we use the results of the back analysis to forecast the formation of a future natural dam. We analyzed two potential scenarios: a partial slope failure of 6.5 * 106 m3 and a worst case where all the unstable volume remaining in the slope fails. The spreading of those potential events shows that a new blockage of the Santa Cruz River is likely to occur. According to their modeled morphometry and the contributing watershed upstream the blockage area, as the one of 2005, the dams would also be unstable. This study shows the importance of back and forward analysis that can be carried out to obtain critical information for land use planning, hazards mitigation, and emergency management.