Acetylation of Cytidine in mRNA Promotes Translation Efficiency.

Generation of the "epitranscriptome" through post-transcriptional ribonucleoside modification embeds a layer of regulatory complexity into RNA structure and function. Here, we describe N4-acetylcytidine (ac4C) as an mRNA modification that is catalyzed by the acetyltransferase NAT10. Transcriptome-wide mapping of ac4C revealed discretely acetylated regions that were enriched within coding sequences. Ablation of NAT10 reduced ac4C detection at the mapped mRNA sites and was globally associated with target mRNA downregulation. Analysis of mRNA half-lives revealed a NAT10-dependent increase in stability in the cohort of acetylated mRNAs. mRNA acetylation was further demonstrated to enhance substrate translation in vitro and in vivo. Codon content analysis within ac4C peaks uncovered a biased representation of cytidine within wobble sites that was empirically determined to influence mRNA decoding efficiency. These findings expand the repertoire of mRNA modifications to include an acetylated residue and establish a role for ac4C in the regulation of mRNA translation.

Detection of ac4C in human mRNA is preserved upon data reassessment.

We recently reported the distribution of N4-acetylcytidine (ac4C) in HeLa mRNA at base resolution through chemical reduction and the induction of C:T mismatches in sequencing (RedaC:T-seq). Our results contradicted an earlier report from Schwartz and colleagues utilizing a similar method termed ac4C-seq. Here, we revisit both datasets and reaffirm our findings. Through RedaC:T-seq reanalysis, we establish a low basal error rate at unmodified nucleotides that is not skewed to any specific mismatch type and a prominent increase in C:T substitutions as the dominant mismatch type in both treated wild-type replicates, with a high degree of reproducibility across replicates. In contrast, through ac4C-seq reanalysis, we uncover significant data quality issues including insufficient depth, with one wild-type replicate yielding 2.7 million reads, inconsistencies in reduction efficiencies between replicates, and an overall increase in mismatches involving thymine that could obscure ac4C detection. These analyses bolster the detection of ac4C in HeLa mRNA through RedaC:T-seq.

Direct epitranscriptomic regulation of mammalian translation initiation through N4-acetylcytidine.

mRNA function is influenced by modifications that modulate canonical nucleobase behavior. We show that a single modification mediates distinct impacts on mRNA translation in a position-dependent manner. Although cytidine acetylation (ac4C) within protein-coding sequences stimulates translation, ac4C within 5' UTRs impacts protein synthesis at the level of initiation. 5' UTR acetylation promotes initiation at upstream sequences, competitively inhibiting annotated start codons. Acetylation further directly impedes initiation at optimal AUG contexts: ac4C within AUG-flanking Kozak sequences reduced initiation in base-resolved transcriptome-wide HeLa results and in vitro utilizing substrates with site-specific ac4C incorporation. Cryo-EM of mammalian 80S initiation complexes revealed that ac4C in the -1 position adjacent to an AUG start codon disrupts an interaction between C and hypermodified t6A at nucleotide 37 of the initiator tRNA. These findings demonstrate the impact of RNA modifications on nucleobase function at a molecular level and introduce mRNA acetylation as a factor regulating translation in a location-specific manner.

Epigenetic Regulation of DNA Repair Pathway Choice by MacroH2A1 Splice Variants Ensures Genome Stability.

The inactive X chromosome (Xi) is inherently susceptible to genomic aberrations. Replication stress (RS) has been proposed as an underlying cause, but the mechanisms that protect from Xi instability remain unknown. Here, we show that macroH2A1.2, an RS-protective histone variant enriched on the Xi, is required for Xi integrity and female survival. Mechanistically, macroH2A1.2 counteracts its structurally distinct and equally Xi-enriched alternative splice variant, macroH2A1.1. Comparative proteomics identified a role for macroH2A1.1 in alternative end joining (alt-EJ), which accounts for Xi anaphase defects in the absence of macroH2A1.2. Genomic instability was rescued by simultaneous depletion of macroH2A1.1 or alt-EJ factors, and mice deficient for both macroH2A1 variants harbor no overt female defects. Notably, macroH2A1 splice variant imbalance affected alt-EJ capacity also in tumor cells. Together, these findings identify macroH2A1 splicing as a modulator of genome maintenance that ensures Xi integrity and may, more broadly, predict DNA repair outcome in malignant cells.

Replication Stress Shapes a Protective Chromatin Environment across Fragile Genomic Regions.

Recent integrative epigenome analyses highlight the importance of functionally distinct chromatin states for accurate cell function. How these states are established and maintained is a matter of intense investigation. Here, we present evidence for DNA damage as an unexpected means to shape a protective chromatin environment at regions of recurrent replication stress (RS). Upon aberrant fork stalling, DNA damage signaling and concomitant H2AX phosphorylation coordinate the FACT-dependent deposition of macroH2A1.2, a histone variant that promotes DNA repair by homologous recombination (HR). MacroH2A1.2, in turn, facilitates the accumulation of the tumor suppressor and HR effector BRCA1 at replication forks to protect from RS-induced DNA damage. Consequently, replicating primary cells steadily accrue macroH2A1.2 at fragile regions, whereas macroH2A1.2 loss in these cells triggers DNA damage signaling-dependent senescence, a hallmark of RS. Altogether, our findings demonstrate that recurrent DNA damage contributes to the chromatin landscape to ensure the epigenomic integrity of dividing cells.

PancrESS - a meta-analysis resource for understanding cell-type specific expression in the human pancreas.

BACKGROUND: The human pancreas is composed of specialized cell types producing hormones and enzymes critical to human health. These specialized functions are the result of cell type-specific transcriptional programs which manifest in cell-specific gene expression. Understanding these programs is essential to developing therapies for pancreatic disorders. Transcription in the human pancreas has been widely studied by single-cell RNA technologies, however the diversity of protocols and analysis methods hinders their interpretability in the aggregate. RESULTS: In this work, we perform a meta-analysis of pancreatic single-cell RNA sequencing data. We present a database for reference transcriptome abundances and cell-type specificity metrics. This database facilitates the identification and definition of marker genes within the pancreas. Additionally, we introduce a versatile tool which is freely available as an R package, and should permit integration into existing workflows. Our tool accepts count data files generated by widely-used single-cell gene expression platforms in their original format, eliminating an additional pre-formatting step. Although we designed it to calculate expression specificity of pancreas cell types, our tool is agnostic to the biological source of count data, extending its applicability to other biological systems. CONCLUSIONS: Our findings enhance the current understanding of expression specificity within the pancreas, surpassing previous work in terms of scope and detail. Furthermore, our database and tool enable researchers to perform similar calculations in diverse biological systems, expanding the applicability of marker gene identification and facilitating comparative analyses.

Activation of Crtc2/Creb1 in skeletal muscle enhances weight loss during intermittent fasting.

The Creb-Regulated Transcriptional Coactivator (Crtc) family of transcriptional coregulators drive Creb1-mediated transcription effects on metabolism in many tissues, but the in vivo effects of Crtc2/Creb1 transcription on skeletal muscle metabolism are not known. Skeletal muscle-specific overexpression of Crtc2 (Crtc2 mice) induced greater mitochondrial activity, metabolic flux capacity for both carbohydrates and fats, improved glucose tolerance and insulin sensitivity, and increased oxidative capacity, supported by upregulation of key metabolic genes. Crtc2 overexpression led to greater weight loss during alternate day fasting (ADF), selective loss of fat rather than lean mass, maintenance of higher energy expenditure during the fast and reduced binge-eating during the feeding period. ADF downregulated most of the mitochondrial electron transport genes, and other regulators of mitochondrial function, that were substantially reversed by Crtc2-driven transcription. Glucocorticoids acted with AMPK to drive atrophy and mitophagy, which was reversed by Crtc2/Creb1 signaling. Crtc2/Creb1-mediated signaling coordinates metabolic adaptations in skeletal muscle that explain how Crtc2/Creb1 regulates metabolism and weight loss.

TET-catalyzed oxidation of intragenic 5-methylcytosine regulates CTCF-dependent alternative splicing.

Intragenic 5-methylcytosine and CTCF mediate opposing effects on pre-mRNA splicing: CTCF promotes inclusion of weak upstream exons through RNA polymerase II pausing, whereas 5-methylcytosine evicts CTCF, leading to exon exclusion. However, the mechanisms governing dynamic DNA methylation at CTCF-binding sites were unclear. Here, we reveal the methylcytosine dioxygenases TET1 and TET2 as active regulators of CTCF-mediated alternative splicing through conversion of 5-methylcytosine to its oxidation derivatives. 5-hydroxymethylcytosine and 5-carboxylcytosine are enriched at an intragenic CTCF-binding sites in the CD45 model gene and are associated with alternative exon inclusion. Reduced TET levels culminate in increased 5-methylcytosine, resulting in CTCF eviction and exon exclusion. In vitro analyses establish the oxidation derivatives are not sufficient to stimulate splicing, but efficiently promote CTCF association. We further show genomewide that reciprocal exchange of 5-hydroxymethylcytosine and 5-methylcytosine at downstream CTCF-binding sites is a general feature of alternative splicing in naïve and activated CD4(+) T cells. These findings significantly expand our current concept of the pre-mRNA "splicing code" to include dynamic intragenic DNA methylation catalyzed by the TET proteins.

Diversity and dynamics of the Drosophila transcriptome.

Animal transcriptomes are dynamic, with each cell type, tissue and organ system expressing an ensemble of transcript isoforms that give rise to substantial diversity. Here we have identified new genes, transcripts and proteins using poly(A)+ RNA sequencing from Drosophila melanogaster in cultured cell lines, dissected organ systems and under environmental perturbations. We found that a small set of mostly neural-specific genes has the potential to encode thousands of transcripts each through extensive alternative promoter usage and RNA splicing. The magnitudes of splicing changes are larger between tissues than between developmental stages, and most sex-specific splicing is gonad-specific. Gonads express hundreds of previously unknown coding and long non-coding RNAs (lncRNAs), some of which are antisense to protein-coding genes and produce short regulatory RNAs. Furthermore, previously identified pervasive intergenic transcription occurs primarily within newly identified introns. The fly transcriptome is substantially more complex than previously recognized, with this complexity arising from combinatorial usage of promoters, splice sites and polyadenylation sites.

Independence between pre-mRNA splicing and DNA methylation in an isogenic minigene resource.

Actively transcribed genes adopt a unique chromatin environment with characteristic patterns of enrichment. Within gene bodies, H3K36me3 and cytosine DNA methylation are elevated at exons of spliced genes and have been implicated in the regulation of pre-mRNA splicing. H3K36me3 is further responsive to splicing, wherein splicing inhibition led to a redistribution and general reduction over gene bodies. In contrast, little is known of the mechanisms supporting elevated DNA methylation at actively spliced genic locations. Recent evidence associating the de novo DNA methyltransferase Dnmt3b with H3K36me3-rich chromatin raises the possibility that genic DNA methylation is influenced by splicing-associated H3K36me3. Here, we report the generation of an isogenic resource to test the direct impact of splicing on chromatin. A panel of minigenes of varying splicing potential were integrated into a single FRT site for inducible expression. Profiling of H3K36me3 confirmed the established relationship to splicing, wherein levels were directly correlated with splicing efficiency. In contrast, DNA methylation was equivalently detected across the minigene panel, irrespective of splicing and H3K36me3 status. In addition to revealing a degree of independence between genic H3K36me3 and DNA methylation, these findings highlight the generated minigene panel as a flexible platform for the query of splicing-dependent chromatin modifications.

Cajal body function in genome organization and transcriptome diversity.

Nuclear bodies contribute to non-random organization of the human genome and nuclear function. Using a major prototypical nuclear body, the Cajal body, as an example, we suggest that these structures assemble at specific gene loci located across the genome as a result of high transcriptional activity. Subsequently, target genes are physically clustered in close proximity in Cajal body-containing cells. However, Cajal bodies are observed in only a limited number of human cell types, including neuronal and cancer cells. Ultimately, Cajal body depletion perturbs splicing kinetics by reducing target small nuclear RNA (snRNA) transcription and limiting the levels of spliceosomal snRNPs, including their modification and turnover following each round of RNA splicing. As such, Cajal bodies are capable of shaping the chromatin interaction landscape and the transcriptome by influencing spliceosome kinetics. Future studies should concentrate on characterizing the direct influence of Cajal bodies upon snRNA gene transcriptional dynamics. Also see the video abstract here.

Controlled DNA double-strand break induction in mice reveals post-damage transcriptome stability.

DNA double-strand breaks (DSBs) and their repair can cause extensive epigenetic changes. As a result, DSBs have been proposed to promote transcriptional and, ultimately, physiological dysfunction via both cell-intrinsic and cell-non-autonomous pathways. Studying the consequences of DSBs in higher organisms has, however, been hindered by a scarcity of tools for controlled DSB induction. Here, we describe a mouse model that allows for both tissue-specific and temporally controlled DSB formation at ∼140 defined genomic loci. Using this model, we show that DSBs promote a DNA damage signaling-dependent decrease in gene expression in primary cells specifically at break-bearing genes, which is reversed upon DSB repair. Importantly, we demonstrate that restoration of gene expression can occur independently of cell cycle progression, underlining its relevance for normal tissue maintenance. Consistent with this, we observe no evidence for persistent transcriptional repression in response to a multi-day course of continuous DSB formation and repair in mouse lymphocytes in vivo Together, our findings reveal an unexpected capacity of primary cells to maintain transcriptome integrity in response to DSBs, pointing to a limited role for DNA damage as a mediator of cell-autonomous epigenetic dysfunction.

Comparative validation of the D. melanogaster modENCODE transcriptome annotation.

Accurate gene model annotation of reference genomes is critical for making them useful. The modENCODE project has improved the D. melanogaster genome annotation by using deep and diverse high-throughput data. Since transcriptional activity that has been evolutionarily conserved is likely to have an advantageous function, we have performed large-scale interspecific comparisons to increase confidence in predicted annotations. To support comparative genomics, we filled in divergence gaps in the Drosophila phylogeny by generating draft genomes for eight new species. For comparative transcriptome analysis, we generated mRNA expression profiles on 81 samples from multiple tissues and developmental stages of 15 Drosophila species, and we performed cap analysis of gene expression in D. melanogaster and D. pseudoobscura. We also describe conservation of four distinct core promoter structures composed of combinations of elements at three positions. Overall, each type of genomic feature shows a characteristic divergence rate relative to neutral models, highlighting the value of multispecies alignment in annotating a target genome that should prove useful in the annotation of other high priority genomes, especially human and other mammalian genomes that are rich in noncoding sequences. We report that the vast majority of elements in the annotation are evolutionarily conserved, indicating that the annotation will be an important springboard for functional genetic testing by the Drosophila community.

The developmental transcriptome of Drosophila melanogaster.

Drosophila melanogaster is one of the most well studied genetic model organisms; nonetheless, its genome still contains unannotated coding and non-coding genes, transcripts, exons and RNA editing sites. Full discovery and annotation are pre-requisites for understanding how the regulation of transcription, splicing and RNA editing directs the development of this complex organism. Here we used RNA-Seq, tiling microarrays and cDNA sequencing to explore the transcriptome in 30 distinct developmental stages. We identified 111,195 new elements, including thousands of genes, coding and non-coding transcripts, exons, splicing and editing events, and inferred protein isoforms that previously eluded discovery using established experimental, prediction and conservation-based approaches. These data substantially expand the number of known transcribed elements in the Drosophila genome and provide a high-resolution view of transcriptome dynamics throughout development.

Design of RNA splicing analysis null models for post hoc filtering of Drosophila head RNA-Seq data with the splicing analysis kit (Spanki).

BACKGROUND: The production of multiple transcript isoforms from one gene is a major source of transcriptome complexity. RNA-Seq experiments, in which transcripts are converted to cDNA and sequenced, allow the resolution and quantification of alternative transcript isoforms. However, methods to analyze splicing are underdeveloped and errors resulting in incorrect splicing calls occur in every experiment. RESULTS: We used RNA-Seq data to develop sequencing and aligner error models. By applying these error models to known input from simulations, we found that errors result from false alignment to minor splice motifs and antisense stands, shifted junction positions, paralog joining, and repeat induced gaps. By using a series of quantitative and qualitative filters, we eliminated diagnosed errors in the simulation, and applied this to RNA-Seq data from Drosophila melanogaster heads. We used high-confidence junction detections to specifically interrogate local splicing differences between transcripts. This method out-performed commonly used RNA-seq methods to identify known alternative splicing events in the Drosophila sex determination pathway. We describe a flexible software package to perform these tasks called Splicing Analysis Kit (Spanki), available at http://www.cbcb.umd.edu/software/spanki. CONCLUSIONS: Splice-junction centric analysis of RNA-Seq data provides advantages in specificity for detection of alternative splicing. Our software provides tools to better understand error profiles in RNA-Seq data and improve inference from this new technology. The splice-junction centric approach that this software enables will provide more accurate estimates of differentially regulated splicing than current tools.

Demasculinization of X chromosomes in the Drosophila genus.

X chromosomes evolve differently from autosomes, but general governing principles have not emerged. For example, genes with male-biased expression are under-represented on the X chromosome of D. melanogaster, but are randomly distributed in the genome of Anopheles gambiae. In direct global profiling experiments using species-specific microarrays, we find a nearly identical paucity of genes with male-biased expression on D. melanogaster, D. simulans, D. yakuba, D. ananassae, D. virilis and D. mojavensis X chromosomes. We observe the same under-representation on the neo-X of D. pseudoobscura. It has been suggested that precocious meiotic silencing of the X chromosome accounts for reduced X chromosome male-biased expression in nematodes, mammals and Drosophila. We show that X chromosome genes with male-biased expression are under-represented in somatic cells and in mitotic male germ cells. These data are incompatible with simple X chromosome inactivation models. Using expression profiling and comparative sequence analysis, we show that selective gene extinction on the X chromosome, creation of new genes on autosomes and changed genomic location of existing genes contribute to the unusual X chromosome gene content.

Constraint and turnover in sex-biased gene expression in the genus Drosophila.

Both genome content and deployment contribute to phenotypic differences between species. Sex is the most important difference between individuals in a species and has long been posited to be rapidly evolving. Indeed, in the Drosophila genus, traits such as sperm length, genitalia, and gonad size are the most obvious differences between species. Comparative analysis of sex-biased expression should deepen our understanding of the relationship between genome content and deployment during evolution. Using existing and newly assembled genomes, we designed species-specific microarrays to examine sex-biased expression of orthologues and species-restricted genes in D. melanogaster, D. simulans, D. yakuba, D. ananassae, D. pseudoobscura, D. virilis and D. mojavensis. We show that averaged sex-biased expression changes accumulate monotonically over time within the genus. However, different genes contribute to expression variance within species groups compared to between groups. We observed greater turnover of species-restricted genes with male-biased expression, indicating that gene formation and extinction may play a significant part in species differences. Genes with male-biased expression also show the greatest expression and DNA sequence divergence. This higher divergence and turnover of genes with male-biased expression may be due to high transcription rates in the male germline, greater functional pleiotropy of genes expressed in females, and/or sexual competition.

Protocol for base resolution mapping of ac4C using RedaC:T-seq.

N4-acetylcytidine (ac4C) is an mRNA modification catalyzed by the enzyme N-acetyltransferase 10 (NAT10), with position-dependent effects on mRNA translation. This protocol details a procedure to map ac4C at base resolution using NaBH4-induced reduction of ac4C and conversion to thymidine followed by sequencing (RedaC:T-seq). Total RNA is ribodepleted and then treated with NaBH4 to reduce ac4C to tetrahydro-ac4C, which specifically alters base pairing during cDNA synthesis, allowing the detection of ac4C at positions called as thymidine following Illumina sequencing. For complete details on the use and execution of this protocol, please refer to Arango et al. (2022).1.

Oncogenic lncRNAs alter epigenetic memory at a fragile chromosomal site in human cancer cells.

Chromosome instability is a critical event in cancer progression. Histone H3 variant CENP-A plays a fundamental role in defining centromere identity, structure, and function but is innately overexpressed in several types of solid cancers. In the cancer background, excess CENP-A is deposited ectopically on chromosome arms, including 8q24/cMYC locus, by invading transcription-coupled H3.3 chaperone pathways. Up-regulation of lncRNAs in many cancers correlates with poor prognosis and recurrence in patients. We report that transcription of 8q24-derived oncogenic lncRNAs plays an unanticipated role in altering the 8q24 chromatin landscape by H3.3 chaperone-mediated deposition of CENP-A-associated complexes. Furthermore, a transgene cassette carrying specific 8q24-derived lncRNA integrated into a naïve chromosome locus recruits CENP-A to the new location in a cis-acting manner. These data provide a plausible mechanistic link between locus-specific oncogenic lncRNAs, aberrant local chromatin structure, and the generation of new epigenetic memory at a fragile site in human cancer cells.