Alphaherpesvirus-mediated remodeling of the cellular transcriptome results in depletion of m6A-containing transcripts.

The mechanisms by which viruses regulate host mRNAs during infection are still poorly understood. Several host transcripts that encode proteins that contribute to the anti-viral response contain the N6-methyladenosine nucleotide (m6A). In this study, we investigated if and how viruses from different (sub) families specifically affect m6A-containing host transcripts. Systematic analysis of host transcriptomes after infection with diverse types of viruses showed that m6A-methylated transcripts are selectively downregulated during infection with Sendai virus, African swine fever virus and the alphaherpesviruses herpes simplex virus 1 (HSV-1) and pseudorabies virus (PRV). Focusing on PRV and HSV-1, we found that downregulation of m6A-methylated transcripts depends on the YTHDF family of m6A-binding proteins, and correlates with localization of these proteins to enlarged P-bodies. Knockdown of YTHDF proteins in primary cells reduced PRV protein expression and increased expression of antiviral interferon-stimulated genes, suggesting that virus-induced depletion of host m6A-containing transcripts constitutes an immune evasion strategy.

Alphaherpesvirus US3 protein-mediated inhibition of the m6A mRNA methyltransferase complex.

Chemical modifications of mRNA, the so-called epitranscriptome, represent an additional layer of post-transcriptional regulation of gene expression. The most common epitranscriptomic modification, N6-methyladenosine (m6A), is generated by a multi-subunit methyltransferase complex. We show that alphaherpesvirus kinases trigger phosphorylation of several components of the m6A methyltransferase complex, including METTL3, METTL14, and WTAP, which correlates with inhibition of the complex and a near complete loss of m6A levels in mRNA of virus-infected cells. Expression of the viral US3 protein is necessary and sufficient for phosphorylation and inhibition of the m6A methyltransferase complex. Although m6A methyltransferase complex inactivation is not essential for virus replication in cell culture, the consensus m6A methylation motif is under-represented in alphaherpesvirus genomes, suggesting evolutionary pressure against methylation of viral transcripts. Together, these findings reveal that phosphorylation can be associated with inactivation of the m6A methyltransferase complex, in this case mediated by the viral US3 protein.

Repurposing an adenine riboswitch into a fluorogenic imaging and sensing tag.

Fluorogenic RNA aptamers are used to genetically encode fluorescent RNA and to construct RNA-based metabolite sensors. Unlike naturally occurring aptamers that efficiently fold and undergo metabolite-induced conformational changes, fluorogenic aptamers can exhibit poor folding, which limits their cellular fluorescence. To overcome this, we evolved a naturally occurring well-folded adenine riboswitch into a fluorogenic aptamer. We generated a library of roughly 1015 adenine aptamer-like RNAs in which the adenine-binding pocket was randomized for both size and sequence, and selected Squash, which binds and activates the fluorescence of green fluorescent protein-like fluorophores. Squash exhibits markedly improved in-cell folding and highly efficient metabolite-dependent folding when fused to a S-adenosylmethionine (SAM)-binding aptamer. A Squash-based ratiometric sensor achieved quantitative SAM measurements, revealed cell-to-cell heterogeneity in SAM levels and revealed metabolic origins of SAM. These studies show that the efficient folding of naturally occurring aptamers can be exploited to engineer well-folded cell-compatible fluorogenic aptamers and devices.

Antibody cross-reactivity accounts for widespread appearance of m1A in 5'UTRs.

N1-methyladenosine (m1A) was proposed to be a highly prevalent modification in mRNA 5'UTRs based on mapping studies using an m1A-binding antibody. We developed a bioinformatic approach to discover m1A and other modifications in mRNA throughout the transcriptome by analyzing preexisting ultra-deep RNA-Seq data for modification-induced misincorporations. Using this approach, we detected appreciable levels of m1A only in one mRNA: the mitochondrial MT-ND5 transcript. As an alternative approach, we also developed an antibody-based m1A-mapping approach to detect m1A at single-nucleotide resolution, and confirmed that the commonly used m1A antibody maps sites to the transcription-start site in mRNA 5'UTRs. However, further analysis revealed that these were false-positives caused by binding of the antibody to the m7G-cap. A different m1A antibody that lacks cap-binding cross-reactivity does not show enriched binding in 5'UTRs. These results demonstrate that high-stoichiometry m1A sites are exceedingly rare in mRNAs and that previous mappings of m1A to 5'UTRs were the result of antibody cross-reactivity to the 5' cap.

m6A enhances the phase separation potential of mRNA.

N6-methyladenosine (m6A) is the most prevalent modified nucleotide in mRNA1,2, with around 25% of mRNAs containing at least one m6A. Methylation of mRNA to form m6A is required for diverse cellular and physiological processes3. Although the presence of m6A in an mRNA can affect its fate in different ways, it is unclear how m6A directs this process and why the effects of m6A can vary in different cellular contexts. Here we show that the cytosolic m6A-binding proteins-YTHDF1, YTHDF2 and YTHDF3-undergo liquid-liquid phase separation in vitro and in cells. This phase separation is markedly enhanced by mRNAs that contain multiple, but not single, m6A residues. Polymethylated mRNAs act as a multivalent scaffold for the binding of YTHDF proteins, juxtaposing their low-complexity domains and thereby leading to phase separation. The resulting mRNA-YTHDF complexes then partition into different endogenous phase-separated compartments, such as P-bodies, stress granules or neuronal RNA granules. m6A-mRNA is subject to compartment-specific regulation, including a reduction in the stability and translation of mRNA. These studies reveal that the number and distribution of m6A sites in cellular mRNAs can regulate and influence the composition of the phase-separated transcriptome, and suggest that the cellular properties of m6A-modified mRNAs are governed by liquid-liquid phase separation principles.

HIF-1α Stabilization Increases miR-210 Eliciting First Trimester Extravillous Trophoblast Mitochondrial Dysfunction.

Preeclampsia is associated with first trimester placental dysfunction. miR-210, a small non-coding RNA, is increased in the preeclamptic placenta. The effects of elevated miR-210 on placental function remain unclear. The objectives of this study were to identify targets of miR-210 in first trimester primary extravillous trophoblasts (EVTs) and to investigate functional pathways altered by elevated placental miR-210 during early pregnancy. EVTs isolated from first trimester placentas were exposed to cobalt chloride (CoCl2), a HIF-1α stabilizer and hypoxia mimetic, and miR-210 expression by qPCR, HIF1α protein levels by western blot and cell invasion were assessed. A custom TruSeq RNA array, including all known/predicted miR-210 targets, was run using miR-210 and miR-negative control transfected EVTs. Mitochondrial function was assessed by high resolution respirometry in transfected EVTs. EVTs exposed to CoCl2 showed a dose and time-dependent increase in miR-210 and HIF1α and reductions in cell invasion. The TruSeq array identified 49 altered genes in miR-210 transfected EVTs with 27 genes repressed and 22 enhanced. Three of the top six significantly repressed genes, NDUFA4, SDHD, and ISCU, are associated with mitochondrial function. miR-210 transfected EVTs had decreased maximal, complex II and complex I+II mitochondrial respiration. This study suggests that miR-210 alters first trimester trophoblast function. miR-210 overexpression alters EVT mitochondrial function in early pregnancy. Mitochondrial dysfunction may lead to increased reactive oxygen species, trophoblast cell damage and likely contributes to the pathogenesis of preeclampsia.

Clock Regulation of Metabolites Reveals Coupling between Transcription and Metabolism.

The intricate connection between the circadian clock and metabolism remains poorly understood. We used high temporal resolution metabolite profiling to explore clock regulation of mouse liver and cell-autonomous metabolism. In liver, ∼50% of metabolites were circadian, with enrichment of nucleotide, amino acid, and methylation pathways. In U2 OS cells, 28% were circadian, including amino acids and NAD biosynthesis metabolites. Eighteen metabolites oscillated in both systems and a subset of these in primary hepatocytes. These 18 metabolites were enriched in methylation and amino acid pathways. To assess clock dependence of these rhythms, we used genetic perturbation. BMAL1 knockdown diminished metabolite rhythms, while CRY1 or CRY2 perturbation generally shortened or lengthened rhythms, respectively. Surprisingly, CRY1 knockdown induced 8 hr rhythms in amino acid, methylation, and vitamin metabolites, decoupling metabolite from transcriptional rhythms, with potential impact on nutrient sensing in vivo. These results provide the first comprehensive views of circadian liver and cell-autonomous metabolism.

Mapping m6A at Individual-Nucleotide Resolution Using Crosslinking and Immunoprecipitation (miCLIP).

N 6 -methyladenosine (m6A) is the most abundant modified base in eukaryotic mRNA and has been linked to diverse effects on mRNA fate. Current m6A mapping approaches localize m6A residues to 100-200 nt-long regions of transcripts. The precise position of m6A in mRNAs cannot be identified on a transcriptome-wide level because there are no chemical methods to distinguish between m6A and adenosine. Here, we describe a method for using anti-m6A antibodies to induce specific mutational signatures at m6A residues after ultraviolet light-induced antibody-RNA crosslinking and reverse transcription. Then, we describe how to use these mutational signatures to map m6A residues at nucleotide resolution. Taken together, our protocol allows for high-throughput detection of individual m6A residues throughout the transcriptome.

MetaPlotR: a Perl/R pipeline for plotting metagenes of nucleotide modifications and other transcriptomic sites.

SUMMARY: An increasing number of studies are mapping protein binding and nucleotide modifications sites throughout the transcriptome. Often, these sites cluster in certain regions of the transcript, giving clues to their function. Hence, it is informative to summarize where in the transcript these sites occur. A metagene is a simple and effective tool for visualizing the distribution of sites along a simplified transcript model. In this work, we introduce MetaPlotR, a Perl/R pipeline for creating metagene plots. AVAILABILITY AND IMPLEMENTATION: The code and associated tutorial are available at https://github.com/olarerin/metaPlotR . CONTACT: srj2003@med.cornell.edu.

Deciphering the Duality of Clock and Growth Metabolism in a Cell Autonomous System Using NMR Profiling of the Secretome.

Oscillations in circadian metabolism are crucial to the well being of organism. Our understanding of metabolic rhythms has been greatly enhanced by recent advances in high-throughput systems biology experimental techniques and data analysis. In an in vitro setting, metabolite rhythms can be measured by time-dependent sampling over an experimental period spanning one or more days at sufficent resolution to elucidate rhythms. We hypothesized that cellular metabolic effects over such a time course would be influenced by both oscillatory and circadian-independent cell metabolic effects. Here we use nuclear magnetic resonance (NMR) spectroscopy-based metabolic profiling of mammalian cell culture media of synchronized U2 OS cells containing an intact transcriptional clock. The experiment was conducted over 48 h, typical for circadian biology studies, and samples collected at 2 h resolution to unravel such non-oscillatory effects. Our data suggest specific metabolic activities exist that change continuously over time in this settting and we demonstrate that the non-oscillatory effects are generally monotonic and possible to model with multivariate regression. Deconvolution of such non-circadian persistent changes are of paramount importance to consider while studying circadian metabolic oscillations.

Single-nucleotide-resolution mapping of m6A and m6Am throughout the transcriptome.

N(6)-methyladenosine (m6A) is the most abundant modified base in eukaryotic mRNA and has been linked to diverse effects on mRNA fate. Current mapping approaches localize m6A residues to transcript regions 100-200 nt long but cannot identify precise m6A positions on a transcriptome-wide level. Here we developed m6A individual-nucleotide-resolution cross-linking and immunoprecipitation (miCLIP) and used it to demonstrate that antibodies to m6A can induce specific mutational signatures at m6A residues after ultraviolet light-induced antibody-RNA cross-linking and reverse transcription. We found that these antibodies similarly induced mutational signatures at N(6),2'-O-dimethyladenosine (m6Am), a modification found at the first nucleotide of certain mRNAs. Using these signatures, we mapped m6A and m6Am at single-nucleotide resolution in human and mouse mRNA and identified small nucleolar RNAs (snoRNAs) as a new class of m6A-containing non-coding RNAs (ncRNAs).

Placental expression of miR-517a/b and miR-517c contributes to trophoblast dysfunction and preeclampsia.

Preeclampsia is a pregnancy specific hypertensive disease that confers significant maternal and fetal risks. While the exact pathophysiology of preeclampsia is unknown, it is widely accepted that placental dysfunction is mechanistically involved. Recent studies reported aberrant expression of placenta-specific microRNAs (miRNAs) in preeclampsia including miR-517a/b and miR-517c. Using placental biopsies from a preeclampsia case-control study, we found increased expression of miR-517a/b in term and preterm preeclampsia vs controls, while, miR-517c was increased only in preterm preeclampsia vs controls. To determine if miR-517a/b and miR-517c are regulated by hypoxia, we treated first trimester primary extravillous trophoblast cells (EVTs) with a hypoxia mimetic and found both were induced. To test for a mechanistic role in placental function, we overexpressed miR-517a/b or miR-517c in EVTs which resulted in decreased trophoblast invasion. Additionally, we found that miR-517a/b and miR-517c overexpression increased expression of the anti-angiogenic protein, sFLT1. The regulation of sFLT1 is mostly unknown, however, TNFSF15, a cytokine involved in FLT1 splicing, was also increased by miR-517a/b and miR-517c in EVTs. In summary, we demonstrate that miR-517a/b and miR-517c contribute to the development of preeclampsia and suggest that these miRNAs play a critical role in regulating trophoblast and placental function.

Assessing the prevalence of mycoplasma contamination in cell culture via a survey of NCBI's RNA-seq archive.

Mycoplasmas are notorious contaminants of cell culture and can have profound effects on host cell biology by depriving cells of nutrients and inducing global changes in gene expression. Over the last two decades, sentinel testing has revealed wide-ranging contamination rates in mammalian culture. To obtain an unbiased assessment from hundreds of labs, we analyzed sequence data from 9395 rodent and primate samples from 884 series in the NCBI Sequence Read Archive. We found 11% of these series were contaminated (defined as ≥100 reads/million mapping to mycoplasma in one or more samples). Ninety percent of mycoplasma-mapped reads aligned to ribosomal RNA. This was unexpected given 37% of contaminated series used poly(A)-selection for mRNA enrichment. Lastly, we examined the relationship between mycoplasma contamination and host gene expression in a single cell RNA-seq dataset and found 61 host genes (P < 0.001) were significantly associated with mycoplasma-mapped read counts. In all, this study suggests mycoplasma contamination is still prevalent today and poses substantial risk to research quality.

Extensive variation in chromatin states across humans.

The majority of disease-associated variants lie outside protein-coding regions, suggesting a link between variation in regulatory regions and disease predisposition. We studied differences in chromatin states using five histone modifications, cohesin, and CTCF in lymphoblastoid lines from 19 individuals of diverse ancestry. We found extensive signal variation in regulatory regions, which often switch between active and repressed states across individuals. Enhancer activity is particularly diverse among individuals, whereas gene expression remains relatively stable. Chromatin variability shows genetic inheritance in trios, correlates with genetic variation and population divergence, and is associated with disruptions of transcription factor binding motifs. Overall, our results provide insights into chromatin variation among humans.

miR-210 inhibits trophoblast invasion and is a serum biomarker for preeclampsia.

Preeclampsia is characterized by hypertension and proteinuria in pregnant women. Its exact cause is unknown. Preeclampsia increases the risk of maternal and fetal morbidity and mortality. Although delivery, often premature, is the only known cure, early targeted interventions may improve maternal and fetal outcomes. Successful intervention requires a better understanding of the molecular etiology of preeclampsia and the development of accurate methods to predict women at risk. To this end, we tested the role of miR-210, a miRNA up-regulated in preeclamptic placentas, in first-trimester extravillous trophoblasts. miR-210 overexpression reduced trophoblast invasion, a process necessary for uteroplacental perfusion, in an extracellular signal-regulated kinase/mitogen-activated protein kinase-dependent manner. Conversely, miR-210 inhibition promoted invasion. Furthermore, given that the placenta secretes miRNAs into the maternal circulation, we tested if serum expression of miR-210 was associated with the disease. We measured miR-210 expression in two clinical studies: a case-control study and a prospective cohort study. Serum miR-210 expression was significantly associated with a diagnosis of preeclampsia (P = 0.007, area under the receiver operator curves = 0.81) and was predictive of the disease, even months before clinical diagnosis (P < 0.0001, area under the receiver operator curve = 0.89). Hence, we conclude that aberrant expression of miR-210 may contribute to trophoblast function and that miR-210 is a novel predictive serum biomarker for preeclampsia that can help in identifying at-risk women for monitoring and treatment.

A functional genomics screen for microRNA regulators of NF-kappaB signaling.

BACKGROUND: The nuclear factor-KappaB (NF-κB) pathway is conserved from fruit flies to humans and is a key mediator of inflammatory signaling. Aberrant regulation of NF-κB is associated with several disorders including autoimmune disease, chronic inflammation, and cancer, making the NF-κB pathway an attractive therapeutic target. Many regulatory components of the NF-κB pathway have been identified, including microRNAs (miRNAs). miRNAs are small non-coding RNAs and are common components of signal transduction pathways. Here we present a cell-based functional genomics screen to systematically identify miRNAs that regulate NF-κB signaling. RESULTS: We screened a library of miRNA mimics using a NF-κB reporter cell line in the presence and absence of tumor necrosis factor (+/- TNF). There were 9 and 15 hits in the -TNF and +TNF screens, respectively. We identified putative functional targets of these hits by integrating computational predictions with NF-κB modulators identified in a previous genome-wide cDNA screen. miR-517a and miR-517c were the top hits, activating the reporter 86- and 126-fold, respectively. Consistent with these results, miR-517a/c induced the expression of endogenous NF-κB targets and promoted the nuclear localization of p65 and the degradation of IκB. We identified TNFAIP3 interacting protein1 (TNIP1) as a target and characterized a functional SNP in the miR-517a/c binding site. Lastly, miR-517a/c induced apoptosis in vitro, which was phenocopied by knockdown of TNIP1. CONCLUSIONS: Our study suggests that miRNAs are common components of NF-κB signaling and miR-517a/c may play an important role in linking NF-κB signaling with cell survival through TNIP1.

Adhesion regulates MAP kinase/ternary complex factor exchange to control a proliferative transcriptional switch.

BACKGROUND: The ternary complex factors (TCFs; Elk1, Net, and Sap-1) are growth factor-responsive transcription cofactors of serum response factor (SRF) and are activated by MAP kinase (MAPK) phosphorylation to regulate immediate early gene transcription. Although cell adhesion also can regulate immediate early genes and proliferation, the mechanism for this effect has remained unexplored. RESULTS: Restricting adhesion and spreading of G(0)-synchronized cells on substrates with decreasing size of micropatterned islands of fibronectin suppressed serum-induced immediate early gene expression and S phase entry. Knockdown of Sap-1 decreased expression of the immediate early genes egr1 and fos and subsequent proliferation normally present with high adhesion, whereas knockdown of Net rescued egr1 and fos expression and proliferation normally suppressed by low adhesion. Chromatin immunoprecipitation studies showed increased occupancy of egr1 and fos promoters by Sap-1 with high adhesion, whereas low adhesion increased Net occupancy. This switch in TCF promoter binding was regulated by an adhesion-mediated switch in MAPK activity. Increasing adhesion enhanced serum-induced JNK activity while suppressing p38 activity, leading to increased Sap-1 phosphorylation and Net dephosphorylation, and switching Net with Sap-1 at egr1 and fos promoters to support proliferation. Microarray studies confirmed this switch in TCF regulation of proliferative genes and uncovered novel gene targets and functions coregulated by Sap-1 and Net. CONCLUSIONS: These data demonstrate a key role for the TCFs in adhesion-induced transcription and proliferation and reveal a novel MAPK/TCF transcriptional switch that controls this process.