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.     

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.     

Evolution and recombination of bovine DNA repeats

The history of the abundant repeat elements in the bovine genome has been studied by comparative hybridization and PCR. The Bov-A and Bov-B SINE elements both emerged just after the divergence of the Camelidae and the true ruminants. A 31-bp subrepeat motif in satellites of the Bovidae species cattle, sheep, and goat is also present in Cervidae (deer) and apparently predates the Bovidae. However, the other components of the bovine satellites were amplified after the divergence of the cattle and the Caprinae (sheep and goat). A 23-bp motif, which as subrepeat of two major satellites occupies 5% of the cattle genome, emerged only after the split of the water buffalo and other cattle species. During the evolution of the Bovidae the satellite repeat units were shaped by recombination events involving subrepeats, other satellite components, and SINE elements. Differences in restriction sites of homologous satellites indicate a continuing rapid horizontal spread of new sequence variants. Correspondence to: J.A. Lenstra     

Detection of Polymorphism and Sequence Characterization of Toll-Like Receptor 7 Gene of Indian Goat Revealing Close Relationship Between Ruminant Species

In this study, approximately 3.4 kb nucleotide sequence of caprine TLR7 (Toll-like receptor 7) gene was generated from twelve different Indian goat breeds belonging to different geographical regions. Goat TLR7 gene ORF (Open Reading Frame) was found to be 3141 nucleotides long coding for 1046 amino acids similar to sheep. The sequence analysis at nucleotide level revealed goat TLR7 having 99.5% homology with sheep, followed by other livestock species. Simple Modular Architecture Research Tool (SMART) was used for the structural analysis of goat TLR7 that showed the presence of 22 leucine rich repeats (LRRs) along with single Toll/interleukin-1 receptor (TIR) domains. TIR domain, when compared, was found to be similar in ruminant species, goat, sheep, cattle, and buffalo. The phylogenetic analysis also revealed grouping of all ruminant species together, goat being closer to sheep followed by cattle and buffalo. A total of 22 polymorphic sites were observed in TLR7 gene of 24 goats representing 12 different breeds, out of which 19 were present within the coding region and three in 3'UTR. Out of the seven nonsynonymous SNPs, two were in ectodomains and one in TIR domain. Overall our results indicate substantial variation within goat TLR7 gene, which could be exploited for association with disease susceptibility.     

Adaptive selection signatures in river buffalo with emphasis on immune and major histocompatibility complex genes

River buffalo is an agriculturally important species with many traits, such as disease tolerance, which promote its use worldwide. Highly contiguous genome assemblies of the river buffalo, goat, pig, human and two cattle subspecies were aligned to study gene gains and losses and signs of positive selection. The gene families that have changed significantly in river buffalo since divergence from cattle play important roles in protein degradation, the olfactory receptor system, detoxification and the immune system. We used the branch site model in PAML to analyse single-copy orthologs to identify positively selected genes that may be involved in skin differentiation, mammary development and bone formation in the river buffalo branch. The high contiguity of the genomes enabled evaluation of differences among species in the major histocompatibility complex. We identified a Babesia-like L1 LINE insertion in the DRB1-like gene in the river buffalo and discuss the implication of this finding.     

Analysis of conserved microsatellite sequences suggests closer relationship between water buffalo Bubalus bubalis and sheep Ovis aries.

The distribution and evolutionary pattern of the conserved microsatellite repeat sequences (CA)n, (TGG)6, and (GGAT)4 were studied to determine the divergence time and phylogenetic position of the water buffalo, Bubalus bubalis. The mean allelic frequencies of these repeat loci showed a high level of heterozygosity among the euartiodactyls (buffalo, cattle, sheep, and goat). Genetic distances calculated from the allelic frequencies of these microsatellites were used to position Bubalus bubalis in the phylogenetic tree. The tree topology revealed a closer proximity of the Bubalus bubalis to the Ovis aries (sheep) genome than to other domestic species. The estimated time of divergence of the water buffalo genome relative to cattle, goat, sheep, pig, rabbit, and horse was found to be 21, 0.5, 0.7, 94, 20.3, and 408 million years (Myr), respectively. Although water buffaloes share morphological and biochemical similarities with cattle, our study using the microsatellite sequences places the bubaline species in an entirely new phylogenetic position. Our results also suggest that with respect to these repeat loci, the water buffalo genome shares a common ancestry with sheep and goat after the divergence of subfamily Bovinae (Bos taurus) from the family Bovidae.     

Pstl repeat: a family of short interspersed nucleotide element (SINE)-like sequences in the genomes of cattle, goat, and buffalo.

The PstI family of elements are short, highly repetitive DNA sequences interspersed throughout the genome of the Bovidae. We have cloned and sequenced some members of the PstI family from cattle, goat, and buffalo. These elements are approximately 500 bp, have a copy number of 2 x 10(5) - 4 x 10(5), and comprise about 4% of the haploid genome. Studies of nucleotide sequence homology indicate that the buffalo and goat PstI repeats (type II) are similar types of short interspersed nucleotide element (SINE) sequences, but the cattle PstI repeat (type I) is considerably more divergent. Additionally, the goat PstI sequence showed significant sequence homology with bovine serine tRNA, and is therefore likely derived from serine tRNA. Interestingly, Southern hybridization suggests that both types of SINEs (I and II) are present in all the species of Bovidae. Dendrogram analysis indicates that cattle PstI SINE is similar to bovine Alu-like SINEs. Goat and buffalo SINEs formed a separate cluster, suggesting that these two types of SINEs evolved separately in the genome of the Bovidae.     

A catalogue of validated single nucleotide polymorphisms in bovine orthologs of mammalian imprinted genes and associations with beef production traits

Genetic (or ‘genomic’) imprinting, a feature of approximately 100 mammalian genes, results in monoallelic expression from one of the two parentally inherited chromosomes. To date, most studies have been directed on imprinted genes in murine or human models; however, there is burgeoning interest in the effects of imprinted genes in domestic livestock species. In particular, attention has focused on imprinted genes that influence foetal growth and development and that are associated with several economically important production traits in cattle, sheep and pigs. We have re-sequenced regions in 20 candidate bovine imprinted genes in order to validate single nucleotide polymorphisms (SNPs) that may influence important production traits in cattle. Putative SNPs detected via re-sequencing were subsequently re-formatted for high-throughput SNP genotyping in 185 cattle samples comprising 138 performance-tested European Bos taurus (all Limousin bulls), 29 African B. taurus and 18 Indian B. indicus samples. Analysis of the resulting genotypic data identified 117 validated SNPs. Preliminary genotype–phenotype association analyses using 83 SNPs that were polymorphic in the Limousin samples with minor allele frequencies ⩾0.05 revealed significant associations between two candidate bovine imprinted genes and a range of important beef production traits: average daily gain, average feed intake, live weight, feed conversion ratio, residual feed intake and residual gain. These genes were the Ras protein-specific guanine nucleotide releasing factor gene (RASGRF1) and the zinc finger, imprinted 2 gene (ZIM2). Despite the relatively small sample size used in these analyses, the observed associations with production traits are supported by the purported biological function of the RASGRF1 and ZIM2 gene products. These results support the hypothesis that imprinted genes contribute significantly to important complex production traits in cattle. Furthermore, these SNPs may be usefully incorporated into future marker-assisted and genomic selection breeding schemes.     

Genomewide identification and annotation of SNPs in Bubalus bubalis

The present study was carried out to identify and annotate the genome wide SNPs in Murrah buffalo genome. A total of 21.2 million raw reads from 4 pooled female Murrah buffalo samples were obtained using restriction enzyme digestion followed by sequencing with Illumina Hiseq 2000. After quality filtration, the reads were aligned to Murrah buffalo genome (ICAR-NBAGR) and Water buffalo genome (UMD_CASPUR_WB_2.0) which resulted in 99.37% and 99.67% of the reads aligning, respectively. A total of 130,688 high quality SNPs along with 35,110 indels were identified versus the Murrah bufffalo genome. Similarly 219,856 high quality SNPs along with 15,201 indels were identified versus the Water buffalo genome. We report 483 SNPs in 66 genes affecting Milk Production, 436 SNPs in 38 genes affecting fertility and 559 SNPs in 72 genes affecting other major traits. The average genome coverage was 13.4% and 14.8% versus the Murrah and Water buffalo genomes, respectively.     

Genome-Wide Association Studies Using Haplotypes and Individual SNPs in Simmental Cattle

Recent advances in high-throughput genotyping technologies have provided the opportunity to map genes using associations between complex traits and markers. Genome-wide association studies (GWAS) based on either a single marker or haplotype have identified genetic variants and underlying genetic mechanisms of quantitative traits. Prompted by the achievements of studies examining economic traits in cattle and to verify the consistency of these two methods using real data, the current study was conducted to construct the haplotype structure in the bovine genome and to detect relevant genes genuinely affecting a carcass trait and a meat quality trait. Using the Illumina BovineHD BeadChip, 942 young bulls with genotyping data were introduced as a reference population to identify the genes in the beef cattle genome significantly associated with foreshank weight and triglyceride levels. In total, 92,553 haplotype blocks were detected in the genome. The regions of high linkage disequilibrium extended up to approximately 200 kb, and the size of haplotype blocks ranged from 22 bp to 199,266 bp. Additionally, the individual SNP analysis and the haplotype-based analysis detected similar regions and common SNPs for these two representative traits. A total of 12 and 7 SNPs in the bovine genome were significantly associated with foreshank weight and triglyceride levels, respectively. By comparison, 4 and 5 haplotype blocks containing the majority of significant SNPs were strongly associated with foreshank weight and triglyceride levels, respectively. In addition, 36 SNPs with high linkage disequilibrium were detected in the GNAQ gene, a potential hotspot that may play a crucial role for regulating carcass trait components.     

Novel SNP Discovery in African Buffalo,Syncerus caffer,Using High-Throughput Sequencing

The African buffalo, Syncerus caffer, is one of the most abundant and ecologically important species of megafauna in the savannah ecosystem. It is an important prey species, as well as a host for a vast array of nematodes, pathogens and infectious diseases, such as bovine tuberculosis and corridor disease. Large-scale SNP discovery in this species would greatly facilitate further research into the area of host genetics and disease susceptibility, as well as provide a wealth of sequence information for other conservation and genomics studies. We sequenced pools of Cape buffalo DNA from a total of 9 animals, on an ABI SOLiD4 sequencer. The resulting short reads were mapped to the UMD3.1 Bos taurus genome assembly using both BWA and Bowtie software packages. A mean depth of 2.7× coverage over the mapped regions was obtained. Btau4 gene annotation was added to all SNPs identified within gene regions. Bowtie and BWA identified a maximum of 2,222,665 and 276,847 SNPs within the buffalo respectively, depending on analysis method. A panel of 173 SNPs was validated by fluorescent genotyping in 87 individuals. 27 SNPs failed to amplify, and of the remaining 146 SNPs, 43–54% of the Bowtie SNPs and 57–58% of the BWA SNPs were confirmed as polymorphic. dN/dS ratios found no evidence of positive selection, and although there were genes that appeared to be under negative selection, these were more likely to be slowly evolving house-keeping genes.     

Draft genome of the river water buffalo

Water buffalo (Bubalus bubalis), a large‐sized member of the Bovidae family, is considered as an important livestock species throughout Southeast Asia. In order to better understand the molecular basis of buffalo improvement and breeding, we sequenced and assembled the genome (2n=50) of a river buffalo species Bubalus bubalis from Bangladesh. Its genome size is 2.77 Gb, with a contig N50 of 25 kb and the scaffold N50 of 6.9 Mbp. Based on the assembled genome, we annotated 24,613 genes for future functional genomics studies. Phylogenetic tree analysis of cattle and water buffalo lineages showed that they diverged about 5.8–9.8 million years ago. Our findings provide an insight into the water buffalo genome which will contribute in further research on buffalo such as molecular breeding, understanding complex traits, conservation, and biodiversity.     

Association of HSL gene E1-c.276C>T and E8-c.51C>T mutation with economical traits of Chinese Simmental cattle

Hormone-sensitive lipase (HSL) is responsible for the decomposition of triglycerides in adipose tissue to release free fatty acids, and it is a key rate-limiting enzyme in the regulation of adipose tissue deposition and decomposition. The objective of this study was to evaluate the association between novel SNPs in the coding region of bovine HSL gene and carcass and meat quality traits of Chinese Simmental-cross steers. Two novel SNPs were genotyped and the 47 traits of carcass and meat quality traits were measured in the population studied. Statistical analysis revealed that the SNPs of HSL gene were associated with the carcass and meat quality traits. The individuals with TT genotypes of E1-276C>T showed significant higher dressing percentage, net meat rate, hind legs circumference, fat coverage rate, mesenteric fat and kidney fat (p < 0.05). E8-51C>T (P17S) also showed a significant association with the pH of beef and fatty acids content in Chinese Simmental cattle (p < 0.01). Our findings indicated that polymorphisms in HSL might be one of important genetic factors that influence carcass yield and meat quality in beef cattle, and it may be a useful marker for meat quality traits in future marker-assisted selection programs in beef cattle breeding and production.     

Novel SNPs in the bovine ADIPOQ and PPARGC1A genes are associated with carcass traits in Hanwoo (Korean cattle)

Adiponectin (ADIPOQ) modulates several biological processes including energy homeostasis, glucose and lipid metabolism. The bovine ADIPOQ gene was located near the QTL affecting marbling, ribeye muscle area and fat thickness on BTA1. The gene encoding peroxisome proliferator-activated receptor-γ coactivator-1α (PPARGC1A) was located within the QTL region of the traits on BTA6. Moreover, its protein product has various biological functions such as cellular energy homeostasis, including adaptive thermogenesis, adipogenesis and gluconeogenesis. Therefore, the ADIPOQ and PPARGC1A genes are a positional and functional candidate gene for carcass traits in beef cattle. The objectives of this study were to identify polymorphisms in the bovine ADIPOQ and PPARGC1A genes, to evaluate their associations with carcass traits in Hanwoo (Korean cattle) population. We identified nine SNPs in the ADIPOQ gene. Two SNPs (DQ156119: g.1436T > C and DQ156119: g.1454A > G) in the promoter region were recognized as new SNPs identified in Hanwoo. Association analysis indicated that the g.1454A > G SNP genotype was significantly associated with effects on LMA (P = 0.004) and BF (P = 0.021). The ADIPOQ haplotype was also found to have significant effect on the LMA. In the PPARGC1A gene, we identified 11 SNPs in the two unexplored regions (intron 3 and 5). Among them, seven SNPs were located in intron 3 and four SNPs were located in intron 5. Of these 11 putative novel SNPs, two SNPs (AY839822: g.292C > T and AY839823: g.1064C > T) with minor allele frequency (MAF) > 0.20 were examined for associations with carcass traits. The association analysis revealed that both SNPs in PPARGC1A gene were significantly associated with LMA (P < 0.05). These findings suggest that the SNPs of bovine ADIPOQ and PPARGC1A genes may be a useful molecular marker for selection of carcass traits in Hanwoo.     

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.     

Intergeneric embryo transfer between water buffalo and domestic cattle

The feasibility of reciprocal embryo transfer between water buffalo (Bubalus bubalis) and domestic cattle (Bos taurus) was studied. Eight mature water buffalo females were superovulated and bred naturally to a male water buffalo. Sixteen water buffalo embryos were recovered nonsurgically using routine bovine embryo transfer procedures. No pregnancies resulted after transfer of thirteen water buffalo embryos to synchronized Holstein heifers. Physiological age of the water buffalo embryos was ahead of chronological age when compared to bovine embryo development.Two cattle embryos were transferred nonsurgically to two synchronized water buffalo recipients. One recipient was diagnosed pregnant; however, she subsequently aborted between 2.5 and 3.0 months of gestation.The significance of successful intergeneric embryo transfer between cattle and water buffaloes would be the provision of a model for the preservation of endangered genera and species for which common recipients are readily available.The reproductive anatomy as well as the reproductive physiology of the estrous cycle of the donor and the recipient species must be reasonably similar. Minor endocrinological incompatibilities might be circumvented by transfer of alien embryos to mated recipients. Problems of immunological incompatibilities might be minimized by modifying the immune response of the recipient and/or the antigenicity of the embryo by microsurgical manipulation.