Shaping Gene Expression Research Articles

Charting the epitranscriptomic landscape across RNA biotypes using native RNA nanopore sequencing.

RNA modifications are conserved chemical features found in all domains of life and across diverse RNA biotypes, shaping gene expression profiles and enabling rapid responses to environmental changes. Their broad chemical diversity and dynamic nature pose significant challenges for studying them comprehensively. These limitations can now be addressed through direct RNA nanopore sequencing (DRS), which allows simultaneous identification of diverse RNA modification types at single-molecule and single-nucleotide resolution. Here, we review recent efforts pioneering the use of DRS to better understand the epitranscriptomic landscape. We highlight how DRS can be applied to investigate different RNA biotypes, emphasizing the use of specialized library preparation protocols and downstream bioinformatic workflows to detect both natural and synthetic RNA modifications. Finally, we provide a perspective on the future role of DRS in epitranscriptomic research, highlighting remaining challenges and emerging opportunities from improved sequencing yields and accuracy enabled by the latest DRS chemistry.

Host and venom evolution in parasitoid wasps: does independently adapting to the same host shape the evolution of the venom gland transcriptome?

BackgroundVenoms have repeatedly evolved over 100 occasions throughout the animal tree of life, making them excellent systems for exploring convergent evolutionary novelty. Growing evidence supports that venom evolution is predominantly driven by prey or host-related selection pressures, and the expression patterns of venom glands reflect adaptive evolution. However, it remains elusive whether the evolution of expression patterns in venom glands is likewise a convergent evolution driven by their prey/host species.ResultsWe utilized parasitoid wasps that had independently adapted to Drosophila hosts as models to investigate the convergent evolution of venom gland transcriptomes in 19 hymenopteran species spanning ~ 200 million years of evolution. Comparative transcriptome analysis reveals that the global expression patterns among the venom glands of Drosophila parasitoid wasps do not achieve higher similarity compared to non-Drosophila parasitoid wasps. Further evolutionary analyses of expression patterns at the single gene, orthogroup, and Gene Ontology (GO) term levels indicate that some orthogroups/GO terms show correlation with the Drosophila parasitoid wasps. However, these groups rarely include genes highly expressed in venom glands or putative venom genes in the Drosophila parasitoid wasps.ConclusionsOur study suggests that convergent evolution may not play a predominant force shaping gene expression levels in the venom gland of the Drosophila parasitoid wasps, offering novel insights into the co-evolution between venom and prey/host.

The Dynamics of Histone Modifications during Mammalian Zygotic Genome Activation.

Mammalian fertilization initiates the reprogramming of oocytes and sperm, forming a totipotent zygote. During this intricate process, the zygotic genome undergoes a maternal-to-zygotic transition (MZT) and subsequent zygotic genome activation (ZGA), marking the initiation of transcriptional control and gene expression post-fertilization. Histone modifications are pivotal in shaping cellular identity and gene expression in many mammals. Recent advances in chromatin analysis have enabled detailed explorations of histone modifications during ZGA. This review delves into conserved and unique regulatory strategies, providing essential insights into the dynamic changes in histone modifications and their variants during ZGA in mammals. The objective is to explore recent advancements in leading mechanisms related to histone modifications governing this embryonic development phase in depth. These considerations will be useful for informing future therapeutic approaches that target epigenetic regulation in diverse biological contexts. It will also contribute to the extensive areas of evolutionary and developmental biology and possibly lay the foundation for future research and discussion on this seminal topic.

Shaping gene expression and its evolution by chromatin architecture and enhancer activity.

Transcriptional regulation plays a pivotal role in orchestrating the intricate genetic programs governing embryonic development. The expression of developmental genes relies on the combined activity of several cis-regulatory elements (CREs), such as enhancers and silencers, which can be located at long linear distances from the genes that they regulate and that interact with them through establishment of chromatin loops. Mutations affecting their activity or interaction with their target genes can lead to developmental disorders and are thought to have importantly contributed to the evolution of the animal body plan. The income of next-generation-sequencing approaches has allowed identifying over a million of sequences with putative regulatory potential in the human genome. Characterizing their function and establishing gene-CREs maps is essential to decode the logic governing developmental gene expression and is one of the major challenges of the post-genomic era. Chromatin 3D organization plays an essential role in determining how CREs specifically contact their target genes while avoiding deleterious off-target interactions. Our understanding of these aspects has greatly advanced with the income of chromatin conformation capture techniques and fluorescence microscopy approaches to visualize the organization of DNA elements in the nucleus. Here we will summarize relevant aspects of how the interplay between CRE activity and chromatin 3D organization regulates developmental gene expression and how it relates to pathological conditions and the evolution of animal body plan.

Hybrid Epigenomes Reveal Extensive Local Genetic Changes to Chromatin Accessibility Contribute to Divergence in Embryonic Gene Expression Between Species.

Chromatin accessibility plays an important role in shaping gene expression, yet little is known about the genetic and molecular mechanisms that influence the evolution of chromatin configuration. Both local (cis) and distant (trans) genetic influences can in principle influence chromatin accessibility and are based on distinct molecular mechanisms. We, therefore, sought to characterize the role that each of these plays in altering chromatin accessibility in 2 closely related sea urchin species. Using hybrids of Heliocidaris erythrogramma and Heliocidaris tuberculata, and adapting a statistical framework previously developed for the analysis of cis and trans influences on the transcriptome, we examined how these mechanisms shape the regulatory landscape at 3 important developmental stages, and compared our results to similar analyses of the transcriptome. We found extensive cis- and trans-based influences on evolutionary changes in chromatin, with cis effects generally larger in effect. Evolutionary changes in accessibility and gene expression are correlated, especially when expression has a local genetic basis. Maternal influences appear to have more of an effect on chromatin accessibility than on gene expression, persisting well past the maternal-to-zygotic transition. Chromatin accessibility near gene regulatory network genes appears to be distinctly regulated, with trans factors appearing to play an outsized role in the configuration of chromatin near these genes. Together, our results represent the first attempt to quantify cis and trans influences on evolutionary divergence in chromatin configuration in an outbred natural study system and suggest that chromatin regulation is more genetically complex than was previously appreciated.

Remodeling of the 3D chromatin architecture in the marine microalga Nannochloropsis oceanica during lipid accumulation

BackgroundGenomic three-dimensional (3D) spatial organization plays a key role in shaping gene expression and associated chromatin modification, and it is highly sensitive to environmental stress conditions. In microalgae, exposure to nitrogen stress can drive lipid accumulation, yet the associated functional alterations in the spatial organization of the microalgal genome have yet to be effectively characterized.ResultsAccordingly, the present study employed RNA-seq, Hi-C, and ChIP-seq approaches to explore the relationship between 3D chromosomal architecture and gene expression during lipid accumulation in the marine microalga Nannochloropsis oceanica in response to nitrogen deprivation (ND). These analyses revealed that ND resulted in various changes in chromosomal organization, including A/B compartment transitions, topologically associating domain (TAD) shifts, and the disruption of short-range interactions. Significantly higher levels of gene expression were evident in A compartments and TAD boundary regions relative to B compartments and TAD interior regions, consistent with observed histone modification enrichment in these areas. ND-induced differentially expressed genes (DEGs) were notably enriched in altered TAD-associated regions and regions exhibiting differential genomic contact. These DEGs were subjected to Gene Ontology (GO) term analyses that indicated they were enriched in the ‘fatty acid metabolism’, ‘response to stress’, ‘carbon fixation’ and ‘photosynthesis’ functional categories, in line with the ND treatment conditions used to conduct this study. These data indicate that Nannochloropsis cells exhibit a clear association between chromatin organization and transcriptional activity under nitrogen stress conditions. Pronounced and extensive histone modifications were evident in response to ND. Observed changes in chromatin architecture were linked to shifts in histone modifications and gene expression.ConclusionsOverall, the reprogramming of many lipid metabolism-associated genes was evident under nitrogen stress conditions with respect to both histone modifications and chromosomal organization. Together these results revealed that higher-order chromatin architecture represents a new layer that can guide efforts to understand the transcriptional regulation of lipid metabolism in nitrogen-deprived microalgae.

Brain Gene Expression of Foraging Behavior and Social Environment in Ceratina calcarata.

Rudimentary social systems have the potential to both advance our understanding of how complex sociality may have evolved and our understanding of how changes in social environment may influence gene expression and cooperation. Recently, studies of primitively social Hymenoptera have greatly expanded empirical evidence for the role of social environment in shaping behavior and gene expression. Here, we compare brain gene expression profiles of foragers across social contexts in the small carpenter bee, Ceratina calcarata. We conducted experimental manipulations of field colonies to examine gene expression profiles among social contexts including foraging mothers, regular daughters, and worker-like dwarf eldest daughters in the presence and absence of mother. Our analysis found significant differences in gene expression associated with female age, reproductive status, and social environment, including circadian clock gene dyw, hexamerin, and genes involved in the regulation of juvenile hormone and chemical communication. We also found that candidate genes differentially expressed in our study were also associated with division of labor, including foraging, in other primitively and advanced eusocial insects. Our results offer evidence for the role of the regulation of key developmental hormones and circadian rhythms in producing cooperative behavior in rudimentary insect societies.

The single-cell chromatin accessibility landscape in mouse perinatal testis development.

Spermatogenesis depends on an orchestrated series of developing events in germ cells and full maturation of the somatic microenvironment. To date, the majority of efforts to study cellular heterogeneity in testis has been focused on single-cell gene expression rather than the chromatin landscape shaping gene expression. To advance our understanding of the regulatory programs underlying testicular cell types, we analyzed single-cell chromatin accessibility profiles in more than 25,000 cells from mouse developing testis. We showed that single-cell sequencing assay for transposase-accessible chromatin (scATAC-Seq) allowed us to deconvolve distinct cell populations and identify cis-regulatory elements (CREs) underlying cell-type specification. We identified sets of transcription factors associated with cell type-specific accessibility, revealing novel regulators of cell fate specification and maintenance. Pseudotime reconstruction revealed detailed regulatory dynamics coordinating the sequential developmental progressions of germ cells and somatic cells. This high-resolution dataset also unveiled previously unreported subpopulations within both the Sertoli and Leydig cell groups. Further, we defined candidate target cell types and genes of several genome-wide association study (GWAS) signals, including those associated with testosterone levels and coronary artery disease. Collectively, our data provide a blueprint of the 'regulon' of the mouse male germline and supporting somatic cells.

FTO mediates LINE1 m6A demethylation and chromatin regulation in mESCs and mouse development.

N6-methyladenosine (m6A) is the most abundant internal modification on mammalian messenger RNA. It is installed by a writer complex and can be reversed by erasers such as the fat mass and obesity-associated protein FTO. Despite extensive research, the primary physiological substrates of FTO in mammalian tissues and development remain elusive. Here, we show that FTO mediates m6A demethylation of long-interspersed element-1 (LINE1) RNA in mouse embryonic stem cells (mESCs), regulating LINE1 RNA abundance and the local chromatin state, which in turn modulates the transcription of LINE1-containing genes. FTO-mediated LINE1 RNA m6A demethylation also plays regulatory roles in shaping chromatin state and gene expression during mouse oocyte and embryonic development. Our results suggest broad effects of LINE1 RNA m6A demethylation by FTO in mammals.

Regnase‐1–related endoribonucleases in health and immunological diseases

Dynamic changes in gene expression are key factors in the development and activation of immune cells. RNA metabolism is one of the critical steps for the control of gene expression. Together with transcriptional regulation, mRNA decay by specific ribonucleases (RNases) plays a vital role in shaping gene expression. In addition to the canonical exoribonuclease-mediated mRNA degradation through the recognition of cis-elements in mRNA 3' untranslated regions by RNA-binding proteins (RBPs), endoribonucleases are involved in the control of mRNAs in immune cells. In this review, we gleam insights on how Regnase-1, an endoribonuclease necessary for regulating immune cell activation and maintenance of immune homeostasis, degrades RNAs involved in immune cell activation. Additionally, we provide insights on recent studies which uncover the role of Regnase-1-related RNases, including Regnase-2, Regnase-3, and Regnase-4, as well as N4BP1 and KHNYN, in immune regulation and antiviral immunity. As the dysregulation of immune mRNA decay leads to pathologies such as autoimmune diseases or impaired activation of immune responses, RNases are deemed as essential components of regulatory feedback mechanisms that modulate inflammation. Given the critical role of RNases in autoimmunity, RNases can be perceived as emerging targets in the development of novel therapeutics.

Current knowledge on genetic variants shaping placental transcriptome and their link to gestational and postnatal health.

Despite the indispensable role of the placenta in the successful course of pregnancy, regulation of its dynamic transcriptome is still underexplored. The purpose of this literature review was to give an overview and draw attention to the contribution of genetic variation in shaping the human placental gene expression. Studies of placental transcriptome shaped by chromosomal variants are limited and may be confounded by cellular mosaicism and somatic genomic rearrangements. Even in relatively simple cases, such as aneuploidies, the placental transcriptome appears to differ from the assumed systematically increased transcript levels of the involved chromosomes. Single nucleotide variants modulating placental gene expression referred to as expression quantitative trait loci (eQTLs) have been analyzed only in ten candidate gene and three genome-wide association studies (GWAS). The latter identified 417 confident placental eGenes, supported by at least two independent studies. Functional profiling of eGenes highlighted biological pathways important in pregnancy, such as immune response or transmembrane transport activity. A fraction of placental eQTLs (1-3%) co-localize with GWAS loci for adult disorders (metabolic, immunological, neurological), suggesting a co-contributory role of the placenta in the developmental programming of health. Some placental eQTLs have been identified as risk factors for adverse pregnancy outcomes, such as rs4769613 (C>T), located near the FLT1 gene and confidently associated with preeclampsia. More studies are needed to map genetic variants shaping gene expression in different placental cell types across three trimesters in normal and complicated gestations and to clarify to what extent these heritable factors contribute to maternal and offspring disease risks.

Extensive fragmentation and re-organization of transcription in Systemic Lupus Erythematosus

Systemic Lupus Erythematosus (SLE) is the prototype of autoimmune diseases, characterized by extensive gene expression perturbations in peripheral blood immune cells. Circumstantial evidence suggests that these perturbations may be due to altered epigenetic profiles and chromatin accessibility but the relationship between transcriptional deregulation and genome organization remains largely unstudied. In this work we propose a genomic approach that leverages patterns of gene coexpression from genome-wide transcriptome profiles in order to identify statistically robust Domains of Co-ordinated gene Expression (DCEs). Application of this method on a large transcriptome profiling dataset of 148 SLE patients and 52 healthy individuals enabled the identification of significant disease-associated alterations in gene co-regulation patterns, which also correlate with SLE activity status. Low disease activity patient genomes are characterized by extensive fragmentation leading to overall fewer DCEs of smaller size. High disease activity genomes display extensive redistribution of co-expression domains with expanded and newly-appearing (emerged) DCEs. The dynamics of domain fragmentation and redistribution are associated with SLE clinical endophenotypes, with genes of the interferon pathway being highly enriched in DCEs that become disrupted and with functions associated to more generalized symptoms, being located in domains that emerge anew in high disease activity genomes. Our results suggest strong links between the SLE phenotype and the underlying genome structure and underline an important role for genome organization in shaping gene expression in SLE.

Metabolic labeling of RNA using multiple ribonucleoside analogs enables the simultaneous evaluation of RNA synthesis and degradation rates.

Gene expression is determined by a balance between RNA synthesis and RNA degradation. To elucidate the underlying regulatory mechanisms and principles of this, simultaneous measurements of RNA synthesis and degradation are required. Here, we report the development of “Dyrec-seq,” which uses 4-thiouridine and 5-bromouridine to simultaneously quantify RNA synthesis and degradation rates. Dyrec-seq enabled the quantification of RNA synthesis and degradation rates of 4702 genes in HeLa cells. Functional enrichment analysis showed that the RNA synthesis and degradation rates of genes are actually determined by the genes’ biological functions. A comparison of theoretical and experimental analyses revealed that the amount of RNA is determined by the ratio of RNA synthesis to degradation rates, whereas the rapidity of responses to external stimuli is determined only by the degradation rate. This study emphasizes that not only RNA synthesis but also RNA degradation is important in shaping gene expression patterns.

Comparative analysis of diet-associated responses in two rice planthopper species

BackgroundHost adaptation is the primary determinant of insect diversification. However, knowledge of different host ranges in closely related species remains scarce. The brown planthopper (Nilaparvata lugens, BPH) and the small brown planthopper (Laodelphax striatellus, SBPH) are the most destructive insect pests within the family Delphacidae. These two species differ in their host range (SBPH can well colonize rice and wheat plants, whereas BPH survives on only rice plants), but the underlying mechanism of this difference remains unknown. High-throughput sequencing provides a powerful approach for analyzing the association between changes in gene expression and the physiological responses of insects. Therefore, gut transcriptomes were performed to elucidate the genes associated with host adaptation in planthoppers. The comparative analysis of planthopper responses to different diets will improve our knowledge of host adaptation regarding herbivorous insects.ResultsIn the present study, we analyzed the change in gene expression of SBPHs that were transferred from rice plants to wheat plants over the short term (rSBPH vs tSBPH) or were colonized on wheat plants over the long term (rSBPH vs wSBPH). The results showed that the majority of differentially expressed genes in SBPH showed similar changes in expression for short-term transfer and long-term colonization. Based on a comparative analysis of BPH and SBPH after transfer, the genes associated with sugar transporters and heat-shock proteins showed similar variation. However, most of the genes were differentially regulated between the two species. The detoxification-related genes were upregulated in SBPH after transfer from the rice plants to the wheat plants, but these genes were downregulated in BPH under the same conditions. In contrast, ribosomal-related genes were downregulated in SBPH after transfer, but these genes were upregulated in BPH under the same conditions.ConclusionsThe results of this study provide evidence that host plants played a dominant role in shaping gene expression and that the low fitness of BPH on wheat plants might be determined within 24 h after transfer. This study deepens our understanding of different host ranges for the two planthopper species, which may provide a potential strategy for pest management.

Homoeologous exchanges occur through intragenic recombination generating novel transcripts and proteins in wheat and other polyploids

Recombination between homeologous chromosomes, also known as homeologous exchange (HE), plays a significant role in shaping genome structure and gene expression in interspecific hybrids and allopolyploids of several plant species. However, the molecular mechanisms that govern HEs are not well understood. Here, we studied HE events in the progeny of a nascent allotetraploid (genome AADD) derived from two diploid progenitors of hexaploid bread wheat using cytological and whole-genome sequence analyses. In total, 37 HEs were identified and HE junctions were mapped precisely. HEs exhibit typical patterns of homologous recombination hotspots, being biased toward low-copy, subtelomeric regions of chromosome arms and showing association with known recombination hotspot motifs. But, strikingly, while homologous recombination preferentially takes place upstream and downstream of coding regions, HEs are highly enriched within gene bodies, giving rise to novel recombinant transcripts, which in turn are predicted to generate new protein fusion variants. To test whether this is a widespread phenomenon, a dataset of high-resolution HE junctions was analyzed for allopolyploid Brassica, rice, Arabidopsis suecica, banana, and peanut. Intragenic recombination and formation of chimeric genes was detected in HEs of all species and was prominent in most of them. HE thus provides a mechanism for evolutionary novelty in transcript and protein sequences in nascent allopolyploids.

Hfq CLASH uncovers sRNA-target interaction networks linked to nutrient availability adaptation.

By shaping gene expression profiles, small RNAs (sRNAs) enable bacteria to efficiently adapt to changes in their environment. To better understand how Escherichia coli acclimatizes to nutrient availability, we performed UV cross-linking, ligation and sequencing of hybrids (CLASH) to uncover Hfq-associated RNA-RNA interactions at specific growth stages. We demonstrate that Hfq CLASH robustly captures bona fide RNA-RNA interactions. We identified hundreds of novel sRNA base-pairing interactions, including many sRNA-sRNA interactions and involving 3'UTR-derived sRNAs. We rediscovered known and identified novel sRNA seed sequences. The sRNA-mRNA interactions identified by CLASH have strong base-pairing potential and are highly enriched for complementary sequence motifs, even those supported by only a few reads. Yet, steady state levels of most mRNA targets were not significantly affected upon over-expression of the sRNA regulator. Our results reinforce the idea that the reproducibility of the interaction, not base-pairing potential, is a stronger predictor for a regulatory outcome.

Species-Specific miRNAs in Human Brain Development and Disease.

Identification of the unique features of human brain development and function can be critical towards the elucidation of intricate processes such as higher cognitive functions and human-specific pathologies like neuropsychiatric and behavioral disorders. The developing primate and human central nervous system (CNS) are distinguished by expanded progenitor zones and a protracted time course of neurogenesis, leading to the expansion in brain size, prominent gyral anatomy, distinctive synaptic properties, and complex neural circuits. Comparative genomic studies have revealed that adaptations of brain capacities may be partly explained by human-specific genetic changes that impact the function of proteins associated with neocortical expansion, synaptic function, and language development. However, the formation of complex gene networks may be most relevant for brain evolution. Indeed, recent studies identified distinct human-specific gene expression patterns across developmental time occurring in brain regions linked to cognition. Interestingly, such modules show species-specific divergence and are enriched in genes associated with neuronal development and synapse formation whilst also being implicated in neuropsychiatric diseases. microRNAs represent a powerful component of gene-regulatory networks by promoting spatiotemporal post-transcriptional control of gene expression in the human and primate brain. It has also been suggested that the divergence in miRNA expression plays an important role in shaping gene expression divergence among species. Primate-specific and human-specific miRNAs are principally involved in progenitor proliferation and neurogenic processes but also associate with human cognition, and neurological disorders. Human embryonic or induced pluripotent stem cells and brain organoids, permitting experimental access to neural cells and differentiation stages that are otherwise difficult or impossible to reach in humans, are an essential means for studying species-specific brain miRNAs. Single-cell sequencing approaches can further decode refined miRNA-mRNA interactions during developmental transitions. Elucidating species-specific miRNA regulation will shed new light into the mechanisms that control spatiotemporal events during human brain development and disease, an important step towards fostering novel, holistic and effective therapeutic approaches for neural disorders. In this review, we discuss species-specific regulation of miRNA function, its contribution to the evolving features of the human brain and in neurological disease, with respect also to future therapeutic approaches.

Integrating distal and proximal information to predict gene expression via a densely connected convolutional neural network.

Interactions among cis-regulatory elements such as enhancers and promoters are main driving forces shaping context-specific chromatin structure and gene expression. Although there have been computational methods for predicting gene expression from genomic and epigenomic information, most of them neglect long-range enhancer-promoter interactions, due to the difficulty in precisely linking regulatory enhancers to target genes. Recently, HiChIP, a novel high-throughput experimental approach, has generated comprehensive data on high-resolution interactions between promoters and distal enhancers. Moreover, plenty of studies suggest that deep learning achieves state-of-the-art performance in epigenomic signal prediction, and thus promoting the understanding of regulatory elements. In consideration of these two factors, we integrate proximal promoter sequences and HiChIP distal enhancer-promoter interactions to accurately predict gene expression. We propose DeepExpression, a densely connected convolutional neural network, to predict gene expression using both promoter sequences and enhancer-promoter interactions. We demonstrate that our model consistently outperforms baseline methods, not only in the classification of binary gene expression status but also in regression of continuous gene expression levels, in both cross-validation experiments and cross-cell line predictions. We show that the sequential promoter information is more informative than the experimental enhancer information; meanwhile, the enhancer-promoter interactions within ±100 kbp around the TSS of a gene are most beneficial. We finally visualize motifs in both promoter and enhancer regions and show the match of identified sequence signatures with known motifs. We expect to see a wide spectrum of applications using HiChIP data in deciphering the mechanism of gene regulation. DeepExpression is freely available at https://github.com/wanwenzeng/DeepExpression. Supplementary data are available at Bioinformatics online.

Genome wide analysis of 3' UTR sequence elements and proteins regulating mRNA stability during maternal-to-zygotic transition in zebrafish.

Posttranscriptional regulation plays a crucial role in shaping gene expression. During the maternal-to-zygotic transition (MZT), thousands of maternal transcripts are regulated. However, how different cis-elements and trans-factors are integrated to determine mRNA stability remains poorly understood. Here, we show that most transcripts are under combinatorial regulation by multiple decay pathways during zebrafish MZT. By using a massively parallel reporter assay, we identified cis-regulatory sequences in the 3′ UTR, including U-rich motifs that are associated with increased mRNA stability. In contrast, miR-430 target sequences, UAUUUAUU AU-rich elements (ARE), CCUC, and CUGC elements emerged as destabilizing motifs, with miR-430 and AREs causing mRNA deadenylation upon genome activation. We identified trans-factors by profiling RNA–protein interactions and found that poly(U)-binding proteins are preferentially associated with 3′ UTR sequences and stabilizing motifs. We show that this activity is antagonized by C-rich motifs and correlated with protein binding. Finally, we integrated these regulatory motifs into a machine learning model that predicts reporter mRNA stability in vivo.

A pan-cancer analysis of driver gene mutations, DNA methylation and gene expressions reveals that chromatin remodeling is a major mechanism inducing global changes in cancer epigenomes

BackgroundRecent large-scale cancer sequencing studies have discovered many novel cancer driver genes (CDGs) in human cancers. Some studies also suggest that CDG mutations contribute to cancer-associated epigenomic and transcriptomic alterations across many cancer types. Here we aim to improve our understanding of the connections between CDG mutations and altered cancer cell epigenomes and transcriptomes on pan-cancer level and how these connections contribute to the known association between epigenome and transcriptome.MethodUsing multi-omics data including somatic mutation, DNA methylation, and gene expression data of 20 cancer types from The Cancer Genome Atlas (TCGA) project, we conducted a pan-cancer analysis to identify CDGs, when mutated, have strong associations with genome-wide methylation or expression changes across cancer types, which we refer as methylation driver genes (MDGs) or expression driver genes (EDGs), respectively.ResultsWe identified 32 MDGs, among which, eight are known chromatin modification or remodeling genes. Many of the remaining 24 MDGs are connected to chromatin regulators through either regulating their transcription or physically interacting with them as potential co-factors. We identified 29 EDGs, 26 of which are also MDGs. Further investigation on target genes’ promoters methylation and expression alteration patterns of these 26 overlapping driver genes shows that hyper-methylation of target genes’ promoters are significantly associated with down-regulation of the same target genes and hypo-methylation of target genes’ promoters are significantly associated with up-regulation of the same target genes.ConclusionThis finding suggests a pivotal role for genetically driven changes in chromatin remodeling in shaping DNA methylation and gene expression patterns during tumor development.