Defining the commonalities between post-transcriptional and post-translational modification communities.

Post-transcriptional modifications of RNA (PRMs) and post-translational modifications of proteins (PTMs) are important regulatory mechanisms in biological processes and have many commonalities. However, the integration of these research areas is lacking. A recent discussion identified the priorities, areas of emphasis, and necessary technologies to advance and integrate these areas of study.

The covalent nucleotide modifications within plant mRNAs: What we know, how we find them, and what should be done in the future.

Although covalent nucleotide modifications were first identified on the bases of transfer RNAs (tRNAs) and ribosomal RNAs (rRNAs), a number of these epitranscriptome marks have also been found to occur on the bases of messenger RNAs (mRNAs). These covalent mRNA features have been demonstrated to have various and significant effects on the processing (e.g. splicing, polyadenylation, etc.) and functionality (e.g. translation, transport, etc.) of these protein-encoding molecules. Here, we focus our attention on the current understanding of the collection of covalent nucleotide modifications known to occur on mRNAs in plants, how they are detected and studied, and the most outstanding future questions of each of these important epitranscriptomic regulatory signals.

Pathogen-induced m6A dynamics affect plant immunity.

Posttranscriptional regulation of mRNA mediated by methylation at the N6 position of adenine (N6-methyladenosine [m6A]) has profound effects on transcriptome regulation in plants. Focused studies across eukaryotes offer glimpses into the processes governed by m6A throughout developmental and disease states. However, we lack an understanding of the dynamics and the regulatory potential of m6A during biotic stress in plants. Here, we provide a comprehensive look into the effects of m6A on both the short-term and long-term responses to pathogen signaling in Arabidopsis (Arabidopsis thaliana). We demonstrate that m6A-deficient plants are more resistant to bacterial and fungal pathogen infections and have altered immune responses. Furthermore, m6A deposition is specifically coordinated on transcripts involved in defense and immunity prior to and proceeding the pathogen signal flagellin. Consequently, the dynamic modulation of m6A on specific stress-responsive transcripts is correlated with changes in abundance and cleavage of these transcripts. Overall, we show that the m6A methylome is regulated prior to and during simulated and active pathogen stress and functions in the coordination and balancing of normal growth and pathogen responses.

The LIN28B-IMP1 post-transcriptional regulon has opposing effects on oncogenic signaling in the intestine.

RNA-binding proteins (RBPs) are expressed broadly during both development and malignant transformation, yet their mechanistic roles in epithelial homeostasis or as drivers of tumor initiation and progression are incompletely understood. Here we describe a novel interplay between RBPs LIN28B and IMP1 in intestinal epithelial cells. Ribosome profiling and RNA sequencing identified IMP1 as a principle node for gene expression regulation downstream from LIN28B In vitro and in vivo data demonstrate that epithelial IMP1 loss increases expression of WNT target genes and enhances LIN28B-mediated intestinal tumorigenesis, which was reversed when we overexpressed IMP1 independently in vivo. Furthermore, IMP1 loss in wild-type or LIN28B-overexpressing mice enhances the regenerative response to irradiation. Together, our data provide new evidence for the opposing effects of the LIN28B-IMP1 axis on post-transcriptional regulation of canonical WNT signaling, with implications in intestinal homeostasis, regeneration and tumorigenesis.

Identification and functional annotation of long intergenic non-coding RNAs in Brassicaceae.

Long intergenic noncoding RNAs (lincRNAs) are a large yet enigmatic class of eukaryotic transcripts that can have critical biological functions. The wealth of RNA-sequencing (RNA-seq) data available for plants provides the opportunity to implement a harmonized identification and annotation effort for lincRNAs that enables cross-species functional and genomic comparisons as well as prioritization of functional candidates. In this study, we processed >24 Tera base pairs of RNA-seq data from >16,000 experiments to identify ∼130,000 lincRNAs in four Brassicaceae: Arabidopsis thaliana, Camelina sativa, Brassica rapa, and Eutrema salsugineum. We used nanopore RNA-seq, transcriptome-wide structural information, peptide data, and epigenomic data to characterize these lincRNAs and identify conserved motifs. We then used comparative genomic and transcriptomic approaches to highlight lincRNAs in our data set with sequence or transcriptional conservation. Finally, we used guilt-by-association analyses to assign putative functions to lincRNAs within our data set. We tested this approach on a subset of lincRNAs associated with germination and seed development, observing germination defects for Arabidopsis lines harboring T-DNA insertions at these loci. LincRNAs with Brassicaceae-conserved putative miRNA binding motifs, small open reading frames, or abiotic-stress modulated expression are a few of the annotations that will guide functional analyses into this cryptic portion of the transcriptome.

Comparative analyses suggest a link between mRNA splicing, stability, and RNA covalent modifications in flowering plants.

BACKGROUND: In recent years, covalent modifications on RNA nucleotides have emerged as pivotal moieties influencing the structure, function, and regulatory processes of RNA Polymerase II transcripts such as mRNAs and lncRNAs. However, our understanding of their biological roles and whether these roles are conserved across eukaryotes remains limited. RESULTS: In this study, we leveraged standard polyadenylation-enriched RNA-sequencing data to identify and characterize RNA modifications that introduce base-pairing errors into cDNA reads. Our investigation incorporated data from three Poaceae (Zea mays, Sorghum bicolor, and Setaria italica), as well as publicly available data from a range of stress and genetic contexts in Sorghum and Arabidopsis thaliana. We uncovered a strong enrichment of RNA covalent modifications (RCMs) deposited on a conserved core set of nuclear mRNAs involved in photosynthesis and translation across these species. However, the cohort of modified transcripts changed based on environmental context and developmental program, a pattern that was also conserved across flowering plants. We determined that RCMs can partly explain accession-level differences in drought tolerance in Sorghum, with stress-associated genes receiving a higher level of RCMs in a drought tolerant accession. To address function, we determined that RCMs are significantly enriched near exon junctions within coding regions, suggesting an association with splicing. Intriguingly, we found that these base-pair disrupting RCMs are associated with stable mRNAs, are highly correlated with protein abundance, and thus likely associated with facilitating translation. CONCLUSIONS: Our data point to a conserved role for RCMs in mRNA stability and translation across the flowering plant lineage.

mRNA N6 -methyladenosine is critical for cold tolerance in Arabidopsis.

Plants respond to low temperatures by altering the mRNA abundance of thousands of genes contributing to numerous physiological and metabolic processes that allow them to adapt. At the post-transcriptional level, these cold stress-responsive transcripts undergo alternative splicing, microRNA-mediated regulation and alternative polyadenylation, amongst others. Recently, m6 A, m5 C and other mRNA modifications that can affect the regulation and stability of RNA were discovered, thus revealing another layer of post-transcriptional regulation that plays an important role in modulating gene expression. The importance of m6 A in plant growth and development has been appreciated, although its significance under stress conditions is still underexplored. To assess the role of m6 A modifications during cold stress responses, methylated RNA immunoprecipitation sequencing was performed in Arabidopsis seedlings esposed to low temperature stress (4°C) for 24 h. This transcriptome-wide m6 A analysis revealed large-scale shifts in this modification in response to low temperature stress. Because m6 A is known to affect transcript stability/degradation and translation, we investigated these possibilities. Interestingly, we found that cold-enriched m6 A-containing transcripts demonstrated the largest increases in transcript abundance coupled with increased ribosome occupancy under cold stress. The significance of the m6 A epitranscriptome on plant cold tolerance was further assessed using the mta mutant in which the major m6 A methyltransferase gene was mutated. Compared to the wild-type, along with the differences in CBFs and COR gene expression levels, the mta mutant exhibited hypersensitivity to cold treatment as determined by primary root growth, biomass, and reactive oxygen species accumulation. Furthermore, and most importantly, both non-acclimated and cold-acclimated mta mutant demonstrated hypersensitivity to freezing tolerance. Taken together, these findings suggest a critical role for the epitranscriptome in cold tolerance of Arabidopsis.

Highly integrated single-base resolution maps of the epigenome in Arabidopsis.

Deciphering the multiple layers of epigenetic regulation that control transcription is critical to understanding how plants develop and respond to their environment. Using sequencing-by-synthesis technology we directly sequenced the cytosine methylome (methylC-seq), transcriptome (mRNA-seq), and small RNA transcriptome (smRNA-seq) to generate highly integrated epigenome maps for wild-type Arabidopsis thaliana and mutants defective in DNA methyltransferase or demethylase activity. At single-base resolution we discovered extensive, previously undetected DNA methylation, identified the context and level of methylation at each site, and observed local sequence effects upon methylation state. Deep sequencing of smRNAs revealed a direct relationship between the location of smRNAs and DNA methylation, perturbation of smRNA biogenesis upon loss of CpG DNA methylation, and a tendency for smRNAs to direct strand-specific DNA methylation in regions of RNA-DNA homology. Finally, strand-specific mRNA-seq revealed altered transcript abundance of hundreds of genes, transposons, and unannotated intergenic transcripts upon modification of the DNA methylation state.

Global analysis of the RNA-protein interaction and RNA secondary structure landscapes of the Arabidopsis nucleus.

Posttranscriptional regulation in eukaryotes requires cis- and trans-acting features and factors including RNA secondary structure and RNA-binding proteins (RBPs). However, a comprehensive view of the structural and RBP interaction landscape of nuclear RNAs has yet to be compiled for any organism. Here, we use our ribonuclease-mediated structure and RBP-binding site mapping approaches to globally profile these features in Arabidopsis seedling nuclei in vivo. We reveal anticorrelated patterns of secondary structure and RBP binding throughout nuclear mRNAs that demarcate sites of alternative splicing and polyadenylation. We also uncover a collection of protein-bound sequence motifs, and identify their structural contexts, co-occurrences in transcripts encoding functionally related proteins, and interactions with putative RBPs. Finally, using these motifs, we find that the chloroplast RBP CP29A also interacts with nuclear mRNAs. In total, we provide a simultaneous view of the RNA secondary structure and RBP interaction landscapes in a eukaryotic nucleus.

mRNA adenosine methylase (MTA) deposits m6A on pri-miRNAs to modulate miRNA biogenesis in Arabidopsis thaliana.

In Arabidopsis thaliana, the METTL3 homolog, mRNA adenosine methylase (MTA) introduces N6-methyladenosine (m6A) into various coding and noncoding RNAs of the plant transcriptome. Here, we show that an MTA-deficient mutant (mta) has decreased levels of microRNAs (miRNAs) but accumulates primary miRNA transcripts (pri-miRNAs). Moreover, pri-miRNAs are methylated by MTA, and RNA structure probing analysis reveals a decrease in secondary structure within stem-loop regions of these transcripts in mta mutant plants. We demonstrate interaction between MTA and both RNA Polymerase II and TOUGH (TGH), a plant protein needed for early steps of miRNA biogenesis. Both MTA and TGH are necessary for efficient colocalization of the Microprocessor components Dicer-like 1 (DCL1) and Hyponastic Leaves 1 (HYL1) with RNA Polymerase II. We propose that secondary structure of miRNA precursors induced by their MTA-dependent m6A methylation status, together with direct interactions between MTA and TGH, influence the recruitment of Microprocessor to plant pri-miRNAs. Therefore, the lack of MTA in mta mutant plants disturbs pri-miRNA processing and leads to the decrease in miRNA accumulation. Furthermore, our findings reveal that reduced miR393b levels likely contributes to the impaired auxin response phenotypes of mta mutant plants.

H3K36 methylation promotes longevity by enhancing transcriptional fidelity.

Epigenetic mechanisms, including histone post-translational modifications, control longevity in diverse organisms. Relatedly, loss of proper transcriptional regulation on a global scale is an emerging phenomenon of shortened life span, but the specific mechanisms linking these observations remain to be uncovered. Here, we describe a life span screen in Saccharomyces cerevisiae that is designed to identify amino acid residues of histones that regulate yeast replicative aging. Our results reveal that lack of sustained histone H3K36 methylation is commensurate with increased cryptic transcription in a subset of genes in old cells and with shorter life span. In contrast, deletion of the K36me2/3 demethylase Rph1 increases H3K36me3 within these genes, suppresses cryptic transcript initiation, and extends life span. We show that this aging phenomenon is conserved, as cryptic transcription also increases in old worms. We propose that epigenetic misregulation in aging cells leads to loss of transcriptional precision that is detrimental to life span, and, importantly, this acceleration in aging can be reversed by restoring transcriptional fidelity.

KMT2A-MAML2 rearrangement emerged and regressed during neuroblastoma therapy without leukemia after 12.8-year follow-up.

Twelvepatients without therapy-related leukemia were studied after completing TOP2 poison chemotherapy in a high-risk neuroblastoma regimen. One patient harbored an inv(11) that was a KMT2A rearrangement. The KMT2A-MAML2 transcript was expressed at low level. The patient was prospectively followed. The inv(11) was undetectable in ensuing samples. Leukemia never developed after a 12.8-year follow-up period. Enriched etoposide-induced TOP2A cleavage in the relevant MAML2 genomic region supports a TOP2A DNA damage mechanism. After completing TOP2 poison chemotherapies, covert KMT2A-R clones may occur in a small minority of patients; however, not all KMT2A rearrangements herald a therapy-related leukemia diagnosis.

A MicroRNA precursor surveillance system in quality control of MicroRNA synthesis.

MicroRNAs (miRNAs) are essential for regulation of gene expression. Though numerous miRNAs have been identified by high-throughput sequencing, few precursor miRNAs (pre-miRNAs) are experimentally validated. Here we report a strategy for constructing high-throughput sequencing libraries enriched for full-length pre-miRNAs. We find widespread and extensive uridylation of Argonaute (Ago)-bound pre-miRNAs, which is primarily catalyzed by two terminal uridylyltransferases: TUT7 and TUT4. Uridylation by TUT7/4 not only polishes pre-miRNA 3' ends, but also facilitates their degradation by the exosome, preventing clogging of Ago with defective species. We show that the exosome exploits distinct substrate preferences of DIS3 and RRP6, its two catalytic subunits, to distinguish productive from defective pre-miRNAs. Furthermore, we identify a positive feedback loop formed by the exosome and TUT7/4 in triggering uridylation and degradation of Ago-bound pre-miRNAs. Our study reveals a pre-miRNA surveillance system that comprises TUT7, TUT4, and the exosome in quality control of miRNA synthesis.

Plant 24-nt reproductive phasiRNAs from intramolecular duplex mRNAs in diverse monocots.

In grasses, two pathways that generate diverse and numerous 21-nt (premeiotic) and 24-nt (meiotic) phased siRNAs are highly enriched in anthers, the male reproductive organs. These "phasiRNAs" are analogous to mammalian piRNAs, yet their functions and evolutionary origins remain largely unknown. The 24-nt meiotic phasiRNAs have only been described in grasses, wherein their biogenesis is dependent on a specialized Dicer (DCL5). To assess how evolution gave rise to this pathway, we examined reproductive phasiRNA pathways in nongrass monocots: garden asparagus, daylily, and lily. The common ancestors of these species diverged approximately 115-117 million years ago (MYA). We found that premeiotic 21-nt and meiotic 24-nt phasiRNAs were abundant in all three species and displayed spatial localization and temporal dynamics similar to grasses. The miR2275-triggered pathway was also present, yielding 24-nt reproductive phasiRNAs, and thus originated more than 117 MYA. In asparagus, unlike in grasses, these siRNAs are largely derived from inverted repeats (IRs); analyses in lily identified thousands of precursor loci, and many were also predicted to form foldback substrates for Dicer processing. Additionally, reproductive phasiRNAs were present in female reproductive organs and thus may function in both male and female germinal development. These data describe several distinct mechanisms of production for 24-nt meiotic phasiRNAs and provide new insights into the evolution of reproductive phasiRNA pathways in monocots.

Nuclear m6A reader YTHDC1 regulates alternative polyadenylation and splicing during mouse oocyte development.

The N6-methyladenosine (m6A) modification is the most prevalent internal RNA modification in eukaryotes. The majority of m6A sites are found in the last exon and 3' UTRs. Here we show that the nuclear m6A reader YTHDC1 is essential for embryo viability and germline development in mouse. Specifically, YTHDC1 is required for spermatogonial development in males and for oocyte growth and maturation in females; Ythdc1-deficient oocytes are blocked at the primary follicle stage. Strikingly, loss of YTHDC1 leads to extensive alternative polyadenylation in oocytes, altering 3' UTR length. Furthermore, YTHDC1 deficiency causes massive alternative splicing defects in oocytes. The majority of splicing defects in mutant oocytes are rescued by introducing wild-type, but not m6A-binding-deficient, YTHDC1. YTHDC1 is associated with the pre-mRNA 3' end processing factors CPSF6, SRSF3, and SRSF7. Thus, YTHDC1 plays a critical role in processing of pre-mRNA transcripts in the oocyte nucleus and may have similar non-redundant roles throughout fetal development.

Lamin B1 depletion in senescent cells triggers large-scale changes in gene expression and the chromatin landscape.

Senescence is a stable proliferation arrest, associated with an altered secretory pathway, thought to promote tumor suppression and tissue aging. While chromatin regulation and lamin B1 down-regulation have been implicated as senescence effectors, functional interactions between them are poorly understood. We compared genome-wide Lys4 trimethylation on histone H3 (H3K4me3) and H3K27me3 distributions between proliferating and senescent human cells and found dramatic differences in senescence, including large-scale domains of H3K4me3- and H3K27me3-enriched "mesas" and H3K27me3-depleted "canyons." Mesas form at lamin B1-associated domains (LADs) in replicative senescence and oncogene-induced senescence and overlap DNA hypomethylation regions in cancer, suggesting that pre-malignant senescent chromatin changes foreshadow epigenetic cancer changes. Hutchinson-Gilford progeria syndrome fibroblasts (mutant lamin A) also show evidence of H3K4me3 mesas, suggesting a link between premature chromatin changes and accelerated cell senescence. Canyons mostly form between LADs and are enriched in genes and enhancers. H3K27me3 loss is correlated with up-regulation of key senescence genes, indicating a link between global chromatin changes and local gene expression regulation. Lamin B1 reduction in proliferating cells triggers senescence and formation of mesas and canyons. Our data illustrate profound chromatin reorganization during senescence and suggest that lamin B1 down-regulation in senescence is a key trigger of global and local chromatin changes that impact gene expression, aging, and cancer.

Genome-wide TOP2A DNA cleavage is biased toward translocated and highly transcribed loci.

Type II topoisomerases orchestrate proper DNA topology, and they are the targets of anti-cancer drugs that cause treatment-related leukemias with balanced translocations. Here, we develop a high-throughput sequencing technology to define TOP2 cleavage sites at single-base precision, and use the technology to characterize TOP2A cleavage genome-wide in the human K562 leukemia cell line. We find that TOP2A cleavage has functionally conserved local sequence preferences, occurs in cleavage cluster regions (CCRs), and is enriched in introns and lincRNA loci. TOP2A CCRs are biased toward the distal regions of gene bodies, and TOP2 poisons cause a proximal shift in their distribution. We find high TOP2A cleavage levels in genes involved in translocations in TOP2 poison-related leukemia. In addition, we find that a large proportion of genes involved in oncogenic translocations overall contain TOP2A CCRs. The TOP2A cleavage of coding and lincRNA genes is independently associated with both length and transcript abundance. Comparisons to ENCODE data reveal distinct TOP2A CCR clusters that overlap with marks of transcription, open chromatin, and enhancers. Our findings implicate TOP2A cleavage as a broad DNA damage mechanism in oncogenic translocations as well as a functional role of TOP2A cleavage in regulating transcription elongation and gene activation.

Electronic cigarette exposure disrupts blood-brain barrier integrity and promotes neuroinflammation.

Electronic cigarette (e-cigarette) use has grown substantially since inception, particularly among adolescents and combustible tobacco users. Several cigarette smoke constituents with known neurovascular effect are present in e-cigarette liquids or formed during the vapor generation. The present study establishes inhaled models of cigarette and e-cigarette use with normalized nicotine delivery, then characterizes the impact on blood-brain barrier (BBB) function. Sequencing of microvessel RNA following exposure revealed downregulation of several genes with critical roles in BBB function. Reduced protein expression of Occludin and Glut1 is also observed at the tight junction in all groups following exposure. Pro-inflammatory changes in leukocyte-endothelial cell interaction are also noted, and mice exposed to nicotine-free e-cigarettes have impaired novel object recognition performance. On this basis, it is concluded that long term e-cigarette use may adversely impact neurovascular health. The observed effects are noted to be partly independent of nicotine content and nicotine may even serve to moderate the effects of non-nicotinic components on the blood-brain barrier.

Global analysis of ribosome-associated noncoding RNAs unveils new modes of translational regulation.

Eukaryotic transcriptomes contain a major non-protein-coding component that includes precursors of small RNAs as well as long noncoding RNA (lncRNAs). Here, we utilized the mapping of ribosome footprints on RNAs to explore translational regulation of coding and noncoding RNAs in roots of Arabidopsis thaliana shifted from replete to deficient phosphorous (Pi) nutrition. Homodirectional changes in steady-state mRNA abundance and translation were observed for all but 265 annotated protein-coding genes. Of the translationally regulated mRNAs, 30% had one or more upstream ORF (uORF) that influenced the number of ribosomes on the principal protein-coding region. Nearly one-half of the 2,382 lncRNAs detected had ribosome footprints, including 56 with significantly altered translation under Pi-limited nutrition. The prediction of translated small ORFs (sORFs) by quantitation of translation termination and peptidic analysis identified lncRNAs that produce peptides, including several deeply evolutionarily conserved and significantly Pi-regulated lncRNAs. Furthermore, we discovered that natural antisense transcripts (NATs) frequently have actively translated sORFs, including five with low-Pi up-regulation that correlated with enhanced translation of the sense protein-coding mRNA. The data also confirmed translation of miRNA target mimics and lncRNAs that produce trans-acting or phased small-interfering RNA (tasiRNA/phasiRNAs). Mutational analyses of the positionally conserved sORF of TAS3a linked its translation with tasiRNA biogenesis. Altogether, this systematic analysis of ribosome-associated mRNAs and lncRNAs demonstrates that nutrient availability and translational regulation controls protein and small peptide-encoding mRNAs as well as a diverse cadre of regulatory RNAs.