Application of SNPs for population genetics of nonmodel organisms:new opportunities and challenges

Recent improvements in the speed, cost and accuracy of next generation sequencing are revolutionizing the discovery of single nucleotide polymorphisms (SNPs). SNPs are increasingly being used as an addition to the molecular ecology toolkit in nonmodel organisms, but their efficient use remains challenging. Here, we discuss common issues when employing SNP markers, including the high numbers of markers typically employed, the effects of ascertainment bias and the inclusion of nonneutral loci in a marker panel. We provide a critique of considerations specifically associated with the application and population genetic analysis of SNPs in nonmodel taxa, focusing specifically on some of the most commonly applied methods.

DiveRsity: An R package for the estimation and exploration of population genetics parameters and their associated errors

Summary: We present a new R package, diveRsity, for the calculation of various diversity statistics, including common diversity partitioning statistics (?, G) and population differentiation statistics (D, GST ', ? test for population heterogeneity), among others. The package calculates these estimators along with their respective bootstrapped confidence intervals for loci, sample population pairwise and global levels. Various plotting tools are also provided for a visual evaluation of estimated values, allowing users to critically assess the validity and significance of statistical tests from a biological perspective. diveRsity has a set of unique features, which facilitate the use of an informed framework for assessing the validity of the use of traditional F-statistics for the inference of demography, with reference to specific marker types, particularly focusing on highly polymorphic microsatellite loci. However, the package can be readily used for other co-dominant marker types (e.g. allozymes, SNPs). Detailed examples of usage and descriptions of package capabilities are provided. The examples demonstrate useful strategies for the exploration of data and interpretation of results generated by diveRsity. Additional online resources for the package are also described, including a GUI web app version intended for those with more limited experience using R for statistical analysis. © 2013 British Ecological Society.

Empirical assessment of non-invasive population genetics in bats:Comparison of DNA quality from faecal and tissue samples

Non-invasive population genetics has become a valuable tool in ecology and conservation biology, allowing genetic studies of wild populations without the need to catch, handle or even observe the study subjects directly. We address some of the concerns regarding the limitations of using non-invasive samples by comparing the quality of population genetic information gained through DNA extracted from faecal samples and biopsy samples of two elusive bat species, Myotis mystacinus and Myotis nattereri. We demonstrate that DNA extracted from faeces and tissue samples gives comparable results for frequency based population genetic analyses, despite the occurrence of genotyping errors when using faecal DNA. We conclude that non-invasive genetic sampling for population genetic analysis in bats is viable, and although more labour-intensive and expensive, it is an alternative to tissue sampling, which is particularly pertinent when specimens are rare, endangered or difficult to capture. © 2012 Museum and Institute of Zoology PAS.

A set of primers for plastid indels and nuclear microsatellites in the invasive plant Heracleum mantegazzianum (Apiaceae) and their transferability to Heracleum sphondylium

This study reports the isolation and polymorphism characterization of four plastid indels and six nuclear microsatellite loci in the invasive plant Heracleum mantegazzianum. These markers were tested in 27 individuals from two distant H. mantegazzianum populations. Plastid indels revealed the presence of five chlorotypes while five nuclear microsatellite loci rendered polymorphism. Applications of these markers include population genetics and phylogeography of H. mantegazzianum. A very good transferability of markers to Heracleum sphondylium was demonstrated.

Reducing the future threat from (liver) fluke: realistic prospect or quixotic fantasy?

The liver fluke remains an economically significant parasite of livestock and is emerging as an important zoonotic infection of humans. The incidence of the disease has increased in the last few years, as a possible consequence of changes to the World's climate. Future predictions suggest that this trend is likely to continue. Allied to the changing pattern of disease, reports of resistance to triclabendazole (TCBZ) have appeared in the literature, although they do not all represent genuine cases of resistance. Nevertheless, any reports of resistance are a concern, because triclabendazole is the only drug that has high activity against the migratory and damaging juvenile stages of infection. How to deal with the twin problems (of increasing incidence and drug resistance) is the overall theme of the session on “Trematodes: Fasciola hepatica epidemiology and control” and of this review to introduce the session.

Greater knowledge of fluke epidemiology and population genetics will highlight those regions where surveillance is most required and indicate how quickly resistant populations of fluke may arise. Models of disease risk are becoming increasingly sophisticated and precise, with more refined data analysis programmes and Geographic Information Systems (GIS) data. Recent improvements have been made in our understanding of the action of triclabendazole and the ways in which flukes have become resistant to it. While microtubules are the most likely target for drug action, tubulin mutations do not seem to be involved in the resistance mechanism. Rather, upregulation of drug uptake and metabolism processes appear to be more important and the data relating to them will be discussed. The information may help in the design of new treatment strategies or pinpoint potential molecular markers for monitoring fluke populations. Advances in the identification of novel targets for drugs and vaccines will be made by the various “-omics” technologies that are now being applied to Fasciola. A major area of concern in the current control of fasciolosis is the lack of reliable tests for the diagnosis of drug (TCBZ) resistance. This has led to inaccurate reports of resistance, which is hindering successful disease management, as farmers may be encouraged to switch to less effective drugs. Progress with the development of a number of new diagnostic tests will be reviewed.

Is the number of species on earth increasing or decreasing? Time, chaos and the origin of species

Darwin's On the Origin of Species has led to a theory of evolution with a mass of empirical detail on population genetics below species level, together with heated debate on the details of macroevolutionary patterns above species level. Most of the main principles are clear and generally accepted, notably that life originated once and has evolved over time by descent with modification. Here, I review the fossil and molecular phylogenetic records of the response of life on Earth to Quaternary climatic changes. I suggest that the record can be best understood in terms of the nonlinear dynamics of the relationship between genotype and phenotype, and between climate and environments. 'The origin of species' is essentially unpredictable, but is nevertheless an inevitable consequence of the way that organisms reproduce through time. The process is 'chaotic', but not 'random'. I suggest that biodiversity is best considered as continuously branching systems of lineages, where 'species' are the branch tips. The Earth's biodiversity should thus (1) be in a state of continuous increase and (2) show continuous discrepancies between genetic and morphological data in time and space. © The Palaeontological Association.

SNP discovery using Next Generation Transcriptomic Sequencing in Atlantic herring (Clupea harengus)

The introduction of Next Generation Sequencing (NGS) has revolutionised population genetics, providing studies of non-model species with unprecedented genomic coverage, allowing evolutionary biologists to address questions previously far beyond the reach of available resources. Furthermore, the simple mutation model of Single Nucleotide Polymorphisms (SNPs) permits cost-effective high-throughput genotyping in thousands of individuals simultaneously. Genomic resources are scarce for the Atlantic herring (Clupea harengus), a small pelagic species that sustains high revenue fisheries. This paper details the development of 578 SNPs using a combined NGS and high-throughput genotyping approach. Eight individuals covering the species distribution in the eastern Atlantic were bar-coded and multiplexed into a single cDNA library and sequenced using the 454 GS FLX platform. SNP discovery was performed by de novo sequence clustering and contig assembly, followed by the mapping of reads against consensus contig sequences. Selection of candidate SNPs for genotyping was conducted using an in silico approach. SNP validation and genotyping were performed simultaneously using an Illumina 1,536 GoldenGate assay. Although the conversion rate of candidate SNPs in the genotyping assay cannot be predicted in advance, this approach has the potential to maximise cost and time efficiencies by avoiding expensive and time-consuming laboratory stages of SNP validation. Additionally, the in silico approach leads to lower ascertainment bias in the resulting SNP panel as marker selection is based only on the ability to design primers and the predicted presence of intron-exon boundaries. Consequently SNPs with a wider spectrum of minor allele frequencies (MAFs) will be genotyped in the final panel. The genomic resources presented here represent a valuable multi-purpose resource for developing informative marker panels for population discrimination, microarray development and for population genomic studies in the wild.

Life history of turbot in Icelandic waters: Intra- and inter-population genetic diversity and otolith tracking of environmental temperatures

The stock structure of turbot was investigated between samples from S-Norway, the Irish Sea and the Kattegat, using 12 microsatellite loci and compared to the turbot caught in Icelandic waters. Highly significant genetic differentiation was observed between samples from Kattegat and other areas. Significant genetic differentiation was also observed between the Irish Sea sample on one hand and Iceland and S-Norway on the other hand. No significant genetic differentiation was observed between Iceland and S-Norway. Otoliths of 25 turbot, age ranging from 3 to 19 years, were subjected to nearly 300 mass spectrometry determinations of stable oxygen and carbon isotopes. Oxygen isotope composition (δ¹⁸O) in the otolith samples was used to estimate ambient temperature at time of otolith accretion, and yielded estimated temperatures experienced by the turbot ranging from 3 to 15°C. Overall, the genetic analysis indicates panmixia between turbot in Icelandic and Norwegian waters. While the extensive migration of larvae between Norway and Iceland is unlikely, passive drift of turbot larva from other areas (e.g. Ireland) cannot be ruled out. 

The occurrence and significance of triploidy in the liver fluke, Fasciola hepatica.

Karyotyping of Fasciola hepatica samples from Britain and Ireland has identified a triploid isolate which is effectively aspermic, rendering it necessarily asexually reproducing. Considering the extensive presence of asexually reproducing diploid and triploid Fasciola in Asia it is suggested that facultative gynogenesis is widespread in this parasite. This has important implications for the population genetics and evolution of Fasciola, especially in relation to the development and spread of drug resistance, and must be considered in the mathematical modelling of this process.

Development and characterization of 11 polymorphic compound tetranucleotide microsatellite loci for the Leisler's bat, Nyctalus leisleri (Vespertilionidae, Chiroptera)

Eleven polymorphic microsatellite marker loci were developed from a Leisler's bat (Nyctalus leisleri) genomic enriched library. Assessment of the usefulness of these markers for population genetics studies of Leisler's bats was carried out by screening 100 specimens sampled from eight locations in Ireland and two in Northeastern France. Both moderately and highly polymorphic marker loci were identified. Five to 28 alleles were found to be segregating per locus with observed heterozygosities values ranging from 28.4 to 94%. Initial evaluation indicates that these microsatellites will be useful for genetic based studies aiming, for instance, at parentage and population structure of Leisler's bats.

Unique genetic variation at a species' rear edge is under threat from global climate change

Global climate change is having a significant effect on the distributions of a wide variety of species, causing both range shifts and population extinctions. To date, however, no consensus has emerged on how these processes will affect the range-wide genetic diversity of impacted species. It has been suggested that species that recolonized from low-latitude refugia might harbour high levels of genetic variation in rear-edge populations, and that loss of these populations could cause a disproportionately large reduction in overall genetic diversity in such taxa. In the present study, we have examined the distribution of genetic diversity across the range of the seaweed Chondrus crispus, a species that has exhibited a northward shift in its southern limit in Europe over the last 40 years. Analysis of 19 populations from both sides of the North Atlantic using mitochondrial single nucleotide polymorphisms (SNPs), sequence data from two singlecopy nuclear regions and allelic variation at eight microsatellite loci revealed unique genetic variation for all marker classes in the rear-edge populations in Iberia, but not in the rear-edge populations in North America. Palaeodistribution modelling and statistical testing of alternative phylogeographic scenarios indicate that the unique genetic diversity in Iberian populations is a result not only of persistence in the region during the last glacial maximum, but also because this refugium did not contribute substantially to the recolonization of Europe after the retreat of the ice. Consequently, loss of these rear-edge populations as a result of ongoing climate change will have a major effect on the overall genetic diversity of the species, particularly in Europe, and this could compromise the adaptive potential of the species as a whole in the face of future global warming.

Further evidence for cryptic north-western refugia in Europe? Mitochondrial phylogeography of the sibling species Pipistrellus pipistrellus and Pipistrellus pygmaeus

The geographic ranges of European plants and animals underwent periods of contraction and re-colonisation during the climatic oscillations of the Pleistocene. The southern Mediterranean peninsulas (Iberian, Italian and Balkan) have been considered the most likely refugia for temperate/warm adapted species. Recent studies however have revealed the existence of extra-Mediterranean refugia, including the existence of cryptic north-west European refugia during the Last Glacial Maxima (24-14.6 kyr BP). In this study we elucidated the phylogeographic history of two sibling bat species, Pipistrellus pipistrellus and P. pygmaeus in their western European range. We sequenced the highly variable mtDNA D-loop for 167 samples of P. pipistrellus (n = 99) and P. pygmaeus (n = 68) and combined our data with published sequences from 331 individuals. Using phylogenetic methodologies we assessed their biogeographic history. Our data support a single eastern European origin for populations of P. pygmaeus s.str., yet multiple splits and origins for populations of P. pipistrellus s.str., including evidence for refugia within refugia and potential cryptic refugia in north western Europe and in the Caucasus. This complex pattern in the distribution of mtDNA haplotypes supports a long history for P. pipistrellus s.str. in Europe, and the hypothesis that species with a broad ecological niche may have adapted and survived outside southern peninsula throughout the LGM.

The paradox of invasion: Reeves’ muntjac deer invade the British Isles from a limited number of founding females

High levels of genetic diversity and high propagule pressure are favoured by conservation biologists as the basis for successful reintroductions and ensuring the persistence of populations. However, invasion ecologists recognize the ‘paradox of invasion’, as successful species introductions may often be characterized by limited numbers of individuals and associated genetic bottlenecks. In the present study, we used a combination of high-resolution nuclear and mitochondrial genetic markers to investigate the invasion history of Reeves' muntjac deer in the British Isles. This invasion has caused severe economic and ecological damage, with secondary spread currently a concern throughout Europe and potentially globally. Microsatellite analysis based on eight loci grouped all 176 introduced individuals studied from across the species' range in the UK into one genetic cluster, and seven mitochondrial D-loop haplotypes were recovered, two of which were present at very low frequency and were related to more common haplotypes. Our results indicate that the entire invasion can be traced to a single founding event involving a low number of females. These findings highlight the fact that even small releases of species may, if ignored, result in irreversible and costly invasion, regardless of initial genetic diversity or continual genetic influx.