Rgtsp: a generalized top scoring pairs package for class prediction.

SUMMARY: A top scoring pair (TSP) classifier consists of a pair of variables whose relative ordering can be used for accurately predicting the class label of a sample. This classification rule has the advantage of being easily interpretable and more robust against technical variations in data, as those due to different microarray platforms. Here we describe a parallel implementation of this classifier which significantly reduces the training time, and a number of extensions, including a multi-class approach, which has the potential of improving the classification performance. AVAILABILITY AND IMPLEMENTATION: Full C++ source code and R package Rgtsp are freely available from http://lausanne.isb-sib.ch/~vpopovic/research/. The implementation relies on existing OpenMP libraries.

On the success of generalized food deception in orchids : the role of sympatric rewarding species and pollinators

Résumé : La production de nectar assure aux plantes entomophiles un important succès reproducteur. Malgré cela, de nombreuses espèces d'orchidées ne produisent pas de nectar. La majorité de ces orchidées dites trompeuses exploitent simplement l'instinct des pollinisateurs généralistes, qui les pousse à chercher du nectar dans les fleurs. Afin d'optimiser la récolte de nectar, les pollinisateurs apprennent à différencier les fleurs trompeuses des nectarifères, et à concentrer leurs visites sur ces dernières, au détriment des plantes trompeuses. Chez les orchidées non autogames, la reproduction est assurée uniquement par les pollinisateurs. L'apprentissage des pollinisateurs a donc un impact négatif sur la reproduction des orchidées trompeuses. Cependant, les caractéristiques d'une espèce trompeuse et des espèces nectarifères au sein d'une communauté végétale peuvent affecter l'apprentissage et le taux de visite des pollinisateurs aux plantes trompeuses. J'ai réalisé des expériences en milieu naturel et en milieu contrôlé, pour déterminer si les caractéristiques florales, spatiales et temporelles des communautés affectent le taux de visite et le succès reproducteur de plantes trompeuses. Une agrégation spatiale élevée des plantes trompeuses et des plantes nectarifères diminue le succès reproducteur des plantes trompeuses. De plus, les pollinisateurs visitent plus souvent l'espèce trompeuse Iorsque ses fleurs sont de couleur similaire à celles de l'espèce nectarifère. Cet effet bénéfique de la similarité pour la couleur des fleurs s'accentue si les deux espèces sont mélangées et proches spatialement, ou si l'espèce trompeuse fleurit après l'espèce nectarifère. Enfin, le comportement des pollinisateurs n'est pas tout de suite affecté lorsque les caractéristiques de la communauté changent. Les caractéristiques des communautés végétales affectent donc la reproduction des espèces trompeuses. Bien que L'absence de coûts associés à la production de nectar, l'exportation efficace de pollen et la production de graines de qualité dont bénéficient les orchidées trompeuses favorisent Ieur maintien, les caractéristiques de la communauté peuvent aussi y contribuer. Mon étude fournit donc une explication alternative et complémentaire au maintien des orchidées trompeuses. Je conclus par une discussion des implications possibles de ces résultats sur le maintien et l'évolution des orchidées trompeuses, en tenant compte de la dynamique des caractéristiques des communautés végétales naturelles. Abstract : Despite the importance of producing food to ensure a high reproductive success, many orchid species lack such rewards. The majority of deceptive orchids simply exploit the instinctive food-foraging behaviour of generalist pollinators. This strategy is termed generalized food deception. To optimize their foraging efficiency, pollinators can learn to discriminate deceptive from rewarding flowers and to focus their visits to the rewarding plants, to the disadvantage of the deceptive plants. Because the reproductive success of non-autogamous orchids entirely relies on pollinator visitation rate, pollinator learning decreases the reproductive success of deceptive orchids. However, the characteristics of deceptive and rewarding plants within a community may affect pollinator learning and visitation rate to a deceptive orchid. Therefore, the biological characteristics of natural plant communities may be crucial to the maintenance of generalized food deceptive orchids. My study focused on the floral, spatial and temporal characteristics of plant communities. I used both in and ex sitar experiments to investigate whether these characteristics influence pollinator visitation rates and the reproductive success of deceptive orchids. A high spatial aggregation of both deceptive and rewarding species decreased the reproductive success of the deceptive species. Also, being of similar flower colour to rewarding sympatric species increased pollinator visitation rates to a deceptive species. The beneficial effect of flower colour similarity was even more pronounced when both species were spatially closely mingled or when the deceptive species flowered after the rewarding species. Finally, pollinator behaviour was unaffected in the short term by a change in the characteristics of plant communities, indicating that pollinators need time to learn under new conditions. Thus, the characteristics of plant communities may crucially affect the reproductive success of deceptive orchids. Although the absence of costs associated with nectar production, the efficient pollen export and the high seed quality of deceptive orchids may favour their maintenance, the characteristics of plant communities may also contribute to it. Therefore, my study provides an alternative yet complementary explanation to the maintenance of generalized food deceptive orchids in natural populations. I discuss the possible implications for the maintenance and the evolution of generalized food deceptive orchids with regards to the floral and temporal dynamics of natural plant communities.

Derivation and validation of a prognostic model for pulmonary embolism.

RATIONALE: An objective and simple prognostic model for patients with pulmonary embolism could be helpful in guiding initial intensity of treatment. OBJECTIVES: To develop a clinical prediction rule that accurately classifies patients with pulmonary embolism into categories of increasing risk of mortality and other adverse medical outcomes. METHODS: We randomly allocated 15,531 inpatient discharges with pulmonary embolism from 186 Pennsylvania hospitals to derivation (67%) and internal validation (33%) samples. We derived our prediction rule using logistic regression with 30-day mortality as the primary outcome, and patient demographic and clinical data routinely available at presentation as potential predictor variables. We externally validated the rule in 221 inpatients with pulmonary embolism from Switzerland and France. MEASUREMENTS: We compared mortality and nonfatal adverse medical outcomes across the derivation and two validation samples. MAIN RESULTS: The prediction rule is based on 11 simple patient characteristics that were independently associated with mortality and stratifies patients with pulmonary embolism into five severity classes, with 30-day mortality rates of 0-1.6% in class I, 1.7-3.5% in class II, 3.2-7.1% in class III, 4.0-11.4% in class IV, and 10.0-24.5% in class V across the derivation and validation samples. Inpatient death and nonfatal complications were <or= 1.1% among patients in class I and <or= 1.9% among patients in class II. CONCLUSIONS: Our rule accurately classifies patients with pulmonary embolism into classes of increasing risk of mortality and other adverse medical outcomes. Further validation of the rule is important before its implementation as a decision aid to guide the initial management of patients with pulmonary embolism.

Fourier transform convolution integrals applied to generalized Born molecular volume.

Generalized Born methods are currently among the solvation models most commonly used for biological applications. We reformulate the generalized Born molecular volume method initially described by (Lee et al, 2003, J Phys Chem, 116, 10606; Lee et al, 2003, J Comp Chem, 24, 1348) using fast Fourier transform convolution integrals. Changes in the initial method are discussed and analyzed. Finally, the method is extensively checked with snapshots from common molecular modeling applications: binding free energy computations and docking. Biologically relevant test systems are chosen, including 855-36091 atoms. It is clearly demonstrated that, precision-wise, the proposed method performs as good as the original, and could better benefit from hardware accelerated boards.

Estimating parameters for generalized mass action models using constraint propagation.

As modern molecular biology moves towards the analysis of biological systems as opposed to their individual components, the need for appropriate mathematical and computational techniques for understanding the dynamics and structure of such systems is becoming more pressing. For example, the modeling of biochemical systems using ordinary differential equations (ODEs) based on high-throughput, time-dense profiles is becoming more common-place, which is necessitating the development of improved techniques to estimate model parameters from such data. Due to the high dimensionality of this estimation problem, straight-forward optimization strategies rarely produce correct parameter values, and hence current methods tend to utilize genetic/evolutionary algorithms to perform non-linear parameter fitting. Here, we describe a completely deterministic approach, which is based on interval analysis. This allows us to examine entire sets of parameters, and thus to exhaust the global search within a finite number of steps. In particular, we show how our method may be applied to a generic class of ODEs used for modeling biochemical systems called Generalized Mass Action Models (GMAs). In addition, we show that for GMAs our method is amenable to the technique in interval arithmetic called constraint propagation, which allows great improvement of its efficiency. To illustrate the applicability of our method we apply it to some networks of biochemical reactions appearing in the literature, showing in particular that, in addition to estimating system parameters in the absence of noise, our method may also be used to recover the topology of these networks.

Primary γδ T cell lymphoma of the lung: report of a case with features suggesting derivation from intraepithelial γδ T lymphocytes [Primary gammadelta T cell lymphoma of the lung: report of a case with features suggesting derivation from intraepithelial gammadelta T lymphocytes].

T cell lymphoma of γδ T cell origin is a rare disease that mainly involves extranodal sites and shows aggressive clinical behavior. Here, we report a case of primary γδ T cell lymphoma of the lungs with epitheliotropism in the respiratory epithelium, a feature somewhat reminiscent of what is observed in enteropathy-associated T cell lymphoma. A 63-year-old man presented with chest pain and dyspnea on exertion, weight loss, and general weakness. On a positron emission tomography (PET) scan, multiple hypermetabolic lesions were found in both lungs. Microscopic examination of the wedge lung biopsy revealed nodular infiltration of monomorphic, medium- to large-sized atypical lymphocytes with round nuclei, coarse chromatin, and a variable amount of clear to eosinophilic cytoplasm. Of note, intraepithelial lymphocytosis by atypical lymphoid cells was observed in the respiratory epithelium within and around the nodule. Immunohistochemically, the tumor cells were CD3+, TCRβF1-, TCRγ+, CD5-, CD7+, CD20-, CD79a-, CD30-, CD4-, CD8-, CD10-, BCL6-, CD21-, CD56+, CD57-, and CD138-, and expressed cytotoxic molecules. Epstein-Barr virus (EBV) was not detected by an in situ hybridization assay for EBV-encoded RNA. Interestingly, CD103 was expressed by a subset of tumor cells, especially those infiltrating the epithelium. T cell clonality was detected by multiplex PCR analysis of TRG and TRD gene rearrangements. After 2 months of systemic chemotherapy, PET scan showed regression of the size and metabolic activity of the lesions. This case represents a unique γδ T cell lymphoma of the lungs showing epitheliotropism by CD103+ γδ T cells that is suggestive of tissue-resident γδ T cells as the cell of origin.

Derivation of traceable and transplantable photoreceptors from mouse embryonic stem cells.

Retinal degenerative diseases resulting in the loss of photoreceptors are one of the major causes of blindness. Photoreceptor replacement therapy is a promising treatment because the transplantation of retina-derived photoreceptors can be applied now to different murine retinopathies to restore visual function. To have an unlimited source of photoreceptors, we derived a transgenic embryonic stem cell (ESC) line in which the Crx-GFP transgene is expressed in photoreceptors and assessed the capacity of a 3D culture protocol to produce integration-competent photoreceptors. This culture system allows the production of a large number of photoreceptors recapitulating the in vivo development. After transplantation, integrated cells showed the typical morphology of mature rods bearing external segments and ribbon synapses. We conclude that a 3D protocol coupled with ESCs provides a safe and renewable source of photoreceptors displaying a development and transplantation competence comparable to photoreceptors from age-matched retinas.

Basin-internal derivation of hydrocarbons in the Witwatersrand Basin, South Africa: evidence from bulk and molecular δ13C data

Gold in the quartz-pebble conglomerates of the late Archean Witwatersrand Basin, South Africa, is often intimately associated with carbonaceous matter of organic/biogenic origin which occurs in the form of stratiform carbon seams and paragenetically late bitumen nodules. Both carbon forms are believed to be formed by solidification of migrating hydrocarbons. This paper presents bulk and molecular chemical and stable carbon isotope data for the carbonaceous matter, all of which are used to provide a clue to the source of the hydrocarbons. These data are compared with those from intra-basinal shales and overlying dolostone of the Transvaal Supergroup. The delta C-13 values of the extracts from the Witwatersrand carbonaceous material show small differences (up to 2.4 parts per thousand) compared to the associated insoluble organic matter. This suggests that the auriferous rocks were stained by mobile hydrocarbons produced by thermal and oxidative alteration of indigenous bitumens, a contribution from hydrocarbons derived from intra-basinal Witwatersrand shales cannot be excluded. Individual aliphatic hydrocarbons of the various carbonaceous materials were subjected to compound specific isotope analysis using on-line gas chromatography/combustion/stable isotope ratio mass spectrometry (GC/C/IRMS). The limited variability of the molecular parameters and uniform delta C-13 values of individual n-alkanes (-31.1 +/- 1.7 parts per thousand) and isoprenoids (-30.7 +/- 1.1 parts per thousand) in the Witwatersrand samples exclude the mixing of oils from different sources. Carbonaceous matter in the dolostones shows distinctly different bulk and molecular isotope characteristics and thus cannot have been the source of the hydrocarbons in the Witwatersrand deposits. All the various forms of Witwatersrand carbon appear indigenous to the Witwatersrand Basin, and the differences between them are explained by variable, in general probably short (centimeter- to meter-scale) hydrocarbon migration during diagenesis and subsequent hydrothermal infiltration. (C) 2001 Elsevier Science B.V. All rights reserved.

Robust Estimators of the Generalized Log-Gamma Distribution

We propose robust estimators of the generalized log-gamma distribution and, more generally, of location-shape-scale families of distributions. A (weighted) Q tau estimator minimizes a tau scale of the differences between empirical and theoretical quantiles. It is n(1/2) consistent; unfortunately, it is not asymptotically normal and, therefore, inconvenient for inference. However, it is a convenient starting point for a one-step weighted likelihood estimator, where the weights are based on a disparity measure between the model density and a kernel density estimate. The one-step weighted likelihood estimator is asymptotically normal and fully efficient under the model. It is also highly robust under outlier contamination. Supplementary materials are available online.

PRPF mutations are associated with generalized defects in spliceosome formation and pre-mRNA splicing in patients with retinitis pigmentosa.

Proteins PRPF31, PRPF3 and PRPF8 (RP-PRPFs) are ubiquitously expressed components of the spliceosome, a macromolecular complex that processes nearly all pre-mRNAs. Although these spliceosomal proteins are conserved in eukaryotes and are essential for survival, heterozygous mutations in human RP-PRPF genes lead to retinitis pigmentosa, a hereditary disease restricted to the eye. Using cells from patients with 10 different mutations, we show that all clinically relevant RP-PRPF defects affect the stoichiometry of spliceosomal small nuclear RNAs (snRNAs), the protein composition of tri-small nuclear ribonucleoproteins and the kinetics of spliceosome assembly. These mutations cause inefficient splicing in vitro and affect constitutive splicing ex-vivo by impairing the removal of at least 9% of endogenously expressed introns. Alternative splicing choices are also affected when RP-PRPF defects are present. Furthermore, we show that the steady-state levels of snRNAs and processed pre-mRNAs are highest in the retina, indicating a particularly elevated splicing activity. Our results suggest a role for PRPFs defects in the etiology of PRPF-linked retinitis pigmentosa, which appears to be a truly systemic splicing disease. Although these mutations cause widespread and important splicing defects, they are likely tolerated by the majority of human tissues but are critical for retinal cell survival.