Intraspecific polymorphism to interspecific divergence: genetics of pigmentation in Drosophila

Genetic Changes to a Transcriptional Silencer Element Confers Phenotypic Diversity within and between Drosophila Species

The modification of transcriptional regulation has become increasingly appreciated as a major contributor to morphological evolution. However, the role of negative-acting control elements (e.g. silencers) in generating morphological diversity has been generally overlooked relative to positive-acting “enhancer” elements. The highly variable body coloration patterns among Drosophilid insects represents a powerful model system in which the molecular alterations that underlie phenotypic diversity can be defined. In a survey of pigment phenotypes among geographically disparate Japanese populations of Drosophila auraria, we discovered a remarkable degree of variation in male-specific abdominal coloration. In testing the expression patterns of the major pigment-producing enzymes, we found that phenotypes uniquely correlated with differences in the expression of ebony, a gene required for yellow-colored cuticle. Assays of ebony’s transcriptional control region indicated that a lightly pigmented strain harbored cis-regulatory mutations that caused correlated changes in its expression. Through a series of chimeric reporter constructs between light and dark strain alleles, we localized function-altering mutations to a conserved silencer that mediates a male-specific pattern of ebony repression. This suggests that the light allele was derived through the loss of this silencer’s activity. Furthermore, examination of the ebony gene of D. serrata, a close relative of D. auraria which secondarily lost male-specific pigmentation revealed the parallel loss of this silencer element. These results demonstrate how loss-of-function mutations in a silencer element resulted in increased gene expression. We propose that the mutational inactivation of silencer elements may represent a favored path to evolve gene expression, impacting morphological traits.

One of the greatest challenges in understanding the relationship between genotype and phenotype is to discern how changes in DNA affect the normal functioning of genes. Mutations may generate a new function for a gene, yet it is frequently observed that they inactivate some aspect of a gene’s normal capacity. Investigations focused on understanding the developmental basis for the evolution of anatomical structures has found a prevalent role for mutations that alter developmental gene regulation. In animals, genes are transcriptionally activated in specific tissues during development by regulatory sequences distributed across their expansive non-protein coding regions. Regulatory elements known as silencers act to prevent genes from being expressed in certain tissues, providing a mechanism for precise control. Here, we show how a silencer that prevents expression of a pigment-producing enzyme in certain Drosophila species has repeatedly been subject to inactivating mutations that increased this gene’s expression. This example illustrates how such negative-acting regulatory sequences can represent a convenient target for increasing gene expression through the loss of a genetic element.

The evolutionary origination and diversification of a dimorphic gene regulatory network through parallel innovations in cis and trans

The origination and diversification of morphological characteristics represents a key problem in understanding the evolution of development. Morphological traits result from gene regulatory networks (GRNs) that form a web of transcription factors, which regulate multiple cis-regulatory element (CRE) sequences to control the coordinated expression of differentiation genes. The formation and modification of GRNs must ultimately be understood at the level of individual regulatory linkages (i.e., transcription factor binding sites within CREs) that constitute the network. Here, we investigate how elements within a network originated and diversified to generate a broad range of abdominal pigmentation phenotypes among Sophophora fruit flies. Our data indicates that the coordinated expression of two melanin synthesis enzymes, Yellow and Tan, recently evolved through novel CRE activities that respond to the spatial patterning inputs of Hox proteins and the sex-specific input of Bric-à-brac transcription factors. Once established, it seems that these newly evolved activities were repeatedly modified by evolutionary changes in the network’s trans-regulators to generate large-scale changes in pigment pattern. By elucidating how yellow and tan are connected to the web of abdominal trans-regulators, we discovered that the yellow and tan abdominal CREs are composed of distinct regulatory inputs that exhibit contrasting responses to the same Hox proteins and Hox cofactors. These results provide an example in which CRE origination underlies a recently evolved novel trait, and highlights how coordinated expression patterns can evolve in parallel through the generation of unique regulatory linkages.

The genomic content of regulatory genes such as transcription factors is surprisingly conserved between diverse animal species, raising the paradox of how new traits emerge, and are subsequently modified and lost. In this study we make a connection between the developmental basis for the formation of a fruit fly trait and the evolutionary basis for that trait’s origin, diversification, and loss. We show how the origin of a novel pigmentation trait is associated with the evolution of two regulatory sequences that control the co-expression of two key pigmentation genes. These sequences interact in unique ways with evolutionarily conserved Hox transcription factors to drive gene co-expression. Once these unique connections evolved, the alteration of this trait appears to have proceeded through changes to regulatory genes rather than regulatory sequences of the pigmentation genes. Thus, our findings support a scenario where regulatory sequence evolution provided new functions to old transcription factors, how co-expression can emerge from different utilizations of the same transcription factors, and that trait diversity was surprisingly shaped by changes in some manner to the deeply conserved regulatory genes.

A single gene causes an interspecific difference in pigmentation in Drosophila

The genetic basis of species differences remains understudied. Studies in insects have contributed significantly to our understanding of morphological evolution. Pigmentation traits in particular have received a great deal of attention and several genes in the insect pigmentation pathway have been implicated in inter- and intraspecific differences. Nonetheless, much remains unknown about many of the genes in this pathway and their potential role in understudied taxa. Here we genetically analyze the puparium color difference between members of the virilis group of Drosophila. The puparium of Drosophila virilis is black, while those of D. americana, D. novamexicana, and D. lummei are brown. We used a series of backcross hybrid populations between D. americana and D. virilis to map the genomic interval responsible for the difference between this species pair. First, we show that the pupal case color difference is caused by a single Mendelizing factor, which we ultimately map to an ∼11-kb region on chromosome 5. The mapped interval includes only the first exon and regulatory region(s) of the dopamine N-acetyltransferase gene (Dat). This gene encodes an enzyme that is known to play a part in the insect pigmentation pathway. Second, we show that this gene is highly expressed at the onset of pupation in light brown taxa (D. americana and D. novamexicana) relative to D. virilis, but not in the dark brown D. lummei. Finally, we examine the role of Dat in adult pigmentation between D. americana (heavily melanized) and D. novamexicana (lightly melanized) and find no discernible effect of this gene in adults. Our results demonstrate that a single gene is entirely or almost entirely responsible for a morphological difference between species.

A new arylalkylamine N-acetyltransferase in silkworm (Bombyx mori) affects integument pigmentation

Dopamine is a precursor for melanin synthesis. Arylalkylamine N-acetyltransferase (AANAT) is involved in the melatonin formation in insects because it could catalyze the transformation from dopamine to dopamine-N-acetyldopamine. In this study, we identified a new AANAT gene in the silkworm (Bombyx mori) and assessed its role in the silkworm. The cDNA of this gene encodes 233 amino acids that shares 57 % amino acid identity with the Bm-iAANAT protein. We thus refer to this gene as Bm-iAANAT2. To investigate the role of Bm-iAANAT2, we constructed a transgenic interference system using a 3xp3 promoter to suppress the expression of Bm-iAANAT2 in the silkworm. We observed that melanin deposition occurs in the head and integument in transgenic lines. To verify the melanism pattern, dopamine content and the enzyme activity of AANAT were determined by high-performance liquid chromatography (HPLC). We found that an increase in dopamine levels affects melanism patterns on the heads of transgenic B. mori. A reduction in the enzyme activity of AANAT leads to changes in dopamine levels. We analyzed the expression of the Bm-iAANAT2 genes by qPCR and found that the expression of Bm-iAANAT2 gene is significantly lower in transgenic lines. Our results lead us to conclude that Bm-iAANAT2 is a new arylalkylamine N-acetyltransferase gene in the silkworm and is involved in the metabolism of the dopamine to avoid the generation of melanin.

Intraspecific polymorphism, interspecific divergence, and the origins of function-altering mutations in deer mouse hemoglobin

Major challenges for illuminating the genetic basis of phenotypic evolution are to identify causative mutations, to quantify their functional effects, to trace their origins as new or preexisting variants, and to assess the manner in which segregating variation is transduced into species differences. Here, we report an experimental analysis of genetic variation in hemoglobin (Hb) function within and among species of Peromyscus mice that are native to different elevations. A multilocus survey of sequence variation in the duplicated HBA and HBB genes in Peromyscus maniculatus revealed that function-altering amino acid variants are widely shared among geographically disparate populations from different elevations, and numerous amino acid polymorphisms are also shared with closely related species. Variation in Hb-O₂ affinity within and among populations of P. maniculatus is attributable to numerous amino acid mutations that have individually small effects. One especially surprising feature of the Hb polymorphism in P. maniculatus is that an appreciable fraction of functional standing variation in the two transcriptionally active HBA paralogs is attributable to recurrent gene conversion from a tandemly linked HBA pseudogene. Moreover, transpecific polymorphism in the duplicated HBA genes is not solely attributable to incomplete lineage sorting or introgressive hybridization; instead, it is mainly attributable to recurrent interparalog gene conversion that has occurred independently in different species. Partly as a result of concerted evolution between tandemly duplicated globin genes, the same amino acid changes that contribute to variation in Hb function within P. maniculatus also contribute to divergence in Hb function among different species of Peromyscus. In the case of function-altering Hb mutations in Peromyscus, there is no qualitative or quantitative distinction between segregating variants within species and fixed differences between species.

Extensive transcriptional response associated with seasonal plasticity of butterfly wing patterns

In the eastern United States, the buckeye butterfly, Junonia coenia, shows seasonal wing colour plasticity where adults emerging in the spring are tan, while those emerging in the autumn are dark red. This variation can be artificially induced in laboratory colonies, thus making J. coenia a useful model system to examine the mechanistic basis of plasticity. To better understand the developmental basis of seasonal plasticity, we used RNA-seq to quantify transcription profiles associated with development of alternative seasonal wing morphs. Depending on the developmental stage, between 547 and 1420 transfrags were significantly differentially expressed between morphs. These extensive differences in gene expression stand in contrast to the much smaller numbers of differentially expressed transcripts identified in previous studies of genetic wing pattern variation in other species and suggest that environmentally induced phenotypic shifts arise from very broad systemic processes. Analyses of candidate endocrine and pigmentation transcripts revealed notable genes upregulated in the red morph, including several ecdysone-associated genes, and cinnabar, an ommochrome pigmentation gene implicated in colour pattern variation in other butterflies. We also found multiple melanin-related transcripts strongly upregulated in the red morph, including tan and yellow-family genes, leading us to speculate that dark red pigmentation in autumn J. coenia may involve nonommochrome pigments. While we identified several endocrine and pigmentation genes as obvious candidates for seasonal colour morph differentiation, we speculate that the majority of observed expression differences were due to thermal stress response. The buckeye transcriptome provides a basis for further developmental studies of phenotypic plasticity.

A flexible Bayesian method for detecting allelic imbalance in RNA-seq data

Background

One method of identifying cis regulatory differences is to analyze allele-specific expression (ASE) and identify cases of allelic imbalance (AI). RNA-seq is the most common way to measure ASE and a binomial test is often applied to determine statistical significance of AI. This implicitly assumes that there is no bias in estimation of AI. However, bias has been found to result from multiple factors including: genome ambiguity, reference quality, the mapping algorithm, and biases in the sequencing process. Two alternative approaches have been developed to handle bias: adjusting for bias using a statistical model and filtering regions of the genome suspected of harboring bias. Existing statistical models which account for bias rely on information from DNA controls, which can be cost prohibitive for large intraspecific studies. In contrast, data filtering is inexpensive and straightforward, but necessarily involves sacrificing a portion of the data.

Results

Here we propose a flexible Bayesian model for analysis of AI, which accounts for bias and can be implemented without DNA controls. In lieu of DNA controls, this Poisson-Gamma (PG) model uses an estimate of bias from simulations. The proposed model always has a lower type I error rate compared to the binomial test. Consistent with prior studies, bias dramatically affects the type I error rate. All of the tested models are sensitive to misspecification of bias. The closer the estimate of bias is to the true underlying bias, the lower the type I error rate. Correct estimates of bias result in a level alpha test.

Conclusions

To improve the assessment of AI, some forms of systematic error (e.g., map bias) can be identified using simulation. The resulting estimates of bias can be used to correct for bias in the PG model, without data filtering. Other sources of bias (e.g., unidentified variant calls) can be easily captured by DNA controls, but are missed by common filtering approaches. Consequently, as variant identification improves, the need for DNA controls will be reduced. Filtering does not significantly improve performance and is not recommended, as information is sacrificed without a measurable gain. The PG model developed here performs well when bias is known, or slightly misspecified. The model is flexible and can accommodate differences in experimental design and bias estimation.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-920) contains supplementary material, which is available to authorized users.

Ancient homology underlies adaptive mimetic diversity across butterflies

Convergent evolution provides a rare, natural experiment with which to test the predictability of adaptation at the molecular level. Little is known about the molecular basis of convergence over macro-evolutionary timescales. Here we use a combination of positional cloning, population genomic resequencing, association mapping and developmental data to demonstrate that positionally orthologous nucleotide variants in the upstream region of the same gene, WntA, are responsible for parallel mimetic variation in two butterfly lineages that diverged >65 million years ago. Furthermore, characterization of spatial patterns of WntA expression during development suggests that alternative regulatory mechanisms underlie wing pattern variation in each system. Taken together, our results reveal a strikingly predictable molecular basis for phenotypic convergence over deep evolutionary time.

Little is known about the genetic basis of convergent evolution in deeply diverged species. Here, the authors show that variation in the WntA gene is associated with parallel wing pattern variation in two butterflies that diverged more than 65 million years ago.

Interspecific tests of allelism reveal the evolutionary timing and pattern of accumulation of reproductive isolation mutations

Despite extensive theory, little is known about the empirical accumulation and evolutionary timing of mutations that contribute to speciation. Here we combined QTL (Quantitative Trait Loci) analyses of reproductive isolation, with information on species evolutionary relationships, to reconstruct the order and timing of mutations contributing to reproductive isolation between three plant (Solanum) species. To evaluate whether reproductive isolation QTL that appear to coincide in more than one species pair are homologous, we used cross-specific tests of allelism and found evidence for both homologous and lineage-specific (non-homologous) alleles at these co-localized loci. These data, along with isolation QTL unique to single species pairs, indicate that >85% of isolation-causing mutations arose later in the history of divergence between species. Phylogenetically explicit analyses of these data support non-linear models of accumulation of hybrid incompatibility, although the specific best-fit model differs between seed (pairwise interactions) and pollen (multi-locus interactions) sterility traits. Our findings corroborate theory that predicts an acceleration ('snowballing') in the accumulation of isolation loci as lineages progressively diverge, and suggest different underlying genetic bases for pollen versus seed sterility. Pollen sterility in particular appears to be due to complex genetic interactions, and we show this is consistent with a snowball model where later arising mutations are more likely to be involved in pairwise or multi-locus interactions that specifically involve ancestral alleles, compared to earlier arising mutations.

Pigmentation in Drosophila melanogaster reaches its maximum in Ethiopia and correlates most strongly with ultra-violet radiation in sub-Saharan Africa

Background

Pigmentation has a long history of investigation in evolutionary biology. In Drosophila melanogaster, latitudinal and altitudinal clines have been found but their underlying causes remain unclear. Moreover, most studies were conducted on cosmopolitan populations which have a relatively low level of genetic structure and diversity compared to sub-Saharan African populations. We investigated: 1) the correlation between pigmentation traits within and between the thorax and the fourth abdominal segment, and 2) their associations with different geographical and ecological variables, using 710 lines belonging to 30 sub-Saharan and cosmopolitan populations.

Results

Pigmentation clines substantially differed between sub-Saharan and cosmopolitan populations. While positive correlations with latitude have previously been described in Europe, India and Australia, in agreement with Bogert's rule or the thermal melanism hypothesis, we found a significant negative correlation in Africa. This correlation persisted even after correction for altitude, which in its turn showed a positive correlation with pigmentation independently from latitude. More importantly, we found that thoracic pigmentation reaches its maximal values in this species in high-altitude populations of Ethiopia (1,600-3,100 m). Ethiopian flies have a diffuse wide thoracic trident making the mesonotum and the head almost black, a phenotype that is absent from all other sub-Saharan or cosmopolitan populations including high-altitude flies from Peru (~3,400 m). Ecological analyses indicated that the variable most predictive of pigmentation in Africa, especially for the thorax, was ultra-violet (UV) intensity, consistent with the so-called Gloger's rule invoking a role of melanin in UV protection.

Conclusion

Our data suggest that different environmental factors may shape clinal variation in tropical and temperate regions, and may lead to the evolution of different degrees of melanism in different high altitude populations in the tropics.

Cardamine hirsuta: a versatile genetic system for comparative studies

A major goal in biology is to identify the genetic basis for phenotypic diversity. This goal underpins research in areas as diverse as evolutionary biology, plant breeding and human genetics. A limitation for this research is no longer the availability of sequence information but the development of functional genetic tools to understand the link between changes in sequence and phenotype. Here we describe Cardamine hirsuta, a close relative of the reference plant Arabidopsis thaliana, as an experimental system in which genetic and transgenic approaches can be deployed effectively for comparative studies. We present high-resolution genetic and cytogenetic maps for C. hirsuta and show that the genome structure of C. hirsuta closely resembles the eight chromosomes of the ancestral crucifer karyotype and provides a good reference point for comparative genome studies across the Brassicaceae. We compared morphological and physiological traits between C. hirsuta and A. thaliana and analysed natural variation in stamen number in which lateral stamen loss is a species characteristic of C. hirsuta. We constructed a set of recombinant inbred lines and detected eight quantitative trait loci that can explain stamen number variation in this population. We found clear phylogeographic structure to the genetic variation in C. hirsuta, thus providing a context within which to address questions about evolutionary changes that link genotype with phenotype and the environment.

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.

Divergence of iron metabolism in wild Malaysian yeast

Comparative genomic studies have reported widespread variation in levels of gene expression within and between species. Using these data to infer organism-level trait divergence has proven to be a key challenge in the field. We have used a wild Malaysian population of S. cerevisiae as a test bed in the search to predict and validate trait differences based on observations of regulatory variation. Malaysian yeast, when cultured in standard medium, activated regulatory programs that protect cells from the toxic effects of high iron. Malaysian yeast also showed a hyperactive regulatory response during culture in the presence of excess iron and had a unique growth defect in conditions of high iron. Molecular validation experiments pinpointed the iron metabolism factors AFT1, CCC1, and YAP5 as contributors to these molecular and cellular phenotypes; in genome-scale sequence analyses, a suite of iron toxicity response genes showed evidence for rapid protein evolution in Malaysian yeast. Our findings support a model in which iron metabolism has diverged in Malaysian yeast as a consequence of a change in selective pressure, with Malaysian alleles shifting the dynamic range of iron response to low-iron concentrations and weakening resistance to extreme iron toxicity. By dissecting the iron scarcity specialist behavior of Malaysian yeast, our work highlights the power of expression divergence as a signpost for biologically and evolutionarily relevant variation at the organismal level. Interpreting the phenotypic relevance of gene expression variation is one of the primary challenges of modern genomics.

The genetic basis of a rare flower color polymorphism in Mimulus lewisii provides insight into the repeatability of evolution

A long-standing question in evolutionary biology asks whether the genetic changes contributing to phenotypic evolution are predictable. Here, we identify a genetic change associated with segregating variation in flower color within a population of Mimulus lewisii. To determine whether these types of changes are predictable, we combined this information with data from other species to investigate whether the spectrum of mutations affecting flower color transitions differs based on the evolutionary time-scale since divergence. We used classic genetic techniques, along with gene expression and population genetic approaches, to identify the putative, loss-of-function mutation that generates rare, white flowers instead of the common, pink color in M. lewisii. We found that a frameshift mutation in an anthocyanin pathway gene is responsible for the white-flowered polymorphism found in this population of M. lewisii. Comparison of our results with data from other species reveals a broader spectrum of flower color mutations segregating within populations relative to those that fix between populations. These results suggest that the genetic basis of fixed differences in flower color may be predictable, but that for segregating variation is not.

Contemporary climate change and terrestrial invertebrates: evolutionary versus plastic changes

To forecast the responses of species to future climate change, an understanding of the ability of species to adapt to long-term shifts in temperature is crucial. We present a review on evolutionary adaptation and phenotypic plasticity of temperature-related traits in terrestrial invertebrates. The evidence for adaptive evolution in melanization is good, but we caution that genetic determination needs to be tested in each individual species, and complex genetic correlations may exist. For phenological traits allochronic data sets provide powerful means to track climate-induced changes; however, rarely are responses deconstructed into evolutionary and plastic responses. Laboratory studies suggest climate change responses in these traits will be driven by both. For stress resistance, the evidence for shifts in traits is poor. Studies leaning heavily on Drosophila have demonstrated potential limits to evolutionary responses in desiccation and heat resistance. Quantifying the capacity for these species to respond plastically and extending this work to other taxa will be an important next step. We also note that, although not strictly speaking a species trait, the response of endosymbionts to heat stress requires further study. Finally, while clearly genetic, and possibly adaptive, the anonymous nature of latitudinal shifts in clines of genetic markers in Drosophila prevents further interpretation.

A perspective on micro-evo-devo: progress and potential

The term "micro-evo-devo" refers to the combined study of the genetic and developmental bases of natural variation in populations and the evolutionary forces that have shaped this variation. It thus represents a synthesis of the fields of evolutionary developmental biology and population genetics. As has been pointed out by several others, this synthesis can provide insights into the evolution of organismal form and function that have not been possible within these individual disciplines separately. Despite a number of important successes in micro-evo-devo, however, it appears that evo devo and population genetics remain largely separate spheres of research, limiting their ability to address evolutionary questions. This also risks pushing contemporary evo devo to the fringes of evolutionary biology because it does not describe the causative molecular changes underlying evolution or the evolutionary forces involved. Here we reemphasize the theoretical and practical importance of micro-evo-devo as a strategy for understanding phenotypic evolution, review the key recent insights that it has provided, and present a perspective on both the potential and the remaining challenges of this exciting interdisciplinary field.

Drosophila americana as a model species for comparative studies on the molecular basis of phenotypic variation

Understanding the molecular basis of within and between species phenotypic variation is one of the main goals of Biology. In Drosophila, most of the work regarding this issue has been performed in D. melanogaster, but other distantly related species must also be studied to verify the generality of the findings obtained for this species. Here, we make the case for D. americana, a species of the virilis group of Drosophila that has been diverging from the model species, D. melanogaster, for approximately 40 Myr. To determine the suitability of this species for such studies, polymorphism and recombination estimates are presented for D. americana based on the largest nucleotide sequence polymorphism data set so far analyzed (more than 100 data sets) for this species. The polymorphism estimates are also compared with those obtained from the comparison of the genome assembly of two D. americana strains (H5 and W11) here reported. As an example of the general utility of these resources, we perform a preliminary study on the molecular basis of lifespan differences in D. americana. First, we show that there are lifespan differences between D. americana populations from different regions of the distribution range. Then, we perform five F2 association experiments using markers for 21 candidate genes previously identified in D. melanogaster. Significant associations are found between polymorphism at two genes (hep and Lim3) and lifespan. For the F2 association study involving the two sequenced strains (H5 and W11), we identify amino acid differences at Lim3 and Hep that could be responsible for the observed changes in lifespan. For both genes, no large gene expression differences were observed between the two strains.

Flower color as a model system for studies of plant evo-devo

Even though pigmentation traits have had substantial impacts on the field of animal evolutionary developmental biology, they have played only relatively minor roles in plant evo-devo. This is surprising given the often direct connection between flower color and fitness variation mediated through the effects of pollinators. At the same time, ecological and evolutionary genetic studies have utilized the molecular resources available for the anthocyanin pathway to generate several examples of the molecular basis of putatively adaptive transitions in flower color. Despite this opportunity to synthesize experimental approaches in ecology, evolution, and developmental biology, the investigation of many fundamental questions in evo-devo using this powerful model is only at its earliest stages. For example, a long-standing question is whether predictable genetic changes accompany the repeated evolution of a trait. Due to the conserved nature of the biochemical and regulatory control of anthocyanin biosynthesis, it has become possible to determine whether, and under what circumstances, different types of mutations responsible for flower color variation are preferentially targeted by natural selection. In addition, because plants use anthocyanin and related compounds in vegetative tissue for other important physiological functions, the identification of naturally occurring transitions from unpigmented to pigmented flowers provides the opportunity to examine the mechanisms by which regulatory networks are co-opted into new developmental domains. Here, we review what is known about the ecological and molecular basis of anthocyanic flower color transitions in natural systems, focusing on the evolutionary and developmental features involved. In doing so, we provide suggestions for future work on this trait and suggest that there is still much to be learned from the evolutionary development of flower color transitions in nature.

Combining genome-wide methods to investigate the genetic complexity of courtship song variation in Drosophila melanogaster

Little is currently known about the genetic complexity of quantitative behavioral variation, the types of genes involved, or their effects on intermediate phenotypes. Here, we conduct a genome-wide association study of Drosophila melanogaster courtship song variation using 168 sequenced inbred lines, and fail to find highly significant associations. However, by combining these data with results from a well-powered Evolve and Resequence (E&R) study on the same trait, we provide statistical evidence that some power to associate genotype and phenotype is available. Genes that are significant in both analyses are enriched for expression in the nervous system, and affect neural development and synaptic growth when perturbed. Quantitative complementation at one of these loci, Syntrophin-like 1, supports a hypothesis that variation at this locus affects variation in the inter-pulse interval of courtship song. These results suggest that experimental evolution may provide an approach for genome-scale replication in Drosophila.

Pupariation site preference within and between Drosophila sibling species

Holometabolous insects pass through a sedentary pupal stage and often choose a location for pupation that is different from the site of larval feeding. We have characterized a difference in pupariation site choice within and between sibling species of Drosophila. We found that, in nature, Drosophila sechellia pupariate within their host fruit, Morinda citrifolia, and that they perform this behavior in laboratory assays. In contrast, in the laboratory, geographically diverse strains of Drosophila simulans vary in their pupariation site preference; D. simulans lines from the ancestral range in southeast Africa pupariate on fruit, or a fruit substitute, whereas populations from Europe or the New World select sites off of fruit. We explored the genetic basis for the evolved preference in puariation site preference by performing quantitative trait locus mapping within and between species. We found that the interspecific difference is controlled largely by loci on chromosomes X and II. In contrast, variation between two strains of D. simulans appears to be highly polygenic, with the majority of phenotypic effects due to loci on chromosome III. These data address the genetic basis of how new traits arise as species diverge and populations disperse.

Human genetic variation within neural crest enhancers: molecular and phenotypic implications

Developmental gene expression programmes are coordinated by the specialized distal cis-regulatory elements called enhancers, which integrate lineage- and signalling-dependent inputs to guide morphogenesis. In previous work, we characterized the genome-wide repertoire of active enhancers in human neural crest cells (hNCC), an embryonic cell population with critical roles in craniofacial development. We showed that in hNCC, co-occupancy of a master regulator TFAP2A with nuclear receptors NR2F1 and NR2F2 correlates with the presence of permissive enhancer chromatin states. Here, we take advantage of pre-existing human genetic variation to further explore potential cooperation between TFAP2A and NR2F1/F2. We demonstrate that isolated single nucleotide polymorphisms affecting NR2F1/F2-binding sites within hNCC enhancers can alter TFAP2A occupancy and overall chromatin features at the same enhancer allele. We propose that a similar strategy can be used to elucidate other cooperative relationships between transcription factors involved in developmental transitions. Using the neural crest and its major contribution to human craniofacial phenotypes as a paradigm, we discuss how genetic variation might modulate the molecular properties and activity of enhancers, and ultimately impact human phenotypic diversity.

The influence of abdominal pigmentation on desiccation and ultraviolet resistance in two species of Drosophila

Drosophila yakuba and D. santomea are sister species that differ in their levels of abdominal pigmentation; D. yakuba shows heavily pigmented posterior abdominal segments in both sexes, whereas D. santomea lacks dark pigment anywhere on its body. Using naturally collected lines, we demonstrate the existence of altitudinal variation in abdominal pigmentation in D. yakuba but not in D. santomea. We use the variation in pigmentation within D. yakuba and two body-color mutants in D. yakuba to elucidate selective advantage of differences in pigmentation. Our results indicate that although differences in abdominal pigmentation have no effect on desiccation resistance, lighter pigmentation confers ultraviolet radiation resistance in this pair of species.

The Loci of repeated evolution: a catalog of genetic hotspots of phenotypic variation

What is the nature of the genetic changes underlying phenotypic evolution? We have catalogued 1008 alleles described in the literature that cause phenotypic differences among animals, plants, and yeasts. Surprisingly, evolution of similar traits in distinct lineages often involves mutations in the same gene ("gene reuse"). This compilation yields three important qualitative implications about repeated evolution. First, the apparent evolution of similar traits by gene reuse can be traced back to two alternatives, either several independent causative mutations or a single original mutational event followed by sorting processes. Second, hotspots of evolution-defined as the repeated occurrence of de novo mutations at orthologous loci and causing similar phenotypic variation-are omnipresent in the literature with more than 100 examples covering various levels of analysis, including numerous gain-of-function events. Finally, several alleles of large effect have been shown to result from the aggregation of multiple small-effect mutations at the same hotspot locus, thus reconciling micromutationist theories of adaptation with the empirical observation of large-effect variants. Although data heterogeneity and experimental biases prevented us from extracting quantitative trends, our synthesis highlights the existence of genetic paths of least resistance leading to viable evolutionary change.

Rapid, habitat-related evolution of land snail colour morphs on reclaimed land

I made use of the known dates of reclamation (and of afforestations) in the IJsselmeerpolders in The Netherlands to assess evolutionary adaptation in Cepaea nemoralis. At 12 localities (three in each polder), I sampled a total of 4390 adult individuals in paired open and shaded habitats, on average 233 m apart, and scored these for genetic shell colour polymorphisms. The results show (highly) significant differentiation at most localities, although the genes involved differed per locality. Overall, though, populations in shaded habitats had evolved towards darker shells than those in adjacent open habitats, whereas a 'Cain & Sheppard' diagram (proportion yellow shells plotted against 'effectively unbanded' shells) failed to reveal a clear pattern. This might suggest that thermal selection is more important than visual selection in generating this pattern. Trait differentiation, regardless of whether they were plotted against polder age or habitat age, showed a linear increase of differentiation with time, corresponding to a mean rate of trait evolution of 15-31 kilodarwin. In conclusion, C. nemoralis is capable of rapid and considerable evolutionary differentiation over 1-25 snail generations, though equilibrium may be reached only at longer time scales.

Dissecting sources of quantitative gene expression pattern divergence between Drosophila species

Gene expression patterns can diverge between species due to changes in a gene's regulatory DNA or changes in the proteins, e.g., transcription factors (TFs), that regulate the gene. We developed a modeling framework to uncover the sources of expression differences in blastoderm embryos of three Drosophila species, focusing on the regulatory circuit controlling expression of the hunchback (hb) posterior stripe. Using this framework and cellular-resolution expression measurements of hb and its regulating TFs, we found that changes in the expression patterns of hb's TFs account for much of the expression divergence. We confirmed our predictions using transgenic D. melanogaster lines, which demonstrate that this set of orthologous cis-regulatory elements (CREs) direct similar, but not identical, expression patterns. We related expression pattern differences to sequence changes in the CRE using a calculation of the CRE's TF binding site content. By applying this calculation in both the transgenic and endogenous contexts, we found that changes in binding site content affect sensitivity to regulating TFs and that compensatory evolution may occur in circuit components other than the CRE.

One quantitative trait locus for intra- and interspecific variation in a sex pheromone

Even though premating isolation is hypothesized to be a major driving force in speciation, its genetic basis is poorly known. In the noctuid moth Heliothis subflexa, one group of sex pheromone components, the acetates, emitted by the female, plays a crucial isolating role in preventing interspecific matings to males of the closely related Heliothis virescens, in which females do not produce acetates and males are repelled by them. We previously found intraspecific variation in acetates in H. subflexa: females in eastern North America contain significantly more acetates than females in Western Mexico. Here we describe the persistence of this intraspecific variation in laboratory-reared strains and the identification of one major quantitative trait locus (QTL), explaining 40% of the variance in acetate amounts. We homologized this intraspecific QTL to our previously identified interspecific QTL using restriction-associated DNA (RAD) tags. We found that a major intraspecific QTL overlaps with one of the two major interspecific QTL. To identify candidate genes underlying the acetate variation, we investigated a number of gene families with known or suspected acetyl- or acyltransferase activity. The most likely candidate genes did not map to our QTL, so that we currently hypothesize that a transcription factor underlies this QTL. Finding a single, large QTL that impacts variation in pheromone blends between and within species is, to our knowledge, the first such example for traits that have been demonstrated to affect premating isolation.

Noncompetitive Gametic Isolation between Sibling Species of Cricket: A Hypothesized Link between Within-Population Incompatibility and Reproductive Isolation between Species

Postmating, prezygotic phenotypes are a common mechanism of reproductive isolation. Here, we describe the dynamics of a noncompetitive gametic isolation phenotype (namely, the ability of a male to induce a female to lay eggs) in a group of recently diverged crickets that are primarily isolated from each other by this phenotype. We not only show that heterospecific males are less able to induce females to lay eggs but that there are male by female incompatibilities in this phenotype that occur within populations. We also identify a protein in the female reproductive tract that correlates with the number of eggs that she was induced to lay. Functional genetic tests using RNAi confirm that the function of this protein is linked to egg-laying induction. Moreover, the dysfunction of this protein appears to underlie both within-population incompatibilities and between-species divergence—thus suggesting a common genetic pathway underlies both. However, this is only correlative evidence and further research is needed to assess whether or not the same mutations in the same genes underlie variation at both levels.

The ontogeny of color: developmental origins of divergent pigmentation in Drosophila americana and D. novamexicana

Pigmentation is a model trait for evolutionary and developmental analysis that is particularly amenable to molecular investigation in the genus Drosophila. To better understand how this phenotype evolves, we examined divergent pigmentation and gene expression over developmental time in the dark-bodied D. americana and its light-bodied sister species D. novamexicana. Prior genetic analysis implicated two enzyme-encoding genes, tan and ebony, in pigmentation divergence between these species, but questions remain about the underlying molecular mechanisms. Here, we describe stages of pupal development in both species and use this staging to determine when pigmentation develops and diverges between D. americana and D. novamexicana. For the developmental stages encompassing pigment divergence, we compare mRNA expression of tan and ebony over time and between species. Finally, we use allele-specific expression assays to determine whether interspecific differences in mRNA abundance have a cis-regulatory basis and find evidence of cis-regulatory divergence for both tan and ebony. cis-regulatory divergence affecting tan had a small effect on mRNA abundance and was limited to a few developmental stages, yet previous data suggests that this divergence is likely to be biologically meaningful. Our study suggests that small and developmentally transient expression changes may contribute to phenotypic diversification more often than commonly appreciated. Recognizing the potential phenotypic impact of such changes is important for a scientific community increasingly focused on dissecting quantitative variation, but detecting these types of changes will be a major challenge to elucidating the molecular basis of complex traits.

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.

The genetic architecture of growth rate in juvenile Takifugu species

Closely related species have often evolved dramatic differences in body size. Takifugu rubripes (fugu) is a large marine pufferfish whose genome has been sequenced, whereas T. niphobles is the smallest species among Takifugu. We show that, unsurprisingly, the juvenile growth rate of T. rubripes is higher than that of T. niphobles in a laboratory setting. We produced F(2) progenies of their F(1) hybrids and found one quantitative trait locus (QTL) significantly associated with variation in juvenile body size. This QTL region (3.5 Mb) contains no known genes directly related to growth phenotype (such as IGFs) except Fgf21, which inhibits growth hormone signaling in mouse. The QTL in Takifugu spp. is distinct from the region previously known to control body size variations in stickleback or tilapia. Our results suggest that in the fish tested herein, genomic regions underlying body size evolution might have different genetic origins. They also suggest that many diverse traits in Takifugu spp. are amenable to genetic mapping.

Effects of altered catecholamine metabolism on pigmentation and physical properties of sclerotized regions in the silkworm melanism mutant

Catecholamine metabolism plays an important role in the determination of insect body color and cuticle sclerotization. To date, limited research has focused on these processes in silkworm. In the current study, we analyzed the interactions between catecholamines and melanin genes and their effects on the pigmentation patterns and physical properties of sclerotized regions in silkworm, using the melanic mutant melanism (mln) silkworm strain as a model. Injection of β-alanine into mln mutant silkworm induced a change in catecholamine metabolism and turned its body color yellow. Further investigation of the catecholamine content and expression levels of the corresponding melanin genes from different developmental stages of Dazao-mln (mutant) and Dazao (wild-type) silkworm revealed that at the larval and adult stages, the expression patterns of melanin genes precipitated dopamine accumulation corresponding to functional loss of Bm-iAANAT, a repressive effect of excess NBAD on ebony, and upregulation of tan in the Dazao-mln strain. During the early pupal stage, dopamine did not accumulate in Dazao-mln, since upregulation of ebony and black genes led to conversion of high amounts of dopamine into NBAD, resulting in deep yellow cuticles. Scanning electron microscope analysis of a cross-section of adult dorsal plates from both wild-type and mutant silkworm disclosed the formation of different layers in Dazao-mln owing to lack of NADA, compared to even and dense layers in Dazao. Analysis of the mechanical properties of the anterior wings revealed higher storage modulus and lower loss tangent in Dazao-mln, which was closely associated with the altered catecholamine metabolism in the mutant strain. Based on these findings, we conclude that catecholamine metabolism is crucial for the color pattern and physical properties of cuticles in silkworm. Our results should provide a significant contribution to Lepidoptera cuticle tanning research.

Extensive compensatory cis-trans regulation in the evolution of mouse gene expression

Gene expression levels are thought to diverge primarily via regulatory mutations in trans within species, and in cis between species. To test this hypothesis in mammals we used RNA-sequencing to measure gene expression divergence between C57BL/6J and CAST/EiJ mouse strains and allele-specific expression in their F1 progeny. We identified 535 genes with parent-of-origin specific expression patterns, although few of these showed full allelic silencing. This suggests that the number of imprinted genes in a typical mouse somatic tissue is relatively small. In the set of nonimprinted genes, 32% showed evidence of divergent expression between the two strains. Of these, 2% could be attributed purely to variants acting in trans, while 43% were attributable only to variants acting in cis. The genes with expression divergence driven by changes in trans showed significantly higher sequence constraint than genes where the divergence was explained by variants acting in cis. The remaining genes with divergent patterns of expression (55%) were regulated by a combination of variants acting in cis and variants acting in trans. Intriguingly, the changes in expression induced by the cis and trans variants were in opposite directions more frequently than expected by chance, implying that compensatory regulation to stabilize gene expression levels is widespread. We propose that expression levels of genes regulated by this mechanism are fine-tuned by cis variants that arise following regulatory changes in trans, suggesting that many cis variants are not the primary targets of natural selection.

Genotype-phenotype mapping in a post-GWAS world

Understanding how metabolic reactions, cell signaling, and developmental pathways translate the genome of an organism into its phenotype is a grand challenge in biology. Genome-wide association studies (GWAS) statistically connect genotypes to phenotypes, without any recourse to known molecular interactions, whereas a molecular biology approach directly ties gene function to phenotype through gene regulatory networks (GRNs). Using natural variation in allele-specific expression, GWAS and GRN approaches can be merged into a single framework via structural equation modeling (SEM). This approach leverages the myriad of polymorphisms in natural populations to elucidate and quantitate the molecular pathways that underlie phenotypic variation. The SEM framework can be used to quantitate a GRN, evaluate its consistency across environments or sexes, identify the differences in GRNs between species, and annotate GRNs de novo in non-model organisms.

Transcriptome analysis reveals novel patterning and pigmentation genes underlying Heliconius butterfly wing pattern variation

Background

Heliconius butterfly wing pattern diversity offers a unique opportunity to investigate how natural genetic variation can drive the evolution of complex adaptive phenotypes. Positional cloning and candidate gene studies have identified a handful of regulatory and pigmentation genes implicated in Heliconius wing pattern variation, but little is known about the greater developmental networks within which these genes interact to pattern a wing. Here we took a large-scale transcriptomic approach to identify the network of genes involved in Heliconius wing pattern development and variation. This included applying over 140 transcriptome microarrays to assay gene expression in dissected wing pattern elements across a range of developmental stages and wing pattern morphs of Heliconius erato.

Results

We identified a number of putative early prepattern genes with color-pattern related expression domains. We also identified 51 genes differentially expressed in association with natural color pattern variation. Of these, the previously identified color pattern “switch gene” optix was recovered as the first transcript to show color-specific differential expression. Most differentially expressed genes were transcribed late in pupal development and have roles in cuticle formation or pigment synthesis. These include previously undescribed transporter genes associated with ommochrome pigmentation. Furthermore, we observed upregulation of melanin-repressing genes such as ebony and Dat1 in non-melanic patterns.

Conclusions

This study identifies many new genes implicated in butterfly wing pattern development and provides a glimpse into the number and types of genes affected by variation in genes that drive color pattern evolution.

Comprehensive microarray-based analysis for stage-specific larval camouflage pattern-associated genes in the swallowtail butterfly, Papilio xuthus

Background

Body coloration is an ecologically important trait that is often involved in prey-predator interactions through mimicry and crypsis. Although this subject has attracted the interest of biologists and the general public, our scientific knowledge on the subject remains fragmentary. In the caterpillar of the swallowtail butterfly Papilio xuthus, spectacular changes in the color pattern are observed; the insect mimics bird droppings (mimetic pattern) as a young larva, and switches to a green camouflage coloration (cryptic pattern) in the final instar. Despite the wide variety and significance of larval color patterns, few studies have been conducted at a molecular level compared with the number of studies on adult butterfly wing patterns.

Results

To obtain a catalog of genes involved in larval mimetic and cryptic pattern formation, we constructed expressed sequence tag (EST) libraries of larval epidermis for P. xuthus, and P. polytes that contained 20,736 and 5,376 clones, respectively, representing one of the largest collections available in butterflies. A comparison with silkworm epidermal EST information revealed the high expression of putative blue and yellow pigment-binding proteins in Papilio species. We also designed a microarray from the EST dataset information, analyzed more than five stages each for six markings, and confirmed spatial expression patterns by whole-mount in situ hybridization. Hence, we succeeded in elucidating many novel marking-specific genes for mimetic and cryptic pattern formation, including pigment-binding protein genes, the melanin-associated gene yellow-h3, the ecdysteroid synthesis enzyme gene 3-dehydroecdysone 3b-reductase, and Papilio-specific genes. We also found many cuticular protein genes with marking specificity that may be associated with the unique surface nanostructure of the markings. Furthermore, we identified two transcription factors, spalt and ecdysteroid signal-related E75, as genes expressed in larval eyespot markings. This finding suggests that E75 is a strong candidate mediator of the hormone-dependent coordination of larval pattern formation.

Conclusions

This study is one of the most comprehensive molecular analyses of complicated morphological features, and it will serve as a new resource for studying insect mimetic and cryptic pattern formation in general. The wide variety of marking-associated genes (both regulatory and structural genes) identified by our screening indicates that a similar strategy will be effective for understanding other complex traits.

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.

Rapid Evolution of Assortative Fertilization between Recently Allopatric Species of Drosophila

The virilis group of Drosophila represents a relatively unexplored but potentially useful model to investigate the genetics of speciation. Good resolution of phylogenetic relationships and the ability to obtain fertile hybrid offspring make the group especially promising for analysis of genetic changes underlying reproductive isolation separate from hybrid sterility and inviability. Phylogenetic analyses reveal a close relationship between the sister species, Drosophila americana and D. novamexicana, yet excepting their contemporary allopatric distributions, factors that contribute to reproductive isolation between this species pair remain uncharacterized. A previous report has shown reduced progeny numbers in laboratory crosses between the two species, especially when female D. novamexicana are crossed with male D. americana. We show that the hatch rate of eggs produced from heterospecific matings is reduced relative to conspecific matings. Failure of eggs to hatch, and consequent reduction in hybrid progeny number, is caused by low fertilization success of heterospecific sperm, thus representing a postmating, prezygotic incompatibility. Following insemination, storage and motility of heterospecific sperm is visibly compromised in female D. novamexicana. Our results provide evidence for a mechanism of reproductive isolation that is seldom reported for Drosophila species, and indicate the rapid evolution of postmating, prezygotic reproductive barriers in allopatry.

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.

A conserved developmental patterning network produces quantitatively different output in multiple species of Drosophila

Differences in the level, timing, or location of gene expression can contribute to alternative phenotypes at the molecular and organismal level. Understanding the origins of expression differences is complicated by the fact that organismal morphology and gene regulatory networks could potentially vary even between closely related species. To assess the scope of such changes, we used high-resolution imaging methods to measure mRNA expression in blastoderm embryos of Drosophila yakuba and Drosophila pseudoobscura and assembled these data into cellular resolution atlases, where expression levels for 13 genes in the segmentation network are averaged into species-specific, cellular resolution morphological frameworks. We demonstrate that the blastoderm embryos of these species differ in their morphology in terms of size, shape, and number of nuclei. We present an approach to compare cellular gene expression patterns between species, while accounting for varying embryo morphology, and apply it to our data and an equivalent dataset for Drosophila melanogaster. Our analysis reveals that all individual genes differ quantitatively in their spatio-temporal expression patterns between these species, primarily in terms of their relative position and dynamics. Despite many small quantitative differences, cellular gene expression profiles for the whole set of genes examined are largely similar. This suggests that cell types at this stage of development are conserved, though they can differ in their relative position by up to 3–4 cell widths and in their relative proportion between species by as much as 5-fold. Quantitative differences in the dynamics and relative level of a subset of genes between corresponding cell types may reflect altered regulatory functions between species. Our results emphasize that transcriptional networks can diverge over short evolutionary timescales and that even small changes can lead to distinct output in terms of the placement and number of equivalent cells.

For a gene to function properly, it must be active in the right place, at the right time, and in the right amount. Changes in any of these features can lead to observable differences between individuals and species and in some cases can lead to disease. We do not currently understand how the position, timing, and amount of gene expression is encoded in DNA sequence. One approach to this problem is to compare how gene expression differs between species and to try to relate changes in DNA sequence to changes in gene expression. Here, we take the first step by comparing gene expression patterns at high spatial and temporal resolution between embryos of three species of fruit flies. We develop methods for comparing gene expression in individual cells, which allow us to control for variation in the size, shape, and number of nuclei between embryos. We find measurable quantitative differences in the patterns for all individual genes that we have examined. However, by considering all genes in our dataset at once, we show that many genes are changing together, leading to largely equivalent types of cells in these three species.

Multiple Quantitative Trait Loci Influence the Shape of a Male-Specific Genital Structure in Drosophila melanogaster

The observation that male genitalia diverge more rapidly than other morphological traits during evolution is taxonomically widespread and likely due to some form of sexual selection. One way to elucidate the evolutionary forces acting on these traits is to detail the genetic architecture of variation both within and between species, a program of research that is considerably more tractable in a model system. Drosophila melanogaster and its sibling species, D. simulans, D. mauritiana, and D. sechellia, are morphologically distinguishable only by the shape of the posterior lobe, a male-specific elaboration of the genital arch. We extend earlier studies identifying quantitative trait loci (QTL) responsible for lobe divergence across species and report the first genetic dissection of lobe shape variation within a species. Using an advanced intercross mapping design, we identify three autosomal QTL contributing to the difference in lobe shape between a pair of D. melanogaster inbred lines. The QTL each contribute 4.6–10.7% to shape variation, and two show a significant epistatic interaction. Interestingly, these intraspecific QTL map to the same locations as interspecific lobe QTL, implying some shared genetic control of the trait within and between species. As a first step toward a mechanistic understanding of natural lobe shape variation, we find an association between our QTL data and a set of genes that show sex-biased expression in the developing genital imaginal disc (the precursor of the adult genitalia). These genes are good candidates to harbor naturally segregating polymorphisms contributing to posterior lobe shape.

Divergent enhancer haplotype of ebony on inversion In(3R)Payne associated with pigmentation variation in a tropical population of Drosophila melanogaster

The pattern and intensity of pigmentation have direct impact on individual fitness through various ecological factors. In a Drosophila melanogaster population from southern Japan, thoracic trident pigmentation intensity of most of the strains could be classified into Dark or Light-type. The expression level variation of the ebony gene correlated well with this phenotype and the allelic differences in expression indicated that the variation is partly due to cis-regulatory changes. In the ∼13 kb gene region, we identified 17 nucleotide sites and 2 indels that were in complete association with the thoracic trident pigmentation intensity. Interestingly, 11 out of 19 sites located within ∼0.5 kb of the core epidermis enhancer. These sites had no obvious association with the abdominal pigmentation intensity in the previously analysed African populations from Uganda and Kenya, which suggested that multiple potential mutational pathways in the cis-regulatory control region of a single gene could lead to similar phenotypic variation within this species. We also found that the Light-type enhancer haplotype is strongly linked to a cosmopolitan inversion, In(3R)Payne, which is predominant in warmer climatic regions in both hemispheres. The sequence pattern suggested that the strong linkage may be due to selective forces related to thermal adaptation. The inferred selection for lighter pigmentation in the Japanese population is in the opposite direction of the previously reported case of selection for darker individuals in African populations. Nevertheless, both adaptive changes involved cis-regulatory changes of ebony, which shows that this gene is likely to be a common target of natural selection.

The genetic basis of rapidly evolving male genital morphology in Drosophila

The external genitalia are some of the most rapidly evolving morphological structures in insects. The posterior lobe of the male genital arch shows striking differences in both size and shape among closely related species of the Drosophila melanogaster species subgroup. Here, we dissect the genetic basis of posterior lobe morphology between D. mauritiana and D. sechellia, two island endemic species that last shared a common ancestor ∼300,000 years ago. We test a large collection of genome-wide homozygous D. mauritiana genetic introgressions, which collectively cover ∼50% of the genome, for their morphological effects when placed in a D. sechellia genetic background. We find several introgressions that have large effects on posterior lobe morphology and that posterior lobe size and posterior lobe shape can be separated genetically for some of the loci that specify morphology. Using next generation sequencing technology, we perform whole transcriptome gene expression analyses of the larval genital imaginal disc of D. mauritiana, D. sechellia, and two D. mauritiana-D. sechellia hybrid introgression genotypes that each have large effects on either posterior lobe size or posterior lobe shape. Many of the genes we identify as differentially expressed are expressed at levels similar to D. mauritiana in one introgression hybrid, but are expressed at levels similar to D. sechellia in the other introgression hybrid. However, we also find that both introgression hybrids express some of the same genes at levels similar to D. mauritiana, and notably, that both introgression hybrids possess genes in the insulin receptor signaling pathway, which are expressed at D. mauritiana expression levels. These results suggest the possibility that the insulin signaling pathway might integrate size and shape genetic inputs to establish differences in overall posterior lobe morphology between D. mauritiana and D. sechellia.

The molecular genetics of clinal variation: a case study of ebony and thoracic trident pigmentation in Drosophila melanogaster from eastern Australia

Widespread pigmentation diversity coupled with a well-defined genetic system of melanin synthesis and patterning in Drosophila provides an excellent opportunity to study phenotypes undergoing evolutionary change. Pigmentation variation is highly correlated with different ecological variables and is thought to reflect adaptations to different environments. Several studies have linked candidate genes from Drosophila melanogaster to intra-population variation and interspecific morphological divergence, but less clearly to variation among populations forming pigmentation clines. We characterized a new thoracic trident pigmentation cline in D. melanogaster populations from eastern Australia, and applied a candidate gene approach to explain the majority of the geographically structured phenotypic variation. More melanized populations from higher latitudes tended to express less ebony than their tropical counterparts, and an independent artificial selection experiment confirmed this association. By partitioning temperature dependent effects, we showed that the genetic differences underlying clinal patterns for trident variation at 25 °C do not explain the patterns observed at 16 °C. Changes in thoracic trident pigmentation could be a common evolutionary response to climatically mediated environmental pressures. On the Australian east coast most of the changes appear to be associated with regulatory divergence of the ebony gene but this depends on temperature.