SNPs in ecological and conservation studies: a test in the Scandinavian wolf population

The IGF1 small dog haplotype is derived from Middle Eastern grey wolves: a closer look at statistics, sampling, and the alleged Middle Eastern origin of small dogs

This paper is a response to Gray MM, Sutter NB, Ostrander EA, Wayne RK: The IGF1 small dog haplotype is derived from Middle Eastern grey wolves. BMC Biology 2010, 8:16.

See research article at http://www.biomedcentral.com/1741-7007/8/16.

A comparison of SNP and STR loci for delineating population structure and performing individual genetic assignment

Background

Technological advances have lead to the rapid increase in availability of single nucleotide polymorphisms (SNPs) in a range of organisms, and there is a general optimism that SNPs will become the marker of choice for a range of evolutionary applications. Here, comparisons between 300 polymorphic SNPs and 14 short tandem repeats (STRs) were conducted on a data set consisting of approximately 500 Atlantic salmon arranged in 10 samples/populations.

Results

Global F_ST ranged from 0.033-0.115 and -0.002-0.316 for the 14 STR and 300 SNP loci respectively. Global F_ST was similar among 28 linkage groups when averaging data from mapped SNPs. With the exception of selecting a panel of SNPs taking the locus displaying the highest global F_ST for each of the 28 linkage groups, which inflated estimation of genetic differentiation among the samples, inferred genetic relationships were highly similar between SNP and STR data sets and variants thereof. The best 15 SNPs (30 alleles) gave a similar level of self-assignment to the best 4 STR loci (83 alleles), however, addition of further STR loci did not lead to a notable increase assignment whereas addition of up to 100 SNP loci increased assignment.

Conclusion

Whilst the optimal combinations of SNPs identified in this study are linked to the samples from which they were selected, this study demonstrates that identification of highly informative SNP loci from larger panels will provide researchers with a powerful approach to delineate genetic relationships at the individual and population levels.

Geographical patterning of sixteen goat breeds from Italy, Albania and Greece assessed by Single Nucleotide Polymorphisms

Background

SNP data of goats of three Mediterranean countries were used for population studies and reconstruction of geographical patterning. 496 individuals belonging to Italian, Albanian and Greek breeds were genotyped to assess the basic population parameters.

Results

A total of 26 SNPs were used, for a total of 12,896 genotypes assayed. Statistical analysis revealed that breeds are not so similar in terms of genetic variability, as reported in studies performed using different markers. The Mantel test showed a strongly significant correlation between genetic and geographic distance. Also, PCA analysis revealed that breeds are grouped according to geographical origin, with the exception of the Greek Skopelos breed.

Conclusion

Our data point out that the use of SNP markers to analyze a wider breed sample could help in understanding the recent evolutionary history of domestic goats. We found correlation between genetic diversity and geographic distance. Also PCA analysis shows that the breeds are well differentiated, with good correspondence to geographical locations, thus confirming the correlation between geographical and genetic distances. This suggests that migration history of the species played a pivotal role in the present-day structure of the breeds and a scenario in which coastal routes were easier for migrating in comparison with inland routes. A westward coastal route to Italy through Greece could have led to gene flow along the Northern Mediterranean.

Comparative performance of single nucleotide polymorphism and microsatellite markers for population genetic analysis

Microsatellite loci are standard genetic markers for population genetic analysis, whereas single nucleotide polymorphisms (SNPs) are more recent tools that require assessment of neutrality and appropriate use in population genetics. Twelve SNP markers were used to describe the genetic structure of Diabrotica virgifera virgifera (LeConte; Coleoptera: Chrysomelidae) in the United States of America and revealed a high mean observed heterozygosity (0.40 +/- 0.059) and low global F(ST) (0.029). Pairwise F(ST) estimates ranged from 0.007 to 0.045, and all but 2 populations showed significant levels of genetic differentiation (P < or = 0.008). Population parameters and conclusions based on SNP markers were analogous to that obtained by use of microsatellite markers from the identical population samples. SNP-based F(ST) estimates were 3-fold higher than corresponding estimates from microsatellites, wherein lower microsatellite F(ST) estimates likely resulted from an overestimate of migration rates between subpopulations due to convergence of allele size (homoplasy). No significant difference was observed in the proportion of SNP or microsatellite markers loci that were nonneutral within populations. SNP markers provided estimates of population genetic parameters consistent with those from microsatellite data, and their low back mutation rates may result in reduced propensity for error in estimation of population parameters.

Linkage disequilibrium and demographic history of wild and domestic canids

Assessing the extent of linkage disequilibrium (LD) in natural populations of a nonmodel species has been difficult due to the lack of available genomic markers. However, with advances in genotyping and genome sequencing, genomic characterization of natural populations has become feasible. Using sequence data and SNP genotypes, we measured LD and modeled the demographic history of wild canid populations and domestic dog breeds. In 11 gray wolf populations and one coyote population, we find that the extent of LD as measured by the distance at which r2=0.2 extends <10 kb in outbred populations to >1.7 Mb in populations that have experienced significant founder events and bottlenecks. This large range in the extent of LD parallels that observed in 18 dog breeds where the r2 value varies from approximately 20 kb to >5 Mb. Furthermore, in modeling demographic history under a composite-likelihood framework, we find that two of five wild canid populations exhibit evidence of a historical population contraction. Five domestic dog breeds display evidence for a minor population contraction during domestication and a more severe contraction during breed formation. Only a 5% reduction in nucleotide diversity was observed as a result of domestication, whereas the loss of nucleotide diversity with breed formation averaged 35%.

Population genomics of the inbred Scandinavian wolf

The Scandinavian wolf population represents one of the genetically most well-characterized examples of a severely bottlenecked natural population (with only two founders), and of how the addition of new genetic material (one immigrant) can at least temporarily provide a 'genetic rescue'. However, inbreeding depression has been observed in this population and in the absence of additional immigrants, its long-term viability is questioned. To study the effects of inbreeding and selection on genomic diversity, we performed a genomic scan with approximately 250 microsatellite markers distributed across all autosomes and the X chromosome. We found linkage disequilibrium (LD) that extended up to distances of 50 Mb, exceeding that of most outbreeding species studied thus far. LD was particularly pronounced on the X chromosome. Overall levels of observed genomic heterozygosity did not deviate significantly from simulations based on known population history, giving no support for a general selection for heterozygotes. However, we found evidence supporting balancing selection at a number of loci and also evidence suggesting directional selection at other loci. For markers on chromosome 23, the signal of selection was particularly strong, indicating that purifying selection against deleterious alleles may have occurred even in this very small population. These data suggest that population genomics allows the exploration of the effects of neutral and non-neutral evolution on a finer scale than what has previously been possible.

SNP genotyping by the 5'-nuclease reaction: advances in high-throughput genotyping with nonmodel organisms

Population genetics studies play an increasingly important role in the management and conservation of nonmodel organisms. Unlike studies with model organisms, a typical population genetics study of a nonmodel organism may be conducted by analyzing thousands or hundreds of thousands of individuals for several dozen single nucleotide polymorphisms (SNPs). The use of robust, robotically mediated TaqMan reactions provides substantial advantages in these types of studies. We describe the methods and laboratory setup for analyzing a sustained high throughput of SNP assays in routine university or natural resource agency laboratories with a handful of thermal cyclers. Agencies sustain rates of nearly 150,000 assays per week using uniplex reactions with the Applied Biosystems 7900HT Fast Real-Time PCR System (AB 7900HT). We further describe the medium-density array run on the BioMark from Fluidigm, which increases this rate to over 500,000 assays per week by multiplexing 96 samples for 96 SNPs.

Genetics in conservation and wildlife management: a revolution since Caughley

In his 1994 review of conservation biology, Graeme Caughley questioned the central role for genetics in that discipline. His central theme was that there was no known case of genetic malfunction leading to the extinction of a population or species, and that driving forces such as overkill, habitat fragmentation and introduced predators as well as environmental and demographic stochasticity of small populations should be considered ahead of genetics in the debate about extinction prevention. At the time, only indirect and theoretical evidence existed for genetic contributions to the declines of wildlife and most of the debate revolved around the impact of genetic variation on fitness and long-term persistence. In addition, the application of DNA technologies to the study of wildlife was in its infancy. Though this was not Caughley’s intention, many within wildlife management took his criticisms of genetic aspects of species decline as the cue to dismiss this branch of science as of minor relevance to conservation biology. Since Caughley’s critique, there has been a revolution in genetic technologies for non-model organisms with the arrival of highly informative hypervariable DNA markers. Perhaps even more importantly, developments in DNA and gene technologies have provided the opportunity to study fundamental life-history traits such as disease resistance in more direct ways than previously possible. In concert with these tools, conservation geneticists have risen to Caughley’s challenge and demonstrated unambiguously a clear role for genetic analysis in conservation biology. Despite these impressive advances, there remains an important gap between the genetic approaches available and their uptake by managers. Bridging this gap will greatly increase the capacity of wildlife managers to generate the data necessary for sound management.

Isolation and characterization of SNP variation at 90 anonymous loci in the banded wren (Thryothorus pleurostictus)

Single nucleotide polymorphisms (SNPs) are becoming more commonly used as molecular markers in conservation studies. However, relatively few studies have employed SNPs for species with little or no existing sequence data, partly due to the practical challenge of locating appropriate SNP loci in these species. Here we describe an application of SNP discovery via shotgun cloning that requires no pre-existing sequence data and is readily applied to all taxa. Using this method, we isolated, cloned and screened for SNP variation at 90 anonymous sequence loci (51kb total) from the banded wren (Thryothorus pleurostictus), a Central American species with minimal pre-existing sequence data. We identified 168 SNPs (a mean of one SNP/305 bp, with SNPs unevenly distributed across loci). Further characterization of variation at 41 of these SNP loci among 256 individuals including 37 parent-offspring families suggests that they provide substantial information for defining the genetic mating system of this species, and that SNPs may be generally useful for this purpose when other markers are problematic.

Molecular ecological approaches to studying the evolutionary impact of selective harvesting in wildlife

Harvesting of wildlife populations by humans is usually targeted by sex, age or phenotypic criteria, and is therefore selective. Selective harvesting has the potential to elicit a genetic response from the target populations in several ways. First, selective harvesting may affect population demographic structure (age structure, sex ratio), which in turn may have consequences for effective population size and hence genetic diversity. Second, wildlife-harvesting regimes that use selective criteria based on phenotypic characteristics (e.g. minimum body size, horn length or antler size) have the potential to impose artificial selection on harvested populations. If there is heritable genetic variation for the target characteristic and harvesting occurs before the age of maturity, then an evolutionary response over time may ensue. Molecular ecological techniques offer ways to predict and detect genetic change in harvested populations, and therefore have great utility for effective wildlife management. Molecular markers can be used to assess the genetic structure of wildlife populations, and thereby assist in the prediction of genetic impacts by delineating evolutionarily meaningful management units. Genetic markers can be used for monitoring genetic diversity and changes in effective population size and breeding systems. Tracking evolutionary change at the phenotypic level in the wild through quantitative genetic analysis can be made possible by genetically determined pedigrees. Finally, advances in genome sequencing and bioinformatics offer the opportunity to study the molecular basis of phenotypic variation through trait mapping and candidate gene approaches. With this understanding, it could be possible to monitor the selective impacts of harvesting at a molecular level in the future. Effective wildlife management practice needs to consider more than the direct impact of harvesting on population dynamics. Programs that utilize molecular genetic tools will be better positioned to assess the long-term evolutionary impact of artificial selection on the evolutionary trajectory and viability of harvested populations.

Insertion-deletion polymorphisms (indels) as genetic markers in natural populations

BACKGROUND

We introduce the use of short insertion-deletion polymorphisms (indels) for genetic analysis of natural populations.

RESULTS

Sequence reads from light shot-gun sequencing efforts of different dog breeds were aligned to the dog genome reference sequence and gaps corresponding to indels were identified. One hundred candidate markers (4-bp indels) were selected and genotyped in unrelated dogs (n = 7) and wolves (n = 18). Eighty-one and 76 out of 94 could be validated as polymorphic loci in the respective sample. Mean indel heterozygosity in a diverse set of wolves was 19%, and 74% of the loci had a minor allele frequency of >10%. Indels found to be polymorphic in wolves were subsequently genotyped in a highly bottlenecked Scandinavian wolf population. Fifty-one loci turned out to be polymorphic, showing their utility even in a population with low genetic diversity. In this population, individual heterozygosity measured at indel and microsatellite loci were highly correlated.

CONCLUSION

With an increasing amount of sequence information gathered from non-model organisms, we suggest that indels will come to form an important source of genetic markers, easy and cheap to genotype, for studies of natural populations.

Estimating the number of founder lineages from haplotypes of closely linked SNPs

We consider an isolated population founded by a small number of individuals randomly chosen from a source population of known genetic composition at a known time in the past. We develop a Monte-Carlo maximum-likelihood method for estimating the number of founding individuals from the haplotype frequencies at several SNP (single nucleotide polymorphism) loci in a sample. We assume the isolated population was founded recently enough that that mutation can be ignored and that haplotype frequencies in the source population have not changed. We apply the method to simulated data and show that it is unbiased. With a reasonable number of individuals sampled, it is possible to estimate the number of founders within a factor of 2. We show that the performance of the method is not degraded substantially if the frequencies of the rare haplotypes in the source are not known precisely and if there is some recombination. We illustrate the use of our method by applying it to a previously published data set from a recently founded population of wolves (Canis lupus) in Scandinavia.

SNP discovery by mismatch-targeting of Mu transposition

Single nucleotide polymorphisms (SNPs) represent a valuable resource for the mapping of human disease genes and induced mutations in model organisms. SNPs may become the markers of choice also for population ecology and evolutionary studies, but their isolation for non-model organisms with unsequenced genomes is often difficult. Here, we describe a rapid and cost-effective strategy to isolate SNPs that exploits the property of the bacteriophage Mu transposition machinery to target mismatched DNA sites and thereby to effectively detect polymorphic loci. To demonstrate the methodology, we isolated 164 SNPs from the unsequenced genome of the Glanville fritillary butterfly (Melitaea cinxia), a much-studied species in population biology, and we validated 24 of them. The strategy involves standard molecular biology techniques as well as undemanding MuA transposase-catalyzed in vitro transposition reactions, and it is applicable to any organism.

Genomics and conservation genetics

In large part, the relevance of genetics to conservation rests on the premise that neutral marker variation in populations reflects levels of detrimental and adaptive genetic variation. Despite its prominence, this tenet has been difficult to evaluate, until now. As we discuss here, genome sequence information and new technological and bioinformatics platforms now enable comprehensive surveys of neutral variation and more direct inferences of detrimental and adaptive variation in species with sequenced genomes and in 'genome-enabled' endangered taxa. Moreover, conservation schemes could begin to consider specific pathological genetic variants. A new conservation genetic agenda would utilize data from enhanced surveys of genomic variation in endangered species to better manage functional genetic variation.

Characterization of 37 Breed-Specific Single-Nucleotide Polymorphisms in Sheep

We identified 37 single-nucleotide polymorphisms (SNPs) in sheep and screened 16 individuals from 8 different sheep breeds selected throughout Europe. Population genetic measures based on the genotyping of about 30 sheep from the same 8 breeds are reported. To date, there are no sheep SNPs documented in the National Center for Biotechnology Information dbSNP database. Therefore, the markers presented here contribute significantly to those currently available.

Single nucleotide polymorphism (SNP) discovery in duplicated genomes: intron-primed exon-crossing (IPEC) as a strategy for avoiding amplification of duplicated loci in Atlantic salmon (Salmo salar) and other salmonid fishes

Background

Single nucleotide polymorphisms (SNPs) represent the most abundant type of DNA variation in the vertebrate genome, and their applications as genetic markers in numerous studies of molecular ecology and conservation of natural populations are emerging. Recent large-scale sequencing projects in several fish species have provided a vast amount of data in public databases, which can be utilized in novel SNP discovery in salmonids. However, the suggested duplicated nature of the salmonid genome may hamper SNP characterization if the primers designed in conserved gene regions amplify multiple loci.

Results

Here we introduce a new intron-primed exon-crossing (IPEC) method in an attempt to overcome this duplication problem, and also evaluate different priming methods for SNP discovery in Atlantic salmon (Salmo salar) and other salmonids. A total of 69 loci with differing priming strategies were screened in S. salar, and 27 of these produced ~13 kb of high-quality sequence data consisting of 19 SNPs or indels (one per 680 bp). The SNP frequency and the overall nucleotide diversity (3.99 × 10^-4) in S. salar was lower than reported in a majority of other organisms, which may suggest a relative young population history for Atlantic salmon. A subset of primers used in cross-species analyses revealed considerable variation in the SNP frequencies and nucleotide diversities in other salmonids.

Conclusion

Sequencing success was significantly higher with the new IPEC primers; thus the total number of loci to screen in order to identify one potential polymorphic site was six times less with this new strategy. Given that duplication may hamper SNP discovery in some species, the IPEC method reported here is an alternative way of identifying novel polymorphisms in such cases.

Genetic diversity, population structure, effective population size and demographic history of the Finnish wolf population

The Finnish wolf population (Canis lupus) was sampled during three different periods (1996-1998, 1999-2001 and 2002-2004), and 118 individuals were genotyped with 10 microsatellite markers. Large genetic variation was found in the population despite a recent demographic bottleneck. No spatial population subdivision was found even though a significant negative relationship between genetic relatedness and geographic distance suggested isolation by distance. Very few individuals did not belong to the local wolf population as determined by assignment analyses, suggesting a low level of immigration in the population. We used the temporal approach and several statistical methods to estimate the variance effective size of the population. All methods gave similar estimates of effective population size, approximately 40 wolves. These estimates were slightly larger than the estimated census size of breeding individuals. A Bayesian model based on Markov chain Monte Carlo simulations indicated strong evidence for a long-term population decline. These results suggest that the contemporary wolf population size is roughly 8% of its historical size, and that the population decline dates back to late 19th century or early 20th century. Despite an increase of over 50% in the census size of the population during the whole study period, there was only weak evidence that the effective population size during the last period was higher than during the first. This may be caused by increased inbreeding, diminished dispersal within the population, and decreased immigration to the population during the last study period.

Use of sequence data from rainbow trout and Atlantic salmon for SNP detection in Pacific salmon

Single nucleotide polymorphisms (SNPs) are a class of genetic markers that are well suited to a broad range of research and management applications. Although advances in genotyping chemistries and analysis methods continue to increase the potential advantages of using SNPs to address molecular ecological questions, the scarcity of available DNA sequence data for most species has limited marker development. As the number and diversity of species being targeted for large-scale sequencing has increased, so has the potential for using sequence from sister taxa for marker development in species of interest. We evaluated the use of Oncorhynchus mykiss and Salmo salar sequence data to identify SNPs in three other species (Oncorhynchus tshawytscha, Oncorhynchus nerka and Oncorhynchus keta). Primers designed based on O. mykiss and S. salar alignments were more successful than primers designed based on Oncorhynchus-only alignments for sequencing target species, presumably due to the much larger number of potential targets available from the former alignments and possibly greater sequence conservation in those targets. In sequencing approximately 89 kb we observed a frequency of 4.30 x 10(-3) SNPs per base pair. Approximately half (53/101) of the subsequently designed validation assays resulted in high-throughput SNP genotyping markers. We speculate that this relatively low conversion rate may reflect the duplicated nature of the salmon genome. Our results suggest that a large number of SNPs could be developed for Pacific salmon using sequence data from other species. While the costs of DNA sequencing are still significant, these must be compared to the costs of using other marker classes for a given application.