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.     

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.     

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.     

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.     

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.     

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.     

Relationship between single nucleotide polymorphism of interleukin-18 and susceptibility to pulmonary tuberculosis in the Chinese Han population

Interleukin-18 (IL-18) is a multi-functional cytokine capable of inducing either Th1 or Th2 polarization depending on the immunologic milieu. IL-18 may influence the host response to Mycobacterium tuberculosis (M.tb) infection. To investigate the relationship between single nucleotide polymorphisms of the IL-18 and susceptibility to pulmonary tuberculosis in the Chinese Han population, the IL-18 gene was sequenced to detect polymorphisms and to examine the genotype frequencies in 300 patients and 702 healthy controls. DNA sequencing revealed three IL-18 variants: rs1946518, rs5744247, and rs549908. It also revealed that allele A of rs1946518 confers a 1.47-fold increased risk of developing tuberculosis (TB) (P = 0.0001, OR [95%CI] = 1.47 [1.21-1.78]), and that the C allele of rs5744247 confers a 0.77-fold decreased risk of disease (P = 0.01, R [95%CI] = 0.77 [0.632-0.937]). The genotypes rs1946518, rs5744247 and rs549908 were found to be significantly associated with TB. Estimation of the frequencies of haplotypes revealed a potential risk haplotype AGA (P = 0.01, OR [95%CI] = 1.41 [1.15-1.72]) and a protective haplotype CCA (P = 0.01, OR [95%CI] = 0.70 [0.57-0.85]) for TB. The present findings suggest that polymorphisms in the IL-18 gene may affect susceptibility to TB and increase the risk of developing the disease in the Chinese Han population.     

Genetic association of single nucleotide polymorphisms in dystrobrevin binding protein 1 gene with schizophrenia in a Malaysian population

Dystrobrevin binding protein 1 (DTNBP1) gene is pivotal in regulating the glutamatergic system. Genetic variants of the DTNBP1 affect cognition and thus may be particularly relevant to schizophrenia. We therefore evaluated the association of six single nucleotide polymorphisms (SNPs) with schizophrenia in a Malaysian population (171 cases; 171 controls). Associations between these six SNPs and schizophrenia were tested in two stages. Association signals with p < 0.05 and minor allele frequency > 0.05 in stage 1 were followed by genotyping the SNPs in a replication phase (stage 2). Genotyping was performed with sequenced specific primer (PCR-SSP) and restriction fragment length polymorphism (PCR-RFLP). In our sample, we found significant associations between rs2619522 (allele p = 0.002, OR = 1.902, 95%CI = 1.266 – 2.859; genotype p = 0.002) and rs2619528 (allele p = 0.008, OR = 1.606, 95%CI = 1.130 – 2.281; genotype p = 6.18 × 10⁻⁵) and schizophrenia. Given that these two SNPs may be associated with the pathophysiology of schizophrenia, further studies on the other DTNBP1 variants are warranted.     

Polymerase chain reaction multiplexing of microsatellites and single nucleotide polymorphism markers for quantitative trait loci mapping of wild house mice

We tested 96 microsatellites and 10 single nucleotide polymorphisms for their allelic distribution in two subspecies of the house mouse, Mus musculus musculus and M. m. domesticus. Sixty‐two microsatellites discriminated strain‐specific differences among nine wild‐derived ‘musculus’ and ‘domesticus’ and three ‘classical’ laboratory strains. For efficient genotyping, we optimized multiplex conditions using five microsatellites per polymerase chain reaction. All 10 single nucleotide polymorphisms were also optimized for simultaneous analysis in one reaction using SNaPshot multiplex. The uniform distribution of markers on autosomes and on the X chromosome makes these panels potentially useful tools for quantitative trait loci mapping of wild house mice.     

Assessment of the population structure of five Finnish dog breeds with microsatellites

Genetic variabilities within and between Finnish populations of Golden Retrievers, German Shepherds, Wirehaired Dachshunds, Pembroke Welsh Corgis and Bedlington Terriers were quantified with microsatellite allele numbers, observed heterozygosities, expected heterozygosities, FIS estimates, FST estimates and DS distances. In a sample of 50 individuals from each, breed and ten polymorphic loci, the highest genetic diversity was exhibited in the Wirehaired Dachshunds and lowest in the Bedlington Terriers. Although statistically significant deviations from the Hardy-Weinberg (H-W) equilibrium were observed, they occurred at an unexpectedly low frequency. Most strikingly, the extremely small Bedlington Terrier population displayed genotypes in H-W proportions in all investigated loci. The H-W deviations always occurred with positive FIS estimates, which, on average, were not larger than values reported for free-living canids. Genetic differentiation between the breeds was very large. As a comparison, present estimates were, on average, over two times higher than previously observed between breeds of sheep, and over two times higher than the highest estimates reported between human populations. Moreover, the highest DS distances were only slightly lower than the lowest values inferred between humans and chimpanzees. Severe bottlenecks in the recent past of the examined breeds were not statistically supported. The presented data imply genetic isolation and intense artificial selection in the history of these breeds of dogs.     

Relationships among calpastatin single nucleotide polymorphisms,calpastatin expression and tenderness in pork longissimus1

Genome scans in the pig have identified a region on chromosome 2 (SSC2) associated with tenderness. Calpastatin is a likely positional candidate gene in this region because of its inhibitory role in the calpain system that is involved in postmortem tenderization. Novel single nucleotide polymorphisms (SNP) in calpastatin were identified and used to genotype a population (n = 1042) of Duroc–Landrace–Yorkshire swine for association with longissimus lumborum slice shear force (SSF) measured at days 7 and 14 postmortem. Three genetic markers residing in the calpastatin gene were significantly associated with SSF (P < 0.0005). Haplotypes constructed from markers in the calpastatin gene were significantly associated with SSF (F‐ratio = 3.93; P‐value = 0.002). The levels of normalized mRNA expression of calpastatin in the longissimus lumborum of 162 animals also were evaluated by real‐time RT‐PCR and were associated with the genotype of the most significant marker for SSF (P < 0.02). This evidence suggests that the causative variation alters expression of calpastatin, thus affecting tenderness. In summary, these data provide evidence of several significant, publicly available SNP markers associated with SSF that may be useful to the swine industry for marker assisted selection of animals that have more tender meat.     

Genotyping of five single nucleotide polymorphisms in the OCA2 and HERC2 genes associated with blue‐brown eye color in the Japanese population

Human eye color is a polymorphic phenotype influenced by multiple genes. It has recently been reported that three single nucleotide polymorphisms (SNPs) within intron 1 of the OCA2 gene (rs7495174, rs4778241, rs4778138) and two SNPs in intron 86 (rs12913832) and the 3′ UTR region (rs1129038) of the HERC2 gene—located in the upstream of the OCA2 locus —have a high statistical association with human eye color. The present study is the first to examine in detail the genotype and haplotype frequencies for these five SNPs in an Asian (Japanese) population (n = 523) comprising solely brown‐eyed individuals. Comparison of the genotype and haplotype distributions in Japanese with those in African and European subjects revealed significant differences between Japanese and other populations. Analysis of haplotypes consisting of four SNPs at the HERC2‐OCA2 locus (rs12913832/rs7495174/rs4778241/rs4778138) showed that the most frequent haplotype in the Japanese population is A‐GAG (0.568), while the frequency of this haplotype is rather low in the European population, even in the brown‐eyed group (0.167). The haplotype distribution in the Japanese population was significantly different from that in the brown‐eyed European group (F_ST = 0.18915). Copyright © 2009 John Wiley & Sons, Ltd.     

Development of new nuclear markers and characterization of single nucleotide polymorphisms in kelp gull (Larus dominicanus)

The present study seeks to develop nuclear markers for the kelp gull (Larus dominicanus). We hereby report the characterization of 12 independent nuclear introns, where 104 single nucleotide polymorphisms (SNPs) in 8138 sequenced base pairs were observed. These SNP markers are the first to be designed for genotyping a gull species. The markers will provide useful tools for understanding which processes act or acted upon kelp gulls to cause their low genetic variability in mitochondrial DNA. In addition, these markers open a new opportunity for population genetic and evolutionary studies in the Laridae group.