Pleiotropic effects of regulatory variation in tan result in correlation of two pigmentation traits in Drosophila melanogaster

Traits with a common genetic basis frequently display correlated phenotypic responses to selection or environmental conditions. In Drosophila melanogaster, pigmentation of the abdomen and a trident‐shaped region on the thorax are genetically correlated. Here, we used a pooled replicated genomewide association approach (Pool‐GWAS) to identify the genetic basis of variation in thoracic trident pigmentation in two Drosophila melanogaster populations. We confirmed the previously reported large effect of ebony and the association of the cosmopolitan inversion In(3R)Payne. For the first time, we identified tan as another major locus contributing to variation in trident pigmentation. Intriguingly, the regulatory variants of tan that were most strongly associated with female abdominal pigmentation also showed a strong association with trident pigmentation. We validated this common genetic basis in transgenic assays and found qualitatively similar effects on trident and abdominal pigmentation. Further work is required to determine whether this genetic correlation is favoured by natural selection or reflects a neutral by‐product of a shared regulatory architecture.     

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.     

Selection for abdominal tergite pigmentation and correlated responses in the trident: a case study in Drosophila melanogaster

In Drosophila melanogaster, abdominal tergite pigmentation and the appearance of a trident‐shaped thoracic pattern exhibit similar biogeographical variation and sensitivity to temperature. These pigmentation traits may be under common selection pressure in natural populations or may be genetically correlated. To investigate the nature of this interaction, replicated populations of D. melanogaster were selected for increased or decreased melanization of the abdominal tergites for 40 generations. Selection for abdominal tergite pigmentation leads to correlated changes in trident formation. Although selection was performed only on female flies, male pigmentation also responded to selection. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 106 , 287–294.     

Complex patterns of cis‐regulatory polymorphisms in ebony underlie standing pigmentation variation in Drosophila melanogaster

Pigmentation traits in adult Drosophila melanogaster were used in this study to investigate how phenotypic variations in continuous ecological traits can be maintained in a natural population. First, pigmentation variation in the adult female was measured at seven different body positions in 20 strains from the Drosophila melanogaster Genetic Reference Panel (DGRP) originating from a natural population in North Carolina. Next, to assess the contributions of cis‐regulatory polymorphisms of the genes involved in the melanin biosynthesis pathway, allele‐specific expression levels of four genes were quantified by amplicon sequencing using a 454 GS Junior. Among those genes, ebony was significantly associated with pigmentation intensity of the thoracic segment. Detailed sequence analysis of the gene regulatory regions of this gene indicated that many different functional cis‐regulatory alleles are segregating in the population and that variations outside the core enhancer element could potentially play important roles in the regulation of gene expression. In addition, a slight enrichment of distantly associated SNP pairs was observed in the ~10 kb cis‐regulatory region of ebony, which suggested the presence of interacting elements scattered across the region. In contrast, sequence analysis in the core cis‐regulatory region of tan indicated that SNPs within the region are significantly associated with allele‐specific expression level of this gene. Collectively, the data suggest that the underlying genetic differences in the cis‐regulatory regions that control intraspecific pigmentation variation can be more complex than those of interspecific pigmentation trait differences, where causal genetic changes are typically confined to modular enhancer elements.     

Annotating ebony on the fly

The distinctive black phenotype of ebony mutants has made it one of the most widely used phenotypic markers in Drosophila genetics. Without doubt, ebony showcases the fruits of the fly community's labours to annotate gene function. As of this writing, FlyBase lists 142 references, 1277 fly stocks, 15 phenotypes and 44 alleles. In addition to its namesake pigmentation phenotype, ebony mutants affect other traits, including phototaxis and courtship. With phenotypic consequences of ebony variants readily apparent in the laboratory, does natural selection also see them in the wild? In this issue of Molecular Ecology, Pool & Aquadro investigate this question and found signs of natural selection on the ebony gene that appear to have resulted from selection for darker pigmentation at higher elevations in sub‐Saharan populations of Drosophila melanogaster. Such findings from population genomic analysis of wild‐derived strains should be included in gene annotations to provide a more holistic view of a gene's function. The evolutionary annotation of ebony added by Pool & Aquadro substantiates that pigmentation can be adaptive and implicates elevation as an important selective factor. This is important progress because the selective factors seem to differ between populations and species. In addition, the study raises issues to consider when extrapolating from selection at the molecular level to selection at the phenotypic level.     

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.     

Phenotypic variability of thoracic pigmentation in Indian populations of Drosophila melanogaster

Fourteen natural populations of Drosophila melanogaster collected along diverse latitudinal and altitudinal ranges, were analysed for thoracic trident pigmentation at two different breeding temperatures (17 and 25°C) for both sexes. Statistical analyses showed significant intra- as well as interpopulational genetic variations. Clinal patterns for thoracic trident pigmentation along latitude as well as altitude are significantly correlated with thermal amplitude ( T _cv). Crosses between dark (Shimla) and light (Kalka) variants of trident pigmentation produced intermediate phenotypes, and a clear maternal effect was observed in the reciprocal F ₁ crosses. Data suggest that the temperature differences have resulted in phenotypic plasticity of thoracic pigmentation across the Indian subcontinent.     

Genetic analysis of differentiation among breeding ponds reveals a candidate gene for local adaptation in Rana arvalis

One of the main questions in evolutionary and conservation biology is how geographical and environmental features of the landscape shape neutral and adaptive genetic variation in natural populations. The identification of genomic polymorphisms that account for adaptive variation can aid in finding candidate loci for local adaptation. Consequently, a comparison of spatial patterns in neutral markers and loci under selection may help disentangle the effects of gene flow, genetic drift and selection at the landscape scale. Many amphibians breed in wetlands, which differ in environmental conditions and in the degree of isolation, enhancing the potential for local adaptation. We used microsatellite markers to measure genetic differentiation among 17 local populations of Rana arvalis breeding in a network of wetlands. We found that locus RC08604 deviated from neutral expectations, suggesting that it is a good candidate for directional selection. We used a genetic network analysis to show that the allele distribution in this locus is correlated with habitat characteristics, whereas this was not the case at neutral markers that displayed a different allele distribution and population network in the study area. The graph approach illustrated the genomic heterogeneity (neutral loci vs. the candidate locus for directional selection) of gene exchange and genetic divergence among populations under directional selection. Limited gene flow between wetlands was only observed at the candidate genomic region under directional selection. RC08604 is partially located inside an up‐regulated thyroid‐hormone receptor (TRβ) gene coordinating the expression of other genes during metamorphosis and appears to be linked with variation in larval life‐history traits found among R. arvalis populations. We suggest that directional selection on genes coding larval life‐history traits is strong enough to maintain the divergence in these genomic regions, reducing the effective recombination of locally adapted alleles but not in other regions of the genome. Integrating this knowledge into conservation plans at the landscape scale will improve the design of management strategies to preserve adaptive genetic diversity in wetland networks.     

Genetic Architecture of Abdominal Pigmentation in Drosophila melanogaster

Pigmentation varies within and between species and is often adaptive. The amount of pigmentation on the abdomen of Drosophila melanogaster is a relatively simple morphological trait, which serves as a model for mapping the genetic basis of variation in complex phenotypes. Here, we assessed natural variation in female abdominal pigmentation in 175 sequenced inbred lines of the Drosophila melanogaster Genetic Reference Panel, derived from the Raleigh, NC population. We quantified the proportion of melanization on the two most posterior abdominal segments, tergites 5 and 6 (T5, T6). We found significant genetic variation in the proportion of melanization and high broad-sense heritabilities for each tergite. Genome-wide association studies identified over 150 DNA variants associated with the proportion of melanization on T5 (84), T6 (34), and the difference between T5 and T6 (35). Several of the top variants associated with variation in pigmentation are in tan, ebony, and bric-a-brac1, genes known to affect D. melanogaster abdominal pigmentation. Mutational analyses and targeted RNAi-knockdown showed that 17 out of 28 (61%) novel candidate genes implicated by the genome-wide association study affected abdominal pigmentation. Several of these genes are involved in developmental and regulatory pathways, chitin production, cuticle structure, and vesicle formation and transport. These findings show that genetic variation may affect multiple steps in pathways involved in tergite development and melanization. Variation in these novel candidates may serve as targets for adaptive evolution and sexual selection in D. melanogaster.     

Morphometric differentiation across House Sparrow Passer domesticus populations in Finland in comparison with the neutral expectation for divergence

To understand the biology of organisms it is important to take into account the evolutionary forces that have acted on their constituent populations. Neutral genetic variation is often assumed to reflect variation in quantitative traits under selection, though with even low neutral divergence there can be substantial differentiation in quantitative genetic variation associated with locally adapted phenotypes. To study the relative roles of natural selection and genetic drift in shaping phenotypic variation, the levels of quantitative divergence based on phenotypes (P_ST) and neutral genetic divergence (F_ST) can be compared. Such a comparison was made between 10 populations of Finnish House Sparrows (n = 238 individuals) collected in 2009 across the whole country. Phenotypic variation in tarsus‐length, wing‐length, bill‐depth, bill‐length and body mass were considered and 13 polymorphic microsatellite loci were analysed to quantify neutral genetic variation. Calculations of P_ST were based on Markov‐Chain Monte Carlo Bayesian estimates of phenotypic variances across and within populations. The robustness of the conclusions of the P_ST–F_ST comparison was evaluated by varying the proportion of variation due to additive genetic effects within and across populations. Our results suggest that body mass is under directional selection, whereas the divergence in other traits does not differ from neutral expectations. These findings suggest candidate traits for considering gene‐based studies of local adaptation. The recognition of locally adapted populations may be of value in the conservation of this declining species.     

Genetic basis of continuous variation in the levels and modular inheritance of pigmentation in cichlid fishes

Variation in pigmentation type and levels is a hallmark of myriad evolutionary radiations, and biologists have long been fascinated by the factors that promote and maintain variation in coloration across populations. Here, we provide insights into the genetic basis of complex and continuous patterns of colour variation in cichlid fishes, which offer a vast diversity of pigmentation patterns that have evolved in response to both natural and sexual selection. Specifically, we crossed two divergent cichlid species to generate an F₂ mapping population that exhibited extensive variation in pigmentation levels and patterns. Our experimental design is robust in that it combines traditional quantitative trait locus (QTL) analysis with population genomics, which has allowed us to move efficiently from QTL interval to candidate gene. In total, we detected 41 QTL and 13 epistatic interactions that underlie melanocyte‐ and xanthophore‐based coloration across the fins and flanks of these fishes. We also identified 2 QTL and 1 interaction for variation in the magnitude of integration among these colour traits. This finding in particular is notable as there are marked differences both within and between species with respect to the complexity of pigmentation patterns. While certain individuals are characterized by more uniform ‘integrated’ colour patterns, others exhibit many more degrees of freedom with respect to the distribution of colour ‘modules’ across the fins and flank. Our data reveal, for the first time, a genetic basis for this difference. Finally, we implicate pax3a as a mediator of continuous variation in the levels of xanthophore‐based colour along the cichlid flank.     

The quantitative genetic basis of clinal divergence in phenotypic plasticity

Phenotypic plasticity is thought to be an important mechanism for adapting to environmental heterogeneity. Nonetheless, the genetic basis of plasticity is still not well understood. In Drosophila melanogaster and D. simulans, body size and thermal stress resistance show clinal patterns along the east coast of Australia, and exhibit plastic responses to different developmental temperatures. The genetic basis of thermal plasticity, and whether the genetic effects underlying clinal variation in traits and their plasticity are similar, remains unknown. Here, we use line‐cross analyses between a tropical and temperate population of Drosophila melanogaster and D. simulans developed at three constant temperatures (18°C, 25°C, and 29°C) to investigate the quantitative genetic basis of clinal divergence in mean thermal response (elevation) and plasticity (slope and curvature) for thermal stress and body size traits. Generally, the genetic effects underlying divergence in mean response and plasticity differed, suggesting that different genetic models may be required to understand the evolution of trait means and plasticity. Furthermore, our results suggest that nonadditive genetic effects, in particular epistasis, may commonly underlie plastic responses, indicating that current models that ignore epistasis may be insufficient to understand and predict evolutionary responses to environmental change.     

NO EFFECT OF ENVIRONMENTAL HETEROGENEITY ON THE MAINTENANCE OF GENETIC VARIATION IN WING SHAPE IN DROSOPHILA MELANOGASTER

Theory suggests that heterogeneous environments should maintain more genetic variation within populations than homogeneous environments, yet experimental evidence for this effect in quantitative traits has been inconsistent. To examine the effect of heterogeneity on quantitative genetic variation, we maintained replicate populations of Drosophila melanogaster under treatments with constant temperatures, temporally variable temperature, or spatially variable temperature with either panmictic or limited migration. Despite observing differences in fitness and divergence in several wing traits between the environments, we did not find any differences in the additive genetic variance for any wing traits among any of the treatments. Although we found an effect of gene flow constraining adaptive divergence between cages in the limited migration treatment, it did not tend to increase within‐population genetic variance relative to any of the other treatments. The lack of any clear and repeatable patterns of response to heterogeneous versus homogeneous environments across several empirical studies suggests that a single general mechanism for the maintenance of standing genetic variation is unlikely; rather, the relative importance of putative mechanisms likely varies considerably from one trait and ecological context to another.     

Müllerian mimicry of a quantitative trait despite contrasting levels of genomic divergence and selection

Hybrid zones, where distinct populations meet and interbreed, give insight into how differences between populations are maintained despite gene flow. Studying clines in genetic loci and adaptive traits across hybrid zones is a powerful method for understanding how selection drives differentiation within a single species, but can also be used to compare parallel divergence in different species responding to a common selective pressure. Here, we study parallel divergence of wing colouration in the butterflies Heliconius erato and H. melpomene, which are distantly related Müllerian mimics which show parallel geographic variation in both discrete variation in pigmentation, and quantitative variation in structural colour. Using geographic cline analysis, we show that clines in these traits are positioned in roughly the same geographic region for both species, which is consistent with direct selection for mimicry. However, the width of the clines varies markedly between species. This difference is explained in part by variation in the strength of selection acting on colour traits within each species, but may also be influenced by differences in the dispersal rate and total strength of selection against hybrids between the species. Genotyping‐by‐sequencing also revealed weaker population structure in H. melpomene, suggesting the hybrid zones may have evolved differently in each species, which may also contribute to the patterns of phenotypic divergence in this system. Overall, we conclude that multiple factors are needed to explain patterns of clinal variation within and between these species, although mimicry has probably played a central role.     

Secondary contact and local adaptation contribute to genome‐wide patterns of clinal variation in Drosophila melanogaster

Populations arrayed along broad latitudinal gradients often show patterns of clinal variation in phenotype and genotype. Such population differentiation can be generated and maintained by both historical demographic events and local adaptation. These evolutionary forces are not mutually exclusive and can in some cases produce nearly identical patterns of genetic differentiation among populations. Here, we investigate the evolutionary forces that generated and maintain clinal variation genome‐wide among populations of Drosophila melanogaster sampled in North America and Australia. We contrast patterns of clinal variation in these continents with patterns of differentiation among ancestral European and African populations. Using established and novel methods we derive here, we show that recently derived North America and Australia populations were likely founded by both European and African lineages and that this hybridization event likely contributed to genome‐wide patterns of parallel clinal variation between continents. The pervasive effects of admixture mean that differentiation at only several hundred loci can be attributed to the operation of spatially varying selection using an F_ST outlier approach. Our results provide novel insight into the well‐studied system of clinal differentiation in D. melanogaster and provide a context for future studies seeking to identify loci contributing to local adaptation in a wide variety of organisms, including other invasive species as well as temperate endemics.     

Adaptive divergence despite strong genetic drift: genomic analysis of the evolutionary mechanisms causing genetic differentiation in the island fox (Urocyon littoralis)

The evolutionary mechanisms generating the tremendous biodiversity of islands have long fascinated evolutionary biologists. Genetic drift and divergent selection are predicted to be strong on islands and both could drive population divergence and speciation. Alternatively, strong genetic drift may preclude adaptation. We conducted a genomic analysis to test the roles of genetic drift and divergent selection in causing genetic differentiation among populations of the island fox (Urocyon littoralis). This species consists of six subspecies, each of which occupies a different California Channel Island. Analysis of 5293 SNP loci generated using Restriction‐site Associated DNA (RAD) sequencing found support for genetic drift as the dominant evolutionary mechanism driving population divergence among island fox populations. In particular, populations had exceptionally low genetic variation, small N_e (range = 2.1–89.7; median = 19.4), and significant genetic signatures of bottlenecks. Moreover, islands with the lowest genetic variation (and, by inference, the strongest historical genetic drift) were most genetically differentiated from mainland grey foxes, and vice versa, indicating genetic drift drives genome‐wide divergence. Nonetheless, outlier tests identified 3.6–6.6% of loci as high F_ST outliers, suggesting that despite strong genetic drift, divergent selection contributes to population divergence. Patterns of similarity among populations based on high F_ST outliers mirrored patterns based on morphology, providing additional evidence that outliers reflect adaptive divergence. Extremely low genetic variation and small N_e in some island fox populations, particularly on San Nicolas Island, suggest that they may be vulnerable to fixation of deleterious alleles, decreased fitness and reduced adaptive potential.     

CRISPR/Cas9‐mediated ebony knockout results in puparium melanism in Spodoptera litura

Insect body pigmentation and coloration are critical to adaption to the environment. To explore the mechanisms that drive pigmentation, we used the clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) genome editing system to target the ebony gene in the non-model insect Spodoptera litura. Ebony is crucial to melanin synthesis in insects. By directly injecting Cas9 messenger RNA and ebony-specific guide RNAs into S. litura embryos, we successfully induced a typical ebony-deficient phenotype of deep coloration of the puparium and induction of melanin formation during the pupal stage. Polymerase chain reaction-based genotype analysis demonstrated that various mutations had occurred at the sites targeted in ebony. Our study clearly demonstrates the function of ebony in the puparium coloration and also provides a potentially useful marker gene for functional studies in S. litura as well as other lepidopteran pests.     

Do different rates of gene flow underlie variation in phenotypic and phenological clines in a montane grasshopper community?

Species responses to environmental change are likely to depend on existing genetic and phenotypic variation, as well as evolutionary potential. A key challenge is to determine whether gene flow might facilitate or impede genomic divergence among populations responding to environmental change, and if emergent phenotypic variation is dependent on gene flow rates. A general expectation is that patterns of genetic differentiation in a set of codistributed species reflect differences in dispersal ability. In less dispersive species, we predict greater genetic divergence and reduced gene flow. This could lead to covariation in life‐history traits due to local adaptation, although plasticity or drift could mirror these patterns. We compare genome‐wide patterns of genetic structure in four phenotypically variable grasshopper species along a steep elevation gradient near Boulder, Colorado, and test the hypothesis that genomic differentiation is greater in short‐winged grasshopper species, and statistically associated with variation in growth, reproductive, and physiological traits along this gradient. In addition, we estimate rates of gene flow under competing demographic models, as well as potential gene flow through surveys of phenological overlap among populations within a species. All species exhibit genetic structure along the elevation gradient and limited gene flow. The most pronounced genetic divergence appears in short‐winged (less dispersive) species, which also exhibit less phenological overlap among populations. A high‐elevation population of the most widespread species, Melanoplus sanguinipes, appears to be a sink population derived from low elevation populations. While dispersal ability has a clear connection to the genetic structure in different species, genetic distance does not predict growth, reproductive, or physiological trait variation in any species, requiring further investigation to clearly link phenotypic divergence to local adaptation.