Assessment of transcript reconstruction methods for RNA-seq

We evaluated 25 protocol variants of 14 independent computational methods for exon identification, transcript reconstruction and expression-level quantification from RNA-seq data. Our results show that most algorithms are able to identify discrete transcript components with high success rates but that assembly of complete isoform structures poses a major challenge even when all constituent elements are identified. Expression-level estimates also varied widely across methods, even when based on similar transcript models. Consequently, the complexity of higher eukaryotic genomes imposes severe limitations on transcript recall and splice product discrimination that are likely to remain limiting factors for the analysis of current-generation RNA-seq data.This work was supported by European Molecular Biology Laboratory, US National Institutes of Health/NHGRI grants U54HG004555 and U54HG004557, Wellcome Trust grant WT098051, and grants BIO2011-26205 and CSD2007-00050 from the Ministerio de Educación y Ciencia

Genome-Wide Survey for Biologically Functional Pseudogenes

According to current estimates there exist about 20,000 pseudogenes in a mammalian genome. The vast majority of these are disabled and nonfunctional copies of protein-coding genes which, therefore, evolve neutrally. Recent findings that a Makorin1 pseudogene, residing on mouse Chromosome 5, is, indeed, in vivo vital and also evolutionarily preserved, encouraged us to conduct a genome-wide survey for other functional pseudogenes in human, mouse, and chimpanzee. We identify to our knowledge the first examples of conserved pseudogenes common to human and mouse, originating from one duplication predating the human–mouse species split and having evolved as pseudogenes since the species split. Functionality is one possible way to explain the apparently contradictory properties of such pseudogene pairs, i.e., high conservation and ancient origin. The hypothesis of functionality is tested by comparing expression evidence and synteny of the candidates with proper test sets. The tests suggest potential biological function. Our candidate set includes a small set of long-lived pseudogenes whose unknown potential function is retained since before the human–mouse species split, and also a larger group of primate-specific ones found from human–chimpanzee searches. Two processed sequences are notable, their conservation since the human–mouse split being as high as most protein-coding genes; one is derived from the protein Ataxin 7-like 3 (ATX7NL3), and one from the Spinocerebellar ataxia type 1 protein (ATX1). Our approach is comparative and can be applied to any pair of species. It is implemented by a semi-automated pipeline based on cross-species BLAST comparisons and maximum-likelihood phylogeny estimations. To separate pseudogenes from protein-coding genes, we use standard methods, utilizing in-frame disablements, as well as a probabilistic filter based on Ka/Ks ratios

Performance of the CMS Cathode Strip Chambers with Cosmic Rays

The Cathode Strip Chambers (CSCs) constitute the primary muon tracking device in the CMS endcaps. Their performance has been evaluated using data taken during a cosmic ray run in fall 2008. Measured noise levels are low, with the number of noisy channels well below 1%. Coordinate resolution was measured for all types of chambers, and fall in the range 47 microns to 243 microns. The efficiencies for local charged track triggers, for hit and for segments reconstruction were measured, and are above 99%. The timing resolution per layer is approximately 5 ns

Analysis of RNA splicing defects in PITX2 mutants supports a gene dosage model of Axenfeld-Rieger syndrome

BACKGROUND: Axenfeld-Rieger syndrome (ARS) is associated with mutations in the PITX2 gene that encodes a homeobox transcription factor. Several intronic PITX2 mutations have been reported in Axenfeld-Rieger patients but their effects on gene expression have not been tested. METHODS: We present two new families with recurrent PITX2 intronic mutations and use PITX2c minigenes and transfected cells to address the hypothesis that intronic mutations effect RNA splicing. Three PITX2 mutations have been analyzed: a G>T mutation within the AG 3' splice site (ss) junction associated with exon 4 (IVS4-1G>T), a G>C mutation at position +5 of the 5' (ss) of exon 4 (IVS4+5G>C), and a previously reported A>G substitution at position -11 of 3'ss of exon 5 (IVS5-11A>G). RESULTS: Mutation IVS4+5G>C showed 71% retention of the intron between exons 4 and 5, and poorly expressed protein. Wild-type protein levels were proportionally expressed from correctly spliced mRNA. The G>T mutation within the exon 4 AG 3'ss junction shifted splicing exclusively to a new AG and resulted in a severely truncated, poorly expressed protein. Finally, the A>G substitution at position -11 of the 3'ss of exon 5 shifted splicing exclusively to a newly created upstream AG and resulted in generation of a protein with a truncated homeodomain. CONCLUSION: This is the first direct evidence to support aberrant RNA splicing as the mechanism underlying the disorder in some patients and suggests that the magnitude of the splicing defect may contribute to the variability of ARS phenotypes, in support of a gene dosage model of Axenfeld-Rieger syndrome

Chemical Defense by the Native Winter Ant (Prenolepis imparis) against the Invasive Argentine Ant (Linepithema humile)

The invasive Argentine ant (Linepithema humile) is established worldwide and displaces native ant species. In northern California, however, the native winter ant (Prenolepis imparis) persists in invaded areas. We found that in aggressive interactions between the two species, P. imparis employs a potent defensive secretion. Field observations were conducted at P. imparis nest sites both in the presence and absence of L. humile. These observations suggested and laboratory assays confirmed that P. imparis workers are more likely to secrete when outnumbered by L. humile. Workers of P. imparis were also more likely to secrete near their nest entrances than when foraging on trees. One-on-one laboratory trials showed that the P. imparis secretion is highly lethal to L. humile, causing 79% mortality. The nonpolar fraction of the secretion was chemically analyzed with gas chromatography/mass spectrometry, and found to be composed of long-chain and cyclic hydrocarbons. Chemical analysis of dissected P. imparis workers showed that the nonpolar fraction is derived from the Dufour's gland. Based on these conclusions, we hypothesize that this chemical defense may help P. imparis to resist displacement by L. humile

Genetic and Chemical Modifiers of a CUG Toxicity Model in Drosophila

Non-coding CUG repeat expansions interfere with the activity of human Muscleblind-like (MBNL) proteins contributing to myotonic dystrophy 1 (DM1). To understand this toxic RNA gain-of-function mechanism we developed a Drosophila model expressing 60 pure and 480 interrupted CUG repeats in the context of a non-translatable RNA. These flies reproduced aspects of the DM1 pathology, most notably nuclear accumulation of CUG transcripts, muscle degeneration, splicing misregulation, and diminished Muscleblind function in vivo. Reduced Muscleblind activity was evident from the sensitivity of CUG-induced phenotypes to a decrease in muscleblind genetic dosage and rescue by MBNL1 expression, and further supported by the co-localization of Muscleblind and CUG repeat RNA in ribonuclear foci. Targeted expression of CUG repeats to the developing eye and brain mushroom bodies was toxic leading to rough eyes and semilethality, respectively. These phenotypes were utilized to identify genetic and chemical modifiers of the CUG-induced toxicity. 15 genetic modifiers of the rough eye phenotype were isolated. These genes identify putative cellular processes unknown to be altered by CUG repeat RNA, and they include mRNA export factor Aly, apoptosis inhibitor Thread, chromatin remodelling factor Nurf-38, and extracellular matrix structural component Viking. Ten chemical compounds suppressed the semilethal phenotype. These compounds significantly improved viability of CUG expressing flies and included non-steroidal anti-inflammatory agents (ketoprofen), muscarinic, cholinergic and histamine receptor inhibitors (orphenadrine), and drugs that can affect sodium and calcium metabolism such as clenbuterol and spironolactone. These findings provide new insights into the DM1 phenotype, and suggest novel candidates for DM1 treatments

Patterns of Sequence Divergence and Evolution of the S1 Orthologous Regions between Asian and African Cultivated Rice Species

A strong postzygotic reproductive barrier separates the recently diverged Asian and African cultivated rice species, Oryza sativa and O. glaberrima. Recently a model of genetic incompatibilities between three adjacent loci: S1A, S1 and S1B (called together the S1 regions) interacting epistatically, was postulated to cause the allelic elimination of female gametes in interspecific hybrids. Two candidate factors for the S1 locus (including a putative F-box gene) were proposed, but candidates for S1A and S1B remained undetermined. Here, to better understand the basis of the evolution of regions involved in reproductive isolation, we studied the genic and structural changes accumulated in the S1 regions between orthologous sequences. First, we established an 813 kb genomic sequence in O. glaberrima, covering completely the S1A, S1 and the majority of the S1B regions, and compared it with the orthologous regions of O. sativa. An overall strong structural conservation was observed, with the exception of three isolated regions of disturbed collinearity: (1) a local invasion of transposable elements around a putative F-box gene within S1, (2) the multiple duplication and subsequent divergence of the same F-box gene within S1A, (3) an interspecific chromosomal inversion in S1B, which restricts recombination in our O. sativa×O. glaberrima crosses. Beside these few structural variations, a uniform conservative pattern of coding sequence divergence was found all along the S1 regions. Hence, the S1 regions have undergone no drastic variation in their recent divergence and evolution between O. sativa and O. glaberrima, suggesting that a small accumulation of genic changes, following a Bateson-Dobzhansky-Muller (BDM) model, might be involved in the establishment of the sterility barrier. In this context, genetic incompatibilities involving the duplicated F-box genes as putative candidates, and a possible strengthening step involving the chromosomal inversion might participate to the reproductive barrier between Asian and African rice species

Transcriptome Analysis Reveals Strain-Specific and Conserved Stemness Genes in Schmidtea mediterranea

The planarian Schmidtea mediterranea is a powerful model organism for studying stem cell biology due to its extraordinary regenerative ability mediated by neoblasts, a population of adult somatic stem cells. Elucidation of the S. mediterranea transcriptome and the dynamics of transcript expression will increase our understanding of the gene regulatory programs that regulate stem cell function and differentiation. Here, we have used RNA-Seq to characterize the S. mediterranea transcriptome in sexual and asexual animals and in purified neoblast and differentiated cell populations. Our analysis identified many uncharacterized genes, transcripts, and alternatively spliced isoforms that are differentially expressed in a strain or cell type-specific manner. Transcriptome profiling of purified neoblasts and differentiated cells identified neoblast-enriched transcripts, many of which likely play important roles in regeneration and stem cell function. Strikingly, many of the neoblast-enriched genes are orthologs of genes whose expression is enriched in human embryonic stem cells, suggesting that a core set of genes that regulate stem cell function are conserved across metazoan species

Identification and Classification of Conserved RNA Secondary Structures in the Human Genome

The discoveries of microRNAs and riboswitches, among others, have shown functional RNAs to be biologically more important and genomically more prevalent than previously anticipated. We have developed a general comparative genomics method based on phylogenetic stochastic context-free grammars for identifying functional RNAs encoded in the human genome and used it to survey an eight-way genome-wide alignment of the human, chimpanzee, mouse, rat, dog, chicken, zebra-fish, and puffer-fish genomes for deeply conserved functional RNAs. At a loose threshold for acceptance, this search resulted in a set of 48,479 candidate RNA structures. This screen finds a large number of known functional RNAs, including 195 miRNAs, 62 histone 3′UTR stem loops, and various types of known genetic recoding elements. Among the highest-scoring new predictions are 169 new miRNA candidates, as well as new candidate selenocysteine insertion sites, RNA editing hairpins, RNAs involved in transcript auto regulation, and many folds that form singletons or small functional RNA families of completely unknown function. While the rate of false positives in the overall set is difficult to estimate and is likely to be substantial, the results nevertheless provide evidence for many new human functional RNAs and present specific predictions to facilitate their further characterization