Enhancer evolution and the origins of morphological novelty

Blooming balloons: Searching for mechanisms of the inflated calyx

Studying morphological novelties offers special insights into developmental biology and evolution. The inflated calyx syndrome (ICS) is a largely unrecognized but fascinating feature of flower development, where sepals form balloon-like husks that encapsulate fruits. Despite its independent emergence in many lineages of flowering plants, the genetic and molecular mechanisms of ICS remain unknown. Early studies in the Solanaceae genus Physalis put forth key roles of MADS-box genes in ICS. However, recent work suggests these classical floral identity transcription factors were false leads. With newfound capabilities that allow rapid development of genetic systems through genomics and genome editing, Physalis has re-emerged as the most tractable model species for dissecting ICS. This review revisits current understanding of ICS and highlights how recent advancements enable a reset in the search for genetic and molecular mechanisms using unbiased, systematic approaches.

Evolution of diapause in the African turquoise killifish by remodeling the ancient gene regulatory landscape

Suspended animation states allow organisms to survive extreme environments. The African turquoise killifish has evolved diapause as a form of suspended development to survive a complete drought. However, the mechanisms underlying the evolution of extreme survival states are unknown. To understand diapause evolution, we performed integrative multi-omics (gene expression, chromatin accessibility, and lipidomics) in the embryos of multiple killifish species. We find that diapause evolved by a recent remodeling of regulatory elements at very ancient gene duplicates (paralogs) present in all vertebrates. CRISPR-Cas9-based perturbations identify the transcription factors REST/NRSF and FOXOs as critical for the diapause gene expression program, including genes involved in lipid metabolism. Indeed, diapause shows a distinct lipid profile, with an increase in triglycerides with very-long-chain fatty acids. Our work suggests a mechanism for the evolution of complex adaptations and offers strategies to promote long-term survival by activating suspended animation programs in other species.

CpG island turnover events predict evolutionary changes in enhancer activity

BACKGROUND

Genetic changes that modify the function of transcriptional enhancers have been linked to the evolution of biological diversity across species. Multiple studies have focused on the role of nucleotide substitutions, transposition, and insertions and deletions in altering enhancer function. CpG islands (CGIs) have recently been shown to influence enhancer activity, and here we test how their turnover across species contributes to enhancer evolution.

RESULTS

We integrate maps of CGIs and enhancer activity-associated histone modifications obtained from multiple tissues in nine mammalian species and find that CGI content in enhancers is strongly associated with increased histone modification levels. CGIs show widespread turnover across species and species-specific CGIs are strongly enriched for enhancers exhibiting species-specific activity across all tissues and species. Genes associated with enhancers with species-specific CGIs show concordant biases in their expression, supporting that CGI turnover contributes to gene regulatory innovation. Our results also implicate CGI turnover in the evolution of Human Gain Enhancers (HGEs), which show increased activity in human embryonic development and may have contributed to the evolution of uniquely human traits. Using a humanized mouse model, we show that a highly conserved HGE with a large CGI absent from the mouse ortholog shows increased activity at the human CGI in the humanized mouse diencephalon.

CONCLUSIONS

Collectively, our results point to CGI turnover as a mechanism driving gene regulatory changes potentially underlying trait evolution in mammals.

The paradise fish, an advanced animal model for behavioral genetics and evolutionary developmental biology

Paradise fish (Macropodus opercularis) is an air-breathing freshwater fish species with a signature labyrinth organ capable of extracting oxygen from the air that helps these fish to survive in hypoxic environments. The appearance of this evolutionary innovation in anabantoids resulted in a rewired circulatory system, but also in the emergence of species-specific behaviors, such as territorial display, courtship and parental care in the case of the paradise fish. Early zoologists were intrigued by the structure and function of the labyrinth apparatus and a series of detailed descriptive histological studies at the beginning of the 20th century revealed the ontogenesis and function of this specialized system. A few decades later, these fish became the subject of numerous ethological studies, and detailed ethograms of their behavior were constructed. These latter studies also demonstrated a strong genetic component underlying their behavior, but due to lack of adequate molecular tools, the fine genetic dissection of the behavior was not possible at the time. The technological breakthroughs that transformed developmental biology and behavioral genetics in the past decades, however, give us now a unique opportunity to revisit these old questions. Building on the classic descriptive studies, the new methodologies will allow us to follow the development of the labyrinth apparatus at a cellular resolution, reveal the genes involved in this process and also the genetic architecture behind the complex behaviors that we can observe in this species.

Gene regulatory network co-option is sufficient to induce a morphological novelty in Drosophila

Identifying the molecular origins by which new morphological structures evolve is one of the long standing problems in evolutionary biology. To date, vanishingly few examples provide a compelling account of how new morphologies were initially formed, thereby limiting our understanding of how diverse forms of life derived their complex features. Here, we provide evidence that the large projections on the Drosophila eugracilis phallus that are implicated in sexual conflict have evolved through co-option of the trichome genetic network. These unicellular apical projections on the phallus postgonal sheath are reminiscent of trichomes that cover the Drosophila body but are up to 20-fold larger in size. During their development, they express the transcription factor Shavenbaby, the master regulator of the trichome network. Consistent with the co-option of the Shavenbaby network during the evolution of the D. eugracilis projections, somatic mosaic CRISPR/Cas9 mutagenesis shows that shavenbaby is necessary for their proper length. Moreover, mis-expression of Shavenbaby in the sheath of D. melanogaster , a naïve species that lacks these extensions, is sufficient to induce small trichomes. These induced extensions rely on a genetic network that is shared to a large extent with the D. eugracilis projections, indicating its co-option but also some genetic rewiring. Thus, by leveraging a genetically tractable evolutionarily novelty, our work shows that the trichome-forming network is flexible enough that it can be co-opted in a new context, and subsequently refined to produce unique apical projections that are barely recognizable compared to their simpler ancestral beginnings.

Sox21b underlies the rapid diversification of a novel male genital structure between Drosophila species

The emergence and diversification of morphological novelties is a major feature of animal evolution.1,2,3,4,5,6,7,8,9 However, relatively little is known about the genetic basis of the evolution of novel structures and the mechanisms underlying their diversification. The epandrial posterior lobes of male genitalia are a novelty of particular Drosophila species.10,11,12,13 The lobes grasp the female ovipositor and insert between her abdominal tergites and, therefore, are important for copulation and species recognition.10,11,12,14,15,16,17 The posterior lobes likely evolved from co-option of a Hox-regulated gene network from the posterior spiracles10 and have since diversified in morphology in the D. simulans clade, in particular, over the last 240,000 years, driven by sexual selection.18,19,20,21 The genetic basis of this diversification is polygenic but, to the best of our knowledge, none of the causative genes have been identified.22,23,24,25,26,27,28,29,30 Identifying the genes underlying the diversification of these secondary sexual structures is essential to understanding the evolutionary impact on copulation and species recognition. Here, we show that Sox21b negatively regulates posterior lobe size. This is consistent with expanded Sox21b expression in D. mauritiana, which develops smaller posterior lobes than D. simulans. We tested this by generating reciprocal hemizygotes and confirmed that changes in Sox21b underlie posterior lobe evolution between these species. Furthermore, we found that posterior lobe size differences caused by the species-specific allele of Sox21b significantly affect copulation duration. Taken together, our study reveals the genetic basis for the sexual-selection-driven diversification of a novel morphological structure and its functional impact on copulatory behavior.

The interaction between lncRNAs and transcription factors regulating autophagy in human cancers: A comprehensive and therapeutical survey

Autophagy, as a highly conserved cellular process, participates in cellular homeostasis by degradation and recycling of damaged organelles and proteins. Besides, autophagy has been evidenced to play a dual role through cancer initiation and progression. In the early stage, it may have a tumor-suppressive function through inducing apoptosis and removing damaged cells and organelles. However, late stages promote tumor progression by maintaining stemness features and induction of chemoresistance. Therefore, identifying and targeting molecular mechanisms involved in autophagy is a potential therapeutic strategy for human cancers. Multiple transcription factors (TFs) are involved in the regulation of autophagy by modulating the expression of autophagy-related genes (ATGs). In addition, a wide array of long noncoding RNAs (lncRNAs), a group of regulatory ncRNAs, have been evidenced to regulate the function of these autophagy-related TFs through tumorigenesis. Subsequently, the lncRNAs/TFs/ATGs axis shows great potential as a therapeutic target for human cancers. Therefore, this review aimed to summarize new findings about the role of lncRNAs in regulating autophagy-related TFs with therapeutic perspectives.

Identification of ancestral gnathostome Gli3 enhancers with activity in mammals

Abnormal expression of the transcriptional regulator and hedgehog (Hh) signaling pathway effector Gli3 is known to trigger congenital disease, most frequently affecting the central nervous system (CNS) and the limbs. Accurate delineation of the genomic cis-regulatory landscape controlling Gli3 transcription during embryonic development is critical for the interpretation of noncoding variants associated with congenital defects. Here, we employed a comparative genomic analysis on fish species with a slow rate of molecular evolution to identify seven previously unknown conserved noncoding elements (CNEs) in Gli3 intronic intervals (CNE15-21). Transgenic assays in zebrafish revealed that most of these elements drive activities in Gli3 expressing tissues, predominantly the fins, CNS, and the heart. Intersection of these CNEs with human disease associated SNPs identified CNE15 as a putative mammalian craniofacial enhancer, with conserved activity in vertebrates and potentially affected by mutation associated with human craniofacial morphology. Finally, comparative functional dissection of an appendage-specific CNE conserved in slowly evolving fish (elephant shark), but not in teleost (CNE14/hs1586) indicates co-option of limb specificity from other tissues prior to the divergence of amniotes and lobe-finned fish. These results uncover a novel subset of intronic Gli3 enhancers that arose in the common ancestor of gnathostomes and whose sequence components were likely gradually modified in other species during the process of evolutionary diversification.

Role of Genetic Variation in Transcriptional Regulatory Elements in Heart Rhythm

Genetic predisposition to cardiac arrhythmias has been a field of intense investigation. Research initially focused on rare hereditary arrhythmias, but over the last two decades, the role of genetic variation (single nucleotide polymorphisms) in heart rate, rhythm, and arrhythmias has been taken into consideration as well. In particular, genome-wide association studies have identified hundreds of genomic loci associated with quantitative electrocardiographic traits, atrial fibrillation, and less common arrhythmias such as Brugada syndrome. A significant number of associated variants have been found to systematically localize in non-coding regulatory elements that control the tissue-specific and temporal transcription of genes encoding transcription factors, ion channels, and other proteins. However, the identification of causal variants and the mechanism underlying their impact on phenotype has proven difficult due to the complex tissue-specific, time-resolved, condition-dependent, and combinatorial function of regulatory elements, as well as their modest conservation across different model species. In this review, we discuss research efforts aimed at identifying and characterizing-trait-associated variant regulatory elements and the molecular mechanisms underlying their impact on heart rate or rhythm.

Rapid and Quantitative Functional Interrogation of Human Enhancer Variant Activity in Live Mice

Functional analysis of non-coding variants associated with human congenital disorders remains challenging due to the lack of efficient in vivo models. Here we introduce dual-enSERT, a robust Cas9-based two-color fluorescent reporter system which enables rapid, quantitative comparison of enhancer allele activities in live mice of any genetic background. We use this new technology to examine and measure the gain- and loss-of-function effects of enhancer variants linked to limb polydactyly, autism, and craniofacial malformation. By combining dual-enSERT with single-cell transcriptomics, we characterize variant enhancer alleles at cellular resolution, thereby implicating candidate molecular pathways in pathogenic enhancer misregulation. We further show that independent, polydactyly-linked enhancer variants lead to ectopic expression in the same cell populations, indicating shared genetic mechanisms underlying non-coding variant pathogenesis. Finally, we streamline dual-enSERT for analysis in F0 animals by placing both reporters on the same transgene separated by a synthetic insulator. Dual-enSERT allows researchers to go from identifying candidate enhancer variants to analysis of comparative enhancer activity in live embryos in under two weeks.

Allele biased transcription factor binding across human brain regions gives mechanistic insight into eQTLs

Transcription Factors (TFs) influence gene expression by facilitating or disrupting the formation of transcription initiation machinery at particular genomic loci. Because genomic localization of TFs is in part driven by TF recognition of DNA sequence, variation in TF binding sites can disrupt TF-DNA associations and affect gene regulation. To identify variants that impact TF binding in human brain tissues, we quantified allele bias for 93 TFs analyzed with ChIP-seq experiments of multiple structural brain regions from two donors. Using graph genomes constructed from phased genomic sequence data, we compared ChIP-seq signal between alleles at heterozygous variants within each tissue sample from each donor. Comparison of results from different brain regions within donors and the same regions between donors provided measures of allele bias reproducibility. We identified thousands of DNA variants that show reproducible bias in ChIP-seq for at least one TF. We found that alleles that are rarer in the general population were more likely than common alleles to exhibit large biases, and more frequently led to reduced TF binding. Combining ChIP-seq with RNA-seq, we identified TF-allele interaction biases with RNA bias in a phased allele linked to 6,709 eQTL variants identified in GTEx data, 3,309 of which were found in neural contexts. Our results provide insights into the effects of both common and rare variation on gene regulation in the brain. These findings can facilitate mechanistic understanding of cis-regulatory variation associated with biological traits, including disease.

A Kuhnian revolution in molecular biology: Most genes in complex organisms express regulatory RNAs

Thomas Kuhn described the progress of science as comprising occasional paradigm shifts separated by interludes of 'normal science'. The paradigm that has held sway since the inception of molecular biology is that genes (mainly) encode proteins. In parallel, theoreticians posited that mutation is random, inferred that most of the genome in complex organisms is non-functional, and asserted that somatic information is not communicated to the germline. However, many anomalies appeared, particularly in plants and animals: the strange genetic phenomena of paramutation and transvection; introns; repetitive sequences; a complex epigenome; lack of scaling of (protein-coding) genes and increase in 'noncoding' sequences with developmental complexity; genetic loci termed 'enhancers' that control spatiotemporal gene expression patterns during development; and a plethora of 'intergenic', overlapping, antisense and intronic transcripts. These observations suggest that the original conception of genetic information was deficient and that most genes in complex organisms specify regulatory RNAs, some of which convey intergenerational information. Also see the video abstract here: https://youtu.be/qxeGwahBANw.

Neural crest cells as a source of microevolutionary variation

Vertebrates have some of the most complex and diverse features in animals, from varied craniofacial morphologies to colorful pigmentation patterns and elaborate social behaviors. All of these traits have their developmental origins in a multipotent embryonic lineage of neural crest cells. This "fourth germ layer" is a vertebrate innovation and the source of a wide range of adult cell types. While others have discussed the role of neural crest cells in human disease and animal domestication, less is known about their role in contributing to adaptive changes in wild populations. Here, we review how variation in the development of neural crest cells and their derivatives generates considerable phenotypic diversity in nature. We focus on the broad span of traits under natural and sexual selection whose variation may originate in the neural crest, with emphasis on behavioral factors such as intraspecies communication that are often overlooked. In all, we encourage the integration of evolutionary ecology with developmental biology and molecular genetics to gain a more complete understanding of the role of this single cell type in trait covariation, evolutionary trajectories, and vertebrate diversity.

Cis-regulatory landscapes in the evolution and development of the mammalian skull

Extensive morphological variation found in mammals reflects the wide spectrum of their ecological adaptations. The highest morphological diversity is present in the craniofacial region, where geometry is mainly dictated by the bony skull. Mammalian craniofacial development represents complex multistep processes governed by numerous conserved genes that require precise spatio-temporal control. A central question in contemporary evolutionary biology is how a defined set of conserved genes can orchestrate formation of fundamentally different structures, and therefore how morphological variability arises. In principle, differential gene expression patterns during development are the source of morphological variation. With the emergence of multicellular organisms, precise regulation of gene expression in time and space is attributed to cis-regulatory elements. These elements contribute to higher-order chromatin structure and together with trans-acting factors control transcriptional landscapes that underlie intricate morphogenetic processes. Consequently, divergence in cis-regulation is believed to rewire existing gene regulatory networks and form the core of morphological evolution. This review outlines the fundamental principles of the genetic code and genomic regulation interplay during development. Recent work that deepened our comprehension of cis-regulatory element origin, divergence and function is presented here to illustrate the state-of-the-art research that uncovered the principles of morphological novelty. This article is part of the theme issue 'The mammalian skull: development, structure and function'.

Shadow enhancers mediate trade-offs between transcriptional noise and fidelity

Enhancers are stretches of regulatory DNA that bind transcription factors (TFs) and regulate the expression of a target gene. Shadow enhancers are two or more enhancers that regulate the same target gene in space and time and are associated with most animal developmental genes. These multi-enhancer systems can drive more consistent transcription than single enhancer systems. Nevertheless, it remains unclear why shadow enhancer TF binding sites are distributed across multiple enhancers rather than within a single large enhancer. Here, we use a computational approach to study systems with varying numbers of TF binding sites and enhancers. We employ chemical reaction networks with stochastic dynamics to determine the trends in transcriptional noise and fidelity, two key performance objectives of enhancers. This reveals that while additive shadow enhancers do not differ in noise and fidelity from their single enhancer counterparts, sub- and superadditive shadow enhancers have noise and fidelity trade-offs not available to single enhancers. We also use our computational approach to compare the duplication and splitting of a single enhancer as mechanisms for the generation of shadow enhancers and find that the duplication of enhancers can decrease noise and increase fidelity, although at the metabolic cost of increased RNA production. A saturation mechanism for enhancer interactions similarly improves on both of these metrics. Taken together, this work highlights that shadow enhancer systems may exist for several reasons: genetic drift or the tuning of key functions of enhancers, including transcription fidelity, noise and output.

CpG island turnover events predict evolutionary changes in enhancer activity

Genetic changes that modify the function of transcriptional enhancers have been linked to the evolution of biological diversity across species. Multiple studies have focused on the role of nucleotide substitutions, transposition, and insertions and deletions in altering enhancer function. Here we show that turnover of CpG islands (CGIs), which contribute to enhancer activation, is broadly associated with changes in enhancer activity across mammals, including humans. We integrated maps of CGIs and enhancer activity-associated histone modifications obtained from multiple tissues in nine mammalian species and found that CGI content in enhancers was strongly associated with increased histone modification levels. CGIs showed widespread turnover across species and species-specific CGIs were strongly enriched for enhancers exhibiting species-specific activity across all tissues and species we examined. Genes associated with enhancers with species-specific CGIs showed concordant biases in their expression, supporting that CGI turnover contributes to gene regulatory innovation. Our results also implicate CGI turnover in the evolution of Human Gain Enhancers (HGEs), which show increased activity in human embryonic development and may have contributed to the evolution of uniquely human traits. Using a humanized mouse model, we show that a highly conserved HGE with a large CGI absent from the mouse ortholog shows increased activity at the human CGI in the humanized mouse diencephalon. Collectively, our results point to CGI turnover as a mechanism driving gene regulatory changes potentially underlying trait evolution in mammals.

A cis-regulatory sequence of the selector gene vestigial drives the evolution of wing scaling in Drosophila species

Scaling between specific organs and overall body size has long fascinated biologists, being a primary mechanism by which organ shapes evolve. Yet, the genetic mechanisms that underlie the evolution of scaling relationships remain elusive. Here, we compared wing and fore tibia lengths (the latter as a proxy of body size) in Drosophila melanogaster, Drosophila simulans, Drosophila ananassae and Drosophila virilis, and show that the first three of these species have roughly a similar wing-to-tibia scaling behavior. In contrast, D. virilis exhibits much smaller wings relative to their body size compared with the other species and this is reflected in the intercept of the wing-to-tibia allometry. We then asked whether the evolution of this relationship could be explained by changes in a specific cis-regulatory region or enhancer that drives expression of the wing selector gene, vestigial (vg), whose function is broadly conserved in insects and contributes to wing size. To test this hypothesis directly, we used CRISPR/Cas9 to replace the DNA sequence of the predicted Quadrant Enhancer (vgQE) from D. virilis for the corresponding vgQE sequence in the genome of D. melanogaster. Strikingly, we discovered that D. melanogaster flies carrying the D. virilis vgQE sequence have wings that are significantly smaller with respect to controls, partially shifting the intercept of the wing-to-tibia scaling relationship towards that observed in D. virilis. We conclude that a single cis-regulatory element in D. virilis contributes to constraining wing size in this species, supporting the hypothesis that scaling could evolve through genetic variations in cis-regulatory elements.

Evolution of modular and pleiotropic enhancers

Cis-regulatory elements (CREs), or enhancers, are segments of noncoding DNA that regulate the spatial and temporal expression of nearby genes. Sometimes, genes are expressed in more than one tissue, and this can be driven by two main types of CREs: tissue-specific "modular" CREs, where different CREs drive expression of the gene in the different tissues, or by "pleiotropic" CREs, where the same CRE drives expression in the different tissues. In this perspective, we will discuss some of the ways (i) modular and pleiotropic CREs might originate; (ii) propose that modular CREs might derive from pleiotropic CREs via a process of duplication, degeneration, and complementation (the CRE-DDC model); and (iii) propose that hotspot loci of evolution are associated with the origin of modular CREs belonging to any gene in a regulatory network.

Resolving between novelty and homology in the rapidly evolving phallus of Drosophila

The genitalia present some of the most rapidly evolving anatomical structures in the animal kingdom, possessing a variety of parts that can distinguish recently diverged species. In the Drosophila melanogaster group, the phallus is adorned with several processes, pointed outgrowths, that are similar in size and shape between species. However, the complex three-dimensional nature of the phallus can obscure the exact connection points of each process. Previous descriptions based upon adult morphology have primarily assigned phallic processes by their approximate positions in the phallus and have remained largely agnostic regarding their homology relationships. In the absence of clearly identified homology, it can be challenging to model when each structure first evolved. Here, we employ a comparative developmental analysis of these processes in eight members of the melanogaster species group to precisely identify the tissue from which each process forms. Our results indicate that adult phallic processes arise from three pupal primordia in all species. We found that in some cases the same primordia generate homologous structures whereas in other cases, different primordia produce phenotypically similar but remarkably non-homologous structures. This suggests that the same gene regulatory network may have been redeployed to different primordia to induce phenotypically similar traits. Our results highlight how traits diversify and can be redeployed, even at short evolutionary scales.

Leveraging massively parallel reporter assays for evolutionary questions

A long-standing goal of evolutionary biology is to decode how gene regulation contributes to organismal diversity. Doing so is challenging because it is hard to predict function from non-coding sequence and to perform molecular research with non-model taxa. Massively parallel reporter assays (MPRAs) enable the testing of thousands to millions of sequences for regulatory activity simultaneously. Here, we discuss the execution, advantages, and limitations of MPRAs, with a focus on evolutionary questions. We propose solutions for extending MPRAs to rare taxa and those with limited genomic resources, and we underscore MPRA's broad potential for driving genome-scale, functional studies across organisms.

Deep cis-regulatory homology of the butterfly wing pattern ground plan

Butterfly wing patterns derive from a deeply conserved developmental ground plan yet are diverse and evolve rapidly. It is poorly understood how gene regulatory architectures can accommodate both deep homology and adaptive change. To address this, we characterized the cis-regulatory evolution of the color pattern gene WntA in nymphalid butterflies. Comparative assay for transposase-accessible chromatin using sequencing (ATAC-seq) and in vivo deletions spanning 46 cis-regulatory elements across five species revealed deep homology of ground plan-determining sequences, except in monarch butterflies. Furthermore, noncoding deletions displayed both positive and negative regulatory effects that were often broad in nature. Our results provide little support for models predicting rapid enhancer turnover and suggest that deeply ancestral, multifunctional noncoding elements can underlie rapidly evolving trait systems.

Site pleiotropy of a stickleback Bmp6 enhancer

Development and regeneration are orchestrated by gene regulatory networks that operate in part through transcriptional enhancers. Although many enhancers are pleiotropic and are active in multiple tissues, little is known about whether enhancer pleiotropy is due to 1) site pleiotropy, in which individual transcription factor binding sites (TFBS) are required for activity in multiple tissues, or 2) multiple distinct sites that regulate expression in different tissues. Here, we investigated the pleiotropy of an intronic enhancer of the stickleback Bone morphogenetic protein 6 (Bmp6) gene. This enhancer was previously shown to regulate evolved changes in tooth number and tooth regeneration, and is highly pleiotropic, with robust activity in both fins and teeth throughout embryonic, larval, and adult life, and in the heart and kidney in adult fish. We tested the hypothesis that the pleiotropy of this enhancer is due to site pleiotropy of an evolutionarily conserved predicted Foxc1 TFBS. Transgenic analysis and site-directed mutagenesis experiments both deleting and scrambling this predicted Foxc1 TFBS revealed that the binding site is required for enhancer activity in both teeth and fins throughout embryonic, larval, and adult development, and in the heart and kidney in adult fish. Collectively these data support a model where the pleiotropy of this Bmp6 enhancer is due to site pleiotropy and this putative binding site is required for enhancer activity in multiple anatomical sites from the embryo to the adult.

Species-specific sensitivity to TGFβ signaling and changes to the Mmp13 promoter underlie avian jaw development and evolution

Precise developmental control of jaw length is critical for survival, but underlying molecular mechanisms remain poorly understood. The jaw skeleton arises from neural crest mesenchyme (NCM), and we previously demonstrated that these progenitor cells express more bone-resorbing enzymes including Matrix metalloproteinase 13 (Mmp13) when they generate shorter jaws in quail embryos versus longer jaws in duck. Moreover, if we inhibit bone resorption or Mmp13, we can increase jaw length. In the current study, we uncover mechanisms establishing species-specific levels of Mmp13 and bone resorption. Quail show greater activation of and sensitivity to transforming growth factor beta (TGFβ) signaling than duck; where intracellular mediators like SMADs and targets like Runt-related transcription factor 2 (Runx2), which bind Mmp13, become elevated. Inhibiting TGFβ signaling decreases bone resorption, and overexpressing Mmp13 in NCM shortens the duck lower jaw. To elucidate the basis for this differential regulation, we examine the Mmp13 promoter. We discover a SMAD-binding element and single nucleotide polymorphisms (SNPs) near a RUNX2-binding element that distinguish quail from duck. Altering the SMAD site and switching the SNPs abolish TGFβ sensitivity in the quail Mmp13 promoter but make the duck promoter responsive. Thus, differential regulation of TGFβ signaling and Mmp13 promoter structure underlie avian jaw development and evolution.

Snake venom gene expression is coordinated by novel regulatory architecture and the integration of multiple co-opted vertebrate pathways

Understanding how regulatory mechanisms evolve is critical for understanding the processes that give rise to novel phenotypes. Snake venom systems represent a valuable and tractable model for testing hypotheses related to the evolution of novel regulatory networks, yet the regulatory mechanisms underlying venom production remain poorly understood. Here, we use functional genomics approaches to investigate venom regulatory architecture in the prairie rattlesnake and identify cis-regulatory sequences (enhancers and promoters), trans-regulatory transcription factors, and integrated signaling cascades involved in the regulation of snake venom genes. We find evidence that two conserved vertebrate pathways, the extracellular signal-regulated kinase and unfolded protein response pathways, were co-opted to regulate snake venom. In one large venom gene family (snake venom serine proteases), this co-option was likely facilitated by the activity of transposable elements. Patterns of snake venom gene enhancer conservation, in some cases spanning 50 million yr of lineage divergence, highlight early origins and subsequent lineage-specific adaptations that have accompanied the evolution of venom regulatory architecture. We also identify features of chromatin structure involved in venom regulation, including topologically associated domains and CTCF loops that underscore the potential importance of novel chromatin structure to coevolve when duplicated genes evolve new regulatory control. Our findings provide a model for understanding how novel regulatory systems may evolve through a combination of genomic processes, including tandem duplication of genes and regulatory sequences, cis-regulatory sequence seeding by transposable elements, and diverse transcriptional regulatory proteins controlled by a co-opted regulatory cascade.

The roles of inter-tissue adhesion in development and morphological evolution

The study of how neighboring tissues physically interact with each other, inter-tissue adhesion, is an emerging field at the interface of cell biology, biophysics and developmental biology. Inter-tissue adhesion can be mediated by either cell-extracellular matrix adhesion or cell-cell adhesion, and both the mechanisms and consequences of inter-tissue adhesion have been studied in vivo in numerous vertebrate and invertebrate species. In this Review, we discuss recent progress in understanding the many functions of inter-tissue adhesion in development and evolution. Inter-tissue adhesion can couple the motion of adjacent tissues, be the source of mechanical resistance that constrains morphogenesis, and transmit tension required for normal development. Tissue-tissue adhesion can also create mechanical instability that leads to tissue folding or looping. Transient inter-tissue adhesion can facilitate tissue invasion, and weak tissue adhesion can generate friction that shapes and positions tissues within the embryo. Lastly, we review studies that reveal how inter-tissue adhesion contributes to the diversification of animal morphologies.

Sex-specific evolution of a Drosophila sensory system via interacting cis- and trans-regulatory changes

The evolution of gene expression via cis-regulatory changes is well established as a major driver of phenotypic evolution. However, relatively little is known about the influence of enhancer architecture and intergenic interactions on regulatory evolution. We address this question by examining chemosensory system evolution in Drosophila. Drosophila prolongata males show a massively increased number of chemosensory bristles compared to females and males of sibling species. This increase is driven by sex-specific transformation of ancestrally mechanosensory organs. Consistent with this phenotype, the Pox neuro transcription factor (Poxn), which specifies chemosensory bristle identity, shows expanded expression in D. prolongata males. Poxn expression is controlled by nonadditive interactions among widely dispersed enhancers. Although some D. prolongata Poxn enhancers show increased activity, the additive component of this increase is slight, suggesting that most changes in Poxn expression are due to epistatic interactions between Poxn enhancers and trans-regulatory factors. Indeed, the expansion of D. prolongata Poxn enhancer activity is only observed in cells that express doublesex (dsx), the gene that controls sexual differentiation in Drosophila and also shows increased expression in D. prolongata males due to cis-regulatory changes. Although expanded dsx expression may contribute to increased activity of D. prolongata Poxn enhancers, this interaction is not sufficient to explain the full expansion of Poxn expression, suggesting that cis-trans interactions between Poxn, dsx, and additional unknown genes are necessary to produce the derived D. prolongata phenotype. Overall, our results demonstrate the importance of epistatic gene interactions for evolution, particularly when pivotal genes have complex regulatory architecture.

Architecture and evolution of the cis-regulatory system of the echinoderm kirrelL gene

The gene regulatory network (GRN) that underlies echinoderm skeletogenesis is a prominent model of GRN architecture and evolution. KirrelL is an essential downstream effector gene in this network and encodes an Ig-superfamily protein required for the fusion of skeletogenic cells and the formation of the skeleton. In this study, we dissected the transcriptional control region of the kirrelL gene of the purple sea urchin, Strongylocentrotus purpuratus. Using plasmid- and bacterial artificial chromosome-based transgenic reporter assays, we identified key cis-regulatory elements (CREs) and transcription factor inputs that regulate Sp-kirrelL, including direct, positive inputs from two key transcription factors in the skeletogenic GRN, Alx1 and Ets1. We next identified kirrelL cis-regulatory regions from seven other echinoderm species that together represent all classes within the phylum. By introducing these heterologous regulatory regions into developing sea urchin embryos we provide evidence of their remarkable conservation across ~500 million years of evolution. We dissected in detail the kirrelL regulatory region of the sea star, Patiria miniata, and demonstrated that it also receives direct inputs from Alx1 and Ets1. Our findings identify kirrelL as a component of the ancestral echinoderm skeletogenic GRN. They support the view that GRN subcircuits, including specific transcription factor–CRE interactions, can remain stable over vast periods of evolutionary history. Lastly, our analysis of kirrelL establishes direct linkages between a developmental GRN and an effector gene that controls a key morphogenetic cell behavior, cell–cell fusion, providing a paradigm for extending the explanatory power of GRNs.

Co-option of the limb patterning program in cephalopod eye development

Background

Across the Metazoa, similar genetic programs are found in the development of analogous, independently evolved, morphological features. The functional significance of this reuse and the underlying mechanisms of co-option remain unclear. Cephalopods have evolved a highly acute visual system with a cup-shaped retina and a novel refractive lens in the anterior, important for a number of sophisticated behaviors including predation, mating, and camouflage. Almost nothing is known about the molecular-genetics of lens development in the cephalopod.

Results

Here we identify the co-option of the canonical bilaterian limb patterning program during cephalopod lens development, a functionally unrelated structure. We show radial expression of transcription factors SP6-9/sp1, Dlx/dll, Pbx/exd, Meis/hth, and a Prdl homolog in the squid Doryteuthis pealeii, similar to expression required in Drosophila limb development. We assess the role of Wnt signaling in the cephalopod lens, a positive regulator in the developing Drosophila limb, and find the regulatory relationship reversed, with ectopic Wnt signaling leading to lens loss.

Conclusion

This regulatory divergence suggests that duplication of SP6-9 in cephalopods may mediate the co-option of the limb patterning program. Thus, our study suggests that this program could perform a more universal developmental function in radial patterning and highlights how canonical genetic programs are repurposed in novel structures.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01182-2.

Evolved Bmp6 enhancer alleles drive spatial shifts in gene expression during tooth development in sticklebacks

Mutations in enhancers have been shown to often underlie natural variation but the evolved differences in enhancer activity can be difficult to identify in vivo. Threespine sticklebacks (Gasterosteus aculeatus) are a robust system for studying enhancer evolution due to abundant natural genetic variation, a diversity of evolved phenotypes between ancestral marine and derived freshwater forms, and the tractability of transgenic techniques. Previous work identified a series of polymorphisms within an intronic enhancer of the Bone morphogenetic protein 6 (Bmp6) gene that are associated with evolved tooth gain, a derived increase in freshwater tooth number that arises late in development. Here, we use a bicistronic reporter construct containing a genetic insulator and a pair of reciprocal two-color transgenic reporter lines to compare enhancer activity of marine and freshwater alleles of this enhancer. In older fish, the two alleles drive partially overlapping expression in both mesenchyme and epithelium of developing teeth, but the freshwater enhancer drives a reduced mesenchymal domain and a larger epithelial domain relative to the marine enhancer. In younger fish, these spatial shifts in enhancer activity are less pronounced. Comparing Bmp6 expression by in situ hybridization in developing teeth of marine and freshwater fish reveals similar evolved spatial shifts in gene expression. Together, these data support a model in which the polymorphisms within this enhancer underlie evolved tooth gain by shifting the spatial expression of Bmp6 during tooth development, and provide a general strategy to identify spatial differences in enhancer activity in vivo.

Repression precedes independent evolutionary gains of a highly specific gene expression pattern

Highly specific expression patterns can be caused by the overlapping activities of activator and repressor sequences in enhancers. However, few studies illuminate how these sequences evolve in the origin of new enhancers. Here, we show that expression of the bond gene in the semicircular wall epithelium (swe) of the Drosophila melanogaster male ejaculatory bulb (EB) is controlled by an enhancer consisting of an activator region that requires Abdominal-B driving expression in the entire EB and a repressor region that restricts this expression to the EB swe. Although this expression pattern is independently gained in the distantly related Scaptodrosophila lebanonensis and does not require Abdominal-B, we show that functionally similar repressor sequences are present in Scaptodrosophila and also in species that do not express bond in the EB. We suggest that during enhancer evolution, repressor sequences can precede the evolution of activator sequences and may lead to similar but independently evolved expression patterns.

Pu et al. show that the independent gain of a highly specific expression pattern across distantly related species may be because of the preexistence of repressor sequences that precedes the diversification of these species. This may reflect a general mechanism underlying the evolution of highly specific enhancers.

Evolution of mouse circadian enhancers from transposable elements

BACKGROUND

Transposable elements are increasingly recognized as a source of cis-regulatory variation. Previous studies have revealed that transposons are often bound by transcription factors and some have been co-opted into functional enhancers regulating host gene expression. However, the process by which transposons mature into complex regulatory elements, like enhancers, remains poorly understood. To investigate this process, we examined the contribution of transposons to the cis-regulatory network controlling circadian gene expression in the mouse liver, a well-characterized network serving an important physiological function.

RESULTS

ChIP-seq analyses reveal that transposons and other repeats contribute ~ 14% of the binding sites for core circadian regulators (CRs) including BMAL1, CLOCK, PER1/2, and CRY1/2, in the mouse liver. RSINE1, an abundant murine-specific SINE, is the only transposon family enriched for CR binding sites across all datasets. Sequence analyses and reporter assays reveal that the circadian regulatory activity of RSINE1 stems from the presence of imperfect CR binding motifs in the ancestral RSINE1 sequence. These motifs matured into canonical motifs through point mutations after transposition. Furthermore, maturation occurred preferentially within elements inserted in the proximity of ancestral CR binding sites. RSINE1 also acquired motifs that recruit nuclear receptors known to cooperate with CRs to regulate circadian gene expression specifically in the liver.

CONCLUSIONS

Our results suggest that the birth of enhancers from transposons is predicated both by the sequence of the transposon and by the cis-regulatory landscape surrounding their genomic integration site.

Deploying Big Data to Crack the Genotype to Phenotype Code

Mechanistically connecting genotypes to phenotypes is a longstanding and central mission of biology. Deciphering these connections will unite questions and datasets across all scales from molecules to ecosystems. Although high-throughput sequencing has provided a rich platform on which to launch this effort, tools for deciphering mechanisms further along the genome to phenome pipeline remain limited. Machine learning approaches and other emerging computational tools hold the promise of augmenting human efforts to overcome these obstacles. This vision paper is the result of a Reintegrating Biology Workshop, bringing together the perspectives of integrative and comparative biologists to survey challenges and opportunities in cracking the genotype to phenotype code and thereby generating predictive frameworks across biological scales. Key recommendations include promoting the development of minimum "best practices" for the experimental design and collection of data; fostering sustained and long-term data repositories; promoting programs that recruit, train, and retain a diversity of talent; and providing funding to effectively support these highly cross-disciplinary efforts. We follow this discussion by highlighting a few specific transformative research opportunities that will be advanced by these efforts.

A cis-regulatory atlas in maize at single-cell resolution

cis-regulatory elements (CREs) encode the genomic blueprints of spatiotemporal gene expression programs enabling highly specialized cell functions. Using single-cell genomics in six maize organs, we determined the cis- and trans-regulatory factors defining diverse cell identities and coordinating chromatin organization by profiling transcription factor (TF) combinatorics, identifying TFs with non-cell-autonomous activity, and uncovering TFs underlying higher-order chromatin interactions. Cell-type-specific CREs were enriched for enhancer activity and within unmethylated long terminal repeat retrotransposons. Moreover, we found cell-type-specific CREs are hotspots for phenotype-associated genetic variants and were targeted by selection during modern maize breeding, highlighting the biological implications of this CRE atlas. Through comparison of maize and Arabidopsis thaliana developmental trajectories, we identified TFs and CREs with conserved and divergent chromatin dynamics, showcasing extensive evolution of gene regulatory networks. In addition to this rich dataset, we developed single-cell analysis software, Socrates, which can be used to understand cis-regulatory variation in any species.

Evolutionary divergence of a Hoxa2b hindbrain enhancer in syngnathids mimics results of functional assays

Hoxa2 genes provide critical patterning signals during development, and their regulation and function have been extensively studied. We report a previously uncharacterized significant sequence divergence of a highly conserved hindbrain hoxa2b enhancer element in the family syngnathidae (pipefishes, seahorses, pipehorses, seadragons). We compared the hox cis-regulatory element variation in the Gulf pipefish and two species of seahorse against eight other species of fish, as well as human and mouse. We annotated the hoxa2b enhancer element binding sites across three species of seahorse, four species of pipefish, and one species of ghost pipefish. Finally, we performed in situ hybridization analysis of hoxa2b expression in Gulf pipefish embryos. We found that all syngnathid fish examined share a modified rhombomere 4 hoxa2b enhancer element, despite the fact that this element has been found to be highly conserved across all vertebrates examined previously. Binding element sequence motifs and spacing between binding elements have been modified for the hoxa2b enhancer in several species of pipefish and seahorse, and that the loss of the Prep/Meis binding site and further space shortening happened after ghost pipefish split from the rest of the syngnathid clade. We showed that expression of this gene in rhombomere 4 is lower relative to the surrounding rhombomeres in developing Gulf pipefish embryos, reflecting previously published functional tests for this enhancer. Our findings highlight the benefits of studying highly derived, diverse taxa for understanding of gene regulatory evolution and support the hypothesis that natural mutations can occur in deeply conserved pathways in ways potentially related to phenotypic diversity.

A cis-regulatory atlas in maize at single-cell resolution

cis-regulatory elements (CREs) encode the genomic blueprints of spatiotemporal gene expression programs enabling highly specialized cell functions. Using single-cell genomics in six maize organs, we determined the cis- and trans-regulatory factors defining diverse cell identities and coordinating chromatin organization by profiling transcription factor (TF) combinatorics, identifying TFs with non-cell-autonomous activity, and uncovering TFs underlying higher-order chromatin interactions. Cell-type-specific CREs were enriched for enhancer activity and within unmethylated long terminal repeat retrotransposons. Moreover, we found cell-type-specific CREs are hotspots for phenotype-associated genetic variants and were targeted by selection during modern maize breeding, highlighting the biological implications of this CRE atlas. Through comparison of maize and Arabidopsis thaliana developmental trajectories, we identified TFs and CREs with conserved and divergent chromatin dynamics, showcasing extensive evolution of gene regulatory networks. In addition to this rich dataset, we developed single-cell analysis software, Socrates, which can be used to understand cis-regulatory variation in any species.

How to build a crop plant: Defining the cis-regulatory landscape of maize

The functional regulatory elements of agronomically important plant genomes have been long sought after. Marand et. al. generate a comprehensive atlas of cis-regulatory elements at single cell resolution in maize, providing a powerful resource for inquiries into the rules of multicellular development and for precision crop engineering.

Conserved pleiotropy of an ancient plant homeobox gene uncovered by cis-regulatory dissection

Divergence of gene function is a hallmark of evolution, but assessing functional divergence over deep time is not trivial. The few alleles available for cross-species studies often fail to expose the entire functional spectrum of genes, potentially obscuring deeply conserved pleiotropic roles. Here, we explore the functional divergence of WUSCHEL HOMEOBOX9 (WOX9), suggested to have species-specific roles in embryo and inflorescence development. Using a cis-regulatory editing drive system, we generate a comprehensive allelic series in tomato, which revealed hidden pleiotropic roles for WOX9. Analysis of accessible chromatin and conserved cis-regulatory sequences identifies the regions responsible for this pleiotropic activity, the functions of which are conserved in groundcherry, a tomato relative. Mimicking these alleles in Arabidopsis, distantly related to tomato and groundcherry, reveals new inflorescence phenotypes, exposing a deeply conserved pleiotropy. We suggest that targeted cis-regulatory mutations can uncover conserved gene functions and reduce undesirable effects in crop improvement.

Developmental hourglass and heterochronic shifts in fin and limb development

How genetic changes are linked to morphological novelties and developmental constraints remains elusive. Here, we investigate genetic apparatuses that distinguish fish fins from tetrapod limbs by analyzing transcriptomes and open-chromatin regions (OCRs). Specifically, we compared mouse forelimb buds with the pectoral fin buds of an elasmobranch, the brown-banded bamboo shark (Chiloscyllium punctatum). A transcriptomic comparison with an accurate orthology map revealed both a mass heterochrony and hourglass-shaped conservation of gene expression between fins and limbs. Furthermore, open-chromatin analysis suggested that access to conserved regulatory sequences is transiently increased during mid-stage limb development. During this stage, stage-specific and tissue-specific OCRs were also enriched. Together, early and late stages of fin/limb development are more permissive to mutations than middle stages, which may have contributed to major morphological changes during the fin-to-limb evolution. We hypothesize that the middle stages are constrained by regulatory complexity that results from dynamic and tissue-specific transcriptional controls.

Animals come in all shapes and sizes. This diversity arose through genetic mutations during evolution, but it is unclear exactly how these variations led to the formation of new shapes. There is increasing evidence to suggest that not all shapes are possible and that variability between animals is limited by a phenomenon known as “developmental constraint”. These limitations direct parts of the body towards a specific shape as they develop in the embryo. Therefore, understanding the mechanisms underlying these developmental constraints could help explain how different body shapes evolved.

The limbs of humans and other mammals evolved from the fins of fish, and this transition is often used to study the role developmental constraints play in evolution. This is an ideal model as there is already a detailed fossil record mapping this evolutionary event, and data pinpointing some of the genes involved in the development of limbs and fins. But this data is incomplete, and a full comparison between the genes activated in the fin and the limb during embryonic development had not been achieved. This is because most fish used for research have undergone recent genetic changes, making it hard to spot which genetic differences are linked to the evolution of the limb.

To overcome this barrier, Onimaru et al. compared genetic data from the developing limbs of mice to the developing fins of the brown-banded bamboo shark, which evolves much slower than other fish. This revealed that although many genes commonly played a role in the development of the fin and the limb in the embryo, the activity of these shared genes was not the same. For example, genes that switched on in the late stages of limb development, switched off in the late stages of fin development. But in the middle of development, those differences were relatively small and both species activated very similar sets of genes. Many of these genes were pleiotropic, which means they have important roles in other tissues and therefore mutate less often. This suggests that the mid-stage of limb development is under the strongest level of constraint.

Darwin’s theory of natural selection explains that mutations drive evolution. But the theory cannot predict what kinds of new body shapes new mutations will produce. Understanding how the activity levels of different genes affect development could help to fill this knowledge gap. This has potential medical applications, for example, understanding why some genetic changes cause more serious problems than others. This work suggests that mutations in genes that are active during the mid-stage of limb development may have the most serious impact.

Conservation of peripheral nervous system formation mechanisms in divergent ascidian embryos

Ascidians with very similar embryos but highly divergent genomes are thought to have undergone extensive developmental system drift. We compared, in four species (Ciona and Phallusia for Phlebobranchia, Molgula and Halocynthia for Stolidobranchia), gene expression and gene regulation for a network of six transcription factors regulating peripheral nervous system (PNS) formation in Ciona. All genes, but one in Molgula, were expressed in the PNS with some differences correlating with phylogenetic distance. Cross-species transgenesis indicated strong levels of conservation, except in Molgula, in gene regulation despite lack of sequence conservation of the enhancers. Developmental system drift in ascidians is thus higher for gene regulation than for gene expression and is impacted not only by phylogenetic distance, but also in a clade-specific manner and unevenly within a network. Finally, considering that Molgula is divergent in our analyses, this suggests deep conservation of developmental mechanisms in ascidians after 390 My of separate evolution.

Active turnover of DNA methylation during cell fate decisions

DNA methylation is a key layer of epigenetic regulation. The deposition of methylation marks relies on the catalytic activity of DNA methyltransferases (DNMTs), and their active removal relies on the activity of ten-eleven translocation (TET) enzymes. Paradoxically, in important biological contexts these antagonistic factors are co-expressed and target overlapping genomic regions. The ensuing cyclic biochemistry of cytosine modifications gives rise to a continuous, out-of-thermal equilibrium transition through different methylation states. But what is the purpose of this intriguing turnover of DNA methylation? Recent evidence demonstrates that methylation turnover is enriched at gene distal regulatory elements, including enhancers, and can give rise to large-scale oscillatory dynamics. We discuss this phenomenon and propose that DNA methylation turnover might facilitate key lineage decisions.

Adaptation of the carbamoyl-phosphate synthetase enzyme in an extremophile fish

Tetrapods and fish have adapted distinct carbamoyl-phosphate synthase (CPS) enzymes to initiate the ornithine urea cycle during the detoxification of nitrogenous wastes. We report evidence that in the ureotelic subgenus of extremophile fish Oreochromis Alcolapia, CPS III has undergone convergent evolution and adapted its substrate affinity to ammonia, which is typical of terrestrial vertebrate CPS I. Unusually, unlike in other vertebrates, the expression of CPS III in Alcolapia is localized to the skeletal muscle and is activated in the myogenic lineage during early embryonic development with expression remaining in mature fish. We propose that adaptation in Alcolapia included both convergent evolution of CPS function to that of terrestrial vertebrates, as well as changes in development mechanisms redirecting CPS III gene expression to the skeletal muscle.

Enhancer evolutionary co-option through shared chromatin accessibility input

Form diversity is fueled by changes in the expression of genes that build organisms. New expression often results from the emergence of new DNA switches, known as transcriptional enhancers. Many enhancers are thought to appear through the recycling of older enhancers, a process called evolutionary co-option. Enhancer co-option is difficult to assess, and the molecular mechanisms explaining its prevalence are elusive. Using state-of-the-art quantification and analyses, we reveal that the sequences of an ancestral and a derived enhancer overlap extensively. They contain specific binding sites for regulators imparting spatial activities. We found that the two enhancers also share a site facilitating access to chromatin in a region where they overlap.

The diversity of forms in multicellular organisms originates largely from the spatial redeployment of developmental genes [S. B. Carroll, Cell 134, 25–36 (2008)]. Several scenarios can explain the emergence of cis-regulatory elements that govern novel aspects of a gene expression pattern [M. Rebeiz, M. Tsiantis, Curr. Opin. Genet. Dev. 45, 115–123 (2017)]. One scenario, enhancer co-option, holds that a DNA sequence producing an ancestral regulatory activity also becomes the template for a new regulatory activity, sharing regulatory information. While enhancer co-option might fuel morphological diversification, it has rarely been documented [W. J. Glassford et al., Dev. Cell 34, 520–531 (2015)]. Moreover, if two regulatory activities are borne from the same sequence, their modularity, considered a defining feature of enhancers [J. Banerji, L. Olson, W. Schaffner, Cell 33, 729–740 (1983)], might be affected by pleiotropy. Sequence overlap may thereby play a determinant role in enhancer function and evolution. Here, we investigated this problem with two regulatory activities of the Drosophila gene yellow, the novel spot enhancer and the ancestral wing blade enhancer. We used precise and comprehensive quantification of each activity in Drosophila wings to systematically map their sequences along the locus. We show that the spot enhancer has co-opted the sequences of the wing blade enhancer. We also identified a pleiotropic site necessary for DNA accessibility of a shared regulatory region. While the evolutionary steps leading to the derived activity are still unknown, such pleiotropy suggests that enhancer accessibility could be one of the molecular mechanisms seeding evolutionary co-option.

Evaluating Enhancer Function and Transcription

Cell-type- and condition-specific profiles of gene expression require coordination between protein-coding gene promoters and cis-regulatory sequences called enhancers. Enhancers can stimulate gene activity at great genomic distances from their targets, raising questions about how enhancers communicate with specific gene promoters and what molecular mechanisms underlie enhancer function. Characterization of enhancer loci has identified the molecular features of active enhancers that accompany the binding of transcription factors and local opening of chromatin. These characteristics include coactivator recruitment, histone modifications, and noncoding RNA transcription. However, it remains unclear which of these features functionally contribute to enhancer activity. Here, we discuss what is known about how enhancers regulate their target genes and how enhancers and promoters communicate. Further, we describe recent data demonstrating many similarities between enhancers and the gene promoters they control, and we highlight unanswered questions in the field, such as the potential roles of transcription at enhancers.

Temporal flexibility of gene regulatory network underlies a novel wing pattern in flies

Developmental genes can be coopted to generate evolutionary novelties by changing their spatial regulation. However, developmental genes seldom act independently, but rather work in a gene regulatory network (GRN). How is it possible to recruit a single gene from a whole GRN? What are the properties that allow parallel cooptions of the same genes during evolution? Here, we show that a novel engrailed gene expression underlies a novel wing color pattern in flies. We show that cooption is facilitated 1) because of GRN flexibility over development and 2) because every single gene of the GRN has its own functional time window. We suggest these two temporal properties could explain why the same gene can be independently recruited several times during evolution.

Organisms have evolved endless morphological, physiological, and behavioral novel traits during the course of evolution. Novel traits were proposed to evolve mainly by orchestration of preexisting genes. Over the past two decades, biologists have shown that cooption of gene regulatory networks (GRNs) indeed underlies numerous evolutionary novelties. However, very little is known about the actual GRN properties that allow such redeployment. Here we have investigated the generation and evolution of the complex wing pattern of the fly Samoaia leonensis. We show that the transcription factor Engrailed is recruited independently from the other players of the anterior–posterior specification network to generate a new wing pattern. We argue that partial cooption is made possible because 1) the anterior–posterior specification GRN is flexible over time in the developing wing and 2) this flexibility results from the fact that every single gene of the GRN possesses its own functional time window. We propose that the temporal flexibility of a GRN is a general prerequisite for its possible cooption during the course of evolution.

On the specificity of gene regulatory networks: How does network co-option affect subsequent evolution?

The process of multicellular organismal development hinges upon the specificity of developmental programs: for different parts of the organism to form unique features, processes must exist to specify each part. This specificity is thought to be hardwired into gene regulatory networks, which activate cohorts of genes in particular tissues at particular times during development. However, the evolution of gene regulatory networks sometimes occurs by mechanisms that sacrifice specificity. One such mechanism is network co-option, in which existing gene networks are redeployed in new developmental contexts. While network co-option may offer an efficient mechanism for generating novel phenotypes, losses of tissue specificity at redeployed network genes could restrict the ability of the affected traits to evolve independently. At present, there has not been a detailed discussion regarding how tissue specificity of network genes might be altered due to gene network co-option at its initiation, as well as how trait independence can be retained or restored after network co-option. A lack of clarity about network co-option makes it more difficult to speculate on the long-term evolutionary implications of this mechanism. In this review, we will discuss the possible initial outcomes of network co-option, outline the mechanisms by which networks may retain or subsequently regain specificity after network co-option, and comment on some of the possible evolutionary consequences of network co-option. We place special emphasis on the need to consider selectively-neutral outcomes of network co-option to improve our understanding of the role of this mechanism in trait evolution.

Disparate expression specificities coded by a shared Hox-C enhancer

Can a single regulatory sequence be shared by two genes undergoing functional divergence? Here we describe a single promiscuous enhancer within the Drosophila Antennapedia Complex, EO053, that directs aspects of the expression of two adjacent genes, pb (a Hox2 ortholog) and zen2 (a divergent Hox3 paralog), with disparate spatial and temporal expression patterns. We were unable to separate the pb-like and zen2-like specificities within EO053, and we identify sequences affecting both expression patterns. Importantly, genomic deletion experiments demonstrate that EO053 cooperates with additional pb- and zen2-specific enhancers to regulate the mRNA expression of both genes. We examine sequence conservation of EO053 within the Schizophora, and show that patterns of synteny between the Hox2 and Hox3 orthologs in Arthropods are consistent with a shared regulatory relationship extending prior to the Hox3/zen divergence. Thus, EO053 represents an example of two genes having evolved disparate outputs while utilizing this shared regulatory region.

Editorial note: This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).

Evolution of wing pigmentation in Drosophila: Diversity, physiological regulation, and cis-regulatory evolution

Fruit flies (Drosophila and its close relatives, or “drosophilids”) are a group that includes an important model organism, Drosophila melanogaster, and also very diverse species distributed worldwide. Many of these species have black or brown pigmentation patterns on their wings, and have been used as material for evo‐devo research. Pigmentation patterns are thought to have evolved rapidly compared with body plans or body shapes; hence they are advantageous model systems for studying evolutionary gains of traits and parallel evolution. Various groups of drosophilids, including genus Idiomyia (Hawaiian Drosophila), have a variety of pigmentations, ranging from simple black pigmentations around crossveins to a single antero‐distal spot and a more complex mottled pattern. Pigmentation patterns are sometimes obviously used for sexual displays; however, in some cases they may have other functions. The process of wing formation in Drosophila, the general mechanism of pigmentation formation, and the transport of substances necessary for pigmentation, including melanin precursors, through wing veins are summarized here. Lastly, the evolution of the expression of genes regulating pigmentation patterns, the role of cis‐regulatory regions, and the conditions required for the evolutionary emergence of pigmentation patterns are discussed. Future prospects for research on the evolution of wing pigmentation pattern formation in drosophilids are presented, particularly from the point of view of how they compare with other studies of the evolution of new traits.

Contrasted evolutionary trajectories of plant transcription factors

Because of their prominent roles in plant development, transcription factors (TF) play central roles as drivers of innovation in the evolution of the green lineage (viridiplantae). The advent of massive sequencing combined with comparative genetics/genomics allows a rigorous investigation of how TF families have contributed to plant diversification from charophyte algae to bryophytes to angiosperms. Here, we review recent progress on TF family reconstruction and the identification of distantly related TFs present throughout the evolutionary timeline from algae to angiosperms. These data provide examples of contrasting evolutionary trajectories of TF families and illustrate how conserved TFs adopt diverse roles over the course of evolution.

Taxonomic Notes on the 'Mahat' (Artocarpus lacucha and A. thailandicus, Moraceae) Species Complex in Thailand

‘Mahat’ is a well-known medicinal plant utilized in Thailand. The Thai name ‘Mahat’ has been used in many scientific articles for years. However, it is, unpredictably, a homonym of two scientific names in Flora of Thailand, i.e., A. lacucha and A. thailandicus. Additionally, both species are complex due to their high morphological variation. This causes difficulties in species identification especially when this Thai name is referred to as the scientific name for research publication, quality control of pharmaceutical raw materials, and registration of pharmaceutical products. In this study, we scrutinized the taxonomy of ‘Mahat’ by detailed examination of its morphology and distribution, including molecular and qualitative phytochemical studies. Leaf surfaces were inspected using scanning electron microscopy. The phylogeny of both species was studied using DNA sequences of nuclear and plastid regions. Chromatographic fingerprints, focusing on the major active compound oxyresveratrol, were identified using high-performance liquid chromatography. According to our current study, phylogenetic evidence showed that some samples of both species were clustered together in the same clade and phytochemical fingerprints were almost identical. These results are valuable data for taxonomic revision in the near future and reveal the possible utilization of A. thailandicus as a new material source of oxyresveratrol in the pharmaceutical industry.