Isolation and characterization of polymorphic microsatellite loci for the three Iberian vipers, Vipera aspis, V. Latastei and V. seoanei by Illumina MiSeq sequencing.

BACKGROUND: European vipers (genus Vipera) are a well-studied taxonomic group, but the low resolution of nuclear sanger-sequenced regions has precluded thorough studies at systematic, ecological, evolutionary and conservation levels. In this study, we developed novel microsatellite markers for the three Iberian vipers, Vipera aspis, V. latastei and V. seoanei, and assessed their polymorphism in north-central Iberian populations. METHODS AND RESULTS: Genomic libraries were developed for each species using an Illumina Miseq sequencing approach. From the 70 primer pairs initially tested, 48 amplified reliably and were polymorphic within species. Cross-species transferability was achieved for 31 microsatellites loci in the three target species and four additional loci that were transferable to one species only. The 48 loci amplified in average seven alleles, and detected average expected and observed heterozygosities of 0.7 and 0.55, in the three genotyped populations/species (26 V. aspis, 20 V. latastei and 10 V. seoanei). CONCLUSIONS: Our study provides a selection of 48 polymorphic microsatellite markers that will contribute significantly to current knowledge on genetic diversity, gene flow, population structure, demographic dynamics, systematics, reproduction and heritability in these species, and potentially in other congeneric taxa.

Development of novel microsatellite markers for population differentiation and detection of natural selection in wild populations of butter catfish, Ompok bimaculatus (Bloch, 1794).

BACKGROUND: Butter catfish (Ompok bimaculatus) is a preferred species in South East Asia, with huge aquaculture potential. However, there is limited information about genetic stock composition due to insufficient markers. The goal of this study was to develop de novo microsatellite markers. METHODS AND RESULTS: For sequencing, genomic SMRT bell libraries (1.5 Kbp size) were prepared for O. bimaculatus. A total of 114 SSR containing sequences were used for primer designing. Polymorphic loci were validated by genotyping 83 individuals from four distant riverine populations, viz., Brahmaputra, Bichiya, Gomti and Kaveri. A total of 30 microsatellite loci were polymorphic, of which five were found to be associated with functional genes and eight (four positive and four negative) loci were found to be under selection pressure. A total of 115 alleles were detected in all loci and PIC ranged from 0.539 to 0.927 and pair-wise FST values from 0.1267 to 0.26002 (p < 0.001), with an overall FST value of 0.17047, indicating the presence of population sub-structure. Cross-species transferability of 29 loci (96.67%) was successful in congener species, Ompok pabda. CONCLUSION: The novel SSR markers developed in this study would facilitate stock characterization of natural populations, to be used in future selection breeding programs and planning conservation strategies in these species. Identified non-neutral markers will give insights into the effect of local adaptation on genetic differentiation in the natural population of this species.

Molecular Characterization of Tinospora cordifolia (Willd.) Miers Using Novel g-SSR Markers and Their Comparison with EST-SSR and SCoT Markers for Genetic Diversity Study.

In the present study, novel genomic-SSR (g-SSR) markers generated in our laboratory were used to characterize Tinospora cordifolia and related species. The g-SSR marker was also compared with EST-SSR and SCoT markers used earlier in our laboratory to assess the genetic diversity of T. cordifolia. A total of 26 accessions of T. cordifolia and 1 accession each of Tinospora rumphii and Tinospora sinensis were characterized using 65 novel g-SSR markers. A total of 125 alleles were detected with 49 polymorphic g-SSR markers. The number of alleles per locus varied from 1-4 with a mean value of 2.55 alleles per locus. Mean PIC, gene diversity, and heterozygosity were estimated to be 0.33, 0.41, and 0.65, respectively. The two species, namely T. rumphii and T. sinensis, showed cross-species transferability of g-SSRs developed in T. cordifolia. The success rate of cross-species transferability in T. rumphii was 95.3% and 93.8% in T. sinensis, proving the usefulness of this marker in genetic diversity studies of related species. The Tinospora accessions were also used for molecular characterization using SCoT and EST-SSR markers and compared for genetic diversity and cross-species transferability. The PIC, gene diversity, heterozygosity, and principal coordinate analysis showed that g-SSR is the better maker for a genetic diversity study of T. cordifolia. Additionally, high cross-species transferability of g-SSRs was found (95.3% and 93.8%) compared to EST-SSRs (68.8% and 67.7%) in T. rumphii and T. sinensis, respectively.

Identification and characterization of chickpea genotypes for early flowering and higher seed germination through molecular markers.

BACKGROUND: Chickpea is the fourth most important legume crop contributing 15.42% to the total legume production and a rich source of proteins, minerals, and vitamins. Determination of genetic diversity of wild and elite cultivars coupled with early flowering and higher seed germination lines are quintessential for variety improvement. METHODS AND RESULTS: In the present study, we have analyzed the genetic diversity, population structure, cross-species transferability, and allelic richness in 50 chickpea collections using 23 Inter simple sequence repeats (ISSR) markers. The observed parameters such as allele number varied from 3 to 16, range of allele size varied from 150 to 1600 bp and polymorphic information content (PIC) range lies in between 0.15 and 0.49. Dendrogram was constructed with ISSR marker genotypic data and classified 50 chickpea germplasms into groups I and II, where the accession P 74 - 1 is in group I and the rest are in group II. Dendrogram, Principal component analysis (PCA), dissimilarity matrix, and Bayesian model-based genetic clustering of 50 chickpea germplasms revealed that P 74 - 1 and P 1883 are very diverse chickpea accessions. CONCLUSION: Based on genetic diversity analysis, 15 chickpea germplasm having been screened for early flowering and higher seed germination and found that the P 1857-1 and P 3971 have early flowering and higher seed germination percentage in comparison to P 1883 and other germplasm. These agronomic traits are essential for crop improvement and imply the potential of ISSR markers in crop improvement.

Comparative Association Mapping Reveals Conservation of Major Gene Resistance to White Pine Blister Rust in Southwestern White Pine (Pinus strobiformis) and Limber Pine (P. flexilis).

All native North American white pines are highly susceptible to white pine blister rust (WPBR) caused by Cronartium ribicola. Understanding genomic diversity and molecular mechanisms underlying genetic resistance to WPBR remains one of the great challenges in improvement of white pines. To compare major gene resistance (MGR) present in two species, southwestern white pine (Pinus strobiformis) Cr3 and limber pine (P. flexilis) Cr4, we performed association analyses of Cr3-controlled resistant traits using single nucleotide polymorphism (SNP) assays designed with Cr4-linked polymorphic genes. We found that ∼70% of P. flexilis SNPs were transferable to P. strobiformis. Furthermore, several Cr4-linked SNPs were significantly associated with the Cr3-controlled traits in P. strobiformis families. The most significantly associated SNP (M326511_1126R) almost colocalized with Cr4 on the Pinus consensus linkage group 8, suggesting that Cr3 and Cr4 might be the same R locus, or have localizations very close to each other in the syntenic region of the P. strobiformis and P. flexilis genomes. M326511_1126R was identified as a nonsynonymous SNP, causing amino acid change (Val376Ile) in a putative pectin acetylesterase, with coding sequences identical between the two species. Moreover, top Cr3-associated SNPs were further developed as TaqMan genotyping assays, suggesting their usefulness as marker-assisted selection (MAS) tools to distinguish genotypes between quantitative resistance and MGR. This work demonstrates the successful transferability of SNP markers between two closely related white pine species in the hybrid zone, and the possibility for deployment of MAS tools to facilitate long-term WPBR management in P. strobiformis breeding and conservation.

Genome-wide microsatellites in amaranth: development, characterization, and cross-species transferability.

Amaranth (Amaranthus spp.) belonging to Amaranthaceae, is known as "the crop of the future" because of its incredible nutritional quality. Amaranthus spp. (> 70) have a huge diversity in terms of their plant morphology, production and nutritional quality; however, these species are not well characterized at molecular level due to unavailability of robust and reproducible molecular markers, which is essential for crop improvement programs. In the present study, 13,051 genome-wide microsatellite motifs were identified and subsequently utilized for marker development using A. hypochondriacus (L.) genome (JPXE01.1). Out of those, 1538 motifs were found with flanking sequences suitable for primer designing. Among designed primers, 225 were utilized for validation of which 119 (52.89%) primers were amplified. Cross-species transferability and evolutionary relatedness among ten species of Amaranthus (A. hypochondriacus, A. caudatus, A. retroflexus, A. cruentus, A. tricolor, A. lividus, A. hybridus, A. viridis, A. edulis, and A. dubius) were also studied using 45 microsatellite motifs. The maximum (86.67%) and minimum (28.89%) cross-species transferability were observed in A. caudatus and A. dubius, respectively, that indicated high variability present across the Amaranthus spp. Total 97 alleles were detected among 10 species of Amaranthus. The averages of major allele frequency, gene diversity, heterozygosity and PIC were 0.733, 0.347, 0.06, and 0.291, respectively. Nei's genetic dissimilarity coefficients ranged from 0.0625 (between A. tricolor and A. hybridus) to 0.7918 (between A. viridis and A. lividus). The phylogenetic tree grouped ten species into three major clusters. Genome-wide development of microsatellite markers and their transferability revealed relationships among amaranth species which ultimately can be useful for species identification, DNA fingerprinting, and QTLs/gene(s) identification. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s13205-021-02930-5.

Hybridization, characterization and transferability of SSRs in the genus Morchella.

Recently, Morchella importuna, M. sextelata, M. eximia, M. exuberans, Mel-13, and Mel-21 have been successfully cultivated in China and planting areas rapidly expanded because of their economic importance. Effective molecular markers are urgently needed for accurately identifying morel cultivars. Microsatellites are widely utilized for strain authentication in many fungal groups. To our knowledge, for the first time we characterized the distribution of microsatellites (simple sequence repeats, SSRs) in the M. importuna genome with 12902 SSRs and reported the first set of SSRs developed for Morchella species. Mono-nucleotides (66.2 %) were the most frequent motifs, followed by tri- (15.4 %), di- (12.1 %), tetra- (3.7 %), penta- (1.3 %) and hexa-nucleotides (1.3 %). We tested the cross-species amplification of 180 SSRs on 24 samples from the six species and high cross-species transferability of SSRs (87.7 %) was found. Among twenty-two microsatellites selected for genetic diversity analysis on 127 samples from the six species, fifteen to twenty polymorphic loci were identified in M. importuna, M. sextelata, M. eximia, M. exuberans, Mel-13 and Mel-21. Interspecific hybridization events were detected among morel species, indicating the potential application of morel crossbreeding. Ninety-one cultivated samples were characterized as new cultivars with different genotypes, but cultivar names used for these by farmers was confusing, with misnaming, synonyms and homonyms. Our results are not only helpful for cultivar identification and morel breeding programs in China, but also provide molecular tools for genetic studies in morels.

The potentiality of rice microsatellite markers in assessment of cross-species transferability and genetic diversity of rice and its wild relatives.

The main aim of this study is to assess the potentiality of SSR markers for the identification of the cross-species transferability frequency in a large set of the diverse genome types of wild relative rice along with cultivated rice. Here, we used 18 different rice genotypes representing nine different genome types with 70 SSR markers to investigate the potentiality of cross-species transferability rate. The overall cross-species transferability of SSR markers across the 18 rice genotypes ranged from 38.9% (RM280 and RM447) to 100% (RM490, RM318, RM279, RM18877 and RM20033, RM19303) with an average of 76.58%. Also, cross-species transferability across chromosome ranged from 54.4% (chromosome 4) to 86.5% (chromosome 2) with an average of 74.35%. The polymorphism information content of the markers varied from 0.198 (RM263) to 0.868 (RM510) with a mean of 0.549 ± 0.153, showing high discriminatory power. The highest rate of cross-transferability was observed in O. rufipogon (97%), The highest rate of cross-species transferability was in O. rufipogon (97.00%), followed by O. glaberrima (94.20%), O. nivara (92.80%), Swarna (92.80%), O. longistaminata (91.40%), O. eichingeri (90%), O. barthii (88.50%), O. alta (82.80%), O. australiensis (77.10%), O. grandiglumis (74.20%), O. officinalis (74.20%), Zizania latifolia (70.00%), O. latifolia (68.50%), O. brachyantha (62.80%), Leersia perrieri (57.10%) and O. ridleyi (41.40%) with least in O. coarctata (28.50%). A total of 341 alleles from 70 loci were detected with the number of alleles per locus ranged from 2 to 12. Based on dendrogram analysis, the AA genome groups was separated as distinct group from the rest of the genome types. Similarly, principal coordinate analysis and structure analysis clearly separated the AA genome type from the rest of the genome types. Through the analysis of molecular variance, more variance (51%) was observed among the individual, whereas less (14%) was observed among the population. Thus, our findings may offer a valuable resource for studying the genetic diversity and relationship to facilitate the understanding of the complex mechanism of the origin and evolutionary processes of different Oryza species and wild relative rice.

Microsatellite marker development by partial sequencing of the sour passion fruit genome (Passiflora edulis Sims).

BACKGROUND: The Passiflora genus comprises hundreds of wild and cultivated species of passion fruit used for food, industrial, ornamental and medicinal purposes. Efforts to develop genomic tools for genetic analysis of P. edulis, the most important commercial Passiflora species, are still incipient. In spite of many recognized applications of microsatellite markers in genetics and breeding, their availability for passion fruit research remains restricted. Microsatellite markers in P. edulis are usually limited in number, show reduced polymorphism, and are mostly based on compound or imperfect repeats. Furthermore, they are confined to only a few Passiflora species. We describe the use of NGS technology to partially assemble the P. edulis genome in order to develop hundreds of new microsatellite markers. RESULTS: A total of 14.11 Gbp of Illumina paired-end sequence reads were analyzed to detect simple sequence repeat sites in the sour passion fruit genome. A sample of 1300 contigs containing perfect repeat microsatellite sequences was selected for PCR primer development. Panels of di- and tri-nucleotide repeat markers were then tested in P. edulis germplasm accessions for validation. DNA polymorphism was detected in 74% of the markers (PIC = 0.16 to 0.77; number of alleles/locus = 2 to 7). A core panel of highly polymorphic markers (PIC = 0.46 to 0.77) was used to cross-amplify PCR products in 79 species of Passiflora (including P. edulis), belonging to four subgenera (Astrophea, Decaloba, Distephana and Passiflora). Approximately 71% of the marker/species combinations resulted in positive amplicons in all species tested. DNA polymorphism was detected in germplasm accessions of six closely related Passiflora species (P. edulis, P. alata, P. maliformis, P. nitida, P. quadrangularis and P. setacea) and the data used for accession discrimination and species assignment. CONCLUSION: A database of P. edulis DNA sequences obtained by NGS technology was examined to identify microsatellite repeats in the sour passion fruit genome. Markers were submitted to evaluation using accessions of cultivated and wild Passiflora species. The new microsatellite markers detected high levels of DNA polymorphism in sour passion fruit and can potentially be used in genetic analysis of P. edulis and other Passiflora species.

Genetic diversity and differentiation in reef-building Millepora species, as revealed by cross-species amplification of fifteen novel microsatellite loci.

Quantifying the genetic diversity in natural populations is crucial to address ecological and evolutionary questions. Despite recent advances in whole-genome sequencing, microsatellite markers have remained one of the most powerful tools for a myriad of population genetic approaches. Here, we used the 454 sequencing technique to develop microsatellite loci in the fire coral Millepora platyphylla, an important reef-builder of Indo-Pacific reefs. We tested the cross-species amplification of these loci in five other species of the genus Millepora and analysed its success in correlation with the genetic distances between species using mitochondrial 16S sequences. We succeeded in discovering fifteen microsatellite loci in our target species M. platyphylla, among which twelve were polymorphic with 2-13 alleles and a mean observed heterozygosity of 0.411. Cross-species amplification in the five other Millepora species revealed a high probability of amplification success (71%) and polymorphism (59%) of the loci. Our results show no evidence of decreased heterozygosity with increasing genetic distance. However, only one locus enabled measures of genetic diversity in the Caribbean species M. complanata due to high proportions of null alleles for most of the microsatellites. This result indicates that our novel markers may only be useful for the Indo-Pacific species of Millepora. Measures of genetic diversity revealed significant linkage disequilibrium, moderate levels of observed heterozygosity (0.323-0.496) and heterozygote deficiencies for the Indo-Pacific species. The accessibility to new polymorphic microsatellite markers for hydrozoan Millepora species creates new opportunities for future research on processes driving the complexity of their colonisation success on many Indo-Pacific reefs.

Twenty microsatellite markers for the endangered Vatica mangachapoi (Dipterocarpaceae).

PREMISE OF THE STUDY: Microsatellite markers were developed for Vatica mangachapoi (Dipterocarpaceae), an endangered species indigenous to Southeast Asia and southern China. METHODS AND RESULTS: Twenty microsatellite markers, including 12 polymorphic markers, were identified from V. mangachapoi using high-throughput sequencing. Polymorphism at each marker was evaluated using 87 individuals from three natural populations. The number of alleles per polymorphic locus ranged from six to 15, and the observed and expected heterozygosity varied from 0.000 to 0.926 and from 0.177 to 0.864, respectively. These markers were transferred successfully to the endangered species V. guangxiensis. CONCLUSIONS: These markers may be used to investigate the genetic diversity and gene flow of V. mangachapoi and V. guangxiensis.

Development and Characterization of Novel EST-SSRs from Larix gmelinii and Their Cross-Species Transferability.

A set of 899 L. gmelinii expression sequence tags (ESTs), available at the National Center of Biotechnology Information (NCBI), was employed to address the feasibility on development of simple sequence repeat (SSR) markers for Larch species. Totally, 634 non-redundant unigenes including 145 contigs and 489 singletons were finally identified and mainly involved in biosynthetic, metabolic processes and response to stress according to BLASTX results, gene ontology (GO) categories and Kyoto Encyclopedia of Genes and Genomes (KEGG) maps. Approximately 11.7% (74) unigenes contained 90 candidate SSRs, which were mainly trinucleotides (29, 32.2%) and dinucleotides (26, 28.9%). A relatively high frequency of SSRs was respectively found in the Open Reading Frame (ORF, about 54.4%) and 5'-untranslated region (5'-UTR, 31.2%), while a low frequency was observed in the 3'-untranslated region (3'-UTR, about 14.4%). Of the 45 novel EST-SSRs markers, nine were found to be polymorphic at two L. gmelinii populations. The number of alleles per locus (Na) ranged from two to four, and the observed (Ho) and expected (He) heterozygosity values were 0.200-0.733 and 0.408-0.604, respectively. The inbreeding coefficients (FIS) for all loci were more than zero except Lg41. Most of these 9EST-SSR markers were transferable to its related species L. kaempferi, L. principis-rupprechtii and L. olgensis. These novel EST-SSRs will be useful for further research on comparative genomics, genetic resources conservation and molecular breeding in larch trees.

Evolutionary factors affecting the cross-species utility of newly developed microsatellite markers in seabirds.

Microsatellite loci are ideal for testing hypotheses relating to genetic segregation at fine spatio-temporal scales. They are also conserved among closely related species, making them potentially useful for clarifying interspecific relationships between recently diverged taxa. However, mutations at primer binding sites may lead to increased nonamplification, or disruptions that may result in decreased polymorphism in nontarget species. Furthermore, high mutation rates and constraints on allele size may also with evolutionary time, promote an increase in convergently evolved allele size classes, biasing measures of interspecific genetic differentiation. Here, we used next-generation sequencing to develop microsatellite markers from a shotgun genome sequence of the sub-Antarctic seabird, the thin-billed prion (Pachyptila belcheri), that we tested for cross-species amplification in other Pachyptila and related sub-Antarctic species. We found that heterozygosity decreased and the proportion of nonamplifying loci increased with phylogenetic distance from the target species. Surprisingly, we found that species trees estimated from interspecific FST provided better approximations of mtDNA relationships among the studied species than those estimated using DC , even though FST was more affected by null alleles. We observed a significantly nonlinear second order polynomial relationship between microsatellite and mtDNA distances. We propose that the loss of linearity with increasing mtDNA distance stems from an increasing proportion of homoplastic allele size classes that are identical in state, but not identical by descent. Therefore, despite high cross-species amplification success and high polymorphism among the closely related Pachyptila species, we caution against the use of microsatellites in phylogenetic inference among distantly related taxa.

Microsatellite Development and Potential Application in Authentication, Conservation, and Genetic Improvement of Chinese Medicinal Plants / 中草药·英文版

Medicinal plants are popular and widely used as a major source of herbal drugs and pharmaceutical compounds. Ever-growing demands make medicinal plants faced to several problems including efficacy and safety to meet business requirements, conservation, and artificially assisted breeding. As a powerful molecular tool, microsatellites offer the great potentials for various purposes in plants. This review provides a scenario of microsatellites in medicinal plants including development from genomic or expressed sequence tag libraries, cross-species transferability, genotyping, and potential applications. We emphasized on the authentication of medicinal plants by microsatellite markers.

Characterization and high cross-species transferability of microsatellite markers from the floral transcriptome of Aspidistra saxicola (Asparagaceae).

Recent studies utilizing transcriptome sequences generated by next-generation sequencing (NGS) technologies have demonstrated the ability to rapidly detect and characterize thousands of gene-based microsatellites from different plants. However, these simple sequence repeats (SSRs) were seldom used directly to test interspecific transferability in the populations of closely related species. Aspidistra Ker-Gawl. is a monocot genus with high species richness and diversity in flower structure, but its fresh floral materials are not easy to obtain. Until now, little is known about genetic background in the species of Aspidistra, quite apart from the fearful reduction of their natural habitats. In this study, the floral transcriptome of Aspidistra saxicola was obtained using NGS. Based on these data, a total of 5527 SSRs were identified in the unigenes. Among these SSRs, the proportions of di- and tri-nucleotide repeats were quite close (49.6% verse 46.8%), and the most tri-nucleotide repeats were AGG/CCT followed by AAG/CTT and AGC/GCT in A. saxicola, showing distinct differences with other angiosperm species. To assess genetic diversity in the species of Aspidistra, 48 SSR loci were tested in four available populations of A. elatior. The results revealed that more than a third of the loci were polymorphic. The majority of these primers could be amplified in 24 species representing the main clades of Aspidistra. The primer subsets from transcriptome data proved highly useful for detecting polymorphisms in the related species, supporting the finding that NGS is an efficient approach to molecular marker development at both intra- and interspecies levels, especially in endangered nonmodel species.

Isolation, via 454 sequencing, and characterization of microsatellites for Vachellia farnesiana (Fabaceae: Mimosoideae).

PREMISE OF THE STUDY: We isolated 15 polymorphic microsatellite markers from Vachellia farnesiana for use in population genetic studies to determine the native range of the species. • METHODS AND RESULTS: Initially, 454 shotgun sequencing was used to identify and design primers for 68 microsatellite loci. Of these, we trialed 47 loci in the target species, and 42 (89%) amplified a product of expected size. Fifteen of the 47 loci were screened for variation in 21 individuals from the native range of V. farnesiana in southern Mexico and 20 from northwestern Australia. Fourteen loci were polymorphic, with observed heterozygosity ranging from 0.026 to 1.00 (mean = 0.515) and two to 12 alleles per locus (average = 5.2). Cross-amplification was successful in four to 11 loci in three other Vachellia species. • CONCLUSIONS: The new microsatellite loci will be useful in understanding genetic variation and investigating the role of human-mediated dispersal in the current distribution of V. farnesiana.

Isolation and characterization of 13 new polymorphic microsatellite markers in the Phaseolus vulgaris L. (Common Bean) genome.

In this study, 13 polymorphic microsatellite markers were isolated from the Phaseolus vulgaris L. (common bean) by using the Fast Isolation by AFLP of Sequence COntaining Repeats (FIASCO) protocol. These markers revealed two to seven alleles, with an average of 3.64 alleles per locus. The polymorphic information content (PIC) values ranged from 0.055 to 0.721 over 13 loci, with a mean value of 0.492, and 7 loci having PIC greater than 0.5. The expected heterozygosity (H(E)) and observed heterozygosity (H(O)) levels ranged from 0.057 to 0.814 and from 0.026 to 0.531, respectively. Cross-species amplification of the 13 prime pairs was performed in its related specie of Vigna unguiculata L. Seven out of all these markers showed cross-species transferability. These markers will be useful for future genetic diversity and population genetics studies for this agricultural specie and its related species.

A global view of 54,001 single nucleotide polymorphisms (SNPs) on the Illumina BovineSNP50 BeadChip and their transferability to water buffalo.

The Illumina BovineSNP50 BeadChip features 54,001 informative single nucleotide polymorphisms (SNPs) that uniformly span the entire bovine genome. Among them, 52,255 SNPs have locations assigned in the current genome assembly (Btau_4.0), including 19,294 (37%) intragenic SNPs (i.e., located within genes) and 32,961 (63%) intergenic SNPs (i.e., located between genes). While the SNPs represented on the Illumina Bovine50K BeadChip are evenly distributed along each bovine chromosome, there are over 14,000 genes that have no SNPs placed on the current BeadChip. Kernel density estimation, a non-parametric method, was used in the present study to identify SNP-poor and SNP-rich regions on each bovine chromosome. With bandwidth = 0.05 Mb, we observed that most regions have SNP densities within 2 standard deviations of the chromosome SNP density mean. The SNP density on chromosome X was the most dynamic, with more than 30 SNP-rich regions and at least 20 regions with no SNPs. Genotyping ten water buffalo using the Illumina BovineSNP50 BeadChip revealed that 41,870 of the 54,001 SNPs are fully scored on all ten water buffalo, but 6,771 SNPs are partially scored on one to nine animals. Both fully scored and partially/no scored SNPs are clearly clustered with various sizes on each chromosome. However, among 43,687 bovine SNPs that were successfully genotyped on nine and ten water buffalo, only 1,159 were polymorphic in the species. These results indicate that the SNPs sites, but not the polymorphisms, are conserved between two species. Overall, our present study provides a solid foundation to further characterize the SNP evolutionary process, thus improving understanding of within- and between-species biodiversity, phylogenetics and adaption to environmental changes.