Redundancy and the evolution of cis-regulatory element multiplicity

Comprehensive in-silico characterization and expression pattern of calmodulin genes under various abiotic and biotic stresses in Indian mustard (Brassica juncea)

Calcium (Ca2⁺) as a secondary messenger has a multidimensional role, including the growth and development of plants and the adaptive response to stress conditions. Calmodulin (CaM), a calcium-binding protein, uniquely binds with these Ca2⁺ ions and transmits Ca2⁺ signals. Calmodulin proteins have been well-reported in various plants for playing a role in abiotic and biotic stress signaling; however, a comprehensive analysis of the CaM genes of Indian mustard (Brassica juncea) has not been studied much. This study reports their chromosome placements, phylogenetic relations, the presence of protein motifs and cis-acting elements, and their expression patterns under stress due to salt, heat, cadmium, Xanthomonas campestris, and Alternaria brassicae. We identified 23 BjCaM genes coding for eight BjCaM proteins possessing the signature EF-hand domains. Chromosome locations, intron-exon structure, and in-silico protein characterization pointed toward genetic diversification. Phylogenetic analysis revealed a close relationship with previously characterized CaM proteins from Arabidopsis and rice. Cis-acting elements in the promoter regions suggested the potential role of BjCaM candidates in hormone signaling and various stress-responsive regulatory mechanisms. qRT-PCR analysis showed differential expression patterns, of which BjCaM17 and BjCaM19 showed higher expression under all stresses. The seven selected BjCaM genes were sensitive to cadmium stress. Interestingly, despite translating to same protein, BjCaM15, BjCaM17, and BjCaM19 showed differential expressions under the same stresses. This research represents the first genome-wide analysis of calmodulin genes in Indian mustard, providing a valuable reference for decoding calcium signaling via calmodulin and its potential exploitation to improve crop resistance to stress conditions.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-025-01561-x.

On the importance of evolutionary constraint for regulatory sequence identification

Regulation of gene expression relies on the activity of specialized genomic elements, enhancers or silencers, distributed over sometimes large distance from their target gene promoters. A significant part of vertebrate genomes consists in such regulatory elements, but their identification and that of their target genes remains challenging, due to the lack of clear signature at the nucleotide level. For many years the main hallmark used for identifying functional elements has been their sequence conservation between genomes of distant species, indicative of purifying selection. More recently, genome-wide biochemical assays have opened new avenues for detecting regulatory regions, shifting attention away from evolutionary constraints. Here, we review the respective contributions of comparative genomics and biochemical assays for the definition of regulatory elements and their targets and advocate that both sequence conservation and preserved synteny, taken as signature of functional constraint, remain essential tools in this task.

Dual Eigen-modules of Cis-Element Regulation Profiles and Selection of Cognition-Language Eigen-direction along Evolution in Hominidae

To understand the genomic basis accounting for the phenotypic differences between human and apes, we compare the matrices consisting of the cis-element frequencies in the proximal regulatory regions of their genomes. One such frequency matrix is represented by a robust singular value decomposition. For each singular value, the negative and positive ends of the sorted motif eigenvector correspond to the dual ends of the sorted gene eigenvector, respectively, comprising a dual eigen-module defined by cis-regulatory element frequencies (CREF). The CREF eigen-modules at levels 1, 2, 3, and 6 are highly conserved across humans, chimpanzees, and orangutans. The key biological processes embedded in the top three CREF eigen-modules are reproduction versus embryogenesis, fetal maturation versus immune system, and stress responses versus mitosis. Although the divergence at the nucleotide level between the chimpanzee and human genome was small, their cis-element frequency matrices crossed a singularity point, at which the fourth and fifth singular values were identical. The CREF eigen-modules corresponding to the fourth and fifth singular values were reorganized along the evolution from apes to human. Interestingly, the fourth sorted gene eigenvector encodes the phenotypes unique to human such as long-term memory, language development, and social behavior. The number of motifs present on Alu elements increases substantially at the fourth level. The motif analysis together with the cases of human-specific Alu insertions suggests that mutations related to Alu elements play a critical role in the evolution of the human-phenotypic gene eigenvector.

The potential application of genome editing by using CRISPR/Cas9, and its engineered and ortholog variants for studying the transcription factors involved in the maintenance of phosphate homeostasis in model plants

Phosphorus (P), an essential macronutrient, is pivotal for growth and development of plants. Availability of phosphate (Pi), the only assimilable P, is often suboptimal in rhizospheres. Pi deficiency triggers an array of spatiotemporal adaptive responses including the differential regulation of several transcription factors (TFs). Studies on MYB TF PHR1 in Arabidopsis thaliana (Arabidopsis) and its orthologs OsPHRs in Oryza sativa (rice) have provided empirical evidence of their significant roles in the maintenance of Pi homeostasis. Since the functional characterization of PHR1 in 2001, several other TFs have now been identified in these model plants. This raised a pertinent question whether there are any likely interactions across these TFs. Clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated protein 9 (Cas9) system has provided an attractive paradigm for editing genome in plants. Here, we review the applications and challenges of this technique for genome editing of the TFs for deciphering the function and plausible interactions across them. This technology could thus provide a much-needed fillip towards engineering TFs for generating Pi use efficient plants for sustainable agriculture. Furthermore, we contemplate whether this technology could be a viable alternative to the controversial genetically modified (GM) rice or it may also eventually embroil into a limbo.

Making sense of transcription networks

When transcription regulatory networks are compared among distantly related eukaryotes, a number of striking similarities are observed: a larger-than-expected number of genes, extensive overlapping connections, and an apparently high degree of functional redundancy. It is often assumed that the complexity of these networks represents optimized solutions, precisely sculpted by natural selection; their common features are often asserted to be adaptive. Here, we discuss support for an alternative hypothesis: the common structural features of transcription networks arise from evolutionary trajectories of "least resistance"--that is, the relative ease with which certain types of network structures are formed during their evolution.

Mechanisms of mutational robustness in transcriptional regulation

Robustness is the invariance of a phenotype in the face of environmental or genetic change. The phenotypes produced by transcriptional regulatory circuits are gene expression patterns that are to some extent robust to mutations. Here we review several causes of this robustness. They include robustness of individual transcription factor binding sites, homotypic clusters of such sites, redundant enhancers, transcription factors, redundant transcription factors, and the wiring of transcriptional regulatory circuits. Such robustness can either be an adaptation by itself, a byproduct of other adaptations, or the result of biophysical principles and non-adaptive forces of genome evolution. The potential consequences of such robustness include complex regulatory network topologies that arise through neutral evolution, as well as cryptic variation, i.e., genotypic divergence without phenotypic divergence. On the longest evolutionary timescales, the robustness of transcriptional regulation has helped shape life as we know it, by facilitating evolutionary innovations that helped organisms such as flowering plants and vertebrates diversify.

Spurious transcription factor binding: non-functional or genetically redundant?

Transcription factor binding sites (TFBSs) on the DNA are generally accepted as the key nodes of gene control. However, the multitudes of TFBSs identified in genome-wide studies, some of them seemingly unconstrained in evolution, have prompted the view that in many cases TF binding may serve no biological function. Yet, insights from transcriptional biochemistry, population genetics and functional genomics suggest that rather than segregating into 'functional' or 'non-functional', TFBS inputs to their target genes may be generally cumulative, with varying degrees of potency and redundancy. As TFBS redundancy can be diminished by mutations and environmental stress, some of the apparently 'spurious' sites may turn out to be important for maintaining adequate transcriptional regulation under these conditions. This has significant implications for interpreting the phenotypic effects of TFBS mutations, particularly in the context of genome-wide association studies for complex traits.

Identifying Cis-Regulatory Elements and Modules Using Conditional Random Fields

Accurate identification of cis-regulatory elements and their correlated modules is essential for analysis of transcriptional regulation, which is a challenging problem in computational biology. Unsupervised learning has the advantage of compensating for missing annotated data, and is thus promising to be effective to identify cis-regulatory elements and modules. We introduced a Conditional Random Fields model, referred to as CRFEM, to integrate sequence features and long-range dependency of genomic sequences such as epigenetic features to identify cis-regulatory elements and modules at the same time. The proposed method is able to automatically learn model parameters with no labeled data and explicitly optimize the predictive probability of cis-regulatory elements and modules. In comparison with existing methods, our method is more accurate and can be used for genome-wide studies of gene regulation.

The evolution of complex gene regulation by low-specificity binding sites

Requirements for gene regulation vary widely both within and among species. Some genes are constitutively expressed, whereas other genes require complex regulatory control. Transcriptional regulation is often controlled by a module of multiple transcription factor binding sites that, in combination, mediate the expression of a target gene. Here, we study how such regulatory modules evolve in response to natural selection. Using a population-genetic model, we show that complex regulatory modules which contain a larger number of binding sites must employ binding motifs that are less specific, on average, compared with smaller regulatory modules. This effect is extremely general, and it holds regardless of the selected binding logic that a module experiences. We attribute this phenomenon to the inability of stabilizing selection to maintain highly specific sites in large regulatory modules. Our analysis helps to explain broad empirical trends in the Saccharomyces cerevisiae regulatory network: those genes with a greater number of distinct transcriptional regulators feature less-specific binding motifs, compared with genes with fewer regulators. Our results also help to explain empirical trends in module size and motif specificity across species, ranging from prokaryotes to single-cellular and multi-cellular eukaryotes.

The underlying molecular and network level mechanisms in the evolution of robustness in gene regulatory networks

Gene regulatory networks show robustness to perturbations. Previous works identified robustness as an emergent property of gene network evolution but the underlying molecular mechanisms are poorly understood. We used a multi-tier modeling approach that integrates molecular sequence and structure information with network architecture and population dynamics. Structural models of transcription factor-DNA complexes are used to estimate relative binding specificities. In this model, mutations in the DNA cause changes on two levels: (a) at the sequence level in individual binding sites (modulating binding specificity), and (b) at the network level (creating and destroying binding sites). We used this model to dissect the underlying mechanisms responsible for the evolution of robustness in gene regulatory networks. Results suggest that in sparse architectures (represented by short promoters), a mixture of local-sequence and network-architecture level changes are exploited. At the local-sequence level, robustness evolves by decreasing the probabilities of both the destruction of existent and generation of new binding sites. Meanwhile, in highly interconnected architectures (represented by long promoters), robustness evolves almost entirely via network level changes, deleting and creating binding sites that modify the network architecture.

Lengthening of 3'UTR increases with morphological complexity in animal evolution

MOTIVATION

Evolutionary expansion of gene regulatory circuits seems to boost morphological complexity. However, the expansion patterns and the quantification relationships have not yet been identified. In this study, we focus on the regulatory circuits at the post-transcriptional level, investigating whether and how this principle may apply.

RESULTS

By analysing the structure of mRNA transcripts in multiple metazoan species, we observed a striking exponential correlation between the length of 3' untranslated regions (3'UTR) and morphological complexity as measured by the number of cell types in each organism. Cellular diversity was similarly associated with the accumulation of microRNA genes and their putative targets. We propose that the lengthening of 3'UTRs together with a commensurate exponential expansion in post-transcriptional regulatory circuits can contribute to the emergence of new cell types during animal evolution.

Coevolution within and between regulatory loci can preserve promoter function despite evolutionary rate acceleration

Phenotypes that appear to be conserved could be maintained not only by strong purifying selection on the underlying genetic systems, but also by stabilizing selection acting via compensatory mutations with balanced effects. Such coevolution has been invoked to explain experimental results, but has rarely been the focus of study. Conserved expression driven by the unc-47 promoters of Caenorhabditis elegans and C. briggsae persists despite divergence within a cis-regulatory element and between this element and the trans-regulatory environment. Compensatory changes in cis and trans are revealed when these promoters are used to drive expression in the other species. Functional changes in the C. briggsae promoter, which has experienced accelerated sequence evolution, did not lead to alteration of gene expression in its endogenous environment. Coevolution among promoter elements suggests that complex epistatic interactions within cis-regulatory elements may facilitate their divergence. Our results offer a detailed picture of regulatory evolution in which subtle, lineage-specific, and compensatory modifications of interacting cis and trans regulators together maintain conserved gene expression patterns.

Some phenotypes, including gene expression patterns, are conserved between distantly related species. However, the molecular bases of those phenotypes are not necessarily conserved. Instead, regulatory DNA sequences and the proteins with which they interact can change over time with balanced effects, preserving expression patterns and concealing regulatory divergence. Coevolution between interacting molecules makes gene regulation highly species-specific, and it can be detected when the cis-regulatory DNA of one species is used to drive expression in another species. In this way, we identified regions of the C. elegans and C. briggsae unc-47 promoters that have coevolved with the lineage-specific trans-regulatory environments of these organisms. The C. briggsae promoter experienced accelerated sequence change relative to related species. All of this evolution occurred without changing the expression pattern driven by the promoter in its endogenous environment.

ncDNA and drift drive binding site accumulation

Background

The amount of transcription factor binding sites (TFBS) in an organism’s genome positively correlates with the complexity of the regulatory network of the organism. However, the manner by which TFBS arise and accumulate in genomes and the effects of regulatory network complexity on the organism’s fitness are far from being known. The availability of TFBS data from many organisms provides an opportunity to explore these issues, particularly from an evolutionary perspective.

Results

We analyzed TFBS data from five model organisms – E. coli K12, S. cerevisiae, C. elegans, D. melanogaster, A. thaliana – and found a positive correlation between the amount of non-coding DNA (ncDNA) in the organism’s genome and regulatory complexity. Based on this finding, we hypothesize that the amount of ncDNA, combined with the population size, can explain the patterns of regulatory complexity across organisms. To test this hypothesis, we devised a genome-based regulatory pathway model and subjected it to the forces of evolution through population genetic simulations. The results support our hypothesis, showing neutral evolutionary forces alone can explain TFBS patterns, and that selection on the regulatory network function does not alter this finding.

Conclusions

The cis-regulome is not a clean functional network crafted by adaptive forces alone, but instead a data source filled with the noise of non-adaptive forces. From a regulatory perspective, this evolutionary noise manifests as complexity on both the binding site and pathway level, which has significant implications on many directions in microbiology, genetics, and synthetic biology.

Escape from Adaptive Conflict follows from weak functional trade-offs and mutational robustness

A fundamental question in molecular evolution is how proteins can adapt to new functions while being conserved for an existing function at the same time. Several theoretical models have been put forward to explain this apparent paradox. The most popular models include neofunctionalization, subfunctionalization (SUBF) by degenerative mutations, and dosage models. All of these models focus on adaptation after gene duplication. A newly proposed model named "Escape from Adaptive Conflict" (EAC) includes adaptive processes before and after gene duplication that lead to multifunctional proteins, and divergence (SUBF). Support for the importance of multifunctionality for the evolution of new protein functions comes from two experimental observations. First, many enzymes have highly evolvable promiscuous side activities. Second, different structural states of the same protein can be associated with different functions. How these observations may be related to the EAC model, under which conditions EAC is possible, and how the different models relate to each other is still unclear. Here, we present a theoretical framework that uses biophysical principles to infer the roles of functional promiscuity, gene dosage, gene duplication, point mutations, and selection pressures in the evolution of proteins. We find that selection pressures can determine whether neofunctionalization or SUBF is the more likely evolutionary process. Multifunctional proteins, arising during EAC evolution, allow rapid adaptation independent of gene duplication. This becomes a crucial advantage when gene duplications are rare. Finally, we propose that an increase in mutational robustness, not necessarily functional optimization, can be the sole driving force behind SUBF.

Evidence for directed evolution of larger size motif in Arabidopsis thaliana genome

Transcription control of gene expression depends on a variety of interactions mediated by the core promoter region, sequence specific DNA-binding proteins, and their cognate promoter elements. The prominent group of cis acting elements in plants contains an ACGT core. The cis element with this core has been shown to be involved in abscisic acid, salicylic acid, and light response. In this study, genome-wide comparison of the frequency of occurrence of two ACGT elements without any spacers as well as those separated by spacers of different length was carried out. In the first step, the frequency of occurrence of the cis element sequences across the whole genome was determined by using BLAST tool. In another approach the spacer sequence was randomized before making the query. As expected, the sequence ACGTACGT had maximum occurrence in Arabidopsis thaliana genome. As we increased the spacer length, one nucleotide at a time, the probability of its occurrence in genome decreased. This trend continued until an unexpectedly sharp rise in frequency of (ACGT)N25(ACGT). The observation of higher probability of bigger size motif suggests its directed evolution in Arabidopsis thaliana genome.

Shadow enhancers: frequently asked questions about distributed cis-regulatory information and enhancer redundancy

This paper, in the form of a frequently asked questions page (FAQ), addresses outstanding questions about "shadow enhancers", quasi-redundant cis-regulatory elements, and their proposed roles in transcriptional control. Questions include: What exactly are shadow enhancers? How many genes have shadow/redundant/distributed enhancers? How redundant are these elements? What is the function of distributed enhancers? How modular are enhancers? Is it useful to study a single enhancer in isolation? In addition, a revised definition of "shadow enhancers" is proposed, and possible mechanisms of shadow enhancer function and evolution are discussed.

Evolutionary origins of transcription factor binding site clusters

Empirical studies have revealed that regulatory DNA sequences such as enhancers or promoters often harbor multiple binding sites for the same transcription factor. Such "homotypic site clustering" has been hypothesized as arising out of functional requirements of the sequences. Here, we propose an alternative explanation of this phenomenon that multisite enhancers are common because they are favored by evolutionary sampling of the genotype-phenotype landscape. To test this hypothesis, we developed a new computational framework specialized for population genetic simulations of enhancer evolution. It uses a thermodynamics-based model of enhancer function, integrating information from strong as well as weak binding sites, to determine the strength of selection. Using this framework, we found that even when simpler genotypes exist for a desired strength of regulation, relatively complex genotypes (enhancers with more sites) are more readily reached by the simulated evolutionary process. We show that there are more ways to "build" a fit genotype with many weak sites than with a few strong sites, and this is why evolution finds complex genotypes more often. Our claims are consistent with an empirical analysis of binding site content in enhancers characterized in Drosophila melanogaster and their orthologs in other Drosophila species. We also characterized a subtle but significant difference between genotypes likely to be sampled by evolution and equally fit genotypes one would obtain by uniform sampling of the fitness landscape, that is, an "evolutionary signature" in enhancer sequences. Finally, we investigated potential effects of other factors, such as rugged fitness landscapes, short local duplications, and noise characteristics of enhancers, on the emergence of homotypic site clustering. Homotypic site clustering is an important contributor to the complexity and function of cis-regulatory sequences. This work provides a simple null hypothesis for its origin, against which alternative adaptationist explanations may be evaluated, and cautions against "evolutionary mirages" present in common features of genomic sequence. The quantitative framework we develop here can be used more generally to understand how mechanisms of enhancer action influence their composition and evolution.

Adult Drosophila melanogaster evolved for antibacterial defense invest in infection-induced expression of both humoral and cellular immunity genes

Background

While the transcription of innate immunity genes in response to bacterial infection has been well-characterised in the Drosophila model, we recently demonstrated the capacity for such transcription to evolve in flies selected for improved antibacterial defense. Here we use this experimental system to examine how insects invest in constitutive versus infection-induced transcription of immunity genes. These two strategies carry with them different consequences with respect to energetic and pleiotropic costs and may be more or less effective in improving defense depending on whether the genes contribute to humoral or cellular aspects of immunity.

Findings

Contrary to expectation we show that selection preferentially increased the infection-induced expression of both cellular and humoral immunity genes. Given their functional roles, infection induced increases in expression were expected for the humoral genes, while increases in constitutive expression were expected for the cellular genes. We also report a restricted ability to improve transcription of immunity genes that is on the order of 2-3 fold regardless of total transcription level of the gene.

Conclusions

The evolved increases in infection-induced expression of the cellular genes may result from specific cross talk with humoral pathways or from generalised strategies for enhancing immunity gene transcription. A failure to see improvements in constitutive expression of the cellular genes suggests either that increases might come at too great a cost or that patterns of expression in adults are decoupled from the larval phase where increases would be most effective. The similarity in fold change increase across all immunity genes may suggest a shared mechanism for the evolution of increased transcription in small, discrete units such as duplication of cis-regulatory elements.