Reconciling Differences in Pool-GWAS Between Populations: A Case Study of Female Abdominal Pigmentation in Drosophila melanogaster

The Genetic Mechanisms Underlying the Concerted Expression of the yellow and tan Genes in Complex Patterns on the Abdomen and Wings of Drosophila guttifera

How complex morphological patterns form is an intriguing question in developmental biology. However, the mechanisms that generate complex patterns remain largely unknown. Here, we sought to identify the genetic mechanisms that regulate the tan (t) gene in a multi-spotted pigmentation pattern on the abdomen and wings of Drosophila guttifera. Previously, we showed that yellow (y) gene expression completely prefigures the abdominal and wing pigment patterns of this species. In the current study, we demonstrate that the t gene is co-expressed with the y gene in nearly identical patterns, both transcripts foreshadowing the adult abdominal and wing melanin spot patterns. We identified cis-regulatory modules (CRMs) of t, one of which drives reporter expression in six longitudinal rows of spots on the developing pupal abdomen, while the second CRM activates the reporter gene in a spotted wing pattern. Comparing the abdominal spot CRMs of y and t, we found a similar composition of putative transcription factor binding sites that are thought to regulate the complex expression patterns of both terminal pigmentation genes y and t. In contrast, the y and t wing spots appear to be regulated by distinct upstream factors. Our results suggest that the D. guttifera abdominal and wing melanin spot patterns have been established through the co-regulation of y and t, shedding light on how complex morphological traits may be regulated through the parallel coordination of downstream target genes.

Genomic architecture of phenotypic extremes in a wild cervid

Identifying the genes underlying fitness-related traits such as body size and male ornamentation can provide tools for conservation and management and are often subject to various selective pressures. Here we performed high-depth whole genome re-sequencing of pools of individuals representing the phenotypic extremes for antler and body size in white-tailed deer (Odocoileus virginianus). Samples were selected from a tissue repository containing phenotypic data for 4,466 male white-tailed deer from Anticosti Island, Quebec, with four pools representing the extreme phenotypes for antler and body size after controlling for age. Our results revealed a largely homogenous population but detected highly divergent windows between pools for both traits, with the mean allele frequency difference of 14% for and 13% for antler and body SNPs in outlier windows, respectively. Genes in outlier antler windows were enriched for pathways associated with cell death and protein metabolism and some of the most differentiated windows included genes associated with oncogenic pathways and reproduction, processes consistent with antler evolution and growth. Genes associated with body size were more nuanced, suggestive of a highly complex trait. Overall, this study revealed the complex genomic make-up of both antler morphology and body size in free-ranging white-tailed deer and identified target loci for additional analyses.

Genome-wide patterns of genetic diversity, population structure and demographic history in mānuka (Leptospermum scoparium) growing on indigenous Māori land

Leptospermum scoparium J. R. Forst et G. Forst, known as mānuka by Māori, the indigenous people of Aotearoa (New Zealand), is a culturally and economically significant shrub species, native to New Zealand and Australia. Chemical, morphological and phylogenetic studies have indicated geographical variation of mānuka across its range in New Zealand, and genetic differentiation between New Zealand and Australia. We used pooled whole genome re-sequencing of 76 L. scoparium and outgroup populations from New Zealand and Australia to compile a dataset totalling ~2.5 million SNPs. We explored the genetic structure and relatedness of L. scoparium across New Zealand, and between populations in New Zealand and Australia, as well as the complex demographic history of this species. Our population genomic investigation suggests there are five geographically distinct mānuka gene pools within New Zealand, with evidence of gene flow occurring between these pools. Demographic modelling suggests three of these gene pools have undergone expansion events, whilst the evolutionary histories of the remaining two have been subjected to contractions. Furthermore, mānuka populations in New Zealand are genetically distinct from populations in Australia, with coalescent modelling suggesting these two clades diverged ~9-12 million years ago. We discuss the evolutionary history of this species and the benefits of using pool-seq for such studies. Our research will support the management and conservation of mānuka by landowners, particularly Māori, and the development of a provenance story for the branding of mānuka based products.

Genomic basis for drought resistance in European beech forests threatened by climate change

In the course of global climate change, Central Europe is experiencing more frequent and prolonged periods of drought. The drought years 2018 and 2019 affected European beeches (Fagus sylvatica L.) differently: even in the same stand, drought-damaged trees neighboured healthy trees, suggesting that the genotype rather than the environment was responsible for this conspicuous pattern. We used this natural experiment to study the genomic basis of drought resistance with Pool-GWAS. Contrasting the extreme phenotypes identified 106 significantly associated single-nucleotide polymorphisms (SNPs) throughout the genome. Most annotated genes with associated SNPs (>70%) were previously implicated in the drought reaction of plants. Non-synonymous substitutions led either to a functional amino acid exchange or premature termination. An SNP assay with 70 loci allowed predicting drought phenotype in 98.6% of a validation sample of 92 trees. Drought resistance in European beech is a moderately polygenic trait that should respond well to natural selection, selective management, and breeding.

Climate change is having a serious impact on many ecosystems. In the summer of 2018 and 2019, around two thirds of European beech trees were damaged or killed by extreme drought. It is critical to keep these beech woods healthy, as they are central to the survival of over 6,000 other species of animals and plants.

The level of damage caused by the drought varied between forests. However, not all the trees in each forest responded in the same way, with severely damaged trees often sitting next to fully healthy ones. This suggests that the genetic make-up of each tree determines how well it can adapt to drought rather than its local environment.

To investigate this further, Pfenninger et al. studied the genome of over 400 European beech trees from the Hesse region in Germany. The samples came from pairs of neighbouring trees that had responded differently to the droughts. The analysis found more than 80 parts of the genome that differed between healthy and damaged trees.

Pfenninger et al. then used this information to create a genetic test which can quickly and inexpensively predict how well an individual beech tree might survive in a drought. Applying this test to another 92 trees revealed that it can reliably detect which ones were healthy and which ones were damaged.

Beech forests are typically managed by private owners, agencies or breeders that could use this genetic test to select and reproduce trees that are better adapted to drought. The goal now is to develop the test so that it can be used more widely to manage European beech trees and potentially other species.

Many ways to make darker flies: Intra- and interspecific variation in Drosophila body pigmentation components

Body pigmentation is an evolutionarily diversified and ecologically relevant trait with substantial variation within and between species, and important roles in animal survival and reproduction. Insect pigmentation, in particular, provides some of the most compelling examples of adaptive evolution, including its ecological significance and genetic bases. Pigmentation includes multiple aspects of color and color pattern that may vary more or less independently, and can be under different selective pressures. We decompose Drosophila thorax and abdominal pigmentation, a valuable eco-evo-devo model, into distinct measurable traits related to color and color pattern. We investigate intra- and interspecific variation for those traits and assess its different sources. For each body part, we measured overall darkness, as well as four other pigmentation properties distinguishing between background color and color of the darker pattern elements that decorate each body part. By focusing on two standard D. melanogaster laboratory populations, we show that pigmentation components vary and covary in distinct manners depending on sex, genetic background, and temperature during development. Studying three natural populations of D. melanogaster along a latitudinal cline and five other Drosophila species, we then show that evolution of lighter or darker bodies can be achieved by changing distinct component traits. Our results paint a much more complex picture of body pigmentation variation than previous studies could uncover, including patterns of sexual dimorphism, thermal plasticity, and interspecific diversity. These findings underscore the value of detailed quantitative phenotyping and analysis of different sources of variation for a better understanding of phenotypic variation and diversification, and the ecological pressures and genetic mechanisms underlying them.

[Phenotypic plasticity in insects]

Insects represent 85% of the animals. They have adapted to many environments and play a major role in ecosystems. Many insect species exhibit phenotypic plasticity. We here report on the mechanisms involved in phenotypic plasticity of different insects (aphids, migratory locust, map butterfly, honeybee) and also on the nutritional size plasticity in Drosophila and the plasticity of the wing eye-spots of the butterfly Bicyclus anynana. We also describe in more detail our work concerning the thermal plasticity of pigmentation in Drosophila. We have shown that the expression of the tan, yellow and Ddc genes, encoding enzymes of the melanin synthesis pathway, is modulated by temperature and that it is a consequence, at least in part, of the temperature-sensitive expression of the bab locus genes that repress them.

Genetic architecture of a body colour cline in Drosophila americana

Phenotypic variation within a species is often structured geographically in clines. In Drosophila americana, a longitudinal cline for body colour exists within North America that appears to be due to local adaptation. The tan and ebony genes have been hypothesized to contribute to this cline, with alleles of both genes that lighten body colour found in D. americana. These alleles are similar in sequence and function to the allele fixed in D. americana's more lightly pigmented sister species, Drosophila novamexicana. Here, we examine the frequency and geographic distribution of these D. novamexicana-like alleles in D. americana. Among alleles from over 100 strains of D. americana isolated from 21 geographic locations, we failed to identify additional alleles of tan or ebony with as much sequence similarity to D. novamexicana as the D. novamexicana-like alleles previously described. However, using genetic analysis of 51 D. americana strains derived from 20 geographic locations, we identified one new allele of ebony and one new allele of tan segregating in D. americana that are functionally equivalent to the D. novamexicana allele. An additional 5 alleles of tan also showed marginal evidence of functional similarity. Given the rarity of these alleles, however, we conclude that they are unlikely to be driving the pigmentation cline. Indeed, phenotypic distributions of the 51 backcross populations analysed indicate a more complex genetic architecture, with diversity in the number and effects of loci altering pigmentation observed both within and among populations of D. americana. This genetic heterogeneity poses a challenge to association studies and genomic scans for clinal variation, but might be common in natural populations.

Shared genetic susceptibilities for irritable bowel syndrome and depressive disorder in Chinese patients uncovered by pooled whole-exome sequencing

Irritable bowel syndrome (IBS) is the most prevalent functional gastrointestinal disorder presenting a high comorbidity with depressive disorder (DD). Many studies have confirmed that these two disease share the similar pathophysiological process, but evidence of the genetic risks is limited. This study aimed to analyze the genetic susceptibilities for IBS and DD in Chinese patients. Pooled whole-exome sequencing (pooled-WES) was performed to identify the candidate variants in the group of diarrhea predominant IBS (IBS-D) patients, DD patients, and healthy controls (HC). Then, targeted sequencing was used to validate the candidate variants in three additional cohorts of IBS-D, DD, and HC. Four variants associated with both IBS-D and DD were identified through pooled-WES, and three of them were validated in targeted sequencing. SYT8 rs3741231 G allele and SSPO rs12536873 TT genotype were associated with both IBS-D and DD. The genes of these variants are important in neurogenesis and neurotransmission. In addition, we found COL6A1 rs13051496, a unique risk variation for IBS-D. It increased the IBS-D risk and had a positive correlation with the scores of abdominal bloating and dissatisfaction of bowel habits. Through the results of this study, it provides a genetic basis for the high comorbidity of IBS-D and DD.

Genomic basis of European ash tree resistance to ash dieback fungus

Populations of European ash trees (Fraxinus excelsior) are being devastated by the invasive alien fungus Hymenoscyphus fraxineus, which causes ash dieback. We sequenced whole genomic DNA from 1,250 ash trees in 31 DNA pools, each pool containing trees with the same ash dieback damage status in a screening trial and from the same seed-source zone. A genome-wide association study identified 3,149 single nucleotide polymorphisms (SNPs) associated with low versus high ash dieback damage. Sixty-one of the 192 most significant SNPs were in, or close to, genes with putative homologues already known to be involved in pathogen responses in other plant species. We also used the pooled sequence data to train a genomic prediction model, cross-validated using individual whole genome sequence data generated for 75 healthy and 75 damaged trees from a single seed source. The model's genomic estimated breeding values (GEBVs) allocated these 150 trees to their observed health statuses with 67% accuracy using 10,000 SNPs. Using the top 20% of GEBVs from just 200 SNPs, we could predict observed tree health with over 90% accuracy. We infer that ash dieback resistance in F. excelsior is a polygenic trait that should respond well to both natural selection and breeding, which could be accelerated using genomic prediction.

Pleiotropic Effects of ebony and tan on Pigmentation and Cuticular Hydrocarbon Composition in Drosophila melanogaster

Pleiotropic genes are genes that affect more than one trait. For example, many genes required for pigmentation in the fruit fly Drosophila melanogaster also affect traits such as circadian rhythms, vision, and mating behavior. Here, we present evidence that two pigmentation genes, ebony and tan, which encode enzymes catalyzing reciprocal reactions in the melanin biosynthesis pathway, also affect cuticular hydrocarbon (CHC) composition in D. melanogaster females. More specifically, we report that ebony loss-of-function mutants have a CHC profile that is biased toward long (>25C) chain CHCs, whereas tan loss-of-function mutants have a CHC profile that is biased toward short (<25C) chain CHCs. Moreover, pharmacological inhibition of dopamine synthesis, a key step in the melanin synthesis pathway, reversed the changes in CHC composition seen in ebony mutants, making the CHC profiles similar to those seen in tan mutants. These observations suggest that genetic variation affecting ebony and/or tan activity might cause correlated changes in pigmentation and CHC composition in natural populations. We tested this possibility using the Drosophila Genetic Reference Panel (DGRP) and found that CHC composition covaried with pigmentation as well as levels of ebony and tan expression in newly eclosed adults in a manner consistent with the ebony and tan mutant phenotypes. These data suggest that the pleiotropic effects of ebony and tan might contribute to covariation of pigmentation and CHC profiles in Drosophila.

bric à brac (bab), a central player in the gene regulatory network that mediates thermal plasticity of pigmentation in Drosophila melanogaster

Drosophila body pigmentation has emerged as a major Evo-Devo model. Using two Drosophila melanogaster lines, Dark and Pale, selected from a natural population, we analyse here the interaction between genetic variation and environmental factors to produce this complex trait. Indeed, pigmentation varies with genotype in natural populations and is sensitive to temperature during development. We demonstrate that the bric à brac (bab) genes, that are differentially expressed between the two lines and whose expression levels vary with temperature, participate in the pigmentation difference between the Dark and Pale lines. The two lines differ in a bab regulatory sequence, the dimorphic element (called here bDE). Both bDE alleles are temperature-sensitive, but the activity of the bDE allele from the Dark line is lower than that of the bDE allele from the Pale line. Our results suggest that this difference could partly be due to differential regulation by AbdB. bab has been previously reported to be a repressor of abdominal pigmentation. We show here that one of its targets in this process is the pigmentation gene tan (t), regulated via the tan abdominal enhancer (t_MSE). Furthermore, t expression is strongly modulated by temperature in the two lines. Thus, temperature sensitivity of t expression is at least partly a consequence of bab thermal transcriptional plasticity. We therefore propose that a gene regulatory network integrating both genetic variation and temperature sensitivity modulates female abdominal pigmentation. Interestingly, both bDE and t_MSE were previously shown to have been recurrently involved in abdominal pigmentation evolution in drosophilids. We propose that the environmental sensitivity of these enhancers has turned them into evolutionary hotspots.

Complex traits such as size or disease susceptibility are typically modulated by both genetic variation and environmental conditions. Model organisms such as fruit flies (Drosophila) are particularly appropriate to analyse the interactions between genetic variation and environmental factors during the development of complex phenotypes. Natural populations carry high genetic variation and can be grown in controlled conditions in the laboratory. Here, we use Drosophila melanogaster female abdominal pigmentation, which is both genetically variable and modulated by the environment (temperature) to dissect this kind of interaction. We show that the pigmentation difference between two inbred fly lines is caused by genetic variation in an enhancer of the bab locus, which encodes two transcription factors controlling abdominal pigmentation. Indeed, this enhancer drives differential expression between the two lines. Interestingly, this enhancer is sensitive to temperature in both lines. We show that the effect of bab on pigmentation is mediated by the pigmentation gene tan (t) that is repressed by bab. Thus, the previously reported temperature-sensitive expression of t is a direct consequence of bab transcriptional plasticity.

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.

Genomics of coloration in natural animal populations

Animal coloration has traditionally been the target of genetic and evolutionary studies. However, until very recently, the study of the genetic basis of animal coloration has been mainly restricted to model species, whereas research on non-model species has been either neglected or mainly based on candidate approaches, and thereby limited by the knowledge obtained in model species. Recent high-throughput sequencing technologies allow us to overcome previous limitations, and open new avenues to study the genetic basis of animal coloration in a broader number of species and colour traits, and to address the general relevance of different genetic structures and their implications for the evolution of colour. In this review, we highlight aspects where genome-wide studies could be of major utility to fill in the gaps in our understanding of the biology and evolution of animal coloration. The new genomic approaches have been promptly adopted to study animal coloration although substantial work is still needed to consider a larger range of species and colour traits, such as those exhibiting continuous variation or based on reflective structures. We argue that a robust advancement in the study of animal coloration will also require large efforts to validate the functional role of the genes and variants discovered using genome-wide tools.This article is part of the themed issue 'Animal coloration: production, perception, function and application'.

The flexible stem hypothesis: evidence from genetic data

Phenotypic plasticity, the ability of a given genotype to produce different phenotypes in response to distinct environmental conditions, is widely observed in the wild. It is believed to facilitate evolution and, under the "flexible stem hypothesis", it is thought that an ancestral plastic species can be at the origin of sister lineages with divergent phenotypes fixed by genetic assimilation of alternative morphs. We review here the genetic mechanisms underlying such phenomenon. We show several examples in which the same gene shows transcriptional plasticity in response to environmental factors and divergence of expression within or between species. Thus, the same gene is involved both in the plasticity of a trait and in the evolution of that trait. In a few cases, it can be traced down to cis-regulatory variation in this gene and, in one case, in the very same regulatory sequence whose activity is modulated by the environment. These data are compatible with the "flexible stem hypothesis" and also suggest that the evolution of the plasticity of a trait and the evolution of the trait are not completely uncoupled as they often involve the same locus. Furthermore, the "flexible stem hypothesis" implies that it is possible to canalize initially plastic phenotypes. Several studies have shown that it was possible through modification of chromatin regulation or hormonal signalling/response. Further studies of phenotypic plasticity in an evolutionary framework are needed to see how much the findings described in this review can be generalized.

Strong epistatic and additive effects of linked candidate SNPs for Drosophila pigmentation have implications for analysis of genome-wide association studies results

Background

The mapping resolution of genome-wide association studies (GWAS) is limited by historic recombination events and effects are often assigned to haplotype blocks rather than individual SNPs. It is not clear how many of the SNPs in the block, and which ones, are causative. Drosophila pigmentation is a powerful model to dissect the genetic basis of intra-specific and inter-specific phenotypic variation. Three tightly linked SNPs in the t-MSE enhancer have been identified in three D. melanogaster populations as major contributors to female abdominal pigmentation. This enhancer controls the expression of the pigmentation gene tan (t) in the abdominal epidermis. Two of the three SNPs were confirmed in an independent study using the D. melanogaster Genetic Reference Panel established from a North American population.

Results

We determined the functional impact of SNP1, SNP2, and SNP3 using transgenic lines to test all possible haplotypes in vivo. We show that all three candidate SNPs contribute to female Drosophila abdominal pigmentation. Interestingly, only two SNPs agree with the effect predicted by GWAS; the third one goes in the opposite direction because of linkage disequilibrium between multiple functional SNPs. Our experimental design uncovered strong additive effects for the three SNPs, but we also found significant epistatic effects explaining up to 11% of the total variation.

Conclusions

Our results suggest that linked causal variants are important for the interpretation of GWAS and functional validation is needed to understand the genetic architecture of traits.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1262-7) contains supplementary material, which is available to authorized users.

The Genomic Architecture of Flowering Time Varies Across Space and Time in Mimulus guttatus

The degree to which genomic architecture varies across space and time is central to the evolution of genomes in response to natural selection. Bulked-segregant mapping combined with pooled sequencing provides an efficient means to estimate the effect of genetic variants on quantitative traits. We develop a novel likelihood framework to identify segregating variation within multiple populations and generations while accommodating estimation error on a sample- and SNP-specific basis. We use this method to map loci for flowering time within natural populations of Mimulus guttatus, collecting the early- and late-flowering plants from each of three neighboring populations and two consecutive generations. Structural variants, such as inversions, and genes from multiple flowering-time pathways exhibit the strongest associations with flowering time. We find appreciable variation in genetic effects on flowering time across both time and space; the greatest differences evident between populations, where numerous factors (environmental variation, genomic background, and private polymorphisms) likely contribute to heterogeneity. However, the changes across years within populations clearly identify genotype-by-environment interactions as an important influence on flowering time variation.

A Variable Genetic Architecture of Melanic Evolution in Drosophila melanogaster

Unraveling the genetic architecture of adaptive phenotypic divergence is a fundamental quest in evolutionary biology. In Drosophila melanogaster, high-altitude melanism has evolved in separate mountain ranges in sub-Saharan Africa, potentially as an adaptation to UV intensity. We investigated the genetic basis of this melanism in three populations using a new bulk segregant analysis mapping method. We identified 19 distinct QTL regions from nine mapping crosses, with several QTL peaks overlapping between two or all populations, and yet different crosses involving the same melanic population commonly yielded distinct QTL. The strongest QTL often overlapped well-known pigmentation genes, but we typically did not find wide signals of genetic differentiation (F_ST) between lightly and darkly pigmented populations at these genes. Instead, we found small numbers of highly differentiated SNPs at the probable causative genes. A simulation analysis showed that these patterns of polymorphism were consistent with selection on standing genetic variation. Overall, our results suggest that, even for potentially simpler traits like pigmentation, the complexity of adaptive trait evolution poses important challenges for QTL mapping and population genetic analysis.

Ancestral population reconstitution from isofemale lines as a tool for experimental evolution

Experimental evolution is a powerful tool to study adaptation under controlled conditions. Laboratory natural selection experiments mimic adaptation in the wild with better‐adapted genotypes having more offspring. Because the selected traits are frequently not known, adaptation is typically measured as fitness increase by comparing evolved populations against an unselected reference population maintained in a laboratory environment. With adaptation to the laboratory conditions and genetic drift, however, it is not clear to what extent such comparisons provide unbiased estimates of adaptation. Alternatively, ancestral variation could be preserved in isofemale lines that can be combined to reconstitute the ancestral population. Here, we assess the impact of selection on alleles segregating in newly established Drosophila isofemale lines. We reconstituted two populations from isofemale lines and compared them to two original ancestral populations (AP) founded from the same lines shortly after collection. No significant allele frequency changes could be detected between both AP and simulations showed that drift had a low impact compared to Pool‐Seq‐associated sampling effects. We conclude that laboratory selection on segregating variation in isofemale lines is too weak to have detectable effects, which validates ancestral population reconstitution from isofemale lines as an unbiased approach for measuring adaptation in evolved populations.

The Genetic Basis of Pigmentation Differences Within and Between Drosophila Species

In Drosophila, as well as in many other plants and animals, pigmentation is highly variable both within and between species. This variability, combined with powerful genetic and transgenic tools as well as knowledge of how pigment patterns are formed biochemically and developmentally, has made Drosophila pigmentation a premier system for investigating the genetic and molecular mechanisms responsible for phenotypic evolution. In this chapter, we review and synthesize findings from a rapidly growing body of case studies examining the genetic basis of pigmentation differences in the abdomen, thorax, wings, and pupal cases within and between Drosophila species. A core set of genes, including genes required for pigment synthesis (eg, yellow, ebony, tan, Dat) as well as developmental regulators of these genes (eg, bab1, bab2, omb, Dll, and wg), emerge as the primary sources of this variation, with most genes having been shown to contribute to pigmentation differences both within and between species. In cases where specific genetic changes contributing to pigmentation divergence were identified in these genes, the changes were always located in noncoding sequences and affected cis-regulatory activity. We conclude this chapter by discussing these and other lessons learned from evolutionary genetic studies of Drosophila pigmentation and identify topics we think should be the focus of future work with this model system.