Comprehensive Analysis of a cis-Regulatory Region Reveals Pleiotropy in Enhancer Function

Beyond modular enhancers: new questions in cis-regulatory evolution

Our understanding of how cis-regulatory elements work has advanced rapidly, outpacing our evolutionary models. In this review, we consider the implications of new mechanistic findings for evolutionary developmental biology. We focus on three different debates: whether evolutionary innovation occurs more often via the modification of old cis-regulatory elements or the emergence of new ones; the extent to which individual elements are specific and autonomous or multifunctional and interdependent; and how the robustness of cis-regulatory architectures influences the rate of trait evolution. These discussions lead us to propose new questions for the evo-devo of cis-regulation.

Notch-Dependent Expression of the Drosophila Hey Gene Is Supported by a Pair of Enhancers with Overlapping Activities

Drosophila Hey is a basic helix-loop-helix-orange (bHLH-O) protein with an important role in the establishment of distinct identities of postmitotic cells. We have previously identified Hey as a transcriptional target and effector of Notch signalling during the asymmetric division of neuronal progenitors, generating neurons of two types, and we have shown that Notch-dependent expression of Hey also marks a subpopulation of the newborn enteroendocrine (EE) cells in the midgut primordium of the embryo. Here, we investigate the transcriptional regulation of Hey in neuronal and intestinal tissues. We isolated two genomic regions upstream of the promoter (HeyUP) and in the second intron (HeyIN2) of the Hey gene, based on the presence of binding motifs for Su(H), the transcription factor that mediates Notch activity. We found that both regions can direct the overlapping expression patterns of reporter transgenes recapitulating endogenous Hey expression. Moreover, we showed that while HeyIN2 represents a Notch-dependent enhancer, HeyUP confers both Notch-dependent and independent transcriptional regulation. We induced mutations that removed the Su(H) binding motifs in either region and then studied the enhancer functionality in the respective Hey mutant lines. Our results provide direct evidence that although both enhancers support Notch-dependent regulation of the Hey gene, their role is redundant, as a Hey loss-of-function lethal phenotype is observed only after deletion of all their Su(H) binding motifs by CRISPR/Cas9.

Validated Negative Regions (VNRs) in the VISTA Database might be Truncated Forms of Bona Fide Enhancers

The VISTA enhancer database is a valuable resource for evaluating predicted enhancers in humans and mice. In addition to thousands of validated positive regions (VPRs) in the human and mouse genomes, the database also contains similar numbers of validated negative regions (VNRs). It is previously shown that the VPRs are on average half as long as predicted overlapping enhancers that are highly conserved and hypothesize that the VPRs may be truncated forms of long bona fide enhancers. Here, it is shown that like the VPRs, the VNRs also are under strong evolutionary constraints and overlap predicted enhancers in the genomes. The VNRs are also on average half as long as predicted overlapping enhancers that are highly conserved. Moreover, the VNRs and the VPRs display similar cell/tissue-specific modification patterns of key epigenetic marks of active enhancers. Furthermore, the VNRs and the VPRs show similar impact score spectra of in silico mutagenesis. These highly similar properties between the VPRs and the VNRs suggest that like the VPRs, the VNRs may also be truncated forms of long bona fide enhancers.

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.

Cell-type-directed design of synthetic enhancers

Transcriptional enhancers act as docking stations for combinations of transcription factors and thereby regulate spatiotemporal activation of their target genes1. It has been a long-standing goal in the field to decode the regulatory logic of an enhancer and to understand the details of how spatiotemporal gene expression is encoded in an enhancer sequence. Here we show that deep learning models2-6, can be used to efficiently design synthetic, cell-type-specific enhancers, starting from random sequences, and that this optimization process allows detailed tracing of enhancer features at single-nucleotide resolution. We evaluate the function of fully synthetic enhancers to specifically target Kenyon cells or glial cells in the fruit fly brain using transgenic animals. We further exploit enhancer design to create 'dual-code' enhancers that target two cell types and minimal enhancers smaller than 50 base pairs that are fully functional. By examining the state space searches towards local optima, we characterize enhancer codes through the strength, combination and arrangement of transcription factor activator and transcription factor repressor motifs. Finally, we apply the same strategies to successfully design human enhancers, which adhere to enhancer rules similar to those of Drosophila enhancers. Enhancer design guided by deep learning leads to better understanding of how enhancers work and shows that their code can be exploited to manipulate cell states.

Shaping gene expression and its evolution by chromatin architecture and enhancer activity

Transcriptional regulation plays a pivotal role in orchestrating the intricate genetic programs governing embryonic development. The expression of developmental genes relies on the combined activity of several cis-regulatory elements (CREs), such as enhancers and silencers, which can be located at long linear distances from the genes that they regulate and that interact with them through establishment of chromatin loops. Mutations affecting their activity or interaction with their target genes can lead to developmental disorders and are thought to have importantly contributed to the evolution of the animal body plan. The income of next-generation-sequencing approaches has allowed identifying over a million of sequences with putative regulatory potential in the human genome. Characterizing their function and establishing gene-CREs maps is essential to decode the logic governing developmental gene expression and is one of the major challenges of the post-genomic era. Chromatin 3D organization plays an essential role in determining how CREs specifically contact their target genes while avoiding deleterious off-target interactions. Our understanding of these aspects has greatly advanced with the income of chromatin conformation capture techniques and fluorescence microscopy approaches to visualize the organization of DNA elements in the nucleus. Here we will summarize relevant aspects of how the interplay between CRE activity and chromatin 3D organization regulates developmental gene expression and how it relates to pathological conditions and the evolution of animal body plan.

The Density of Regulatory Information Is a Major Determinant of Evolutionary Constraint on Noncoding DNA in Drosophila

Evolutionary analyses have estimated that ∼60% of nucleotides in intergenic regions of the Drosophila melanogaster genome are functionally relevant, suggesting that regulatory information may be encoded more densely in intergenic regions than has been revealed by most functional dissections of regulatory DNA. Here, we approached this issue through a functional dissection of the regulatory region of the gene shavenbaby (svb). Most of the ∼90 kb of this large regulatory region is highly conserved in the genus Drosophila, though characterized enhancers occupy a small fraction of this region. By analyzing the regulation of svb in different contexts of Drosophila development, we found that the regulatory information that drives svb expression in the abdominal pupal epidermis is organized in a different way than the elements that drive svb expression in the embryonic epidermis. While in the embryonic epidermis svb is activated by compact enhancers separated by large inactive DNA regions, svb expression in the pupal epidermis is driven by regulatory information distributed over broader regions of svb cis-regulatory DNA. In the same vein, we observed that other developmental genes also display a dense distribution of putative regulatory elements in their regulatory regions. Furthermore, we found that a large percentage of conserved noncoding DNA of the Drosophila genome is contained within regions of open chromatin. These results suggest that part of the evolutionary constraint on noncoding DNA of Drosophila is explained by the density of regulatory information, which may be greater than previously appreciated.

Genetics and the Evolution of Prezygotic Isolation

The significance of prezygotic isolation for speciation has been recognized at least since the Modern Synthesis. However, fundamental questions remain. For example, how are genetic associations between traits that contribute to prezygotic isolation maintained? What is the source of genetic variation underlying the evolution of these traits? And how do prezygotic barriers affect patterns of gene flow? We address these questions by reviewing genetic features shared across plants and animals that influence prezygotic isolation. Emerging technologies increasingly enable the identification and functional characterization of the genes involved, allowing us to test established theoretical expectations. Embedding these genes in their developmental context will allow further predictions about what constrains the evolution of prezygotic isolation. Ongoing improvements in statistical and computational tools will reveal how pre- and postzygotic isolation may differ in how they influence gene flow across the genome. Finally, we highlight opportunities for progress by combining theory with appropriate data.

Shifting landscapes: the role of 3D genomic organizations in gene regulatory strategies

3D genome folding enables the physical storage of chromosomes into the compact volume of a cell's nucleus, allows for the accurate segregation of chromatin to daughter cells, and has been shown to be tightly coupled to the way in which genetic information is converted into transcriptional programs [1-3]. Importantly, this link between chromatin architecture and gene regulation is a selectable feature in which modifications to chromatin organization accompany, or perhaps even drive the establishment of new regulatory strategies with enduring impacts on animal body plan complexity. Here, we discuss the nature of different 3D genome folding systems found across the tree of life, with particular emphasis on metazoans, and the relative influence of these systems on gene regulation. We suggest how the properties of these folding systems have influenced regulatory strategies employed by different lineages and may have catalyzed the partitioning and specialization of genetic programs that enabled multicellularity and organ-grade body plan complexity.

Genome-Wide Association Mapping for Yield and Yield-Related Traits in Rice (Oryza Sativa L.) Using SNPs Markers

Rice (Oryza sativa L.) is a staple food for more than 50% of the world's population. Rice cultivar improvement is critical in order to feed the world's growing population. Improving yield is one of the main aims of rice breeders. However, yield is a complex quantitative trait controlled by many genes. The presence of genetic diversity is the key factor to improve the yield hence, the presence of diversity in any germplasm is important for yield improvement. In the current study, the rice germplasm was collected from Pakistan and the United States of America and a panel of 100 diverse genotypes was utilized to identify important yield and yield-related traits. For this, a genome-wide association study (GWAS) was performed to identify the genetic loci related to yield. The GWAS on the diverse germplasm will lead to the identification of new genes which can be utilized in the breeding program for improvement of yield. For this reason, firstly, the germplasm was phenotypically evaluated in two growing seasons for yield and yield-related traits. The analysis of variance results showed significant differences among traits which showed the presence of diversity in the current germplasm. Secondly, the germplasm was also genotypically evaluated using 10K SNP. Genetic structure analysis showed the presence of four groups which showed that enough genetic diversity was present in the rice germplasm to be used for association mapping analysis. The results of GWAS identified 201 significant marker trait associations (MTAs. 16 MTAs were identified for plant height, 49 for days to flowering, three for days to maturity, four for tillers per plant, four for panicle length, eight for grains per panicle, 20 unfilled grains per panicle, 81 for seed setting %, four for thousand-grain weight, five for yield per plot and seven for yield per hectare. Apart from this, some pleiotropic loci were also identified. The results showed that panicle length (PL) and thousand-grain weight (TGW) were controlled by a pleiotropic locus OsGRb23906 on chromosome 1 at 10,116,371 cM. The loci OsGRb25803 and OsGRb15974 on chromosomes 4 and 8 at the position of 14,321,111 cM and 6,205,816 cM respectively, showed pleiotropic effects for seed setting % (SS) and unfilled grain per panicle (UG/P). A locus OsGRb09180 on chromosome 4 at 19,850,601 cM was significantly linked with SS and yield/ha. Furthermore, gene annotation was performed, and results indicated that the 190 candidate genes or QTLs that closely linked with studied traits. These candidate genes and novel significant markers could be useful in marker-assisted gene selection and QTL pyramiding to improve rice yield and the selection of potential parents, recombinants and MTAs which could be used in rice breeding programs to develop high-yielding rice varieties for sustainable food security.

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.

Enhancer architecture and chromatin accessibility constrain phenotypic space during Drosophila development

Developmental enhancers bind transcription factors and dictate patterns of gene expression during development. Their molecular evolution can underlie phenotypical evolution, but the contributions of the evolutionary pathways involved remain little understood. Here, using mutation libraries in Drosophila melanogaster embryos, we observed that most point mutations in developmental enhancers led to changes in gene expression levels but rarely resulted in novel expression outside of the native pattern. In contrast, random sequences, often acting as developmental enhancers, drove expression across a range of cell types; random sequences including motifs for transcription factors with pioneer activity acted as enhancers even more frequently. Our findings suggest that the phenotypic landscapes of developmental enhancers are constrained by enhancer architecture and chromatin accessibility. We propose that the evolution of existing enhancers is limited in its capacity to generate novel phenotypes, whereas the activity of de novo elements is a primary source of phenotypic novelty.

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.

Evolutionary Changes in the Chromatin Landscape Contribute to Reorganization of a Developmental Gene Network During Rapid Life History Evolution in Sea Urchins

Chromatin configuration is highly dynamic during embryonic development in animals, exerting an important point of control in transcriptional regulation. Yet there exists remarkably little information about the role of evolutionary changes in chromatin configuration to the evolution of gene expression and organismal traits. Genome-wide assays of chromatin configuration, coupled with whole-genome alignments, can help address this gap in knowledge in several ways. In this study we present a comparative analysis of regulatory element sequences and accessibility throughout embryogenesis in three sea urchin species with divergent life histories: a lecithotroph Heliocidaris erythrogramma, a closely related planktotroph H. tuberculata, and a distantly related planktotroph Lytechinus variegatus. We identified distinct epigenetic and mutational signatures of evolutionary modifications to the function of putative cis-regulatory elements in H. erythrogramma that have accumulated nonuniformly throughout the genome, suggesting selection, rather than drift, underlies many modifications associated with the derived life history. Specifically, regulatory elements composing the sea urchin developmental gene regulatory network are enriched for signatures of positive selection and accessibility changes which may function to alter binding affinity and access of developmental transcription factors to these sites. Furthermore, regulatory element changes often correlate with divergent expression patterns of genes involved in cell type specification, morphogenesis, and development of other derived traits, suggesting these evolutionary modifications have been consequential for phenotypic evolution in H. erythrogramma. Collectively, our results demonstrate that selective pressures imposed by changes in developmental life history rapidly reshape the cis-regulatory landscape of core developmental genes to generate novel traits and embryonic programs.

Temporal analysis of enhancers during mouse cerebellar development reveals dynamic and novel regulatory functions

We have identified active enhancers in the mouse cerebellum at embryonic and postnatal stages which provides a view of novel enhancers active during cerebellar development. The majority of cerebellar enhancers have dynamic activity between embryonic and postnatal development. Cerebellar enhancers were enriched for neural transcription factor binding sites with temporally specific expression. Putative gene targets displayed spatially restricted expression patterns, indicating cell-type specific expression regulation. Functional analysis of target genes indicated that enhancers regulate processes spanning several developmental epochs such as specification, differentiation and maturation. We use these analyses to discover one novel regulator and one novel marker of cerebellar development: Bhlhe22 and Pax3, respectively. We identified an enrichment of de novo mutations and variants associated with autism spectrum disorder in cerebellar enhancers. Furthermore, by comparing our data with relevant brain development ENCODE histone profiles and cerebellar single-cell datasets we have been able to generalize and expand on the presented analyses, respectively. We have made the results of our analyses available online in the Developing Mouse Cerebellum Enhancer Atlas (https://goldowitzlab.shinyapps.io/developing_mouse_cerebellum_enhancer_atlas/), where our dataset can be efficiently queried, curated and exported by the scientific community to facilitate future research efforts. Our study provides a valuable resource for studying the dynamics of gene expression regulation by enhancers in the developing cerebellum and delivers a rich dataset of novel gene-enhancer associations providing a basis for future in-depth studies in the cerebellum.

Pleiotropic Enhancers are Ubiquitous Regulatory Elements in the Human Genome

Enhancers are regulatory elements of genomes that determine spatio-temporal patterns of gene expression. The human genome contains a vast number of enhancers, which largely outnumber protein-coding genes. Historically, enhancers have been regarded as highly tissue-specific. However, recent evidence has demonstrated that many enhancers are pleiotropic, with activity in multiple developmental contexts. Yet, the extent and impact of pleiotropy remain largely unexplored. In this study we analyzed active enhancers across human organs based on the analysis of both eRNA transcription (FANTOM5 consortium data sets) and chromatin architecture (ENCODE consortium data sets). We show that pleiotropic enhancers are pervasive in the human genome and that most enhancers active in a particular organ are also active in other organs. In addition, our analysis suggests that the proportion of context-specific enhancers of a given organ is explained, at least in part, by the proportion of context-specific genes in that same organ. The notion that such a high proportion of human enhancers can be pleiotropic suggests that small regions of regulatory DNA contain abundant regulatory information and that these regions evolve under important evolutionary constraints.

The role of the epidermis enhancer element in positive and negative transcriptional regulation of ebony in Drosophila melanogaster

The spatiotemporal regulation of gene expression is essential to ensure robust phenotypic outcomes. Pigmentation patterns in Drosophila are determined by pigments biosynthesized in the developing epidermis and the cis-regulatory elements of the genes involved in this process are well-characterized. Here, we report that the known primary epidermal enhancer is dispensable for the transcriptional activation of ebony (involved in light-colored pigment synthesis) in the developing epidermis of Drosophila melanogaster. The evidence was obtained by introducing an approximately 1 kbp deletion at the primary epidermal enhancer by genome editing. The effect of the primary epidermal enhancer deletion on pigmentation and on the endogenous expression pattern of a mCherry-fused ebony allele was examined in the abdomen. The expression levels of the mCherry-fused ebony in the primary epidermal enhancer-deleted strains were slightly higher than that of the control strain, indicating that the sequences outside the primary epidermal enhancer have an ability to drive an expression of this gene in the epidermis. Interestingly, the primary epidermal enhancer deletion resulted in a derepression of this gene in the dorsal midline of the abdominal tergites, where dark pigmentation is present in the wild-type individuals. This indicated that the primary epidermal enhancer fragment contains a silencer. Furthermore, the endogenous expression pattern of ebony in the 2 additional strains with partially deleted primary epidermal enhancer revealed that the silencer resides within a 351-bp fragment in the 5' portion of the primary epidermal enhancer. These results demonstrated that deletion assays combined with reporter assays are highly effective in detecting the presence of positively and negatively regulating sequences within and outside the focal cis-regulatory elements.

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.

The Non-continuum Nature of Eukaryotic Transcriptional Regulation

Eukaryotic transcription factors are versatile mediators of specificity in gene regulation. This versatility is achieved through mutual specification by context-specific DNA binding on the one hand, and identity-specific protein-protein partnerships on the other. This interactivity, known as combinatorial control, enables a repertoire of complex transcriptional outputs that are qualitatively disjoint, or non-continuum, with respect to binding affinity. This feature contrasts starkly with prokaryotic gene regulators, whose activities in general vary quantitatively in step with binding affinity. Biophysical studies on prokaryotic model systems and more recent investigations on transcription factors highlight an important role for folded state dynamics and molecular hydration in protein/DNA recognition. Analysis of molecular models of combinatorial control and recent literature in low-affinity gene regulation suggest that transcription factors harbor unique conformational dynamics that are inaccessible or unused by prokaryotic DNA-binding proteins. Thus, understanding the intrinsic dynamics involved in DNA binding and co-regulator recruitment appears to be a key to understanding how transcription factors mediate non-continuum outcomes in eukaryotic gene expression, and how such capability might have evolved from ancient, structurally conserved counterparts.

Delilah, prospero, and D-Pax2 constitute a gene regulatory network essential for the development of functional proprioceptors

Coordinated animal locomotion depends on the development of functional proprioceptors. While early cell-fate determination processes are well characterized, little is known about the terminal differentiation of cells within the proprioceptive lineage and the genetic networks that control them. In this work we describe a gene regulatory network consisting of three transcription factors–Prospero (Pros), D-Pax2, and Delilah (Dei)–that dictates two alternative differentiation programs within the proprioceptive lineage in Drosophila. We show that D-Pax2 and Pros control the differentiation of cap versus scolopale cells in the chordotonal organ lineage by, respectively, activating and repressing the transcription of dei. Normally, D-Pax2 activates the expression of dei in the cap cell but is unable to do so in the scolopale cell where Pros is co-expressed. We further show that D-Pax2 and Pros exert their effects on dei transcription via a 262 bp chordotonal-specific enhancer in which two D-Pax2- and three Pros-binding sites were identified experimentally. When this enhancer was removed from the fly genome, the cap- and ligament-specific expression of dei was lost, resulting in loss of chordotonal organ functionality and defective larval locomotion. Thus, coordinated larval locomotion depends on the activity of a dei enhancer that integrates both activating and repressive inputs for the generation of a functional proprioceptive organ.

Hearing loss genes reveal patterns of adaptive evolution at the coding and non-coding levels in mammals

Background

Mammals possess unique hearing capacities that differ significantly from those of the rest of the amniotes. In order to gain insights into the evolution of the mammalian inner ear, we aim to identify the set of genetic changes and the evolutionary forces that underlie this process. We hypothesize that genes that impair hearing when mutated in humans or in mice (hearing loss (HL) genes) must play important roles in the development and physiology of the inner ear and may have been targets of selective forces across the evolution of mammals. Additionally, we investigated if these HL genes underwent a human-specific evolutionary process that could underlie the evolution of phenotypic traits that characterize human hearing.

Results

We compiled a dataset of HL genes including non-syndromic deafness genes identified by genetic screenings in humans and mice. We found that many genes including those required for the normal function of the inner ear such as LOXHD1, TMC1, OTOF, CDH23, and PCDH15 show strong signatures of positive selection. We also found numerous noncoding accelerated regions in HL genes, and among them, we identified active transcriptional enhancers through functional enhancer assays in transgenic zebrafish.

Conclusions

Our results indicate that the key inner ear genes and regulatory regions underwent adaptive evolution in the basal branch of mammals and along the human-specific branch, suggesting that they could have played an important role in the functional remodeling of the cochlea. Altogether, our data suggest that morphological and functional evolution could be attained through molecular changes affecting both coding and noncoding regulatory regions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01170-6.

Pri smORF Peptides Are Wide Mediators of Ecdysone Signaling, Contributing to Shape Spatiotemporal Responses

There is growing evidence that peptides encoded by small open-reading frames (sORF or smORF) can fulfill various cellular functions and define a novel class regulatory molecules. To which extend transcripts encoding only smORF peptides compare with canonical protein-coding genes, yet remain poorly understood. In particular, little is known on whether and how smORF-encoding RNAs might need tightly regulated expression within a given tissue, at a given time during development. We addressed these questions through the analysis of Drosophila polished rice (pri, a.k.a. tarsal less or mille pattes), which encodes four smORF peptides (11-32 amino acids in length) required at several stages of development. Previous work has shown that the expression of pri during epidermal development is regulated in the response to ecdysone, the major steroid hormone in insects. Here, we show that pri transcription is strongly upregulated by ecdysone across a large panel of cell types, suggesting that pri is a core component of ecdysone response. Although pri is produced as an intron-less short transcript (1.5 kb), genetic assays reveal that the developmental functions of pri require an unexpectedly large array of enhancers (spanning over 50 kb), driving a variety of spatiotemporal patterns of pri expression across developing tissues. Furthermore, we found that separate pri enhancers are directly activated by the ecdysone nuclear receptor (EcR) and display distinct regulatory modes between developmental tissues and/or stages. Alike major developmental genes, the expression of pri in a given tissue often involves several enhancers driving apparently redundant (or shadow) expression, while individual pri enhancers can harbor pleiotropic functions across tissues. Taken together, these data reveal the broad role of Pri smORF peptides in ecdysone signaling and show that the cis-regulatory architecture of the pri gene contributes to shape distinct spatial and temporal patterns of ecdysone response throughout development.

Favorable pleiotropic loci for fiber yield and quality in upland cotton (Gossypium hirsutum)

Upland cotton (Gossypium hirsutum L.) is an important economic crop for renewable textile fibers. However, the simultaneous improvement of yield and fiber quality in cotton is difficult as the linkage drag. Compared with breaking the linkage drag, identification of the favorable pleiotropic loci on the genome level by genome-wide association study (GWAS) provides a new way to improve the yield and fiber quality simultaneously. In our study restriction-site-associated DNA sequencing (RAD-seq) was used to genotype 316 cotton accessions. Eight major traits in three categories including yield, fiber quality and maturation were investigated in nine environments (3 sites × 3 years). 231 SNPs associated with these eight traits (− log₁₀(P) > 5.27) were identified, located in 27 genomic regions respectively by linkage disequilibrium analysis. Further analysis showed that four genomic regions (the region 1, 6, 8 and 23) held favorable pleiotropic loci and 6 candidate genes were identified. Through genotyping, 14 elite accessions carrying the favorable loci on four pleiotropic regions were identified. These favorable pleiotropic loci and elite genotypes identified in this study will be utilized to improve the yield and fiber quality simultaneously in future cotton breeding.

Enhancer pleiotropy, gene expression, and the architecture of human enhancer-gene interactions

Enhancers are often studied as noncoding regulatory elements that modulate the precise spatiotemporal expression of genes in a highly tissue-specific manner. This paradigm has been challenged by recent evidence of individual enhancers acting in multiple tissues or developmental contexts. However, the frequency of these enhancers with high degrees of “pleiotropy” out of all putative enhancers is not well understood. Consequently, it is unclear how the variation of enhancer pleiotropy corresponds to the variation in expression breadth of target genes. Here, we use multi-tissue chromatin maps from diverse human tissues to investigate the enhancer–gene interaction architecture while accounting for 1) the distribution of enhancer pleiotropy, 2) the variations of regulatory links from enhancers to target genes, and 3) the expression breadth of target genes. We show that most enhancers are tissue-specific and that highly pleiotropy enhancers account for <1% of all putative regulatory sequences in the human genome. Notably, several genomic features are indicative of increasing enhancer pleiotropy, including longer sequence length, greater number of links to genes, increasing abundance and diversity of encoded transcription factor motifs, and stronger evolutionary conservation. Intriguingly, the number of enhancers per gene remains remarkably consistent for all genes (∼14). However, enhancer pleiotropy does not directly translate to the expression breadth of target genes. We further present a series of Gaussian Mixture Models to represent this organization architecture. Consequently, we demonstrate that a modest trend of more pleiotropic enhancers targeting more broadly expressed genes can generate the observed diversity of expression breadths in the human genome.

Steroid-dependent switch of OvoL/Shavenbaby controls self-renewal versus differentiation of intestinal stem cells

Adult stem cells must continuously fine-tune their behavior to regenerate damaged organs and avoid tumors. While several signaling pathways are well known to regulate somatic stem cells, the underlying mechanisms remain largely unexplored. Here, we demonstrate a cell-intrinsic role for the OvoL family transcription factor, Shavenbaby (Svb), in balancing self-renewal and differentiation of Drosophila intestinal stem cells. We find that svb is a downstream target of Wnt and EGFR pathways, mediating their activity for stem cell survival and proliferation. This requires post-translational processing of Svb into a transcriptional activator, whose upregulation induces tumor-like stem cell hyperproliferation. In contrast, the unprocessed form of Svb acts as a repressor that imposes differentiation into enterocytes, and suppresses tumors induced by altered signaling. We show that the switch between Svb repressor and activator is triggered in response to systemic steroid hormone, which is produced by ovaries. Therefore, the Svb axis allows intrinsic integration of local signaling cues and inter-organ communication to adjust stem cell proliferation versus differentiation, suggesting a broad role of OvoL/Svb in adult and cancer stem cells.

A complex gene regulatory architecture underlies the development and evolution of cuticle morphology in Drosophila

The cuticle of insects is decorated with non-sensory hairs called trichomes. A few Drosophila species independently lost most of the dorso-lateral trichomes on first instar larvae. Genetic experiments revealed that this naked cuticle phenotype was caused by the evolution of enhancer function at the ovo/shavenbaby (ovo/svb) locus. Here we explore how this discovery catalyzed major new insights into morphological evolution in different developmental contexts, enhancer pleiotropy in gene regulation and the functionality and evolution of the Svb gene regulatory network (GRN). Taken together this highlights the importance of understanding the architecture and evolution of gene regulatory networks in detail and the great potential for further study of the Svb GRN.

Cis-regulatory variation in the shavenbaby gene underlies intraspecific phenotypic variation, mirroring interspecific divergence in the same trait

Despite considerable progress in recent decades in dissecting the genetic causes of natural morphological variation, there is limited understanding of how variation within species ultimately contributes to species differences. We have studied patterning of the non-sensory hairs, commonly known as "trichomes," on the dorsal cuticle of first-instar larvae of Drosophila. Most Drosophila species produce a dense lawn of dorsal trichomes, but a subset of these trichomes were lost in D. sechellia and D. ezoana due entirely to regulatory evolution of the shavenbaby (svb) gene. Here, we describe intraspecific variation in dorsal trichome patterns of first-instar larvae of D. virilis that is similar to the trichome pattern variation identified previously between species. We found that a single large effect QTL, which includes svb, explains most of the trichome number difference between two D. virilis strains and that svb expression correlates with the trichome difference between strains. This QTL does not explain the entire difference between strains, implying that additional loci contribute to variation in trichome numbers. Thus, the genetic architecture of intraspecific variation exhibits similarities and differences with interspecific variation that may reflect differences in long-term and short-term evolutionary processes.

Three distinct mechanisms, Notch instructive, permissive, and independent, regulate the expression of two different pericardial genes to specify cardiac cell subtypes

The development of a complex organ involves the specification and differentiation of diverse cell types constituting that organ. Two major cell subtypes, contractile cardial cells (CCs) and nephrocytic pericardial cells (PCs), comprise the Drosophila heart. Binding sites for Suppressor of Hairless [Su(H)], an integral transcription factor in the Notch signaling pathway, are enriched in the enhancers of PC-specific genes. Here we show three distinct mechanisms regulating the expression of two different PC-specific genes, Holes in muscle (Him), and Zn finger homeodomain 1 (zfh1). Him transcription is activated in PCs in a permissive manner by Notch signaling: in the absence of Notch signaling, Su(H) forms a repressor complex with co-repressors and binds to the Him enhancer, repressing its transcription; upon alleviation of this repression by Notch signaling, Him transcription is activated. In contrast, zfh1 is transcribed by a Notch-instructive mechanism in most PCs, where mere alleviation of repression by preventing the binding of Su(H)-co-repressor complex is not sufficient to activate transcription. Our results suggest that upon activation of Notch signaling, the Notch intracellular domain associates with Su(H) to form an activator complex that binds to the zfh1 enhancer, and that this activator complex is necessary for bringing about zfh1 transcription in these PCs. Finally, a third, Notch-independent mechanism activates zfh1 transcription in the remaining, even skipped-expressing, PCs. Collectively, our data show how the same feature, enrichment of Su(H) binding sites in PC-specific gene enhancers, is utilized by two very distinct mechanisms, one permissive, the other instructive, to contribute to the same overall goal: the specification and differentiation of a cardiac cell subtype by activation of the pericardial gene program. Furthermore, our results demonstrate that the zfh1 enhancer drives expression in two different domains using distinct Notch-instructive and Notch-independent mechanisms.

Dense and pleiotropic regulatory information in a developmental enhancer

Changes in gene regulation underlie much of phenotypic evolution¹. However, our understanding of the potential for regulatory evolution is biased, because most evidence comes from either natural variation or limited experimental perturbations². Using an automated robotics pipeline, we surveyed an unbiased mutation library for a developmental enhancer in Drosophila melanogaster. We found that almost all mutations altered gene expression and that parameters of gene expression-levels, location, and state-were convolved. The widespread pleiotropic effects of most mutations may constrain the evolvability of developmental enhancers. Consistent with these observations, comparisons of diverse Drosophila larvae revealed apparent biases in the phenotypes influenced by the enhancer. Developmental enhancers may encode a higher density of regulatory information than has been appreciated previously, imposing constraints on regulatory evolution.

Systemic evolutionary changes in mammalian gene expression

Changes in gene expression play an important role in evolution and can be relevant to evolutionary medicine. In this work, a strong relationship was found between the statistical significance of evolutionary changes in the expression of orthologous genes in the five or six homologous mammalian tissues and the across-tissues unidirectionality of changes (i.e., they occur in the same direction in different tissues -- all upward or all downward). In the area of highly significant changes, the fraction of unidirectionally changed genes (UCG) was above 0.9 (random expectation is 0.03). This observation indicates that the most pronounced evolutionary changes in mammalian gene expression are systemic (i.e., they operate at the whole-organism level). The UCG are strongly enriched in the housekeeping genes. More specifically, in the human-chimpanzee comparison, the UCG are enriched in the pathways belonging to gene expression (translation is prominent), cell cycle control, ubiquitin-dependent protein degradation (mostly related to cell cycle control), apoptosis, and Parkinson's disease. In the human-macaque comparison, the two other neurodegenerative diseases (Alzheimer's and Huntington's) are added to the enriched pathways. The consolidation of gene expression changes at the level of pathways indicates that they are not neutral but functional. The systemic expression changes probably maintain the across-tissues balance of basic physiological processes in the course of evolution (e.g., during the movement along the fast-slow life axis). These results can be useful for understanding the variation in longevity and susceptibility to cancer and widespread neurodegenerative diseases. This approach can also guide the choice of prospective genes for studies aiming to decipher cis-regulatory code (the gene list is provided).

Visual mate preference evolution during butterfly speciation is linked to neural processing genes

Many animal species remain separate not because their individuals fail to produce viable hybrids but because they "choose" not to mate. However, we still know very little of the genetic mechanisms underlying changes in these mate preference behaviours. Heliconius butterflies display bright warning patterns, which they also use to recognize conspecifics. Here, we couple QTL for divergence in visual preference behaviours with population genomic and gene expression analyses of neural tissue (central brain, optic lobes and ommatidia) across development in two sympatric Heliconius species. Within a region containing 200 genes, we identify five genes that are strongly associated with divergent visual preferences. Three of these have previously been implicated in key components of neural signalling (specifically an ionotropic glutamate receptor and two regucalcins), and overall our candidates suggest shifts in behaviour involve changes in visual integration or processing. This would allow preference evolution without altering perception of the wider environment.

Transcriptional Enhancers in Drosophila

Key discoveries in Drosophila have shaped our understanding of cellular "enhancers." With a special focus on the fly, this chapter surveys properties of these adaptable cis-regulatory elements, whose actions are critical for the complex spatial/temporal transcriptional regulation of gene expression in metazoa. The powerful combination of genetics, molecular biology, and genomics available in Drosophila has provided an arena in which the developmental role of enhancers can be explored. Enhancers are characterized by diverse low- or high-throughput assays, which are challenging to interpret, as not all of these methods of identifying enhancers produce concordant results. As a model metazoan, the fly offers important advantages to comprehensive analysis of the central functions that enhancers play in gene expression, and their critical role in mediating the production of phenotypes from genotype and environmental inputs. A major challenge moving forward will be obtaining a quantitative understanding of how these cis-regulatory elements operate in development and disease.

Shadow enhancers can suppress input transcription factor noise through distinct regulatory logic

Shadow enhancers, groups of seemingly redundant enhancers, are found in a wide range of organisms and are critical for robust developmental patterning. However, their mechanism of action is unknown. We hypothesized that shadow enhancers drive consistent expression levels by buffering upstream noise through a separation of transcription factor (TF) inputs at the individual enhancers. By measuring the transcriptional dynamics of several Kruppel shadow enhancer configurations in live Drosophila embryos, we showed that individual member enhancers act largely independently. We found that TF fluctuations are an appreciable source of noise that the shadow enhancer pair can better buffer than duplicated enhancers. The shadow enhancer pair is also uniquely able to maintain low levels of expression noise across a wide range of temperatures. A stochastic model demonstrated the separation of TF inputs is sufficient to explain these findings. Our results suggest the widespread use of shadow enhancers is partially due to their noise suppressing ability.

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.

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).

Developmental and Evolutionary Allometry of the Mammalian Limb Skeleton

The variety of limb skeletal proportions enables a remarkable diversity of behaviors that include powered flight in bats and flipper-propelled swimming in whales using extremes of a range of homologous limb architectures. Even within human limbs, bone lengths span more than an order of magnitude from the short finger and toe bones to the long arm and leg bones. Yet all of this diversity arises from embryonic skeletal elements that are each a very similar size at formation. In this review article, I survey what is and is not yet known of the development and evolution of skeletal proportion at multiple hierarchical levels of biological organization. These include the cellular parameters of skeletal elongation in the cartilage growth plate, genes associated with differential growth, and putative gene regulatory mechanisms that would allow both covariant and independent evolution of the forelimbs and hindlimbs and of individual limb segments. Although the genetic mechanisms that shape skeletal proportion are still largely unknown, and most of what is known is limited to mammals, it is becoming increasingly apparent that the diversity of bone lengths is an emergent property of a complex system that controls elongation of individual skeletal elements using a genetic toolkit shared by all.

Trajectories of Serum Sodium on In-Hospital and 1-Year Survival among Hospitalized Patients

BACKGROUND AND OBJECTIVES

This study aimed to investigate the association between in-hospital trajectories of serum sodium and risk of in-hospital and 1-year mortality in patients in hospital.

DESIGN, SETTING, PARTICIPANTS, & MEASUREMENTS

This is a single-center cohort study. All adult patients who were hospitalized from years 2011 through 2013 who had available admission serum sodium and at least three serum sodium measurements during hospitalization were included. The trend of serum sodium during hospitalization was analyzed using group-based trajectory modeling; the five main trajectories were grouped as follows: (1) stable normonatremia, (2) uncorrected hyponatremia, (3) borderline high serum sodium, (4) corrected hyponatremia, and (5) fluctuating serum sodium. The outcome of interest was in-hospital mortality and 1-year mortality. Stable normonatremia was used as the reference group for outcome comparison.

RESULTS

A total of 43,539 patients were analyzed. Of these, 47% had stable normonatremia, 15% had uncorrected hyponatremia, 31% had borderline high serum sodium, 3% had corrected hyponatremia, and 5% had fluctuating serum sodium trajectory. In adjusted analysis, there was a higher in-hospital mortality among those with uncorrected hyponatremia (odds ratio [OR], 1.33; 95% CI, 1.06 to 1.67), borderline high serum sodium (OR, 1.66; 95% CI, 1.38 to 2.00), corrected hyponatremia (OR, 1.50; 95% CI, 1.02 to 2.20), and fluctuating serum sodium (OR, 4.61; 95% CI, 3.61 to 5.88), compared with those with the normonatremia trajectory. One-year mortality was higher among those with uncorrected hyponatremia (hazard ratio [HR], 1.28; 95% CI, 1.19 to 1.38), borderline high serum sodium (HR, 1.18; 95% CI, 1.11 to 1.26), corrected hyponatremia (HR, 1.24; 95% CI, 1.08 to 1.42), and fluctuating serum sodium (HR, 2.10; 95% CI, 1.89 to 2.33) compared with those with the normonatremia trajectory.

CONCLUSIONS

More than half of patients who had been hospitalized had an abnormal serum sodium trajectory during hospitalization. This study demonstrated that not only the absolute serum sodium levels but also their in-hospital trajectories were significantly associated with in-hospital and 1-year mortality. The highest in-hospital and 1-year mortality risk was associated with the fluctuating serum sodium trajectory.

PODCAST

This article contains a podcast at https://www.asn-online.org/media/podcast/CJASN/2020_03_25_CJN.12281019.mp3.

An Atlas of Transcription Factors Expressed in Male Pupal Terminalia of Drosophila melanogaster

During development, transcription factors and signaling molecules govern gene regulatory networks to direct the formation of unique morphologies. As changes in gene regulatory networks are often implicated in morphological evolution, mapping transcription factor landscapes is important, especially in tissues that undergo rapid evolutionary change. The terminalia (genital and anal structures) of Drosophila melanogaster and its close relatives exhibit dramatic changes in morphology between species. While previous studies have identified network components important for patterning the larval genital disc, the networks governing adult structures during pupal development have remained uncharted. Here, we performed RNA-seq in whole Drosophila melanogaster male terminalia followed by in situ hybridization for 100 highly expressed transcription factors during pupal development. We find that the male terminalia are highly patterned during pupal stages and that specific transcription factors mark separate structures and substructures. Our results are housed online in a searchable database (https://flyterminalia.pitt.edu/) as a resource for the community. This work lays a foundation for future investigations into the gene regulatory networks governing the development and evolution of Drosophila terminalia.

Neurobiological functions of transcriptional enhancers

Transcriptional enhancers are regulatory DNA elements that underlie the specificity and dynamic patterns of gene expression. Over the past decade, large-scale functional genomics projects have driven transformative progress in our understanding of enhancers. These data have relevance for identifying mechanisms of gene regulation in the CNS, elucidating the function of non-coding regulatory sequences in neurobiology and linking sequence variation within enhancers to genetic risk for neurological and psychiatric disorders. However, the sheer volume and complexity of genomic data presents a challenge to interpreting enhancer function in normal and pathogenic neurobiological processes. Here, to advance the application of genome-scale enhancer data, we offer a primer on current models of enhancer function in the CNS, we review how enhancers regulate gene expression across the neuronal lifespan, and we suggest how emerging findings regarding the role of non-coding sequence variation offer opportunities for understanding brain disorders and developing new technologies for neuroscience.

Parallel evolution of ancient, pleiotropic enhancers underlies butterfly wing pattern mimicry

Color pattern mimicry in Heliconius butterflies is a classic case study of complex trait adaptation via selection on a few large effect genes. Association studies have linked color pattern variation to a handful of noncoding regions, yet the presumptive cis-regulatory elements (CREs) that control color patterning remain unknown. Here we combine chromatin assays, DNA sequence associations, and genome editing to functionally characterize 5 cis-regulatory elements of the color pattern gene optix We were surprised to find that the cis-regulatory architecture of optix is characterized by pleiotropy and regulatory fragility, where deletion of individual cis-regulatory elements has broad effects on both color pattern and wing vein development. Remarkably, we found orthologous cis-regulatory elements associate with wing pattern convergence of distantly related comimics, suggesting that parallel coevolution of ancestral elements facilitated pattern mimicry. Our results support a model of color pattern evolution in Heliconius where changes to ancient, multifunctional cis-regulatory elements underlie adaptive radiation.

Influence of genetic polymorphism on transcriptional enhancer activity in the malaria vector Anopheles coluzzii

Enhancers are cis-regulatory elements that control most of the developmental and spatial gene expression in eukaryotes. Genetic variation of enhancer sequences is known to influence phenotypes, but the effect of enhancer variation upon enhancer functional activity and downstream phenotypes has barely been examined in any species. In the African malaria vector, Anopheles coluzzii, we identified candidate enhancers in the proximity of genes relevant for immunity, insecticide resistance, and development. The candidate enhancers were functionally validated using luciferase reporter assays, and their activity was found to be essentially independent of their physical orientation, a typical property of enhancers. All of the enhancers segregated genetically polymorphic alleles, which displayed significantly different levels of functional activity. Deletion mutagenesis and functional testing revealed a fine structure of positive and negative regulatory elements that modulate activity of the enhancer core. Enhancer polymorphisms segregate in wild A. coluzzii populations in West Africa. Thus, enhancer variants that modify target gene expression leading to likely phenotypic consequences are frequent in nature. These results demonstrate the existence of naturally polymorphic A. coluzzii enhancers, which may help explain important differences between individuals or populations for malaria transmission efficiency and vector adaptation to the environment.

Reconstruction of the Global Neural Crest Gene Regulatory Network In Vivo

Precise control of developmental processes is encoded in the genome in the form of gene regulatory networks (GRNs). Such multi-factorial systems are difficult to decode in vertebrates owing to their complex gene hierarchies and dynamic molecular interactions. Here we present a genome-wide in vivo reconstruction of the GRN underlying development of the multipotent neural crest (NC) embryonic cell population. By coupling NC-specific epigenomic and transcriptional profiling at population and single-cell levels with genome/epigenome engineering in vivo, we identify multiple regulatory layers governing NC ontogeny, including NC-specific enhancers and super-enhancers, novel trans-factors, and cis-signatures allowing reverse engineering of the NC-GRN at unprecedented resolution. Furthermore, identification and dissection of divergent upstream combinatorial regulatory codes has afforded new insights into opposing gene circuits that define canonical and neural NC fates early during NC ontogeny. Our integrated approach, allowing dissection of cell-type-specific regulatory circuits in vivo, has broad implications for GRN discovery and investigation.

Evolution: How Many Phenotypes Do Regulatory Mutations Affect?

Mutations in gene regulatory regions are thought to play an important role in the evolution of morphological structures. This is largely due to their minimal pleiotropic effects, limiting their impact to one particular body part. A recent study finds that one such regulatory mutation may affect two particular morphological structures.

Modular tissue-specific regulation of doublesex underpins sexually dimorphic development in Drosophila

The ability of a single genome to produce distinct and often dramatically different male and female forms is one of the wonders of animal development. In Drosophila melanogaster, most sexually dimorphic traits are controlled by sex-specific isoforms of the doublesex (dsx) transcription factor, and dsx expression is mostly limited to cells that give rise to sexually dimorphic traits. However, it is unknown how this mosaic of sexually dimorphic and monomorphic organs arises. Here, we characterize the cis-regulatory sequences that control dsx expression in the foreleg, which contains multiple types of sex-specific sensory organs. We find that separate modular enhancers are responsible for dsx expression in each sexually dimorphic organ. Expression of dsx in the sex comb is co-regulated by two enhancers with distinct spatial and temporal specificities that are separated by a genitalia-specific enhancer. The sex comb-specific enhancer from D. willistoni, a species that primitively lacks sex combs, is not active in the foreleg. Thus, the mosaic of sexually dimorphic and monomorphic organs depends on modular regulation of dsx transcription by dedicated cell type-specific enhancers.

Multi-enhancer transcriptional hubs confer phenotypic robustness

We previously showed in Drosophila melanogaster embryos that low-affinity Ultrabithorax (Ubx)-responsive shavenbaby (svb) enhancers drive expression using localized transcriptional environments and that active svb enhancers on different chromosomes tended to colocalize (Tsai et al., 2017). Here, we test the hypothesis that these multi-enhancer 'hubs' improve phenotypic resilience to stress by buffering against decreases in transcription factor concentrations and transcriptional output. Deleting a redundant enhancer from the svb locus led to reduced trichome numbers in embryos raised at elevated temperatures. Using high-resolution fluorescence microscopy, we observed lower Ubx concentration and transcriptional output in this deletion allele. Transcription sites of the full svb cis-regulatory region inserted into a different chromosome colocalized with the svb locus, increasing Ubx concentration, the transcriptional output of svb, and partially rescuing the phenotype. Thus, multiple enhancers could reinforce a local transcriptional hub to buffer against environmental stresses and genetic perturbations, providing a mechanism for phenotypical robustness.

Cloudy with a Chance of Insights: Context Dependent Gene Regulation and Implications for Evolutionary Studies

Research in various fields of evolutionary biology has shown that divergence in gene expression is a key driver for phenotypic evolution. An exceptional contribution of cis-regulatory divergence has been found to contribute to morphological diversification. In the light of these findings, the analysis of genome-wide expression data has become one of the central tools to link genotype and phenotype information on a more mechanistic level. However, in many studies, especially if general conclusions are drawn from such data, a key feature of gene regulation is often neglected. With our article, we want to raise awareness that gene regulation and thus gene expression is highly context dependent. Genes show tissue- and stage-specific expression. We argue that the regulatory context must be considered in comparative expression studies.