Evolution of cis-regulatory sequence and function in Diptera

Genome-Wide Identification, Expression Profiling, and Characterization of Cyclin-like Genes Reveal Their Role in the Fertility of the Diamondback Moth

Simple Summary

Cyclin genes are known as cell cycle regulators and play a significant role in the fertility of different organisms, including mice and insects. Until now, no study has been performed on the complete identification of the cyclin genes in insects. Here, we identified 21 cyclin genes in the diamondback moth (DBM) genome through a comprehensive genome-wide analysis and evaluated the gene structure, genomic location, and evolutionary relationship. Cis-regulatory elements and potential miRNA targeting the cyclin genes were also assessed. By analyzing the transcriptomic and RT-qPCR based expression profiling at different stages and tissues, we found that the majority of the cyclin genes were significantly expressed in the reproductive tissues. Moreover, RNAi-mediated characterization of PxCyc B1 showed its role in female fertility. The current study provides a basis for further evaluation of the cyclin genes, which may be used as a potential target for pest management programs.

Abstract

Cyclin-like genes are primarily considered as cell cycle regulators and have shown to be crucial for insect growth, development, reproduction, and fertility. However, no research has been performed on the cyclin-like genes in the diamondback moth (Plutella xylostella). Here, we identified the 21 cyclin genes in the genome of P. xylostella and clustered them into four groups. Most cyclin genes showed a well-maintained gene structure and motif distribution within the same group. The putative promoter regions of cyclin genes contained several transcription binding factors related to reproduction, along with growth and development. Furthermore, 16 miRNAs were identified targeting the 13 cyclin genes. Transcriptome and quantitative real-time PCR (qRT-PCR)-based expression profiling of cyclin-like genes at different stages and tissues were evaluated, revealing that 16 out of 21 cyclin genes were highly expressed in reproductive tissues of adult females and males. The Cyclin B1 gene (PxCyc B1) was only expressed in the ovary of the adult female and selected for the subsequent analysis. RNAi-mediated suppression of PxCyc B1 interrupted the external genitalia and length of the ovariole of female adults. Furthermore, the egg-laying capacity and hatching rate were also significantly decreased by suppressing the PxCyc B1, indicating the importance of cyclin genes in the reproduction and fertility of P. xylostella. The current study explained the detailed genome-wide analysis of cyclin-like genes in P. xylostella, which provided a basis for subsequent research to assess the roles of cyclin genes in reproduction, and the cyclin gene may be considered an effective target site to control this pest.

Expression Profile of Sorghum Genes and Cis-Regulatory Elements under Salt-Stress Conditions

Salinity stress is one of the most important abiotic stresses that causes great losses in crop production worldwide. Identifying the molecular mechanisms of salt resistance in sorghum will help develop salt-tolerant crops with high yields. Sorghum (Sorghum bicolor (L.) Moench) is one of the world’s four major grains and is known as a plant with excellent adaptability to salt stress. Among the various genotypes of sorghum, a Korean cultivar Nampungchal is also highly tolerant to salt. However, little is known about how Nampungchal responds to salt stress. In this study, we measured various physiological parameters, including Na⁺ and K⁺ contents, in leaves grown under saline conditions and investigated the expression patterns of differentially expressed genes (DEGs) using QuantSeq analysis. These DEG analyses revealed that genes up-regulated in a 150 mM NaCl treatment have various functions related to abiotic stresses, such as ERF and DREB. In addition, transcription factors such as ABA, WRKY, MYB, and bZip bind to the CREs region of sorghum and are involved in the regulation of various abiotic stress-responsive transcriptions, including salt stress. These findings may deepen our understanding of the mechanisms of salt tolerance in sorghum and other crops.

Distinct genetic architectures underlie divergent thorax, leg, and wing pigmentation between Drosophila elegans and D. gunungcola

Pigmentation divergence between Drosophila species has emerged as a model trait for studying the genetic basis of phenotypic evolution, with genetic changes contributing to pigmentation differences often mapping to genes in the pigment synthesis pathway and their regulators. These studies of Drosophila pigmentation have tended to focus on pigmentation changes in one body part for a particular pair of species, but changes in pigmentation are often observed in multiple body parts between the same pair of species. The similarities and differences of genetic changes responsible for divergent pigmentation in different body parts of the same species thus remain largely unknown. Here we compare the genetic basis of pigmentation divergence between Drosophila elegans and D. gunungcola in the wing, legs, and thorax. Prior work has shown that regions of the genome containing the pigmentation genes yellow and ebony influence the size of divergent male-specific wing spots between these two species. We find that these same two regions of the genome underlie differences in leg and thorax pigmentation; however, divergent alleles in these regions show differences in allelic dominance and epistasis among the three body parts. These complex patterns of inheritance can be explained by a model of evolution involving tissue-specific changes in the expression of Yellow and Ebony between D. elegans and D. gunungcola.

Intra-Strain Genetic Variation of Platyfish (Xiphophorus maculatus) Strains Determines Tumorigenic Trajectory

Xiphophorus interspecies hybrids represent a valuable model system to study heritable tumorigenesis, and the only model system that exhibits both spontaneous and inducible tumors. Types of tumorigenesis depend on the specific pedigree of the parental species, X. maculatus, utilized to produce interspecies hybrids. Although the ancestors of the two currently used X. maculatus parental lines, Jp163 A and Jp163 B, were originally siblings produced by the same mother, backcross interspecies hybrid progeny between X. hellerii and X. maculatus Jp163 A develop spontaneous melanoma initiating at the dorsal fin due to segregation of an oncogene and a regulator encoded by the X. maculatus genome, while the backcross hybrid progeny with X. hellerii or X. couchianus and Jp163 B exhibit melanoma on the flanks of their bodies, especially after treatment with ultraviolet light. Therefore, dissecting the genetic differences between these two closely related lines may lead to better understanding of functional molecular differences associated with tumorigenic mechanisms. For this purpose, comparative genomic analyses were undertaken to establish genetic variants between these two X. maculatus lines. Surprisingly, given the heritage of these two fish lines, we found genetic variants are clustered together in select chromosomal regions. Among these variants are non-synonymous mutations located in 381 genes. The non-random distribution of genetic variants between these two may highlight ancestral chromosomal recombination patterns that became fixed during subsequent inbreeding. Employing comparative transcriptomics, we also determined differences in the skin transcriptional landscape between the two lines. The genetic differences observed are associated with pathways highlighting fundamental cellular functions including inter-cellular and microenvironment-cellular interactions, and DNA repair. These results collectively lead to the conclusion that diverged functional genetic baselines are present between Jp163 A and B strains. Further, disruption of these fixed genetic baselines in the hybrids may give rise to spontaneous or inducible mechanisms of tumorigenesis.

Modular cis-regulatory logic of yellow gene expression in silkmoth larvae

Colour patterns in butterflies and moths are crucial traits for adaptation. Previous investigations have highlighted genes responsible for pigmentation (ie yellow and ebony). However, the mechanisms by which these genes are regulated in lepidopteran insects remain poorly understood. To elucidate this, molecular studies involving dipterans have largely analysed the cis-regulatory regions of pigmentation genes and have revealed cis-regulatory modularity. Here, we used well-developed transgenic techniques in Bombyx mori and demonstrated that cis-regulatory modularity controls tissue-specific expression of the yellow gene. We first identified which body parts are regulated by the yellow gene via black pigmentation. We then isolated three discrete regulatory elements driving tissue-specific gene expression in three regions of B. mori larvae. Finally, we found that there is no apparent sequence conservation of cis-regulatory regions between B. mori and Drosophila melanogaster, and no expression driven by the regulatory regions of one species when introduced into the other species. Therefore, the trans-regulatory landscapes of the yellow gene differ significantly between the two taxa. The results of this study confirm that lepidopteran species use cis-regulatory modules to control gene expression related to pigmentation, and represent a powerful cadre of transgenic tools for studying evolutionary developmental mechanisms.

Different Evolutionary Strategies To Conserve Chromatin Boundary Function in the Bithorax Complex

Chromatin boundary elements subdivide chromosomes in multicellular organisms into physically independent domains. In addition to this architectural function, these elements also play a critical role in gene regulation. Here we investigated the evolution of a Drosophila Bithorax complex boundary element called Fab-7, which is required for the proper parasegment specific expression of the homeotic Abd-B gene. Using a “gene” replacement strategy, we show that Fab-7 boundaries from two closely related species, D. erecta and D. yakuba, and a more distant species, D. pseudoobscura, are able to substitute for the melanogaster boundary. Consistent with this functional conservation, the two known Fab-7 boundary factors, Elba and LBC, have recognition sequences in the boundaries from all species. However, the strategies used for maintaining binding and function in the face of sequence divergence is different. The first is conventional, and depends upon conservation of the 8 bp Elba recognition sequence. The second is unconventional, and takes advantage of the unusually large and flexible sequence recognition properties of the LBC boundary factor, and the deployment of multiple LBC recognition elements in each boundary. In the former case, binding is lost when the recognition sequence is altered. In the latter case, sequence divergence is accompanied by changes in the number, relative affinity, and location of the LBC recognition elements.

Sensitivity of Allelic Divergence to Genomic Position: Lessons from the Drosophila tan Gene

To identify genetic variants underlying changes in phenotypes within and between species, researchers often utilize transgenic animals to compare the function of alleles in different genetic backgrounds. In Drosophila, targeted integration mediated by the ΦC31 integrase allows activity of alternative alleles to be compared at the same genomic location. By using the same insertion site for each transgene, position effects are generally assumed to be controlled for because both alleles are surrounded by the same genomic context. Here, we test this assumption by comparing the activity of tan alleles from two Drosophila species, D. americana and D. novamexicana, at five different genomic locations in D. melanogaster. We found that the relative effects of these alleles varied among insertion sites, with no difference in activity observed between them at two sites. One of these sites simply silenced both transgenes, but the other allowed expression of both alleles that was sufficient to rescue a mutant phenotype yet failed to reveal the functional differences between the two alleles. These results suggest that more than one insertion site should be used when comparing the activity of transgenes because failing to do so could cause functional differences between alleles to go undetected.

Functionally conserved enhancers with divergent sequences in distant vertebrates

BACKGROUND

To examine the contributions of sequence and function conservation in the evolution of enhancers, we systematically identified enhancers whose sequences are not conserved among distant groups of vertebrate species, but have homologous function and are likely to be derived from a common ancestral sequence. Our approach combined comparative genomics and epigenomics to identify potential enhancer sequences in the genomes of three groups of distantly related vertebrate species.

RESULTS

We searched for sequences that were conserved within groups of closely related species but not between groups of more distant species, and were associated with an epigenetic mark of enhancer activity. To facilitate inferring orthology between non-conserved sequences, we limited our search to introns whose orthology could be unambiguously established by mapping the bracketing exons. We show that a subset of these non-conserved but syntenic sequences from the mouse and zebrafish genomes have homologous functions in a zebrafish transgenic enhancer assay. The conserved expression patterns driven by these enhancers are probably associated with short transcription factor-binding motifs present in the divergent sequences.

CONCLUSIONS

We have identified numerous potential enhancers with divergent sequences but a conserved function. These results indicate that selection on function, rather than sequence, may be a common mode of enhancer evolution; evidence for selection at the sequence level is not a necessary criterion to define a gene regulatory element.

cis-regulatory analysis of the Drosophila pdm locus reveals a diversity of neural enhancers

Background

One of the major challenges in developmental biology is to understand the regulatory events that generate neuronal diversity. During Drosophila embryonic neural lineage development, cellular temporal identity is established in part by a transcription factor (TF) regulatory network that mediates a cascade of cellular identity decisions. Two of the regulators essential to this network are the POU-domain TFs Nubbin and Pdm-2, encoded by adjacent genes collectively known as pdm. The focus of this study is the discovery and characterization of cis-regulatory DNA that governs their expression.

Results

Phylogenetic footprinting analysis of a 125 kb genomic region that spans the pdm locus identified 116 conserved sequence clusters. To determine which of these regions function as cis-regulatory enhancers that regulate the dynamics of pdm gene expression, we tested each for in vivo enhancer activity during embryonic development and postembryonic neurogenesis. Our screen revealed 77 unique enhancers positioned throughout the noncoding region of the pdm locus. Many of these activated neural-specific gene expression during different developmental stages and many drove expression in overlapping patterns. Sequence comparisons of functionally related enhancers that activate overlapping expression patterns revealed that they share conserved elements that can be predictive of enhancer behavior. To facilitate data accessibility, the results of our analysis are catalogued in cisPatterns, an online database of the structure and function of these and other Drosophila enhancers.

Conclusions

These studies reveal a diversity of modular enhancers that most likely regulate pdm gene expression during embryonic and adult development, highlighting a high level of temporal and spatial expression specificity. In addition, we discovered clusters of functionally related enhancers throughout the pdm locus. A subset of these enhancers share conserved elements including sequences that correspond to known TF DNA binding sites. Although comparative analysis of the nubbin and pdm-2 encoding sequences indicate that these two genes most likely arose from a duplication event, we found only partial evidence of sequence duplication between their enhancers, suggesting that after the putative duplication their cis-regulatory DNA diverged at a higher rate than their coding sequences.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1897-2) contains supplementary material, which is available to authorized users.

The achaete-scute complex in Diptera: patterns of noncoding sequence evolution

The achaete‐scute complex (AS‐C) has been a useful paradigm for the study of pattern formation and its evolution. achaete‐scute genes have duplicated and evolved distinct expression patterns during the evolution of cyclorraphous Diptera. Are the expression patterns in different species driven by conserved regulatory elements? If so, when did such regulatory elements arise? Here, we have sequenced most of the AS‐C of the fly Calliphora vicina (including the genes achaete, scute and lethal of scute) to compare noncoding sequences with known cis‐regulatory sequences in Drosophila. The organization of the complex is conserved with respect to Drosophila species. There are numerous small stretches of conserved noncoding sequence that, in spite of high sequence turnover, display binding sites for known transcription factors. Synteny of the blocks of conserved noncoding sequences is maintained suggesting not only conservation of the position of regulatory elements but also an origin prior to the divergence between these two species. We propose that some of these enhancers originated by duplication with their target genes.

Gain of cis-regulatory activities underlies novel domains of wingless gene expression in Drosophila

Changes in gene expression during animal development are largely responsible for the evolution of morphological diversity. However, the genetic and molecular mechanisms responsible for the origins of new gene-expression domains have been difficult to elucidate. Here, we sought to identify molecular events underlying the origins of three novel features of wingless (wg) gene expression that are associated with distinct pigmentation patterns in Drosophila guttifera. We compared the activity of cis-regulatory sequences (enhancers) across the wg locus in D. guttifera and Drosophila melanogaster and found strong functional conservation among the enhancers that control similar patterns of wg expression in larval imaginal discs that are essential for appendage development. For pupal tissues, however, we found three novel wg enhancer activities in D. guttifera associated with novel domains of wg expression, including two enhancers located surprisingly far away in an intron of the distant Wnt10 gene. Detailed analysis of one enhancer (the vein-tip enhancer) revealed that it overlapped with a region controlling wg expression in wing crossveins (crossvein enhancer) in D. guttifera and other species. Our results indicate that one novel domain of wg expression in D. guttifera wings evolved by co-opting pre-existing regulatory sequences governing gene activity in the developing wing. We suggest that the modification of existing enhancers is a common path to the evolution of new gene-expression domains and enhancers.

The role of self-organization in developmental evolution

In developmental and evolutionary biology, particular emphasis has been given to the relationship between transcription factors and the cognate cis-regulatory elements of their target genes. These constitute the gene regulatory networks that control expression and are assumed to causally determine the formation of structures and body plans. Comparative analysis has, however, established a broad sequence homology among species that nonetheless display quite different anatomies. Transgenic experiments have also confirmed that many developmentally important elements are, in fact, functionally interchangeable. Although dependent upon the appropriate degree of gene expression, the actual construction of specific structures appears not directly linked to the functions of gene products alone. Instead, the self-formation of complex patterns, due in large part to epigenetic and non-genetic determinants, remains a persisting theme in the study of ontogeny and regenerative medicine. Recent evidence indeed points to the existence of a self-organizing process, operating through a set of intrinsic rules and forces, which imposes coordination and a holistic order upon cells and tissue. This has been repeatedly demonstrated in experiments on regeneration as well as in the autonomous formation of structures in vitro. The process cannot be wholly attributed to the functional outcome of protein-protein interactions or to concentration gradients of diffusible chemicals. This phenomenon is examined here along with some of the methodological and theoretical approaches that are now used in understanding the causal basis for self-organization in development and its evolution.

Conserved boundary elements from the Hox complex of mosquito, Anopheles gambiae

The conservation of hox genes as well as their genomic organization across the phyla suggests that this system of anterior–posterior axis formation arose early during evolution and has come under strong selection pressure. Studies in the split Hox cluster of Drosophila have shown that proper expression of hox genes is dependent on chromatin domain boundaries that prevent inappropriate interactions among different types of cis-regulatory elements. To investigate whether boundary function and their role in regulation of hox genes is conserved in insects with intact Hox clusters, we used an algorithm to locate potential boundary elements in the Hox complex of mosquito, Anopheles gambiae. Several potential boundary elements were identified that could be tested for their functional conservation. Comparative analysis revealed that like Drosophila, the bithorax region in A. gambiae contains an extensive array of boundaries and enhancers organized into domains. We analysed a subset of candidate boundary elements and show that they function as enhancer blockers in Drosophila. The functional conservation of boundary elements from mosquito in fly suggests that regulation of hox genes involving chromatin domain boundaries is an evolutionary conserved mechanism and points to an important role of such elements in key developmentally regulated loci.

Novel genes from formation to function

The study of the evolution of novel genes generally focuses on the formation of new coding sequences. However, equally important in the evolution of novel functional genes are the formation of regulatory regions that allow the expression of the genes and the effects of the new genes in the organism as well. Herein, we discuss the current knowledge on the evolution of novel functional genes, and we examine in more detail the youngest genes discovered. We examine the existing data on a very recent and rapidly evolving cluster of duplicated genes, the Sdic gene cluster. This cluster of genes is an excellent model for the evolution of novel genes, as it is very recent and may still be in the process of evolving.

Rapid evolutionary rewiring of a structurally constrained eye enhancer

BACKGROUND

Enhancers are genomic cis-regulatory sequences that integrate spatiotemporal signals to control gene expression. Enhancer activity depends on the combination of bound transcription factors as well as-in some cases-the arrangement and spacing of binding sites for these factors. Here, we examine evolutionary changes to the sequence and structure of sparkling, a Notch/EGFR/Runx-regulated enhancer that activates the dPax2 gene in cone cells of the developing Drosophila eye.

RESULTS

Despite functional and structural constraints on its sequence, sparkling has undergone major reorganization in its recent evolutionary history. Our data suggest that the relative strengths of the various regulatory inputs into sparkling change rapidly over evolutionary time, such that reduced input from some factors is compensated by increased input from different regulators. These gains and losses are at least partly responsible for the changes in enhancer structure that we observe. Furthermore, stereotypical spatial relationships between certain binding sites ("grammar elements") can be identified in all sparkling orthologs-although the sites themselves are often recently derived. We also find that low binding affinity for the Notch-regulated transcription factor Su(H), a conserved property of sparkling, is required to prevent ectopic responses to Notch in noncone cells.

CONCLUSIONS

Rapid DNA sequence turnover does not imply either the absence of critical cis-regulatory information or the absence of structural rules. Our findings demonstrate that even a severely constrained cis-regulatory sequence can be significantly rewired over a short evolutionary timescale.

Genome-wide identification of conserved regulatory function in diverged sequences

Plasticity of gene regulatory encryption can permit DNA sequence divergence without loss of function. Functional information is preserved through conservation of the composition of transcription factor binding sites (TFBS) in a regulatory element. We have developed a method that can accurately identify pairs of functional noncoding orthologs at evolutionarily diverged loci by searching for conserved TFBS arrangements. With an estimated 5% false-positive rate (FPR) in approximately 3000 human and zebrafish syntenic loci, we detected approximately 300 pairs of diverged elements that are likely to share common ancestry and have similar regulatory activity. By analyzing a pool of experimentally validated human enhancers, we demonstrated that 7/8 (88%) of their predicted functional orthologs retained in vivo regulatory control. Moreover, in 5/7 (71%) of assayed enhancer pairs, we observed concordant expression patterns. We argue that TFBS composition is often necessary to retain and sufficient to predict regulatory function in the absence of overt sequence conservation, revealing an entire class of functionally conserved, evolutionarily diverged regulatory elements that we term "covert."

Morphological evolution caused by many subtle-effect substitutions in regulatory DNA

Morphology evolves often through changes in developmental genes, but the causal mutations, and their effects, remain largely unknown. The evolution of naked cuticle—rather than trichomes—on larvae of Drosophila sechellia resulted from changes in five transcriptional enhancers of shavenbaby, a gene encoding a transcription factor that governs trichome morphogenesis. Here we show that the function of one of these enhancers evolved through multiple single nucleotide substitutions that altered both the timing and level of shavenbaby expression. The consequences of these nucleotide substitutions on larval morphology were quantified with a novel functional assay. We found that each substitution had a relatively small phenotypic effect, and that many nucleotide changes account for this large morphological difference. In addition, we observed that the substitutions displayed non-additive effects to generate a large phenotypic change. These data provide unprecedented resolution of the phenotypic effects of substitutions and show how individual nucleotide changes in a transcriptional enhancer have caused morphological evolution.

Ongoing purifying selection on intergenic spacers in group A streptococcus

Bacterial intergenic spacers are non-coding genomic regions enriched with cis-regulatory elements for gene expression. A population genetics approach was used to investigate the evolutionary force shaping the genetic diversity of intergenic spacers among 13 genomes of group A streptococcus (GAS). Analysis of 590 genes and their linked 5' intergenic spacers showed reduced nucleotide diversity in spacers compared to synonymous nucleotide diversity in protein-coding regions, suggestive of past purifying selection on spacers. Certain spacers showed elevated nucleotide diversity indicative of past homologous recombination with divergent genotypes. In addition, analysis of the difference between mean nucleotide difference and number of segregating sites showed evidence of an excess of rare variants both at nonsynonymous sites in genes and at sites in spacers, which is evidence that there are numerous slightly deleterious variants in GAS populations with potential effects on both protein sequences and gene expression.

Nomadic enhancers: tissue-specific cis-regulatory elements of yellow have divergent genomic positions among Drosophila species

cis-regulatory DNA sequences known as enhancers control gene expression in space and time. They are central to metazoan development and are often responsible for changes in gene regulation that contribute to phenotypic evolution. Here, we examine the sequence, function, and genomic location of enhancers controlling tissue- and cell-type specific expression of the yellow gene in six Drosophila species. yellow is required for the production of dark pigment, and its expression has evolved largely in concert with divergent pigment patterns. Using Drosophila melanogaster as a transgenic host, we examined the expression of reporter genes in which either 5' intergenic or intronic sequences of yellow from each species controlled the expression of Green Fluorescent Protein. Surprisingly, we found that sequences controlling expression in the wing veins, as well as sequences controlling expression in epidermal cells of the abdomen, thorax, and wing, were located in different genomic regions in different species. By contrast, sequences controlling expression in bristle-associated cells were located in the intron of all species. Differences in the precise pattern of spatial expression within the developing epidermis of D. melanogaster transformants usually correlated with adult pigmentation in the species from which the cis-regulatory sequences were derived, which is consistent with cis-regulatory evolution affecting yellow expression playing a central role in Drosophila pigmentation divergence. Sequence comparisons among species favored a model in which sequential nucleotide substitutions were responsible for the observed changes in cis-regulatory architecture. Taken together, these data demonstrate frequent changes in yellow cis-regulatory architecture among Drosophila species. Similar analyses of other genes, combining in vivo functional tests of enhancer activity with in silico comparative genomics, are needed to determine whether the pattern of regulatory evolution we observed for yellow is characteristic of genes with rapidly evolving expression patterns.

The multiple origins of aluminium resistance in hexaploid wheat include Aegilops tauschii and more recent cis mutations to TaALMT1

Acid soils limit plant production worldwide because their high concentrations of soluble aluminium cations (Al(3+) ) inhibit root growth. Major food crops such as wheat (Triticum aestivum L.) have evolved mechanisms to resist Al(3+) toxicity, thus enabling wider distribution. The origins of Al(3+) resistance in wheat are perplexing because all progenitors of this hexaploid species are reportedly sensitive to Al(3+) stress. The large genotypic variation for Al(3+) resistance in wheat is largely controlled by expression of an anion channel, TaALMT1, which releases malate anions from the root apices. A current hypothesis proposes that the malate anions protect this sensitive growth zone by binding to Al(3+) in the apoplasm. We investigated the evolution of this trait in wheat, and demonstrated that it has multiple independent origins that enhance Al(3+) resistance by increasing TaALMT1 expression. One origin is likely to be Aegilops tauschii while other origins occurred more recently from a series of cis mutations that have generated tandemly repeated elements in the TaALMT1 promoter. We generated transgenic plants to directly compare these promoter alleles and demonstrate that the tandemly repeated elements act to enhance gene expression. This study provides an example from higher eukaryotes in which perfect tandem repeats are linked with transcriptional regulation and phenotypic change in the context of evolutionary adaptation to a major abiotic stress.

The genetic basis of evolutionary change in gene expression levels

The regulation of gene expression is an important determinant of organismal phenotype and evolution. However, the widespread recognition of this fact occurred long after the synthesis of evolution and genetics. Here, we give a brief sketch of thoughts regarding gene regulation in the history of evolution and genetics. We then review the development of genome-wide studies of gene regulatory variation in the context of the location and mode of action of the causative genetic changes. In particular, we review mapping of the genetic basis of expression variation through expression quantitative trait locus studies and measuring the cis/trans component of expression variation in allele-specific expression studies. We conclude by proposing a systematic integration of ideas that combines global mapping studies, cis/trans tests and modern population genetics methodologies, in order to directly estimate the forces acting on regulatory variation within and between species.

Wingless and aristaless2 define a developmental ground plan for moth and butterfly wing pattern evolution

Butterfly wing patterns have long been a favorite system for studying the evolutionary radiation of complex morphologies. One of the key characteristics of the system is that wing patterns are based on a highly conserved ground plan of pattern homologies. In fact, the evolution of lepidopteran wing patterns is proposed to have occurred through the repeated gain, loss, and modification of only a handful of serially repeated elements. In this study, we examine the evolution and development of stripe wing pattern elements. We show that expression of the developmental morphogen wingless (wg) is associated with early determination of the major basal (B), discal (DI and DII), and marginal (EI) stripe patterns in a broad sampling of Lepidoptera, suggesting homology of these pattern elements across moths and butterflies. We describe for the first time a novel Lepidoptera-specific homeobox gene, aristaless2 (al2), which precedes wg expression during the early determination of DII stripe patterns. We show that al2 was derived from a tandem duplication of the aristaless gene, whereupon it underwent a rapid coding and cis-regulatory divergence relative to its more conserved paralog aristaless1 (al1), which retained an ancestral expression pattern. The al2 stripe expression domain evolutionarily preceded the appearance of the DII pattern elements in multiple lineages, leading us to speculate that al2 represented preexisting positional information that may have facilitated DII evolution via a developmental drive mechanism. In contrast to butterfly eyespot patterns, which are often cited as a key example of developmental co-option of preexisting developmental genes, this study provides an example where the origin of a major color pattern element is associated with the evolution of a novel lepidopteran homeobox gene.

Phylogeny disambiguates the evolution of heat-shock cis-regulatory elements in Drosophila

Heat-shock genes have a well-studied control mechanism for their expression that is mediated through cis-regulatory motifs known as heat-shock elements (HSEs). The evolution of important features of this control mechanism has not been investigated in detail, however. Here we exploit the genome sequencing of multiple Drosophila species, combined with a wealth of available information on the structure and function of HSEs in D. melanogaster, to undertake this investigation. We find that in single-copy heat shock genes, entire HSEs have evolved or disappeared 14 times, and the phylogenetic approach bounds the timing and direction of these evolutionary events in relation to speciation. In contrast, in the multi-copy gene Hsp70, the number of HSEs is nearly constant across species. HSEs evolve in size, position, and sequence within heat-shock promoters. In turn, functional significance of certain features is implicated by preservation despite this evolutionary change; these features include tail-to-tail arrangements of HSEs, gapped HSEs, and the presence or absence of entire HSEs. The variation among Drosophila species indicates that the cis-regulatory encoding of responsiveness to heat and other stresses is diverse. The broad dimensions of variation uncovered are particularly important as they suggest a substantial challenge for functional studies.

Structural rules and complex regulatory circuitry constrain expression of a Notch- and EGFR-regulated eye enhancer

Enhancers integrate spatiotemporal information to generate precise patterns of gene expression. How complex is the regulatory logic of a typical developmental enhancer, and how important is its internal organization? Here, we examine in detail the structure and function of sparkling, a Notch- and EGFR/MAPK-regulated, cone cell-specific enhancer of the Drosophila Pax2 gene, in vivo. In addition to its 12 previously identified protein-binding sites, sparkling is densely populated with previously unmapped regulatory sequences, which interact in complex ways to control gene expression. One segment is essential for activation at a distance, yet dispensable for other activation functions and for cell type patterning. Unexpectedly, rearranging sparkling's regulatory sites converts it into a robust photoreceptor-specific enhancer. Our results show that a single combination of regulatory inputs can encode multiple outputs, and suggest that the enhancer's organization determines the correct expression pattern by facilitating certain short-range regulatory interactions at the expense of others.

Evolution of cis-regulatory sequences in Drosophila

Cross-species comparison is an emerging paradigm for identifying cis-regulatory sequences and understanding their function and evolution. In this chapter, we review probabilistic models of evolution of transcription factor binding sites, which provide the theoretical basis for a number of new bioinformatics tools for comparative sequence analysis. We illustrate how important functional and evolutionary insights on binding site gain and loss can be acquired through sequence comparison. This includes the observation that binding site turnover follows a molecular clock and that its rate correlates with the strength of binding sites and the presence of other sites in the neighborhood. We also comment on emerging trends that go beyond individual binding sites to a more holistic study of regulatory evolution. We point out common technical challenges, such as reliable sequence alignment and binding site prediction, when doing comparative regulatory sequence analysis and note some potential solutions thereof.

Functional evolution of cis-regulatory modules at a homeotic gene in Drosophila

It is a long-held belief in evolutionary biology that the rate of molecular evolution for a given DNA sequence is inversely related to the level of functional constraint. This belief holds true for the protein-coding homeotic (Hox) genes originally discovered in Drosophila melanogaster. Expression of the Hox genes in Drosophila embryos is essential for body patterning and is controlled by an extensive array of cis-regulatory modules (CRMs). How the regulatory modules functionally evolve in different species is not clear. A comparison of the CRMs for the Abdominal-B gene from different Drosophila species reveals relatively low levels of overall sequence conservation. However, embryonic enhancer CRMs from other Drosophila species direct transgenic reporter gene expression in the same spatial and temporal patterns during development as their D. melanogaster orthologs. Bioinformatic analysis reveals the presence of short conserved sequences within defined CRMs, representing gap and pair-rule transcription factor binding sites. One predicted binding site for the gap transcription factor KRUPPEL in the IAB5 CRM was found to be altered in Superabdominal (Sab) mutations. In Sab mutant flies, the third abdominal segment is transformed into a copy of the fifth abdominal segment. A model for KRUPPEL-mediated repression at this binding site is presented. These findings challenge our current understanding of the relationship between sequence evolution at the molecular level and functional activity of a CRM. While the overall sequence conservation at Drosophila CRMs is not distinctive from neighboring genomic regions, functionally critical transcription factor binding sites within embryonic enhancer CRMs are highly conserved. These results have implications for understanding mechanisms of gene expression during embryonic development, enhancer function, and the molecular evolution of eukaryotic regulatory modules.

Unravelling cis-regulatory elements in the genome of the smallest photosynthetic eukaryote: phylogenetic footprinting in Ostreococcus

We used a phylogenetic footprinting approach, adapted to high levels of divergence, to estimate the level of constraint in intergenic regions of the extremely gene dense Ostreococcus algae genomes (Chlorophyta, Prasinophyceae). We first benchmarked our method against the Saccharomyces sensu stricto genome data and found that the proportion of conserved non-coding sites was consistent with those obtained with methods using calibration by the neutral substitution rate. We then applied our method to the complete genomes of Ostreococcus tauri and O. lucimarinus, which are the most divergent species from the same genus sequenced so far. We found that 77% of intergenic regions in Ostreococcus still contain some phylogenetic footprints, as compared to 88% for Saccharomyces, corresponding to an average rate of constraint on intergenic region of 17% and 30%, respectively. A comparison with some known functional cis-regulatory elements enabled us to investigate whether some transcriptional regulatory pathways were conserved throughout the green lineage. Strikingly, the size of the phylogenetic footprints depends on gene orientation of neighboring genes, and appears to be genus-specific. In Ostreococcus, 5' intergenic regions contain four times more conserved sites than 3' intergenic regions, whereas in yeast a higher frequency of constrained sites in intergenic regions between genes on the same DNA strand suggests a higher frequency of bidirectional regulatory elements. The phylogenetic footprinting approach can be used despite high levels of divergence in the ultrasmall Ostreococcus algae, to decipher structure of constrained regulatory motifs, and identify putative regulatory pathways conserved within the green lineage.

Homology of dipteran bristles and lepidopteran scales: requirement for the Bombyx mori achaete-scute homologue ASH2

Lepidopteran wing scales and Drosophila bristles are considered homologous structures on the basis of the similarities in their cell lineages. However, the molecular mechanisms underlying scale development are essentially unknown as analysis of gene function in Lepidoptera is sorely limited. In this study, we used the Bombyx mori mutant scaleless (sl), which displays a nearly complete loss of wing scales, to explore the mechanism of lepidopteran wing-scale formation. We found that Bm-ASH2, one of four Bombyx achaete-scute homologs, is highly expressed in early pupal wings of wild-type silkworms, but its expression is severely reduced in sl pupal wings. Through molecular characterization of the mutant locus using luciferase and gel shift assays, genetic analysis of recombining populations, and in vivo rescue experiments, we provide evidence that a 26-bp deletion within the Bm-ASH2 promoter is closely linked to the sl locus and leads to loss of Bm-ASH2 expression and the scaleless-wings phenotype. Thus, the Bm-ASH2 appears to play a critical role in scale formation in B. mori. This finding supports the proposed homology of lepidopteran scales and dipteran bristles and provides evidence for conservation of the genetic pathway in scale/bristle development at the level of gene function.

Rapid evolution of sex pheromone-producing enzyme expression in Drosophila

A wide range of organisms use sex pheromones to communicate with each other and to identify appropriate mating partners. While the evolution of chemical communication has been suggested to cause sexual isolation and speciation, the mechanisms that govern evolutionary transitions in sex pheromone production are poorly understood. Here, we decipher the molecular mechanisms underlying the rapid evolution in the expression of a gene involved in sex pheromone production in Drosophilid flies. Long-chain cuticular hydrocarbons (e.g., dienes) are produced female-specifically, notably via the activity of the desaturase DESAT-F, and are potent pheromones for male courtship behavior in Drosophila melanogaster. We show that across the genus Drosophila, the expression of this enzyme is correlated with long-chain diene production and has undergone an extraordinary number of evolutionary transitions, including six independent gene inactivations, three losses of expression without gene loss, and two transitions in sex-specificity. Furthermore, we show that evolutionary transitions from monomorphism to dimorphism (and its reversion) in desatF expression involved the gain (and the inactivation) of a binding-site for the sex-determination transcription factor, DOUBLESEX. In addition, we documented a surprising example of the gain of particular cis-regulatory motifs of the desatF locus via a set of small deletions. Together, our results suggest that frequent changes in the expression of pheromone-producing enzymes underlie evolutionary transitions in chemical communication, and reflect changing regimes of sexual selection, which may have contributed to speciation among Drosophila.

Reverse-engineering a transcriptional enhancer: a case study in Drosophila

Enhancers, or cis-regulatory elements, are the principal determinants of spatiotemporal patterning of gene expression. For reasons of clinical and research utility, it is desirable to build customized enhancers that drive novel gene expression patterns, but currently, we largely rely on "found" genomic elements. Synthetic enhancers, assembled from transcription factor binding sites taken from natural signal-regulated enhancers, generally fail to behave like their wild-type counterparts when placed in transgenic animals, suggesting that important aspects of enhancer function are still unexplored. As a step toward the creation of a truly synthetic regulatory element, we have undertaken an extensive structure-function study of an enhancer of the Drosophila decapentaplegic (dpp) gene that drives expression in the developing visceral mesoderm (VM). Although considerable past efforts have been made to dissect the dppVM enhancer, transgenic experiments presented here indicate that its activity cannot be explained by the known regulators alone. dppVM contains multiple, previously uncharacterized, regulatory sites, some of which exhibit functional redundancy. The results presented here suggest that even the best-studied enhancers must be further dissected before they can be fully understood, and before faithful synthetic elements based on them can be created. Implications for developmental genetics, mathematical modeling, and therapeutic applications are discussed.

Spatial regulation of Corazonin neuropeptide expression requires multiple cis-acting elements in Drosophila melanogaster

Although most invertebrate neuropeptide-encoding genes display distinct expression patterns in the central nervous system (CNS), the molecular mechanisms underlying spatial regulation of the neuropeptide genes are largely unknown. Expression of the neuropeptide Corazonin (Crz) is limited to only 24 neurons in the larval CNS of Drosophila melanogaster, and these neurons have been categorized into three groups, namely, DL, DM, and vCrz. To identify cis-regulatory elements that control transcription of Crz in each neuronal group, reporter gene expression patterns driven by various 5' flanking sequences of Crz were analyzed to assess their promoter activities in the CNS. We show that the 504-bp 5' upstream sequence is the shortest promoter directing reporter activities in all Crz neurons. Further dissection of this sequence revealed two important regions responsible for group specificity: -504::-419 for DM expression and -380::-241 for DL and vCrz expression. The latter region is further subdivided into three sites (proximal, center, and distal), in which any combinations of the two are sufficient for DL expression, whereas both proximal and distal sites are required for vCrz expression. Interestingly, the TATA box does not play a role in Crz transcription in most neurons. We also show that a 434-bp 5' upstream sequence of the D. virilis Crz gene, when introduced into the D. melanogaster genome, drives reporter expression in the DL and vCrz neurons, suggesting that regulatory mechanisms for Crz expression in at least two such neuronal groups are conserved between the two species.

Functional diversification of sonic hedgehog paralog enhancers identified by phylogenomic reconstruction

BACKGROUND

Cis-regulatory modules of developmental genes are targets of evolutionary changes that underlie the morphologic diversity of animals. Little is known about the 'grammar' of interactions between transcription factors and cis-regulatory modules and therefore about the molecular mechanisms that underlie changes in these modules, particularly after gene and genome duplications. We investigated the ar-C midline enhancer of sonic hedgehog (shh) orthologs and paralogs from distantly related vertebrate lineages, from fish to human, including the basal vertebrate Latimeria menadoensis.

RESULTS

We demonstrate that the sonic hedgehog a (shha) paralogs sonic hedgehog b (tiggy winkle hedgehog; shhb) genes of fishes have a modified ar-C enhancer, which specifies a diverged function at the embryonic midline. We have identified several conserved motifs that are indicative of putative transcription factor binding sites by local alignment of ar-C enhancers of numerous vertebrate sequences. To trace the evolutionary changes among paralog enhancers, phylogenomic reconstruction was carried out and lineage-specific motif changes were identified. The relation between motif composition and observed developmental differences was evaluated through transgenic functional analyses. Altering and exchanging motifs between paralog enhancers resulted in reversal of enhancer specificity in the floor plate and notochord. A model reconstructing enhancer divergence during vertebrate evolution was developed.

CONCLUSION

Our model suggests that the identified motifs of the ar-C enhancer function as binary switches that are responsible for specific activity between midline tissues, and that these motifs are adjusted during functional diversification of paralogs. The unraveled motif changes can also account for the complex interpretation of activator and repressor input signals within a single enhancer.

Evidence of spatially varying selection acting on four chromatin-remodeling loci in Drosophila melanogaster

The packaging of DNA into proper chromatin structure contributes to transcriptional regulation. This packaging is environment sensitive, yet its role in adaptation to novel environmental conditions is completely unknown. We set out to identify candidate chromatin-remodeling loci that are differentiated between tropical and temperate populations in Drosophila melanogaster, an ancestrally equatorial African species that has recently colonized temperate environments around the world. Here we describe sequence variation at seven such chromatin-remodeling loci, four of which (chd1, ssrp, chm, and glu) exhibit strong differentiation between tropical and temperate populations. An in-depth analysis of chm revealed sequence differentiation restricted to a small portion of the gene, as well as evidence of clinal variation along the east coasts of both the United States and Australia. The functions of chd1, chm, ssrp, and glu point to several novel hypotheses for the role of chromatin-based transcriptional regulation in adaptation to a novel environment. Specifically, both stress-induced transcription and developmental homeostasis emerge as potential functional targets of environment-dependent selection.

Evolution of genes and genomes on the Drosophila phylogeny

Comparative analysis of multiple genomes in a phylogenetic framework dramatically improves the precision and sensitivity of evolutionary inference, producing more robust results than single-genome analyses can provide. The genomes of 12 Drosophila species, ten of which are presented here for the first time (sechellia, simulans, yakuba, erecta, ananassae, persimilis, willistoni, mojavensis, virilis and grimshawi), illustrate how rates and patterns of sequence divergence across taxa can illuminate evolutionary processes on a genomic scale. These genome sequences augment the formidable genetic tools that have made Drosophila melanogaster a pre-eminent model for animal genetics, and will further catalyse fundamental research on mechanisms of development, cell biology, genetics, disease, neurobiology, behaviour, physiology and evolution. Despite remarkable similarities among these Drosophila species, we identified many putatively non-neutral changes in protein-coding genes, non-coding RNA genes, and cis-regulatory regions. These may prove to underlie differences in the ecology and behaviour of these diverse species.

Incorporating evolution of transcription factor binding sites into annotated alignments

Identifying transcription factor binding sites (TFBSs) is essential to elucidate putative regulatory mechanisms. A common strategy is to combine cross-species conservation with single sequence TFBS annotation to yield "conserved TFBSs". Most current methods in this field adopt a multi-step approach that segregates the two aspects. Again, it is widely accepted that the evolutionary dynamics of binding sites differ from those of the surrounding sequence. Hence, it is desirable to have an approach that explicitly takes this factor into account. Although a plethora of approaches have been proposed for the prediction of conserved TFBSs, very few explicitly model TFBS evolutionary properties, while additionally being multi-step. Recently, we introduced a novel approach to simultaneously align and annotate conserved TFBSs in a pair of sequences. Building upon the standard Smith-Waterman algorithm for local alignments, SimAnn introduces additional states for profiles to output extended alignments or annotated alignments. That is, alignments with parts annotated as gaplessly aligned TFBSs (pair-profile hits)are generated. Moreover,the pair- profile related parameters are derived in a sound statistical framework. In this article, we extend this approach to explicitly incorporate evolution of binding sites in the SimAnn framework. We demonstrate the extension in the theoretical derivations through two position-specific evolutionary models, previously used for modelling TFBS evolution. In a simulated setting, we provide a proof of concept that the approach works given the underlying assumptions,as compared to the original work. Finally, using a real dataset of experimentally verified binding sites in human-mouse sequence pairs,we compare the new approach (eSimAnn) to an existing multi-step tool that also considers TFBS evolution. Although it is widely accepted that binding sites evolve differently from the surrounding sequences, most comparative TFBS identification methods do not explicitly consider this.Additionally, prediction of conserved binding sites is carried out in a multi-step approach that segregates alignment from TFBS annotation. In this paper, we demonstrate how the simultaneous alignment and annotation approach of SimAnn can be further extended to incorporate TFBS evolutionary relationships. We study how alignments and binding site predictions interplay at varying evolutionary distances and for various profile qualities.

Evaluating the role of natural selection in the evolution of gene regulation

Surveys of gene expression reveal extensive variability both within and between a wide range of species. Compelling cases have been made for adaptive changes in gene regulation, but the proportion of expression divergence attributable to natural selection remains unclear. Distinguishing adaptive changes driven by positive selection from neutral divergence resulting from mutation and genetic drift is critical for understanding the evolution of gene expression. Here, we review the various methods that have been used to test for signs of selection in genomic expression data. We also discuss properties of regulatory systems relevant to neutral models of gene expression. Despite some potential caveats, published studies provide considerable evidence for adaptive changes in gene expression. Future challenges for studies of regulatory evolution will be to quantify the frequency of adaptive changes, identify the genetic basis of expression divergence and associate changes in gene expression with specific organismal phenotypes.

The stars and stripes of animal bodies: evolution of regulatory elements mediating pigment and bristle patterns in Drosophila

Evolution has generated enormous morphological diversity in animals and one of the genetic processes that might have contributed to this is evolution of the cis-regulatory sequences responsible for the temporal and spatial expression of genes regulating embryonic development. This could be particularly relevant to pleiotropic genes with multiple independently acting regulatory modules. Loss or gain of modules enables altered expression without loss of other functions. Here I focus on recent studies correlating differences in morphological traits between related species of Drosophila to changes in cis-regulatory sequences. They show that ancestral regulatory modules have evolved to mediate different transcriptional outputs and suggest that evolution of cis-regulatory sequences might reflect a general mechanism driving evolutionary change.

Emerging principles of regulatory evolution

Understanding the genetic and molecular mechanisms governing the evolution of morphology is a major challenge in biology. Because most animals share a conserved repertoire of body-building and -patterning genes, morphological diversity appears to evolve primarily through changes in the deployment of these genes during development. The complex expression patterns of developmentally regulated genes are typically controlled by numerous independent cis-regulatory elements (CREs). It has been proposed that morphological evolution relies predominantly on changes in the architecture of gene regulatory networks and in particular on functional changes within CREs. Here, we discuss recent experimental studies that support this hypothesis and reveal some unanticipated features of how regulatory evolution occurs. From this growing body of evidence, we identify three key operating principles underlying regulatory evolution, that is, how regulatory evolution: (i) uses available genetic components in the form of preexisting and active transcription factors and CREs to generate novelty; (ii) minimizes the penalty to overall fitness by introducing discrete changes in gene expression; and (iii) allows interactions to arise among any transcription factor and downstream CRE. These principles endow regulatory evolution with a vast creative potential that accounts for both relatively modest morphological differences among closely related species and more profound anatomical divergences among groups at higher taxonomical levels.

Tissue-specific activity of the blind mole rat and the two nucleotide-mutated mouse alphaB-crystallin promoter in transgenic mice

The alphaB-crystallin and HspB2 genes are located approximately 0.9 kb apart in a head-to-head arrangement in mammals. Previous experiments have shown that a truncated -668/+45 alphaB-crystallin enhancer/promoter fragment from blind mole rats (Spalax ehrenbergi), which have nonfunctional lenses, lacks lens activity and has enhanced muscle activity in transgenic mice. Here we show that the full-length mole rat alphaB-crystallin intergenic region behaves similarly in transgenic mice. A two-nucleotide mutation ((-273)CA-->G) in the mouse alphaB-crystallin enhancer/promoter fragment mimicking the wild-type mole rat sequence functionally converted the mouse promoter fragment to that of the wild-type mole rat promoter when tested in transgenic mice. The reciprocal mutation in the mole rat promoter fragment ((-272)G-->CA) did not affect its activity. Oligonucleotides from the wild-type mouse and mole rat alphaB-crystallin promoter region under study formed distinct complexes with nuclear proteins from cultured cells. The mouse mutant sequence lost binding ability, whereas the mutated mole rat sequence gained the ability to form a complex similar in size to that of the wild-type mouse oligonucleotide. Our data support the idea that blind mole rats' alphaB-crystallin promoter activity was modified during the evolution of subterranean life and shows that tissue-specific promoter activity can be modulated by changing as few as two apparently neutral nucleotides in the mouse alphaB-crystallin enhancer region, implying the importance of the context of regulatory sequences for promoter activity.