A genome-wide, fine-scale map of natural pigmentation variation in Drosophila melanogaster

Molecular genetic foundation of a sex-linked tailless trait in Hongshan chicken by whole genome data analysis

As a Chinese local chicken breed, Hongshan chickens have 2 kinds of tail feather phenotypes, normal and taillessness. Our previous studies showed that taillessness was a sex-linked dominant trait. Abnormal development of the tail vertebrae could be explained this phenomenon in some chicken breeds. However, the number of caudal vertebrae in rumpless Hongshan chickens was normal, so rumplessness in Hongshan chicken was not related to the development of the caudal vertebrae. Afterwards, we found that rumplessness in Hongshan was due to abnormal development of tail feather rather than abnormal development of caudal vertebrae. In order to understand the genetic foundation of the rumplessness of Hongshan chickens, we compared and reanalyzed 2 sets of data in normal and rumpless Hongshan chickens from our previous studies. By joint analysis of genome-wide selection signature analysis and genome-wide association approach, we found that 1 overlapping gene (EDIL3) and 16 peak genes (ENSGALG00000051843, ENSGALG00000053498, ENSGALG00000054800, KIF27, PTPRD, ENSGALG00000047579, ENSGALG00000041052, ARHGEF28, CAMK4, SERINC5, ENSGALG00000050776, ERCC8, MCC, ADAMTS19, ENSGALG00000053322, CHRNA8) located on the Z chromosome was associated with the rumpless trait. The results of this study furtherly revealed the molecular mechanism of the rumpless trait in Hongshan chickens, and identified the candidate genes associated with this trait. Our results will help to improve the shape of chicken tail feathers and to rise individual economic value in some specific market in China.

Identification of a Unique Genomic Region in Sweet Chestnut (Castanea sativa Mill.) That Controls Resistance to Asian Chestnut Gall Wasp Dryocosmus kuriphilus Yasumatsu

The Asian chestnut gall wasp (ACGW) (Hymenoptera Dryocosmus kuriphilus Yasumatsu) is a severe pest of sweet chestnut (Castanea sativa Mill.) with a strong impact on growth and nut production. A comparative field trial in Central Italy, including provenances from Spain, Italy, and Greece, was screened for ACGW infestation over consecutive years. The Greek provenance Hortiatis expressed a high proportion of immune plants and was used to perform a genome-wide association study based on DNA pool sequencing (Pool-GWAS) by comparing two DNA pools from 25 susceptible and 25 resistant plants. DNA pools were sequenced with 50X coverage depth. Sequence reads were aligned to a C. mollissima reference genome and the pools were compared to identify SNPs associated with resistance. Twenty-one significant SNPs were identified and highlighted a small genomic region on pseudochromosome 3 (Chr 3), containing 12 candidate genes of three gene families: Cytochrome P450, UDP-glycosyltransferase, and Rac-like GTP-binding protein. Functional analyses revealed a putative metabolic gene cluster related to saccharide biosynthesis in the genomic regions associated with resistance that could be involved in the production of a toxic metabolite against parasites. The comparison with previous genetic studies confirmed the involvement of Chr 3 in the control of resistance to ACGW.

Genetic basis of variation in thermal developmental plasticity for Drosophila melanogaster body pigmentation

Seasonal differences in insect pigmentation are attributed to the influence of ambient temperature on pigmentation development. This thermal plasticity is adaptive and heritable, and thereby capable of evolving. However, the specific genes contributing to the variation in plasticity that can drive its evolution remain largely unknown. To address this, we analysed pigmentation and pigmentation plasticity in Drosophila melanogaster. We measured two components of pigmentation in the thorax and abdomen: overall darkness and the proportion of length covered by darker pattern elements (a trident in the thorax and bands in the abdomen) in females from two developmental temperatures (17 or 28°C) and 191 genotypes. Using a GWAS approach to identify the genetic basis of variation in pigmentation and its response to temperature, we identified numerous dispersed QTLs, including some mapping to melanogenesis genes (yellow, ebony, and tan). Remarkably, we observed limited overlap between QTLs for variation within specific temperatures and those influencing thermal plasticity, as well as minimal overlap between plasticity QTLs across pigmentation components and across body parts. For most traits, consistent with selection favouring the retention of plasticity, we found that lower plasticity alleles were often at lower frequencies. The functional analysis of selected candidate QTLs and pigmentation genes largely confirmed their contributions to variation in pigmentation and/or pigmentation plasticity. Overall, our study reveals the existence and underlying basis of extensive and trait-specific genetic variation for pigmentation and pigmentation plasticity, offering a rich reservoir of raw material for natural selection to shape the evolution of these traits independently.

Regulatory evolution tuning pigmentation intensity quantitatively in Drosophila

Quantitative variation in attributes such as color, texture, or stiffness dominates morphological diversification. It results from combinations of alleles at many Mendelian loci. Here, we identify an additional source of quantitative variation among species, continuous evolution in a gene regulatory region. Specifically, we examined the modulation of wing pigmentation in a group of fly species and showed that inter-species variation correlated with the quantitative expression of the pigmentation gene yellow. This variation results from an enhancer of yellow determining darkness through species-specific activity. We mapped the divergent activities between two sister species and found the changes to be broadly distributed along the enhancer. Our results demonstrate that enhancers can act as dials fueling quantitative morphological diversification by modulating trait properties.

The Ancestral Modulation Hypothesis: Predicting Mechanistic Control of Sexually Heteromorphic Traits Using Evolutionary History

AbstractAcross the animal kingdom there are myriad forms within a sex across, and even within, species, rendering concepts of universal sex traits moot. The mechanisms that regulate the development of these trait differences are varied, although in vertebrates, common pathways involve gonadal steroid hormones. Gonadal steroids are often associated with heteromorphic trait development, where the steroid found at higher circulating levels is the one involved in trait development for that sex. Occasionally, there are situations in which a gonadal steroid associated with heteromorphic trait development in one sex is involved in heteromorphic or monomorphic trait development in another sex. We propose a verbal hypothesis, the ancestral modulation hypothesis (AMH), that uses the evolutionary history of the trait-particularly which sex ancestrally possessed higher trait values-to predict the regulatory pathway that governs trait expression. The AMH predicts that the genomic architecture appears first to resolve sexual conflict in an initially monomorphic trait. This architecture takes advantage of existing sex-biased signals, the gonadal steroid pathway, to generate trait heteromorphism. In cases where the other sex experiences evolutionary pressure for the new phenotype, that sex will co-opt the existing architecture by altering its signal to match that of the original high-trait-value sex. We describe the integrated levels needed to produce this pattern and what the expected outcomes will be given the evolutionary history of the trait. We present this framework as a testable hypothesis for the scientific community to investigate and to create further engagement and analysis of both ultimate and proximate approaches to sexual heteromorphism.

Complexities of recapitulating polygenic effects in natural populations: replication of genetic effects on wing shape in artificially selected and wild-caught populations of Drosophila melanogaster

Identifying the genetic architecture of complex traits is important to many geneticists, including those interested in human disease, plant and animal breeding, and evolutionary genetics. Advances in sequencing technology and statistical methods for genome-wide association studies have allowed for the identification of more variants with smaller effect sizes, however, many of these identified polymorphisms fail to be replicated in subsequent studies. In addition to sampling variation, this failure to replicate reflects the complexities introduced by factors including environmental variation, genetic background, and differences in allele frequencies among populations. Using Drosophila melanogaster wing shape, we ask if we can replicate allelic effects of polymorphisms first identified in a genome-wide association studies in three genes: dachsous, extra-macrochaete, and neuralized, using artificial selection in the lab, and bulk segregant mapping in natural populations. We demonstrate that multivariate wing shape changes associated with these genes are aligned with major axes of phenotypic and genetic variation in natural populations. Following seven generations of artificial selection along the dachsous shape change vector, we observe genetic differentiation of variants in dachsous and genomic regions containing other genes in the hippo signaling pathway. This suggests a shared direction of effects within a developmental network. We also performed artificial selection with the extra-macrochaete shape change vector, which is not a part of the hippo signaling network, but showed a largely shared direction of effects. The response to selection along the emc vector was similar to that of dachsous, suggesting that the available genetic diversity of a population, summarized by the genetic (co)variance matrix (G), influenced alleles captured by selection. Despite the success with artificial selection, bulk segregant analysis using natural populations did not detect these same variants, likely due to the contribution of environmental variation and low minor allele frequencies, coupled with small effect sizes of the contributing variants.

Drosophilids with darker cuticle have higher body temperature under light

Cuticle pigmentation was shown to be associated with body temperature for several relatively large species of insects, but it was questioned for small insects. Here we used a thermal camera to assess the association between drosophilid cuticle pigmentation and body temperature increase when individuals are exposed to light. We compared mutants of large effects within species (Drosophila melanogaster ebony and yellow mutants). Then we analyzed the impact of naturally occurring pigmentation variation within species complexes (Drosophila americana/Drosophila novamexicana and Drosophila yakuba/Drosophila santomea). Finally we analyzed lines of D. melanogaster with moderate differences in pigmentation. We found significant differences in temperatures for each of the four pairs we analyzed. The temperature differences appeared to be proportional to the differently pigmented area: between Drosophila melanogaster ebony and yellow mutants or between Drosophila americana and Drosophila novamexicana, for which the whole body is differently pigmented, the temperature difference was around 0.6 °C ± 0.2 °C. By contrast, between D. yakuba and D. santomea or between Drosophila melanogaster Dark and Pale lines, for which only the posterior abdomen is differentially pigmented, we detected a temperature difference of about 0.14 °C ± 0.10 °C. This strongly suggests that cuticle pigmentation has ecological implications in drosophilids regarding adaptation to environmental temperature.

Identifying genes associated with genetic control of color polymorphism in the pearl oyster Pinctada margaritifera var. cumingii (Linnaeus 1758) using a comparative whole genome pool-sequencing approach

For hundreds of years, the color diversity of Mollusca shells has been a topic of interest for humanity. However, the genetic control underlying color expression is still poorly understood in mollusks. The pearl oyster Pinctada margaritifera is increasingly becoming a biological model to study this process due to its ability to produce a large range of colors. Previous breeding experiments demonstrated that color phenotypes were partly under genetic control, and while a few genes were found in comparative transcriptomics and epigenetic experiments, genetic variants associated with the phenotypes have not yet been investigated. Here, we used a pooled-sequencing approach on 172 individuals to investigate color-associated variants on three color phenotypes of economic interest for pearl farming, in three wild and one hatchery populations. While our results uncovered SNPs targeting pigment-related genes already identified in previous studies, such as PBGD, tyrosinases, GST, or FECH, we also identified new color-related genes occurring in the same pathways, like CYP4F8, CYP3A4, and CYP2R1. Moreover, we identified new genes involved in novel pathways unknown to be involved in shell coloration for P. margaritifera, like the carotenoid pathway, BCO1. These findings are essential to possibly implement future breeding programs focused on individual selection for specific color production in pearl oysters and improve the footprint of perliculture on the Polynesian lagoon by producing less but with a better quality.

Bait-ER: A Bayesian method to detect targets of selection in Evolve-and-Resequence experiments

For over a decade, experimental evolution has been combined with high-throughput sequencing techniques. In so-called Evolve-and-Resequence (E&R) experiments, populations are kept in the laboratory under controlled experimental conditions where their genomes are sampled and allele frequencies monitored. However, identifying signatures of adaptation in E&R datasets is far from trivial, and it is still necessary to develop more efficient and statistically sound methods for detecting selection in genome-wide data. Here, we present Bait-ER - a fully Bayesian approach based on the Moran model of allele evolution to estimate selection coefficients from E&R experiments. The model has overlapping generations, a feature that describes several experimental designs found in the literature. We tested our method under several different demographic and experimental conditions to assess its accuracy and precision, and it performs well in most scenarios. Nevertheless, some care must be taken when analysing trajectories where drift largely dominates and starting frequencies are low. We compare our method with other available software and report that ours has generally high accuracy even for trajectories whose complexity goes beyond a classical sweep model. Furthermore, our approach avoids the computational burden of simulating an empirical null distribution, outperforming available software in terms of computational time and facilitating its use on genome-wide data. We implemented and released our method in a new open-source software package that can be accessed at https://doi.org/10.5281/zenodo.7351736.

Evolution and genomic basis of the plant-penetrating ovipositor: a key morphological trait in herbivorous Drosophilidae

Herbivorous insects are extraordinarily diverse, yet are found in only one-third of insect orders. This skew may result from barriers to plant colonization, coupled with phylogenetic constraint on plant-colonizing adaptations. The plant-penetrating ovipositor, however, is one trait that surmounts host plant physical defences and may be evolutionarily labile. Ovipositors densely lined with hard bristles have evolved repeatedly in herbivorous lineages, including within the Drosophilidae. However, the evolution and genetic basis of this innovation has not been well studied. Here, we focused on the evolution of this trait in Scaptomyza, a genus sister to Hawaiian Drosophila, that contains a herbivorous clade. Our phylogenetic approach revealed that ovipositor bristle number increased as herbivory evolved in the Scaptomyza lineage. Through a genome-wide association study, we then dissected the genomic architecture of variation in ovipositor bristle number within S. flava. Top-associated variants were enriched for transcriptional repressors, and the strongest associations included genes contributing to peripheral nervous system development. Individual genotyping supported the association at a variant upstream of Gαi, a neural development gene, contributing to a gain of 0.58 bristles/major allele. These results suggest that regulatory variation involving conserved developmental genes contributes to this key morphological trait involved in plant colonization.

Transgenerational effects of gamma radiation dose and dose rate on Drosophila flies irradiated at an early embryonal stage

Ionizing radiation (IR) kills cells mainly through induction of DNA damages and the surviving cells may suffer from mutations. Transgenerational effects of IR are well documented, but the exact mechanisms underlying them are less well understood; they include induction of mutations in germ cells and epigenetic inheritance. Previously, effects in the offspring of mice and zebrafish exposed to IR have been reported. A few studies also showed indications of transgenerational effects of radiation in humans, particularly in nuclear power workers. In the present project, short- and long-term effects of low-dose-rate (LDR; 50 and 97 mGy/h) and high-dose-rate (HDR; 23.4, 47.1 and 495 Gy/h) IR in Drosophila embryos were investigated. The embryos were irradiated at different doses and dose rates and radiosensitivity at different developmental stages was investigated. Also, the survival of larvae, pupae and adults developed from embryos irradiated at an early stage (30 min after egg laying) were studied. The larval crawling and pupation height assays were applied to investigate radiation effects on larval locomotion and pupation behavior, respectively. In parallel, the offspring from 3 Gy irradiated early-stage embryos were followed up to 12 generations and abnormal phenotypes were studied. Acute exposure of embryos at different stages of development showed that the early stage embryo is the most sensitive. The effects on larval locomotion showed no significant differences between the dose rates but a significant decrease of locomotion activity above 7 Gy was observed. The results indicate that embryos exposed to the low dose rates have shorter eclosion times. At the same cumulative dose (1 up to 7 Gy), HDR is more embryotoxic than LDR. We also found a radiation-induced depigmentation on males (A5 segment of the dorsal abdomen, A5pig^-) that can be transmitted up to 12 generations. The phenomenon does not follow the classical Mendelian laws of segregation.

Multi-omics analysis identifies a CYP9K1 haplotype conferring pyrethroid resistance in the malaria vector Anopheles funestus in East Africa

Metabolic resistance to pyrethroids is a menace to the continued effectiveness of malaria vector controls. Its molecular basis is complex and varies geographically across Africa. Here, we used a multi‐omics approach, followed‐up with functional validation to show that a directionally selected haplotype of a cytochrome P450, CYP9K1 is a major driver of resistance in Anopheles funestus. A PoolSeq GWAS using mosquitoes alive and dead after permethrin exposure, from Malawi and Cameroon, detected candidate genomic regions, but lacked consistency across replicates. Targeted sequencing of candidate resistance genes detected several SNPs associated with known pyrethroid resistance QTLs. The most significant SNPs were in the cytochrome P450 CYP304B1 (Cameroon), CYP315A1 (Uganda) and the ABC transporter gene ABCG4 (Malawi). However, when comparing field resistant mosquitoes to laboratory susceptible, the pyrethroid resistance locus rp1 and SNPs around the ABC transporter ABCG4 were consistently significant, except for Uganda where SNPs in the P450 CYP9K1 was markedly significant. In vitro heterologous metabolism assays with recombinant CYP9K1 revealed that it metabolises type II pyrethroid (deltamethrin; 64% depletion) but not type I (permethrin; 0%), while moderately metabolising DDT (17%). CYP9K1 exhibited reduced genetic diversity in Uganda underlying an extensive selective sweep. Furthermore, a glycine to alanine (G454A) amino acid change in CYP9K1 was fixed in Ugandan mosquitoes but not in other An. funestus populations. This study sheds further light on the evolution of metabolic resistance in a major malaria vector by implicating more genes and variants that can be used to design field‐applicable markers to better track resistance Africa‐wide.

The role of the epidermis enhancer element in positive and negative transcriptional regulation of ebony in Drosophila melanogaster

The spatiotemporal regulation of gene expression is essential to ensure robust phenotypic outcomes. Pigmentation patterns in Drosophila are determined by pigments biosynthesized in the developing epidermis and the cis-regulatory elements of the genes involved in this process are well-characterized. Here, we report that the known primary epidermal enhancer is dispensable for the transcriptional activation of ebony (involved in light-colored pigment synthesis) in the developing epidermis of Drosophila melanogaster. The evidence was obtained by introducing an approximately 1 kbp deletion at the primary epidermal enhancer by genome editing. The effect of the primary epidermal enhancer deletion on pigmentation and on the endogenous expression pattern of a mCherry-fused ebony allele was examined in the abdomen. The expression levels of the mCherry-fused ebony in the primary epidermal enhancer-deleted strains were slightly higher than that of the control strain, indicating that the sequences outside the primary epidermal enhancer have an ability to drive an expression of this gene in the epidermis. Interestingly, the primary epidermal enhancer deletion resulted in a derepression of this gene in the dorsal midline of the abdominal tergites, where dark pigmentation is present in the wild-type individuals. This indicated that the primary epidermal enhancer fragment contains a silencer. Furthermore, the endogenous expression pattern of ebony in the 2 additional strains with partially deleted primary epidermal enhancer revealed that the silencer resides within a 351-bp fragment in the 5' portion of the primary epidermal enhancer. These results demonstrated that deletion assays combined with reporter assays are highly effective in detecting the presence of positively and negatively regulating sequences within and outside the focal cis-regulatory elements.

Pool-GWAS on reproductive dormancy in Drosophila simulans suggests a polygenic architecture

The genetic basis of adaptation to different environments has been of long-standing interest to evolutionary biologists. Dormancy is a well-studied adaptation to facilitate overwintering. In Drosophila melanogaster, a moderate number of genes with large effects have been described, which suggests a simple genetic basis of dormancy. On the other hand, genome-wide scans for dormancy suggest a polygenic architecture in insects. In D. melanogaster, the analysis of the genetic architecture of dormancy is complicated by the presence of cosmopolitan inversions. Here, we performed a genome-wide scan to characterize the genetic basis of this ecologically extremely important trait in the sibling species of D. melanogaster, D. simulans that lacks cosmopolitan inversions. We performed Pool-GWAS in a South African D. simulans population for dormancy incidence at 2 temperature regimes (10 and 12°C, LD 10:14). We identified several genes with SNPs that showed a significant association with dormancy (P-value < 1e-13), but the overall modest response suggests that dormancy is a polygenic trait with many loci of small effect. Our results shed light on controversies on reproductive dormancy in Drosophila and have important implications for the characterization of the genetic basis of this trait.

Powerful, efficient QTL mapping in Drosophila melanogaster using bulked phenotyping and pooled sequencing

Despite the value of recombinant inbred lines for the dissection of complex traits, large panels can be difficult to maintain, distribute, and phenotype. An attractive alternative to recombinant inbred lines for many traits leverages selecting phenotypically extreme individuals from a segregating population, and subjecting pools of selected and control individuals to sequencing. Under a bulked or extreme segregant analysis paradigm, genomic regions contributing to trait variation are revealed as frequency differences between pools. Here, we describe such an extreme quantitative trait locus, or extreme quantitative trait loci, mapping strategy that builds on an existing multiparental population, the Drosophila Synthetic Population Resource, and involves phenotyping and genotyping a population derived by mixing hundreds of Drosophila Synthetic Population Resource recombinant inbred lines. Simulations demonstrate that challenging, yet experimentally tractable extreme quantitative trait loci designs (≥4 replicates, ≥5,000 individuals/replicate, and selecting the 5-10% most extreme animals) yield at least the same power as traditional recombinant inbred line-based quantitative trait loci mapping and can localize variants with sub-centimorgan resolution. We empirically demonstrate the effectiveness of the approach using a 4-fold replicated extreme quantitative trait loci experiment that identifies 7 quantitative trait loci for caffeine resistance. Two mapped extreme quantitative trait loci factors replicate loci previously identified in recombinant inbred lines, 6/7 are associated with excellent candidate genes, and RNAi knock-downs support the involvement of 4 genes in the genetic control of trait variation. For many traits of interest to drosophilists, a bulked phenotyping/genotyping extreme quantitative trait loci design has considerable advantages.

Bulk segregation analysis in the NGS era: a review of its teenage years

Bulk segregation analysis (BSA) utilizes a strategy of pooling individuals with extreme phenotypes to conduct economical and rapidly linked marker screening or quantitative trait locus (QTL) mapping. With the development of next-generation sequencing (NGS) technology in the past 10 years, BSA methods and technical systems have been gradually developed and improved. At the same time, the ever-decreasing costs of sequencing accelerate NGS-based BSA application in different species, including eukaryotic yeast, grain crops, economic crops, horticultural crops, trees, aquatic animals, and insects. This paper provides a landscape of BSA methods and reviews the BSA development process in the past decade, including the sequencing method for BSA, different populations, different mapping algorithms, associated region threshold determination, and factors affecting BSA mapping. Finally, we summarize related strategies in QTL fine mapping combining BSA.

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.

Population Structure and Diversity in European Honey Bees (Apis mellifera L.)—An Empirical Comparison of Pool and Individual Whole-Genome Sequencing

Background: Whole-genome sequencing has become routine for population genetic studies. Sequencing of individuals provides maximal data but is rather expensive and fewer samples can be studied. In contrast, sequencing a pool of samples (pool-seq) can provide sufficient data, while presenting less of an economic challenge. Few studies have compared the two approaches to infer population genetic structure and diversity in real datasets. Here, we apply individual sequencing (ind-seq) and pool-seq to the study of Western honey bees (Apis mellifera). Methods: We collected honey bee workers that belonged to 14 populations, including 13 subspecies, totaling 1347 colonies, who were individually (139 individuals) and pool-sequenced (14 pools). We compared allele frequencies, genetic diversity estimates, and population structure as inferred by the two approaches. Results: Pool-seq and ind-seq revealed near identical population structure and genetic diversities, albeit at different costs. While pool-seq provides genome-wide polymorphism data at considerably lower costs, ind-seq can provide additional information, including the identification of population substructures, hybridization, or individual outliers. Conclusions: If costs are not the limiting factor, we recommend using ind-seq, as population genetic structure can be inferred similarly well, with the advantage gained from individual genetic information. Not least, it also significantly reduces the effort required for the collection of numerous samples and their further processing in the laboratory.

Drosophila Evolution over Space and Time (DEST): A New Population Genomics Resource

Drosophila melanogaster is a leading model in population genetics and genomics, and a growing number of whole-genome data sets from natural populations of this species have been published over the last years. A major challenge is the integration of disparate data sets, often generated using different sequencing technologies and bioinformatic pipelines, which hampers our ability to address questions about the evolution of this species. Here we address these issues by developing a bioinformatics pipeline that maps pooled sequencing (Pool-Seq) reads from D. melanogaster to a hologenome consisting of fly and symbiont genomes and estimates allele frequencies using either a heuristic (PoolSNP) or a probabilistic variant caller (SNAPE-pooled). We use this pipeline to generate the largest data repository of genomic data available for D. melanogaster to date, encompassing 271 previously published and unpublished population samples from over 100 locations in >20 countries on four continents. Several of these locations have been sampled at different seasons across multiple years. This data set, which we call Drosophila Evolution over Space and Time (DEST), is coupled with sampling and environmental metadata. A web-based genome browser and web portal provide easy access to the SNP data set. We further provide guidelines on how to use Pool-Seq data for model-based demographic inference. Our aim is to provide this scalable platform as a community resource which can be easily extended via future efforts for an even more extensive cosmopolitan data set. Our resource will enable population geneticists to analyze spatiotemporal genetic patterns and evolutionary dynamics of D. melanogaster populations in unprecedented detail.

Next-Generation Sequencing and Its Impacts on Entomological Research in Ecology and Evolution

The advent of NGS-based methods has been profoundly transforming entomological research. Through continual development and improvement of different methods and sequencing platforms, NGS has promoted mass elucidation of partial or whole genetic materials associated with beneficial insects, pests (of agriculture, forestry and animal, and human health), and species of conservation concern, helping to unravel ecological and evolutionary mechanisms and characterizing survival, trophic interactions, and dispersal. It is shifting the scale of biodiversity and environmental analyses from individuals and biodiversity indicator species to the large-scale study of communities and ecosystems using bulk samples of species or a mixed "soup" of environmental DNA. As the NGS-based methods have become more affordable, complexity demystified, and specificity and sensitivity proven, their use in entomological research has spread widely. This article presents several examples on how NGS-based methods have been used in entomology to provide incentives to apply them when appropriate and to open our minds to the expected advances in entomology that are yet to come.

Optimizing sampling design and sequencing strategy for the genomic analysis of quantitative traits in natural populations

Mapping the genes underlying ecologically relevant traits in natural populations is fundamental to develop a molecular understanding of species adaptation. Current sequencing technologies enable the characterization of a species' genetic diversity across the landscape or even over its whole range. The relevant capture of the genetic diversity across the landscape is critical for a successful genetic mapping of traits and there are no clear guidelines on how to achieve an optimal sampling and which sequencing strategy to implement. Here we determine, through simulation, the sampling scheme that maximizes the power to map the genetic basis of a complex trait in an outbreeding species across an idealized landscape and draw genomic predictions for the trait, comparing individual and pool sequencing strategies. Our results show that quantitative trait locus detection power and prediction accuracy are higher when more populations over the landscape are sampled and this is more cost-effectively done with pool sequencing than with individual sequencing. Additionally, we recommend sampling populations from areas of high genetic diversity. As progress in sequencing enables the integration of trait-based functional ecology into landscape genomics studies, these findings will guide study designs allowing direct measures of genetic effects in natural populations across the environment.

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.

bric à brac controls sex pheromone choice by male European corn borer moths

The sex pheromone system of ~160,000 moth species acts as a powerful form of assortative mating whereby females attract conspecific males with a species-specific blend of volatile compounds. Understanding how female pheromone production and male preference coevolve to produce this diversity requires knowledge of the genes underlying change in both traits. In the European corn borer moth, pheromone blend variation is controlled by two alleles of an autosomal fatty-acyl reductase gene expressed in the female pheromone gland (pgFAR). Here we show that asymmetric male preference is controlled by cis-acting variation in a sex-linked transcription factor expressed in the developing male antenna, bric à brac (bab). A genome-wide association study of preference using pheromone-trapped males implicates variation in the 293 kb bab intron 1, rather than the coding sequence. Linkage disequilibrium between bab intron 1 and pgFAR further validates bab as the preference locus, and demonstrates that the two genes interact to contribute to assortative mating. Thus, lack of physical linkage is not a constraint for coevolutionary divergence of female pheromone production and male behavioral response genes, in contrast to what is often predicted by evolutionary theory.

Identification of major QTL for waterlogging tolerance in maize using genome-wide association study and bulked sample analysis

Waterlogging has increasingly become one of the major constraints to maize (Zea mays L.) production in some maize growing areas as it seriously decreases the yield. Waterlogging tolerance in maize germplasm provides a basis for maize waterlogging improvement. In this study, nine seedling traits, plant height (PH), root length (RL), shoot dry weight (SDW), root dry weight (RDW), adventitious root number (ARN), node number of brace root (BRNN), brace root number (BRN), brace root dry weigh (BRDW), survival rate (SR), and the secondary traits that were defined as relative phenotypic value of seedling traits under waterlogging and control treatments were used in a natural population that contain 365 inbred lines to evaluate the waterlogging tolerance of tropical maize. The result showed that maize waterlogging tolerance was genetically controlled and seedling traits were significantly different between the control and waterlogging treatments. PH, RL, SDW, and RDW are important seedling traits for waterlogging tolerance identification. Some tropical maize inbred lines were identified with extreme waterlogging tolerance that can provide an important germplasm resource for breeding. Population structure analysis showed that two major phylogenetic subgroups in tropical maize could be identified. Genome-wide association study (GWAS) using 39,266 single nucleotide polymorphisms (SNPs) across the whole genome identified 49 trait-SNPs distributed on over all 10 chromosomes excluding chromosome 10. Seventy-one significant SNPs, distributed on all 10 chromosomes excluding chromosome 5, were identified by extend bulked sample analysis (Ext-BSA) based on the inbred lines with extreme phenotypes. GWAS and Ext-BSA identified the same loci on bin1.07, bin6.01, bin2.09, bin6.04, bin7.02, and bin7.03. Nine genes were proposed as potential candidate genes. Cloning and functional validation of these genes would be helpful for understanding the molecular mechanism of waterlogging tolerance in maize.

Population Genomics on the Fly: Recent Advances in Drosophila

Drosophila melanogaster, a small dipteran of African origin, represents one of the best-studied model organisms. Early work in this system has uniquely shed light on the basic principles of genetics and resulted in a versatile collection of genetic tools that allow to uncover mechanistic links between genotype and phenotype. Moreover, given its worldwide distribution in diverse habitats and its moderate genome-size, Drosophila has proven very powerful for population genetics inference and was one of the first eukaryotes whose genome was fully sequenced. In this book chapter, we provide a brief historical overview of research in Drosophila and then focus on recent advances during the genomic era. After describing different types and sources of genomic data, we discuss mechanisms of neutral evolution including the demographic history of Drosophila and the effects of recombination and biased gene conversion. Then, we review recent advances in detecting genome-wide signals of selection, such as soft and hard selective sweeps. We further provide a brief introduction to background selection, selection of noncoding DNA and codon usage and focus on the role of structural variants, such as transposable elements and chromosomal inversions, during the adaptive process. Finally, we discuss how genomic data helps to dissect neutral and adaptive evolutionary mechanisms that shape genetic and phenotypic variation in natural populations along environmental gradients. In summary, this book chapter serves as a starting point to Drosophila population genomics and provides an introduction to the system and an overview to data sources, important population genetic concepts and recent advances in the field.

T-lex3: an accurate tool to genotype and estimate population frequencies of transposable elements using the latest short-read whole genome sequencing data

Motivation

Transposable elements (TEs) constitute a significant proportion of the majority of genomes sequenced to date. TEs are responsible for a considerable fraction of the genetic variation within and among species. Accurate genotyping of TEs in genomes is therefore crucial for a complete identification of the genetic differences among individuals, populations and species.

Results

In this work, we present a new version of T-lex, a computational pipeline that accurately genotypes and estimates the population frequencies of reference TE insertions using short-read high-throughput sequencing data. In this new version, we have re-designed the T-lex algorithm to integrate the BWA-MEM short-read aligner, which is one of the most accurate short-read mappers and can be launched on longer short-reads (e.g. reads >150 bp). We have added new filtering steps to increase the accuracy of the genotyping, and new parameters that allow the user to control both the minimum and maximum number of reads, and the minimum number of strains to genotype a TE insertion. We also showed for the first time that T-lex3 provides accurate TE calls in a plant genome.

Availability and implementation

To test the accuracy of T-lex3, we called 1630 individual TE insertions in Drosophila melanogaster, 1600 individual TE insertions in humans, and 3067 individual TE insertions in the rice genome. We showed that this new version of T-lex is a broadly applicable and accurate tool for genotyping and estimating TE frequencies in organisms with different genome sizes and different TE contents. T-lex3 is available at Github: https://github.com/GonzalezLab/T-lex3.

Supplementary information

Supplementary data are available at Bioinformatics online.

High-Throughput Genome-Wide Genotyping To Optimize the Use of Natural Genetic Resources in the Grassland Species Perennial Ryegrass (Lolium perenne L.)

The natural genetic diversity of agricultural species is an essential genetic resource for breeding programs aiming to improve their ecosystem and production services. A large natural ecotype diversity is usually available for most grassland species. This could be used to recombine natural climatic adaptations and agronomic value to create improved populations of grassland species adapted to future regional climates. However describing natural genetic resources can be long and costly. Molecular markers may provide useful information to help this task. This opportunity was investigated for Lolium perenne L., using a set of 385 accessions from the natural diversity of this species collected right across Europe and provided by genebanks of several countries. For each of these populations, genotyping provided the allele frequencies of 189,781 SNP markers. GWAS were implemented for over 30 agronomic and/or putatively adaptive traits recorded in three climatically contrasted locations (France, Belgium, Germany). Significant associations were detected for hundreds of markers despite a strong confounding effect of the genetic background; most of them pertained to phenology traits. It is likely that genetic variability in these traits has had an important contribution to environmental adaptation and ecotype differentiation. Genomic prediction models calibrated using natural diversity were found to be highly effective to describe natural populations for almost all traits as well as commercial synthetic populations for some important traits such as disease resistance, spring growth or phenological traits. These results will certainly be valuable information to help the use of natural genetic resources of other species.

Next-generation biological control: the need for integrating genetics and genomics

Biological control is widely successful at controlling pests, but effective biocontrol agents are now more difficult to import from countries of origin due to more restrictive international trade laws (the Nagoya Protocol). Coupled with increasing demand, the efficacy of existing and new biocontrol agents needs to be improved with genetic and genomic approaches. Although they have been underutilised in the past, application of genetic and genomic techniques is becoming more feasible from both technological and economic perspectives. We review current methods and provide a framework for using them. First, it is necessary to identify which biocontrol trait to select and in what direction. Next, the genes or markers linked to these traits need be determined, including how to implement this information into a selective breeding program. Choosing a trait can be assisted by modelling to account for the proper agro‐ecological context, and by knowing which traits have sufficiently high heritability values. We provide guidelines for designing genomic strategies in biocontrol programs, which depend on the organism, budget, and desired objective. Genomic approaches start with genome sequencing and assembly. We provide a guide for deciding the most successful sequencing strategy for biocontrol agents. Gene discovery involves quantitative trait loci analyses, transcriptomic and proteomic studies, and gene editing. Improving biocontrol practices includes marker‐assisted selection, genomic selection and microbiome manipulation of biocontrol agents, and monitoring for genetic variation during rearing and post‐release. We conclude by identifying the most promising applications of genetic and genomic methods to improve biological control efficacy.

[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.

BSA-Seq Discovery and Functional Analysis of Candidate Hessian Fly (Mayetiola destructor) Avirulence Genes

The Hessian fly (HF, Mayetiola destructor) is a plant-galling parasite of wheat (Triticum spp.). Seven percent of its genome is composed of highly diversified signal-peptide-encoding genes that are transcribed in HF larval salivary glands. These observations suggest that they encode effector proteins that are injected into wheat cells to suppress basal wheat immunity and redirect wheat development towards gall formation. Genetic mapping has determined that mutations in four of these genes are associated with HF larval survival (virulence) on plants carrying four different resistance (R) genes. Here, this line of investigation was pursued further using bulked-segregant analysis combined with whole genome resequencing (BSA-seq). Virulence to wheat R genes H6, Hdic, and H5 was examined. Mutations associated with H6 virulence had been mapped previously. Therefore, we used H6 to test the capacity of BSA-seq to map virulence using a field-derived HF population. This was the first time a non-structured HF population had been used to map HF virulence. Hdic virulence had not been mapped previously. Using a structured laboratory population, BSA-seq associated Hdic virulence with mutations in two candidate effector-encoding genes. Using a laboratory population, H5 virulence was previously positioned in a region spanning the centromere of HF autosome 2. BSA-seq resolved H5 virulence to a 1.3 Mb fragment on the same chromosome but failed to identify candidate mutations. Map-based candidate effectors were then delivered to Nicotiana plant cells via the type III secretion system of Burkholderia glumae bacteria. These experiments demonstrated that the genes associated with virulence to wheat R genes H6 and H13 are capable of suppressing plant immunity. Results are consistent with the hypothesis that effector proteins underlie the ability of HFs to survive on wheat.

Complex Genetic System Involved in Fusarium Ear Rot Resistance in Maize as Revealed by GWAS, Bulked Sample Analysis, and Genomic Prediction

Fusarium ear rot (FER) caused by Fusarium verticillioides is one of the most prevalent maize diseases in China and worldwide. Resistance to FER is a complex trait controlled by multiple genes highly affected by environment. In this paper, genome-wide association study (GWAS), bulked sample analysis (BSA), and genomic prediction were performed for understanding FER resistance using 509 diverse inbred lines, which were genotyped by 37,801 high-quality single-nucleotide polymorphisms (SNPs). Ear rot evaluation was performed using artificial inoculation in four environments in China: Xinxiang, Henan, and Shunyi, Beijing, during 2017 and 2018. Significant phenotypic and genetic variation for FER severity was observed, and FER resistance was significantly correlated among the four environments with a generalized heritability of 0.78. GWAS identified 23 SNPs that were associated with FER resistance, 2 of which (1_226233417 on chromosome 1 and 10_14501044 on chromosome 10) were associated at threshold of 2.65 × 10^-7 [-log(0.01/37,801)]. Using BSA, resistance quantitative trait loci were identified on chromosomes 3, 4, 7, 9, and 10 at the 90% confidence level and on chromosomes 3 and 10 at the 95% confidence level. A key region, bin 10.03, was detected by both GWAS and BSA. Genomic prediction for FER resistance showed that the prediction accuracy by trait-related markers was higher than that by randomly selected markers under different levels of marker density. Marker-assisted selection using genomic prediction could be an efficient strategy for genetic improvement for complex traits like FER resistance.

Sex determination in the GIFT strain of tilapia is controlled by a locus in linkage group 23

Background

Tilapias (Family Cichlidae) are the second most important group of aquaculture species in the world. They have been the subject of much research on sex determination due to problems caused by early maturation in culture and their complex sex-determining systems. Different sex-determining loci (linkage group 1, 20 and 23) have been detected in various tilapia stocks. The ‘genetically improved farmed tilapia’ (GIFT) stock, founded from multiple Nile tilapia (Oreochromis niloticus) populations, with some likely to have been introgressed with O. mossambicus, is a key resource for tilapia aquaculture. The sex-determining mechanism in the GIFT stock was unknown, but potentially complicated due to its multiple origins.

Results

A bulk segregant analysis (BSA) version of double-digest restriction-site associated DNA sequencing (BSA-ddRADseq) was developed and used to detect and position sex-linked single nucleotide polymorphism (SNP) markers in 19 families from the GIFT strain breeding nucleus and two Stirling families as controls (a single XY locus had been previously mapped to LG1 in the latter). About 1500 SNPs per family were detected across the genome. Phenotypic sex in Stirling families showed strong association with LG1, whereas only SNPs located in LG23 showed clear association with sex in the majority of the GIFT families. No other genomic regions linked to sex determination were apparent. This region was validated using a series of LG23-specific DNA markers (five SNPs with highest association to sex from this study, the LG23 sex-associated microsatellite UNH898 and ARO172, and the recently isolated amhy marker for individual fish (n = 284).

Conclusions

Perhaps surprisingly given its multiple origins, sex determination in the GIFT strain breeding nucleus was associated only with a locus in LG23. BSA-ddRADseq allowed cost-effective analysis of multiple families, strengthening this conclusion. This technique has potential to be applied to other complex traits. The sex-linked SNP markers identified will be useful for potential marker-assisted selection (MAS) to control sex-ratio in GIFT tilapia to suppress unwanted reproduction during growout.

Comparative study of population genomic approaches for mapping colony-level traits

Social insect colonies exhibit colony-level phenotypes such as social immunity and task coordination, which are produced by the individual phenotypes. Mapping the genetic basis of such phenotypes requires associating the colony-level phenotype with the genotypes in the colony. In this paper, we examine alternative approaches to DNA extraction, library construction, and sequencing for genome wide association studies (GWAS) of colony-level traits using a population sample of Cataglyphis niger ants. We evaluate the accuracy of allele frequency estimation from sequencing a pool of individuals (pool-seq) from each colony using either whole-genome sequencing or reduced representation genomic sequencing. Based on empirical measurement of the experimental noise in sequenced DNA pools, we show that reduced representation pool-seq is drastically less accurate than whole-genome pool-seq. Surprisingly, normalized pooling of samples did not result in greater accuracy than un-normalized pooling. Subsequently, we evaluate the power of the alternative approaches for detecting quantitative trait loci (QTL) of colony-level traits by using simulations that account for an environmental effect on the phenotype. Our results can inform experimental designs and enable optimizing the power of GWAS depending on budget, availability of samples and research goals. We conclude that for a given budget, sequencing un-normalized pools of individuals from each colony provides optimal QTL detection power.

Genomic mapping techniques are used to map phenotypes to genotypes. Mapping is of general interest in any biological system, including fundamental studies of biological traits, clinical studies of genetic predisposition to disease, and agro- and bio-technological studies of domesticated plants and animals. Typically, such studies associate phenotypic measurements of individuals with their genotypes. Here we evaluate methodological approaches for genomic mapping of phenotypes that are expressed at the level of a group rather than that of individuals. We demonstrate that genomic sequencing of a DNA pool from multiple samples provides increased statistical power within a limited budget. Our results facilitate more efficient use of resources in genomic mapping studies that investigate group-level phenotypes.

Trait mapping in diverse arthropods by bulked segregant analysis

Bulked segregant analysis (BSA) is a cross-based method for genetic mapping in sexually reproducing organisms. The method's use of bulked (pooled) samples markedly reduces the genotyping effort associated with traditional linkage mapping studies. Further, it can be applied to species with life histories or physical attributes (as for micro-insects) that render genetic mapping with other methods impractical. Recent studies in both insects and mites have revealed that advanced BSA experimental designs can resolve causal loci to narrow genomic intervals, facilitating follow-up investigations. As high-quality genomes become more widely available, BSA methods are poised to become an increasingly important tool for the rapid mapping of both monogenic and polygenic traits in diverse arthropod species.

Optimizing the Power to Identify the Genetic Basis of Complex Traits with Evolve and Resequence Studies

Evolve and resequence (E&R) studies are frequently used to dissect the genetic basis of quantitative traits. By subjecting a population to truncating selection for several generations and estimating the allele frequency differences between selected and nonselected populations using next-generation sequencing (NGS), the loci contributing to the selected trait may be identified. The role of different parameters, such as, the population size or the number of replicate populations has been examined in previous works. However, the influence of the selection regime, that is the strength of truncating selection during the experiment, remains little explored. Using whole genome, individual based forward simulations of E&R studies, we found that the power to identify the causative alleles may be maximized by gradually increasing the strength of truncating selection during the experiment. Notably, such an optimal selection regime comes at no or little additional cost in terms of sequencing effort and experimental time. Interestingly, we also found that a selection regime which optimizes the power to identify the causative loci is not necessarily identical to a regime that maximizes the phenotypic response. Finally, our simulations suggest that an E&R study with an optimized selection regime may have a higher power to identify the genetic basis of quantitative traits than a genome-wide association study, highlighting that E&R is a powerful approach for finding the loci underlying complex traits.

Altering the Temporal Regulation of One Transcription Factor Drives Evolutionary Trade-Offs between Head Sensory Organs

Size trade-offs of visual versus olfactory organs is a pervasive feature of animal evolution. This could result from genetic or functional constraints. We demonstrate that head sensory organ size trade-offs in Drosophila are genetically encoded and arise through differential subdivision of the head primordium into visual versus non-visual fields. We discover that changes in the temporal regulation of the highly conserved eyeless/Pax6 gene expression during development is a conserved mechanism for sensory trade-offs within and between Drosophila species. We identify a natural single nucleotide polymorphism in the cis-regulatory region of eyeless in a binding site of its repressor Cut that is sufficient to alter its temporal regulation and eye size. Because eyeless/Pax6 is a conserved regulator of head sensory placode subdivision, we propose that its temporal regulation is key to define the relative size of head sensory organs.

Benchmarking software tools for detecting and quantifying selection in evolve and resequencing studies

BACKGROUND

The combination of experimental evolution with whole-genome resequencing of pooled individuals, also called evolve and resequence (E&R) is a powerful approach to study the selection processes and to infer the architecture of adaptive variation. Given the large potential of this method, a range of software tools were developed to identify selected SNPs and to measure their selection coefficients.

RESULTS

In this benchmarking study, we compare 15 test statistics implemented in 10 software tools using three different scenarios. We demonstrate that the power of the methods differs among the scenarios, but some consistently outperform others. LRT-1, CLEAR, and the CMH test perform best despite LRT-1 and the CMH test not requiring time series data. CLEAR provides the most accurate estimates of selection coefficients.

CONCLUSION

This benchmark study will not only facilitate the analysis of already existing data, but also affect the design of future data collections.

Extreme-QTL mapping of monepantel resistance in Haemonchus contortus

Background

Haemonchus contortus, a gastrointestinal nematode parasite of sheep, is mainly controlled by anthelmintics; the occurrence of anthelmintic resistance leads to treatment failures and increases economic burden. Because molecular mechanisms involved in drug resistance can be elucidated by genomic studies, an extreme quantitative trait locus (X-QTL) mapping approach was used to identify co-segregation of the resistance phenotype with genetic markers to detect the genome-wide variants associated with monepantel resistance in H. contortus.

Methods

A cross between H. contortus isolates using parental susceptible (Par-S) males and monepantel resistant (Par-R) females resulted in SR progeny, while reciprocal cross resulted in RS progeny. Pools (n = 30,000) of infective larvae (L3) recovered from Par-R, and from SR and RS populations in the F3 generation, collected both before (unselected group) and 7 days after (selected group) selection with monepantel treatment in sheep hosts, were subjected to genome sequencing (Pool-Seq). Pairwise comparisons of allele frequencies between unselected and selected groups were performed for each population by Fisher’s exact test (FET) and for both populations combined by a Cochran-Mantel-Haenszel (CMH) test.

Results

Mapping rates varied from 80.29 to 81.77% at a 90.4X mean coverage of aligned reads. After correction for multiple testing, significant (P < 0.05) changes in allele frequencies were detected by FET for 6 and 57 single nucleotide polymorphisms (SNPs) in the SR and RS populations, respectively, and by the CMH test for 124 SNPs in both populations. The significant variants located on chromosome 2 generated a selection signal in a genomic region harboring the mptl-1, deg-3 and des-2 genes, previously reported as candidates for monepantel resistance. In addition, three new variants were identified in the mptl-1 gene.

Conclusions

This study expands knowledge on genome-wide molecular events underlying H. contortus resistance to monepantel. The identification of a genome region harboring major genes previously associated with monepantel resistance supports the results of the employed X-QTL approach. In addition, a deletion in exon 11 of the mptl-1 gene should be further investigated as the putative causal mutation leading to monepantel resistance.

Electronic supplementary material

The online version of this article (10.1186/s13071-019-3663-9) contains supplementary material, which is available to authorized users.

Comparing Pool-seq, Rapture, and GBS genotyping for inferring weak population structure: The American lobster (Homarus americanus) as a case study

Unraveling genetic population structure is challenging in species potentially characterized by large population size and high dispersal rates, often resulting in weak genetic differentiation. Genotyping a large number of samples can improve the detection of subtle genetic structure, but this may substantially increase sequencing cost and downstream bioinformatics computational time. To overcome this challenge, alternative, cost-effective sequencing approaches, namely Pool-seq and Rapture, have been developed. We empirically measured the power of resolution and congruence of these two methods in documenting weak population structure in nonmodel species with high gene flow comparatively to a conventional genotyping-by-sequencing (GBS) approach. For this, we used the American lobster (Homarus americanus) as a case study. First, we found that GBS, Rapture, and Pool-seq approaches gave similar allele frequency estimates (i.e., correlation coefficient over 0.90) and all three revealed the same weak pattern of population structure. Yet, Pool-seq data showed F ST estimates three to five times higher than GBS and Rapture, while the latter two methods returned similar F ST estimates, indicating that individual-based approaches provided more congruent results than Pool-seq. We conclude that despite higher costs, GBS and Rapture are more convenient approaches to use in the case of species exhibiting very weak differentiation. While both GBS and Rapture approaches provided similar results with regard to estimates of population genetic parameters, GBS remains more cost-effective in project involving a relatively small numbers of genotyped individuals (e.g., <1,000). Overall, this study illustrates the complexity of estimating genetic differentiation and other summary statistics in complex biological systems characterized by large population size and migration rates.

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.

Stress response, behavior, and development are shaped by transposable element-induced mutations in Drosophila

Most of the current knowledge on the genetic basis of adaptive evolution is based on the analysis of single nucleotide polymorphisms (SNPs). Despite increasing evidence for their causal role, the contribution of structural variants to adaptive evolution remains largely unexplored. In this work, we analyzed the population frequencies of 1,615 Transposable Element (TE) insertions annotated in the reference genome of Drosophila melanogaster, in 91 samples from 60 worldwide natural populations. We identified a set of 300 polymorphic TEs that are present at high population frequencies, and located in genomic regions with high recombination rate, where the efficiency of natural selection is high. The age and the length of these 300 TEs are consistent with relatively young and long insertions reaching high frequencies due to the action of positive selection. Besides, we identified a set of 21 fixed TEs also likely to be adaptive. Indeed, we, and others, found evidence of selection for 84 of these reference TE insertions. The analysis of the genes located nearby these 84 candidate adaptive insertions suggested that the functional response to selection is related with the GO categories of response to stimulus, behavior, and development. We further showed that a subset of the candidate adaptive TEs affects expression of nearby genes, and five of them have already been linked to an ecologically relevant phenotypic effect. Our results provide a more complete understanding of the genetic variation and the fitness-related traits relevant for adaptive evolution. Similar studies should help uncover the importance of TE-induced adaptive mutations in other species as well.

Evolutionary insights from large scale resequencing datasets in Drosophila melanogaster

Drosophila melanogaster has long been used as an evolutionary model system. Its small genome size, well-annotated genome, and ease of sampling, also makes it a choice species for genome resequencing studies. Hundreds of genomic samples from populations worldwide are available and are currently being used to tackle a wide range of evolutionary questions. In this review, we focused on three insights that have increased our understanding of the evolutionary history of this species, and that have implications for the study of evolutionary processes in other species as well. Because of technical limitations, most of the studies so far have focused on SNP variants. However, long-read sequencing techniques should allow us in the near future to include other type of genomic variants that also influence genome evolution.

Regulatory and sequence evolution in response to selection for improved associative learning ability in Nasonia vitripennis

Background

Selection acts on the phenotype, yet only the genotype is inherited. While both the phenotypic and genotypic response to short-term selection can be measured, the link between these is a major unsolved problem in evolutionary biology, in particular for complex behavioural phenotypes.

Results

Here we characterize the genomic and the transcriptomic basis of associative learning ability in the parasitic wasp Nasonia vitripennis and use gene network analysis to link the two. We artificially selected for improved associative learning ability in four independent pairs of lines and identified signatures of selection across the genome. Allele frequency diverged consistently between the selected and control lines in 118 single nucleotide polymorphisms (SNPs), clustering in 51 distinct genomic regions containing 128 genes. The majority of SNPs were found in regulatory regions, suggesting a potential role for gene expression evolution. We therefore sequenced the transcriptomes of selected and control lines and identified 36 consistently differentially expressed transcripts with large changes in expression. None of the differentially expressed genes also showed sequence divergence as a result of selection. Instead, gene network analysis showed many of the genes with consistent allele frequency differences and all of the differentially expressed genes to cluster in a single co-expression network. At a functional level, both genomic and transcriptomic analyses implicated members of gene networks known to be involved in neural plasticity and cognitive processes.

Conclusions

Taken together, our results reveal how specific cognitive abilities can readily respond to selection via a complex interplay between regulatory and sequence evolution.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5310-9) contains supplementary material, which is available to authorized users.

Association mapping desiccation resistance within chromosomal inversions in the African malaria vector Anopheles gambiae

Inversion polymorphisms are responsible for many ecologically important phenotypes and are often found under balancing selection. However, the same features that ensure their large role in local adaptation-especially reduced recombination between alternate arrangements-mean that uncovering the precise loci within inversions that control these phenotypes is unachievable using standard mapping approaches. Here, we take advantage of long-term balancing selection on a pair of inversions in the mosquito Anopheles gambiae to map desiccation tolerance via pool-GWAS. Two polymorphic inversions on chromosome 2 of this species (denoted 2La and 2Rb) are associated with arid and hot conditions in Africa and are maintained in spatially and temporally heterogeneous environments. After measuring thousands of wild-caught individuals for survival under desiccation stress, we used phenotypically extreme individuals homozygous for alternative arrangements at the 2La inversion to construct pools for whole-genome sequencing. Genomewide association mapping using these pools revealed dozens of significant SNPs within both 2La and 2Rb, many of which neighboured genes controlling ion channels or related functions. Our results point to the promise of similar approaches in systems with inversions maintained by balancing selection and provide a list of candidate genes underlying the specific phenotypes controlled by the two inversions studied here.

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.