Methylation at the N6 position of adenosine (m6A) is a highly prevalent and reversible modification within eukaryotic mRNAs that has been linked to many stages of RNA processing and fate. Recent studies suggest that m6A deposition and proteins involved in the m6A pathway play a diverse set of roles in either restricting or modulating the lifecycles of select viruses. Here, we report that m6A levels are significantly increased in cells infected with the oncogenic human DNA virus Kaposi's sarcoma-associated herpesvirus (KSHV). Transcriptome-wide m6A-sequencing of the KSHV-positive renal carcinoma cell line iSLK.219 during lytic reactivation revealed the presence of m6A across multiple kinetic classes of viral transcripts, and a concomitant decrease in m6A levels across much of the host transcriptome. However, we found that depletion of the m6A machinery had differential pro- and anti-viral impacts on viral gene expression depending on the cell-type analyzed. In iSLK.219 and iSLK.BAC16 cells the pathway functioned in a pro-viral manner, as depletion of the m6A writer METTL3 and the reader YTHDF2 significantly impaired virion production. In iSLK.219 cells the defect was linked to their roles in the post-transcriptional accumulation of the major viral lytic transactivator ORF50, which is m6A modified. In contrast, although the ORF50 mRNA was also m6A modified in KSHV infected B cells, ORF50 protein expression was instead increased upon depletion of METTL3, or, to a lesser extent, YTHDF2. These results highlight that the m6A pathway is centrally involved in regulating KSHV gene expression, and underscore how the outcome of this dynamically regulated modification can vary significantly between cell types.

Infection with the flavivirus Zika (ZIKV) causes neurological, immunological, and developmental defects through incompletely understood mechanisms. We report that ZIKV infection affects viral and human RNAs by altering the topology and function of N⁶-adenosine methylation (m⁶A), a modification affecting RNA structure and function. m⁶A nucleosides are abundant in ZIKV RNA, with twelve m⁶A peaks identified across full-length ZIKV RNA. m⁶A in ZIKV RNA is controlled by host methyltransferases METTL3 and METTL14 and demethylases ALKBH5 and FTO, and knockdown of methyltransferases increases, while silencing demethylases decreases, ZIKV production. YTHDF family proteins, which regulate the stability of m⁶A-modified RNA, bind to ZIKV RNA, and their silencing increases ZIKV replication. Profiling of the m⁶A methylome of host mRNAs reveals that ZIKV infection alters m⁶A location in mRNAs, methylation motifs, and target genes modified by methyltransferases. Our results identify a mechanism by which ZIKV interacts with and alters host cell functions.

Determining the spatial organization of chromatin in cells mainly relies on crosslinking-based chromosome conformation capture techniques, but resolution and signal-to-noise ratio of these approaches is limited by interference from DNA-bound proteins. Here we introduce chemical-crosslinking assisted proximity capture (CAP-C), a method that uses multifunctional chemical crosslinkers with defined sizes to capture chromatin contacts. CAP-C generates chromatin contact maps at subkilobase (sub-kb) resolution with low background noise. We applied CAP-C to formaldehyde prefixed mouse embryonic stem cells (mESCs) and investigated loop domains (median size of 200 kb) and nonloop domains (median size of 9 kb). Transcription inhibition caused a greater loss of contacts in nonloop domains than loop domains. We uncovered conserved, transcription-state-dependent chromatin compartmentalization at high resolution that is shared from Drosophila to human, and a transcription-initiation-dependent nuclear subcompartment that brings multiple nonloop domains in close proximity. We also showed that CAP-C could be used to detect native chromatin conformation without formaldehyde prefixing.

Despite explosive growth in genomic datasets, the methods for studying epigenomic mechanisms of gene regulation remain primitive. Here we present a model-based approach to systematically analyze the epigenomic functions in modulating transcription factor-DNA binding. Based on the first principles of statistical mechanics, this model considers the interactions between epigenomic modifications and a cis-regulatory module, which contains multiple binding sites arranged in any configurations. We compiled a comprehensive epigenomic dataset in mouse embryonic stem (mES) cells, including DNA methylation (MeDIP-seq and MRE-seq), DNA hydroxymethylation (5-hmC-seq), and histone modifications (ChIP-seq). We discovered correlations of transcription factors (TFs) for specific combinations of epigenomic modifications, which we term epigenomic motifs. Epigenomic motifs explained why some TFs appeared to have different DNA binding motifs derived from in vivo (ChIP-seq) and in vitro experiments. Theoretical analyses suggested that the epigenome can modulate transcriptional noise and boost the cooperativity of weak TF binding sites. ChIP-seq data suggested that epigenomic boost of binding affinities in weak TF binding sites can function in mES cells. We showed in theory that the epigenome should suppress the TF binding differences on SNP-containing binding sites in two people. Using personal data, we identified strong associations between H3K4me2/H3K9ac and the degree of personal differences in NFκかっぱB binding in SNP-containing binding sites, which may explain why some SNPs introduce much smaller personal variations on TF binding than other SNPs. In summary, this model presents a powerful approach to analyze the functions of epigenomic modifications. This model was implemented into an open source program APEG (Affinity Prediction by Epigenome and Genome, http://systemsbio.ucsd.edu/apeg).

DNA methylation at selective cytosine residues (5-methylcytosine (5mC)) and their removal by TET-mediated DNA demethylation are critical for setting up pluripotent states in early embryonic development. TET enzymes successively convert 5mC to 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC), and 5-carboxylcytosine (5caC), with 5fC and 5caC subject to removal by thymine DNA glycosylase (TDG) in conjunction with base excision repair. Early reports indicate that 5fC and 5caC could be stably detected on enhancers, promoters and gene bodies, with distinct effects on gene expression, but the mechanisms have remained elusive. Here we determined the X-ray crystal structure of yeast elongating RNA polymerase II (Pol II) in complex with a DNA template containing oxidized 5mCs, revealing specific hydrogen bonds between the 5-carboxyl group of 5caC and the conserved epi-DNA recognition loop in the polymerase. This causes a positional shift for incoming nucleoside 5'-triphosphate (NTP), thus compromising nucleotide addition. To test the implication of this structural insight in vivo, we determined the global effect of increased 5fC/5caC levels on transcription, finding that such DNA modifications indeed retarded Pol II elongation on gene bodies. These results demonstrate the functional impact of oxidized 5mCs on gene expression and suggest a novel role for Pol II as a specific and direct epigenetic sensor during transcription elongation.

In mammals, cytosine methylation (5mC) is widely distributed throughout the genome but is notably depleted from active promoters and enhancers. While the role of DNA methylation in promoter silencing has been well documented, the function of this epigenetic mark at enhancers remains unclear. Recent experiments have demonstrated that enhancers are enriched for 5-hydroxymethylcytosine (5hmC), an oxidization product of the Tet family of 5mC dioxygenases and an intermediate of DNA demethylation. These results support the involvement of Tet proteins in the regulation of dynamic DNA methylation at enhancers. By mapping DNA methylation and hydroxymethylation at base resolution, we find that deletion of Tet2 causes extensive loss of 5hmC at enhancers, accompanied by enhancer hypermethylation, reduction of enhancer activity, and delayed gene induction in the early steps of differentiation. Our results reveal that DNA demethylation modulates enhancer activity, and its disruption influences the timing of transcriptome reprogramming during cellular differentiation.

The blood-brain barrier (BBB) remains a formidable obstacle in medicine, preventing efficient penetration of chemotherapeutic and diagnostic agents to malignant gliomas. Here, a transactivator of transcription (TAT) peptide-modified gold nanoparticle platform (TAT-Au NP) with a 5 nm core size is demonstrated to be capable of crossing the BBB efficiently and delivering cargoes such as the anticancer drug doxorubicin (Dox) and Gd(3+) contrast agents to brain tumor tissues. Treatment of mice bearing intracranial glioma xenografts with pH-sensitive Dox-conjugated TAT-Au NPs via a single intravenous administration leads to significant survival benefit when compared to the free Dox. Furthermore, it is demonstrated that TAT-Au NPs are capable of delivering Gd(3+) chelates for enhanced brain tumor imaging with a prolonged retention time of Gd(3+) when compared to the free Gd(3+) chelates. Collectively, these results show promising applications of the TAT-Au NPs for enhanced malignant brain tumor therapy and non-invasive imaging.

Long noncoding RNAs (lncRNAs) represent a multidimensional class of regulatory molecules that are involved in many aspects of brain function. Emerging evidence indicates that lncRNAs are localized to the synapse; however, a direct role for their activity in this subcellular compartment in memory formation has yet to be demonstrated. Using lncRNA capture-seq, we identified a specific set of lncRNAs that accumulate in the synaptic compartment within the infralimbic prefrontal cortex of adult male C57/Bl6 mice. Among these was a splice variant related to the stress-associated lncRNA, Gas5. RNA immunoprecipitation followed by mass spectrometry and single-molecule imaging revealed that this Gas5 isoform, in association with the RNA binding proteins G3BP2 and CAPRIN1, regulates the activity-dependent trafficking and clustering of RNA granules. In addition, we found that cell-type-specific, activity-dependent, and synapse-specific knockdown of the Gas5 variant led to impaired fear extinction memory. These findings identify a new mechanism of fear extinction that involves the dynamic interaction between local lncRNA activity and RNA condensates in the synaptic compartment.

Tibetan frogs, Nanorana parkeri, are differentiated genetically but not morphologically along geographical and elevational gradients in a challenging environment, presenting a unique opportunity to investigate processes leading to speciation. Analyses of whole genomes of 63 frogs reveal population structuring and historical demography, characterized by highly restricted gene flow in a narrow geographic zone lying between matrilines West (W) and East (E). A population found only along a single tributary of the Yalu Zangbu River has the mitogenome only of E, whereas nuclear genes of W comprise 89-95% of the nuclear genome. Selection accounts for 579 broadly scattered, highly divergent regions (HDRs) of the genome, which involve 365 genes. These genes fall into 51 gene ontology (GO) functional classes, 14 of which are likely to be important in driving reproductive isolation. GO enrichment analyses of E reveal many overrepresented functional categories associated with adaptation to high elevations, including blood circulation, response to hypoxia, and UV radiation. Four genes, including DNAJC8 in the brain, TNNC1 and ADORA1 in the heart, and LAMB3 in the lung, differ in levels of expression between low- and high-elevation populations. High-altitude adaptation plays an important role in maintaining and driving continuing divergence and reproductive isolation. Use of total genomes enabled recognition of selection and adaptation in and between populations, as well as documentation of evolution along a stepped cline toward speciation.