Lee E. Vandivier

Regulatory Impact of RNA Secondary Structure across the Arabidopsis Transcriptome

This study employs a high-throughput, sequencing-based, structure-mapping approach to investigate RNA secondary structure throughout the Arabidopsis transcriptome. By combining these structural insights with a number of other RNA sequencing–based approaches, it provides a global assessment of RNA folding and its significant regulatory effects in a plant transcriptome. The secondary structure of an RNA molecule plays an integral role in its maturation, regulation, and function. However, the global influence of this feature on plant gene expression is still largely unclear. Here, we use a high-throughput, sequencing-based, structure-mapping approach in conjunction with transcriptome-wide sequencing of rRNA-depleted (RNA sequencing), small RNA, and ribosome-bound RNA populations to investigate the impact of RNA secondary structure on gene expression regulation in Arabidopsis thaliana. From this analysis, we find that highly unpaired and paired RNAs are strongly correlated with euchromatic and heterochromatic epigenetic histone modifications, respectively, providing evidence that secondary structure is necessary for these RNA-mediated posttranscriptional regulatory pathways. Additionally, we uncover key structural patterns across protein-coding transcripts that indicate RNA folding demarcates regions of protein translation and likely affects microRNA-mediated regulation of mRNAs in this model plant. We further reveal that RNA folding is significantly anticorrelated with overall transcript abundance, which is often due to the increased propensity of highly structured mRNAs to be degraded and/or processed into small RNAs. Finally, we find that secondary structure affects mRNA translation, suggesting that this feature regulates plant gene expression at multiple levels. These findings provide a global assessment of RNA folding and its significant regulatory effects in a plant transcriptome.

Fungi of the murine gut: episodic variation and proliferation during antibiotic treatment.

Antibiotic use in humans has been associated with outgrowth of fungi. Here we used a murine model to investigate the gut microbiome over 76 days of treatment with vancomycin, ampicillin, neomycin, and metronidazole and subsequent recovery. Mouse stool was studied as a surrogate for the microbiota of the lower gastrointestinal tract. The abundance of fungi and bacteria was measured using quantitative PCR, and the proportional composition of the communities quantified using 454/Roche pyrosequencing of rRNA gene tags. Prior to treatment, bacteria outnumbered fungi by >3 orders of magnitude. Upon antibiotic treatment, bacteria dropped in abundance >3 orders of magnitude, so that the predominant 16S sequences detected became transients derived from food. Upon cessation of treatment, bacterial communities mostly returned to their previous numbers and types after 8 weeks, though communities remained detectably different from untreated controls. Fungal communities varied substantially over time, even in the untreated controls. Separate cages within the same treatment group showed radical differences, but mice within a cage generally behaved similarly. Fungi increased ∼40-fold in abundance upon antibiotic treatment but declined back to their original abundance after cessation of treatment. At the last time point, Candida remained more abundant than prior to treatment. These data show that 1) gut fungal populations change radically during normal mouse husbandry, 2) fungi grow out in the gut upon suppression of bacterial communities with antibiotics, and 3) perturbations due to antibiotics persist long term in both the fungal and bacterial microbiota.

Succession in the gut microbiome following antibiotic and antibody therapies for Clostridium difficile.

Antibiotic disruption of the intestinal microbiota may cause susceptibility to pathogens that is resolved by progressive bacterial outgrowth and colonization. Succession is central to ecological theory but not widely documented in studies of the vertebrate microbiome. Here, we study succession in the hamster gut after treatment with antibiotics and exposure to Clostridium difficile. C. difficile infection is typically lethal in hamsters, but protection can be conferred with neutralizing antibodies against the A and B toxins. We compare treatment with neutralizing monoclonal antibodies (mAb) to treatment with vancomycin, which prolongs the lives of animals but ultimately fails to protect them from death. We carried out longitudinal deep sequencing analysis and found distinctive waves of succession associated with each form of treatment. Clindamycin sensitization prior to infection was associated with the temporary suppression of the previously dominant Bacteroidales and the fungus Saccinobaculus in favor of Proteobacteria. In mAb-treated animals, C. difficile proliferated before joining Proteobacteria in giving way to re-expanding Bacteroidales and the fungus Wickerhamomyces. However, the Bacteroidales lineages returning by day 7 were different from those that were present initially, and they persisted for the duration of the experiment. Animals treated with vancomycin showed a different set of late-stage lineages that were dominated by Proteobacteria as well as increased disparity between the tissue-associated and luminal cecal communities. The control animals showed no change in their gut microbiota. These data thus suggest different patterns of ecological succession following antibiotic treatment and C. difficile infection.

Essential Role for Endogenous siRNAs during Meiosis in Mouse Oocytes

The RNase III enzyme DICER generates both microRNAs (miRNAs) and endogenous short interfering RNAs (endo-siRNAs). Both small RNA species silence gene expression post-transcriptionally in association with the ARGONAUTE (AGO) family of proteins. In mammals, there are four AGO proteins (AGO1-4), of which only AGO2 possesses endonucleolytic activity. siRNAs trigger endonucleolytic cleavage of target mRNAs, mediated by AGO2, whereas miRNAs cause translational repression and mRNA decay through association with any of the four AGO proteins. Dicer deletion in mouse oocytes leads to female infertility due to defects during meiosis I. Because mouse oocytes express both miRNAs and endo-siRNAs, this phenotype could be due to the absence of either class of small RNA, or both. However, we and others demonstrated that miRNA function is suppressed in mouse oocytes, which suggested that endo-siRNAs, not miRNAs, are essential for female meiosis. To determine if this was the case we generated mice that express a catalytically inactive knock-in allele of Ago2 (Ago2ADH) exclusively in oocytes and thereby disrupted the function of siRNAs. Oogenesis and hormonal response are normal in Ago2ADH oocytes, but meiotic maturation is impaired, with severe defects in spindle formation and chromosome alignment that lead to meiotic catastrophe. The transcriptome of these oocytes is widely perturbed and shows a highly significant correlation with the transcriptome of Dicer null and Ago2 null oocytes. Expression of the mouse transcript (MT), the most abundant transposable element in mouse oocytes, is increased. This study reveals that endo-siRNAs are essential during meiosis I in mouse females, demonstrating a role for endo-siRNAs in mammals.

The Msi Family of RNA-Binding Proteins Function Redundantly as Intestinal Oncoproteins

Members of the Msi family of RNA-binding proteins have recently emerged as potent oncoproteins in a range of malignancies. MSI2 is highly expressed in hematopoietic cancers, where it is required for disease maintenance. In contrast to the hematopoietic system, colorectal cancers can express both Msi family members, MSI1 and MSI2. Here, we demonstrate that, in the intestinal epithelium, Msi1 and Msi2 have analogous oncogenic effects. Further, comparison of Msi1/2-induced gene expression programs and transcriptome-wide analyses of Msi1/2-RNA-binding targets reveal significant functional overlap, including induction of the PDK-Akt-mTORC1 axis. Ultimately, we demonstrate that concomitant loss of function of both MSI family members is sufficient to abrogate the growth of human colorectal cancer cells, and Msi gene deletion inhibits tumorigenesis in several mouse models of intestinal cancer. Our findings demonstrate that MSI1 and MSI2 act as functionally redundant oncoproteins required for the ontogeny of intestinal cancers.

The Conservation and Function of RNA Secondary Structure in Plants.

RNA transcripts fold into secondary structures via intricate patterns of base pairing. These secondary structures impart catalytic, ligand binding, and scaffolding functions to a wide array of RNAs, forming a critical node of biological regulation. Among their many functions, RNA structural elements modulate epigenetic marks, alter mRNA stability and translation, regulate alternative splicing, transduce signals, and scaffold large macromolecular complexes. Thus, the study of RNA secondary structure is critical to understanding the function and regulation of RNA transcripts. Here, we review the origins, form, and function of RNA secondary structure, focusing on plants. We then provide an overview of methods for probing secondary structure, from physical methods such as X-ray crystallography and nuclear magnetic resonance (NMR) imaging to chemical and nuclease probing methods. Combining these latter methods with high-throughput sequencing has enabled them to scale across whole transcriptomes, yielding tremendous new insights into the form and function of RNA secondary structure.

Chemical Modifications Mark Alternatively Spliced and Uncapped Messenger RNAs in Arabidopsis

Global identification of RNA modifications that affect the Watson-Crick base-pairing edge across the Arabidopsis transcriptome provides insights into the many functions of these important additions. Posttranscriptional chemical modification of RNA bases is a widespread and physiologically relevant regulator of RNA maturation, stability, and function. While modifications are best characterized in short, noncoding RNAs such as tRNAs, growing evidence indicates that mRNAs and long noncoding RNAs (lncRNAs) are likewise modified. Here, we apply our high-throughput annotation of modified ribonucleotides (HAMR) pipeline to identify and classify modifications that affect Watson-Crick base pairing at three different levels of the Arabidopsis thaliana transcriptome (polyadenylated, small, and degrading RNAs). We find this type of modifications primarily within uncapped, degrading mRNAs and lncRNAs, suggesting they are the cause or consequence of RNA turnover. Additionally, modifications within stable mRNAs tend to occur in alternatively spliced introns, suggesting they regulate splicing. Furthermore, these modifications target mRNAs with coherent functions, including stress responses. Thus, our comprehensive analysis across multiple RNA classes yields insights into the functions of covalent RNA modifications in plant transcriptomes.

Transcriptome-wide measurement of plant RNA secondary structure.

RNAs fold into intricate and precise secondary structures. These structural patterns regulate multiple steps of the RNA lifecycle, while also conferring catalytic and scaffolding functions to certain transcripts. Therefore, a full understanding of RNA posttranscriptional regulation requires a comprehensive picture of secondary structure. Here, we review several high throughput sequencing-based methods to globally survey plant RNA secondary structure. These methods are more accurate than computational prediction, and more scalable than physical techniques such as crystallography. We note hurdles to reliably measuring secondary structure, including RNA-binding proteins, RNA base modifications, and intramolecular duplexes. Finally, we survey the functional knowledge that has been gleaned from each of these methods, and identify some unanswered questions that remain.

In Silico Identification of RNA Modifications from High-Throughput Sequencing Data Using HAMR

RNA molecules are often altered post-transcriptionally by the covalent modification of their nucleotides. These modifications are known to modulate the structure, function, and activity of RNAs. When reverse transcribed into cDNA during RNA sequencing library preparation, atypical (modified) ribonucleotides that affect Watson-Crick base pairing will interfere with reverse transcriptase (RT), resulting in cDNA products with mis-incorporated bases or prematurely terminated RNA products. These interactions with RT can therefore be inferred from mismatch patterns in the sequencing reads, and are distinguishable from simple base-calling errors, single-nucleotide polymorphisms (SNPs), or RNA editing sites. Here, we describe a computational protocol for the in silico identification of modified ribonucleotides from RT-based RNA-seq read-out using the High-throughput Analysis of Modified Ribonucleotides (HAMR) software. HAMR can identify these modifications transcriptome-wide with single nucleotide resolution, and also differentiate between different types of modifications to predict modification identity. Researchers can use HAMR to identify and characterize RNA modifications using RNA-seq data from a variety of common RT-based sequencing protocols such as Poly(A), total RNA-seq, and small RNA-seq.

High-Throughput Nuclease-Mediated Probing of RNA Secondary Structure in Plant Transcriptomes

Empirical measurement of RNA secondary structure is an invaluable tool that has provided a more complete understanding of the RNA life cycle and functionality of this extremely important molecule. In general, methods for probing structural information involve treating RNA with either a chemical or an enzyme that preferentially targets regions of the RNA in a single- or double-stranded conformation (ssRNA and dsRNA, respectively). Here, we describe an approach that utilizes a combination of ssRNA- and dsRNA-specific nuclease (ss- and dsRNase, respectively) treatments along with high-throughput sequencing technology to provide comprehensive and robust measurements of RNA secondary structure across entire plant transcriptomes.