Repeated morphological evolution through cis-regulatory changes in a pleiotropic gene

Pigmentation and behavior: potential association through pleiotropic genes in Drosophila

The molecular basis of pigmentation variation within and among Drosophila species is largely attributed to genes in melanin biosynthesis pathway, which involves dopamine metabolism. Most of the genetic changes underlying pigmentation variations reported to date are changes at the expression levels of the structural genes in the pathway. Within D. melanogaster, changes in cis-regulatory regions of a gene, ebony, are responsible for the naturally occurring variation of the body pigmentation intensity. This gene is also known to be expressed in glia, and many visual and behavioral abnormalities of its mutants have been reported. This implies that the gene has pleiotropic functions in the nervous systems. In this review, current knowledge on pigmentation variation and melanin biosynthesis pathway are summarized, with some focus on pleiotropic features of ebony and other genes in the pathway. A potential association between pigmentation and behavior through such pleiotropic genes is discussed in light of cis-regulatory structure and pleiotropic mutations.

The making of a pest: the evolution of a fruit-penetrating ovipositor in Drosophila suzukii and related species

Evolutionary innovation can allow a species access to a new ecological niche, potentially reducing competition with closely related species. While the vast majority of Drosophila flies feed on rotting fruit and other decaying matter, and are harmless to human activity, Drosophila suzukii, which has a morphologically modified ovipositor, is capable of colonizing live fruit that is still in the process of ripening, causing massive agricultural damage. Here, we conducted the first comparative analysis of this species and its close relatives, analysing both ovipositor structure and fruit susceptibility. We found that the ovipositor of the species most closely related to D. suzukii, Drosophila subpulchrella, has a similar number of enlarged, evolutionarily derived bristles, but a notably different overall shape. Like D. suzukii, D. subpulchrella flies are capable of puncturing the skin of raspberries and cherries, but we found no evidence that they could penetrate the thicker skin of two varieties of grapes. More distantly related species, one of which has previously been mistaken for D. suzukii, have blunt ovipositors with small bristles. While they did not penetrate fruit skin in any of the assays, they readily colonized fruit interiors where the skin was broken. Our results suggest that considering evolutionary context may be beneficial to the management of invasive species.

A genome-wide scan for evidence of selection in a maize population under long-term artificial selection for ear number

A genome-wide scan to detect evidence of selection was conducted in the Golden Glow maize long-term selection population. The population had been subjected to selection for increased number of ears per plant for 30 generations, with an empirically estimated effective population size ranging from 384 to 667 individuals and an increase of more than threefold in the number of ears per plant. Allele frequencies at >1.2 million single-nucleotide polymorphism loci were estimated from pooled whole-genome resequencing data, and FST values across sliding windows were employed to assess divergence between the population preselection and the population postselection. Twenty-eight highly divergent regions were identified, with half of these regions providing gene-level resolution on potentially selected variants. Approximately 93% of the divergent regions do not demonstrate a significant decrease in heterozygosity, which suggests that they are not approaching fixation. Also, most regions display a pattern consistent with a soft-sweep model as opposed to a hard-sweep model, suggesting that selection mostly operated on standing genetic variation. For at least 25% of the regions, results suggest that selection operated on variants located outside of currently annotated coding regions. These results provide insights into the underlying genetic effects of long-term artificial selection and identification of putative genetic elements underlying number of ears per plant in maize.

Whole genome phylogeny for 21 Drosophila species using predicted 2b-RAD fragments

Type IIB restriction endonucleases are site-specific endonucleases that cut both strands of double-stranded DNA upstream and downstream of their recognition sequences. These restriction enzymes have recognition sequences that are generally interrupted and range from 5 to 7 bases long. They produce DNA fragments which are uniformly small, ranging from 21 to 33 base pairs in length (without cohesive ends). The fragments are generated from throughout the entire length of a genomic DNA providing an excellent fractional representation of the genome. In this study we simulated restriction enzyme digestions on 21 sequenced genomes of various Drosophila species using the predicted targets of 16 Type IIB restriction enzymes to effectively produce a large and arbitrary selection of loci from these genomes. The fragments were then used to compare organisms and to calculate the distance between genomes in pair-wise combination by counting the number of shared fragments between the two genomes. Phylogenetic trees were then generated for each enzyme using this distance measure and the consensus was calculated. The consensus tree obtained agrees well with the currently accepted tree for the Drosophila species. We conclude that multi-locus sub-genomic representation combined with next generation sequencing, especially for individuals and species without previous genome characterization, can accelerate studies of comparative genomics and the building of accurate phylogenetic trees.

Two genomic regions together cause dark abdominal pigmentation in Drosophila tenebrosa

Pigmentation is a rapidly evolving trait that is under both natural and sexual selection in many organisms. In the quinaria group of Drosophila, nearly all of the 30 species have an abdomen that is light in color with distinct markings; D. tenebrosa is the exception in that it has a completely melanic abdomen with no visible markings. In this study, we use a combination of quantitative genetic and candidate gene approaches to investigate the genetic basis of abdominal pigmentation in D. tenebrosa. We find that abdominal pigmentation is invariant across wild-caught lines of D. tenebrosa and is not sexually dimorphic. Quantitative genetic mapping utilizing crosses between D. tenebrosa and the light-colored D. suboccidentalis indicates that two genomic regions together underlie abdominal pigmentation, including the X-chromosome and an autosome (Muller Element C/E). Further support for their central importance in pigmentation is that experimental introgression of one phenotype into the other species, in either direction, results in introgression of these two genomic regions. Finally, the expression of the X-linked gene yellow in the pupae exactly foreshadows the adult melanization pattern in the abdomen of both species, suggesting that changes in the regulation of yellow are important for the phenotypic divergence of D. tenebrosa from the rest of the quinaria group. These results contribute to a body of work that demonstrates how changes in expression of highly conserved genes can cause substantial phenotypic differences even between closely related species.

Genomics and the Evolution of Phenotypic Traits

Evolutionary genetics has entered an unprecedented era of discovery, catalyzed in large part by the development of technologies that provide information about genome sequence and function. An important benefit is the ability to move beyond a handful of model organisms in lab settings to identify the genetic basis for evolutionarily interesting traits in many organisms in natural settings. Other benefits are the abilities to identify causal mutations and validate their phenotypic consequences more readily and in many more species. Genomic technologies have reinvigorated interest in some of the most fundamental and persistent questions in evolutionary genetics, revealed previously unsuspected evolutionary phenomena, and opened the door to a wide range of new questions.

Substrate vibrations during courtship in three Drosophila species

While a plethora of studies have focused on the role of visual, chemical and near-field airborne signals in courtship of Drosophila fruit flies, the existence of substrate-borne vibrational signals has been almost completely overlooked. Here we describe substrate vibrations generated during courtship in three species of the D. melanogaster group, from the allegedly mute species D. suzukii, its sister species D. biarmipes, and from D. melanogaster. In all species, we recorded several types of substrate vibrations which were generated by locomotion, abdominal vibrations and most likely through the activity of thoracic wing muscles. In D. melanogaster and D. suzukii, all substrate vibrations described in intact males were also recorded in males with amputated wings. Evidence suggests that vibrational signalling may be widespread among Drosophila species, and fruit flies may provide an ideal model to study various aspects of this widespread form of animal communication.

Evolution of transcriptional enhancers and animal diversity

Deciphering the genetic bases that drive animal diversity is one of the major challenges of modern biology. Although four decades ago it was proposed that animal evolution was mainly driven by changes in cis-regulatory DNA elements controlling gene expression rather than in protein-coding sequences, only now are powerful bioinformatics and experimental approaches available to accelerate studies into how the evolution of transcriptional enhancers contributes to novel forms and functions. In the introduction to this Theme Issue, we start by defining the general properties of transcriptional enhancers, such as modularity and the coexistence of tight sequence conservation with transcription factor-binding site shuffling as different mechanisms that maintain the enhancer grammar over evolutionary time. We discuss past and current methods used to identify cell-type-specific enhancers and provide examples of how enhancers originate de novo, change and are lost in particular lineages. We then focus in the central part of this Theme Issue on analysing examples of how the molecular evolution of enhancers may change form and function. Throughout this introduction, we present the main findings of the articles, reviews and perspectives contributed to this Theme Issue that together illustrate some of the great advances and current frontiers in the field.

Evolutionary change within a bipotential switch shaped the sperm/oocyte decision in hermaphroditic nematodes

A subset of transcription factors like Gli2 and Oct1 are bipotential--they can activate or repress the same target, in response to changing signals from upstream genes. Some previous studies implied that the sex-determination protein TRA-1 might also be bipotential; here we confirm this hypothesis by identifying a co-factor, and use it to explore how the structure of a bipotential switch changes during evolution. First, null mutants reveal that C. briggsae TRR-1 is required for spermatogenesis, RNA interference implies that it works as part of the Tip60 Histone Acetyl Transferase complex, and RT-PCR data show that it promotes the expression of Cbr-fog-3, a gene needed for spermatogenesis. Second, epistasis tests reveal that TRR-1 works through TRA-1, both to activate Cbr-fog-3 and to control the sperm/oocyte decision. Since previous studies showed that TRA-1 can repress fog-3 as well, these observations demonstrate that it is bipotential. Third, TRR-1 also regulates the development of the male tail. Since Cbr-tra-2 Cbr-trr-1 double mutants resemble Cbr-tra-1 null mutants, these two regulatory branches control all tra-1 activity. Fourth, striking differences in the relationship between these two branches of the switch have arisen during recent evolution. C. briggsae trr-1 null mutants prevent hermaphrodite spermatogenesis, but not Cbr-fem null mutants, which disrupt the other half of the switch. On the other hand, C. elegans fem null mutants prevent spermatogenesis, but not Cel-trr-1 mutants. However, synthetic interactions confirm that both halves of the switch exist in each species. Thus, the relationship between the two halves of a bipotential switch can shift rapidly during evolution, so that the same phenotype is produce by alternative, complementary mechanisms.

Evolution of eye development in the darkness of caves: adaptation, drift, or both?

Animals inhabiting the darkness of caves are generally blind and de-pigmented, regardless of the phylum they belong to. Survival in this environment is an enormous challenge, the most obvious being to find food and mates without the help of vision, and the loss of eyes in cave animals is often accompanied by an enhancement of other sensory apparatuses. Here we review the recent literature describing developmental biology and molecular evolution studies in order to discuss the evolutionary mechanisms underlying adaptation to life in the dark. We conclude that both genetic drift (neutral hypothesis) and direct and indirect selection (selective hypothesis) occurred together during the loss of eyes in cave animals. We also identify some future directions of research to better understand adaptation to total darkness, for which integrative analyses relying on evo-devo approaches associated with thorough ecological and population genomic studies should shed some light.

Linking genomics and ecology to investigate the complex evolution of an invasive Drosophila pest

Drosophilid fruit flies have provided science with striking cases of behavioral adaptation and genetic innovation. A recent example is the invasive pest Drosophila suzukii, which, unlike most other Drosophila, lays eggs and feeds on undamaged, ripening fruits. This not only poses a serious threat for fruit cultivation but also offers an interesting model to study evolution of behavioral innovation. We developed genome and transcriptome resources for D. suzukii. Coupling analyses of these data with field observations, we propose a hypothesis of the origin of its peculiar ecology. Using nuclear and mitochondrial phylogenetic analyses, we confirm its Asian origin and reveal a surprising sister relationship between the eugracilis and the melanogaster subgroups. Although the D. suzukii genome is comparable in size and repeat content to other Drosophila species, it has the lowest nucleotide substitution rate among the species analyzed in this study. This finding is compatible with the overwintering diapause of D. suzukii, which results in a reduced number of generations per year compared with its sister species. Genome-scale relaxed clock analyses support a late Miocene origin of D. suzukii, concomitant with paleogeological and climatic conditions that suggest an adaptation to temperate montane forests, a hypothesis confirmed by field trapping. We propose a causal link between the ecological adaptations of D. suzukii in its native habitat and its invasive success in Europe and North America.

Recurrent modification of a conserved cis-regulatory element underlies fruit fly pigmentation diversity

The development of morphological traits occurs through the collective action of networks of genes connected at the level of gene expression. As any node in a network may be a target of evolutionary change, the recurrent targeting of the same node would indicate that the path of evolution is biased for the relevant trait and network. Although examples of parallel evolution have implicated recurrent modification of the same gene and cis-regulatory element (CRE), little is known about the mutational and molecular paths of parallel CRE evolution. In Drosophila melanogaster fruit flies, the Bric-à-brac (Bab) transcription factors control the development of a suite of sexually dimorphic traits on the posterior abdomen. Female-specific Bab expression is regulated by the dimorphic element, a CRE that possesses direct inputs from body plan (ABD-B) and sex-determination (DSX) transcription factors. Here, we find that the recurrent evolutionary modification of this CRE underlies both intraspecific and interspecific variation in female pigmentation in the melanogaster species group. By reconstructing the sequence and regulatory activity of the ancestral Drosophila melanogaster dimorphic element, we demonstrate that a handful of mutations were sufficient to create independent CRE alleles with differing activities. Moreover, intraspecific and interspecific dimorphic element evolution proceeded with little to no alterations to the known body plan and sex-determination regulatory linkages. Collectively, our findings represent an example where the paths of evolution appear biased to a specific CRE, and drastic changes in function were accompanied by deep conservation of key regulatory linkages.

Critiquing blind dating: the dangers of over-confident date estimates in comparative genomics

Phylogenomic advances provide more rigorous estimates for the timing of evolutionary divergences than previously available (e.g., Bayesian relaxed-clock estimates with soft fossil constraints). However, because many family-level clades and higher, as well as model species within those clades, have not been included in phylogenomic studies, the literature presents temporal estimates likely harboring substantial errors. Blindly using such dates can substantially retard scientific advancement. We suggest a way forward by conducting analyses that minimize prior assumptions and use large datasets, and demonstrate how using such a phylogenomic approach can lead to significantly more parsimonious conclusions without a good fossil record. We suggest that such an approach calls for research into the biological causes of conflict between molecular and fossil signatures.

Divergent genetic mechanisms underlie reversals to radial floral symmetry from diverse zygomorphic flowered ancestors

Malpighiaceae possess flowers with a unique bilateral symmetry (zygomorphy), which is a hypothesized adaptation associated with specialization on neotropical oil bee pollinators. Gene expression of two representatives of the CYC2 lineage of floral symmetry TCP genes, CYC2A and CYC2B, demarcate the adaxial (dorsal) region of the flower in the characteristic zygomorphic flowers of most Malpighiaceae. Several clades within the family, however, have independently lost their specialized oil bee pollinators and reverted to radial flowers (actinomorphy) like their ancestors. Here, we investigate CYC2 expression associated with four independent reversals to actinomorphy. We demonstrate that these reversals are always associated with alteration of the highly conserved CYC2 expression pattern observed in most New World (NW) Malpighiaceae. In NW Lasiocarpus and Old World (OW) Microsteria, the expression of CYC2-like genes has expanded to include the ventral region of the corolla. Thus, the pattern of gene expression in these species has become radialized, which is comparable to what has been reported in the radial flowered legume clade Cadia. In striking contrast, in NW Psychopterys and OW Sphedamnocarpus, CYC2-like expression is entirely absent or at barely detectable levels. This is more similar to the pattern of CYC2 expression observed in radial flowered Arabidopsis. These results collectively indicate that, regardless of geographic distribution, reversals to similar floral phenotypes in this large tropical angiosperm clade have evolved via different genetic changes from an otherwise highly conserved developmental program.

Employing BAC-reporter constructs in the sea anemone Nematostella vectensis

Changes in the expression and function of genes drive evolutionary change. Comparing how genes are regulated in different species is therefore becoming an important part of evo-devo studies. A key tool for investigating the regulation of genes is represented by bacterial artificial chromosomes (BAC)-reporter constructs. BACs are large insert libraries, often >100 kb, which thus capture the genomic sequences surrounding a gene of interest, including all, or nearly all, of the elements underpinning regulation. Recombinant BACs, containing a reporter gene in place of the endogenous coding sequence of genes, can be utilized to drive the expression of reporter genes under the regulatory control of the gene of interest while still embedded within its genomic context. Systematic deletions within the BAC-reporter construct can be used to identify the minimal reporter in an unbiased way, avoiding the risk of overlooking regulatory elements that may be many kilobases away from the transcription start-site. Nematostella vectensis (Edwardsiidae, Anthozoa, Cnidaria) has become an important model in regenerative biology, ecology, and especially in studies of evo-devo and gene-regulatory networks due to its interesting phylogenetic position and amenability to molecular techniques. The increasing interest in this rising model system also led to a demand for methods that can be used to study the regulation of genes in Nematostella. Here, we present our progress in employing BAC-reporter constructs to visualize gene-expression in Nematostella. Using a new Nematostella-specific recombination cassette, we made nine different BAC-reporter constructs. Although five BAC recombinants gave variable effects, three constructs, namely Nv-bra:eGFP::L10 BAC, Nv-dpp:eGFP::L10 BAC, and Nv-grm:eGFP::L10 BAC delivered promising results. We show that these three constructs express the reporter gene eGFP in 10.4-17.2% of all analyzed larvae, out of which 26.2-41.9% express GFP in a mosaic fashion within the expected domain. In addition to the expression within the known domains, we also observed cases of misexpression of eGFP and examples that could represent actual expression outside the described domain. Furthermore, we deep-sequenced and assembled five different BACs containing Nv-chordin, Nv-foxa, Nv-dpp, Nv-wnta, and Nv-wnt1, to improve assembly around these genes. The use of BAC-reporter constructs will foster cis-regulatory analyses in Nematostella and thus help to improve our understanding of the regulatory network in this cnidarian system. Ultimately, this will advance the comparison of gene-regulation across species and lead to a much better understanding of evolutionary changes and novelties.

Stepwise evolution of essential centromere function in a Drosophila neogene

Evolutionarily young genes that serve essential functions represent a paradox; they must perform a function that either was not required until after their birth or was redundant with another gene. How young genes rapidly acquire essential function is largely unknown. We traced the evolutionary steps by which the Drosophila gene Umbrea acquired an essential role in chromosome segregation in D. melanogaster since the gene's origin less than 15 million years ago. Umbrea neofunctionalization occurred via loss of an ancestral heterochromatin-localizing domain, followed by alterations that rewired its protein interaction network and led to species-specific centromere localization. Our evolutionary cell biology approach provides temporal and mechanistic detail about how young genes gain essential function. Such innovations may constantly alter the repertoire of centromeric proteins in eukaryotes.

From remote enhancers to gene regulation: charting the genome's regulatory landscapes

Vertebrate genes are characterized by the presence of cis-regulatory elements located at great distances from the genes they control. Alterations of these elements have been implicated in human diseases and evolution, yet little is known about how these elements interact with their surrounding sequences. A recent survey of the mouse genome with a regulatory sensor showed that the regulatory activities of these elements are not organized in a gene-centric manner, but instead are broadly distributed along chromosomes, forming large regulatory landscapes with distinct tissue-specific activities. A large genome-wide collection of expression data from this regulatory sensor revealed some basic principles of this complex genome regulatory architecture, including a substantial interplay between enhancers and other types of activities to modulate gene expression. We discuss the implications of these findings for our understanding of non-coding transcription, and of the possible consequences of structural genomic variations in disease and evolution.

Population genomics and transcriptional consequences of regulatory motif variation in globally diverse Saccharomyces cerevisiae strains

Noncoding genetic variation is known to significantly influence gene expression levels in a growing number of specific cases; however, the patterns of genome-wide noncoding variation present within populations, the evolutionary forces acting on noncoding variants, and the relative effects of regulatory polymorphisms on transcript abundance are not well characterized. Here, we address these questions by analyzing patterns of regulatory variation in motifs for 177 DNA binding proteins in 37 strains of Saccharomyces cerevisiae. Between S. cerevisiae strains, we found considerable polymorphism in regulatory motifs across strains (mean π = 0.005) as well as diversity in regulatory motifs (mean 0.91 motifs differences per regulatory region). Population genetics analyses reveal that motifs are under purifying selection, and there is considerable heterogeneity in the magnitude of selection across different motifs. Finally, we obtained RNA-Seq data in 22 strains and identified 49 polymorphic DNA sequence motifs in 30 distinct genes that are significantly associated with transcriptional differences between strains. In 22 of these genes, there was a single polymorphic motif associated with expression in the upstream region. Our results provide comprehensive insights into the evolutionary trajectory of regulatory variation in yeast and the characteristics of a compendium of regulatory alleles.

Emergence and diversification of fly pigmentation through evolution of a gene regulatory module

The typical pattern of morphological evolution associated with the radiation of a group of related species is the emergence of a novel trait and its subsequent diversification. Yet the genetic mechanisms associated with these two evolutionary steps are poorly characterized. Here, we show that a spot of dark pigment on fly wings emerged from the assembly of a novel gene regulatory module in which a set of pigmentation genes evolved to respond to a common transcriptional regulator determining their spatial distribution. The primitive wing spot pattern subsequently diversified through changes in the expression pattern of this regulator. These results suggest that the genetic changes underlying the emergence and diversification of wing pigmentation patterns are partitioned within genetic networks.

Highly conserved elements discovered in vertebrates are present in non-syntenic loci of tunicates, act as enhancers and can be transcribed during development

Co-option of cis-regulatory modules has been suggested as a mechanism for the evolution of expression sites during development. However, the extent and mechanisms involved in mobilization of cis-regulatory modules remains elusive. To trace the history of non-coding elements, which may represent candidate ancestral cis-regulatory modules affirmed during chordate evolution, we have searched for conserved elements in tunicate and vertebrate (Olfactores) genomes. We identified, for the first time, 183 non-coding sequences that are highly conserved between the two groups. Our results show that all but one element are conserved in non-syntenic regions between vertebrate and tunicate genomes, while being syntenic among vertebrates. Nevertheless, in all the groups, they are significantly associated with transcription factors showing specific functions fundamental to animal development, such as multicellular organism development and sequence-specific DNA binding. The majority of these regions map onto ultraconserved elements and we demonstrate that they can act as functional enhancers within the organism of origin, as well as in cross-transgenesis experiments, and that they are transcribed in extant species of Olfactores. We refer to the elements as 'Olfactores conserved non-coding elements'.

Proteomic profiling of cereal aphid saliva reveals both ubiquitous and adaptive secreted proteins

The secreted salivary proteins from two cereal aphid species, Sitobion avenae and Metopolophium dirhodum, were collected from artificial diets and analysed by tandem mass spectrometry. Protein identification was performed by searching MS data against the official protein set from the current pea aphid (Acyrthosiphon pisum) genome assembly and revealed 12 and 7 proteins in the saliva of S. avenae and M. dirhodum, respectively. When combined with a comparable dataset from A. pisum, only three individual proteins were common to all the aphid species; two paralogues of the GMC oxidoreductase family (glucose dehydrogenase; GLD) and ACYPI009881, an aphid specific protein previously identified as a putative component of the salivary sheath. Antibodies were designed from translated protein sequences obtained from partial cDNA sequences for ACYPI009881 and both saliva associated GLDs. The antibodies detected all parent proteins in secreted saliva from the three aphid species, but could only detect ACYPI009881, and not saliva associated GLDs, in protein extractions from the salivary glands. This result was confirmed by immunohistochemistry using whole and sectioned salivary glands, and in addition, localised ACYPI009881 to specific cell types within the principal salivary gland. The implications of these findings for the origin of salivary components and the putative role of the proteins identified are discussed in the context of our limited understanding of the functional relationship between aphid saliva and the plants they feed on. The mass spectrometry data have been deposited to the ProteomeXchange and can be accessed under the identifier PXD000113.

Drosophila sex combs as a model of evolutionary innovations

The diversity of animal and plant forms is shaped by nested evolutionary innovations. Understanding the genetic and molecular changes responsible for these innovations is therefore one of the key goals of evolutionary biology. From the genetic point of view, the origin of novel traits implies the origin of new regulatory pathways to control their development. To understand how these new pathways are assembled in the course of evolution, we need model systems that combine relatively recent innovations with a powerful set of genetic and molecular tools. One such model is provided by the Drosophila sex comb-a male-specific morphological structure that evolved in a relatively small lineage related to the model species D. melanogaster. Our extensive knowledge of sex comb development in D. melanogaster provides the basis for investigating the genetic changes responsible for sex comb origin and diversification. At the same time, sex combs can change on microevolutionary timescales and differ spectacularly among closely related species, providing opportunities for direct genetic analysis and for integrating developmental and population-genetic approaches. Sex comb evolution is associated with the origin of novel interactions between Hox and sex determination genes. Activity of the sex determination pathway was brought under the control of the Hox code to become segment-specific, while Hox gene expression became sexually dimorphic. At the same time, both Hox and sex determination genes were integrated into the intrasegmental spatial patterning network, and acquired new joint downstream targets. Phylogenetic analysis shows that similar sex comb morphologies evolved independently in different lineages. Convergent evolution at the phenotypic level reflects convergent changes in the expression of Hox and sex determination genes, involving both independent gains and losses of regulatory interactions. However, the downstream cell-differentiation programs have diverged between species, and in some lineages, similar adult morphologies are produced by different morphogenetic mechanisms. These features make the sex comb an excellent model for examining not only the genetic changes responsible for its evolution, but also the cellular processes that translate DNA sequence changes into morphological diversity. The origin and diversification of sex combs provides insights into the roles of modularity, cooption, and regulatory changes in evolutionary innovations, and can serve as a model for understanding the origin of the more drastic novelties that define higher order taxa.

RNA-seq analysis of allele-specific expression, hybrid effects, and regulatory divergence in hybrids compared with their parents from natural populations

Hybridization is a prominent process among natural plant populations that can result in phenotypic novelty, heterosis, and changes in gene expression. The effects of intraspecific hybridization on F1 hybrid gene expression were investigated using parents from divergent, natural populations of Cirsium arvense, an invasive Compositae weed. Using an RNA-seq approach, the expression of 68,746 unigenes was quantified in parents and hybrids. The expression levels of 51% of transcripts differed between parents, a majority of which had less than 1.25× fold-changes. More unigenes had higher expression in the invasive parent (P1) than the noninvasive parent (P2). Of those that were divergently expressed between parents, 10% showed additive and 81% showed nonadditive (transgressive or dominant) modes of gene action in the hybrids. A majority of the dominant cases had P2-like expression patterns in the hybrids. Comparisons of allele-specific expression also enabled a survey of cis- and trans-regulatory effects. Cis- and trans-regulatory divergence was found at 70% and 68% of 62,281 informative single-nucleotide polymorphism sites, respectively. Of the 17% of sites exhibiting both cis- and trans-effects, a majority (70%) had antagonistic regulatory interactions (cis x trans); trans-divergence tended to drive higher expression of the P1 allele, whereas cis-divergence tended to increase P2 transcript abundance. Trans-effects correlated more highly than cis with parental expression divergence and accounted for a greater proportion of the regulatory divergence at sites with additive compared with nonadditive inheritance patterns. This study explores the nature of, and types of mechanisms underlying, expression changes that occur in upon intraspecific hybridization in natural populations.

Deciphering the evolutionary history and developmental mechanisms of a complex sexual ornament: the abdominal appendages of Sepsidae (Diptera)

Male abdomen appendages are a novel trait found within Sepsidae (Diptera). Here we demonstrate that they are likely to have evolved once, were lost three times, and then secondarily gained in one lineage. The developmental basis of these appendages was investigated by counting the number of histoblast cells in each abdominal segment in four species: two that represented the initial instance of appendage evolution, one that has secondarily gained appendages, and one species that did not have appendages. Males of all species with appendages have elevated cell counts for the fourth segment, which gives rise to the appendages. In Perochaeta dikowi, which reacquired the trait, the females also have elevated cell count on the fourth segment despite the fact that females do not develop appendages. The species without appendages has similar cell counts in all segments regardless of sex. These results suggest that the basis for appendage development is shared in males across all species, but the sexual dimorphism is regulated differently in P. dikowi.

Conserved regulatory architecture underlies parallel genetic changes and convergent phenotypic evolution

Similar morphological, physiological, and behavioral features have evolved independently in different species, a pattern known as convergence. It is known that morphological convergence can occur through changes in orthologous genes. In some cases of convergence, cis-regulatory changes generate parallel modifications in the expression patterns of orthologous genes. Our understanding of how changes in cis-regulatory regions contribute to convergence is hampered, usually, by a limited understanding of the global cis-regulatory structure of the evolving genes. Here we examine the genetic causes of a case of precise phenotypic convergence between Drosophila sechellia and Drosophila ezoana, species that diverged ~40 Mya. Previous studies revealed that changes in multiple transcriptional enhancers of shavenbaby (svb, a transcript of the ovo locus) caused phenotypic evolution in the D. sechellia lineage. It has also been shown that the convergent phenotype of D. ezoana was likely caused by cis-regulatory evolution of svb. Here we show that the large-scale cis-regulatory architecture of svb is conserved between these Drosophila species. Furthermore, we show that the D. ezoana orthologs of the evolved D. sechellia enhancers have also evolved expression patterns that correlate precisely with the changes in the phenotype. Our results suggest that phenotypic convergence resulted from multiple noncoding changes that occurred in parallel in the D. sechellia and D. ezoana lineages.

Unraveling the thread of nature's tapestry: the genetics of diversity and convergence in animal pigmentation

Animals display incredibly diverse color patterns yet little is known about the underlying genetic basis of these phenotypes. However, emerging results are reshaping our view of how the process of phenotypic evolution occurs. Here, we outline recent research from three particularly active areas of investigation: melanin pigmentation in Drosophila, wing patterning in butterflies, and pigment variation in lizards. For each system, we highlight (i) the function and evolution of color variation, (ii) various approaches that have been used to explore the genetic basis of pigment variation, and (iii) conclusions regarding the genetic basis of convergent evolution which have emerged from comparative analyses. Results from these studies indicate that natural variation in pigmentation is a particularly powerful tool to examine the molecular basis of evolution, especially with regard to convergent or parallel evolution. Comparison of these systems also reveals that the molecular basis of convergent evolution is heterogeneous, sometimes involving conserved mechanisms and sometimes not. In the near future, additional work in other emerging systems will substantially expand the scope of available comparisons.

Hundreds of conserved non-coding genomic regions are independently lost in mammals

Conserved non-protein-coding DNA elements (CNEs) often encode cis-regulatory elements and are rarely lost during evolution. However, CNE losses that do occur can be associated with phenotypic changes, exemplified by pelvic spine loss in sticklebacks. Using a computational strategy to detect complete loss of CNEs in mammalian genomes while strictly controlling for artifacts, we find >600 CNEs that are independently lost in at least two mammalian lineages, including a spinal cord enhancer near GDF11. We observed several genomic regions where multiple independent CNE loss events happened; the most extreme is the DIAPH2 locus. We show that CNE losses often involve deletions and that CNE loss frequencies are non-uniform. Similar to less pleiotropic enhancers, we find that independently lost CNEs are shorter, slightly less constrained and evolutionarily younger than CNEs without detected losses. This suggests that independently lost CNEs are less pleiotropic and that pleiotropic constraints contribute to non-uniform CNE loss frequencies. We also detected 35 CNEs that are independently lost in the human lineage and in other mammals. Our study uncovers an interesting aspect of the evolution of functional DNA in mammalian genomes. Experiments are necessary to test if these independently lost CNEs are associated with parallel phenotype changes in mammals.

Parallel molecular evolution in an herbivore community

Numerous insects have independently evolved the ability to feed on plants that produce toxic secondary compounds called cardenolides and can sequester these compounds for use in their defense. We surveyed the protein target for cardenolides, the alpha subunit of the sodium pump, Na(+),K(+)-ATPase (ATPα), in 14 species that feed on cardenolide-producing plants and 15 outgroups spanning three insect orders. Despite the large number of potential targets for modulating cardenolide sensitivity, amino acid substitutions associated with host-plant specialization are highly clustered, with many parallel substitutions. Additionally, we document four independent duplications of ATPα with convergent tissue-specific expression patterns. We find that unique substitutions are disproportionately associated with recent duplications relative to parallel substitutions. Together, these findings support the hypothesis that adaptation tends to take evolutionary paths that minimize negative pleiotropy.

Convergent evolution of a reproductive trait through distinct developmental mechanisms in Drosophila

Convergent morphologies often arise due to similar selective pressures in independent lineages. It is poorly understood whether the same or different developmental genetic mechanisms underlie such convergence. Here we show that independent evolution of a reproductive trait, ovariole number, has resulted from changes in distinct developmental mechanisms, each of which may have a different underlying genetic basis in Drosophila. Ovariole number in Drosophila is species-specific, highly variable, and largely under genetic control. Convergent changes in Drosophila ovariole number have evolved independently within and between species. We previously showed that the number of a specific ovarian cell type, terminal filament (TF) cells, determines ovariole number. Here we examine TF cell development in different Drosophila lineages that independently evolved a significantly lower ovariole number than the D. melanogaster Oregon R strain. We show that in these Drosophila lineages, reduction in ovariole number occurs primarily through variations in one of two different developmental mechanisms: (1) reduced number of somatic gonad precursors (SGP cells) specified during embryogenesis; or (2) alterations of somatic gonad cell morphogenesis and differentiation in larval life. Mutations in the D. melanogaster Insulin Receptor (InR) alter SGP cell number but not ovarian morphogenesis, while targeted loss of function of bric-à-brac 2 (bab2) affects morphogenesis without changing SGP cell number. Thus, evolution can produce similar ovariole numbers through distinct developmental mechanisms, likely controlled by different genetic mechanisms.

Cellular hyperproliferation and cancer as evolutionary variables

Technological advances in biology have begun to dramatically change the way we think about evolution, development, health and disease. The ability to sequence the genomes of many individuals within a population, and across multiple species, has opened the door to the possibility of answering some long-standing and perplexing questions about our own genetic heritage. One such question revolves around the nature of cellular hyperproliferation. This cellular behavior is used to effect wound healing in most animals, as well as, in some animals, the regeneration of lost body parts. Yet at the same time, cellular hyperproliferation is the fundamental pathological condition responsible for cancers in humans. Here, I will discuss why microevolution, macroevolution and developmental biology all have to be taken into consideration when interpreting studies of both normal and malignant hyperproliferation. I will also illustrate how a synthesis of evolutionary sciences and developmental biology through the study of diverse model organisms can inform our understanding of both health and disease.

The population genetics of cooperative gene regulation

Background

Changes in gene regulatory networks drive the evolution of phenotypic diversity both within and between species. Rewiring of transcriptional networks is achieved either by changes to transcription factor binding sites or by changes to the physical interactions among transcription factor proteins. It has been suggested that the evolution of cooperative binding among factors can facilitate the adaptive rewiring of a regulatory network.

Results

We use a population-genetic model to explore when cooperative binding of transcription factors is favored by evolution, and what effects cooperativity then has on the adaptive re-writing of regulatory networks. We consider a pair of transcription factors that regulate multiple targets and overlap in the sets of target genes they regulate. We show that, under stabilising selection, cooperative binding between the transcription factors is favoured provided the amount of overlap between their target genes exceeds a threshold. The value of this threshold depends on several population-genetic factors: strength of selection on binding sites, cost of pleiotropy associated with protein-protein interactions, rates of mutation and population size. Once it is established, we find that cooperative binding of transcription factors significantly accelerates the adaptive rewiring of transcriptional networks under positive selection. We compare our qualitative predictions to systematic data on Saccharomyces cerevisiae transcription factors, their binding sites, and their protein-protein interactions.

Conclusions

Our study reveals a rich set of evolutionary dynamics driven by a tradeoff between the beneficial effects of cooperative binding at targets shared by a pair of factors, and the detrimental effects of cooperative binding for non-shared targets. We find that cooperative regulation will evolve when transcription factors share a sufficient proportion of their target genes. These findings help to explain empirical pattens in datasets of transcription factors in Saccharomyces cerevisiae and, they suggest that changes to physical interactions between transcription factors can play a critical role in the evolution of gene regulatory networks.

Evolution of multiple additive loci caused divergence between Drosophila yakuba and D. santomea in wing rowing during male courtship

In Drosophila, male flies perform innate, stereotyped courtship behavior. This innate behavior evolves rapidly between fly species, and is likely to have contributed to reproductive isolation and species divergence. We currently understand little about the neurobiological and genetic mechanisms that contributed to the evolution of courtship behavior. Here we describe a novel behavioral difference between the two closely related species D. yakuba and D. santomea: the frequency of wing rowing during courtship. During courtship, D. santomea males repeatedly rotate their wing blades to face forward and then back (rowing), while D. yakuba males rarely row their wings. We found little intraspecific variation in the frequency of wing rowing for both species. We exploited multiplexed shotgun genotyping (MSG) to genotype two backcross populations with a single lane of Illumina sequencing. We performed quantitative trait locus (QTL) mapping using the ancestry information estimated by MSG and found that the species difference in wing rowing mapped to four or five genetically separable regions. We found no evidence that these loci display epistasis. The identified loci all act in the same direction and can account for most of the species difference.

Why transcription factor binding sites are ten nucleotides long

Gene expression is controlled primarily by transcription factors, whose DNA binding sites are typically 10 nt long. We develop a population-genetic model to understand how the length and information content of such binding sites evolve. Our analysis is based on an inherent trade-off between specificity, which is greater in long binding sites, and robustness to mutation, which is greater in short binding sites. The evolutionary stable distribution of binding site lengths predicted by the model agrees with the empirical distribution (5-31 nt, with mean 9.9 nt for eukaryotes), and it is remarkably robust to variation in the underlying parameters of population size, mutation rate, number of transcription factor targets, and strength of selection for proper binding and selection against improper binding. In a systematic data set of eukaryotic and prokaryotic transcription factors we also uncover strong relationships between the length of a binding site and its information content per nucleotide, as well as between the number of targets a transcription factor regulates and the information content in its binding sites. Our analysis explains these features as well as the remarkable conservation of binding site characteristics across diverse taxa.

Specialized motor-driven dusp1 expression in the song systems of multiple lineages of vocal learning birds

Mechanisms for the evolution of convergent behavioral traits are largely unknown. Vocal learning is one such trait that evolved multiple times and is necessary in humans for the acquisition of spoken language. Among birds, vocal learning is evolved in songbirds, parrots, and hummingbirds. Each time similar forebrain song nuclei specialized for vocal learning and production have evolved. This finding led to the hypothesis that the behavioral and neuroanatomical convergences for vocal learning could be associated with molecular convergence. We previously found that the neural activity-induced gene dual specificity phosphatase 1 (dusp1) was up-regulated in non-vocal circuits, specifically in sensory-input neurons of the thalamus and telencephalon; however, dusp1 was not up-regulated in higher order sensory neurons or motor circuits. Here we show that song motor nuclei are an exception to this pattern. The song nuclei of species from all known vocal learning avian lineages showed motor-driven up-regulation of dusp1 expression induced by singing. There was no detectable motor-driven dusp1 expression throughout the rest of the forebrain after non-vocal motor performance. This pattern contrasts with expression of the commonly studied activity-induced gene egr1, which shows motor-driven expression in song nuclei induced by singing, but also motor-driven expression in adjacent brain regions after non-vocal motor behaviors. In the vocal non-learning avian species, we found no detectable vocalizing-driven dusp1 expression in the forebrain. These findings suggest that independent evolutions of neural systems for vocal learning were accompanied by selection for specialized motor-driven expression of the dusp1 gene in those circuits. This specialized expression of dusp1 could potentially lead to differential regulation of dusp1-modulated molecular cascades in vocal learning circuits.

The evolution of root hairs and rhizoids

BACKGROUND

Almost all land plants develop tip-growing filamentous cells at the interface between the plant and substrate (the soil). Root hairs form on the surface of roots of sporophytes (the multicellular diploid phase of the life cycle) in vascular plants. Rhizoids develop on the free-living gametophytes of vascular and non-vascular plants and on both gametophytes and sporophytes of the extinct rhyniophytes. Extant lycophytes (clubmosses and quillworts) and monilophytes (ferns and horsetails) develop both free-living gametophytes and free-living sporophytes. These gametophytes and sporophytes grow in close contact with the soil and develop rhizoids and root hairs, respectively.

SCOPE

Here we review the development and function of rhizoids and root hairs in extant groups of land plants. Root hairs are important for the uptake of nutrients with limited mobility in the soil such as phosphate. Rhizoids have a variety of functions including water transport and adhesion to surfaces in some mosses and liverworts.

CONCLUSIONS

A similar gene regulatory network controls the development of rhizoids in moss gametophytes and root hairs on the roots of vascular plant sporophytes. It is likely that this gene regulatory network first operated in the gametophyte of the earliest land plants. We propose that later it functioned in sporophytes as the diploid phase evolved a free-living habit and developed an interface with the soil. This transference of gene function from gametophyte to sporophyte could provide a mechanism that, at least in part, explains the increase in morphological diversity of sporophytes that occurred during the radiation of land plants in the Devonian Period.

Phylogenomic analysis reveals dynamic evolutionary history of the Drosophila heterochromatin protein 1 (HP1) gene family

Heterochromatin is the gene-poor, satellite-rich eukaryotic genome compartment that supports many essential cellular processes. The functional diversity of proteins that bind and often epigenetically define heterochromatic DNA sequence reflects the diverse functions supported by this enigmatic genome compartment. Moreover, heterogeneous signatures of selection at chromosomal proteins often mirror the heterogeneity of evolutionary forces that act on heterochromatic DNA. To identify new such surrogates for dissecting heterochromatin function and evolution, we conducted a comprehensive phylogenomic analysis of the Heterochromatin Protein 1 gene family across 40 million years of Drosophila evolution. Our study expands this gene family from 5 genes to at least 26 genes, including several uncharacterized genes in Drosophila melanogaster. The 21 newly defined HP1s introduce unprecedented structural diversity, lineage-restriction, and germline-biased expression patterns into the HP1 family. We find little evidence of positive selection at these HP1 genes in both population genetic and molecular evolution analyses. Instead, we find that dynamic evolution occurs via prolific gene gains and losses. Despite this dynamic gene turnover, the number of HP1 genes is relatively constant across species. We propose that karyotype evolution drives at least some HP1 gene turnover. For example, the loss of the male germline-restricted HP1E in the obscura group coincides with one episode of dramatic karyotypic evolution, including the gain of a neo-Y in this lineage. This expanded compendium of ovary- and testis-restricted HP1 genes revealed by our study, together with correlated gain/loss dynamics and chromosome fission/fusion events, will guide functional analyses of novel roles supported by germline chromatin.

Our genome is comprised of two compartments. The euchromatin harbors abundant genes and regulatory information, while heterochromatin harbors few genes and abundant repetitive DNA. These characteristic features of heterochromatin challenge traditional methods of sequence assembly and molecular dissection. The analysis, instead, of proteins that localize to and often functionally define heterochromatic sequence has illuminated numerous heterochromatin-dependent, essential cellular processes, including chromosome segregation, telomere stability, and gene regulation. With the aim of increasing our sample of heterochromatin-localizing proteins, we performed a comprehensive search for new members of Heterochromatin Protein 1 gene family over 40 million years of Drosophila evolution. Our report expands this family from a modest five genes to 26 genes. Unlike the founding family members, the HP1s we describe are structurally diverse, largely restricted to male reproductive tissue, and highly dynamic over evolutionary time. Despite recurrent HP1 gene birth and death, gene numbers per species are relatively constant. These gene “replacements” likely support a dynamic biological process. We propose, and present evidence for, the hypothesis that recurrent chromosomal rearrangements drive at least some HP1 gene family dynamics observed. We anticipate that these HP1 genes will help define new heterochromatin-dependent processes in the male germline.

Parallel reduction in expression of the eye development gene hedgehog in separately derived cave populations of the amphipod Gammarus minus

Caves provide excellent settings to examine evolutionary questions. Subterranean environments are characterized by similar and consistent conditions. Cave-adapted species often share characteristics such as diminished pigmentation, elongated limbs and reduced or absent eyes. Relatively little is known about the evolution and development of troglomorphic traits in invertebrates. In this study, we compare expression of the eye development genes hedgehog, pax6, sine oculis and dachshund in individuals from multiple independently derived cave populations of the amphipod Gammarus minus. hedgehog expression was significantly reduced in cave populations, compared to genetically related surface populations. Interestingly, no differences were found in pax6, sine oculis or dachshund expression. Because hedgehog-related genes are also involved in eye reduced in Astyanax mexicanus, these genes may be consistent targets of evolution during cave adaptation. These results provide support for the hypothesis of genomic 'hotspots' of evolution and allow comparison of adaptive mechanisms among diverse animals in subterranean environments.

Genome features of "Dark-fly", a Drosophila line reared long-term in a dark environment

Organisms are remarkably adapted to diverse environments by specialized metabolisms, morphology, or behaviors. To address the molecular mechanisms underlying environmental adaptation, we have utilized a Drosophila melanogaster line, termed “Dark-fly”, which has been maintained in constant dark conditions for 57 years (1400 generations). We found that Dark-fly exhibited higher fecundity in dark than in light conditions, indicating that Dark-fly possesses some traits advantageous in darkness. Using next-generation sequencing technology, we determined the whole genome sequence of Dark-fly and identified approximately 220,000 single nucleotide polymorphisms (SNPs) and 4,700 insertions or deletions (InDels) in the Dark-fly genome compared to the genome of the Oregon-R-S strain, a control strain. 1.8% of SNPs were classified as non-synonymous SNPs (nsSNPs: i.e., they alter the amino acid sequence of gene products). Among them, we detected 28 nonsense mutations (i.e., they produce a stop codon in the protein sequence) in the Dark-fly genome. These included genes encoding an olfactory receptor and a light receptor. We also searched runs of homozygosity (ROH) regions as putative regions selected during the population history, and found 21 ROH regions in the Dark-fly genome. We identified 241 genes carrying nsSNPs or InDels in the ROH regions. These include a cluster of alpha-esterase genes that are involved in detoxification processes. Furthermore, analysis of structural variants in the Dark-fly genome showed the deletion of a gene related to fatty acid metabolism. Our results revealed unique features of the Dark-fly genome and provided a list of potential candidate genes involved in environmental adaptation.

Tempo and mode in evolution of transcriptional regulation

Perennial questions of evolutionary biology can be applied to gene regulatory systems using the abundance of experimental data addressing gene regulation in a comparative context. What is the tempo (frequency, rate) and mode (way, mechanism) of transcriptional regulatory evolution? Here we synthesize the results of 230 experiments performed on insects and nematodes in which regulatory DNA from one species was used to drive gene expression in another species. General principles of regulatory evolution emerge. Gene regulatory evolution is widespread and accumulates with genetic divergence in both insects and nematodes. Divergence in cis is more common than divergence in trans. Coevolution between cis and trans shows a particular increase over greater evolutionary timespans, especially in sex-specific gene regulation. Despite these generalities, the evolution of gene regulation is gene- and taxon-specific. The congruence of these conclusions with evidence from other types of experiments suggests that general principles are discoverable, and a unified view of the tempo and mode of regulatory evolution may be achievable.

Gene regulatory networks reused to build novel traits: co-option of an eye-related gene regulatory network in eye-like organs and red wing patches on insect wings is suggested by optix expression

Co-option of the eye developmental gene regulatory network may have led to the appearance of novel functional traits on the wings of flies and butterflies. The first trait is a recently described wing organ in a species of extinct midge resembling the outer layers of the midge's own compound eye. The second trait is red pigment patches on Heliconius butterfly wings connected to the expression of an eye selector gene, optix. These examples, as well as others, are discussed regarding the type of empirical evidence and burden of proof that have been used to infer gene network co-option underlying the origin of novel traits. A conceptual framework describing increasing confidence in inference of network co-option is proposed. Novel research directions to facilitate inference of network co-option are also highlighted, especially in cases where the pre-existent and novel traits do not resemble each other.

Quantification of transcript levels with quantitative RT-PCR

Differential gene expression is a key factor driving phenotypic divergence. Determining when and where gene expression has diverged between organisms requires a quantitative method. While large-scale approaches such as microarrays or high-throughput mRNA sequencing can identify candidates, quantitative RT-PCR is the definitive method for confirming gene expression differences. Here, we describe the steps for performing qRT-PCR including extracting total RNA, reverse-transcribing it to make a pool of cDNA, and then quantifying relative expression of a few candidate genes using real-time or quantitative PCR.

Ancestral and conserved cis-regulatory architectures in developmental control genes

Among developmental control genes, transcription factor-target gene "linkages"--the direct connections between target genes and the factors that control their patterns of expression--can show remarkable evolutionary stability. However, the specific binding sites that mediate and define these regulatory connections are themselves often subject to rapid turnover. Here we describe several instances in which particular transcription factor binding motif combinations have evidently been conserved upstream of orthologous target genes for extraordinarily long evolutionary periods. This occurs against a backdrop in which other binding sites for the same factors are coming and going rapidly. Our examples include a particular Dpp Silencer Element upstream of insect brinker genes, in combination with a novel motif we refer to as the Downstream Element; combinations of a Suppressor of Hairless Paired Site (SPS) and a specific proneural protein binding site associated with arthropod Notch pathway target genes; and a three-motif combination, also including an SPS, upstream of deuterostome Hes repressor genes, which are also Notch targets. We propose that these stable motif architectures have been conserved intact from a deep ancestor, in part because they mediate a special mode of regulation that cannot be supplied by the other, unstable motif instances.

Convergent, modular expression of ebony and tan in the mimetic wing patterns of Heliconius butterflies

The evolution of pigmentation in vertebrates and flies has involved repeated divergence at a small number of genes related to melanin synthesis. Here, we study insect melanin synthesis genes in Heliconius butterflies, a group characterised by its diversity of wing patterns consisting of black (melanin), and yellow and red (ommochrome) pigmented scales. Consistent with their respective biochemical roles in Drosophila melanogaster, ebony is upregulated in non-melanic wing regions destined to be pigmented red whilst tan is upregulated in melanic regions. Wing regions destined to be pigmented yellow, however, are downregulated for both genes. This pattern is conserved across multiple divergent and convergent phenotypes within the Heliconii, suggesting a conserved mechanism for the development of black, red and yellow pattern elements across the genus. Linkage mapping of five melanin biosynthesis genes showed that, in contrast to other organisms, these genes do not control pattern polymorphism. Thus, the pigmentation genes themselves are not the locus of evolutionary change but lie downstream of a wing pattern regulatory factor. The results suggest a modular system in which particular combinations of genes are switched on whenever red, yellow or black pattern elements are favoured by natural selection for diverse and mimetic wing patterns.