Viola Nolte

PoPoolation: A Toolbox for Population Genetic Analysis of Next Generation Sequencing Data from Pooled Individuals

Recent statistical analyses suggest that sequencing of pooled samples provides a cost effective approach to determine genome-wide population genetic parameters. Here we introduce PoPoolation, a toolbox specifically designed for the population genetic analysis of sequence data from pooled individuals. PoPoolation calculates estimates of θ Watterson, θ π, and Tajimas D that account for the bias introduced by pooling and sequencing errors, as well as divergence between species. Results of genome-wide analyses can be graphically displayed in a sliding window plot. PoPoolation is written in Perl and R and it builds on commonly used data formats. Its source code can be downloaded from http://code.google.com/p/popoolation/. Furthermore, we evaluate the influence of mapping algorithms, sequencing errors, and read coverage on the accuracy of population genetic parameter estimates from pooled data.

Sequencing pools of individuals — mining genome-wide polymorphism data without big funding

The analysis of polymorphism data is becoming increasingly important as a complementary tool to classical genetic analyses. Nevertheless, despite plunging sequencing costs, genomic sequencing of individuals at the population scale is still restricted to a few model species. Whole-genome sequencing of pools of individuals (Pool-seq) provides a cost-effective alternative to sequencing individuals separately. With the availability of custom-tailored software tools, Pool-seq is being increasingly used for population genomic research on both model and non-model organisms. In this Review, we not only demonstrate the breadth of questions that are being addressed by Pool-seq but also discuss its limitations and provide guidelines for users.

Diversity in a hidden world: potential and limitation of next-generation sequencing for surveys of molecular diversity of eukaryotic microorganisms.

With the delivery of millions of sequence reads in a single experiment, next‐generation sequencing (NGS) is currently revolutionizing surveys of microorganism diversity. In particular, when applied to Eukaryotes, we are still lacking a rigorous comparison of morphological and NGS‐based diversity estimates. In this report, we studied the diversity and the seasonal community turnover of alveolates (Ciliophora and Dinophyceae) in an oligotrophic freshwater lake by SSU amplicon sequencing with NGS as well as by classical morphological analysis. We complemented the morphological analysis by single‐cell PCR followed by Sanger sequencing to provide an unambiguous link to the NGS data. We show that NGS and morphological analyses generally capture frequency shifts of abundant taxa over our seasonal samples. The observed incongruencies are probably largely due to rDNA copy number variation among taxa and heterogeneity in the efficiency of cell lysis. Overall, NGS‐based amplicon sequencing was superior in detecting rare species. We propose that in the absence of other nuclear markers less susceptible to copy number variation, rDNA‐based diversity studies need to be adjusted for confounding effects of copy number variation.

Genome-wide patterns of latitudinal differentiation among populations of Drosophila melanogaster from North America

Understanding the genetic underpinnings of adaptive change is a fundamental but largely unresolved problem in evolutionary biology. Drosophila melanogaster, an ancestrally tropical insect that has spread to temperate regions and become cosmopolitan, offers a powerful opportunity for identifying the molecular polymorphisms underlying clinal adaptation. Here, we use genome‐wide next‐generation sequencing of DNA pools (‘pool‐seq’) from three populations collected along the North American east coast to examine patterns of latitudinal differentiation. Comparing the genomes of these populations is particularly interesting since they exhibit clinal variation in a number of important life history traits. We find extensive latitudinal differentiation, with many of the most strongly differentiated genes involved in major functional pathways such as the insulin/TOR, ecdysone, torso, EGFR, TGFβ/BMP, JAK/STAT, immunity and circadian rhythm pathways. We observe particularly strong differentiation on chromosome 3R, especially within the cosmopolitan inversion In(3R)Payne, which contains a large number of clinally varying genes. While much of the differentiation might be driven by clinal differences in the frequency of In(3R)P, we also identify genes that are likely independent of this inversion. Our results provide genome‐wide evidence consistent with pervasive spatially variable selection acting on numerous loci and pathways along the well‐known North American cline, with many candidates implicated in life history regulation and exhibiting parallel differentiation along the previously investigated Australian cline.

Contrasting seasonal niche separation between rare and abundant taxa conceals the extent of protist diversity

With the advent of molecular methods, it became clear that microbial biodiversity had been vastly underestimated. Since then, species abundance patterns were determined for several environments, but temporal changes in species composition were not studied to the same level of resolution. Using massively parallel sequencing on the 454 GS FLX platform we identified a highly dynamic turnover of the seasonal abundance of protists in the Austrian lake Fuschlsee. We show that seasonal abundance patterns of protists closely match their biogeographic distribution. The stable predominance of few highly abundant taxa, which previously led to the suggestion of a low global protist species richness, is contrasted by a highly dynamic turnover of rare species. We suggest that differential seasonality of rare and abundant protist taxa explains the—so far—conflicting evidence in the ‘everything is everywhere’ dispute. Consequently temporal sampling is basic for adequate diversity and species richness estimates.

Massive habitat-specific genomic response in D. melanogaster populations during experimental evolution in hot and cold environments

Experimental evolution in combination with whole-genome sequencing (evolve and resequence [E&R]) is a promising approach to define the genotype–phenotype map and to understand adaptation in evolving populations. Many previous studies have identified a large number of putative selected sites (i.e., candidate loci), but it remains unclear to what extent these loci are genuine targets of selection or experimental noise. To address this question, we exposed the same founder population to two different selection regimes—a hot environment and a cold environment—and quantified the genomic response in each. We detected large numbers of putative selected loci in both environments, albeit with little overlap between the two sets of candidates, indicating that most resulted from habitat-specific selection. By quantifying changes across multiple independent biological replicates, we demonstrate that most of the candidate SNPs were false positives that were linked to selected sites over distances much larger than the typical linkage disequilibrium range of Drosophila melanogaster. We show that many of these mid- to long-range associations were attributable to large segregating inversions and confirm by computer simulations that such patterns could be readily replicated when strong selection acts on rare haplotypes. In light of our findings, we outline recommendations to improve the performance of future Drosophila E&R studies which include using species with negligible inversion loads, such as D. mauritiana and D. simulans, instead of D. melanogaster.

A Genome-Wide, Fine-Scale Map of Natural Pigmentation Variation in Drosophila melanogaster

Various approaches can be applied to uncover the genetic basis of natural phenotypic variation, each with their specific strengths and limitations. Here, we use a replicated genome-wide association approach (Pool-GWAS) to fine-scale map genomic regions contributing to natural variation in female abdominal pigmentation in Drosophila melanogaster, a trait that is highly variable in natural populations and highly heritable in the laboratory. We examined abdominal pigmentation phenotypes in approximately 8000 female European D. melanogaster, isolating 1000 individuals with extreme phenotypes. We then used whole-genome Illumina sequencing to identify single nucleotide polymorphisms (SNPs) segregating in our sample, and tested these for associations with pigmentation by contrasting allele frequencies between replicate pools of light and dark individuals. We identify two small regions near the pigmentation genes tan and bric-à-brac 1, both corresponding to known cis-regulatory regions, which contain SNPs showing significant associations with pigmentation variation. While the Pool-GWAS approach suffers some limitations, its cost advantage facilitates replication and it can be applied to any non-model system with an available reference genome.

Detecting Selective Sweeps from Pooled Next-Generation Sequencing Samples

Due to its cost effectiveness, next-generation sequencing of pools of individuals (Pool-Seq) is becoming a popular strategy for characterizing variation in population samples. Because Pool-Seq provides genome-wide SNP frequency data, it is possible to use them for demographic inference and/or the identification of selective sweeps. Here, we introduce a statistical method that is designed to detect selective sweeps from pooled data by accounting for statistical challenges associated with Pool-Seq, namely sequencing errors and random sampling among chromosomes. This allows for an efficient use of the information: all base calls are included in the analysis, but the higher credibility of regions with higher coverage and base calls with better quality scores is accounted for. Computer simulations show that our method efficiently detects sweeps even at very low coverage (0.5× per chromosome). Indeed, the power of detecting sweeps is similar to what we could expect from sequences of individual chromosomes. Since the inference of selective sweeps is based on the allele frequency spectrum (AFS), we also provide a method to accurately estimate the AFS provided that the quality scores for the sequence reads are reliable. Applying our approach to Pool-Seq data from Drosophila melanogaster, we identify several selective sweep signatures on chromosome X that include some previously well-characterized sweeps like the wapl region.

Host adaptation to viruses relies on few genes with different cross-resistance properties.

Significance Despite ample knowledge of the genetics and physiology of host responses to parasites, little is known about the genetic basis of host adaptation to parasites. Moreover, adaptation to one parasite is likely to impact the outcome of different infections. Yet these correlated responses, seminal to the understanding of host evolution in multiparasite environments, remain poorly studied. We determined the genetic and phenotypic changes underlying adaptation upon experimental evolution of a Drosophila melanogaster population under viral infection [Drosophila C virus (DCV)]. After 20 generations, selected flies showed increased survival upon infection with DCV and two other viruses. Using whole-genome sequencing and through RNAi, we identified and functionally validated three genes underlying the adaptive process and revealed their differential roles in the correlated responses observed. Host adaptation to one parasite may affect its response to others. However, the genetics of these direct and correlated responses remains poorly studied. The overlap between these responses is instrumental for the understanding of host evolution in multiparasite environments. We determined the genetic and phenotypic changes underlying adaptation of Drosophila melanogaster to Drosophila C virus (DCV). Within 20 generations, flies selected with DCV showed increased survival after DCV infection, but also after cricket paralysis virus (CrPV) and flock house virus (FHV) infection. Whole-genome sequencing identified two regions of significant differentiation among treatments, from which candidate genes were functionally tested with RNAi. Three genes were validated—pastrel, a known DCV-response gene, and two other loci, Ubc-E2H and CG8492. Knockdown of Ubc-E2H and pastrel also led to increased sensitivity to CrPV, whereas knockdown of CG8492 increased susceptibility to FHV infection. Therefore, Drosophila adaptation to DCV relies on few major genes, each with different cross-resistance properties, conferring host resistance to several parasites.

Genome-wide patterns of natural variation reveal strong selective sweeps and ongoing genomic conflict in Drosophila mauritiana.

Although it is well understood that selection shapes the polymorphism pattern in Drosophila, signatures of classic selective sweeps are scarce. Here, we focus on Drosophila mauritiana, an island endemic, which is closely related to Drosophila melanogaster. Based on a new, annotated genome sequence, we characterized the genome-wide polymorphism by sequencing pooled individuals (Pool-seq). We show that the interplay between selection and recombination results in a genome-wide polymorphism pattern characteristic for D. mauritiana. Two large genomic regions (>500 kb) showed the signature of almost complete selective sweeps. We propose that the absence of population structure and limited geographic distribution could explain why such pronounced sweep patterns are restricted to D. mauritiana. Further evidence for strong adaptive evolution was detected for several nucleoporin genes, some of which were not previously identified as genes involved in genomic conflict. Since this adaptive evolution is continuing after the split of D. mauritiana and Drosophila simulans, we conclude that genomic conflict is not restricted to short episodes, but rather an ongoing process in Drosophila.