Abdominal pigmentation variation in drosophila polymorpha: geographic variation in the trait, and underlying phylogeography

Drosophila polymorpha is a widespread species that exhibits abdominal pigmentation variation throughout its range. To gain insight into this variation we combined phenotypic and genotypic data to test a series of nested hypotheses. First, we tested the null hypothesis that geographic variation in pigmentation is due to neutral factors. We used nested clade analysis to examine the distribution of haplotypes from a nuclear and a mitochrondrial locus. Restricted gene flow via isolation by distance, the primary inference of this phylogeographic analysis, was then used to generate and test the hypothesis of increasing average abdominal pigmentation difference with increasing geographic distance. We found no correlation between geographic distance and phenotypic distance. We then tested the hypothesis that pigmentation is affected by environmental differences among localities. We found a significant effect of habitat type on the average abdominal pigmentation phenotype of different localities. Finally, we tested the hypothesis that pigmentation in D. polymorpha is associated with desiccation resistance. We found that dark individuals of both sexes survived significantly longer in a desiccating environment than light individuals. These patterns combined lead us to hypothesize that abdominal pigmentation variation in D. polymorpha is important in mediating the organism's interactions with local ecological factors.     

Evolution of abdominal pigmentation differences across species in the Drosophila dunni subgroup

The Drosophila dunni subgroup displays a nearly perfect latitudinal cline in abdominal pigmentation that likely resulted from selective forces acting in the habitat of each species during speciation. Here we characterize the nature of this clinal variation by developing a quantitative measure to assess variation in abdominal pigmentation within and between the D. dunni subgroup species. Using discriminant analysis, we confirm the existence of a cline and find that our quantitative measure of pigmentation distinguishes each of the species with singular efficacy. We then combine our quantitative phenotypic analysis of pigmentation with the phylogeny of the D. dunni subgroup species and map the species relationships into the three-dimensional morphological space defined by our pigmentation measures. In this manner, we can visualize how the species have traversed the morphological pigmentation space during the course of speciation. Our analysis reveals that natural selection has caused overall intensity of pigmentation among the northernmost species of the cline to converge. Along with this convergence in phenotype has been a relaxation in expression of sexual dimorphism in these species, indicating a possible shift in the relative intensity of natural and sexual selection. Our analysis indicates an accelerated rate of change in pigmentation for the darkest species in addition to this species evolving a novel abdominal pigmentation trait.     

Contrasting modes of natural selection acting on pigmentation genes in the Drosophila dunni subgroup

Genes that encode for divergent adaptive traits may have genealogies that contrast with those from loci that are not functionally involved in differentiation. Here, we examine DNA sequence variation among the species of the eastern Caribbean Drosophila dunni subgroup at two loci, yellow and dopa decaboxylase (Ddc), which both play integral roles in pigmentation patterning of adult Drosophila. Phylogenetic analyses of these loci produce gene genealogies with topologies that mirror those described for other nuclear genes: the six morphologically distinct species within the subgroup are divided into only three lineages, with one lineage containing four species that share extensive ancestral polymorphism. At the Ddc locus these major lineages are delineated only by silent site variation. We observe a significantly higher rate of synonymous site divergence than non-synonymous divergence, consistent with strong purifying selection acting on the locus. In contrast, the yellow locus exhibits patterns of amino acid divergence and nucleotide diversity that are consistent with recent diversifying selection acting in two different lineages. This selection appears to be targeting amino acid variants in the signal sequence of the Yellow protein, a region which is tightly constrained among members of the larger D. cardini radiation. This result highlights not only the potential importance of yellow in the evolution of divergent pigmentation patterns among members of the D. dunni subgroup, but also hints that variation in signal peptide sequences may play a role in phenotypic diversification.     

Ecological genetics of abdominal pigmentation in Drosophila falleni: a pleiotropic link to nematode parasitism

Drosophila falleni belongs to the quinaria species group, whose species vary considerably in patterns of wing and abdominal pigmentation. Drosophila falleni itself exhibits substantial variation among wild flies in abdominal spotting patterns. A selection experiment revealed that natural populations of D. falleni harbor high levels of genetic variation for spot number: in 10 generations of selection modal spot number within populations declined from 18 (the modal number in wild-caught females) to as low as zero. Rearing flies at different temperatures shows that some of the variation among wild flies is likely to reflect variation in the environmental conditions under which they developed. Fitness assays did not reveal any cost of reduced spot number with respect to development time, adult survival, or female fecundity. However, spotless flies were almost twice as susceptible to infection by the nematode parasite Howardula aoronymphium. Thus, selection exerted by nematode parasites may influence pigmentation patterns and other, genetically correlated traits in natural populations D. falleni.     

Molecular variation at the In(2L)t proximal breakpoint site in natural populations of Drosophila melanogaster and D. simulans

A previous study of nucleotide polymorphism in a Costa Rican population of Drosophila melanogaster found evidence for a nonneutral deficiency in the number of haplotypes near the proximal breakpoint of In(2L)t, a common inversion polymorphism in this species. Another striking feature of the data was a window of unusually high nucleotide diversity spanning the breakpoint site. To distinguish between selective and neutral demographic explanations for the observed patterns in the data, we sample alleles from three additional populations of D. melanogaster and one population of D. simulans. We find that the strength of associations among sites found at the breakpoint varies between populations of D. melanogaster. In D. simulans, analysis of the homologous region reveals unusually elevated levels of nucleotide polymorphism spanning the breakpoint site. As with American populations of D. melanogaster, our D. simulans sample shows a marked reduction in the number of haplotypes but not in nucleotide diversity. Haplotype tests reveal a significant deficiency in the number of haplotypes relative to the neutral expectation in the D. simulans sample and some populations of D. melanogaster. At the breakpoint site, the level of divergence between haplotype classes is comparable to interspecific divergence. The observation of interspecific polymorphisms that differentiate major haplotype classes in both species suggests that haplotype classes at this locus are considerably old. When considered in the context of other studies on patterns of variation within and between populations of D. melanogaster and D. simulans, our data appear more consistent with the operation of selection than with simple demographic explanations.     

Unusual pattern of single nucleotide polymorphism at the exuperantia2 locus of Drosophila pseudoobscura

We have investigated the pattern of DNA sequence variation at the exuperantia2 locus in Drosophila pseudoobscura. This adds to the increasing dataset of genetic variation in D. pseudoobscura, a useful model species for evolutionary genetic studies. The level of silent site nucleotide diversity and the divergence from an outgroup Drosophila miranda are comparable with those for other X-linked loci. One peculiar pattern at the exu2 locus of D. pseudoobscura is a complete linkage disequilibrium between two SNPs, one of which is a replacement site. As a result, there are two distinct haplotype groups in our dataset. Based upon the comparisons with the outgroup sequences from D. miranda and Drosophila persimilis, we show that the newly derived haplotype group has lower diversity than the ancestral haplotype group. The pattern of protein evolution at exu2 shows some deviation from the neutral model. Together, these and other characteristics of the exu2 locus suggest the action of selection on the pattern of SNP variation, consistent with a partial selective sweep associated with the newly derived haplotype.     

Molecular Evolution between Drosophila Melanogaster and D. Simulans: Reduced Codon Bias, Faster Rates of Amino Acid Substitution, and Larger Proteins in D. Melanogaster

Both natural selection and mutational biases contribute to variation in codon usage bias within Drosophila species. This study addresses the cause of codon bias differences between the sibling species, Drosophila melanogaster and D. simulans. Under a model of mutation-selection-drift, variation in mutational processes between species predicts greater base composition differences in neutrally evolving regions than in highly biased genes. Variation in selection intensity, however, predicts larger base composition differences in highly biased loci. Greater differences in the G+C content of 34 coding regions than 46 intron sequences between D. melanogaster and D. simulans suggest that D. melanogaster has undergone a reduction in selection intensity for codon bias. Computer simulations suggest at least a fivefold reduction in N(e)s at silent sites in this lineage. Other classes of molecular change show lineage effects between these species. Rates of amino acid substitution are higher in the D. melanogaster lineage than in D. simulans in 14 genes for which outgroup sequences are available. Surprisingly, protein sizes are larger in D. melanogaster than in D. simulans in the 34 genes compared between the two species. A substantial fraction of silent, replacement, and insertion/deletion mutations in coding regions may be weakly selected in Drosophila.     

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.     

Phylogenetic analysis of the Cardini group of Drosophila with respect to changes in pigmentation

Phenotypic variability is the engine that drives future diversification with the expectation that polymorphic ancestors give rise to descendants harboring a subset of the ancestral variation. Here we examine evolutionary transitions from polymorphism to monomorphism in a visually striking New World radiation of fruit flies, the Drosophila cardini group. This group is distributed across the Americas and the Caribbean islands and exhibits a wide spectrum of abdominal pigmentation variation. Specifically, the D. dunni subgroup consists of Caribbean island endemics, each of which is monomorphic for its pigmentation pattern, with an interspecific cline of pigmentation across the islands. The D. cardini subgroup consists of American continental species with wide-ranging distributions and intraspecifically variable abdominal pigmentation. We determined the phylogeny of 18 species and subspecies using three nuclear and three mitochondrial regions analyzed with maximum parsimony, maximum likelihood, and Bayesian methods. The topology produced from a combined dataset exhibited high support values at all nodes, and differed from earlier phylogenetic hypotheses based on polytene chromosome inversion patterns and isozyme data. We find that the D. dunni subgroup species, with the exception of D. belladunni, are derived from a single source not of direct South American origin and their dispersal across the islands of the Caribbean does not follow a simple stepping-stone model. Morphological changes in pigmentation across the island species are incongruent with the colonization history of the group indicating that natural selection may have played a role in the determination of this character. Finally, we demonstrate that monomorphic species have arisen independently from polymorphic ancestors two to three times.     

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.     

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.     

Clines and adaptive evolution in the methuselah gene region in Drosophila melanogaster

In an effort to characterize further the patterns of selection and adaptive evolution at the methuselah locus in Drosophila species, we extended an analysis of geographical variation to include single nucleotide polymorphisms (SNPs) in adjacent genes on either side of the mth locus, and examined the molecular variation in a neighbouring methuselah paralogue (mth2). An analysis of 13 SNPs spanning a region of nearly 19 kilobases surrounding the mth locus demonstrated that a clinal pattern associated with the most common mth haplotype does not extend to adjacent gene loci, providing compelling evidence that the clinal pattern results from selection on as yet unidentified sites associated with the functional mth locus. mth2 exhibited a significant pattern of adaptive divergence among D. melanogaster, D. simulans and D. yakuba similar to that seen at mth. However, Ka : Ks ratios indicate a difference in levels of functional constraint at the two methuselah, loci with mth2 exhibiting a five- to six-fold reduction in levels of amino acid divergence relative to mth.     

Genetic basis of sex-specific color pattern variation in Drosophila malerkotliana

Pigmentation is a rapidly evolving trait that can play important roles in mimicry, sexual selection, thermoregulation, and other adaptive processes in many groups of animals. In Drosophila, pigmentation can differ dramatically among closely related taxa, presenting a good opportunity to dissect the genetic changes underlying species divergence. In this report, we investigate the genetic basis of color pattern variation between two allopatric subspecies of Drosophila malerkotliana, a widespread member of the ananassae species subgroup. In D. malerkotliana malerkotliana, the last three abdominal segments are darkly pigmented in males but not in females, while in D. malerkotliana pallens both sexes lack dark pigmentation. Composite interval mapping in F(2) hybrid progeny shows that this difference is largely controlled by three quantitative trait loci (QTL) located on the 2L chromosome arm, which is homologous to the 3R of D. melanogaster (Muller element E). Using highly recombinant introgression strains produced by repeated backcrossing and phenotypic selection, we show that these QTL do not correspond to any of the candidate genes known to be involved in pigment patterning and synthesis in Drosophila. These results, in combination with similar analyses in other Drosophila species, indicate that different genetic and molecular changes are responsible for the evolution of similar phenotypic traits in different lineages. This feature makes Drosophila color patterns a powerful model for investigating how the genetic basis of trait evolution is influenced by the intrinsic organization of regulatory pathways controlling the development of these traits.     

Alu insertion polymorphisms for the study of human genomic diversity

We analyze genetic variation at fused1, a locus that is close to the centromere of the X chromosome-autosome (X/4) fusion in Drosophila americana. In contrast to other X-linked and autosomal genes, for which a lack of population subdivision in D. americana has been observed at the DNA level, we find strong haplotype structure associated with the alternative chromosomal arrangements. There are several derived fixed differences at fused1 (including one amino acid replacement) between two haplotype classes of this locus. From these results, we obtain an estimate of an age of approximately 0.61 million years for the origin of the two haplotypes of the fused1 gene. Haplotypes associated with the X/4 fusion have less DNA sequence variation at fused1 than haplotypes associated with the ancestral chromosome arrangement. The X/4 haplotypes also exhibit clinal variation for the allele frequencies of the three most common amino acid replacement polymorphisms, but not for adjacent silent polymorphisms. These patterns of variation are best explained as a result of selection acting on amino acid substitutions, with geographic variation in selection pressures.     

Drosophila duplication hotspots are associated with late-replicating regions of the genome

Duplications play a significant role in both extremes of the phenotypic spectrum of newly arising mutations: they can have severe deleterious effects (e.g. duplications underlie a variety of diseases) but can also be highly advantageous. The phenotypic potential of newly arisen duplications has stimulated wide interest in both the mutational and selective processes shaping these variants in the genome. Here we take advantage of the Drosophila simulans-Drosophila melanogaster genetic system to further our understanding of both processes. Regarding mutational processes, the study of two closely related species allows investigation of the potential existence of shared duplication hotspots, and the similarities and differences between the two genomes can be used to dissect its underlying causes. Regarding selection, the difference in the effective population size between the two species can be leveraged to ask questions about the strength of selection acting on different classes of duplications. In this study, we conducted a survey of duplication polymorphisms in 14 different lines of D. simulans using tiling microarrays and combined it with an analogous survey for the D. melanogaster genome. By integrating the two datasets, we identified duplication hotspots conserved between the two species. However, unlike the duplication hotspots identified in mammalian genomes, Drosophila duplication hotspots are not associated with sequences of high sequence identity capable of mediating non-allelic homologous recombination. Instead, Drosophila duplication hotspots are associated with late-replicating regions of the genome, suggesting a link between DNA replication and duplication rates. We also found evidence supporting a higher effectiveness of selection on duplications in D. simulans than in D. melanogaster. This is also true for duplications segregating at high frequency, where we find evidence in D. simulans that a sizeable fraction of these mutations is being driven to fixation by positive selection.     

Patterns of mutation and selection at synonymous sites in Drosophila

That natural selection affects molecular evolution at synonymous sites in protein-coding sequences is well established and is thought to predominantly reflect selection for translational efficiency/accuracy mediated through codon bias. However, a recently developed maximum likelihood framework, when applied to 18 coding sequences in 3 species of Drosophila, confirmed an earlier report that the Notch gene in Drosophila melanogaster was evolving under selection in favor of those codons defined as unpreferred in this species. This finding opened the possibility that synonymous sites may be subject to a variety of selective pressures beyond weak selection for increased frequencies of the codons currently defined as "preferred" in D. melanogaster. To further explore patterns of synonymous site evolution in Drosophila in a lineage-specific manner, we expanded the application of the maximum likelihood framework to 8,452 protein coding sequences with well-defined orthology in D. melanogaster, Drosophila sechellia, and Drosophila yakuba. Our analyses reveal intragenomic and interspecific variation in mutational patterns as well as in patterns and intensity of selection on synonymous sites. In D. melanogaster, our results provide little statistical evidence for recent selection on synonymous sites, and Notch remains an outlier. In contrast, in D. sechellia our findings provide evidence in support of selection predominantly in favor of preferred codons. However, there is a small subset of genes in this species that appear to be evolving under selection in favor of unpreferred codons, which indicates that selection on synonymous sites is not limited to the preferential fixation of mutations that enhance the speed or accuracy of translation in this species.     

Genetics of a difference in pigmentation between Drosophila yakuba and Drosophila santomea

Abstract.— Drosophila yakuba is a species widespread in Africa, whereas D. santomea, its newly discovered sister species, is endemic to the volcanic island of São Tomé in the Gulf of Guinea. Drosophila santomea probably formed after colonization of the island by its common ancestor with D. yakuba. The two species differ strikingly in pigmentation: D. santomea, unlike the other eight species in the D. melanogaster subgroup, almost completely lacks dark abdominal pigmentation. D. yakuba shows the sexually dimorphic pigmentation typical of the group: both sexes have melanic patterns on the abdomen, but males are much darker than females. A genetic analysis of this species difference using morphological markers shows that the X chromosome accounts for nearly 90% of the species difference in the area of abdomen that is pigmented and that at least three genes (one on each major chromosome) are involved in each sex. The order of chromosome effects on pigmentation area are the same in males and females, suggesting that loss of pigmentation in D. santomea may have involved the same genes in both sexes. Further genetic analysis of the interspecific difference between males in pigmentation area and intensity using molecular markers shows that at least five genes are responsible, with no single locus having an overwhelming effect on the trait. The species difference is thus oligogenic or polygenic. Different chromosomal regions from each of the two species influenced pigmentation in the same direction, suggesting that the species difference (at least in males) is due to natural or sexual selection and not genetic drift. Measurements of sexual isolation between the species in both light and dark conditions show no difference, suggesting that the pigmentation difference is not an important cue for interspecific mate discrimination. Using DNA sequence differences in nine noncoding regions, we estimate that D. santomea and D. yakuba diverged about 400,000 years ago, a time similar to the divergences between two other well‐studied pair of species in the subgroup, both of which also involved island colonization.     

Hidden Electrophoretic Variation at the Xanthine Dehydrogenase Locus in a Natural Population of DROSOPHILA MELANOGASTER

Sixty-two isochromosomal lines of D. melanogaster were screened for cryptic electrophoretic variation at the xanthine dehydrogenase (XDH) locus. Sequential polyacrylamide vertical slab gel electrophoresis was performed using four electrophoretic criteria. A total of 15 classes of electromorphs were revealed. D. melanogaster appears to exhibit as much polymorphism at this locus as other extensively studied Drosophila species.--No evidence for loci on the X or second chromosomes which modified XDH mobility was found. Six of the electromorphs were mapped to the Xdh (ry) structural locus. Eight of the remaining nine classes exhibited mobility variation consistent with structural variation at the Xdh locus. The final class exhibited aberrant patterns and is under further study.     

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.     

Genetic and developmental analysis of abdominal pigmentation differences across species in the Drosophila dunni subgroup

Abdominal pigmentation pattern varies dramatically among the species of the Drosophila dunni subgroup across the islands of the Caribbean. Previously, we developed a quantitative measure of abdominal pigmentation to assess phenotypic variation within and between species of this group. In this paper, we use this quantitative measure in an interspecific genetic analysis to decipher the underlying genetic basis of pigmentation differences between one of the lightest and the darkest species in the group. Our analysis shows that pigmentation expression in different areas of the abdomen is under separate genetic control. For these different abdominal regions, we detected a wide range of genetic effects, including X-linked, autosomal additive, near single-gene dominant, and sex-specific effects. Combining these genetic results with our earlier phenotypic and phylogenetic analyses, we present a simple conceptual model to explain how change in the control of expression of pigmentation has evolved throughout the D. dunni subgroup.