A principal component regression based genome wide analysis approach reveals the presence of a novel QTL on BTA7 for MAP resistance in holstein cattle

Genome-wide analysis reveals the ancient and recent admixture history of East African Shorthorn Zebu from Western Kenya

The Kenyan East African zebu cattle are valuable and widely used genetic resources. Previous studies using microsatellite loci revealed the complex history of these populations with the presence of taurine and zebu genetic backgrounds. Here, we estimate at genome-wide level the genetic composition and population structure of the East African Shorthorn Zebu (EASZ) of western Kenya. A total of 548 EASZ from 20 sub-locations were genotyped using the Illumina BovineSNP50 v. 1 beadchip. STRUCTURE analysis reveals admixture with Asian zebu, African and European taurine cattle. The EASZ were separated into three categories: substantial (⩾12.5%), moderate (1.56%<X<12.5%) and non-introgressed (⩽1.56%) according to the European taurine genetic proportion. The non-European taurine introgressed animals (n=425) show an unfluctuating zebu and taurine ancestry of 0.84±0.009 s.d. and 0.16±0.009 s.d., respectively, with significant differences in African taurine (AT) and Asian zebu backgrounds across chromosomes (P<0.0001). In contrast, no such differences are observed for the European taurine ancestry (P=0.1357). Excluding European introgressed animals, low and nonsignificant genetic differentiation and isolation by distance are observed among sub-locations (F_st=0.0033, P=0.09; r=0.155, P=0.07). Following a short population expansion, a major reduction in effective population size (N_e) is observed from approximately 240 years ago to present time. Our results support ancient zebu × AT admixture in the EASZ population, subsequently shaped by selection and/or genetic drift, followed by a more recent exotic European cattle introgression.

Impact of host genetics on susceptibility and resistance to Mycobacterium avium subspecies Paratuberculosis infection in domestic ruminants

Johne's disease or Paratuberculosis has emerged as major infectious disease of animals in general and domestic livestock in particular on global basis. There have been major initiatives in developed countries for the control of this incurable malady of animals and human beings alike (inflammatory bowel disease or Crohn's disease). Disease has not received similar attention due to inherent complexities of disease, diagnosis and control, in resource poor counties around the world. However, the rich genetic diverstiy of the otherwise low productive animal population offers opportunity for the control of Johne's disease and improve per animal productivity. Present review aims to gather and compile information available on genetics or resistance to Johne's disease and its future exploitation by resource poor countries rich in animal diversity. This review will also help to create awareness and share knowledge and experience on prevalence and opportunities for control of Johne's disease in the livestock population to boost per animal productivity among developing and poor countries of the world. Breeding of animals for disease resistance provides good, safe, effective and cheaper way of controlling Johne's disease in animals, with especial reference to domestic livestock of developing and poor countries. Study will help to establish better understanding of the correlation between host cell factors and resistance to MAP infection which may have ultimately help in the control of Johne's disease in future.

Genome-wide association analysis and genomic prediction of Mycobacterium avium subspecies paratuberculosis infection in US Jersey cattle

Paratuberculosis (Johne’s disease), an enteric disorder in ruminants caused by Mycobacterium avium subspecies paratuberculosis (MAP), causes economic losses in excess of $200 million annually to the US dairy industry. To identify genomic regions underlying susceptibility to MAP infection in Jersey cattle, a case-control genome-wide association study (GWAS) was performed. Blood and fecal samples were collected from ∼5,000 mature cows in 30 commercial Jersey herds from across the US. Discovery data consisted of 450 cases and 439 controls genotyped with the Illumina BovineSNP50 BeadChip. Cases were animals with positive ELISA and fecal culture (FC) results. Controls were animals negative to both ELISA and FC tests that matched cases on birth date and herd. Validation data consisted of 180 animals including 90 cases (positive to FC) and 90 controls (negative to ELISA and FC), selected from discovery herds and genotyped by Illumina BovineLD BeadChip (∼7K SNPs). Two analytical approaches were used: single-marker GWAS using the GRAMMAR-GC method and Bayesian variable selection (Bayes C) using GenSel software. GRAMMAR-GC identified one SNP on BTA7 at 68 megabases (Mb) surpassing a significance threshold of 5×10⁻⁵. ARS-BFGL-NGS-11887 on BTA23 (27.7 Mb) accounted for the highest percentage of genetic variance (3.3%) in the Bayes C analysis. SNPs identified in common by GRAMMAR-GC and Bayes C in both discovery and combined data were mapped to BTA23 (27, 29 and 44 Mb), 3 (100, 101, 106 and 107 Mb) and 17 (57 Mb). Correspondence between results of GRAMMAR-GC and Bayes C was high (70–80% of most significant SNPs in common). These SNPs could potentially be associated with causal variants underlying susceptibility to MAP infection in Jersey cattle. Predictive performance of the model developed by Bayes C for prediction of infection status of animals in validation set was low (55% probability of correct ranking of paired case and control samples).

Single nucleotide polymorphisms for feed efficiency and performance in crossbred beef cattle

Background

This study was conducted to: (1) identify new SNPs for residual feed intake (RFI) and performance traits within candidate genes identified in a genome wide association study (GWAS); (2) estimate the proportion of variation in RFI explained by the detected SNPs; (3) estimate the effects of detected SNPs on carcass traits to avoid undesirable correlated effects on these economically important traits when selecting for feed efficiency; and (4) map the genes to biological mechanisms and pathways. A total number of 339 SNPs corresponding to 180 genes were tested for association with phenotypes using a single locus regression (SLRM) and genotypic model on 726 and 990 crossbred animals for feed efficiency and carcass traits, respectively.

Results

Strong evidence of associations for RFI were located on chromosomes 8, 15, 16, 18, 19, 21, and 28. The strongest association with RFI (P = 0.0017) was found with a newly discovered SNP located on BTA 8 within the ELP3 gene. SNPs rs41820824 and rs41821600 on BTA 16 within the gene HMCN1 were strongly associated with RFI (P = 0.0064 and P = 0.0033, respectively). A SNP located on BTA 18 within the ZNF423 gene provided strong evidence for association with RFI (P = 0.0028). Genomic estimated breeding values (GEBV) from 98 significant SNPs were moderately correlated (0.47) to the estimated breeding values (EBVs) from a mixed animal model. The significant (P < 0.05) SNPs (98) explained 26% of the genetic variance for RFI. In silico functional analysis for the genes suggested 35 and 39 biological processes and pathways, respectively for feed efficiency traits.

Conclusions

This study identified several positional and functional candidate genes involved in important biological mechanisms associated with feed efficiency and performance. Significant SNPs should be validated in other populations to establish their potential utilization in genetic improvement programs.

Bovine CLEC7A genetic variants and their association with seropositivity in Johne's disease ELISA

Mycobacterium avium ssp. paratuberculosis (MAP) infection in cattle causes significant economic losses to the dairy and beef industries resulting from reduced productivity, premature culling and mortality. Bovine Dectin-1, an important pattern recognition molecule that is able to generate a proinflammatory response by acting alongside Toll like receptor (TLR) signaling, is known to co-operate with TLR2 to specifically activate a macrophage proinflammatory response against mycobacterial infections. Therefore, the goal of this study was to identify single nucleotide polymorphisms (SNPs) in the gene encoding bovine Dectin-1 (CLEC7A) and to assess their association with susceptibility to MAP infection in dairy cattle. Blood and milk samples, collected from commercial dairy operations, were tested for MAP infection using blood and milk ELISAs and a resource population consisting of 197 infected and 242 healthy cattle was constructed. Pooled DNA was used for sequencing and eight single nucleotide polymorphisms (SNPs) were identified. Identified SNPs were genotyped on the resource population using the iPLEX MassARRAY system and statistical analysis was performed using logistic regression fitting the additive and dominance effects of each SNP in the model. Out of a total of eight identified SNPs, five were successfully genotyped, and three out of these five SNPs were found to be in complete linkage. Statistical analysis revealed a strong association between a non-synonymous SNP c.589A>G (p=0.008), and MAP infection status of the resource population inferred by seropositivity in MAP antibody specific ELISAs. This SNP c.589A>G was located in the geneic region that encodes the carbohydrate recognition domain of bovine Dectin-1. Therefore, further investigation of its functional relevance is warranted.

Genome wide association studies for body conformation traits in the Chinese Holstein cattle population

BACKGROUND

Genome-wide association study (GWAS) is a powerful tool for revealing the genetic basis of quantitative traits. However, studies using GWAS for conformation traits of cattle is comparatively less. This study aims to use GWAS to find the candidates genes for body conformation traits.

RESULTS

The Illumina BovineSNP50 BeadChip was used to identify single nucleotide polymorphisms (SNPs) that are associated with body conformation traits. A least absolute shrinkage and selection operator (LASSO) was applied to detect multiple SNPs simultaneously for 29 body conformation traits with 1,314 Chinese Holstein cattle and 52,166 SNPs. Totally, 59 genome-wide significant SNPs associated with 26 conformation traits were detected by genome-wide association analysis; five SNPs were within previously reported QTL regions (Animal Quantitative Trait Loci (QTL) database) and 11 were very close to the reported SNPs. Twenty-two SNPs were located within annotated gene regions, while the remainder were 0.6-826 kb away from known genes. Some of the genes had clear biological functions related to conformation traits. By combining information about the previously reported QTL regions and the biological functions of the genes, we identified DARC, GAS1, MTPN, HTR2A, ZNF521, PDIA6, and TMEM130 as the most promising candidate genes for capacity and body depth, chest width, foot angle, angularity, rear leg side view, teat length, and animal size traits, respectively. We also found four SNPs that affected four pairs of traits, and the genetic correlation between each pair of traits ranged from 0.35 to 0.86, suggesting that these SNPs may have a pleiotropic effect on each pair of traits.

CONCLUSIONS

A total of 59 significant SNPs associated with 26 conformation traits were identified in the Chinese Holstein population. Six promising candidate genes were suggested, and four SNPs showed genetic correlation for four pairs of traits.

Association analysis identifies Melampsora ×columbiana poplar leaf rust resistance SNPs

Populus species are currently being domesticated through intensive time- and resource-dependent programs for utilization in phytoremediation, wood and paper products, and conversion to biofuels. Poplar leaf rust disease can greatly reduce wood volume. Genetic resistance is effective in reducing economic losses but major resistance loci have been race-specific and can be readily defeated by the pathogen. Developing durable disease resistance requires the identification of non-race-specific loci. In the presented study, area under the disease progress curve was calculated from natural infection of Melampsora ×columbiana in three consecutive years. Association analysis was performed using 412 P. trichocarpa clones genotyped with 29,355 SNPs covering 3,543 genes. We found 40 SNPs within 26 unique genes significantly associated (permutated P<0.05) with poplar rust severity. Moreover, two SNPs were repeated in all three years suggesting non-race-specificity and three additional SNPs were differentially expressed in other poplar rust interactions. These five SNPs were found in genes that have orthologs in Arabidopsis with functionality in pathogen induced transcriptome reprogramming, Ca²⁺/calmodulin and salicylic acid signaling, and tolerance to reactive oxygen species. The additive effect of non-R gene functional variants may constitute high levels of durable poplar leaf rust resistance. Therefore, these findings are of significance for speeding the genetic improvement of this long-lived, economically important organism.

Genome-wide association mapping for wood characteristics in Populus identifies an array of candidate single nucleotide polymorphisms

Establishing links between phenotypes and molecular variants is of central importance to accelerate genetic improvement of economically important plant species. Our work represents the first genome-wide association study to the inherently complex and currently poorly understood genetic architecture of industrially relevant wood traits. Here, we employed an Illumina Infinium 34K single nucleotide polymorphism (SNP) genotyping array that generated 29,233 high-quality SNPs in c. 3500 broad-based candidate genes within a population of 334 unrelated Populus trichocarpa individuals to establish genome-wide associations. The analysis revealed 141 significant SNPs (α ≤ 0.05) associated with 16 wood chemistry/ultrastructure traits, individually explaining 3-7% of the phenotypic variance. A large set of associations (41% of all hits) occurred in candidate genes preselected for their suggested a priori involvement with secondary growth. For example, an allelic variant in the FRA8 ortholog explained 21% of the total genetic variance in fiber length, when the trait's heritability estimate was considered. The remaining associations identified SNPs in genes not previously implicated in wood or secondary wall formation. Our findings provide unique insights into wood trait architecture and support efforts for population improvement based on desirable allelic variants.

Genome-wide association study to identify chromosomal regions associated with antibody response to Mycobacterium avium subspecies paratuberculosis in milk of Dutch Holstein-Friesians

Heritability of susceptibility to Johne's disease in cattle has been shown to vary from 0.041 to 0.159. Although the presence of genetic variation involved in susceptibility to Johne's disease has been demonstrated, the understanding of genes contributing to the genetic variance is far from complete. The objective of this study was to contribute to further understanding of genetic variation involved in susceptibility to Johne's disease by identifying associated chromosomal regions using a genome-wide association approach. Log-transformed ELISA test results of 265,290 individual Holstein-Friesian cows from 3,927 herds from the Netherlands were analyzed to obtain sire estimated breeding values for Mycobacterium avium subspecies paratuberculosis (MAP)-specific antibody response in milk using a sire-maternal grandsire model with fixed effects for parity, year of birth, lactation stage, and herd; a covariate for milk yield on test day; and random effects for sire, maternal grandsire, and error. For 192 sires with estimated breeding values with a minimum reliability of 70%, single nucleotide polymorphism (SNP) typing was conducted by a multiple SNP analysis with a random polygenic effect fitting 37,869 SNP simultaneously. Five SNP associated with MAP-specific antibody response in milk were identified distributed over 4 chromosomal regions (chromosome 4, 15, 18, and 28). Thirteen putative SNP associated with MAP-specific antibody response in milk were identified distributed over 10 chromosomes (chromosome 4, 14, 16, 18, 19, 20, 21, 26, 27, and 29). This knowledge contributes to the current understanding of genetic variation involved in Johne's disease susceptibility and facilitates control of Johne's disease and improvement of health status by breeding.

[Advances on genome-wide association study for economically important traits in dairy cattle]

Milk performance traits have been considered as the most important production traits. With the development of balance breeding sense, much attention has been paid to the functional traits such as reproduction, type, health, and longevity traits, which is put into breeding programs. Identification of major genes or genetic markers for milk production and functional traits and their applications in breeding programs of dairy cattle are expected to improve the genetic progress of these traits. Occurrence of millions of SNPs in whole genome and high throughout genotyping techniques has made genome-wide association study to be an important strategy to identify genes responsible for economic traits in domestic animals. This paper reviewed the advances on genome-wide association study for milk production traits and functional traits in dairy cattle.

Genes controlling vaccine responses and disease resistance to respiratory viral pathogens in cattle

Farm animals remain at risk of endemic, exotic and newly emerging viruses. Vaccination is often promoted as the best possible solution, and yet for many pathogens, either there are no appropriate vaccines or those that are available are far from ideal. A complementary approach to disease control may be to identify genes and chromosomal regions that underlie genetic variation in disease resistance and response to vaccination. However, identification of the causal polymorphisms is not straightforward as it generally requires large numbers of animals with linked phenotypes and genotypes. Investigation of genes underlying complex traits such as resistance or response to viral pathogens requires several genetic approaches including candidate genes deduced from knowledge about the cellular pathways leading to protection or pathology, or unbiased whole genome scans using markers spread across the genome. Evidence for host genetic variation exists for a number of viral diseases in cattle including bovine respiratory disease and anecdotally, foot and mouth disease virus (FMDV). We immunised and vaccinated a cattle cross herd with a 40-mer peptide derived from FMDV and a vaccine against bovine respiratory syncytial virus (BRSV). Genetic variation has been quantified. A candidate gene approach has grouped high and low antibody and T cell responders by common motifs in the peptide binding pockets of the bovine major histocompatibility complex (BoLA) DRB3 gene. This suggests that vaccines with a minimal number of epitopes that are recognised by most cattle could be designed. Whole genome scans using microsatellite and single nucleotide polymorphism (SNP) markers has revealed many novel quantitative trait loci (QTL) and SNP markers controlling both humoral and cell-mediated immunity, some of which are in genes of known immunological relevance including the toll-like receptors (TLRs). The sequencing, assembly and annotation of livestock genomes and is continuing apace. In addition, provision of high-density SNP chips should make it possible to link phenotypes with genotypes in field populations without the need for structured populations or pedigree information. This will hopefully enable fine mapping of QTL and ultimate identification of the causal gene(s). The research could lead to selection of animals that are more resistant to disease and new ways to improve vaccine efficacy.

Proteomic analysis of plasma from Holstein cows testing positive for Mycobacterium avium subsp. paratuberculosis (MAP)

Johne's disease (JD) is a widespread and economically important chronic inflammatory disease of the small intestine of ruminants caused by Mycobacterium avium subsp. paratuberculosis (MAP). Although there are several techniques available for diagnosis of JD, their sensitivity is questionable. New proteome profiling methods, such as serum/plasma protein fingerprinting by 2-Dimensional Fluorescence Difference Gel Electrophoresis (2D-DIGE), may therefore be useful for identifying novel protein biomarkers of MAP infection. In this study, plasma samples were collected from 380 Holstein cows and screened for the presence of MAP infection using the M.pt. Johne's antibody Kit (IDEXX). Five negative (MAP-), and 5 strongly positive (MAP+) cows were selected for proteomic analysis. Highly abundant proteins were depleted from the plasma samples using the ProteoMiner technology (Bio-Rad) to enhance the resolution of low abundance proteins. Plasma samples from MAP-, MAP+, and a pooled internal control were labelled with different fluorescent dyes and separated based on their isoelectrical point (IP) and then their molecular weight. Gel images of the fluorescent plasma protein maps were acquired using a Typhoon scanner and analyzed using the DeCyder software. Proteins that were differentially expressed were excised from the gels, trypsin digested, and subjected to MS/MS analysis for identification. Six proteins were identified as being up-regulated at least 2-fold in MAP+ cows including: transferrin, gelsolin isoforms α & β (actin binding protein - ABP), complement subcomponent C1r, complement component C3, amine oxidase - copper containing 3 (AOC3), and coagulation factor II (thrombin) (p<0.05). Two proteins that were down-regulated approximately 2-fold in the MAP+ cows included coagulation factor XIII -B polypeptide (COAFXIII), and fibrinogen γ chain (FGG) and its precursor.

Meta-analysis of two genome-wide association studies of bovine paratuberculosis

Background

Bovine paratuberculosis (ParaTB) also known as Johne's disease, is a contagious fatal disease resulting from infection by Mycobacterium avium subspecies paratuberculosis (MAP). Previous studies have identified loci associated with ParaTB using different measurements to define cases and controls. The objective of this study was to combine the data from two recent studies to identify genetic loci associated with MAP tissue infection and humoral immune response, defined by MAP ELISA-positive cattle, by comparing cases and control animals for one or both measures of infection.

Methodology/Principal Findings

The two populations used for the association analyses were a cohort of MAP tissue infected animals and control Holstein cows from the USA and the second cohort composed of ELISA-positive and ELISA-negative Holstein cows from Italy. Altogether 1190 cattle were genotyped with the Illumina BovineSNP50 BeadChip. SNP markers were removed if the minor allele frequency <0.01 or genotyping failure was >5%. Animals were removed with >5% genotyping failure. Whole genome association analyses were conducted with the GRAMMAR-CG method using two different definitions of control populations.

Conclusion/Significance

The analyses identified several loci (P<5 e-05) associated with ParaTB, defined by positive ELISA and presence of bacteria in tissue compared to ELISA and tissue negative animals, on chromosomes 1, 12 and 15 and one unassigned SNP. These results confirmed associations on chromosome 12 and the unassigned SNP with ParaTB which had been found in the Italian population alone. Furthermore, several additional genomic regions were found associated with ParaTB when ELISA and tissue positive animals were compared with tissue negative samples. These loci were on chromosomes 1, 6, 7, 13, 16, 21,23 and 25 (P<5 e-05). The results clearly indicate the importance of the phenotype definition when seeking to identify markers associated with different disease responses.

Genetic susceptibility to paratuberculosis

Multiple studies indicate that host animal genetics play a role in susceptibility to Mycobacterium avium subsp. paratuberculosis (MAP) infection. However, due to differences in methods used to define MAP-infected animals and controls and differences in methods of genetic analysis, there is as yet no clear consensus on the genes or markers to reliably define the MAP infection susceptibility of any animal species. Meta-analysis of combined studies and larger studies will help resolve the situation in the coming years.

Genetic variation of toll-like receptor genes and infection by Mycobacterium avium ssp. paratuberculosis in Holstein-Friesian cattle

Toll-like receptors (TLR) are membrane proteins that play a key role in innate immunity, by recognizing pathogens and subsequently activating appropriate responses. Mutations in TLR genes are associated with susceptibility to inflammatory and infectious diseases in humans. In cattle, 3 members of the TLR family, TLR1, TLR2, and TLR4, are associated with Mycobacterium avium ssp. paratuberculosis infection, although the extent of this association for the TLR1 and TLR4 receptors has not yet been determined. Moreover, the causal variant in the TLR2 gene has not yet been unequivocally established. In this study, 24 single nucleotide polymorphisms (SNP) in the bovine TLR1, TLR2, and TLR4 genes were selected from the literature, databases, and in silico searches, for a population-based genetic association study of a Spanish Holstein-Friesian sample. Whereas previous results regarding the TLR1 gene were not corroborated, a risk haplotype was detected in TLR2; however, its low frequency indicates that this detected association should be interpreted with caution. In the case of the TLR4 gene, 3 tightly linked SNP were found to be associated with susceptibility to M. avium ssp. paratuberculosis infection. Moreover, one of these SNP, the SNP c.-226G>C, which is localized in the 5'UTR region of the TLR4 gene, has been reported to be able to alter TLR4 expression, raising the possibility that this mutation may contribute to the response of the individual to infection.

Association analyses of the MAS-QTL data set using grammar, principal components and Bayesian network methodologies

BACKGROUND

It has been shown that if genetic relationships among individuals are not taken into account for genome wide association studies, this may lead to false positives. To address this problem, we used Genome-wide Rapid Association using Mixed Model and Regression and principal component stratification analyses. To account for linkage disequilibrium among the significant markers, principal components loadings obtained from top markers can be included as covariates. Estimation of Bayesian networks may also be useful to investigate linkage disequilibrium among SNPs and their relation with environmental variables.For the quantitative trait we first estimated residuals while taking polygenic effects into account. We then used a single SNP approach to detect the most significant SNPs based on the residuals and applied principal component regression to take linkage disequilibrium among these SNPs into account. For the categorical trait we used principal component stratification methodology to account for background effects. For correction of linkage disequilibrium we used principal component logit regression. Bayesian networks were estimated to investigate relationship among SNPs.

RESULTS

Using the Genome-wide Rapid Association using Mixed Model and Regression and principal component stratification approach we detected around 100 significant SNPs for the quantitative trait (p<0.05 with 1000 permutations) and 109 significant (p<0.0006 with local FDR correction) SNPs for the categorical trait. With additional principal component regression we reduced the list to 16 and 50 SNPs for the quantitative and categorical trait, respectively.

CONCLUSIONS

GRAMMAR could efficiently incorporate the information regarding random genetic effects. Principal component stratification should be cautiously used with stringent multiple hypothesis testing correction to correct for ancestral stratification and association analyses for binary traits when there are systematic genetic effects such as half sib family structures. Bayesian networks are useful to investigate relationships among SNPs and environmental variables.

Candidate gene and genome-wide association studies of Mycobacterium avium subsp. paratuberculosis infection in cattle and sheep: a review

Paratuberculosis (Johne's disease), caused by Mycobacterium avium subspecies paratuberculosis, is responsible for significant economic losses in livestock industries worldwide. This organism is also of public health concern due to an unconfirmed link to Crohn's disease. Susceptibility to paratuberculosis has been suggested to have a genetic component. In livestock, a number of candidate genes have been studied, selected on their association to susceptibility in other mycobacterial diseases, their known role in disease pathogenesis or links to susceptibility of humans to Crohn's disease. These genes include solute carrier family 11 member 1 (SLC11A1, formerly NRAMP1), toll-like receptors, caspase associated recruitment domain 15 (CARD15, formerly NOD2), major histocompatibility complex (MHC) and cytokines (interleukin-10 and interferon-gamma) and their receptors. Genome wide association studies have attempted to confirm associations found and identify new genes involved in pathogenesis and susceptibility. There are a number of limitations and difficulties in these approaches, some peculiar to paratuberculosis but others generally applicable to identification of genetic associations for complex traits. The technical approaches and available information for paratuberculosis have expanded rapidly, particularly relating to sheep and cattle. Here we review the current published evidence for a genetic association with paratuberculosis susceptibility, technological advances that have progressed the field and potential avenues for future research.

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.