Assessing statistical power of SNPs for population structure and conservation studies

Single nucleotide polymorphisms (SNPs) have been proposed by some as the new frontier for population studies, and several papers have presented theoretical and empirical evidence reporting the advantages and limitations of SNPs. As a practical matter, however, it remains unclear how many SNP markers will be required or what the optimal characteristics of those markers should be in order to obtain sufficient statistical power to detect different levels of population differentiation. We use a hypothetical case to illustrate the process of designing a population genetics project, and present results from simulations that address several issues for maximizing statistical power to detect differentiation while minimizing the amount of effort in developing SNPs. Results indicate that (i) while ~30 SNPs should be sufficient to detect moderate (F_ST = 0.01) levels of differentiation, studies aimed at detecting demographic independence (e.g. F_ST < 0.005) may require 80 or more SNPs and large sample sizes; (ii) different SNP allele frequencies have little affect on power, and thus, selection of SNPs can be relatively unbiased; (iii) increasing the sample size has a strong effect on power, so that the number of loci can be minimized when sample number is known, and increasing sample size is almost always beneficial; and (iv) power is increased by including multiple SNPs within loci and inferring haplotypes, rather than trying to use only unlinked SNPs. This also has the practical benefit of reducing the SNP ascertainment effort, and may influence the decision of whether to seek SNPs in coding or noncoding regions.     

Quantifying ascertainment bias and species-specific length differences in human and chimpanzee microsatellites using genome sequences

Surveys of variability of homologous microsatellite loci among species reveal an ascertainment bias for microsatellite length where microsatellite loci isolated in one species tend to be longer than homologous loci in related species. Here, we take advantage of the availability of aligned human and chimpanzee genome sequences to compare length difference of homologous microsatellites for loci identified in humans to length difference for loci identified in chimpanzees. We are able to quantify ascertainment bias for a range of motifs and microsatellite lengths. Because ascertainment bias should not exist if a microsatellite selected in one species is as likely to be longer as it is to be shorter than its homologue, we propose that the nature of ascertainment bias can provide evidence for understanding how microsatellites evolve. We show that bias is greater for longer microsatellites but also that many long microsatellites have short homologues. These results are consistent with the notion that growth of long microsatellites is constrained by an upper length boundary that, when reached, sometimes results in large deletions. By evaluating ascertainment bias separately for interrupted and uninterrupted repeats we also show that long microsatellites tend to become interrupted, thereby contributing a second component of ascertainment bias. Having accounted for ascertainment bias, in agreement with results published elsewhere, we find that microsatellites in humans are longer on average than those in chimpanzees. This length difference is similar among repeat motifs but surprisingly comprises two roughly equal components, one associated with the repeats themselves and one with the flanking sequences. The differences we find can only be explained if microsatellites are both evolving directionally under a biased mutation process and are doing so at different rates in different closely related species.     

Global lack of flyway structure in a cosmopolitan bird revealed by a genome wide survey of single nucleotide polymorphisms

Knowledge about population structure and connectivity of waterfowl species, especially mallards (Anas platyrhynchos), is a priority because of recent outbreaks of avian influenza. Ringing studies that trace large‐scale movement patterns have to date been unable to detect clearly delineated mallard populations. We employed 363 single nucleotide polymorphism markers in combination with population genetics and phylogeographical approaches to conduct a population genomic test of panmixia in 801 mallards from 45 locations worldwide. Basic population genetic and phylogenetic methods suggest no or very little population structure on continental scales. Nor could individual‐based structuring algorithms discern geographical structuring. Model‐based coalescent analyses for testing models of population structure pointed to strong genetic connectivity among the world's mallard population. These diverse approaches all support the conclusion that there is a lack of clear population structure, suggesting that the world's mallards, perhaps with minor exceptions, form a single large, mainly interbreeding population.     

Analysis of population differentiation in North Eurasian cattle (Bos taurus) using single nucleotide polymorphisms in three genes associated with production traits

Single nucleotide polymorphisms (SNPs) in growth hormone 1 (GH1), insulin-like growth factor 1 (IGF1) and leptin (LEP), all candidates for production traits in cattle, were characterized in North Eurasian cattle breeds. Allele frequencies of IGF1 exhibited significant (P < 0.05) deviation from neutral expectation and therefore, might be associated with divergence in North Eurasian cattle because of genetic selection. Allele frequencies and lower heterozygosity of LEP may indicate a recent introduction of an alternative allele in this geographic region. Locus F(ST) estimates were highest for IGF1 (0.151, sigma = 0.042) and lowest for GH (0.062, sigma = 0.020). Our results suggest a slightly higher population differentiation across the candidate genes (FST = 0.108) than across microsatellites (FST = 0.095), possibly because of selection and stochastic effects.     

Functional impacts of non-synonymous single nucleotide polymorphisms: selective constraint and structural environments

In this work, we studied the correlations between selective constraint, structural environments and functional impacts of non-synonymous single nucleotide polymorphisms (nsSNPs). We found that the relation between solvent accessibility and functional impacts of nsSNPs is not as simple as generally thought. Finer structural classifications need to be taken into account to reveal the complex relations between the characteristics of a structure environment and its influence on the functional impacts of nsSNPs. We introduced two parameters for each structural environment, consensus residue percentage and residue distribution distance, to characterize the selective constraint imposed by the environment. Both parameters significantly correlate with the functional bias of nsSNPs across the structural environments. This result shows that selective constraint underlies the bias of a structural environment towards a certain type of nsSNPs (disease-associated or benign).     

Application of single nucleotide polymorphisms to non-model species: a technical review

Single nucleotide polymorphisms (SNPs) have gained wide use in humans and model species and are becoming the marker of choice for applications in other species. Technology that was developed for work in model species may provide useful tools for SNP discovery and genotyping in non-model organisms. However, SNP discovery can be expensive, labour intensive, and introduce ascertainment bias. In addition, the most efficient approaches to SNP discovery will depend on the research questions that the markers are to resolve as well as the focal species. We discuss advantages and disadvantages of several past and recent technologies for SNP discovery and genotyping and summarize a variety of SNP discovery and genotyping studies in ecology and evolution.     

Pig identification and meat traceability by multiallelic amplification fragments with multiple single nucleotide polymorphisms

Compared with conventional identification methods, DNA-based genetic approaches such as single nucleotide polymorphisms (SNPs) and satellites are much more reliable for pig identification and meat traceability. In this study, multiallelic amplification fragments with multiple SNPs, incorporating the advantages of both SNPs and microsatellites, were explored for the first time for pig identification and meat traceability. Primer pairs for multiallelic fragments and their optimal SNPs were successfully selected and used for identification of individuals from Suzhong and Duroc populations. Meanwhile, the combined panel of the above mentioned primer pairs together with their optimal SNPs for Suzhong and/or Duroc pigs were validated for identification of the hybrids (Suzhong×Duroc). Therefore, we have successfully selected multiallelic amplification fragments with multiple SNPs to identify pigs and their meat samples from Suzhong, Duroc or their hybrids. Our study demonstrates that our method is more powerful for pig identification or meat traceability than SNPs or microsatellites.     

Discovery,validation and characterization of 1039 cattle single nucleotide polymorphisms

We identified ～13 000 putative single nucleotide polymorphisms (SNPs) by comparison of repeat‐masked BAC‐end sequences from the cattle RPCI‐42 BAC library with whole‐genome shotgun contigs of cattle genome assembly Btau 1.0. Genotyping of a subset of these SNPs was performed on a panel containing 186 DNA samples from 18 cattle breeds including 43 trios. Of 1039 SNPs confirmed as polymorphic in the panel, 998 had minor allele frequency ≥0.25 among unrelated individuals of at least one breed. When Btau 4.0 became available, 974 of these validated SNPs were assigned in silico to known cattle chromosomes, while 41 SNPs were mapped to unassigned sequence scaffolds, yielding one SNP every ～3 Mbp on average. Twenty‐four SNPs identified in Btau 1.0 were not mapped to Btau 4.0. Of the 1015 SNPs mapped to Btau 4.0, 959 SNPs had nucleotide bases identical in Btau 4.0 and Btau 1.0 contigs, whereas 56 bases were changed, resulting in the loss of the in silico SNP in Btau 4.0. Because these 1039 SNPs were all directly confirmed by genotyping on the multi‐breed panel, it is likely that the original polymorphisms were correctly identified. The 1039 validated SNPs identified in this study represent a new and useful resource for genome‐wide association studies and applications in animal breeding.     

Using blocks of linked single nucleotide polymorphisms as highly polymorphic genetic markers for parentage analysis

Single nucleotide polymorphisms (SNPs) are plentiful in most genomes and amenable to high throughput genotyping, but they are not yet popular for parentage or paternity analysis. The markers are bi-allelic, so individually they contain little information about parentage, and in nonmodel organisms the process of identifying large numbers of unlinked SNPs can be daunting. We explore the possibility of using blocks of between three and 26 linked SNPs as highly polymorphic molecular markers for reconstructing male genotypes in polyandrous organisms with moderate (five offspring) to large (25 offspring) clutches of offspring. Haplotypes are inferred for each block of linked SNPs using the programs Haplore and Phase 2.1. Each multi-SNP haplotype is then treated as a separate allele, producing a highly polymorphic, 'microsatellite-like' marker. A simulation study is performed using haplotype frequencies derived from empirical data sets from Drosophila melanogaster and Mus musculus populations. We find that the markers produced are competitive with microsatellite loci in terms of single parent exclusion probabilities, particularly when using six or more linked SNPs to form a haplotype. These markers contain only modest rates of missing data and genotyping or phasing errors and thus should be seriously considered as molecular markers for parentage analysis, particularly when the study is interested in the functional significance of polymorphisms across the genome.     

Characterization of 59 canine single nucleotide polymorphisms in the Italian wolf (Canis lupus) population

We characterized 59 canine single nucleotide polymorphisms (SNPs) in the endangered Italian wolf (Canis lupus) population, which were discovered by resequencing sequence‐tagged‐site (STS) DNA sequences that are known to contain SNPs in domestic dogs. Dog SNPs were usually found also in wolves. Additional SNPs unique in dogs or wolves were discovered, which is important for detecting hybrids between dogs and wolves. We developed new primer sets and analysed 15 SNPs by Pyrosequencing. The characterized SNPs will provide an important addition to the genetic markers that are currently available for studying wild populations of canids.     

Characterization of 12 single nucleotide polymorphisms in weathervane scallop

We describe 27 single nucleotide polymorphisms (SNPs) in a commercially important bivalve, the weathervane scallop (Patinopecten caurinus), identified using a targeted‐gene approach. We further characterize 12 of these using 5′‐nuclease and allele‐specific PCR assays. Polymorphisms were identified in both mitochondrial and nuclear genes. These are the first SNPs developed for delineating population structure in the weathervane scallop and will provide a useful complement to currently available genetic markers.     

Comparison of SNPs and microsatellites for assessing the genetic structure of chicken populations

Many studies in human genetics compare informativeness of single‐nucleotide polymorphisms (SNPs) and microsatellites (single sequence repeats; SSR) in genome scans, but it is difficult to transfer the results directly to livestock because of different population structures. The aim of this study was to determine the number of SNPs needed to obtain the same differentiation power as with a given standard set of microsatellites. Eight chicken breeds were genotyped for 29 SSRs and 9216 SNPs. After filtering, only 2931 SNPs remained. The differentiation power was evaluated using two methods: partitioning of the Euclidean distance matrix based on a principal component analysis (PCA) and a Bayesian model‐based clustering approach. Generally, with PCA‐based partitioning, 70 SNPs provide a comparable resolution to 29 SSRs. In model‐based clustering, the similarity coefficient showed significantly higher values between repeated runs for SNPs compared to SSRs. For the membership coefficients, reflecting the proportion to which a fraction segment of the genome belongs to the ith cluster, the highest values were obtained for 29 SSRs and 100 SNPs respectively. With a low number of loci (29 SSRs or ≤100 SNPs), neither marker types could detect the admixture in the Gödöllö Nhx population. Using more than 250 SNPs allowed a more detailed insight into the genetic architecture. Thus, the admixed population could be detected. It is concluded that breed differentiation studies will substantially gain power even with moderate numbers of SNPs.     

宁夏回族原发性膝骨性关节炎与瘦素受体基因多态性的相关性

为探讨宁夏回族原发性膝骨性关节炎(Osteoarthritis, OA)与瘦素受体基因(Leptin receptor, LEPR)A668G位点单核苷酸多态性(SNPs)之间的关系, 文章运用病例-对照研究, 通过聚合酶链反应-限制性片段长度多态性(PCR-RFLP)技术, 对148例兼具原发症状和影象证据的宁夏回族膝OA患者以及155名年龄、性别相匹配的对照群体进行LEPR A668G SNPs检测, 并进行测序验证, 分析LEPR基因多态性与膝OA的易感关联。研究表明, 膝OA组瘦素(Leptin, LEP)水平显著高于对照组(P<0.001), 血浆可溶性瘦素受体(sLEPR)水平较对照组明显降低(P<0.001), 膝OA组LEPR A668G位点AG/GA+GG基因型和G等位基因的分布频率和对照组存在差异(P=0.008和P=0.024)。研究结果提示, LEPR A668G位点的多态性可能与宁夏回族人群中膝OA易感性相关, 可以作为预测宁夏回族膝OA发病危险的遗传标记及早期防治的候选基因之一。     

Genome‐wide single nucleotide polymorphisms reveal population history and adaptive divergence in wild guppies

Adaptation of guppies (Poecilia reticulata) to contrasting upland and lowland habitats has been extensively studied with respect to behaviour, morphology and life history traits. Yet population history has not been studied at the whole‐genome level. Although single nucleotide polymorphisms (SNPs) are the most abundant form of variation in many genomes and consequently very informative for a genome‐wide picture of standing natural variation in populations, genome‐wide SNP data are rarely available for wild vertebrates. Here we use genetically mapped SNP markers to comprehensively survey genetic variation within and among naturally occurring guppy populations from a wide geographic range in Trinidad and Venezuela. Results from three different clustering methods, Neighbor‐net, principal component analysis (PCA) and Bayesian analysis show that the population substructure agrees with geographic separation and largely with previously hypothesized patterns of historical colonization. Within major drainages (Caroni, Oropouche and Northern), populations are genetically similar, but those in different geographic regions are highly divergent from one another, with some indications of ancient shared polymorphisms. Clear genomic signatures of a previous introduction experiment were seen, and we detected additional potential admixture events. Headwater populations were significantly less heterozygous than downstream populations. Pairwise F_ST values revealed marked differences in allele frequencies among populations from different regions, and also among populations within the same region. F_ST outlier methods indicated some regions of the genome as being under directional selection. Overall, this study demonstrates the power of a genome‐wide SNP data set to inform for studies on natural variation, adaptation and evolution of wild populations     

Analysis of public single nucleotide polymorphisms in commercial pig populations

More than 5500 pig single nucleotide polymorphisms (SNPs) were recently identified and deposited in the public domain. To test the usefulness of these public SNPs, 109 SNPs were analysed for polymorphism within six commercial pig populations. A functional polymerase chain reaction (PCR) assay was obtained for 103 SNPs and it was possible to validate c. 59% by PCR-restriction fragment length polymorphism. Furthermore, polymorphism was found using a relatively limited number of genomic DNA samples, indicating that these polymorphisms are segregating at a useful frequency in these populations. The high percentage of validated markers demonstrates the utility of these public pig SNPs to identify loci responsible for economically important traits in commercial pig populations.     

Comparison of information content for microsatellites and SNPs in poultry and cattle

Data were available for 12 poultry microsatellites and 29 poultry single nucleotide polymorphisms (SNPs), and for 34 cattle microsatellites and 36 cattle SNPs. Stochastic permutation was used to determine the number of SNPs needed to obtain the same average information content as a given number of microsatellites. For poultry, the information content averaged 0.71 for the 12 microsatellites compared to 0.72 for the 29 SNPs. For cattle, the information content averaged 0.92 for the 34 microsatellites compared with 0.79 for the 36 SNPs. This study shows that, for each microsatellite, three SNPs are needed to obtain the same average information content.     

Association between single nucleotide polymorphisms of TPH1 and TPH2 genes,and depressive disorders

Tryptophan catabolites pathway disorders are observed in patients with depression. Moreover, single nucleotide polymorphisms of tryptophan hydroxylase genes may modulate the risk of depression occurrence. The objective of our study was to confirm the association between the presence of polymorphic variants of TPH1 and TPH2 genes, and the development of depressive disorders. Six polymorphisms were selected: c.804‐7C>A (rs10488682), c.‐1668T>A (rs623580), c.803+221C>A (rs1800532), c.‐173A>T (rs1799913)—TPH1, c.‐1449C>A (rs7963803), and c.‐844G>T (rs4570625)—TPH2. A total of 510 DNA samples (230 controls and 280 patients) were genotyped using TaqMan probes. Among the studied polymoorphisms, the G/G genotype and G allele of c.804‐7C>A—TPH1, the T/T homozygote of c.803+221C>A—TPH1, the A/A genotype and A allele of c.1668T>A—TPH1, the G/G homozygote and G allele of c.‐844G>T—TPH2, and the C/A heterozygote and A allele of c.‐1449C>A—TPH2 were associated with the occurrence of depression. However, the T/T homozygote of c.‐1668T>A—TPH1, the G/T heterozygote and T allele of c.‐844G>T—TPH2, and the C/C homozygote and C allele of c.‐1449C>A—TPH2 decreased the risk of development of depressive disorders . Each of the studied polymorphisms modulated the risk of depression for selected genotypes and alleles. These results support the hypothesis regarding the involvement of the pathway in the pathogenesis of depression.     

Discovery and evaluation of single nucleotide polymorphisms (SNPs) for Haliotis midae: a targeted EST approach

In this study, we describe the first set of SNP markers for the South African abalone, Haliotis midae. A cDNA library was constructed from which ESTs were selected for the screening of SNPs. The observed frequency of SNPs in this species was estimated at one every 185 bp. When characterized in wild-caught abalone, the minor allele frequencies and F(ST) estimates for every SNP indicated that these markers may potentially be useful for population analysis, parentage assignment and linkage mapping in Haliotis midae. No linkage disequilibrium was observed between SNPs originating from different EST sequences. These SNPs, together with additional SNPs currently being developed, will provide a useful complementary set of markers to the currently available genetic markers in abalone.