Cis-regulatory Changes Research Articles

Copy Number Variation and Transposable Elements Feature in Recent, Ongoing Adaptation at the Cyp6g1 Locus

The increased transcription of the Cyp6g1 gene of Drosophila melanogaster, and consequent resistance to insecticides such as DDT, is a widely cited example of adaptation mediated by cis-regulatory change. A fragment of an Accord transposable element inserted upstream of the Cyp6g1 gene is causally associated with resistance and has spread to high frequencies in populations around the world since the 1940s. Here we report the existence of a natural allelic series at this locus of D. melanogaster, involving copy number variation of Cyp6g1, and two additional transposable element insertions (a P and an HMS-Beagle). We provide evidence that this genetic variation underpins phenotypic variation, as the more derived the allele, the greater the level of DDT resistance. Tracking the spatial and temporal patterns of allele frequency changes indicates that the multiple steps of the allelic series are adaptive. Further, a DDT association study shows that the most resistant allele, Cyp6g1-[BP], is greatly enriched in the top 5% of the phenotypic distribution and accounts for ∼16% of the underlying phenotypic variation in resistance to DDT. In contrast, copy number variation for another candidate resistance gene, Cyp12d1, is not associated with resistance. Thus the Cyp6g1 locus is a major contributor to DDT resistance in field populations, and evolution at this locus features multiple adaptive steps occurring in rapid succession.

The Importance of Being Cis: Evolution of Orthologous Fish and Mammalian Enhancer Activity

Conserved noncoding elements (CNEs) in vertebrate genomes often act as developmental enhancers, but a critical issue is how well orthologous CNE sequences retain the same activity in their respective species, a characteristic important for generalization of model organism studies. To quantify how well CNE enhancer activity has been preserved, we compared the anatomy-specific activities of 41 zebra fish CNEs in zebra fish embryos with the activities of orthologous human CNEs in mouse embryos. We found that 13/41 (∼30%) of the orthologous CNE pairs exhibit conserved positive activity in zebra fish and mouse. Conserved positive activity is only weakly associated with either sequence conservation or the absence of bases undergoing accelerated evolution. A stronger effect is that disparate activity is associated with transcription factor binding site divergence. To distinguish the contributions of cis- versus trans-regulatory changes, we analyzed 13 CNEs in a three-way experimental comparison: human CNE tested in zebra fish, human CNE tested in mouse, and orthologous zebra fish CNE tested in zebra fish. Both cis- and trans-changes affect a significant fraction of CNEs, although human and zebra fish sequences exhibit disparate activity in zebra fish (indicating cis regulatory changes) twice as often as human sequences show disparate activity when tested in mouse and zebra fish (indicating trans regulatory changes). In all four cases where the zebra fish and human CNE display a similar expression pattern in zebra fish, the human CNE also displays a similar expression pattern in mouse. This suggests that the endogenous enhancer activity of ∼30% of human CNEs can be determined from experiments in zebra fish alone, and to identify these CNEs, both the zebra fish and the human sequences should be tested.

Species‐specific coordinated gene expression and trans‐regulation of larval color pattern in three swallowtail butterflies

The diversity of butterfly larval color pattern has been attracted to people since Darwin's time; however, its molecular mechanisms still remain largely unknown. Larval body markings often differ completely between closely related species under natural selection. The final instar larvae of the swallowtail butterflies Papilio xuthus and Papilio polytes show a green camouflage pattern, whereas those of Papilio machaon show a warning color pattern, although P. xuthus and P. machaon are closely related species. To identify the genes that contribute to species divergence, we compared the expression pattern of eight pigment-associated genes between three Papilio species. The spatial expression pattern of melanin-related genes coincided with the species-specific cuticular markings. We newly found that the combination of bilin-binding protein and yellow-related gene (YRG) correlated perfectly with larval blue, yellow, and green coloration. To distinguish whether the interspecific differences in pigment-associated genes are caused by cis-regulatory changes or distribution differences in trans-regulatory proteins, we compared species-specific mRNA expression in an F1 hybrid specimen. Px-YRG and Pp-YRG showed a similar expression pattern, suggesting that the change in expression of YRG is caused mainly by changes in the distribution of trans-regulatory proteins. Our findings shed light on the gene regulatory networks for butterfly larval color pattern.

Identical Skin Toxins by Convergent Molecular Adaptation in Frogs

The Tree of Life is rife with adaptive convergences at all scales and biological levels of complexity. However, natural selection is not likely to result in the independent evolution of identical gene products. Here we report such a striking example of evolutionary convergence in the toxic skin secretions of two distantly related frog lineages. Caeruleins are important decapeptides in pharmacological and clinical research [1] and are commonly believed to represent a single evolutionary class of peptides [2-4]. Instead, our phylogenetic analyses combining transcriptome and genome data reveal that independently evolved precursor genes encode identical caeruleins in Xenopus and Litoria frogs. The former arose by duplication from the cholecystokinin (cck) gene, whereas the latter was derived from the gastrin gene. These hormone genes that are involved in many physiological processes diverged early in vertebrate evolution, after a segmental duplication during the Cambrian period. Besides implicating convergent mutations of the peptide-encoding sequence, recurrent caerulein origins entail parallel shifts of expression from the gut-brain axis to skin secretory glands. These results highlight extreme structural convergence in anciently diverged genes as an evolutionary mechanism through which recurrent adaptation is attained across large phylogenetic distances.

Correlated changes between regulatory cis elements and condition-specific expression in paralogous gene families

Gene duplication is integral to evolution, providing novel opportunities for organisms to diversify in function. One fundamental pathway of functional diversification among initially redundant gene copies, or paralogs, is via alterations in their expression patterns. Although the mechanisms underlying expression divergence are not completely understood, transcription factor binding sites and nucleosome occupancy are known to play a significant role in the process. Previous attempts to detect genomic variations mediating expression divergence in orthologs have had limited success for two primary reasons. First, it is inherently challenging to compare expressions among orthologs due to variable trans-acting effects and second, previous studies have quantified expression divergence in terms of an overall similarity of expression profiles across multiple samples, thereby obscuring condition-specific expression changes. Moreover, the inherently inter-correlated expressions among homologs present statistical challenges, not adequately addressed in many previous studies. Using rigorous statistical tests, here we characterize the relationship between cis element divergence and condition-specific expression divergence among paralogous genes in Saccharomyces cerevisiae. In particular, among all combinations of gene family and TFs analyzed, we found a significant correlation between TF binding and the condition-specific expression patterns in over 20% of the cases. In addition, incorporating nucleosome occupancy reveals several additional correlations. For instance, our results suggest that GAL4 binding plays a major role in the expression divergence of the genes in the sugar transporter family. Our work presents a novel means of investigating the cis regulatory changes potentially mediating expression divergence in paralogous gene families under specific conditions.

A Single Amino Acid Replacement in ETC2 Shapes Trichome Patterning in Natural Arabidopsis Populations

Our understanding of the evolution of organismal diversity is restricted by the current resolution of the genotype-phenotype map. In particular, the genetic basis of environmentally relevant phenotypic variation among natural populations remains poorly understood. Trichomes are single-cell outgrowths on the surface of plant leaves and other above-ground organs. Consistent with trichomes' suggested function to protect plants from predators and abiotic stressors [1-3], trichome density is strikingly variable among natural populations (e.g., [2, 4]). Despite substantial progress in the genetic dissection of trichome development [5], how trichome number is modulated in natural populations remains enigmatic. Here, we show that the ENHANCER OF TRY AND CPC 2 (ETC2) from the single-repeat R3 MYB family is the major locus determining trichome patterning in natural Arabidopsis populations. Our study identifies a single amino acid substitution in ETC2 (K19E) as the causal quantitative trait nucleotide (QTN). We suggest that this amino acid replacement might affect the stability of the ETC2 repressor, which results in a reduced trichome number. This is consistent with the view that morphology can evolve by coding changes that can subtly modulate protein activity as well as cis-regulatory changes that alter expression patterns.

Cis-regulatory Changes atFLOWERING LOCUS TMediate Natural Variation in Flowering Responses ofArabidopsis thaliana

Flowering time, a critical adaptive trait, is modulated by several environmental cues. These external signals converge on a small set of genes that in turn mediate the flowering response. Mutant analysis and subsequent molecular studies have revealed that one of these integrator genes, FLOWERING LOCUS T (FT), responds to photoperiod and temperature cues, two environmental parameters that greatly influence flowering time. As the central player in the transition to flowering, the protein coding sequence of FT and its function are highly conserved across species. Using QTL mapping with a new advanced intercross-recombinant inbred line (AI-RIL) population, we show that a QTL tightly linked to FT contributes to natural variation in the flowering response to the combined effects of photoperiod and ambient temperature. Using heterogeneous inbred families (HIF) and introgression lines, we fine map the QTL to a 6.7 kb fragment in the FT promoter. We confirm by quantitative complementation that FT has differential activity in the two parental strains. Further support for FT underlying the QTL comes from a new approach, quantitative knockdown with artificial microRNAs (amiRNAs). Consistent with the causal sequence polymorphism being in the promoter, we find that the QTL affects FT expression. Taken together, these results indicate that allelic variation at pathway integrator genes such as FT can underlie phenotypic variability and that this may be achieved through cis-regulatory changes.

Defending Evo-Devo: A Response to Hoekstra and Coyne

The study of evolutionary developmental biology (“evo-devo”) has recently experienced a dramatic surge in popularity among researchers and theorists concerned with evolution. However, some biologists and philosophers remain skeptical of the claims of evo-devo. This paper discusses and responds to the recent high profile criticisms of evo-devo presented by biologists Hopi E. Hoekstra and Jerry A. Coyne. I argue that their objections are unconvincing. Indeed, empirical research supports the main tenets of evo-devo, including the claim that morphological evolution is the result of cis-regulatory change and the distinction that evo-devo draws between morphological and physiological traits.

Candida glabrata PHO4 Is Necessary and Sufficient for Pho2-Independent Transcription of Phosphate Starvation Genes

Comparative genomic analyses of Candida glabrata and Saccharomyces cerevisiae suggest many signal transduction pathways are highly conserved. Focusing on the phosphate signal transduction (PHO) pathway of C. glabrata, we demonstrate that components of the pathway are conserved and confirm the role of CgPHO81, CgPHO80, CgPHO4, and CgMSN5 in the PHO pathway through deletion analysis. Unlike S. cerevisiae, C. glabrata shows little dependence on the transcription factor, Pho2, for induction of phosphate-regulated genes during phosphate limitation. We show that the CgPho4 protein is necessary and sufficient for Pho2-independent gene expression; CgPho4 is capable of driving expression of PHO promoters in S. cerevisiae in the absence of ScPHO2. On the basis of the sequences of PHO4 in the hemiascomycetes and complementation analysis, we suggest that Pho2 dependence is a trait only observed in species closely related to S. cerevisiae. Our data are consistent with trans-regulatory changes in the PHO pathway via the transcription factor Pho4 as opposed to cis-regulatory changes (the promoter).

Statistical tests for natural selection on regulatory regions based on the strength of transcription factor binding sites

BackgroundAlthough cis-regulatory changes play an important role in evolution, it remains difficult to establish the contribution of natural selection to regulatory differences between species. For protein coding regions, powerful tests of natural selection have been developed based on comparisons of synonymous and non-synonymous substitutions, and analogous tests for regulatory regions would be of great utility.ResultsHere, tests for natural selection on regulatory regions are proposed based on nucleotide substitutions that occur in characterized transcription factor binding sites (an important type functional element within regulatory regions). In the absence of selection, these substitutions will tend to reduce the strength of existing binding sites. On the other hand, purifying selection will act to preserve the binding sites in regulatory regions, while positive selection can act to create or destroy binding sites, as well as change their strength. Using standard models of binding site strength and molecular evolution in the absence of selection, this intuition can be used to develop statistical tests for natural selection. Application of these tests to two well-characterized regulatory regions in Drosophila provides evidence for purifying selection.ConclusionThis demonstrates that it is possible to develop tests for selection on regulatory regions based on the specific functional constrains on these sequences.

THE LOCI OF EVOLUTION: HOW PREDICTABLE IS GENETIC EVOLUTION?

Is genetic evolution predictable? Evolutionary developmental biologists have argued that, at least for morphological traits, the answer is a resounding yes. Most mutations causing morphological variation are expected to reside in the cis-regulatory, rather than the coding, regions of developmental genes. This “cis-regulatory hypothesis” has recently come under attack. In this review, we first describe and critique the arguments that have been proposed in support of the cis-regulatory hypothesis. We then test the empirical support for the cis-regulatory hypothesis with a comprehensive survey of mutations responsible for phenotypic evolution in multicellular organisms. Cis-regulatory mutations currently represent approximately 22% of 331 identified genetic changes although the number of cis-regulatory changes published annually is rapidly increasing. Above the species level, cis-regulatory mutations altering morphology are more common than coding changes. Also, above the species level cis-regulatory mutations predominate for genes not involved in terminal differentiation. These patterns imply that the simple question “Do coding or cis-regulatory mutations cause more phenotypic evolution?” hides more interesting phenomena. Evolution in different kinds of populations and over different durations may result in selection of different kinds of mutations. Predicting the genetic basis of evolution requires a comprehensive synthesis of molecular developmental biology and population genetics.

RESURRECTING THE ROLE OF TRANSCRIPTION FACTOR CHANGE IN DEVELOPMENTAL EVOLUTION

A long-standing question in evolutionary and developmental biology concerns the relative contribution of cis-regulatory and protein changes to developmental evolution. Central to this argument is which mutations generate evolutionarily relevant phenotypic variation? A review of the growing body of evolutionary and developmental literature supports the notion that many developmentally relevant differences occur in the cis-regulatory regions of protein-coding genes, generally to the exclusion of changes in the protein-coding region of genes. However, accumulating experimental evidence demonstrates that many of the arguments against a role for proteins in the evolution of gene regulation, and the developmental evolution in general, are no longer supported and there is an increasing number of cases in which transcription factor protein changes have been demonstrated in evolution. Here, we review the evidence that cis-regulatory evolution is an important driver of phenotypic evolution and provide examples of protein-mediated developmental evolution. Finally, we present an argument that the evolution of proteins may play a more substantial, but thus far underestimated, role in developmental evolution.

Evolution of metal hyperaccumulation required cis-regulatory changes and triplication of HMA4

Little is known about the types of mutations underlying the evolution of species-specific traits. The metal hyperaccumulator Arabidopsis halleri has the rare ability to colonize heavy-metal-polluted soils, and, as an extremophile sister species of Arabidopsis thaliana, it is a powerful model for research on adaptation. A. halleri naturally accumulates and tolerates leaf concentrations as high as 2.2% zinc and 0.28% cadmium in dry biomass. On the basis of transcriptomics studies, metal hyperaccumulation in A. halleri has been associated with more than 30 candidate genes that are expressed at higher levels in A. halleri than in A. thaliana. Some of these genes have been genetically mapped to broad chromosomal segments of between 4 and 24 cM co-segregating with Zn and Cd hypertolerance. However, the in planta loss-of-function approaches required to demonstrate the contribution of a given candidate gene to metal hyperaccumulation or hypertolerance have not been pursued to date. Using RNA interference to downregulate HMA4 (HEAVY METAL ATPASE 4) expression, we show here that Zn hyperaccumulation and full hypertolerance to Cd and Zn in A. halleri depend on the metal pump HMA4. Contrary to a postulated global trans regulatory factor governing high expression of numerous metal hyperaccumulation genes, we demonstrate that enhanced expression of HMA4 in A. halleri is attributable to a combination of modified cis-regulatory sequences and copy number expansion, in comparison to A. thaliana. Transfer of an A. halleri HMA4 gene to A. thaliana recapitulates Zn partitioning into xylem vessels and the constitutive transcriptional upregulation of Zn deficiency response genes characteristic of Zn hyperaccumulators. Our results demonstrate the importance of cis-regulatory mutations and gene copy number expansion in the evolution of a complex naturally selected extreme trait. The elucidation of a natural strategy for metal hyperaccumulation enables the rational design of technologies for the clean-up of metal-contaminated soils and for bio-fortification.

The Evolution of Gene Regulation Underlies a Morphological Difference between Two Drosophila Sister Species

Understanding the mechanisms underlying the morphological divergence of species is one of the central goals of evolutionary biology. Here, we analyze the genetic and molecular bases of the divergence of body pigmentation patterns between Drosophila yakuba and its sister species Drosophila santomea. We found that loss of pigmentation in D. santomea involved the selective loss of expression of the tan and yellow pigmentation genes. We demonstrate that tan gene expression was eliminated through the mutational inactivation of one specific tan cis-regulatory element (CRE) whereas the Tan protein sequence remained unchanged. Surprisingly, we identify three independent loss-of-function alleles of the tan CRE in the young D. santomea lineage. We submit that there is sufficient empirical evidence to support the general prediction that functional evolutionary changes at pleiotropic loci will most often involve mutations in their discrete, modular cis-regulatory elements.

Regulatory changes underlying expression differences within and between Drosophila species

Differences in gene expression are an important source of phenotypic variation, and can be caused by changes in cis and/or trans regulation. cis-regulatory variants alter allele-specific expression, whereas trans-regulatory variants influence expression of both alleles in a diploid cell. Because of this difference, we hypothesize that natural selection may favor one type of change over the other. Here, we investigate contributions of cis- and trans-regulatory changes to variable intra- and interspecific gene expression using four strains of Drosophila melanogaster, three strains of D. simulans and a total of 78 genes. We show that cis-regulatory changes account for a greater proportion of the expression differences observed between rather than within species. These data are inconsistent with a neutral model assuming equal probabilities of fixation for cis- and trans-regulatory polymorphisms, suggesting that natural selection influences the molecular mechanisms underlying divergent gene expression. Specifically, cis-regulatory changes seem to accumulate preferentially over time.

Analysis of mammalian gene batteries reveals both stable ancestral cores and highly dynamic regulatory sequences

Analysis of the evolutionary dynamics of target gene batteries controlled by 16 different transcription factors reveals stable ancestral cores and highly dynamic regulatory sequences

Cis-Regulatory Changes in Kit Ligand Expression and Parallel Evolution of Pigmentation in Sticklebacks and Humans

Dramatic pigmentation changes have evolved within most vertebrate groups, including fish and humans. Here we use genetic crosses in sticklebacks to investigate the parallel origin of pigmentation changes in natural populations. High-resolution mapping and expression experiments show that light gills and light ventrums map to a divergent regulatory allele of the Kit ligand (Kitlg) gene. The divergent allele reduces expression in gill and skin tissue and is shared by multiple derived freshwater populations with reduced pigmentation. In humans, Europeans and East Asians also share derived alleles at the KITLG locus. Strong signatures of selection map to regulatory regions surrounding the gene, and admixture mapping shows that the KITLG genomic region has a significant effect on human skin color. These experiments suggest that regulatory changes in Kitlg contribute to natural variation in vertebrate pigmentation, and that similar genetic mechanisms may underlie rapid evolutionary change in fish and humans.

An Embarrassment of Switches

What makes a human different from a chimpanzee or a mouse? Of course, we know the answer in broad outline. Mutations in the genome, sifted by natural selection, cause changes in appearance, physiology, and behavior—what geneticists call the phenotype. But we have only a vague picture of a more detailed answer. Precisely which mutations generate phenotypic evolution? It's not that we can't find the mutations. Today's DNA sequencing technology readily identifies all differences between two genomes. There are simply too many differences—tens of millions between human and chimp, for example (1). An unknown fraction of these mutations alter the phenotype. Nonetheless, the molecular effects of mutations provide a rough guide to their phenotypic effects. Some mutations change the amino acid sequence of proteins, thereby altering their functions, and some change so-called cis-regulatory regions, altering when and where proteins are produced. We know a lot about the first class, but much less about the second. Several recent papers, including one by Borneman et al. on page 815 of this issue (2), demonstrate a surprising abundance of cis-regulatory changes between closely related species.

Morphological evolution through multiple cis-regulatory mutations at a single gene

One central, and yet unsolved, question in evolutionary biology is the relationship between the genetic variants segregating within species and the causes of morphological differences between species. The classic neo-darwinian view postulates that species differences result from the accumulation of small-effect changes at multiple loci. However, many examples support the possible role of larger abrupt changes in the expression of developmental genes in morphological evolution. Although this evidence might be considered a challenge to a neo-darwinian micromutationist view of evolution, there are currently few examples of the actual genes causing morphological differences between species. Here we examine the genetic basis of a trichome pattern difference between Drosophila species, previously shown to result from the evolution of a single gene, shavenbaby (svb), probably through cis-regulatory changes. We first identified three distinct svb enhancers from D. melanogaster driving reporter gene expression in partly overlapping patterns that together recapitulate endogenous svb expression. All three homologous enhancers from D. sechellia drive expression in modified patterns, in a direction consistent with the evolved svb expression pattern. To test the influence of these enhancers on the actual phenotypic difference, we conducted interspecific genetic mapping at a resolution sufficient to recover multiple intragenic recombinants. This functional analysis revealed that independent genetic regions upstream of svb that overlap the three identified enhancers are collectively required to generate the D. sechellia trichome pattern. Our results demonstrate that the accumulation of multiple small-effect changes at a single locus underlies the evolution of a morphological difference between species. These data support the view that alleles of large effect that distinguish species may sometimes reflect the accumulation of multiple mutations of small effect at select genes.

Accelerated Evolution of Conserved Noncoding Sequences in Humans

Changes in gene regulation likely influenced the profound phenotypic divergence of humans from other mammals, but the extent of adaptive substitution in human regulatory sequences remains unknown. We identified 992 conserved noncoding sequences (CNSs) with a significant excess of human-specific substitutions. These accelerated elements were disproportionately found near genes involved in neuronal cell adhesion. To assess the uniqueness of human noncoding evolution, we examined CNSs accelerated in chimpanzee and mouse. Although we observed a similar enrichment near neuronal adhesion genes in chimpanzee, the accelerated CNSs themselves exhibited almost no overlap with those in human, suggesting independent evolution toward different neuronal phenotypes in each species. CNSs accelerated in mouse showed no bias toward neuronal cell adhesion. Our results indicate that widespread cis-regulatory changes in human evolution may have contributed to uniquely human features of brain development and function.