A comparison of methods for differential expression analysis of RNA-seq data.

BACKGROUND: Finding genes that are differentially expressed between conditions is an integral part of understanding the molecular basis of phenotypic variation. In the past decades, DNA microarrays have been used extensively to quantify the abundance of mRNA corresponding to different genes, and more recently high-throughput sequencing of cDNA (RNA-seq) has emerged as a powerful competitor. As the cost of sequencing decreases, it is conceivable that the use of RNA-seq for differential expression analysis will increase rapidly. To exploit the possibilities and address the challenges posed by this relatively new type of data, a number of software packages have been developed especially for differential expression analysis of RNA-seq data. RESULTS: We conducted an extensive comparison of eleven methods for differential expression analysis of RNA-seq data. All methods are freely available within the R framework and take as input a matrix of counts, i.e. the number of reads mapping to each genomic feature of interest in each of a number of samples. We evaluate the methods based on both simulated data and real RNA-seq data. CONCLUSIONS: Very small sample sizes, which are still common in RNA-seq experiments, impose problems for all evaluated methods and any results obtained under such conditions should be interpreted with caution. For larger sample sizes, the methods combining a variance-stabilizing transformation with the 'limma' method for differential expression analysis perform well under many different conditions, as does the nonparametric SAMseq method.

Transmission of HIV-1 drug resistance in Switzerland: a 10-year molecular epidemiology survey.

OBJECTIVES: Representative prevalence data of transmitted drug-resistant HIV-1 are essential to establish accurate guidelines addressing resistance testing and first-line treatments. METHODS: Systematic resistance testing was carried out in individuals in Switzerland with documented HIV-1 seroconversion during 1996-2005 and available samples with RNA > 1000 copies/ml obtained within 1 year of estimated seroconversion. Resistance interpretation used the Stanford list of mutations for surveillance of transmitted drug resistance and the French National Agency for AIDS Research algorithm. RESULTS: Viral sequences from 822 individuals were available. Risk groups were men having sex with men (42%), heterosexual contacts (32%) and intravenous drug users (20%); 30% were infected with non-B subtype viruses. Overall, prevalence of transmitted resistance was 7.7% [95% confidence interval (CI), 5.9-9.5] for any drug, 5.5% (95% CI, 3.9-7.1) for nucleoside reverse transcriptase inhibitors, 1.9% (95% CI, 1.0-2.8) for non-nucleoside reverse transcriptase inhibitors and 2.7% (95% CI, 1.6-3.8) for protease inhibitors. Dual- or triple-class resistance was observed in 2% (95% CI, 0.8-2.5). No significant trend in prevalence of transmitted resistance was observed over years. There were no differences according to ethnicity, risk groups or gender, but prevalence of transmitted resistance was highest among individuals infected with subtype B virus. CONCLUSIONS: The transmission rate of drug-resistant HIV-1 has been stable since 1996, with very rare transmission of dual- or triple-class resistance. These data suggest that transmission of drug resistance in the setting of easy access to antiretroviral treatment can remain stable and be kept at a low level.

A system to measure the kinematics during the entire ski jump sequence using inertial sensors.

Three-dimensional analysis of the entire sequence in ski jumping is recommended when studying the kinematics or evaluating performance. Camera-based systems which allow three-dimensional kinematics measurement are complex to set-up and require extensive post-processing, usually limiting ski jumping analyses to small numbers of jumps. In this study, a simple method using a wearable inertial sensors-based system is described to measure the orientation of the lower-body segments (sacrum, thighs, shanks) and skis during the entire jump sequence. This new method combines the fusion of inertial signals and biomechanical constraints of ski jumping. Its performance was evaluated in terms of validity and sensitivity to different performances based on 22 athletes monitored during daily training. The validity of the method was assessed by comparing the inclination of the ski and the slope at landing point and reported an error of -0.2±4.8°. The validity was also assessed by comparison of characteristic angles obtained with the proposed system and reference values in the literature; the differences were smaller than 6° for 75% of the angles and smaller than 15° for 90% of the angles. The sensitivity to different performances was evaluated by comparing the angles between two groups of athletes with different jump lengths and by assessing the association between angles and jump lengths. The differences of technique observed between athletes and the associations with jumps length agreed with the literature. In conclusion, these results suggest that this system is a promising tool for a generalization of three-dimensional kinematics analysis in ski jumping.

The murine interleukin-2 receptor gamma chain gene: organization, chromosomal localization and expression in the adult thymus.

Defects in the interleukin-2 receptor gamma (IL-2R gamma) chain in the man result in an X-linked severe combined immunodeficiency, SCIDX1, characterized by an absence of T-cell differentiation. This phenotype may result from pertubations in IL-2, IL-4-, IL-7- or IL-15-mediated signaling, as the IL-2R gamma chain forms an integral component of these receptor systems. We have isolated and characterized cDNA and genomic clones for the murine IL-2R gamma. The gene (Il2rg) is well conserved between mouse and man with respect to overall structure and size, and contains regions of high conservation in the promoter region as well. Il2rg maps to mouse X chromosome region 40, in a region of synteny with human Xq12-13.1. We have also explored the expression of the IL-2R gamma during thymocyte development. IL-2R gamma transcripts are detected in the earliest thymocyte precursor cells and persist throughout intrathymic development into the mature peripheral compartment. Genomic clones for the murine IL-2R gamma will allow for further studies on the regulation and function of this gene in vivo.

uvbyHbeta photometry of main sequence A type stars.

We present Stroemgren uvby and Hbeta_ photometry for a set of 575 northern main sequence A type stars, most of them belonging to the Hipparcos Input Catalogue, with V from 5mag to 10mag and with known radial velocities. These observations enlarge the catalogue we began to compile some years ago to more than 1500 stars. Our catalogue includes kinematic and astrophysical data for each star. Our future goal is to perform an accurate analysis of the kinematical behaviour of these stars in the solar neighbourhood.

Early onset collagen VI myopathies: Genetic and clinical correlations.

OBJECTIVE: Mutations in the genes encoding the extracellular matrix protein collagen VI (ColVI) cause a spectrum of disorders with variable inheritance including Ullrich congenital muscular dystrophy, Bethlem myopathy, and intermediate phenotypes. We extensively characterized, at the clinical, cellular, and molecular levels, 49 patients with onset in the first 2 years of life to investigate genotype-phenotype correlations. METHODS: Patients were classified into 3 groups: early-severe (18%), moderate-progressive (53%), and mild (29%). ColVI secretion was analyzed in patient-derived skin fibroblasts. Chain-specific transcript levels were quantified by quantitative reverse transcriptase polymerase chain reaction (qRT-PCR), and mutation identification was performed by sequencing of complementary DNA. RESULTS: ColVI secretion was altered in all fibroblast cultures studied. We identified 56 mutations, mostly novel and private. Dominant de novo mutations were detected in 61% of the cases. Importantly, mutations causing premature termination codons (PTCs) or in-frame insertions strikingly destabilized the corresponding transcripts. Homozygous PTC-causing mutations in the triple helix domains led to the most severe phenotypes (ambulation never achieved), whereas dominant de novo in-frame exon skipping and glycine missense mutations were identified in patients of the moderate-progressive group (loss of ambulation). INTERPRETATION: This work emphasizes that the diagnosis of early onset ColVI myopathies is arduous and time-consuming, and demonstrates that quantitative RT-PCR is a helpful tool for the identification of some mutation-bearing genes. Moreover, the clinical classification proposed allowed genotype-phenotype relationships to be explored, and may be useful in the design of future clinical trials.

Identification of a mammalian H(+)-myo-inositol symporter expressed predominantly in the brain.

Inositol and its phosphorylated derivatives play a major role in brain function, either as osmolytes, second messengers or regulators of vesicle endo- and exocytosis. Here we describe the identification and functional characterization of a novel H(+)-myo- inositol co-transporter, HMIT, expressed predominantly in the brain. HMIT cDNA encodes a 618 amino acid polypeptide with 12 predicted transmembrane domains. Functional expression of HMIT in Xenopus oocytes showed that transport activity was specific for myo-inositol and related stereoisomers with a Michaelis-Menten constant of approximately 100 microM, and that transport activity was strongly stimulated by decreasing pH. Electrophysiological measurements revealed that transport was electrogenic with a maximal transport activity reached at pH 5.0. In rat brain membrane preparations, HMIT appeared as a 75-90 kDa protein that could be converted to a 67 kDa band upon enzymatic deglycosylation. Immunofluorescence microscopy analysis showed HMIT expression in glial cells and some neurons. These data provide the first characterization of a mammalian H(+)-coupled myo- inositol transporter. Predominant central expression of HMIT suggests that it has a key role in the control of myo-inositol brain metabolism.

RNA-based gene duplication: mechanistic and evolutionary insights.

Gene copies that stem from the mRNAs of parental source genes have long been viewed as evolutionary dead-ends with little biological relevance. Here we review a range of recent studies that have unveiled a significant number of functional retroposed gene copies in both mammalian and some non-mammalian genomes. These studies have not only revealed previously unknown mechanisms for the emergence of new genes and their functions but have also provided fascinating general insights into molecular and evolutionary processes that have shaped genomes. For example, analyses of chromosomal gene movement patterns via RNA-based gene duplication have shed fresh light on the evolutionary origin and biology of our sex chromosomes.

Leishmania RNA virus controls the severity of mucocutaneous leishmaniasis.

Mucocutaneous leishmaniasis is caused by infections with intracellular parasites of the Leishmania Viannia subgenus, including Leishmania guyanensis. The pathology develops after parasite dissemination to nasopharyngeal tissues, where destructive metastatic lesions form with chronic inflammation. Currently, the mechanisms involved in lesion development are poorly understood. Here we show that metastasizing parasites have a high Leishmania RNA virus-1 (LRV1) burden that is recognized by the host Toll-like receptor 3 (TLR3) to induce proinflammatory cytokines and chemokines. Paradoxically, these TLR3-mediated immune responses rendered mice more susceptible to infection, and the animals developed an increased footpad swelling and parasitemia. Thus, LRV1 in the metastasizing parasites subverted the host immune response to Leishmania and promoted parasite persistence.

Induction of the alternative NF-κB pathway by lymphotoxin αβ (LTαβ) relies on internalization of LTβ receptor.

Several tumor necrosis factor receptor (TNFR) family members activate both the classical and the alternative NF-κB pathways. However, how a single receptor engages these two distinct pathways is still poorly understood. Using lymphotoxin β receptor (LTβR) as a prototype, we showed that activation of the alternative, but not the classical, NF-κB pathway relied on internalization of the receptor. Further molecular analyses revealed a specific cytosolic region of LTβR essential for its internalization, TRAF3 recruitment, and p100 processing. Interestingly, we found that dynamin-dependent, but clathrin-independent, internalization of LTβR appeared to be required for the activation of the alternative, but not the classical, NF-κB pathway. In vivo, ligand-induced internalization of LTβR in mesenteric lymph node stromal cells correlated with induction of alternative NF-κB target genes. Thus, our data shed light on LTβR cellular trafficking as a process required for specific biological functions of NF-κB.

Tandem repeat sequence variation as causative cis-eQTLs for protein-coding gene expression variation: the case of CSTB.

Association studies have revealed expression quantitative trait loci (eQTLs) for a large number of genes. However, the causative variants that regulate gene expression levels are generally unknown. We hypothesized that copy-number variation of sequence repeats contribute to the expression variation of some genes. Our laboratory has previously identified that the rare expansion of a repeat c.-174CGGGGCGGGGCG in the promoter region of the CSTB gene causes a silencing of the gene, resulting in progressive myoclonus epilepsy. Here, we genotyped the repeat length and quantified CSTB expression by quantitative real-time polymerase chain reaction in 173 lymphoblastoid cell lines (LCLs) and fibroblast samples from the GenCord collection. The majority of alleles contain either two or three copies of this repeat. Independent analysis revealed that the c.-174CGGGGCGGGGCG repeat length is strongly associated with CSTB expression (P = 3.14 × 10(-11)) in LCLs only. Examination of both genotyped and imputed single-nucleotide polymorphisms (SNPs) within 2 Mb of CSTB revealed that the dodecamer repeat represents the strongest cis-eQTL for CSTB in LCLs. We conclude that the common two or three copy variation is likely the causative cis-eQTL for CSTB expression variation. More broadly, we propose that polymorphic tandem repeats may represent the causative variation of a fraction of cis-eQTLs in the genome.

Postglacial recolonization at a snail's pace (Trochulus villosus): confronting competing refugia hypotheses using model selection.

The localization of Last Glacial Maximum (LGM) refugia is crucial information to understand a species' history and predict its reaction to future climate changes. However, many phylogeographical studies often lack sampling designs intensive enough to precisely localize these refugia. The hairy land snail Trochulus villosus has a small range centred on Switzerland, which could be intensively covered by sampling 455 individuals from 52 populations. Based on mitochondrial DNA sequences (COI and 16S), we identified two divergent lineages with distinct geographical distributions. Bayesian skyline plots suggested that both lineages expanded at the end of the LGM. To find where the origin populations were located, we applied the principles of ancestral character reconstruction and identified a candidate refugium for each mtDNA lineage: the French Jura and Central Switzerland, both ice-free during the LGM. Additional refugia, however, could not be excluded, as suggested by the microsatellite analysis of a population subset. Modelling the LGM niche of T. villosus, we showed that suitable climatic conditions were expected in the inferred refugia, but potentially also in the nunataks of the alpine ice shield. In a model selection approach, we compared several alternative recolonization scenarios by estimating the Akaike information criterion for their respective maximum-likelihood migration rates. The 'two refugia' scenario received by far the best support given the distribution of genetic diversity in T. villosus populations. Provided that fine-scale sampling designs and various analytical approaches are combined, it is possible to refine our necessary understanding of species responses to environmental changes.

The tomato genome sequence providies insights into fleshy fruit evolution

Tomato (Solanum lycopersicum) is a major crop plant and a model system for fruit development. Solanum is one of the largest angiosperm genera1 and includes annual and perennial plants from diverse habitats. Here we present a high-quality genome sequence of domesticated tomato, a draft sequence of its closest wild relative, Solanum pimpinellifolium2, and compare them to each other and to the potato genome (Solanum tuberosum). The two tomato genomes show only 0.6% nucleotide divergence and signs of recent admixture, but show more than 8% divergence from potato, with nine large and several smaller inversions. In contrast to Arabidopsis, but similar to soybean, tomato and potato small RNAs map predominantly to gene-rich chromosomal regions, including gene promoters. The Solanum lineage has experienced two consecutive genome triplications: one that is ancient and shared with rosids, and a more recent one. These triplications set the stage for the neofunctionalization of genes controlling fruit characteristics, such as colour and fleshiness.

Establishment of the epithelial-specific transcriptome of normal and malignant human breast cells based on MPSS and array expression data.

INTRODUCTION: Diverse microarray and sequencing technologies have been widely used to characterise the molecular changes in malignant epithelial cells in breast cancers. Such gene expression studies to identify markers and targets in tumour cells are, however, compromised by the cellular heterogeneity of solid breast tumours and by the lack of appropriate counterparts representing normal breast epithelial cells. METHODS: Malignant neoplastic epithelial cells from primary breast cancers and luminal and myoepithelial cells isolated from normal human breast tissue were isolated by immunomagnetic separation methods. Pools of RNA from highly enriched preparations of these cell types were subjected to expression profiling using massively parallel signature sequencing (MPSS) and four different genome wide microarray platforms. Functional related transcripts of the differential tumour epithelial transcriptome were used for gene set enrichment analysis to identify enrichment of luminal and myoepithelial type genes. Clinical pathological validation of a small number of genes was performed on tissue microarrays. RESULTS: MPSS identified 6,553 differentially expressed genes between the pool of normal luminal cells and that of primary tumours substantially enriched for epithelial cells, of which 98% were represented and 60% were confirmed by microarray profiling. Significant expression level changes between these two samples detected only by microarray technology were shown by 4,149 transcripts, resulting in a combined differential tumour epithelial transcriptome of 8,051 genes. Microarray gene signatures identified a comprehensive list of 907 and 955 transcripts whose expression differed between luminal epithelial cells and myoepithelial cells, respectively. Functional annotation and gene set enrichment analysis highlighted a group of genes related to skeletal development that were associated with the myoepithelial/basal cells and upregulated in the tumour sample. One of the most highly overexpressed genes in this category, that encoding periostin, was analysed immunohistochemically on breast cancer tissue microarrays and its expression in neoplastic cells correlated with poor outcome in a cohort of poor prognosis estrogen receptor-positive tumours. CONCLUSION: Using highly enriched cell populations in combination with multiplatform gene expression profiling studies, a comprehensive analysis of molecular changes between the normal and malignant breast tissue was established. This study provides a basis for the identification of novel and potentially important targets for diagnosis, prognosis and therapy in breast cancer.