Resolving noise-control conflict by gene duplication

Trans-eQTL hotspots shape complex traits by modulating cellular states

Regulatory genetic variation shapes gene expression, providing an important mechanism connecting DNA variation and complex traits. The causal relationships between gene expression and complex traits remain poorly understood. Here, we integrated transcriptomes and 46 genetically complex growth traits in a large cross between two strains of the yeast Saccharomyces cerevisiae. We discovered thousands of genetic correlations between gene expression and growth, suggesting potential functional connections. Local regulatory variation was a minor source of these genetic correlations. Instead, genetic correlations tended to arise from multiple independent trans-acting regulatory loci. Trans-acting hotspots that affect the expression of numerous genes accounted for particularly large fractions of genetic growth variation and of genetic correlations between gene expression and growth. Genes with genetic correlations were enriched for similar biological processes across traits but with heterogeneous direction of effect. Our results reveal how trans-acting regulatory hotspots shape complex traits by altering cellular states.

Plasticity and environment-specific relationships between gene expression and fitness in Saccharomyces cerevisiae

The environment influences how an organism's genotype determines its phenotype and how this phenotype affects its fitness. Here, to better understand this dual role of environment in the production and selection of phenotypic variation, we determined genotype-phenotype-fitness relationships for mutant strains of Saccharomyces cerevisiae in four environments. Specifically, we measured how promoter mutations of the metabolic gene TDH3 modified expression level and affected growth for four different carbon sources. In each environment, we observed a clear relationship between TDH3 expression level and fitness, but this relationship differed among environments. Mutations with similar effects on expression in different environments often had different effects on fitness and vice versa. Such environment-specific relationships between phenotype and fitness can shape the evolution of phenotypic plasticity. We also found that mutations disrupting binding sites for transcription factors had more variable effects on expression among environments than those disrupting the TATA box, which is part of the core promoter. However, mutations with the most environmentally variable effects on fitness were located in the TATA box, because of both the lack of plasticity of TATA box mutations and environment-specific fitness functions. This observation suggests that mutations affecting different molecular mechanisms contribute unequally to regulatory sequence evolution in changing environments.

Beyond RNA-binding domains: determinants of protein-RNA binding

RNA-binding proteins (RBPs) are composed of RNA-binding domains (RBDs) often linked via intrinsically disordered regions (IDRs). Structural and biochemical analyses have shown that disordered linkers contribute to RNA binding by orienting the adjacent RBDs and also characterized certain disordered repeats that directly contact the RNA. However, the relative contribution of IDRs and predicted RBDs to the in vivo binding pattern is poorly explored. Here, we upscaled the RNA-tagging method to map the transcriptome-wide binding of 16 RBPs in budding yeast. We then performed extensive sequence mutations to distinguish binding determinants within predicted RBDs and the surrounding IDRs in eight of these. The majority of the predicted RBDs tested were not individually essential for mRNA binding. However, multiple IDRs that lacked predicted RNA-binding potential appeared essential for binding affinity or specificity. Our results provide new insights into the function of poorly studied RBPs and emphasize the complex and distributed encoding of RBP-RNA interaction in vivo.

Revisiting the model for coactivator recruitment: Med15 can select its target sites independent of promoter-bound transcription factors

Activation domains (ADs) within transcription factors (TFs) induce gene expression by recruiting coactivators such as the Mediator complex. Coactivators lack DNA binding domains (DBDs) and are assumed to passively follow their recruiting TFs. This is supported by direct AD-coactivator interactions seen in vitro but has not yet been tested in living cells. To examine that, we targeted two Med15-recruiting ADs to a range of budding yeast promoters through fusion with different DBDs. The DBD-AD fusions localized to hundreds of genomic sites but recruited Med15 and induced transcription in only a subset of bound promoters, characterized by a fuzzy-nucleosome architecture. Direct DBD-Med15 fusions shifted DBD localization towards fuzzy-nucleosome promoters, including promoters devoid of the endogenous Mediator. We propose that Med15, and perhaps other coactivators, possess inherent promoter preference and thus actively contribute to the selection of TF-induced genes.

Gene networks governing the response of a calcareous sponge to future ocean conditions reveal lineage-specific XBP1 regulation of the unfolded protein response

Marine sponges are predicted to be winners in the future ocean due to their exemplary adaptive capacity. However, while many sponge groups exhibit tolerance to a wide range of environmental insults, calcifying sponges may be more susceptible to thermo-acidic stress. To describe the gene regulatory networks that govern the stress response of the calcareous sponge, Leucetta chagosensis (class Calcarea, order Clathrinida), individuals were subjected to warming and acidification conditions based on the climate models for 2100. Transcriptome analysis and gene co-expression network reconstruction revealed that the unfolded protein response (UPR) was activated under thermo-acidic stress. Among the upregulated genes were two lineage-specific homologs of X-box binding protein 1 (XBP1), a transcription factor that activates the UPR. Alternative dimerization between these XBP1 gene products suggests a clathrinid-specific mechanism to reversibly sequester the transcription factor into an inactive form, enabling the rapid regulation of pathways linked to the UPR in clathrinid calcareous sponges. Our findings support the idea that transcription factor duplication events may refine evolutionarily conserved molecular pathways and contribute to ecological success.

Ordered and disordered regions of the Origin Recognition Complex direct differential in vivo binding at distinct motif sequences

The Origin Recognition Complex (ORC) seeds replication-fork formation by binding to DNA replication origins, which in budding yeast contain a 17bp DNA motif. High resolution structure of the ORC-DNA complex revealed two base-interacting elements: a disordered basic patch (Orc1-BP4) and an insertion helix (Orc4-IH). To define the ORC elements guiding its DNA binding in vivo, we mapped genomic locations of 38 designed ORC mutants, revealing that different ORC elements guide binding at different sites. At silencing-associated sites lacking the motif, ORC binding and activity were fully explained by a BAH domain. Within replication origins, we reveal two dominating motif variants showing differential binding modes and symmetry: a non-repetitive motif whose binding requires Orc1-BP4 and Orc4-IH, and a repetitive one where another basic patch, Orc1-BP3, can replace Orc4-IH. Disordered basic patches are therefore key for ORC-motif binding in vivo, and we discuss how these conserved, minor-groove interacting elements can guide specific ORC-DNA recognition.

Intrinsically disordered regions of the Msn2 transcription factor encode multiple functions using interwoven sequence grammars

Intrinsically disordered regions (IDRs) are abundant in eukaryotic proteins, but their sequence-function relationship remains poorly understood. IDRs of transcription factors (TFs) can direct promoter selection and recruit coactivators, as shown for the budding yeast TF Msn2. To examine how IDRs encode both these functions, we compared genomic binding specificity, coactivator recruitment, and gene induction amongst a large set of designed Msn2-IDR mutants. We find that both functions depend on multiple regions across the > 600AA IDR. Yet, transcription activity was readily disrupted by mutations that showed no effect on the Msn2 binding specificity. Our data attribute this differential sensitivity to the integration of a relaxed, composition-based code directing binding specificity with a more stringent, motif-based code controlling the recruitment of coactivators and transcription activity. Therefore, Msn2 utilizes interwoven sequence grammars for encoding multiple functions, suggesting a new IDR design paradigm of potentially general use.

Divergence of TORC1-mediated stress response leads to novel acquired stress resistance in a pathogenic yeast

Acquired stress resistance (ASR) enables organisms to prepare for environmental changes that occur after an initial stressor. However, the genetic basis for ASR and how the underlying network evolved remain poorly understood. In this study, we discovered that a short phosphate starvation induces oxidative stress response (OSR) genes in the pathogenic yeast C. glabrata and protects it against a severe H2O2 stress; the same treatment, however, provides little benefit in the low pathogenic-potential relative, S. cerevisiae. This ASR involves the same transcription factors (TFs) as the OSR, but with different combinatorial logics. We show that Target-of-Rapamycin Complex 1 (TORC1) is differentially inhibited by phosphate starvation in the two species and contributes to the ASR via its proximal effector, Sch9. Therefore, evolution of the phosphate starvation-induced ASR involves the rewiring of TORC1's response to phosphate limitation and the repurposing of TF-target gene networks for the OSR using new regulatory logics.

The architecture of binding cooperativity between densely bound transcription factors

The binding of transcription factors (TFs) along genomes is restricted to a subset of sites containing their preferred motifs. TF-binding specificity is often attributed to the co-binding of interacting TFs; however, apart from specific examples, this model remains untested. Here, we define dependencies among budding yeast TFs that localize to overlapping promoters by profiling the genome-wide consequences of co-depleting multiple TFs. We describe unidirectional interactions, revealing Msn2 as a central factor allowing TF binding at its target promoters. By contrast, no case of mutual cooperation was observed. Particularly, Msn2 retained binding at its preferred promoters upon co-depletion of fourteen similarly bound TFs. Overall, the consequences of TF co-depletions were moderate, limited to a subset of promoters, and failed to explain the role of regions outside the DNA-binding domain in directing TF-binding preferences. Our results call for re-evaluating the role of cooperative interactions in directing TF-binding preferences.

Divergence of TORC1-mediated Stress Response Leads to Novel Acquired Stress Resistance in a Pathogenic Yeast

Acquired stress resistance (ASR) enables organisms to prepare for environmental changes that occur after an initial stressor. However, the genetic basis for ASR and how the underlying network evolved remain poorly understood. In this study, we discovered that a short phosphate starvation induces oxidative stress response (OSR) genes in the pathogenic yeast C. glabrata and protects it against a severe H2O2 stress; the same treatment, however, provides little benefit in the low pathogenic-potential relative, S. cerevisiae. This ASR involves the same transcription factors (TFs) as the OSR, but with different combinatorial logics. We show that Target-of-Rapamycin Complex 1 (TORC1) is differentially inhibited by phosphate starvation in the two species and contributes to the ASR via its proximal effector, Sch9. Therefore, evolution of the phosphate starvation-induced ASR involves the rewiring of TORC1's response to phosphate limitation and the repurposing of TF-target gene networks for the OSR using new regulatory logics.

The yeast RNA methylation complex consists of conserved yet reconfigured components with m6A-dependent and independent roles

N6-methyladenosine (m6A), the most abundant mRNA modification, is deposited in mammals/insects/plants by m6A methyltransferase complexes (MTC) comprising a catalytic subunit and at least five additional proteins. The yeast MTC is critical for meiosis and was known to comprise three proteins, of which two were conserved. We uncover three novel MTC components (Kar4/Ygl036w-Vir1/Dyn2). All MTC subunits, except for Dyn2, are essential for m6A deposition and have corresponding mammalian MTC orthologues. Unlike the mammalian bipartite MTC, the yeast MTC is unipartite, yet multifunctional. The mRNA interacting module, comprising Ime4, Mum2, Vir1, and Kar4, exerts the MTC's m6A-independent function, while Slz1 enables the MTC catalytic function in m6A deposition. Both functions are critical for meiotic progression. Kar4 also has a mechanistically separate role from the MTC during mating. The yeast MTC constituents play distinguishable m6A-dependent, MTC-dependent, and MTC-independent functions, highlighting their complexity and paving the path towards dissecting multi-layered MTC functions in mammals.

Paralog-based synthetic lethality: rationales and applications

Tumor cells can result from gene mutations and over-expression. Synthetic lethality (SL) offers a desirable setting where cancer cells bearing one mutated gene of an SL gene pair can be specifically targeted by disrupting the function of the other genes, while leaving wide-type normal cells unharmed. Paralogs, a set of homologous genes that have diverged from each other as a consequence of gene duplication, make the concept of SL feasible as the loss of one gene does not affect the cell's survival. Furthermore, homozygous loss of paralogs in tumor cells is more frequent than singletons, making them ideal SL targets. Although high-throughput CRISPR-Cas9 screenings have uncovered numerous paralog-based SL pairs, the unclear mechanisms of targeting these gene pairs and the difficulty in finding specific inhibitors that exclusively target a single but not both paralogs hinder further clinical development. Here, we review the potential mechanisms of paralog-based SL given their function and genetic combination, and discuss the challenge and application prospects of paralog-based SL in cancer therapeutic discovery.

Evolutionary trade-off and mutational bias could favor transcriptional over translational divergence within paralog pairs

How changes in the different steps of protein synthesis-transcription, translation and degradation-contribute to differences of protein abundance among genes is not fully understood. There is however accumulating evidence that transcriptional divergence might have a prominent role. Here, we show that yeast paralogous genes are more divergent in transcription than in translation. We explore two causal mechanisms for this predominance of transcriptional divergence: an evolutionary trade-off between the precision and economy of gene expression and a larger mutational target size for transcription. Performing simulations within a minimal model of post-duplication evolution, we find that both mechanisms are consistent with the observed divergence patterns. We also investigate how additional properties of the effects of mutations on gene expression, such as their asymmetry and correlation across levels of regulation, can shape the evolution of paralogs. Our results highlight the importance of fully characterizing the distributions of mutational effects on transcription and translation. They also show how general trade-offs in cellular processes and mutation bias can have far-reaching evolutionary impacts.

Molecular evolution of the ependymin-related gene epdl2 in African weakly electric fish

Gene duplication and subsequent molecular evolution can give rise to taxon-specific gene specializations. In previous work, we found evidence that African weakly electric fish (Mormyridae) may have as many as three copies of the epdl2 gene, and the expression of two epdl2 genes is correlated with electric signal divergence. Epdl2 belongs to the ependymin-related family (EPDR), a functionally diverse family of secretory glycoproteins. In this study, we first describe vertebrate EPDR evolution and then present a detailed evolutionary history of epdl2 in Mormyridae with emphasis on the speciose genus Paramormyrops. Using Sanger sequencing, we confirm three apparently functional epdl2 genes in Paramormyrops kingsleyae. Next, we developed a nanopore-based amplicon sequencing strategy and bioinformatics pipeline to obtain and classify full-length epdl2 gene sequences (N = 34) across Mormyridae. Our phylogenetic analysis proposes three or four epdl2 paralogs dating from early Paramormyrops evolution. Finally, we conducted selection tests which detected positive selection around the duplication events and identified ten sites likely targeted by selection in the resulting paralogs. These sites' locations in our modeled 3D protein structure involve four sites in ligand binding and six sites in homodimer formation. Together, these findings strongly imply an evolutionary mechanism whereby epdl2 genes underwent selection-driven functional specialization after tandem duplications in the rapidly speciating Paramormyrops. Considering previous evidence, we propose that epdl2 may contribute to electric signal diversification in mormyrids, an important aspect of species recognition during mating.

Heterosis of fitness and phenotypic variance in the evolution of a diploid gene regulatory network

Heterosis describes the phenomenon, whereby a hybrid population has higher fitness than an inbred population, which has previously been explained by either Mendelian dominance or overdominance under the general assumption of a simple genotype-phenotype relationship. However, recent studies have demonstrated that genes interact through a complex gene regulatory network (GRN). Furthermore, phenotypic variance is reportedly lower for heterozygotes, and the origin of such variance-related heterosis remains elusive. Therefore, a theoretical analysis linking heterosis to GRN evolution and stochastic gene expression dynamics is required. Here, we investigated heterosis related to fitness and phenotypic variance in a system with interacting genes by numerically evolving diploid GRNs. According to the results, the heterozygote population exhibited higher fitness than the homozygote population, indicating fitness-related heterosis resulting from evolution. In addition, the heterozygote population exhibited lower noise-related phenotypic variance in expression levels than the homozygous population, implying that the heterozygote population is more robust to noise. Furthermore, the distribution of the ratio of heterozygote phenotypic variance to homozygote phenotypic variance exhibited quantitative similarity with previous experimental results. By applying dominance and differential gene expression rather than only a single gene expression model, we confirmed the correlation between heterosis and differential gene expression. We explain our results by proposing that the convex high-fitness region is evolutionarily shaped in the genetic space to gain noise robustness under genetic mixing through sexual reproduction. These results provide new insights into the effects of GRNs on variance-related heterosis and differential gene expression.

Intron-mediated induction of phenotypic heterogeneity

Intragenic regions that are removed during maturation of the RNA transcript-introns-are universally present in the nuclear genomes of eukaryotes1. The budding yeast, an otherwise intron-poor species, preserves two sets of ribosomal protein genes that differ primarily in their introns2,3. Although studies have shed light on the role of ribosomal protein introns under stress and starvation4-6, understanding the contribution of introns to ribosome regulation remains challenging. Here, by combining isogrowth profiling7 with single-cell protein measurements8, we show that introns can mediate inducible phenotypic heterogeneity that confers a clear fitness advantage. Osmotic stress leads to bimodal expression of the small ribosomal subunit protein Rps22B, which is mediated by an intron in the 5' untranslated region of its transcript. The two resulting yeast subpopulations differ in their ability to cope with starvation. Low levels of Rps22B protein result in prolonged survival under sustained starvation, whereas high levels of Rps22B enable cells to grow faster after transient starvation. Furthermore, yeasts growing at high concentrations of sugar, similar to those in ripe grapes, exhibit bimodal expression of Rps22B when approaching the stationary phase. Differential intron-mediated regulation of ribosomal protein genes thus provides a way to diversify the population when starvation threatens in natural environments. Our findings reveal a role for introns in inducing phenotypic heterogeneity in changing environments, and suggest that duplicated ribosomal protein genes in yeast contribute to resolving the evolutionary conflict between precise expression control and environmental responsiveness9.

Evolution of binding preferences among whole-genome duplicated transcription factors

Throughout evolution, new transcription factors (TFs) emerge by gene duplication, promoting growth and rewiring of transcriptional networks. How TF duplicates diverge was studied in a few cases only. To provide a genome-scale view, we considered the set of budding yeast TFs classified as whole-genome duplication (WGD)-retained paralogs (~35% of all specific TFs). Using high-resolution profiling, we find that ~60% of paralogs evolved differential binding preferences. We show that this divergence results primarily from variations outside the DNA-binding domains (DBDs), while DBD preferences remain largely conserved. Analysis of non-WGD orthologs revealed uneven splitting of ancestral preferences between duplicates, and the preferential acquiring of new targets by the least conserved paralog (biased neo/sub-functionalization). Interactions between paralogs were rare, and, when present, occurred through weak competition for DNA-binding or dependency between dimer-forming paralogs. We discuss the implications of our findings for the evolutionary design of transcriptional networks.

Coenzyme B12 -dependent and independent photoregulation of carotenogenesis across Myxococcales

Light-induced carotenogenesis in Myxococcus xanthus is controlled by the B₁₂ -based CarH repressor and photoreceptor, and by a separate intricate pathway involving singlet oxygen, the B₁₂ -independent CarH paralog CarA and various other proteins, some eukaryotic-like. Whether other myxobacteria conserve these pathways and undergo photoregulated carotenogenesis is unknown. Here, comparative analyses across 27 Myxococcales genomes identified carotenogenic genes, albeit arranged differently, with carH often in their genomic vicinity, in all three Myxococcales suborders. However, CarA and its associated factors were found exclusively in suborder Cystobacterineae, with carA-carH invariably in tandem in a syntenic carotenogenic operon, except for Cystobacter/Melittangium, which lack CarA but retain all other factors. We experimentally show B₁₂ -mediated photoregulated carotenogenesis in representative myxobacteria, and a remarkably plastic CarH operator design and DNA binding across Myxococcales. Unlike the two characterized CarH from other phyla, which are tetrameric, Cystobacter CarH (the first myxobacterial homolog amenable to analysis in vitro) is a dimer that combines direct CarH-like B₁₂ -based photoregulation with CarA-like DNA-binding and inhibition by an antirepressor. This study provides new molecular insights into B₁₂ -dependent photoreceptors. It further establishes the B₁₂ -dependent pathway for photoregulated carotenogenesis as broadly prevalent across myxobacteria and its evolution, exclusively in one suborder, into a parallel complex B₁₂ -independent circuit. This article is protected by copyright. All rights reserved.

Yeast Crf1p: An activator in need is an activator indeed

Ribosome biogenesis is an energetically costly process, and tight regulation is required for stoichiometric balance between components. This requires coordination of RNA polymerases I, II, and III. Lack of nutrients or the presence of stress leads to downregulation of ribosome biogenesis, a process for which mechanistic target of rapamycin complex I (mTORC1) is key. mTORC1 activity is communicated by means of specific transcription factors, and in yeast, which is a primary model system in which transcriptional coordination has been delineated, transcription factors involved in regulation of ribosomal protein genes include Fhl1p and its cofactors, Ifh1p and Crf1p. Ifh1p is an activator, whereas Crf1p has been implicated in maintaining the repressed state upon mTORC1 inhibition. Computational analyses of evolutionary relationships have indicated that Ifh1p and Crf1p descend from a common ancestor. Here, we discuss recent evidence, which suggests that Crf1p also functions as an activator. We propose a model that consolidates available experimental evidence, which posits that Crf1p functions as an alternate activator to prevent the stronger activator Ifh1p from re-binding gene promoters upon mTORC1 inhibition. The correlation between retention of Crf1p in related yeast strains and duplication of ribosomal protein genes suggests that this backup activation may be important to ensure gene expression when Ifh1p is limiting. With ribosome biogenesis as a hallmark of cell growth, failure to control assembly of ribosomal components leads to several human pathologies. A comprehensive understanding of mechanisms underlying this process is therefore of the essence.

Retention of duplicated genes in evolution

Gene duplication is a prevalent phenomenon across the tree of life. The processes that lead to the retention of duplicated genes are not well understood. Functional genomics approaches in model organisms, such as yeast, provide useful tools to test the mechanisms underlying retention with functional redundancy and divergence of duplicated genes, including fates associated with neofunctionalization, subfunctionalization, back-up compensation, and dosage amplification. Duplicated genes may also be retained as a consequence of structural and functional entanglement. Advances in human gene editing have enabled the interrogation of duplicated genes in the human genome, providing new tools to evaluate the relative contributions of each of these factors to duplicate gene retention and the evolution of genome structure.

The impact of whole genome duplications on the human gene regulatory networks

This work studies the effects of the two rounds of Whole Genome Duplication (WGD) at the origin of the vertebrate lineage on the architecture of the human gene regulatory networks. We integrate information on transcriptional regulation, miRNA regulation, and protein-protein interactions to comparatively analyse the role of WGD and Small Scale Duplications (SSD) in the structural properties of the resulting multilayer network. We show that complex network motifs, such as combinations of feed-forward loops and bifan arrays, deriving from WGD events are specifically enriched in the network. Pairs of WGD-derived proteins display a strong tendency to interact both with each other and with common partners and WGD-derived transcription factors play a prominent role in the retention of a strong regulatory redundancy. Combinatorial regulation and synergy between different regulatory layers are in general enhanced by duplication events, but the two types of duplications contribute in different ways. Overall, our findings suggest that the two WGD events played a substantial role in increasing the multi-layer complexity of the vertebrate regulatory network by enhancing its combinatorial organization, with potential consequences on its overall robustness and ability to perform high-level functions like signal integration and noise control. Lastly, we discuss in detail the RAR/RXR pathway as an illustrative example of the evolutionary impact of WGD duplications in human.

Gene duplication is one of the main mechanisms driving genome evolution. The duplication of a genomic segment can be the result of a local event, involving only a small portion of the genome, or of a dramatic duplication of the whole genome, which is however only rarely retained. All vertebrates descend from two rounds of Whole-Genome Duplication (WGD) that occurred approximately 500 Mya. We show that these events influenced in unique ways the evolution of different human gene regulatory networks, with sizeable effects on their current structure. We find that WGDs statistically increased the presence of specific classes of simple genetic circuits, considered to be fundamental building blocks of more sophisticated circuitry and commonly associated to complex functions. Our findings support the hypothesis that these rare, large-scale events have played a substantial role in the emergence of complex traits in vertebrates.

Independent control of mean and noise by convolution of gene expression distributions

Gene expression noise can reduce cellular fitness or facilitate processes such as alternative metabolism, antibiotic resistance, and differentiation. Unfortunately, efforts to study the impacts of noise have been hampered by a scaling relationship between noise and expression level from individual promoters. Here, we use theory to demonstrate that mean and noise can be controlled independently by expressing two copies of a gene from separate inducible promoters in the same cell. We engineer low and high noise inducible promoters to validate this result in Escherichia coli, and develop a model that predicts the experimental distributions. Finally, we use our method to reveal that the response of a promoter to a repressor is less sensitive with higher repressor noise and explain this result using a law from probability theory. Our approach can be applied to investigate the effects of noise on diverse biological pathways or program cellular heterogeneity for synthetic biology applications.

Nuclear compartmentalization as a mechanism of quantitative control of gene expression

Gene regulation requires the dynamic coordination of hundreds of regulatory factors at precise genomic and RNA targets. Although many regulatory factors have specific affinity for their nucleic acid targets, molecular diffusion and affinity models alone cannot explain many of the quantitative features of gene regulation in the nucleus. One emerging explanation for these quantitative properties is that DNA, RNA and proteins organize within precise, 3D compartments in the nucleus to concentrate groups of functionally related molecules. Recently, nucleic acids and proteins involved in many important nuclear processes have been shown to engage in cooperative interactions, which lead to the formation of condensates that partition the nucleus. In this Review, we discuss an emerging perspective of gene regulation, which moves away from classic models of stoichiometric interactions towards an understanding of how spatial compartmentalization can lead to non-stoichiometric molecular interactions and non-linear regulatory behaviours. We describe key mechanisms of nuclear compartment formation, including emerging roles for non-coding RNAs in facilitating their formation, and discuss the functional role of nuclear compartments in transcription regulation, co-transcriptional and post-transcriptional RNA processing, and higher-order chromatin regulation. More generally, we discuss how compartmentalization may explain important quantitative aspects of gene regulation.

The YTHDF proteins ECT2 and ECT3 bind largely overlapping target sets and influence target mRNA abundance, not alternative polyadenylation

Gene regulation via N6-methyladenosine (m⁶A) in mRNA involves RNA-binding proteins that recognize m⁶A via a YT521-B homology (YTH) domain. The plant YTH domain proteins ECT2 and ECT3 act genetically redundantly in stimulating cell proliferation during organogenesis, but several fundamental questions regarding their mode of action remain unclear. Here, we use HyperTRIBE (targets of RNA-binding proteins identified by editing) to show that most ECT2 and ECT3 targets overlap, with only a few examples of preferential targeting by either of the two proteins. HyperTRIBE in different mutant backgrounds also provides direct views of redundant, ectopic, and specific target interactions of the two proteins. We also show that contrary to conclusions of previous reports, ECT2 does not accumulate in the nucleus. Accordingly, inactivation of ECT2, ECT3, and their surrogate ECT4 does not change patterns of polyadenylation site choice in ECT2/3 target mRNAs, but does lead to lower steady-state accumulation of target mRNAs. In addition, mRNA and microRNA expression profiles show indications of stress response activation in ect2/ect3/ect4 mutants, likely via indirect effects. Thus, previous suggestions of control of alternative polyadenylation by ECT2 are not supported by evidence, and ECT2 and ECT3 act largely redundantly to regulate target mRNA, including its abundance, in the cytoplasm.

Measurement of histone replacement dynamics with genetically encoded exchange timers in yeast

Histone exchange between histones carrying position-specific marks and histones bearing general marks is important for gene regulation, but understanding of histone exchange remains incomplete. To overcome the poor time resolution of conventional pulse-chase histone labeling, we present a genetically encoded histone exchange timer sensitive to the duration that two tagged histone subunits co-reside at an individual genomic locus. We apply these sensors to map genome-wide patterns of histone exchange in yeast using single samples. Comparing H3 exchange in cycling and G1-arrested cells suggests that replication-independent H3 exchange occurs at several hundred nucleosomes (<1% of all nucleosomes) per minute, with a maximal rate at histone promoters. We observed substantial differences between the two nucleosome core subcomplexes: H2A-H2B subcomplexes undergo rapid transcription-dependent replacement within coding regions, whereas H3-H4 replacement occurs predominantly within promoter nucleosomes, in association with gene activation or repression. Our timers allow the in vivo study of histone exchange dynamics with minute time scale resolution.

Yeast cell fate control by temporal redundancy modulation of transcription factor paralogs

Recent single-cell studies have revealed that yeast stress response involves transcription factors that are activated in pulses. However, it remains unclear whether and how these dynamic transcription factors temporally interact to regulate stress survival. Here we show that budding yeast cells can exploit the temporal relationship between paralogous general stress regulators, Msn2 and Msn4, during stress response. We find that individual pulses of Msn2 and Msn4 are largely redundant, and cells can enhance the expression of their shared targets by increasing their temporal divergence. Thus, functional redundancy between these two paralogs is modulated in a dynamic manner to confer fitness advantages for yeast cells, which might feed back to promote the preservation of their redundancy. This evolutionary implication is supported by evidence from Msn2/Msn4 orthologs and analyses of other transcription factor paralogs. Together, we show a cell fate control mechanism through temporal redundancy modulation in yeast, which may represent an evolutionarily important strategy for maintaining functional redundancy between gene duplicates.

The evolution of red color vision is linked to coordinated rhodopsin tuning in lycaenid butterflies

Color vision has evolved multiple times in both vertebrates and invertebrates and is largely determined by the number and variation in spectral sensitivities of distinct opsin subclasses. However, because of the difficulty of expressing long-wavelength (LW) invertebrate opsins in vitro, our understanding of the molecular basis of functional shifts in opsin spectral sensitivities has been biased toward research primarily in vertebrates. This has restricted our ability to address whether invertebrate Gq protein-coupled opsins function in a novel or convergent way compared to vertebrate Gt opsins. Here we develop a robust heterologous expression system to purify invertebrate rhodopsins, identify specific amino acid changes responsible for adaptive spectral tuning, and pinpoint how molecular variation in invertebrate opsins underlie wavelength sensitivity shifts that enhance visual perception. By combining functional and optophysiological approaches, we disentangle the relative contributions of lateral filtering pigments from red-shifted LW and blue short-wavelength opsins expressed in distinct photoreceptor cells of individual ommatidia. We use in situ hybridization to visualize six ommatidial classes in the compound eye of a lycaenid butterfly with a four-opsin visual system. We show experimentally that certain key tuning residues underlying green spectral shifts in blue opsin paralogs have evolved repeatedly among short-wavelength opsin lineages. Taken together, our results demonstrate the interplay between regulatory and adaptive evolution at multiple Gq opsin loci, as well as how coordinated spectral shifts in LW and blue opsins can act together to enhance insect spectral sensitivity at blue and red wavelengths for visual performance adaptation.

Gene regulatory network inference as relaxed graph matching

Bipartite network inference is a ubiquitous problem across disciplines. One important example in the field molecular biology is gene regulatory network inference. Gene regulatory networks are an instrumental tool aiding in the discovery of the molecular mechanisms driving diverse diseases, including cancer. However, only noisy observations of the projections of these regulatory networks are typically assayed. In an effort to better estimate regulatory networks from their noisy projections, we formulate a non-convex but analytically tractable optimization problem called OTTER. This problem can be interpreted as relaxed graph matching between the two projections of the bipartite network. OTTER's solutions can be derived explicitly and inspire a spectral algorithm, for which we provide network recovery guarantees. We also provide an alternative approach based on gradient descent that is more robust to noise compared to the spectral algorithm. Interestingly, this gradient descent approach resembles the message passing equations of an established gene regulatory network inference method, PANDA. Using three cancer-related data sets, we show that OTTER outperforms state-of-the-art inference methods in predicting transcription factor binding to gene regulatory regions. To encourage new graph matching applications to this problem, we have made all networks and validation data publicly available.

Speed-Specificity Trade-Offs in the Transcription Factors Search for Their Genomic Binding Sites

Transcription factors (TFs) regulate gene expression by binding DNA sequences recognized by their DNA-binding domains (DBDs). DBD-recognized motifs are short and highly abundant in genomes. The ability of TFs to bind a specific subset of motif-containing sites, and to do so rapidly upon activation, is fundamental for gene expression in all eukaryotes. Despite extensive interest, our understanding of the TF-target search process is fragmented; although binding specificity and detection speed are two facets of this same process, trade-offs between them are rarely addressed. In this opinion article, we discuss potential speed-specificity trade-offs in the context of existing models. We further discuss the recently described 'distributed specificity' paradigm, suggesting that intrinsically disordered regions (IDRs) promote specificity while reducing the TF-target search time.

Genome-Wide Profiling of Protein-DNA Interactions with Chromatin Endogenous Cleavage and High-Throughput Sequencing (ChEC-Seq )

Interactions between regulatory proteins and specific genomic regions are critical for all chromatin-based processes, including transcription, DNA replication, and DNA repair. Genome-wide mapping of such interactions is most commonly performed with chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-Seq), but a number of orthogonal methods employing targeted enzymatic activity have also been introduced. We previously described a genome-wide implementation of chromatin endogenous cleavage (ChEC-Seq), wherein a protein of interest is fused to micrococcal nuclease (MNase) to enable targeted, calcium-dependent genomic cleavage. Here, we describe the ChEC-Seq protocol for use in budding yeast though it can be used in other organisms in conjunction with appropriate methods for introduction of an MNase fusion protein.

Molecular and evolutionary processes generating variation in gene expression

Heritable variation in gene expression is common within and between species. This variation arises from mutations that alter the form or function of molecular gene regulatory networks that are then filtered by natural selection. High-throughput methods for introducing mutations and characterizing their cis- and trans-regulatory effects on gene expression (particularly, transcription) are revealing how different molecular mechanisms generate regulatory variation, and studies comparing these mutational effects with variation seen in the wild are teasing apart the role of neutral and non-neutral evolutionary processes. This integration of molecular and evolutionary biology allows us to understand how the variation in gene expression we see today came to be and to predict how it is most likely to evolve in the future.

Expression Partitioning of Duplicate Genes at Single Cell Resolution in Arabidopsis Roots

Gene duplication is a key evolutionary phenomenon, prevalent in all organisms but particularly so in plants, where whole genome duplication (WGD; polyploidy) is a major force in genome evolution. Much effort has been expended in attempting to understand the evolution of duplicate genes, addressing such questions as why some paralog pairs rapidly return to single copy status whereas, in other pairs, both paralogs are retained and may diverge in expression pattern or function. The effect of a gene - its site of expression and thus the initial locus of its function - occurs at the level of a cell comprising a single cell type at a given state of the cell's development. Using Arabidopsis thaliana single cell transcriptomic data we categorized patterns of expression for 11,470 duplicate gene pairs across 36 cell clusters comprising nine cell types and their developmental states. Among these 11,470 pairs, 10,187 (88.8%) had at least one copy expressed in at least one of the 36 cell clusters. Pairs produced by WGD more often had both paralogs expressed in root cells than did pairs produced by small scale duplications. Three quarters of gene pairs expressed in the 36 cell clusters (7,608/10,187) showed extreme expression bias in at least one cluster, including 352 cases of reciprocal bias, a pattern consistent with expression subfunctionalization. More than twice as many pairs showed reciprocal expression bias between cell states than between cell types or between roots and leaves. A group of 33 gene pairs with reciprocal expression bias showed evidence of concerted divergence of gene networks in stele vs. epidermis. Pairs with both paralogs expressed without bias were less likely to have paralogs with divergent mutant phenotypes; such bias-free pairs showed evidence of preservation by maintenance of dosage balance. Overall, we found considerable evidence of shifts in gene expression following duplication, including in >80% of pairs encoding 7,653 genes expressed ubiquitously in all root cell types and states for which we inferred the polarity of change.

Beyond Trees: Regulons and Regulatory Motif Characterization

Trees and their seeds regulate their germination, growth, and reproduction in response to environmental stimuli. These stimuli, through signal transduction, trigger transcription factors that alter the expression of various genes leading to the unfolding of the genetic program. A regulon is conceptually defined as a set of target genes regulated by a transcription factor by physically binding to regulatory motifs to accomplish a specific biological function, such as the CO-FT regulon for flowering timing and fall growth cessation in trees. Only with a clear characterization of regulatory motifs, can candidate target genes be experimentally validated, but motif characterization represents the weakest feature of regulon research, especially in tree genetics. I review here relevant experimental and bioinformatics approaches in characterizing transcription factors and their binding sites, outline problems in tree regulon research, and demonstrate how transcription factor databases can be effectively used to aid the characterization of tree regulons.

Transcriptional Rewiring, Adaptation, and the Role of Gene Duplication in the Metabolism of Ethanol of Saccharomyces cerevisiae

Ethanol is the main by-product of yeast sugar fermentation that affects microbial growth parameters, being considered a dual molecule, a nutrient and a stressor. Previous works demonstrated that the budding yeast arose after an ancient hybridization process resulted in a tier of duplicated genes within its genome, many of them with implications in this ethanol "produce-accumulate-consume" strategy. The evolutionary link between ethanol production, consumption, and tolerance versus ploidy and stability of the hybrids is an ongoing debatable issue. The implication of ancestral duplicates in this metabolic rewiring, and how these duplicates differ transcriptionally, remains unsolved. Here, we study the transcriptomic adaptive signatures to ethanol as a nonfermentative carbon source to sustain clonal yeast growth by experimental evolution, emphasizing the role of duplicated genes in the adaptive process. As expected, ethanol was able to sustain growth but at a lower rate than glucose. Our results demonstrate that in asexual populations a complete transcriptomic rewiring was produced, strikingly by downregulation of duplicated genes, mainly whole-genome duplicates, whereas small-scale duplicates exhibited significant transcriptional divergence between copies. Overall, this study contributes to the understanding of evolution after gene duplication, linking transcriptional divergence with duplicates' fate in a multigene trait as ethanol tolerance.IMPORTANCE Gene duplication events have been related with increasing biological complexity through the tree of life, but also with illnesses, including cancer. Early evolutionary theories indicated that duplicated genes could explore alternative functions due to relaxation of selective constraints in one of the copies, as the other remains as ancestral-function backup. In unicellular eukaryotes like yeasts, it has been demonstrated that the fate and persistence of duplicates depend on duplication mechanism (whole-genome or small-scale events), shaping their actual genomes. Although it has been shown that small-scale duplicates tend to innovate and whole-genome duplicates specialize in ancestral functions, the implication of duplicates' transcriptional plasticity and transcriptional divergence on environmental and metabolic responses remains largely obscure. Here, by experimental adaptive evolution, we show that Saccharomyces cerevisiae is able to respond to metabolic stress (ethanol as nonfermentative carbon source) due to the persistence of duplicated genes. These duplicates respond by transcriptional rewiring, depending on their transcriptional background. Our results shed light on the mechanisms that determine the role of duplicates, and on their evolvability.

The origin and evolution of vertebrate neural crest cells

The neural crest is a vertebrate-specific migratory stem cell population that generates a remarkably diverse set of cell types and structures. Because many of the morphological, physiological and behavioural novelties of vertebrates are derived from neural crest cells, it is thought that the origin of this cell population was an important milestone in early vertebrate history. An outstanding question in the field of vertebrate evolutionary-developmental biology (evo-devo) is how this cell type evolved in ancestral vertebrates. In this review, we briefly summarize neural crest developmental genetics in vertebrates, focusing in particular on the gene regulatory interactions instructing their early formation within and migration from the dorsal neural tube. We then discuss how studies searching for homologues of neural crest cells in invertebrate chordates led to the discovery of neural crest-like cells in tunicates and the potential implications this has for tracing the pre-vertebrate origins of the neural crest population. Finally, we synthesize this information to propose a model to explain the origin of neural crest cells. We suggest that at least some of the regulatory components of early stages of neural crest development long pre-date vertebrate origins, perhaps dating back to the last common bilaterian ancestor. These components, originally directing neuroectodermal patterning and cell migration, served as a gene regulatory 'scaffold' upon which neural crest-like cells with limited migration and potency evolved in the last common ancestor of tunicates and vertebrates. Finally, the acquisition of regulatory programmes controlling multipotency and long-range, directed migration led to the transition from neural crest-like cells in invertebrate chordates to multipotent migratory neural crest in the first vertebrates.

Regulation plays a multifaceted role in the retention of gene duplicates

A gene duplication can lead to all sorts of problems in a cell. However, it can also lead to all sorts of benefits. Beneficial or not, the gene duplicates might be kept in the genome because of several different reasons. For instance, if natural selection works towards optimizing one function of a gene at the expense of another, then gene duplication could resolve this conflict by separating the functions in two genes. Here, we outline evolutionary incentives to keep a duplicated gene in the genome, focusing on divergence in expression and trade-off resolution as featured in a new and exciting paper published in this edition of PLOS Biology.