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.     

Color vision and niche partitioning in a diverse neotropical primate community in lowland Amazonian Ecuador

A recent focus in community ecology has been on how within‐species variability shapes interspecific niche partitioning. Primate color vision offers a rich system in which to explore this issue. Most neotropical primates exhibit intraspecific variation in color vision due to allelic variation at the middle‐to‐long‐wavelength opsin gene on the X chromosome. Studies of opsin polymorphisms have typically sampled primates from different sites, limiting the ability to relate this genetic diversity to niche partitioning. We surveyed genetic variation in color vision of five primate species, belonging to all three families of the primate infraorder Platyrrhini, found in the Yasuní Biosphere Reserve in Ecuador. The frugivorous spider monkeys and woolly monkeys (Ateles belzebuth and Lagothrix lagotricha poeppigii, family Atelidae) each had two opsin alleles, and more than 75% of individuals carried the longest‐wavelength (553–556 nm) allele. Among the other species, Saimiri sciureus macrodon (family Cebidae) and Pithecia aequatorialis (family Pitheciidae) had three alleles, while Plecturocebus discolor (family Pitheciidae) had four alleles—the largest number yet identified in a wild population of titi monkeys. For all three non‐atelid species, the middle‐wavelength (545 nm) allele was the most common. Overall, we identified genetic evidence of fourteen different visual phenotypes—seven types of dichromats and seven trichromats—among the five sympatric taxa. The differences we found suggest that interspecific competition among primates may influence intraspecific frequencies of opsin alleles. The diversity we describe invites detailed study of foraging behavior of different vision phenotypes to learn how they may contribute to niche partitioning.     

Melanic pigmentation and light preference within and between two Drosophila species

Environmental adaptation and species divergence often involve suites of co‐evolving traits. Pigmentation in insects presents a variable, adaptive, and well‐characterized class of phenotypes for which correlations with multiple other traits have been demonstrated. In Drosophila, the pigmentation genes ebony and tan have pleiotropic effects on flies'' response to light, creating the potential for correlated evolution of pigmentation and vision. Here, we investigate differences in light preference within and between two sister species, Drosophila americana and D. novamexicana, which differ in pigmentation in part because of evolution at ebony and tan and occupy environments that differ in many variables including solar radiation. We hypothesized that lighter pigmentation would be correlated with a greater preference for environmental light and tested this hypothesis using a habitat choice experiment. In a first set of experiments, using males of D. novamexicana line N14 and D. americana line A00, the light‐bodied D. novamexicana was found slightly but significantly more often than D. americana in the light habitat. A second experiment, which included additional lines and females as well as males, failed to find any significant difference between D. novamexicana‐N14 and D. americana‐A00. Additionally, the other dark line of D. americana (A04) was found in the light habitat more often than the light‐bodied D. novamexicana‐N14, in contrast to our predictions. However, the lightest line of D. americana, A01, was found substantially and significantly more often in the light habitat than the two darker lines of D. americana, thus providing partial support for our hypothesis. Finally, across all four lines, females were found more often in the light habitat than their more darkly pigmented male counterparts. Additional replication is needed to corroborate these findings and evaluate conflicting results, with the consistent effect of sex within and between species providing an especially intriguing avenue for further research.     

RAD sequencing of common whelk,Buccinum undatum,reveals fine‐scale population structuring in Europe and cryptic speciation within the North Atlantic

Buccinum undatum is a subtidal gastropod that exhibits clear spatial variation in several phenotypic shell traits (color, shape, and thickness) across its North Atlantic distribution. Studies of spatial phenotypic variation exist for the species; however, population genetic studies have thus far relied on a limited set of mitochondrial and microsatellite markers. Here, we greatly expand on previous work by characterizing population genetic structure in B. undatum across the North Atlantic from SNP variation obtained by RAD sequencing. There was a high degree of genetic differentiation between Canadian and European populations (Iceland, Faroe Islands, and England) consistent with the divergence of populations in allopatry (F _ST > 0.57 for all pairwise comparisons). In addition, B. undatum populations within Iceland, the Faroe Islands, and England are typified by weak but significant genetic structuring following an isolation‐by‐distance model. Finally, we established a significant correlation between genetic structuring in Iceland and two phenotypic traits: shell shape and color frequency. The works detailed here enhance our understanding of genetic structuring in B. undatum and establish the species as an intriguing model for future genome‐wide association studies.     

Proteomic analysis of Drosophila mojavensis male accessory glands suggests novel classes of seminal fluid proteins

Fruit-flies of the genus Drosophila are characterized by overwhelming variation in fertilization traits such as copulatory plug formation, sperm storage organ use, and nutritional ejaculatory donation. Despite extensive research on the genetic model Drosophila melanogaster, little is known about the molecular underpinnings of these interspecific differences. This study employs a proteomic approach to pin-point candidate seminal fluid proteins in Drosophila mojavensis, a cactophilic fruit-fly that exhibits divergent reproductive biology when compared to D. melanogaster. We identify several classes of candidate seminal fluid proteins not previously documented in the D. melanogaster male ejaculate, including metabolic enzymes, nutrient transport proteins, and clotting factors. Conversely, we also define 29 SFPs that are conserved despite >40 million years of Drosophila evolution. We discuss our results in terms of universal processes in insect reproduction, as well as the specialized reproductive biology of D. mojavensis.     

Coadapted changes in energy metabolites and body color phenotypes for resistance to starvation and desiccation in latitudinal populations of D. melanogaster

In D. melanogaster, resistance to starvation and desiccation vary in opposite directions across a geographical gradient in India but there is lack of such clinal variation on other continents. However, it is not clear whether these resistance traits or other correlated traits are the target of natural selection. For resistance to starvation or desiccation in D. melanogaster, we tested the hypothesis whether body color phenotypes and energy metabolites show correlated selection response. Our results are interesting in several respects. First, based on within population analysis, assorted darker and lighter flies from a given population showed that darker flies store higher amount of trehalose and confer greater desiccation resistance as compared with lighter flies. By contrast, lighter flies store higher lipids content and confer increased starvation tolerance. Thus, there is a trade-off for energy metabolites as well as body color phenotypes for starvation and desiccation stress. Further, trait associations within populations reflect similar patterns in geographical populations. Second, we found opposite clines for trehalose and body lipids. Third, coadapated phenotypes have evolved under contrasting climatic conditions i.e. drier and colder northern localities select darker flies with higher trehalose as well as desiccation resistance while hot and humid localities favor lighter flies with higher lipids level and greater starvation tolerance. Thus, the evolution of coadapated phenotypes associated with starvation and desiccation resistance might have resulted due to specific ecological conditions i.e. humidity changes on the Indian subcontinent.     

Uncovering the genetic diversity of yams (Dioscorea spp.) in China by combining phenotypic trait and molecular marker analyses

Yam is an important edible tuber and root plant worldwide; China as one of the native places of yams has many diverse local resources. The goal of this study was to clarify the genetic diversity of the commonly cultivated yam landraces and the genetic relationship between the main yam species in China. In this study, 26 phenotypic traits of 112 yam accessions from 21 provinces in China were evaluated, and 24 simple sequence repeat (SSR) and 29 sequence‐related amplified polymorphism (SRAP) markers were used for the genetic diversity analysis. Phenotypic traits revealed that Dioscorea opposita had the highest genetic diversity, followed by D. alata, D. persimilis, D. fordii, and D. esculenta. Among the 26 phenotypic traits, the Shannon diversity indexes of leaf shape, petiole color, and stem color were high, and the range in the variation of tuber‐related traits in the underground part was higher than that in the aboveground part. All accessions were divided into six groups by phenotypic trait clustering, which was also supported by principal component analysis (PCA). Molecular marker analysis showed that SSR and SRAP markers had good amplification effects and could effectively and accurately evaluate the genetic variation of yam. The unweighted pair‐group method with arithmetic means analysis based on SSR‐SRAP marker data showed that the 112 accessions were also divided into six groups, similar to the phenotypic trait results. The results of PCA and population structure analysis based on SSR‐SRAP data also produced similar results. In addition, the analysis of the origin and genetic relationship of yam indicated that the species D. opposita may have originated from China. These results demonstrate the genetic diversity and distinctness among the widely cultivated species of Chinese yam and provide a theoretical reference for the classification, breeding, germplasm innovation, utilization, and variety protection of Chinese yam resources.     

Genetic basis of natural variation in body pigmentation in Drosophila melanogaster

Body pigmentation in insects and other organisms is typically variable within and between species and is often associated with fitness. Regulatory variants with large effects at bab1, t and e affect variation in abdominal pigmentation in several populations of Drosophila melanogaster. Recently, we performed a genome wide association (GWA) analysis of variation in abdominal pigmentation using the inbred, sequenced lines of the Drosophila Genetic Reference Panel (DGRP). We confirmed the large effects of regulatory variants in bab1, t and e; identified 81 additional candidate genes; and validated 17 candidate genes (out of 28 tested) using RNAi knockdown of gene expression and mutant alleles. However, these analyses are imperfect proxies for the effects of segregating variants. Here, we describe the results of an extreme quantitative trait locus (xQTL) GWA analysis of female body pigmentation in an outbred population derived from light and dark DGRP lines. We replicated the effects on pigmentation of 28 genes implicated by the DGRP GWA study, including bab1, t and e and 7 genes previously validated by RNAi and/or mutant analyses. We also identified many additional loci. The genetic architecture of Drosophila pigmentation is complex, with a few major genes and many other loci with smaller effects.     

Evolution of wing pigmentation in Drosophila: Diversity,physiological regulation,and cis-regulatory evolution

Fruit flies (Drosophila and its close relatives, or “drosophilids”) are a group that includes an important model organism, Drosophila melanogaster, and also very diverse species distributed worldwide. Many of these species have black or brown pigmentation patterns on their wings, and have been used as material for evo-devo research. Pigmentation patterns are thought to have evolved rapidly compared with body plans or body shapes; hence they are advantageous model systems for studying evolutionary gains of traits and parallel evolution. Various groups of drosophilids, including genus Idiomyia (Hawaiian Drosophila), have a variety of pigmentations, ranging from simple black pigmentations around crossveins to a single antero-distal spot and a more complex mottled pattern. Pigmentation patterns are sometimes obviously used for sexual displays; however, in some cases they may have other functions. The process of wing formation in Drosophila, the general mechanism of pigmentation formation, and the transport of substances necessary for pigmentation, including melanin precursors, through wing veins are summarized here. Lastly, the evolution of the expression of genes regulating pigmentation patterns, the role of cis-regulatory regions, and the conditions required for the evolutionary emergence of pigmentation patterns are discussed. Future prospects for research on the evolution of wing pigmentation pattern formation in drosophilids are presented, particularly from the point of view of how they compare with other studies of the evolution of new traits.     

Plasticity and artificial selection for developmental mode in a poecilogonous sea slug

The contribution of phenotypically plastic traits to evolution depends on the degree of environmental influence on the target of selection (the phenotype) as well as the underlying genetic structure of the trait and plastic response. Likewise, maternal effects can help or hinder evolution through affects to the response to selection. The sacoglossan sea slug Alderia willowi exhibits intraspecific variation for developmental mode (= poecilogony) that is environmentally modulated with populations producing more yolk‐feeding (lecithotrophic) larvae during the summer, and more planktonic‐feeding (planktotrophic) larvae in the winter. I found significant family‐level variation in the reaction norms between 17 maternal families of A. willowi when reared in a split‐brood design in low (16 ppt) versus high (32 ppt) salinity, conditions which mimic seasonal variation in salinity of natural populations. I documented a significant response to selection for lecithotrophic larvae in high and low salinity. The slope of the reaction norm was maintained following one generation of selection for lecithotrophy. When the maternal environment was controlled in the laboratory, I found significant maternal effects, which reduced the response to selection. These results suggest there is standing genetic variation for egg‐mass type in A. willowi, but the ability of selection to act on that variation may depend on the environment in which the phenotype is expressed in preceding generations.     

Realized evolvability: quantifying phenotypic evolution in a Drosophila clade

Understanding the evolutionary potential of morphological evolution is still a major problem in evolutionary biology. In this study, we tried to quantify the amount of variation of different traits among species of a Drosophila clade reared under standard conditions. Nineteen different traits have been measured on nine species of the same clade, the Neotropical saltans group of Drosophila. Measured traits can be distributed into five categories: size‐traits (wing and thorax), shape indices (ratios among the size traits), number of sternopleural bristles on the thorax, number of abdominal bristles on successive sternites, and dorsal pigmentation of abdomen. All species are of medium size with a generally dark pigmentation. A remarkable feature is the presence of numerous bristles on the 6th sternite of the males, while this segment is bare in other Drosophila species. A multivariate analysis revealed that it was possible to discriminate all the investigated species by using the complete set of measured traits. For each trait, phenotypic variability was investigated at the within‐ and between‐species levels. As a rule, the interspecific coefficient of variation (CV) was much greater than the within species CV, and it is proposed to call it realized evolvability. All possible correlations were calculated between traits across species, providing many unexpected results. Size traits were highly correlated among them, but not correlated with shape indices. Abdominal traits (bristles and pigmentation) were correlated, but often in opposite directions, with thoracic shape indices. Tergite pigmentation was negatively correlated with bristle number on sternite. For the moment, most of the correlations cannot be explained by developmental processes or parallel selective pressures. Nonetheless, mapping the evolution of the two characters on a molecular phylogeny of the studied species revealed two opposite phylogenetic trends for abdominal pigmentation and setation, respectively. Our data suggest a need for similar studies in other well‐defined Drosophila clades.     

Genome-Wide Association Studies of Quantitatively Measured Skin,Hair, and Eye Pigmentation in Four European Populations

Pigmentation of the skin, hair, and eyes varies both within and between human populations. Identifying the genes and alleles underlying this variation has been the goal of many candidate gene and several genome-wide association studies (GWAS). Most GWAS for pigmentary traits to date have been based on subjective phenotypes using categorical scales. But skin, hair, and eye pigmentation vary continuously. Here, we seek to characterize quantitative variation in these traits objectively and accurately and to determine their genetic basis. Objective and quantitative measures of skin, hair, and eye color were made using reflectance or digital spectroscopy in Europeans from Ireland, Poland, Italy, and Portugal. A GWAS was conducted for the three quantitative pigmentation phenotypes in 176 women across 313,763 SNP loci, and replication of the most significant associations was attempted in a sample of 294 European men and women from the same countries. We find that the pigmentation phenotypes are highly stratified along axes of European genetic differentiation. The country of sampling explains approximately 35% of the variation in skin pigmentation, 31% of the variation in hair pigmentation, and 40% of the variation in eye pigmentation. All three quantitative phenotypes are correlated with each other. In our two-stage association study, we reproduce the association of rs1667394 at the OCA2/HERC2 locus with eye color but we do not identify new genetic determinants of skin and hair pigmentation supporting the lack of major genes affecting skin and hair color variation within Europe and suggesting that not only careful phenotyping but also larger cohorts are required to understand the genetic architecture of these complex quantitative traits. Interestingly, we also see that in each of these four populations, men are more lightly pigmented in the unexposed skin of the inner arm than women, a fact that is underappreciated and may vary across the world.     

Gudgeon fish with and without genetically determined countershading coexist in heterogeneous littoral environments of an ancient lake

Countershading, characterized by a darker dorsal surface and lighter ventral surface, is common among many animals. This dorsoventral pigment polarity is often thought to be adaptive coloration for camouflage. By contrast, noncountershaded (melanistic) morphs often occur within a species due to genetic color polymorphism in terrestrial animals. However, the polymorphism with either countershaded or melanistic morphs is poorly known in wild aquatic animals. This study explored the genetic nature of diverged color morphs of a lineage of gudgeon fish (genus Sarcocheilichthys) in the ancient Lake Biwa and propose this system as a novel model for testing hypotheses of functional aspects of countershading and its loss in aquatic environments. This system harbors two color morphs that have been treated taxonomically as separate species; Sarcocheilichthys variegatus microoculus which occurs throughout the littoral zone and Sarcocheilichthys biwaensis which occurs in and around rocky areas. First, we confirmed that the divergence of dorsoventral color patterns between the two morphs is under strict genetic control at the levels of chromatophore distribution and melanin‐related gene expression under common garden rearing. The former morph displayed sharp countershading coloration, whereas the latter morph exhibited a strong tendency toward its loss. The crossing results indicated that this divergence was likely controlled by a single locus in a two‐allele Mendelian inheritance pattern. Furthermore, our population genomic and genome‐wide association study analyses detected no genome‐wide divergence between the two morphs, except for one region near a locus that may be associated with the color divergence. Thus, these morphs are either in a state of intraspecific color polymorphism or two incipient species. Evolutionary forces underlying this polymorphism appear to be associated with heterogeneous littoral environments in this lake. Future ecological genomic research will provide insight into adaptive functions of this widespread coloration, including the eco‐evolutionary drivers of its loss, in the aquatic world.     

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.     

Genotype and sex-based host variation in behaviour and susceptibility drives population disease dynamics

Host heterogeneity in pathogen transmission is widespread and presents a major hurdle to predicting and minimizing disease outbreaks. Using Drosophila melanogaster infected with Drosophila C virus as a model system, we integrated experimental measurements of social aggregation, virus shedding, and disease-induced mortality from different genetic lines and sexes into a disease modelling framework. The experimentally measured host heterogeneity produced substantial differences in simulated disease outbreaks, providing evidence for genetic and sex-specific effects on disease dynamics at a population level. While this was true for homogeneous populations of single sex/genetic line, the genetic background or sex of the index case did not alter outbreak dynamics in simulated, heterogeneous populations. Finally, to explore the relative effects of social aggregation, viral shedding and mortality, we compared simulations where we allowed these traits to vary, as measured experimentally, to simulations where we constrained variation in these traits to the population mean. In this context, variation in infectiousness, followed by social aggregation, was the most influential component of transmission. Overall, we show that host heterogeneity in three host traits dramatically affects population-level transmission, but the relative impact of this variation depends on both the susceptible population diversity and the distribution of population-level variation.     

Many ways to be small: different environmental regulators of size generate distinct scaling relationships in Drosophila melanogaster

Static allometries, the scaling relationship between body and trait size, describe the shape of animals in a population or species, and are generated in response to variation in genetic or environmental regulators of size. In principle, allometries may vary with the different size regulators that generate them, which can be problematic since allometric differences are also used to infer patterns of selection on morphology. We test this hypothesis by examining the patterns of scaling in Drosophila melanogaster subjected to variation in three environmental regulators of size: nutrition, temperature and rearing density. Our data indicate that different environmental regulators of size do indeed generate different patterns of scaling. Consequently, flies that are ostensibly the same size may have very different body proportions. These data indicate that trait size is not simply a read-out of body size, but that different environmental factors may regulate body and trait size, and the relationship between the two, through different developmental mechanisms. It may therefore be difficult to infer selective pressures that shape scaling relationships in a wild population without first elucidating the environmental and genetic factors that generate size variation among members of the population.