Generation and analysis of a mouse intestinal metatranscriptome through Illumina based RNA-sequencing

With the advent of high through-put sequencing (HTS), the emerging science of metagenomics is transforming our understanding of the relationships of microbial communities with their environments. While metagenomics aims to catalogue the genes present in a sample through assessing which genes are actively expressed, metatranscriptomics can provide a mechanistic understanding of community inter-relationships. To achieve these goals, several challenges need to be addressed from sample preparation to sequence processing, statistical analysis and functional annotation. Here we use an inbred non-obese diabetic (NOD) mouse model in which germ-free animals were colonized with a defined mixture of eight commensal bacteria, to explore methods of RNA extraction and to develop a pipeline for the generation and analysis of metatranscriptomic data. Applying the Illumina HTS platform, we sequenced 12 NOD cecal samples prepared using multiple RNA-extraction protocols. The absence of a complete set of reference genomes necessitated a peptide-based search strategy. Up to 16% of sequence reads could be matched to a known bacterial gene. Phylogenetic analysis of the mapped ORFs revealed a distribution consistent with ribosomal RNA, the majority from Bacteroides or Clostridium species. To place these HTS data within a systems context, we mapped the relative abundance of corresponding Escherichia coli homologs onto metabolic and protein-protein interaction networks. These maps identified bacterial processes with components that were well-represented in the datasets. In summary this study highlights the potential of exploiting the economy of HTS platforms for metatranscriptomics.

Crystal structure of the yeast eIF4A-eIF4G complex: an RNA-helicase controlled by protein-protein interactions

Translation initiation factors eIF4A and eIF4G form, together with the cap-binding factor eIF4E, the eIF4F complex, which is crucial for recruiting the small ribosomal subunit to the mRNA 5' end and for subsequent scanning and searching for the start codon. eIF4A is an ATP-dependent RNA helicase whose activity is stimulated by binding to eIF4G. We report here the structure of the complex formed by yeast eIF4G's middle domain and full-length eIF4A at 2.6-A resolution. eIF4A shows an extended conformation where eIF4G holds its crucial DEAD-box sequence motifs in a productive conformation, thus explaining the stimulation of eIF4A's activity. A hitherto undescribed interaction involves the amino acid Trp-579 of eIF4G. Mutation to alanine results in decreased binding to eIF4A and a temperature-sensitive phenotype of yeast cells that carry a Trp579Ala mutation as its sole source for eIF4G. Conformational changes between eIF4A's closed and open state provide a model for its RNA-helicase activity.

THE DOMAINS OF THE CATALYTIC SUBUNIT OF THE EUKARYOTIC RNA DEGRADING EXOSOME, RRP44P, HAVE DISTINCT FUNCTIONS

The exosome is a 3’ to 5’ exoribonuclease complex that consists of ten essential subunits. In the cytoplasm, the exosome degrades mRNA in a general mRNA turnover pathway and in several mRNA surveillance pathways. In the nucleus, the exosome processes RNA precursors to form small, stable, mature RNA species, including rRNA, snRNA, and snoRNA. In addition to processing these RNAs, the nuclear exosome is also involved in degrading aberrantly processed forms of these RNAs, and others, including mRNA. The 3’ to 5’ exoribonuclease activity of the exosome is contributed by the RNB domain of the only catalytically active subunit, Rrp44p, a member of the RNase II family of enzymes. In addition to the RNB domain, Rrp44p consists of three putative RNA binding domains and has an uncharacterized N-terminus, which includes a CR3 region and PIN domain. In an effort to characterize the cellular functions of the domains of Rrp44p, this study identified a second nuclease active site in the PIN domain. Specifically, the PIN domain exhibits endoribonuclease activity in vitro and is essential for exosome function. Further analysis of the nuclease activities of Rrp44p indicate a role for the exoribonuclease activity of Rrp44p in the cytoplasmic and nuclear exosome. This work has also characterized the CR3 region of Rrp44p, a region that has not yet been characterized in any other protein. This region is needed for the majority, if not all, of the cytoplasmic exosome functions as well as for interaction with the exosome. The CR3 region, along with a histidine residue in the N-terminus of Rrp44p, may coordinate a zinc atom. Preliminary evidence supports a role for this coordination in exosome function. Further investigation, however, is needed to determine the molecular dependence of the exosome on the CR3 region of Rrp44p. Despite its initial discovery thirteen years ago, the essential function of Rrp44p, and the exosome, is not yet known. The studies presented here, however, indicate that the essential function of Rrp44p and the exosome is in the nucleus and depends on the nuclease activities of Rrp44p.

Accessing novel developmental mechanisms in the mouse by gene trapping

Genetic analysis is a powerful method for analyzing the function of specific genes in development. I sought to identify novel genes in the mouse using a genetic analysis relying on the expression pattern and phenotype of mutated genes. To this end, I have conducted a gene trap screen using the vector $\rm SA\beta geo,$ a promoterless DNA construct that encodes a fusion protein with lacZ and neomycin resistance activities. Productive integration and expression of the $\beta$geo protein in embryonic stem (ES) cells requires integration into an active transcription unit. The endogenous regulatory elements direct reporter gene expression which reflects the expression of the endogenous gene. Of eight mouse lines generated from gene trap ES cell clones, four showed differential regulation of $\beta$geo activity during embryogenesis. These four were analyzed in more detail.^ Three of the lines RNA 1, RNA2 and RNA 3 had similar expression patterns, within subsets of cells in sites of embryonic hematopoiesis. Cloning of the trapped genes revealed that all three integrations had occurred within 45S rRNA precursor transcription units. These results imply that there exists in these cells some mechanism responsible for the efficient production of the $\beta$geo protein from an RNA polymerase I transcript that is not present in most of the cells in the embryo.^ The fourth line, GT-2, showed widespread, dynamic expression. Many of the sites of expression were important classic embryonic induction systems. Cloning of the sequences fused to the $5\sp\prime$ end of the $\beta$geo sequence revealed that the trapped gene contained significant sequence homology with a previously identified human sequence HumORF5. An open reading frame of this sequence is homologous to a group of eukaryotic proteins that are members of the RNA helicase superfamily I.^ Analysis of the gene trap lines suggests that potentially novel developmental mechanisms have been uncovered. In the case of RNA 1, 2 and 3, the differential production of ribosomal RNAs may be required for differentiation or function of the $\beta$geo positive hematopoietic cells. In the GT-2 line, a previously unsuspected temporal and spatial regulation of a putative RNA helicase implies a role for this activity during specific aspects of mouse development. ^

The role of specific regions of ribosomal RNAs in translation termination

The ribosome is a molecular machine that produces proteins in a cell. It consists of RNAs (rRNAs) and proteins. The rRNAs have been implicated in various aspects of protein biosynthesis supporting the idea that they function directly in translation. In this study the direct involvement of rRNA in translation termination was hypothesized and both genetic and biochemical strategies were designed to test this hypothesis. As a result, several regions of rRNAs from both ribosomal subunits were implicated in termination. More specifically, the mutation G1093A in an RNA of the large subunit (23S rRNA) and the mutation C1054A in the small subunit RNA (16S rRNA) of the Escherichia coli ribosome, were shown to affect the binding of the proteins that drive termination, RF1 and RF2. These mutations also caused defects in catalysis of peptidyl-tRNA hydrolysis, the last step of termination. Furthermore, the mutations affected the function of RF2 to a greater extent than that of RF1, a striking result considering the similarity of the RFs. The major defect in RF2 function was consistent with in vivo characteristics of the mutants and can be explained by the inability of the mutant rRNA sites to activate the hydrolytic center, that is the catalytic site for peptidyl-tRNA hydrolysis. Consistent with this explanation is the possibility of a direct interaction between the G1093-region (domain II of 23S rRNA) and the hydrolytic center (most likely domains IV–VI of 23S rRNA). To test that interaction hypothesis selections were performed for mutations in domains IV–VI that compensated for the growth defects caused by G1093A. Several compensatory mutations were isolated which not only restored growth in the presence of G1093A but also appeared to compensate for the termination defects caused by the G1093A. Therefore these results provided genetic evidence for an intramolecular interaction that might lead to peptidyl-tRNA hydrolysis. Finally, a new approach to the study of rRNA involvement in termination was designed. By screening a library of rRNA fragments, a fragment of the 23S rRNA (nt 74-136) was identified that caused readthrough of UGA. The antisense RNA fragment produced a similar effect. The data implicated the corresponding segment of intact 23S rRNA in termination. ^

Analysis of primary Wilms tumors for WT1 RNA processing alterations

Wilms tumor (WT) is an embryonal renal tumor with a heterogeneous genetic etiology that serves as a valuable model for studying tumorigenesis. Biallelic inactivation of the tumor suppressor gene WT1, a zinc-finger transcriptional regulator located at 11p13, is critical for the development of some Wilms tumors. Interestingly, WT1 genomic analysis has demonstrated mutations in less than 20% of WT cases. This suggests either other genes play a more major role in Wilms tumorigenesis or WT1 is functionally altered by mechanisms other than DNA mutation. Previous observations in rat and in WT xenograft cell lines have suggested that abnormal WT1 RNA processing (exon 6 RNA editing and aberrant exon 2 splicing, respectively) is a potential mechanism of altering WT1 function in the absence of a WT1 DNA mutation. However, the role of this abnormal RNA processing has not previously been assessed in primary Wilms tumors. ^ To test the hypothesis that abnormal WT1 RNA processing is a mechanism of WT1alteration during tumor development, WT1 RNA from 85 primary tumors was analyzed using reverse transcription and polymerase chain reaction amplification (RT-PCR). Although no evidence for WT1 RNA editing was observed, variable levels (5% to 50%) of aberrant WT1 exon 2 splicing were detected for 11 tumors in the absence of a detectable WT1 DNA mutation. Also, alteration of normal WT1 alternative splicing, observed as RNA isoform loss, was detected in five tumors with no apparent WT1 genomic alteration, although no consistent pattern of RNA isoform loss was detected. This abnormal WT1 splicing, detected by either loss of exon 2 from some of the transcripts or loss of RNA isoforms, is statistically correlated with relapse (p = 0.005). These studies demonstrate that abnormal WT1 RNA processing is not a common mechanism of abrogating normal WT1 function in primary tumors. However, in those cases in which abnormal WTI splicing is present, these data indicate that it may serve as a useful prognostic marker for relapse in WT patients. ^

Atomic mutagenesis of the ribosomal Sarcin-Ricin-Loop to study EF-G GTPase activation

Translocation factor EF-G, possesses a low basal GTPase activity, which is stimulated by the ribosome. One potential region of the ribosome that triggers GTPase activity of EF-G is the Sarcin-Ricin-Loop (SRL) (helix 95) in domain VI of the 23S rRNA. Structural data showed that the tip of the SRL closely approaches GTP in the active center of EF-G, structural probing data confirmed that EF-G interacts with nucleotides G2655, A2660, G2661 and A2662.1-3 The exocyclic group of adenine at A2660 is required for stimulation of EF-G GTPase activity by the ribosome as demonstrated using atomic mutagenesis.4 Recent crystal structures of EF-G on the ribosome, gave more insights into the molecular mechanism of EF-G GTPase activity.5 Based on the structure of EF-Tu on the ribosome1, the following mechanism of GTPase activation was proposed: upon binding of EF-G to the ribosome, the conserved His92 (E.coli) changes its position, pointing to the γ-phosphate of GTP. In this activated state, the phosphate of residue A2662 of the SRL positions the catalytic His in its active conformation. It was further proposed that the phosphate oxygen of A2662 is involved in a charge-relay system, enabling GTP hydrolysis. In order to test this mechanism, we use the atomic mutagenesis approach, which allows introducing non-natural modifications in the SRL, in the context of the complete 70S ribosome. Therefore, we replaced one of the non-bridging oxygens of A2662 by a methyl group. A methylphosphonat is not able to position or activate a histidine, as it has no free electrons and therefore no proton acceptor function. These modified ribosomes were then tested for stimulation of EF-G GTPase activity. First experiments show that one of the two stereoisomers incorporated into ribosomes does not stimulate GTPase activity of EF-G, whereas the other is active. From this we conclude that indeed the non-bridging phosphate oxygen of A2662 is involved in EF-G GTPase activation by the ribosome. Ongoing experiments aim at revealing the contribution of this non-bridging oxygen at A2662 to the mechanism of EF-G GTPase activation at the atomic level.

tRNA-derived fragments target the small ribosomal subunit to fine-tune translation

Post-transcriptional cleavage of RNA molecules to generate smaller fragments is a widespread mechanism that enlarges the structural and functional complexity of cellular RNomes. In particular, fragments deriving from both precursor and mature tRNAs represent one of the rapidly growing classes of post-transcriptional RNA pieces. Importantly, these tRNA-derived fragments (tRFs) possess distinct expression patterns, abundance, cellular localizations, or biological roles compared with their parental tRNA molecules (1). Here we present evidence that tRFs from the archaeon Haloferax volcanii directly bind to ribosomes. In a previous genomic screen for ribosome-associated small RNAs we have identified a 26 residue long fragment originating from the 5’ part of valine tRNA (Val-tRF) to be by far the most abundant tRF in H. volcanii (2). The Val-tRF is processed in a stress- dependent manner and was found to primarily target the small ribosomal subunit in vitro and in vivo. Translational activity was markedly reduced in the presence of Val-tRF, while control RNA fragments of similar length did not show inhibition of protein biosynthesis. Crosslinking experiments and subsequent primer extension analyses revealed the Val-tRF interaction site to surround the mRNA path in the 30S subunit. In support of this, binding experiments demonstrated that Val-tRF does compete with mRNAs for ribosome binding. Therefore this tRF represents a ribosome-bound non-protein-coding RNA (ncRNA) capable of regulating gene expression in H. volcanii under environmental stress conditions probably by fine-tuning the rate of protein production (1). (1) Gebetsberger J. and Polacek N. (2013), RNA Biol. 10:1798-1808 (2) Gebetsberger J. et. al. (2012), Archaea, Article ID 260909

Ribosomal and Immune Transcripts Associate with Relapse in Acquired ADAMTS13-Deficient Thrombotic Thrombocytopenic Purpura.

Approximately 40% of patients who survive acute episodes of thrombotic thrombocytopenic purpura (TTP) associated with severe acquired ADAMTS13 deficiency experience one or more relapses. Risk factors for relapse other than severe ADAMTS13 deficiency and ADAMTS13 autoantibodies are unknown. ADAMTS13 autoantibodies, TTP episodes following infection or type I interferon treatment and reported ensuing systemic lupus erythematosus in some patients suggest immune dysregulation. This cross-sectional study asked whether autoantibodies against RNA-binding proteins or peripheral blood gene expression profiles measured during remission are associated with history of prior relapse in acquired ADAMTS13-deficient TTP. Peripheral blood from 38 well-characterized patients with autoimmune ADAMTS13-deficient TTP in remission was examined for autoantibodies and global gene expression. A subset of TTP patients (9 patients, 24%) exhibited a peripheral blood gene signature composed of elevated ribosomal transcripts that associated with prior relapse. A non-overlapping subset of TTP patients (9 patients, 24%) displayed a peripheral blood type I interferon gene signature that associated with autoantibodies to RNA-binding proteins but not with history of relapse. Patients who had relapsed bimodally expressed higher HLA transcript levels independently of ribosomal transcripts. Presence of any one potential risk factor (ribosomal gene signature, elevated HLA-DRB1, elevated HLA-DRB5) associated with relapse (OR = 38.4; p = 0.0002) more closely than any factor alone or all factors together. Levels of immune transcripts typical of natural killer (NK) and T lymphocytes positively correlated with ribosomal gene expression and number of prior episodes but not with time since the most recent episode. Flow cytometry confirmed elevated expression of cell surface markers encoded by these transcripts on T and/or NK cell subsets of patients who had relapsed. These data associate elevated ribosomal and immune transcripts with relapse history in acquired, ADAMTS13-deficient TTP.

The integrity of the G2421-C2395 base pair in the ribosomal E-site is crucial for protein synthesis.

During the elongation cycle of protein biosynthesis, tRNAs traverse through the ribosome by consecutive binding to the 3 ribosomal binding sites (A-, P-, and E- sites). While the ribosomal A- and P-sites have been functionally well characterized in the past, the contribution of the E-site to protein biosynthesis is still poorly understood in molecular terms. Previous studies suggested an important functional interaction of the terminal residue A76 of E-tRNA with the nucleobase of the universally conserved 23S rRNA residue C2394. Using an atomic mutagenesis approach to introduce non-natural nucleoside analogs into the 23S rRNA, we could show that removal of the nucleobase or the ribose 2'-OH at C2394 had no effect on protein synthesis. On the other hand, our data disclose the importance of the highly conserved E-site base pair G2421-C2395 for effective translation. Ribosomes with a disrupted G2421-C2395 base pair are defective in tRNA binding to the E-site. This results in an impaired translation of genuine mRNAs, while homo-polymeric templates are not affected. Cumulatively our data emphasize the importance of E-site tRNA occupancy and in particular the intactness of the 23S rRNA base pair G2421-C2395 for productive protein biosynthesis.

Trypanosoma brucei RRM1 is a nuclear RNA-binding protein and modulator of chromatin structure

TbRRM1 of Trypanosoma brucei is a nucleoprotein that was previously identified in a search for splicing factors in T. brucei. We show that TbRRM1 associates with mRNAs and with the auxiliary splicing factor polypyrimidine tract-binding protein 2, but not with components of the core spliceosome. TbRRM1 also interacts with several retrotransposon hot spot (RHS) proteins and histones. RNA immunoprecipitation of a tagged form of TbRRM1 from procyclic (insect) form trypanosomes identified ca. 1,500 transcripts that were enriched and 3,000 transcripts that were underrepresented compared to cellular mRNA. Enriched transcripts encoded RNA-binding proteins, including TbRRM1 itself, several RHS transcripts, mRNAs with long coding regions, and a high proportion of stage-regulated mRNAs that are more highly expressed in bloodstream forms. Transcripts encoding ribosomal proteins, other factors involved in translation, and procyclic-specific transcripts were underrepresented. Knockdown of TbRRM1 by RNA interference caused widespread changes in mRNA abundance, but these changes did not correlate with the binding of the protein to transcripts, and most splice sites were unchanged, negating a general role for TbRRM1 in splice site selection. When changes in mRNA abundance were mapped across the genome, regions with many downregulated mRNAs were identified. Two regions were analyzed by chromatin immunoprecipitation, both of which exhibited increases in nucleosome occupancy upon TbRRM1 depletion. In addition, subjecting cells to heat shock resulted in translocation of TbRRM1 to the cytoplasm and compaction of chromatin, consistent with a second role for TbRRM1 in modulating chromatin structure. IMPORTANCE: Trypanosoma brucei, the parasite that causes human sleeping sickness, is transmitted by tsetse flies. The parasite progresses through different life cycle stages in its two hosts, altering its pattern of gene expression in the process. In trypanosomes, protein-coding genes are organized as polycistronic units that are processed into monocistronic mRNAs. Since genes in the same unit can be regulated independently of each other, it is believed that gene regulation is essentially posttranscriptional. In this study, we investigated the role of a nuclear RNA-binding protein, TbRRM1, in the insect stage of the parasite. We found that TbRRM1 binds nuclear mRNAs and also affects chromatin status. Reduction of nuclear TbRRM1 by RNA interference or heat shock resulted in chromatin compaction. We propose that TbRRM1 regulates RNA polymerase II-driven gene expression both cotranscriptionally, by facilitating transcription and efficient splicing, and posttranscriptionally, via its interaction with nuclear mRNAs.

The effect of RNA base lesions on mRNA translation

The biological effect of oxidatively damaged RNA, unlike oxidatively damaged DNA, has rarely been investigated, although it poses a threat to any living cell. Here we report on the effect of the commonly known RNA base-lesions 8-oxo-rG, 8-oxo-rA, ε-rC, ε-rA, 5-HO-rC, 5-HO-rU and the RNA abasic site (rAS) on ribosomal translation. To this end we have developed an in vitro translation assay based on the mRNA display methodology. A short synthetic mRNA construct containing the base lesion in a predefined position of the open reading frame was 32P-labeled at the 5′-end and equipped with a puromycin unit at the 3′-end. Upon in vitro translation in rabbit reticulocyte lysates, the encoded peptide chain is transferred to the puromycin unit and the products analyzed by gel electrophoresis. Alternatively, the unlabeled mRNA construct was used and incubated with 35S-methionine to prove peptide elongation of the message. We find that all base-lesions interfere substantially with ribosomal translation. We identified two classes, the first containing modifications at the base coding edge (ε-rC, ε-rA and rAS) which completely abolish peptide synthesis at the site of modification, and the second consisting of 8-oxo-rG, 8-oxo-rA, 5-HO-rC and 5-HO-rU that significantly retard full-length peptide synthesis, leading to some abortive peptides at the site of modification.