Eukaryotic Gene Regulation Research Articles

Sumoylation of the human histone H4 tail inhibits p300-mediated transcription by RNA polymerase II in cellular extracts

The post-translational modification of histones by the small ubiquitin-like modifier (SUMO) protein has been associated with gene regulation, centromeric localization, and double-strand break repair in eukaryotes. Although sumoylation of histone H4 was specifically associated with gene repression, this could not be proven due to the challenge of site-specifically sumoylating H4 in cells. Biochemical crosstalk between SUMO and other histone modifications, such as H4 acetylation and H3 methylation, that are associated with active genes also remains unclear. We addressed these challenges in mechanistic studies using an H4 chemically modified at Lys12 by SUMO-3 (H4K12su) and incorporated into mononucleosomes and chromatinized plasmids for functional studies. Mononucleosome-based assays revealed that H4K12su inhibits transcription-activating H4 tail acetylation by the histone acetyltransferase p300, as well as transcription-associated H3K4 methylation by the extended catalytic module of the Set1/COMPASS (complex of proteins associated with Set1) histone methyltransferase complex. Activator- and p300-dependent in vitro transcription assays with chromatinized plasmids revealed that H4K12su inhibits both H4 tail acetylation and RNA polymerase II-mediated transcription. Finally, cell-based assays with a SUMO-H4 fusion that mimics H4 tail sumoylation confirmed the negative crosstalk between histone sumoylation and acetylation/methylation. Thus, our studies establish the key role for histone sumoylation in gene silencing and its negative biochemical crosstalk with active transcription-associated marks in human cells.

Identification and Analysis of Small Interfering RNAs Associated With Heat Stress in Flowering Chinese Cabbage Using High-Throughput Sequencing.

Endogenous small interfering RNAs (siRNAs) are substantial gene regulators in eukaryotes and play key functions in plant development and stress tolerance. Among environmental factors, heat is serious abiotic stress that severely influences the productivity and quality of flowering Chinese cabbage (Brassica campestris L. ssp. chinensis var. utilis Tsen et Lee). However, how siRNAs are involved in regulating gene expression during heat stress is not fully understood in flowering Chinese cabbage. Combining bioinformatical and next-generation sequencing approaches, we identified heat-responsive siRNAs in four small RNA libraries of flowering Chinese cabbage using leaves collected at 0, 1, 6, and 12 h after a 38°C heat-stress treatment; 536, 816, and 829 siRNAs exhibited substantial differential expression at 1, 6, and 12 h, respectively. Seventy-five upregulated and 69 downregulated differentially expressed siRNAs (DE-siRNAs) were common for the three time points of heat stress. We identified 795 target genes of DE-siRNAs, including serine/threonine-protein kinase SRK2I, CTR1-like, disease resistance protein RML1A-like, and RPP1, which may play a role in regulating heat tolerance. Gene ontology showed that predictive targets of DE-siRNAs may have key roles in the positive regulation of biological processes, organismal processes, responses to temperature stimulus, signaling, and growth and development. These novel results contribute to further understanding how siRNAs modulate the expression of their target genes to control heat tolerance in flowering Chinese cabbage.

Mapping the dynamic transfer functions of eukaryotic gene regulation.

Biological information can be encoded within the dynamics of signaling components, which has been implicated in a broad range of physiological processes including stress response, oncogenesis, and stem cell differentiation. To study the complexity of information transfer across the eukaryotic promoter, we screened 119 dynamic conditions-modulating the pulse frequency, amplitude, and pulse width of light-regulating the binding of an epigenome editor to a fluorescent reporter. This system revealed tunable gene expression and filtering behaviors and provided a quantification of the limit to the amount of information that can be reliably transferred across a single promoter as ∼1.7 bits. Using a library of over 100 orthogonal chromatin regulators, we further determined that chromatin state could be used to tune mutual information and expression levels, as well as completely alter the input-output transfer function of the promoter. This system unlocks the information-rich content of eukaryotic gene regulation.

Evaluation of Lee etal.: Clarity and interpretation of mutual information in promoter transfer functions.

One snapshot of the peer review process for "Mapping the dynamic transfer functions of eukaryotic gene regulation" (Lee etal., 2021) appears below.

Proximal-end bias from in-vitro reconstituted nucleosomes and the result on downstream data analysis.

The most basic level of eukaryotic gene regulation is the presence or absence of nucleosomes on DNA regulatory elements. In an effort to elucidate in vivo nucleosome patterns, in vitro studies are frequently used. In vitro, short DNA fragments are more favorable for nucleosome formation, increasing the likelihood of nucleosome occupancy. This may in part result from the fact that nucleosomes prefer to form on the terminal ends of linear DNA. This phenomenon has the potential to bias in vitro reconstituted nucleosomes and skew results. If the ends of DNA fragments are known, the reads falling close to the ends are typically discarded. In this study we confirm the phenomenon of end bias of in vitro nucleosomes. We describe a method in which nearly identical libraries, with different known ends, are used to recover nucleosomes which form towards the terminal ends of fragmented DNA. Finally, we illustrate that although nucleosomes prefer to form on DNA ends, it does not appear to skew results or the interpretation thereof.

ScAPAatlas: an atlas of alternative polyadenylation across cell types in human and mouse.

Alternative polyadenylation (APA) has been widely recognized as a crucial step during the post-transcriptional regulation of eukaryotic genes. Recent studies have demonstrated that APA exerts key regulatory roles in many biological processes and often occurs in a tissue- and cell-type-specific manner. However, to our knowledge, there is no database incorporating information about APA at the cell-type level. Single-cell RNA-seq is a rapidly evolving and powerful tool that enable APA analysis at the cell-type level. Here, we present a comprehensive resource, scAPAatlas (http://www.bioailab.com:3838/scAPAatlas), for exploring APA across different cell types, and interpreting potential biological functions. Based on the curated scRNA-seq data from 24 human and 25 mouse normal tissues, we systematically identified cell-type-specific APA events for different cell types and examined the correlations between APA and gene expression level. We also estimated the crosstalk between cell-type-specific APA events and microRNAs or RNA-binding proteins. A user-friendly web interface has been constructed to support browsing, searching and visualizing multi-layer information of cell-type-specific APA events. Overall, scAPAatlas, incorporating a rich resource for exploration of APA at the cell-type level, will greatly help researchers chart cell type with APA and elucidate the biological functions of APA.

Principles of mRNA targeting via the Arabidopsis m6A-binding protein ECT2.

Specific recognition of N6-methyladenosine (m6A) in mRNA by RNA-binding proteins containing a YT521-B homology (YTH) domain is important in eukaryotic gene regulation. The Arabidopsis YTH domain protein ECT2 is thought to bind to mRNA at URU(m6A)Y sites, yet RR(m6A)CH is the canonical m6A consensus site in all eukaryotes and ECT2 functions require m6A-binding activity. Here, we apply iCLIP (individual nucleotide resolution crosslinking and immunoprecipitation) and HyperTRIBE (targets of RNA-binding proteins identified by editing) to define high-quality target sets of ECT2 and analyze the patterns of enriched sequence motifs around ECT2 crosslink sites. Our analyses show that ECT2 does in fact bind to RR(m6A)CH. Pyrimidine-rich motifs are enriched around, but not at m6A sites, reflecting a preference for N6-adenosine methylation of RRACH/GGAU islands in pyrimidine-rich regions. Such motifs, particularly oligo-U and UNUNU upstream of m6A sites, are also implicated in ECT2 binding via its intrinsically disordered region (IDR). Finally, URUAY-type motifs are enriched at ECT2 crosslink sites, but their distinct properties suggest function as sites of competition between binding of ECT2 and as yet unidentified RNA-binding proteins. Our study provides coherence between genetic and molecular studies of m6A-YTH function in plants and reveals new insight into the mode of RNA recognition by YTH domain-containing proteins.

Position Specific Alternative Splicing and Gene Expression Profiles Along the Tonotopic Axis of Chick Cochlea.

Alternative splicing (AS) refers to the production of multiple mRNA isoforms from a single gene due to alternative selection of exons or splice sites during pre-mRNA splicing. It is a primary mechanism of gene regulation in higher eukaryotes and significantly expands the functional complexity of eukaryotic organisms, contributing to animal development and disease. Recent studies have shown that AS also influences functional diversity by affecting the transcriptomic and proteomic profiles in a position-dependent manner in a single organ. The peripheral hearing organ, the cochlea, is organized to detect sounds at different frequencies depending on its location along the longitudinal axis. This unique functional configuration, the tonotopy, is known to be facilitated by differential gene expression along the cochlear duct. We profiled transcriptome-wide gene expression and AS changes that occur within the different positions of chick cochlea. These analyses revealed distinct gene expression profiles and AS, including a splicing program that is unique to tonotopy. Changes in the expression of splicing factors PTBP3, ESRP1, and ESRP2 were demonstrated to contribute to position-specific AS. RNA-binding motif enrichment analysis near alternatively spliced exons provided further insight into the combinatorial regulation of AS at different positions by different RNA-binding proteins. These data, along with gene ontology (GO) analysis, represent a comprehensive analysis of the dynamic regulation of AS at different positions in chick cochlea.

Genome-wide identification of DCL, AGO and RDR gene families and their associated functional regulatory elements analyses in banana (Musa acuminata).

RNA silencing is mediated through RNA interference (RNAi) pathway gene families, i.e., Dicer-Like (DCL), Argonaute (AGO), and RNA-dependent RNA polymerase (RDR) and their cis-acting regulatory elements. The RNAi pathway is also directly connected with the post-transcriptional gene silencing (PTGS) mechanism, and the pathway controls eukaryotic gene regulation during growth, development, and stress response. Nevertheless, genome-wide identification of RNAi pathway gene families such as DCL, AGO, and RDR and their regulatory network analyses related to transcription factors have not been studied in many fruit crop species, including banana (Musa acuminata). In this study, we studied in silico genome-wide identification and characterization of DCL, AGO, and RDR genes in bananas thoroughly via integrated bioinformatics approaches. A genome-wide analysis identified 3 MaDCL, 13 MaAGO, and 5 MaRDR candidate genes based on multiple sequence alignment and phylogenetic tree related to the RNAi pathway in banana genomes. These genes correspond to the Arabidopsis thaliana RNAi silencing genes. The analysis of the conserved domain, motif, and gene structure (exon-intron numbers) for MaDCL, MaAGO, and MaRDR genes showed higher homogeneity within the same gene family. The Gene Ontology (GO) enrichment analysis exhibited that the identified RNAi genes could be involved in RNA silencing and associated metabolic pathways. A number of important transcription factors (TFs), e.g., ERF, Dof, C2H2, TCP, GATA and MIKC_MADS families, were identified by network and sub-network analyses between TFs and candidate RNAi gene families. Furthermore, the cis-acting regulatory elements related to light-responsive (LR), stress-responsive (SR), hormone-responsive (HR), and other activities (OT) functions were identified in candidate MaDCL, MaAGO, and MaRDR genes. These genome-wide analyses of these RNAi gene families provide valuable information related to RNA silencing, which would shed light on further characterization of RNAi genes, their regulatory elements, and functional roles, which might be helpful for banana improvement in the breeding program.

Allosteric conformational ensembles have unlimited capacity for integrating information.

Integration of binding information by macromolecular entities is fundamental to cellular functionality. Recent work has shown that such integration cannot be explained by pairwise cooperativities, in which binding is modulated by binding at another site. Higher-order cooperativities (HOCs), in which binding is collectively modulated by multiple other binding events, appear to be necessary but an appropriate mechanism has been lacking. We show here that HOCs arise through allostery, in which effective cooperativity emerges indirectly from an ensemble of dynamically interchanging conformations. Conformational ensembles play important roles in many cellular processes but their integrative capabilities remain poorly understood. We show that sufficiently complex ensembles can implement any form of information integration achievable without energy expenditure, including all patterns of HOCs. Our results provide a rigorous biophysical foundation for analysing the integration of binding information through allostery. We discuss the implications for eukaryotic gene regulation, where complex conformational dynamics accompanies widespread information integration.

Quantifying splice-site usage: a simple yet powerful approach to analyze splicing.

RNA splicing, and variations in this process referred to as alternative splicing, are critical aspects of gene regulation in eukaryotes. From environmental responses in plants to being a primary link between genetic variation and disease in humans, splicing differences confer extensive phenotypic changes across diverse organisms (1–3). Regulation of splicing occurs through differential selection of splice sites in a splicing reaction, which results in variation in the abundance of isoforms and/or splicing events. However, genomic determinants that influence splice-site selection remain largely unknown. While traditional approaches for analyzing splicing rely on quantifying variant transcripts (i.e. isoforms) or splicing events (i.e. intron retention, exon skipping etc.) (4), recent approaches focus on analyzing complex/mutually exclusive splicing patterns (5–8). However, none of these approaches explicitly measure individual splice-site usage, which can provide valuable information about splice-site choice and its regulation. Here, we present a simple approach to quantify the empirical usage of individual splice sites reflecting their strength, which determines their selection in a splicing reaction. Splice-site strength/usage, as a quantitative phenotype, allows us to directly link genetic variation with usage of individual splice-sites. We demonstrate the power of this approach in defining the genomic determinants of splice-site choice through GWAS. Our pilot analysis with more than a thousand splice sites hints that sequence divergence in cis rather than trans is associated with variations in splicing among accessions of Arabidopsis thaliana. This approach allows deciphering principles of splicing and has broad implications from agriculture to medicine.

Conserved dual-mode gene regulation programs in higher eukaryotes.

Recent genomic data analyses have revealed important underlying logics in eukaryotic gene regulation, such as CpG islands (CGIs)-dependent dual-mode gene regulation. In mammals, genes lacking CGIs at their promoters are generally regulated by interconversion between euchromatin and heterochromatin, while genes associated with CGIs constitutively remain as euchromatin. Whether a similar mode of gene regulation exists in non-mammalian species has been unknown. Here, through comparative epigenomic analyses, we demonstrate that the dual-mode gene regulation program is common in various eukaryotes, even in the species lacking CGIs. In cases of vertebrates or plants, we find that genes associated with high methylation level promoters are inactivated by forming heterochromatin and expressed in a context-dependent manner. In contrast, the genes with low methylation level promoters are broadly expressed and remain as euchromatin even when repressed by Polycomb proteins. Furthermore, we show that invertebrate animals lacking DNA methylation, such as fruit flies and nematodes, also have divergence in gene types: some genes are regulated by Polycomb proteins, while others are regulated by heterochromatin formation. Altogether, our study establishes gene type divergence and the resulting dual-mode gene regulation as fundamental features shared in a broad range of higher eukaryotic species.

Slaying the last unicorn: discovery of histones in the microalga Nanochlorum eucaryotum.

Histones are the principal constituents of eukaryotic chromatin. The four core histones (H2A, H2B, H3 and H4) are conserved across sequenced eukaryotic genomes and therefore thought to be universal to eukaryotes. In the early 1980s, however, a series of biochemical investigations failed to find evidence for histones or nucleosomal structures in the microscopic green alga Nanochlorum eucaryotum. If true, derived histone loss in this lineage would constitute an exceptional case that might help us further understand the principles governing eukaryotic gene regulation. To substantiate these earlier reports of histone loss in N. eucaryotum, we sequenced, assembled and quantified its transcriptome. Following a systematic search for histone-fold domains in the assembled transcriptome, we detect orthologues to all four core histones. We also find histone mRNAs to be highly expressed, comparable to the situation in other eukaryotes. Finally, we obtain characteristic protection patterns when N. eucaryotum chromatin is subjected to micrococcal nuclease digestion, indicating widespread formation of nucleosomal complexes in vivo. We conclude that previous reports of missing histones in N. eucaryotum were mistaken. By all indications, Nanochlorum eucaryotum has histone-based chromatin characteristic of most eukaryotes.

Fine gene expression regulation by minor sequence variations downstream of the polyadenylation signal.

The termination of transcription is a complex process that substantially contributes to gene regulation in eukaryotes. Previously, it was noted that a single cytosine deletion at the position + 32bp relative to the single polyadenylation signal AAUAAA (hereafter the dC mutation) causes a 2-fold increase in the transcription level of the upstream eGFP reporter in mouse embryonic stem cells. Here, we analyzed the conservation of this phenomenon in immortalized mouse, human and drosophila cell lines and the influence of the dC mutation on the choice of the pre-mRNA cleavage sites. We have constructed dual-reporter plasmids to accurately measure the effect of the dC and other nearby located mutations on eGFP mRNA level by RT-qPCR. In this way, we found that the dC mutation leads to a 2-fold increase in the expression level of the upstream eGFP reporter gene in cultured mouse and human, but not in drosophila cells. In addition, 3' RACE analysis demonstrated that eGFP pre-mRNAs are cut at multiple positions between + 14 to + 31, and that the most proximal cleavage site becomes almost exclusively utilized in the presence of the dC mutation. We alsoidentified new short sequence variations located within positions + 25.. + 40 and + 33.. + 48 that increase eGFP expression up to ~2-4-fold. Altogether, the positive effect of the dC mutation seems to be conserved in mouse embryonic stem cells, mouse embryonic 3T3 fibroblasts and human HEK293T cells. In the latter cells, the dC mutation appears to be involved in regulating pre-mRNA cleavage site selection. Finally, a multiplexed approach is proposed to identify motifs located downstream of cleavage site(s) that are essential for transcription termination.

Functional and transcriptional profiling of non-coding RNAs in yeast reveal context-dependent phenotypes and in trans effects on the protein regulatory network.

Non-coding RNAs (ncRNAs), including the more recently identified Stable Unannotated Transcripts (SUTs) and Cryptic Unstable Transcripts (CUTs), are increasingly being shown to play pivotal roles in the transcriptional and post-transcriptional regulation of genes in eukaryotes. Here, we carried out a large-scale screening of ncRNAs in Saccharomyces cerevisiae, and provide evidence for SUT and CUT function. Phenotypic data on 372 ncRNA deletion strains in 23 different growth conditions were collected, identifying ncRNAs responsible for significant cellular fitness changes. Transcriptome profiles were assembled for 18 haploid ncRNA deletion mutants and 2 essential ncRNA heterozygous deletants. Guided by the resulting RNA-seq data we analysed the genome-wide dysregulation of protein coding genes and non-coding transcripts. Novel functional ncRNAs, SUT125, SUT126, SUT035 and SUT532 that act in trans by modulating transcription factors were identified. Furthermore, we described the impact of SUTs and CUTs in modulating coding gene expression in response to different environmental conditions, regulating important biological process such as respiration (SUT125, SUT126, SUT035, SUT432), steroid biosynthesis (CUT494, SUT053, SUT468) or rRNA processing (SUT075 and snR30). Overall, these data capture and integrate the regulatory and phenotypic network of ncRNAs and protein-coding genes, providing genome-wide evidence of the impact of ncRNAs on cellular homeostasis.

Advances of functional consequences and regulation mechanisms of alternative cleavage and polyadenylation.

During the maturation of pre-mRNAs and some lncRNAs, their 3'ends are cleaved and polyadenylated. The cleavage and polyadenylation (C/P) require the presence of a polyadenylation signal (PAS) at the RNA 3?end. Most eukaryotic genes have multiple PASs, resulting in alternative cleavage and polyadenylation (APA). APA leads to transcript isoforms with different coding potentials and/or variable 3?UTRs. The 3'UTR affects mRNA stability, translation, transportation, and cellular localization. Therefore, APA is an important mechanism of posttranscriptional gene regulation in eukaryotes. In recent years, whole genome sequencing of animals, plants and yeast has revealed that APA is pervasive in eukaryotes, and the functional consequences and regulation of APA have been studied. To date, many cis-acting regulatory elements and trans-acting factors for APA regulation have been identified. In this review, we summarize the recent advances in the functional consequences and regulation of APA and discuss the future directions, aiming to provide clues and references for future APA study.

BET inhibition disrupts transcription but retains enhancer-promoter contact

Enhancers are DNA sequences that enable complex temporal and tissue-specific regulation of genes in higher eukaryotes. Although it is not entirely clear how enhancer-promoter interactions can increase gene expression, this proximity has been observed in multiple systems at multiple loci and is thought to be essential for the maintenance of gene expression. Bromodomain and Extra-Terminal domain (BET) and Mediator proteins have been shown capable of forming phase condensates and are thought to be essential for super-enhancer function. Here, we show that targeting of cells with inhibitors of BET proteins or pharmacological degradation of BET protein Bromodomain-containing protein 4 (BRD4) has a strong impact on transcription but very little impact on enhancer-promoter interactions. Dissolving phase condensates reduces BRD4 and Mediator binding at enhancers and can also strongly affect gene transcription, without disrupting enhancer-promoter interactions. These results suggest that activation of transcription and maintenance of enhancer-promoter interactions are separable events. Our findings further indicate that enhancer-promoter interactions are not dependent on high levels of BRD4 and Mediator, and are likely maintained by a complex set of factors including additional activator complexes and, at some sites, CTCF and cohesin.

Detecting Internal N7-Methylguanosine mRNA Modifications by Differential Enzymatic Digestion Coupled with Mass Spectrometry Analysis.

The recent discovery of reversible chemical modifications on mRNA has opened a new era of post-transcriptional gene regulation in eukaryotes. Among these modifications identified in eukaryotic mRNA, N7-methylguanosine (m7G) is unique owing to its presence in the 5' cap structure. Recently, it has been reported that m7G also exists internally in mRNA. Here, we describe a protocol of combining differential enzymatic digestion with liquid chromatography-electrospray ionization-tandem mass spectrometry (LC-ESI-MS/MS) analysis to detect internal m7G modification in mRNA. This protocol can also be used to quantify the level of m7G at both the 5' cap and internal positions of mRNA.

Chemical Decorations of "MARs" Residents in Orchestrating Eukaryotic Gene Regulation.

Genome organization plays a crucial role in gene regulation, orchestrating multiple cellular functions. A meshwork of proteins constituting a three-dimensional (3D) matrix helps in maintaining the genomic architecture. Sequences of DNA that are involved in tethering the chromatin to the matrix are called scaffold/matrix attachment regions (S/MARs), and the proteins that bind to these sequences and mediate tethering are termed S/MAR-binding proteins (S/MARBPs). The regulation of S/MARBPs is important for cellular functions and is altered under different conditions. Limited information is available presently to understand the structure–function relationship conclusively. Although all S/MARBPs bind to DNA, their context- and tissue-specific regulatory roles cannot be justified solely based on the available information on their structures. Conformational changes in a protein lead to changes in protein–protein interactions (PPIs) that essentially would regulate functional outcomes. A well-studied form of protein regulation is post-translational modification (PTM). It involves disulfide bond formation, cleavage of precursor proteins, and addition or removal of low-molecular-weight groups, leading to modifications like phosphorylation, methylation, SUMOylation, acetylation, PARylation, and ubiquitination. These chemical modifications lead to varied functional outcomes by mechanisms like modifying DNA–protein interactions and PPIs, altering protein function, stability, and crosstalk with other PTMs regulating subcellular localizations. S/MARBPs are reported to be regulated by PTMs, thereby contributing to gene regulation. In this review, we discuss the current understanding, scope, disease implications, and future perspectives of the diverse PTMs regulating functions of S/MARBPs.

The biochemical and genetic discovery of the SAGA complex

One of the major advances in our understanding of gene regulation in eukaryotes was the discovery of factors that regulate transcription by controlling chromatin structure. Prominent among these discoveries was the demonstration that Gcn5 is a histone acetyltransferase, establishing a direct connection between transcriptional activation and histone acetylation. This breakthrough was soon followed by the purification of a protein complex that contains Gcn5, the SAGA complex. In this article, we review the early genetic and biochemical experiments that led to the discovery of SAGA and the elucidation of its multiple activities.