Integrating molecular tools into freshwater ecology: developments and opportunities

相似文献   

Exploring a Pool‐seq‐only approach for gaining population genomic insights in nonmodel species

Developing genomic insights is challenging in nonmodel species for which resources are often scarce and prohibitively costly. Here, we explore the potential of a recently established approach using Pool‐seq data to generate a de novo genome assembly for mining exons, upon which Pool‐seq data are used to estimate population divergence and diversity. We do this for two pairs of sympatric populations of brown trout (Salmo trutta): one naturally sympatric set of populations and another pair of populations introduced to a common environment. We validate our approach by comparing the results to those from markers previously used to describe the populations (allozymes and individual‐based single nucleotide polymorphisms [SNPs]) and from mapping the Pool‐seq data to a reference genome of the closely related Atlantic salmon (Salmo salar). We find that genomic differentiation (F_ST) between the two introduced populations exceeds that of the naturally sympatric populations (F_ST = 0.13 and 0.03 between the introduced and the naturally sympatric populations, respectively), in concordance with estimates from the previously used SNPs. The same level of population divergence is found for the two genome assemblies, but estimates of average nucleotide diversity differ ( ≈ 0.002 and  ≈ 0.001 when mapping to S. trutta and S. salar, respectively), although the relationships between population values are largely consistent. This discrepancy might be attributed to biases when mapping to a haploid condensed assembly made of highly fragmented read data compared to using a high‐quality reference assembly from a divergent species. We conclude that the Pool‐seq‐only approach can be suitable for detecting and quantifying genome‐wide population differentiation, and for comparing genomic diversity in populations of nonmodel species where reference genomes are lacking.     

Genotyping‐by‐sequencing for estimating relatedness in nonmodel organisms: Avoiding the trap of precise bias

There has been remarkably little attention to using the high resolution provided by genotyping‐by‐sequencing (i.e., RADseq and similar methods) for assessing relatedness in wildlife populations. A major hurdle is the genotyping error, especially allelic dropout, often found in this type of data that could lead to downward‐biased, yet precise, estimates of relatedness. Here, we assess the applicability of genotyping‐by‐sequencing for relatedness inferences given its relatively high genotyping error rate. Individuals of known relatedness were simulated under genotyping error, allelic dropout and missing data scenarios based on an empirical ddRAD data set, and their true relatedness was compared to that estimated by seven relatedness estimators. We found that an estimator chosen through such analyses can circumvent the influence of genotyping error, with the estimator of Ritland (Genetics Research, 67, 175) shown to be unaffected by allelic dropout and to be the most accurate when there is genotyping error. We also found that the choice of estimator should not rely solely on the strength of correlation between estimated and true relatedness as a strong correlation does not necessarily mean estimates are close to true relatedness. We also demonstrated how even a large SNP data set with genotyping error (allelic dropout or otherwise) or missing data still performs better than a perfectly genotyped microsatellite data set of tens of markers. The simulation‐based approach used here can be easily implemented by others on their own genotyping‐by‐sequencing data sets to confirm the most appropriate and powerful estimator for their data.     

Estimation of population allele frequencies from next‐generation sequencing data: pool‐versus individual‐based genotyping

Molecular markers produced by next‐generation sequencing (NGS) technologies are revolutionizing genetic research. However, the costs of analysing large numbers of individual genomes remain prohibitive for most population genetics studies. Here, we present results based on mathematical derivations showing that, under many realistic experimental designs, NGS of DNA pools from diploid individuals allows to estimate the allele frequencies at single nucleotide polymorphisms (SNPs) with at least the same accuracy as individual‐based analyses, for considerably lower library construction and sequencing efforts. These findings remain true when taking into account the possibility of substantially unequal contributions of each individual to the final pool of sequence reads. We propose the intuitive notion of effective pool size to account for unequal pooling and derive a Bayesian hierarchical model to estimate this parameter directly from the data. We provide a user‐friendly application assessing the accuracy of allele frequency estimation from both pool‐ and individual‐based NGS population data under various sampling, sequencing depth and experimental error designs. We illustrate our findings with theoretical examples and real data sets corresponding to SNP loci obtained using restriction site–associated DNA (RAD) sequencing in pool‐ and individual‐based experiments carried out on the same population of the pine processionary moth (Thaumetopoea pityocampa). NGS of DNA pools might not be optimal for all types of studies but provides a cost‐effective approach for estimating allele frequencies for very large numbers of SNPs. It thus allows comparison of genome‐wide patterns of genetic variation for large numbers of individuals in multiple populations.     

The design and cross-population application of a genome-wide SNP chip for the great tit Parus major

The vast amount of phenotypic information collected in some wild animal populations makes them extremely valuable for unravelling the genetics of ecologically important traits and understanding how populations adapt to changes in their environment. Next generation sequencing has revolutionized the development of large marker panels in species previously lacking genomic resources. In this study, a unique genomics toolkit was developed for the great tit (Parus major), a model species in ecology and behavioural biology. This toolkit consists of nearly 100,000 SNPs, over 250 million nucleotides of assembled genomic DNA and more than 80 million nucleotides of assembled expressed sequences. A SNP chip with 9193 SNP markers expected to be spaced evenly along the great tit genome was used to genotype 4702 birds from two of the most intensively studied natural vertebrate populations [Wytham Woods/Bagley Woods (United Kingdom) and de Hoge Veluwe/Westerheide (The Netherlands)]. We show that (i) SNPs identified in either of the two populations have a high genotyping success in the other population, (ii) the minor allele frequencies of the SNPs are highly correlated between the two populations and (iii) despite this high correlation, a large number of SNPs display significant differentiation (F(ST) ) between the populations, with an overrepresentation of genes involved in cardiovascular development close to these SNPs. The developed resources provide the basis for unravelling the genetics of important traits in many long-term studies of great tits. More generally, the protocols and pitfalls encountered will be of use for those developing similar resources.     

Recent progress and challenges in population genetics of polyploid organisms: an overview of current state‐of‐the‐art molecular and statistical tools

Despite the importance of polyploidy and the increasing availability of new genomic data, there remain important gaps in our knowledge of polyploid population genetics. These gaps arise from the complex nature of polyploid data (e.g. multiple alleles and loci, mixed inheritance patterns, association between ploidy and mating system variation). Furthermore, many of the standard tools for population genetics that have been developed for diploids are often not feasible for polyploids. This review aims to provide an overview of the state‐of‐the‐art in polyploid population genetics and to identify the main areas where further development of molecular techniques and statistical theory is required. We review commonly used molecular tools (amplified fragment length polymorphism, microsatellites, Sanger sequencing, next‐generation sequencing and derived technologies) and their challenges associated with their use in polyploid populations: that is, allele dosage determination, null alleles, difficulty of distinguishing orthologues from paralogues and copy number variation. In addition, we review the approaches that have been used for population genetic analysis in polyploids and their specific problems. These problems are in most cases directly associated with dosage uncertainty and the problem of inferring allele frequencies and assumptions regarding inheritance. This leads us to conclude that for advancing the field of polyploid population genetics, most priority should be given to development of new molecular approaches that allow efficient dosage determination, and to further development of analytical approaches to circumvent dosage uncertainty and to accommodate ‘flexible’ modes of inheritance. In addition, there is a need for more simulation‐based studies that test what kinds of biases could result from both existing and novel approaches.     

The easy road to genome‐wide medium density SNP screening in a non‐model species: development and application of a 10 K SNP‐chip for the house sparrow (Passer domesticus)

With the advent of next generation sequencing, new avenues have opened to study genomics in wild populations of non‐model species. Here, we describe a successful approach to a genome‐wide medium density Single Nucleotide Polymorphism (SNP) panel in a non‐model species, the house sparrow (Passer domesticus), through the development of a 10 K Illumina iSelect HD BeadChip. Genomic DNA and cDNA derived from six individuals were sequenced on a 454 GS FLX system and generated a total of 1.2 million sequences, in which SNPs were detected. As no reference genome exists for the house sparrow, we used the zebra finch (Taeniopygia guttata) reference genome to determine the most likely position of each SNP. The 10 000 SNPs on the SNP‐chip were selected to be distributed evenly across 31 chromosomes, giving on average one SNP per 100 000 bp. The SNP‐chip was screened across 1968 individual house sparrows from four island populations. Of the original 10 000 SNPs, 7413 were found to be variable, and 99% of these SNPs were successfully called in at least 93% of all individuals. We used the SNP‐chip to demonstrate the ability of such genome‐wide marker data to detect population sub‐division, and compared these results to similar analyses using microsatellites. The SNP‐chip will be used to map Quantitative Trait Loci (QTL) for fitness‐related phenotypic traits in natural populations.     

Discovery and development of exome‐based,co‐dominant single nucleotide polymorphism markers in hexaploid wheat (Triticum aestivum L.)

Globally, wheat is the most widely grown crop and one of the three most important crops for human and livestock feed. However, the complex nature of the wheat genome has, until recently, resulted in a lack of single nucleotide polymorphism (SNP)‐based molecular markers of practical use to wheat breeders. Recently, large numbers of SNP‐based wheat markers have been made available via the use of next‐generation sequencing combined with a variety of genotyping platforms. However, many of these markers and platforms have difficulty distinguishing between heterozygote and homozygote individuals and are therefore of limited use to wheat breeders carrying out commercial‐scale breeding programmes. To identify exome‐based co‐dominant SNP‐based assays, which are capable of distinguishing between heterozygotes and homozygotes, we have used targeted re‐sequencing of the wheat exome to generate large amounts of genomic sequences from eight varieties. Using a bioinformatics approach, these sequences have been used to identify 95 266 putative single nucleotide polymorphisms, of which 10 251 were classified as being putatively co‐dominant. Validation of a subset of these putative co‐dominant markers confirmed that 96% were true polymorphisms and 65% were co‐dominant SNP assays. The new co‐dominant markers described here are capable of genotypic classification of a segregating locus in polyploid wheat and can be used on a variety of genotyping platforms; as such, they represent a powerful tool for wheat breeders. These markers and related information have been made publically available on an interactive web‐based database to facilitate their use on genotyping programmes worldwide.     

Phylogeographic and phenotypic outcomes of brown anole colonization across the Caribbean provide insight into the beginning stages of an adaptive radiation

Some of the most important insights into the ecological and evolutionary processes of diversification and speciation have come from studies of island adaptive radiations, yet relatively little research has examined how these radiations initiate. We suggest that Anolis sagrei is a candidate for understanding the origins of the Caribbean Anolis adaptive radiation and how a colonizing anole species begins to undergo allopatric diversification, phenotypic divergence and, potentially, speciation. We undertook a genomic and morphological analysis of representative populations across the entire native range of A. sagrei, finding that the species originated in the early Pliocene, with the deepest divergence occurring between western and eastern Cuba. Lineages from these two regions subsequently colonized the northern Caribbean. We find that at the broadest scale, populations colonizing areas with fewer closely related competitors tend to evolve larger body size and more lamellae on their toepads. This trend follows expectations for post‐colonization divergence from progenitors and convergence in allopatry, whereby populations freed from competition with close relatives evolve towards common morphological and ecological optima. Taken together, our results show a complex history of ancient and recent Cuban diaspora with populations on competitor‐poor islands evolving away from their ancestral Cuban populations regardless of their phylogenetic relationships, thus providing insight into the original diversification of colonist anoles at the beginning of the radiation. Our research also supplies an evolutionary framework for the many studies of this increasingly important species in ecological and evolutionary research.