How many SNPs are enough?

Since the days of allozyme analysis, we have been enamored with the idea that if we just had enough polymorphic mendelian loci, we could gauge the inbreeding level of individuals by measuring heterozygosity and simultaneously measure the degree of genetic relatedness between pairs of individuals. Given Mendel’s Laws, we have always known that we would need numerous independently segregating loci to achieve any reasonable degree of accuracy. Santure et al. (2010, this issue) use a 771 marker SNP panel to assess heterozygosity levels and to assess pairwise relatedness, and compare both with theoretical expectations obtained from a carefully recorded pedigree of a zebra finch breeding colony, as a function of increasing numbers of SNP markers. They also compare the SNP results with those from a 20‐locus microsatellite panel, showing that adding SNPs to a fairly large microsatellite panel improves accuracy, but given an existing panel of 125 SNPs, little is to be gained by adding microsatellites. They show that the accuracy available for estimating individual levels of inbreeding is somewhat limited. They also show that the average pairwise relatedness measures bracket pedigree relationship very nicely, but the variances for individual pairs remain substantial, even with a very large panel.     

Isolation and characterization of microsatellite loci from forest musk deer (Moschus berezovskii)

This study reported the isolation and characterization of eight polymorphic microsatellite loci in endangered forest musk deer Moschus berezovskii. An improved enrichment protocol was used to isolate microsatellites, and polymorphism was explored with samples from wild musk deer population collected in Miyalo of Sichuan Province in China. Approximately 70% of clones from the genomic library constructed in current study contained dinucleotide (AC) repeats. Eight microsatellite loci amplified were highly polymorphic within forest musk deer population. The number of alleles per locus ranged from 6 to 14, and the observed and expected heterozygosities ranged from 0.41 approximately 1.0 and from 0.8 approximately 0.9, respectively. The average polymorphic information content (PIC) value for these markers was 0.82. This demonstrated that the eight microsatellite loci developed here are highly polymorphic, and can be used as genetic markers for further investigation of musk deer. Also, the results showed that the musk deer distributed in Miyalo had a relatively higher level of genetic variation.     

Genetic introgression between an introduced babbler, the Chinese hwamei Leucodioptron c. canorum, and the endemic Taiwan hwamei L. taewanus: a multiple marker systems analysis

Identification of interspecific hybrids is often a subject of primary concern in the development of conservation strategies. Here we performed a genetic analysis combining mitochondrial DNA (mtDNA), microsatellites and single nucleotide polymorphic sites (SNPs) to assay the level of hybridization and introgression between an introduced babbler, Chinese hwamei Leucodioptron canorus, and its close relative, the endemic Taiwan hwamei L. taewanus in Taiwan. Fifty‐five Chinese hwameis from the Asian mainland and 69 Taiwan hwameis, including nine morphological hybrids, were sampled and analyzed. Evidences of mitochondrial introgression were found in three hybrids and one Taiwan hwamei. Five unlinked interspecific SNPs were identified at nine anonymous nuclear loci with interspecific differentiation (total Fst=0.77) that was much higher than that at seven highly polymorphic microsatellite loci combined (total Fst=0.1). Bayesian cluster analysis based on five interspecific SNP loci and two highly differentiated microsatellite loci (Fst>0.08) suggested that twelve individuals sampled in Taiwan were likely F2 or backcross hybrids, among which eight were morphological intermediates. A total of 20.3% (14/69) individuals sampled in Taiwan were suggested to be hybrids, suggesting that fitness reduction in hybrids might be negligible. These results imply that without an effective management strategy, the entire Taiwan hwamei population could easily become an admixed with Chinese hwamei and loose its evolutionary integrity. To reduce introgressive hybridization, illegal trade of Chinese hwamei should be strictly regulated and only the expensive male Chinese hwameis should be legally imported to minimize the chance for Chinese hwameis being released into the field. In our study we also found interspecific SNP markers to outperform microsatellite loci in detecting hybridization and introgression between two closely related species, which may be ascribed to the lower level of homoplasy of SNP loci.     

Characterization of 151 SNPs for population structure analysis of the endangered Tatra chamois (Rupicapra rupicapra tatrica) and its relative, the Alpine chamois (R. r. rupicapra)

We present 151 SNPs (single nucleotide polymorphisms) identified in Tatra and Alpine chamois individuals after genotyping with the Bovine SNP50 Genotyping BeadChip (Illumina^®). Population structure analysis based on the identified 151 SNPs as well as a subset of 48 SNPs were able to give information about geographic origin of each individual making the markers suitable for future evaluation and monitoring of the genetic status of the Tatra and Alpine chamois. Our study demonstrates how it is possible to quickly identify informative SNPs in non-model organisms based on a SNP marker panel created for a related domestic species.     

Nine novel microsatellite markers for the hog deer (Axis porcinus)

Hog deer (Axis porcinus) is listed as an endangered species in many countries and its taxonomic status remains ambiguous. In this study, we developed nine novel microsatellite markers using enrichment methods. These polymorphic microsatellite loci showed allele numbers of 2–3 with an average of 2.22 in a group of 26 individuals and observed and expected heterozygosity ranging from 0.044 to 0.619 (average 0.397) and 0.194 to 0.632 (average 0.433) respectively. Linkage disequilibrium was not detected. These markers could be used in the genetic study and taxonomic identification of this species.     

Analysis of microsatellite loci in tree of heaven (<Emphasis Type="Italic">Ailanthus altissima</Emphasis> (Mill.) Swingle) using SSR-GBS

Microsatellite markers are still the marker of choice for many research questions in the field of forest genetics. However, the number of available markers is often low for species that have not been studied intensively like the tree of heaven (Ailanthus altissima). During the last decade, next-generation sequencing (NGS) has offered advanced techniques for efficiently identifying microsatellite markers and accurately genotyping samples. Here, we identify new microsatellite markers for the tree of heaven by applying an NGS-based method using the Illumina MiSeq platform. NGS technology was proved to be an effective method for fast and cost-efficient identification of microsatellite markers by implementing a genotyping-by-sequencing approach based on Illumina amplicon sequencing (SSR-GBS). We screened three populations from Eastern Austria for genetic variation at 19 newly identified microsatellite loci. We tested two different genotyping approaches: (1) considering only allele lengths (forming a so-called “allele length dataset”), (2) taking also single nucleotide polymorphisms (SNPs) within the amplified fragments into account (forming a so-called “SNP dataset”). The results revealed higher values for all genetic diversity parameters, as well as a better resolution of genetic assignment, when the latter approach was followed. Thus, by taking advantage of sequence information which is provided by SSR-GBS, one may achieve considerable gains in performance using the same marker set. The developed markers provide a cost-efficient tool for genotyping populations of tree of heaven and the approach presented here promises to be of high value for medium throughput genotyping applications in non-model forest tree species. We will use this method to widen the perspectives for further population genetic investigations of the tree of heaven.     

Genomic heterogeneity in the density of noncoding single-nucleotide and microsatellite polymorphisms in Plasmodium falciparum

The density and distribution of single-nucleotide polymorphisms (SNPs) across the genome has important implications for linkage disequilibrium mapping and association studies, and the level of simple-sequence microsatellite polymorphisms has important implications for the use of oligonucleotide hybridization methods to genotype SNPs. To assess the density of these types of polymorphisms in P. falciparum, we sampled introns and noncoding DNA upstream and downstream of coding regions among a variety of geographically diverse parasites. Across 36,229 base pairs of noncoding sequence representing 41 genetic loci, a total of 307 polymorphisms including 248 polymorphic microsatellites and 39 SNPs were identified. We found a significant excess of microsatellite polymorphisms having a repeat unit length of one or two, compared to those with longer repeat lengths, as well as a nonrandom distribution of SNP polymorphisms. Almost half of the SNPs localized to only three of the 41 genetic loci sampled. Furthermore, we find significant differences in the frequency of polymorphisms across the two chromosomes (2 and 3) examined most extensively, with an excess of SNPs and a surplus of polymorphic microsatellites on chromosome 3 as compared to chromosome 2 (P=0.0001). Furthermore, at some individual genetic loci we also find a nonrandom distribution of polymorphisms between coding and flanking noncoding sequences, where completely monomorphic regions may flank highly polymorphic genes. These data, combined with our previous findings of nonrandom distribution of SNPs across chromosome 2, suggest that the Plasmodium falciparum genome may be a mosaic with regard to genetic diversity, containing chromosomal regions that are highly polymorphic interspersed with regions that are much less polymorphic.     

Cross-amplification tests of ungulate primers in the endangered Neotropical pampas deer (Ozotoceros bezoarticus)

In cross-species amplification tests of 15 ungulate primers in pampas deer, five were retained to form a small panel of highly polymorphic loci that could be used to efficiently screen populations of this endangered species. The polymerase chain reactions were performed incorporating the universal fluorescent labeled M13 (-21) primer. In 69 pampas deer, average allelic diversity was 15, expected heterozygosity was 0.869 and the mean polymorphic information content value was 0.847. Paternity exclusion probabilities over loci were NE-1P = 0.01336 and NE-2P = 0.00135, and combined non-exclusion probability of identity was P(ID) = 3 x 10(-8).     

Isolation and characterization of SNP variation at 90 anonymous loci in the banded wren (<Emphasis Type="Italic">Thryothorus pleurostictus)</Emphasis>

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 (51 kb total) from the banded wren (Thryothorus pleurostictus), a Central American species with minimal pre-existing sequence data and a documented paucity of microsatellite allelic variation. 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.     

Characterization of sequence polymorphisms from microsatellite flanking regions in <Emphasis Type="Italic">Vitis</Emphasis> spp

Single nucleotide polymorphisms or SNPs are the most abundant form of genetic variation in the genome of plants and animals. Microsatellites are hypervariable regions of genome, while their flanking regions are assumed to be as conserved as the average of the genome. In the present study, flanking sequences of 10 microsatellite loci were compared in different cultivars of Vitis to determine the existing polymorphism. For every microsatellite, about 8 homozygous cultivars (regarding the microsatellite genotype) were chosen for sequencing. A total of 45 different varieties of Vitis and 91 sequences were analysed. Sequence polymorphisms were detected for all the microsatellite flanking regions studied, including single nucleotide polymorphisms (SNPs), insertions and deletions. The number of identified changes varied considerably among the loci with a frequency of one polymorphism every 41 nucleotides, being VVMD5 the most polymorphic one. A number of SNPs were used to design SNP markers, which were scored by dideoxy single base primer extension and capillary electrophoresis methodology. These SNP markers were employed to genotype 21 cultivars of Vitis vinifera and 4 varieties of other Vitis species. The utility of the markers developed as well as their utility for varietal identification and pedigree studies is discussed, using a similar study carried out with the 10 microsatellites as a reference.     

Nuclear microsatellite loci for Arabidopsis halleri (Brassicaceae), a model species to study plant adaptation to heavy metals

? Premise of the study: Arabidopsis halleri is a model species to study the adaptation of plants to soils contaminated by zinc, cadmium, and lead. To provide a neutral genetic background with which adaptive genetic markers could be compared, we developed highly polymorphic neutral microsatellite markers. ? Methods and Results: Using a microsatellite-enriched library method, we identified 120 microsatellite loci for quantitative trait locus (QTL) mapping analysis, of which eight primer pairs were developed in a single multiplex for population genetic studies. Analyses were performed on 508 individuals from 26 populations. All loci were polymorphic with six to 23 alleles per locus. Genetic diversity varied between 0.56 and 0.76. ? Conclusions: Our results demonstrated the value of these eight microsatellite markers to investigate neutral population genetic structure in A. halleri. To increase the resolution of population genetic analyses, we suggest adding them to the 11 markers previously developed independently.     

A panel of microsatellite loci for population studies of Sakhalin taimen Parahucho perryi (Brevoort)

Using DNA samples of Sakhalin taimen and a set of microsatellite loci, earlier reported for other salmonid fishes (Salmonidae), successful cross-species amplification was performed. A total of 56 Sakhalin taimen (Parahucho perryi) samples from the Daga, Nabil, Poronai, and Agnevo rivers (Sakhalin Island) were examined at 36 microsatellite loci, most of which were described for other species and first tested in taimen. Among the 21 loci first tested in taimen, two loci produced no amplification products. The remaining 19 loci were successfully amplified (for some loci, new primers were generated). Thirteen of these loci were monomorphic, while six loci were polymorphic and used in further population genetic analysis. In addition, with the purpose of modification of the allele sizes and optimization of the research technique, new primers for the already known 12 loci of Sakhalin taimen were designed. Three more loci were included in analysis without changes. As a result, a joint panel consisting of 21 polymorphic microsatellite markers was suggested for analysis of Sakhalin taimen. This panel was tested with four population samples from the rivers of Sakhalin Island.The results showed that this panel of markers could be used in detailed population studies for evaluation of the level of genetic differentiation, inbreeding, and migrations in Sakhalin taimen, an endangered species with a fragmented range. Using this approach in further studies will make it possible to isolate basic populations, which is necessary for conservation of this rare species.     

Mitochondrial DNA and microsatellite analyses of the genetic status of the presumed subspecies Cervus elaphus montanus (Carpathian red deer)

The possibly distinct Carpathian red deer was compared genetically to other European populations. We screened 120 red deer specimens from Serbia, the Romanian lowland and the Romanian Carpathians for genetic variability using 582 bp of the mitochondrial control region and nine polymorphic nuclear microsatellite loci. The study aimed at a population genetic characterization of the Carpathian red deer, which are often treated as a distinct subspecies (Cervus elaphus montanus). The genetic integrity of the Carpathian populations was confirmed through the haplotype distribution, private alleles and genetic distances. The Carpathian red deer are thus identified as one of the few remaining natural populations of this species, deserving special attention among game and conservation biologists. The history of the populations studied, in particular the introduction of Carpathian red deer into Romanian lowland areas in the 20th century, was reflected by the genetic data.     

Genome-wide polymorphism and comparative analyses in the white-tailed deer (Odocoileus virginianus): a model for conservation genomics

The white-tailed deer (Odocoileus virginianus) represents one of the most successful and widely distributed large mammal species within North America, yet very little nucleotide sequence information is available. We utilized massively parallel pyrosequencing of a reduced representation library (RRL) and a random shotgun library (RSL) to generate a complete mitochondrial genome sequence and identify a large number of putative single nucleotide polymorphisms (SNPs) distributed throughout the white-tailed deer nuclear and mitochondrial genomes. A SNP validation study designed to test specific classes of putative SNPs provides evidence for as many as 10,476 genome-wide SNPs in the current dataset. Based on cytogenetic evidence for homology between cow (Bos taurus) and white-tailed deer chromosomes, we demonstrate that a divergent genome may be used for estimating the relative distribution and density of de novo sequence contigs as well as putative SNPs for species without draft genome assemblies. Our approach demonstrates that bioinformatic tools developed for model or agriculturally important species may be leveraged to support next-generation research programs for species of biological, ecological and evolutionary importance. We also provide a functional annotation analysis for the de novo sequence contigs assembled from white-tailed deer pyrosequencing reads, a mitochondrial phylogeny involving 13,722 nucleotide positions for 10 unique species of Cervidae, and a median joining haplotype network as a putative representation of mitochondrial evolution in O. virginianus. The results of this study are expected to provide a detailed template enabling genome-wide sequence-based studies of threatened, endangered or conservationally important non-model organisms.     

SNPs by AFLP (SBA): a rapid SNP isolation strategy for non-model organisms

Despite the great potential of single nucleotide polymorphism (SNP) markers in evolutionary studies, in particular for inferring population genetic parameters, SNP analysis has almost exclusively been limited to humans and ‘genomic model’ organisms, due to the lack of available sequence data in non-model organisms. Here, we describe a rapid and cost effective method to isolate candidate SNPs in non-model organisms. This SNP isolation strategy consists basically in the direct sequencing of amplified fragment length polymorphism bands. In a first application of this method, 10 unique DNA fragments that contained 24 SNPs were discovered in 11.11 kb of sequenced genomic DNA of a non-model species, the brown trout (Salmo trutta).     

Evaluation of linkage disequilibrium and its effect on non-parametric multipoint linkage analysis using two high density single-nucleotide polymorphism mapping panels

Genotype data from the Illumina Linkage III SNP panel (n = 4,720 SNPs) and the Affymetrix 10 k mapping array (n = 11,120 SNPs) were used to test the effects of linkage disequilibrium (LD) between SNPs in a linkage analysis in the Collaborative Study on the Genetics of Alcoholism pedigree collection (143 pedigrees; 1,614 individuals). The average r2 between adjacent markers across the genetic map was 0.099 +/- 0.003 in the Illumina III panel and 0.17 +/- 0.003 in the Affymetrix 10 k array. In order to determine the effect of LD between marker loci in a nonparametric multipoint linkage analysis, markers in strong LD with another marker (r2 > 0.40) were removed (n = 471 loci in the Illumina panel; n = 1,804 loci in the Affymetrix panel) and the linkage analysis results were compared to the results using the entire marker sets. In all analyses using the ALDX1 phenotype, 8 linkage regions on 5 chromosomes (2, 7, 10, 11, X) were detected (peak markers p < 0.01), and the Illumina panel detected an additional region on chromosome 6. Analysis of the same pedigree set and ALDX1 phenotype using short tandem repeat markers (STRs) resulted in 3 linkage regions on 3 chromosomes (peak markers p < 0.01). These results suggest that in this pedigree set, LD between loci with spacing similar to the SNP panels tested may not significantly affect the overall detection of linkage regions in a genome scan. Moreover, since the data quality and information content are greatly improved in the SNP panels over STR genotyping methods, new linkage regions may be identified due to higher information content and data quality in a dense SNP linkage panel.     

Genetic Variations Revealed by Microsatellite Markers ina Small Population of the Sika Deer (Cervus nippon)on Kinkazan Island, Northern Japan

Genetic variations within a population of the Japanese sika deer, Cervus nippon, on Kinkazan Island were studied by microsatellite analysis. Seventeen pairs of polymerase chain reaction primers designed for several species of ungulates successfully amplified polymorphic microsatellite DNA in sika deer. About 20% of the Kinkazan population was sampled and genotyped for nine diagnostic microsatellite loci. Alleles at those loci in the Kinkazan population were found to be under the Hardy-Weinberg equilibrium. To determine whether the Kinkazan deer have a reduced level of genetic variability, an average heterozygosity in the population was calculated and compared with the values determined for other populations from Hyogo, Yamaguchi, Shimane, Tsushima, and Nagasaki. Neither the observed nor the expected heterozygosity in the Kinkazan deer significantly differed from that in the other populations. Our result indicated that, despite its small population size, the Kinkazan deer preserve extensive microsatellite variations.     

Bovine microsatellite loci are highly conserved in red deer (Cervus elaphus), sika deer (Cervus nippon) and Soay sheep (Ovis aries)

We tested 174 bovine microsatellite primer pairs for use in a primitive breed of sheep and two species of deer. Of 173 markers, 127 (73·4%) gave a product in Soay sheep ( Ovis aries ) of which 54 (42·5%) were polymorphic. One hundred and twenty-nine of 174 (74·1%) markers gave a product in red deer ( Cervus elaphus ) of which 72 (55·8%) were polymorphic. In sika deer ( Cervus nippon ) 126 of 171 (73·7%) microsatellite primers gave a product with 47 (37·3%) polymorphic. The proportion of bovine microsatellite loci conserved across artiodactyl species was significantly greater in this study than previously reported. Reasons for this high degree of microsatellite conservation are discussed. We suggest that a high resolution comparative map of the artiodactyls can be constructed using microsatellites.     

Application of high-resolution DNA melting for genotyping in lepidopteran non-model species: Ostrinia furnacalis (Crambidae)

Development of an ideal marker system facilitates a better understanding of the genetic diversity in lepidopteran non-model organisms, which have abundant species, but relatively limited genomic resources. Single nucleotide polymorphisms (SNPs) discovered within single-copy genes have proved to be desired markers, but SNP genotyping by current techniques remain laborious and expensive. High resolution melting (HRM) curve analysis represents a simple, rapid and inexpensive genotyping method that is primarily confined to clinical and diagnostic studies. In this study, we evaluated the potential of HRM analysis for SNP genotyping in the lepidopteran non-model species Ostrinia furnacalis (Crambidae). Small amplicon and unlabeled probe assays were developed for the SNPs, which were identified in 30 females of O. furnacalis from 3 different populations by our direct sequencing. Both assays were then applied to genotype 90 unknown female DNA by prior mixing with known wild-type DNA. The genotyping results were compared with those that were obtained using bi-directional sequencing analysis. Our results demonstrated the efficiency and reliability of the HRM assays. HRM has the potential to provide simple, cost-effective genotyping assays and facilitates genotyping studies in any non-model lepidopteran species of interest.