Eukaryotic Gene Regulation Research Articles

From bacterial operons to gene therapy: 50 years of the journal Cell.

Recombinant DNA technology has profoundly advanced virtually every aspect of biological and medical sciences, from basic research to biotechnology. Here, I discuss conceptual connections linking fundamental discoveries that were enabled by the technology, advances in the understanding of gene regulation in both prokaryotes and eukaryotes, and the recent FDA-approved CRISPR-based gene therapy for sickle cell anemia and β-thalassemia based on transcriptional derepression.

A fine kinetic balance of interactions directs transcription factor hubs to genes.

Eukaryotic gene regulation relies on the binding of sequence-specific transcription factors (TFs). TFs bind chromatin transiently yet occupy their target sites by forming high-local concentration microenvironments (hubs and condensates) that increase the frequency of binding. Despite their ubiquity, such microenvironments are difficult to study in endogenous contexts due to technical limitations. Here, we use live embryo light-sheet imaging, single-molecule tracking, and genomics to overcome these limitations and investigate how hubs are localized to target genes to drive TF occupancy and transcription. By examining mutants of a hub-forming TF, Zelda, in Drosophila embryos, we find that hub formation propensity, spatial distributions, and temporal stabilities are differentially regulated by DNA binding and disordered protein domains. We show that hub localization to genomic targets is driven by a finely-tuned kinetic balance of interactions between proteins and chromatin, and hubs can be redirected to new genomic sites when this balance is perturbed.

Evolution of a Eukaryotic Transcription Factor's co-TF Dependence Involves Multiple Intrinsically Disordered Regions Affecting Activation and Autoinhibition.

Combinatorial control by multiple transcription factors (TFs) is a hallmark of eukaryotic gene regulation. Despite its prevalence and crucial roles in enhancing specificity and integrating information, the mechanisms behind why eukaryotic TFs depend on one another, and whether such interdependence evolves, are not well understood. We exploit natural variation in co-TF dependence in the yeast phosphate starvation (PHO) response to address this question. In the model yeast Saccharomyces cerevisiae, the main TF, Pho4, relies on the co-TF Pho2 to regulate ~28 genes. In a related yeast pathogen, Candida glabrata, its Pho4 exhibits significantly reduced Pho2 dependence and has an expanded target set of ~70 genes. Biochemical analyses showed C. glabrata Pho4 (CgPho4) binds to the same consensus motif with 3-4-fold higher affinity than ScPho4 does. A machine-learning-based prediction and yeast one-hybrid assay identified two Intrinsically Disordered Regions (IDRs) in CgPho4 that boost the activity of the main activation domain but showed little to no activity on their own. We also found evidence for autoinhibition behind the co-TF dependence in ScPho4. An IDR in ScPho4 next to its DNA binding domain was found to act as a double-edged sword: it both allows for enhanced activity with Pho2, and inhibits Pho4's activity without Pho2. This study provides a detailed molecular picture of how co-TF dependence is mediated and how its evolution, mainly driven by IDR divergence, can lead to significant rewiring of the regulatory network.

Actively Expressed Intergenic Genes Generated by Transposable Element Insertions in Gossypium hirsutum Cotton.

The genomes and annotated genes of allotetraploid cotton Gossypium hirsutum have been extensively studied in recent years. However, the expression, regulation, and evolution of intergenic genes (ITGs) have not been completely deciphered. In this study, we identified a novel set of actively expressed ITGs in G. hirsutum cotton, through transcriptome profiling based on deep sequencing data, as well as chromatin immunoprecipitation, followed by sequencing (ChIP-seq) of histone modifications and how the ITGs evolved. Totals of 17,567 and 8249 ITGs were identified in G. hirsutum and Gossypium arboreum, respectively. The expression of ITGs in G. hirsutum was significantly higher than that in G. arboreum. Moreover, longer exons were observed in G. hirsutum ITGs. Notably, 42.3% of the ITGs from G. hirsutum were generated by the long terminal repeat (LTR) insertions, while their proportion in genic genes was 19.9%. The H3K27ac and H3K4me3 modification proportions and intensities of ITGs were equivalent to genic genes. The H3K4me1 modifications were lower in ITGs. Additionally, evolution analyses revealed that the ITGs from G. hirsutum were mainly produced around 6.6 and 1.6 million years ago (Mya), later than the pegged time for genic genes, which is 7.0 Mya. The characterization of ITGs helps to elucidate the evolution of cotton genomes and shed more light on their biological functions in the transcriptional regulation of eukaryotic genes, along with the roles of histone modifications in speciation and diversification.

A single-cell multimodal view on gene regulatory network inference from transcriptomics and chromatin accessibility data.

Eukaryotic gene regulation is a combinatorial, dynamic, and quantitative process that plays a vital role in development and disease and can be modeled at a systems level in gene regulatory networks (GRNs). The wealth of multi-omics data measured on the same samples and even on the same cells has lifted the field of GRN inference to the next stage. Combinations of (single-cell) transcriptomics and chromatin accessibility allow the prediction of fine-grained regulatory programs that go beyond mere correlation of transcription factor and target gene expression, with enhancer GRNs (eGRNs) modeling molecular interactions between transcription factors, regulatory elements, and target genes. In this review, we highlight the key components for successful (e)GRN inference from (sc)RNA-seq and (sc)ATAC-seq data exemplified by state-of-the-art methods as well as open challenges and future developments. Moreover, we address preprocessing strategies, metacell generation and computational omics pairing, transcription factor binding site detection, and linear and three-dimensional approaches to identify chromatin interactions as well as dynamic and causal eGRN inference. We believe that the integration of transcriptomics together with epigenomics data at a single-cell level is the new standard for mechanistic network inference, and that it can be further advanced with integrating additional omics layers and spatiotemporal data, as well as with shifting the focus towards more quantitative and causal modeling strategies.

Diverse Gene Regulatory Mechanisms Alter Rattlesnake Venom Gene Expression at Fine Evolutionary Scales.

Understanding and predicting the relationships between genotype and phenotype is often challenging, largely due to the complex nature of eukaryotic gene regulation. A step towards this goal is to map how phenotypic diversity evolves through genomic changes that modify gene regulatory interactions. Using the Prairie Rattlesnake (Crotalus viridis) and related species, we integrate mRNA-seq, proteomic, ATAC-seq and whole-genome resequencing data to understand how specific evolutionary modifications to gene regulatory network components produce differences in venom gene expression. Through comparisons within and between species, we find a remarkably high degree of gene expression and regulatory network variation across even a shallow level of evolutionary divergence. We use these data to test hypotheses about the roles of specific trans-factors and cis-regulatory elements, how these roles may vary across venom genes and gene families, and how variation in regulatory systems drive diversity in venom phenotypes. Our results illustrate that differences in chromatin and genotype at regulatory elements play major roles in modulating expression. However, we also find that enhancer deletions, differences in transcription factor expression, and variation in activity of the insulator protein CTCF also likely impact venom phenotypes. Our findings provide insight into the diversity and gene-specificity of gene regulatory features and highlight the value of comparative studies to link gene regulatory network variation to phenotypic variation.

Co-transcriptional gene regulationin eukaryotes and prokaryotes.

Many steps of RNA processing occur during transcription by RNA polymerases. Co-transcriptional activities are deemed commonplace in prokaryotes, in which the lack of membrane barriers allows mixing of all gene expression steps, from transcription to translation. In the past decade, an extraordinary level of coordination between transcription and RNA processing has emerged in eukaryotes. In this Review, we discuss recent developments in our understanding of co-transcriptional gene regulation in both eukaryotes and prokaryotes, comparing methodologies and mechanisms, and highlight striking parallels in how RNA polymerases interact with the machineries that act on nascent RNA. The development of RNA sequencing and imaging techniques that detect transient transcription and RNA processing intermediates has facilitated discoveries of transcription coordination with splicing, 3'-end cleavage and dynamic RNA folding and revealed physical contacts between processing machineries and RNA polymerases. Such studies indicate that intron retention in a given nascent transcript can prevent 3'-end cleavage and cause transcriptional readthrough, which is a hallmark of eukaryotic cellular stress responses. We also discuss how coordination between nascentRNA biogenesis and transcription drives fundamental aspects of gene expression in both prokaryotes and eukaryotes.

Messenger RNA Surveillance: Current Understanding, Regulatory Mechanisms, and Future Implications.

Nonsense-mediated mRNA decay (NMD) is an evolutionarily conserved surveillance mechanism in eukaryotes primarily deployed to ensure RNA quality control by eliminating aberrant transcripts and also involved in modulating the expression of several physiological transcripts. NMD, the mRNA surveillance pathway, is a major form of gene regulation in eukaryotes. NMD serves as one of the most significant quality control mechanisms as it primarily scans the newly synthesized transcripts and differentiates the aberrant and non-aberrant transcripts. The synthesis of truncated proteins is restricted, which would otherwise lead to cellular dysfunctions. The up-frameshift factors (UPFs) play a central role in executing the NMD event, largely by recognizing and recruiting multiple protein factors that result in the decay of non-physiological mRNAs. NMD exhibits astounding variability in its ability across eukaryotes in an array of pathological and physiological contexts. The detailed understanding of NMD and the underlying molecular mechanisms remains blurred. This review outlines our current understanding of NMD, in regulating multifaceted cellular events during development and disease. It also attempts to identify unanswered questions that deserve further investigation.

Molecular basis of chromatin remodelling by DDM1 involved in plant DNA methylation.

Eukaryotic gene regulation occurs at the chromatin level, which requires changing the chromatin structure by a group of ATP-dependent DNA translocases-namely, the chromatin remodellers1. In plants, chromatin remodellers function in various biological processes and possess both conserved and plant-specific components2-5. DECREASE IN DNA METHYLATION 1 (DDM1) is a plant chromatin remodeller that plays a key role in the maintenance DNA methylation6-11. Here we determined the structures of Arabidopsis DDM1 in complex with nucleosome in ADP-BeFx-bound, ADP-bound and nucleotide-free conformations. We show that DDM1 specifically recognizes the H4 tail and nucleosomal DNA. The conformational differences between ADP-BeFx-bound, ADP-bound and nucleotide-free DDM1 suggest a chromatin remodelling cycle coupled to ATP binding, hydrolysis and ADP release. This, in turn, triggers conformational changes in the DDM1-bound nucleosomal DNA, which alters the nucleosome structure and promotes DNA sliding. Together, our data reveal the molecular basis of chromatin remodelling by DDM1.

The Giardia lamblia ribosome structure reveals divergence in several biological pathways and the mode of emetine function

Giardia lamblia is a deeply branching protist and a human pathogen. Its unusual biology presents the opportunity to explore conserved and fundamental molecular mechanisms. We determined the structure of the G.lamblia 80S ribosome bound to tRNA, mRNA, and the antibiotic emetine by cryo-electron microscopy, to an overall resolution of 2.49Å. The structure reveals rapidly evolving protein and nucleotide regions, differences in the peptide exit tunnel, and likely altered ribosome quality control pathways. Examination of translation initiation factor binding sites suggests these interactions are conserved despite a divergent initiation mechanism. Highlighting the potential of G. lamblia to resolve conserved biological principles; our structure reveals the interactions of the translation inhibitor emetine with the ribosome and mRNA, thus providing insight into the mechanism of action for this widely used antibiotic. Our work defines key questions in G. lamblia and motivates future experiments to explore the diversity of eukaryotic gene regulation.

A covalent compound selectively inhibits RNA demethylase ALKBH5 rather than FTO.

N 6-Methyladenosine (m6A) is the most prevalent mRNA modification and is required for gene regulation in eukaryotes. ALKBH5, an m6A demethylase, is a promising target, particularly for anticancer drug discovery. However, the development of selective and potent inhibitors of ALKBH5 rather than FTO remains challenging. Herein, we used a targeted covalent inhibition strategy and identified a covalent inhibitor, TD19, which selectively inhibits ALKBH5 compared with FTO demethylase in protein-based and tumor cell-based assays. TD19 irreversibly modifies the residues C100 and C267, preventing ALKBH5 from binding to m6A-containing RNA. Moreover, TD19 displays good anticancer efficacy in acute myeloid leukemia and glioblastoma multiforme cell lines. Thus, the ALKBH5 inhibitor developed in this study, which selectively targets ALKBH5 compared with FTO, can potentially be used as a probe for investigating the biological functions of RNA demethylase and as a lead compound in anticancer research.

Development of a rapid mass spectrometric method for the analysis of ten-eleven translocation enzymes.

In DNA, methylation at the fifth position of cytosine (5mC) by DNA methyltransferases is essential for eukaryotic gene regulation. Methylation patterns are dynamically controlled by epigenetic machinery. Erasure of 5mC by Fe2+ and 2-ketoglutarate (2KG) dependent dioxygenases in the ten-eleven translocation family (TET1-3), plays a key role in nuclear processes. Through the event of active demethylation, TET proteins iteratively oxidize 5mC to 5-hydroxymethyl cytosine (5hmC), 5-formylcytosine (5fC) and 5-carboxycytosine (5caC), each of which has been implicated in numerous diseases when aberrantly generated. A wide range of biochemical assays have been developed to characterize TET activity, many of which require multi-step processing to detect and quantify the 5mC oxidized products. Herein, we describe the development and optimization of a sensitive MALDI mass spectrometry-based technique that directly measures TET activity and eliminates tedious processing steps. Employing optimized assay conditions, we report the steady-state activity of wild type TET2 enzymes to furnish 5hmC, 5fC and 5caC. We next determine IC50 values of several small-molecule inhibitors of TETs. The utility of this assay is further demonstrated by analyzing the activity of V1395A which is an activating mutant of TET2 that primarily generates 5caC. Lastly, we describe the development of a secondary assay that utilizes bisulfite chemistry to further examine the activity of wildtype TET2 and V1395A in a base-resolution manner. The combined results demonstrate that the activity of TET proteins can be gauged, and their products accurately quantified using our methods.

P300 is an obligate integrator of combinatorial transcription factor inputs

Transcription coactivators are proteins or protein complexes that mediate transcription factor (TF) function. However, they lack DNA-binding capacity, prompting the question of how they engage target loci. Three non-exclusive hypotheses have been posited: coactivators are recruited by complexing with TFs, by binding histones through epigenetic reader domains, or by partitioning into condensates through their extensive intrinsically disordered regions. Using p300 as a prototypical coactivator, we systematically mutated its annotated domains and show by single-molecule tracking in live U2OS cells that coactivator-chromatin binding depends entirely on combinatorial binding of multiple TF-interaction domains. Furthermore, we demonstrate that acetyltransferase activity opposes p300-chromatin association and that the N-terminal TF-interaction domains regulate that activity. Single TF-interaction domains are insufficient for chromatin binding and regulation of catalytic activity, implying a principle that we speculate could broadly apply to eukaryotic gene regulation: a TF must act in coordination with other TFs to recruit coactivator activity.

A YTHDF-PABP interaction is required for m6 A-mediated organogenesis in plants.

N6-methyladenosine (m6 A) in mRNA is key to eukaryotic gene regulation. Many m6 A functions involve RNA-binding proteins that recognize m6 A via a YT521-B Homology (YTH) domain. YTH domain proteins contain long intrinsically disordered regions (IDRs) that may mediate phase separation and interaction with protein partners, but whose precise biochemical functions remain largely unknown. The Arabidopsis thaliana YTH domain proteins ECT2, ECT3, and ECT4 accelerate organogenesis through stimulation of cell division in organ primordia. Here, we use ECT2 to reveal molecular underpinnings of this function. We show that stimulation of leaf formation requires the long N-terminal IDR, and we identify two short IDR elements required for ECT2-mediated organogenesis. Of these two, a 19-amino acid region containing a tyrosine-rich motif conserved in both plant and metazoan YTHDF proteins is necessary for binding to the major cytoplasmic poly(A)-binding proteins PAB2, PAB4, and PAB8. Remarkably, overexpression of PAB4 in leaf primordia partially rescues the delayed leaf formation in ect2 ect3 ect4 mutants, suggesting that the ECT2-PAB2/4/8 interaction on target mRNAs of organogenesis-related genes may overcome limiting PAB concentrations in primordial cells.

Identification and Profiling Analysis of microRNAs in Guava Fruit (Psidium guajava L.) and Their Role during Ripening.

The guava (Psidium guajava L.) is a climacteric fruit with an accelerated post-harvest overripening. miRNAs are small RNA sequences that function as gene regulators in eukaryotes and are essential for their survival and development. In this study, miRNA libraries were constructed, sequenced and analyzed from the breaker and ripe stages of guava fruit cv. Siglo XXI. One hundred and seventy-four mature miRNA sequences from 28 miRNA families were identified. The taxonomic distribution of the guava miRNAs showed a high level of conservation among the dicotyledonous plants. Most of the predicted miRNA target genes were transcription factors and genes involved in the metabolism of phytohormones such as abscisic acid, auxins, and ethylene, as revealed through an ontology enrichment analysis. The miRNA families miR168, miR169, miR396, miR397, and miR482 were classified as being directly associated with maturation, whereas the miRNA families miR160, miR165, miR167, miR3930, miR395, miR398, and miR535 were classified as being indirectly associated. With this study, we intended to increase our knowledge and understanding of the regulatory process involved in the ripening process, thereby providing valuable information for future research on the ripening of guava fruit.

High-throughput data and modeling reveal insights into the mechanisms of cooperative DNA-binding by transcription factor proteins.

Cooperative DNA-binding by transcription factor (TF) proteins is critical for eukaryotic gene regulation. In the human genome, many regulatory regions contain TF-binding sites in close proximity to each other, which can facilitate cooperative interactions. However, binding site proximity does not necessarily imply cooperative binding, as TFs can also bind independently to each of their neighboring target sites. Currently, the rules that drive cooperative TF binding are not well understood. In addition, it is oftentimes difficult to infer direct TF-TF cooperativity from existing DNA-binding data. Here, we show that in vitro binding assays using DNA libraries of a few thousand genomic sequences with putative cooperative TF-binding events can be used to develop accurate models of cooperativity and to gain insights into cooperative binding mechanisms. Using factors ETS1 and RUNX1 as our case study, we show that the distance and orientation between ETS1 sites are critical determinants of cooperative ETS1-ETS1 binding, while cooperative ETS1-RUNX1 interactions show more flexibility in distance and orientation and can be accurately predicted based on the affinity and sequence/shape features of the binding sites. The approach described here, combining custom experimental design with machine-learning modeling, can be easily applied to study the cooperative DNA-binding patterns of any TFs.

Reading and erasing of histone crotonyllysine mimics by the AF9 YEATS domain and SIRT2 deacylase

Lysine acylations on histones and their recognition by chromatin-binding reader domains and removal by histone deacylases function as an important mechanism for eukaryotic gene regulation. Histone lysine crotonylation (Kcr) is an epigenetic mark associated with active transcription, and its installation and removal are dynamically regulated by cellular epigenetic enzymes. Here, we report binding studies and enzyme assays with histone H3K9 peptides bearing simplest Kcr analogs with varying hydrocarbon chain length, bulkiness, rigidity and polarity. We demonstrate that the AF9 YEATS domain displays selectivity for binding of different acylation modifications on histone H3K9 peptides and exhibits preference for bulkier cinnamoylated lysine over crotonylated lysine and its mimics. SIRT2 shows deacylase activity against most of acylated H3K9 peptides bearing different crotonyllysine mimics, however, it displays a poor ability for the removal of cinnamoyl and trifluorocrotonyl groups. These results demonstrate different substrate selectivities of epigenetic proteins acting on crotonyllysine and pave the way for rational design and development of AF9 YEATS and SIRT2 inhibitors for treatment of human diseases, including cancer.

Role of RNA Alternative Splicing in T Cell Function and Disease.

Alternative RNA splicing, a ubiquitous mechanism of gene regulation in eukaryotes, expands genome coding capacity and proteomic diversity. It has essential roles in all aspects of human physiology, including immunity. This review highlights the importance of RNA alternative splicing in regulating immune T cell function. We discuss how mutations that affect the alternative splicing of T cell factors can contribute to abnormal T cell function and ultimately lead to autoimmune diseases. We also explore the potential applications of strategies that target the alternative splicing changes of T cell factors. These strategies could help design therapeutic approaches to treat autoimmune disorders and improve immunotherapy.

Driving forces behind phase separation of the carboxy-terminal domain of RNA polymerase II

Eukaryotic gene regulation and pre-mRNA transcription depend on the carboxy-terminal domain (CTD) of RNA polymerase (Pol) II. Due to its highly repetitive, intrinsically disordered sequence, the CTD enables clustering and phase separation of Pol II. The molecular interactions that drive CTD phase separation and Pol II clustering are unclear. Here, we show that multivalent interactions involving tyrosine impart temperature- and concentration-dependent self-coacervation of the CTD. NMR spectroscopy, molecular ensemble calculations and all-atom molecular dynamics simulations demonstrate the presence of diverse tyrosine-engaging interactions, including tyrosine-proline contacts, in condensed states of human CTD and other low-complexity proteins. We further show that the network of multivalent interactions involving tyrosine is responsible for the co-recruitment of the human Mediator complex and CTD during phase separation. Our work advances the understanding of the driving forces of CTD phase separation and thus provides the basis to better understand CTD-mediated Pol II clustering in eukaryotic gene transcription.

Single molecule analysis of CENP-A chromatin by high-speed atomic force microscopy.

Chromatin accessibility is modulated in a variety of ways to create open and closed chromatin states, both of which are critical for eukaryotic gene regulation. At the single molecule level, how accessibility is regulated of the chromatin fiber composed of canonical or variant nucleosomes is a fundamental question in the field. Here, we developed a single-molecule tracking method where we could analyze thousands of canonical H3 and centromeric variant nucleosomes imaged by high-speed atomic force microscopy. This approach allowed us to investigate how changes in nucleosome dynamics in vitro inform us about transcriptional potential in vivo. By high-speed atomic force microscopy, we tracked chromatin dynamics in real time and determined the mean square displacement and diffusion constant for the variant centromeric CENP-A nucleosome. Furthermore, we found that an essential kinetochore protein CENP-C reduces the diffusion constant and mobility of centromeric nucleosomes along the chromatin fiber. We subsequently interrogated how CENP-C modulates CENP-A chromatin dynamics in vivo. Overexpressing CENP-C resulted in reduced centromeric transcription and impaired loading of new CENP-A molecules. From these data, we speculate that factors altering nucleosome mobility in vitro, also correspondingly alter transcription in vivo. Subsequently, we propose a model in which variant nucleosomes encode their own diffusion kinetics and mobility, and where binding partners can suppress or enhance nucleosome mobility.