Recent radiation of endemic Caribbean Drosophila of the dunni subgroup inferred from multilocus DNA sequence variation

Studies of island endemism provide a unique opportunity to elucidate fundamental mechanisms of speciation. Here we examine intra- and interspecific DNA sequence variation at four unlinked genetic loci among populations of the Drosophila dunni subgroup to provide a detailed genealogical portrait of the process of speciation among these island endemic species. Our data indicate two major rounds of diversification that have shaped the D. dunni subgroup. The first occurred 1.6-2.6 million years ago and separated three major lineages, one in Puerto Rico and the Virgin Islands, a second in the northern Lesser Antilles and Barbados, and a third in St. Vincent and Grenada. A second round of diversification occurred in the last 96,000 years in the northern Lesser Antilles and Barbados. The four distinct species that resulted from this recent round of diversification maintain relatively high amounts of genetic variation, similar to that of a closely related mainland species, and share extensive ancestral polymorphism. These data suggest a minimal role for population bottlenecks linked to founder events in the history of the D. dunni subgroup. Further, the recent divergence of these island populations highlights the extremely rapid development of reproductive isolation and distinct patterns of abdominal pigmentation that has occurred in these species.     

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.     

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.     

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.     

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.     

Comparative population genetics of the immunity gene, Relish: is adaptive evolution idiosyncratic?

The frequency of adaptive evolution acting on common loci in distant lineages remains an outstanding question in evolutionary biology. We asked whether the immunity factor, Relish, a gene with a history of directional selection in Drosophila simulans, shows evidence of a similar selective history in other Drosophila species. We found only weak evidence of recurrent adaptive protein evolution at the Relish locus in three sister species pairs, suggesting that this key component of the insect immune system has an idiosyncratic evolutionary history in Drosophila.     

Divergence in coloration and ecological speciation in the Anolis marmoratus species complex

Adaptive divergence in coloration is expected to produce reproductive isolation in species that use colourful signals in mate choice and species recognition. Indeed, many adaptive radiations are characterized by differentiation in colourful signals, suggesting that divergent selection acting on coloration may be an important component of speciation. Populations in the Anolis marmoratus species complex from the Caribbean island of Guadeloupe display striking divergence in the colour and pattern of adult males that occurs over small geographic distances, suggesting strong divergent selection. Here we test the hypothesis that divergence in coloration results in reduced gene flow among populations. We quantify variation in adult male coloration across a habitat gradient between mesic and xeric habitats, use a multilocus coalescent approach to infer historical demographic parameters of divergence, and examine gene flow and population structure using microsatellite variation. We find that colour variation evolved without geographic isolation and in the face of gene flow, consistent with strong divergent selection and that both ecological and sexual selection are implicated. However, we find no significant differentiation at microsatellite loci across populations, suggesting little reproductive isolation and high levels of contemporary gene exchange. Strong divergent selection on loci affecting coloration probably maintains clinal phenotypic variation despite high gene flow at neutral loci, supporting the notion of a porous genome in which adaptive portions of the genome remain fixed whereas neutral portions are homogenized by gene flow and recombination. We discuss the impact of these findings for studies of colour evolution and ecological speciation.     

Antagonistic selection factors induce a continuous population divergence in a polymorphism

Understanding the relative importance of selection and stochastic factors in population divergence of adaptive traits is a classical topic in evolutionary biology. However, it is difficult to separate these factors and detect the effects of selection when two or more contrasting selective factors are simultaneously acting on a single locus. In the damselfly Ischnura senegalensis, females exhibit color dimorphism and morph frequencies change geographically. We here evaluated the role of selection and stochastic factors in population divergence of morph frequencies by comparing the divergences in color locus and neutral loci. Comparisons between population pairwise F_ST for neutral loci and for the color locus did not detect any stochastic factors affecting color locus. Although comparison between population divergence in color and neutral loci using all populations detected only divergent selection, we detected two antagonistic selective factors acting on the color locus, that is, balancing and divergent selection, when considering geographical distance between populations. Our results suggest that a combination of two antagonistic selective factors, rather than stochastic factors, establishes the geographic cline in morph frequency in this system.     

Mito‐nuclear discord in six congeneric lineages of Holarctic ducks (genus Anas)

Many species have Holarctic distributions that extend across Europe, Asia and North America. Most genetics research on these species has examined only mitochondrial (mt) DNA, which has revealed wide variance in divergence between Old World (OW) and New World (NW) populations, ranging from shallow, unstructured genealogies to deeply divergent lineages. In this study, we sequenced 20 nuclear introns to test for concordant patterns of OW–NW differentiation between mtDNA and nuclear (nu) DNA for six lineages of Holarctic ducks (genus Anas). Genetic differentiation for both marker types varied widely among these lineages (idiosyncratic population histories), but mtDNA and nuDNA divergence within lineages was not significantly correlated. Moreover, compared with the association between mtDNA and nuDNA divergence observed among different species, OW–NW nuDNA differentiation was generally lower than mtDNA divergence, at least for lineages with deeply divergent mtDNA. Furthermore, coalescent estimates indicated significantly higher rates of gene flow for nuDNA than mtDNA for four of the six lineages. Thus, Holarctic ducks show prominent mito‐nuclear discord between OW and NW populations, and we reject differences in sorting rates as the sole cause of the within‐species discord. Male‐mediated intercontinental gene flow is likely a leading contributor to this discord, although selection could also cause increased mtDNA divergence relative to weak nuDNA differentiation. The population genetics of these ducks contribute to growing evidence that mtDNA can be an unreliable indicator of stage of speciation and that more holistic approaches are needed for species delimitation.     

Local changes in GC/AT substitution biases and in crossover frequencies on Drosophila chromosomes

I present here evidence of remarkable local changes in GC/AT substitution biases and in crossover frequencies on Drosophila chromosomes. The substitution pattern at 10 loci in the telomeric region of the X chromosome was studied for four species of the Drosophila melanogaster species subgroup. Drosophila orena and Drosophila erecta are clearly the most closely related species pair (the erecta complex) among the four species studied; however, the overall data at the 10 loci revealed a clear dichotomy in the silent substitution patterns between the AT-biased- substitution melanogaster and erecta lineages and the GC-biased-substitution yakuba and orena lineages, suggesting two or more independent changes in GC/AT substitution biases. More importantly, the results indicated a between- loci heterogeneity in GC/AT substitution bias in this small region independently in the yakuba and orena lineages. Indeed, silent substitutions in the orena lineage were significantly biased toward G and C at the consecutive yellow, lethal of scute, and asense loci, but they were significantly biased toward A and T at sta. The substitution bias toward G and C was centered in different areas in yakuba (significantly biased at EG:165H7.3, EG:171D11.2, and suppressor of sable). The similar silent substitution patterns in coding and noncoding regions, furthermore, suggested mutational biases as a cause of the substitution biases. On the other hand, previous study reveals that Drosophila yakuba has about 20-fold higher crossover frequencies in the telomeric region of the X chromosome than does D. melanogaster; this study revealed that the total genetic map length of the yakuba X chromosome was only about 1.5 times as large as that of melanogaster and that the map length of the X-telomeric y-sta region did not differ between Drosophila yakuba and D. erecta. Taken together, the data strongly suggested that an approximately 20- fold reduction in the X-telomeric crossover frequencies occurred in the ancestral population of D. melanogaster after the melanogaster-yakuba divergence but before the melanogaster-simulans divergence.     

Quantitative trait loci affecting the difference in pigmentation between Drosophila yakuba and D. santomea

Using quantitative trait locus (QTL) mapping, we studied the genetic basis of the difference in pigmentation between two sister species of Drosophila: Drosophila yakuba, which, like other members of the D. melanogaster subgroup, shows heavy black pigmentation on the abdomen of males and females, and D. santomea, an endemic to the African island of S?o Tomé, which has virtually no pigmentation. Here we mapped four QTL with large effects on this interspecific difference in pigmentation: two on the X chromosome and one each on the second and third chromosomes. The same four QTL were detected in male hybrids in the backcrosses to both D. santomea and D. yakuba and in the female D. yakuba backcross hybrids. All four QTL exhibited strong epistatic interactions in male backcross hybrids, but only one pair of QTL interacted in females from the backcross to D. yabuka. All QTL from each species affected pigmentation in the same direction, consistent with adaptive evolution driven by directional natural selection. The regions delimited by the QTL included many positional candidate loci in the pigmentation pathway, including genes affecting catecholamine biosynthesis, melanization of the cuticle, and many additional pleiotropic effects.     

Selection from parasites favours immunogenetic diversity but not divergence among locally adapted host populations

The unprecedented polymorphism in the major histocompatibility complex (MHC) genes is thought to be maintained by balancing selection from parasites. However, do parasites also drive divergence at MHC loci between host populations, or do the effects of balancing selection maintain similarities among populations? We examined MHC variation in populations of the livebearing fish Poecilia mexicana and characterized their parasite communities. Poecilia mexicana populations in the Cueva del Azufre system are locally adapted to darkness and the presence of toxic hydrogen sulphide, representing highly divergent ecotypes or incipient species. Parasite communities differed significantly across populations, and populations with higher parasite loads had higher levels of diversity at class II MHC genes. However, despite different parasite communities, marked divergence in adaptive traits and in neutral genetic markers, we found MHC alleles to be remarkably similar among host populations. Our findings indicate that balancing selection from parasites maintains immunogenetic diversity of hosts, but this process does not promote MHC divergence in this system. On the contrary, we suggest that balancing selection on immunogenetic loci may outweigh divergent selection causing divergence, thereby hindering host divergence and speciation. Our findings support the hypothesis that balancing selection maintains MHC similarities among lineages during and after speciation (trans‐species evolution).     

Genomic variation in rice: genesis of highly polymorphic linkage blocks during domestication

Genomic regions that are unusually divergent between closely related species or racial groups can be particularly informative about the process of speciation or the operation of natural selection. The two sequenced genomes of cultivated Asian rice, Oryza sativa, reveal that at least 6% of the genomes are unusually divergent. Sequencing of ten unlinked loci from the highly divergent regions consistently identified two highly divergent haplotypes with each locus in nearly complete linkage disequilibrium among 25 O. sativa cultivars and 35 lines from six wild species. The existence of two highly divergent haplotypes in high divergence regions in species from all geographical areas (Africa, Asia, and Oceania) was in contrast to the low polymorphism and low linkage disequilibrium that were observed in other parts of the genome, represented by ten reference loci. While several natural processes are likely to contribute to this pattern of genomic variation, domestication may have greatly exaggerated the trend. In this hypothesis, divergent haplotypes that were adapted to different geographical and ecological environments migrated along with humans during the development of domesticated varieties. If true, these high divergence regions of the genome would be enriched for loci that contribute to the enormous range of phenotypic variation observed among domesticated breeds.     

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.     

The evolution of gene regulation underlies a morphological difference between two Drosophila sister species

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

Digest: The role of linkage in mimicking “magic traits”

Can divergence in a mating trait increase local adaption by increasing ecological divergence? Servedio and Bürger propose that “pseudomagic traits,” tightly linked complexes consisting of an ecological locus under divergent selection and a locus acting as a mating cue, can effectively mimic pleiotropy. Such pseudomagic traits can form even when linkage between ecological and mating loci is limited.     

Lineage pattern, trans-species polymorphism, and selection pressure among the major lineages of feline Mhc-DRB peptide-binding region

The long-term evolution of major histocompatibility complex (MHC) involves the birth-and-death process and independent divergence of loci during episodes punctuated by natural selection. Here, we investigated the molecular signatures of natural selection at exon-2 of MHC class II DRB gene which includes a part of the peptide-binding region (PBR) in seven of eight putative extant Felidae lineages. The DRB alleles in felids can be mainly divided into five lineages. Signatures of trans-species polymorphism among major allelic lineages indicate that balancing selection has maintained the MHC polymorphism for a long evolutionary time. Analysis based on maximum likelihood models of codon substitution revealed overall purifying selection acting on the feline DRB. Sites that have undergone positive selection and those that are under divergent selective pressure among lineages were detected and found to fall within the putative PBR. This study increased our understanding of the nature of selective forces acting on DRB during feline radiation.     

Nuclear gene genealogies reveal historical,demographic and selective factors associated with speciation in field crickets

Population genetics theory predicts that genetic drift should eliminate shared polymorphism, leading to monophyly or exclusivity of populations, when the elapsed time between lineage-splitting events is large relative to effective population size. We examined patterns of nucleotide variation in introns at four nuclear loci to relate processes affecting the history of genes to patterns of divergence among natural populations and species. Ancestral polymorphisms were shared among three recognized species, Gryllus firmus, G. pennsylvanicus, and G. ovisopis, and genealogical patterns suggest that successive speciation events occurred recently and rapidly relative to effective population size. High levels of shared polymorphism among these morphologically, behaviorally, and ecologically distinct species indicate that only a small fraction of the genome needs to become differentiated for speciation to occur. Among the four nuclear gene loci there was a 10-fold range in nucleotide diversity, and patterns of polymorphism and divergence suggest that natural selection has acted to maintain or eliminate variation at some loci. While nuclear gene genealogies may have limited applications in phylogeography or other approaches dependent on population monophyly, they provide important insights into the historical, demographic, and selective forces that shape speciation.     

Genetic diversity of MHC class I loci in six non-model frogs is shaped by positive selection and gene duplication

Comparative studies of major histocompatibility complex (MHC) genes across vertebrate species can reveal the evolutionary processes that shape the structure and function of immune regulatory proteins. In this study, we characterized MHC class I sequences from six frog species representing three anuran families (Hylidae, Centrolenidae and Ranidae). Using cDNA from our focal species, we amplified a total of 79 unique sequences spanning exons 2-4 that encode the extracellular domains of the functional alpha chain protein. We compared intra- and interspecific nucleotide and amino-acid divergence, tested for recombination, and identified codon sites under selection by estimating the rate of non-synonymous to synonymous substitutions with multiple codon-based maximum likelihood methods. We determined that positive (diversifying) selection was acting on specific amino-acid sites located within the domains that bind pathogen-derived peptides. We also found significant signals of recombination across the physical distance of the genes. Finally, we determined that all the six species expressed two or three putative classical class I loci, in contrast to the single locus condition of Xenopus laevis. Our results suggest that MHC evolution in anurans is a dynamic process and that variation in numbers of loci and genetic diversity can exist among taxa. Thus, the accumulation of genetic data for more species will be useful in further characterizing the relative importance of processes such as selection, recombination and gene duplication in shaping MHC loci among amphibian lineages.