Investigation of X- and Y-specific single nucleotide polymorphisms in taurine (Bos taurus) and indicine (Bos indicus) cattle

Initially, domesticated African cattle were of taurine type. Today, we find both African Bos taurus and Bos indicus cattle, as well as their crossbreeds, on the continent of Africa and they all share the same set of African taurine mitochondrial DNA haplogroups. In this study, we report genetic variation as substitutions and insertions/deletions (indels) on both the X and Y chromosomes, and use the variation to assess hybridization between taurine and indicine cattle. Six African cattle breeds (four Sanga breeds, including Raya Azebu, Danakil, Caprivi, Nguni; and two Zebu breeds, including Kilimanjaro Zebu and South Kavirondo Zebu) were screened for six new X-chromosomal markers, specifically three single nucleotide polymorphisms and three indels in the DDX3X (previously DBX ) and ZFX genes, and five previously identified Y-chromosomal markers in the DDX3Y (previously DBY ) and ZFY genes. In total, 90 (57 bulls and 33 cows) samples from the African breeds were analysed. We identify five diagnostic haplotypes of indicine and taurine origins on both the X and Y chromosomes. For each breed, the level of indicine introgression varies; in addition to pure taurine, indicine and hybrid X-chromosome individuals, recombinant X-chromosome variants were also detected. These markers are useful molecular tools for assessing the level of indicine admixture in African cattle breeds.     

Association between Sscp Haplotypes at the Bovine Growth Hormone Gene and Milk Protein Percentage

The bovine Growth Hormone gene (bGH) is an attractive candidate gene for milk production in cattle. Single-strand conformation polymorphisms at bGH were identified and used to define haplotype configurations at this gene in the Israeli Holstein dairy cattle population (Bos taurus) and in the parent animals of the International Bovine Reference Family Panel (a collection of B. taurus and B. indicus crosses). B. taurus and B. indicus haplotypes at the bGH gene differed qualitatively, confirming the previously proposed long evolutionary separation of these cattle subraces. Only a small number of bGH haplotypes were present in the Israel Holstein population. One of the haplotypes, apparently of B. indicus origin, was found to have a highly significant positive effect on milk protein percentage. This illustrates the utility of the haplotype approach for uncovering candidate gene involvement in quantitative genetic variation in agricultural populations. The strong effect of an indicine haplotype in a taurine background raises the possibility that indicine alleles at other candidate genes may comprise a genetic resource for improvement of taurine populations. It is proposed that haplotype analysis may be a useful adjunct to measures of genetic distance for evaluating rare breeds with respect to gene conservation.     

Genome sequence and assembly of Bos indicus

Cattle are divided into 2 groups referred to as taurine and indicine, both of which have been under strong artificial selection due to their importance for human nutrition. A side effect of this domestication includes a loss of genetic diversity within each specialized breed. Recently, the first taurine genome was sequenced and assembled, allowing for a better understanding of this ruminant species. However, genetic information from indicine breeds has been limited. Here, we present the first genome sequence of an indicine breed (Nellore) generated with 52X coverage by SOLiD sequencing platform. As expected, both genomes share high similarity at the nucleotide level for all autosomes and the X chromosome. Regarding the Y chromosome, the homology was considerably lower, most likely due to uncompleted assembly of the taurine Y chromosome. We were also able to cover 97% of the annotated taurine protein-coding genes.     

Geographic distribution and frequency of a taurine Bos taurus and an indicine Bos indicus Y specific allele amongst sub-saharan African cattle breeds

We report for the first time, and for the whole of sub-Saharan Africa, the geographical distribution and the frequency of an indicine and a taurine Y specific allele amongst African cattle breeds. A total of 984 males from 69 indigenous African populations from 22 countries were analysed at the microsatellite locus INRA 124. The taurine allele is probably the oldest one on the continent. However, the taurine and the indicine alleles were present in 291 males (30%), and 693 males (70%), respectively. More particularly, 96% of zebu males (n = 470), 50% of taurine males (n = 263), 29% of sanga males (crossbreed Bos taurus x Bos indicus, n = 263) and 95% of zebu x sanga crossbred males (n = 56) had the indicine allele. The Borgou, a breed classified as zebu x taurine cross showed only the zebu allele (n = 12). The indicine allele dominates today in the Abyssinian region, a large part of the Lake Victoria region and the sahelian belt of West Africa. All the sanga males (n = 64) but only one from the Abyssinian region had the indicine allele. The taurine allele is the commonest only among the sanga breeds of the southern African region and the trypanotolerant taurine breeds of West Africa. In West Africa and in the southern Africa regions, zones of introgression were detected with breeds showing both Y chromosome alleles. Our data also reveal a pattern of male zebu introgression in Mozambique and Zimbabwe, probably originating from the Mozambique coast. The sanga cattle from the Lake Victoria region and the Kuri cattle of Lake Chad, cattle populations surrounded by zebu breeds were, surprisingly, completely devoid of the indicine allele. Human migration, phenotypic preferences by the pastoralists, adaptation to specific habitats and to specific diseases are the main factors explaining the present-day distribution of the alleles in sub-Saharan Africa.     

Genome-wide identification of runs of homozygosity islands and associated genes in local dairy cattle breeds

Runs of homozygosity (ROH) are widely used as predictors of whole-genome inbreeding levels in cattle. They identify regions that have an unfavorable effect on a phenotype when homozygous, but also identify the genes associated with traits of economic interest present in these regions. Here, the distribution of ROH islands and enriched genes within these regions in four dairy cattle breeds were investigated. Cinisara (71), Modicana (72), Reggiana (168) and Italian Holstein (96) individuals were genotyped using the 50K v2 Illumina BeadChip. The genomic regions most commonly associated with ROHs were identified by selecting the top 1% of the single nucleotide polymorphisms (SNPs) most commonly observed in the ROH of each breed. In total, 11 genomic regions were identified in Cinisara and Italian Holstein, and eight in Modicana and Reggiana, indicating an increased ROH frequency level. Generally, ROH islands differed between breeds. The most homozygous region (>45% of individuals with ROH) was found in Modicana on chromosome 6 within a quantitative trail locus affecting milk fat and protein concentrations. We identified between 126 and 347 genes within ROH islands, which are involved in multiple signaling and signal transduction pathways in a wide variety of biological processes. The gene ontology enrichment provided information on possible molecular functions, biological processes and cellular components under selection related to milk production, reproduction, immune response and resistance/susceptibility to infection and diseases. Thus, scanning the genome for ROH could be an alternative strategy to detect genomic regions and genes related to important economic traits.     

Genetic diversity among Canadienne, Brown Swiss, Holstein, and Jersey cattle of Canada based on 15 bovine microsatellite markers.

The genetic diversity among Canadienne, Brown Swiss, Holstein, and Jersey cattle was estimated from relationships determined by genotyping 20 distantly related animals in each breed for 15 microsatellites located on separate chromosomes. The Canadienne, Holstein, and Jersey cattle had an average of six alleles per loci compared with five alleles for Brown Swiss. Furthermore, a number of potentially breed-specific alleles were identified. The allele size variance among breeds was similar, but varied considerably among loci. All of the loci studied were equally heterozygous, as were Brown Swiss, Canadienne, and Holstein cattle (0.68-0.69) whereas Jersey cattle showed lower heterozygosity (0.59). The within-breed estimates of genetic distance were greater than zero and significant. The genetic distance between Canadienne and Holstein (0.156), Brown Swiss (0.243), and Jersey (0.235) was negligible, suggesting close relationship. Concurrently, Brown Swiss and Holstein (0.211) cattle also demonstrated close relationship. In contrast, the Jersey breed was genetically distant from the Brown Swiss and Holstein cattle (0.427 and 0.320, respectively). The characterization of Canadienne cattle, as part of the genetic resource conservation effort currently underway in Canada, underscores the difficulty in scientifically establishing unique breeds. Therefore, the need to consider all relevant morphological characteristics and production performance in combination with available cultural, historical, pedigree, and molecular information becomes relevant when identifying breeds for conservation.     

Dynamics of genomic architecture during composite breed development in cattle

Some livestock breeds face the challenge of reduced genetic variation, increased inbreeding depression owing to genetic drift and selection. Hybridization can reverse these processes and increase levels of productivity and adaptation to various environmental stressors. Samples from American Brangus were used to evaluate the indicine/taurine composition through nine generations (~45 years) after the hybridization process was completed. The purpose was to determine how hybridization alters allelic combinations of a breed over time when genetic factors such as selection and drift are operating. Furthermore, we explored genomic regions with deviations from the expected composition from the progenitor breeds and related these regions to traits under selection. The Brangus composition deviated from the theoretical expectation, defined by the breed association, of 62.5% taurine, showing taurine composition to be 70.4 ± 0.6%. Taurine and indicine proportion were not consistent across chromosomes. Furthermore, these non-uniform areas were found to be associated with traits that were probably under selection such as intermuscular fat and average daily gain. Interestingly, the sex chromosomes were predominantly taurine, which could be due to the composite being formed particularly in the final cross that resulted in progeny designated as purebred Brangus. This work demonstrated the process of new breed formation on a genomic level. It suggests that factors like genetic drift, selection and complementarity shift the genetic architecture into a uniquely different population. These findings are important to better understand how hybridization and crossbreeding systems shape the genetic architecture of composite populations.     

Use of breed-specific single nucleotide polymorphisms to discriminate between Holstein and Jersey dairy cattle breeds

Emphasis on livestock genetic improvement in the past decades has led to commercialization of different breeds of livestock species. Breed validation has become increasingly important to assess the safety and authenticity of livestock products in global and domestic markets. The objective of this study was to evaluate the use of breed-specific single nucleotide polymorphisms (SNPs) in discriminating between Holstein and Jersey dairy cattle breeds. Two separate resource populations were used, including a reference population consisting of 498 Holstein and 83 Jersey bull DNA samples, and a validation population consisting of 260 Holstein and 34 Jersey cow DNA samples. Five Jersey-specific and four Holstein-specific SNPs were identified and genotyped on the reference and validation resource populations. The reference population was used to validate the breed-specific SNPs used in this study and to predict the allocation efficiencies and misclassification probabilities of different combinations of SNPs. Individual animals in the validation population were allocated to either breed based on the presence of breed-specific alleles. It was found that any combination of three breed-specific SNPs had, on average, high breed allocation efficiency of >95% and low misclassification probability of <5%. In conclusion, this study demonstrates a simple, yet effective, method of using breed-specific SNPs to discriminate between Jersey and Holstein cattle breeds.     

The extent of linkage disequilibrium in beef cattle breeds using high-density SNP genotypes

Background

The extent of linkage disequilibrium (LD) between molecular markers impacts genome-wide association studies and implementation of genomic selection. The availability of high-density single nucleotide polymorphism (SNP) genotyping platforms makes it possible to investigate LD at an unprecedented resolution. In this work, we characterised LD decay in breeds of beef cattle of taurine, indicine and composite origins and explored its variation across autosomes and the X chromosome.

Findings

In each breed, LD decayed rapidly and r² was less than 0.2 for marker pairs separated by 50 kb. The LD decay curves clustered into three groups of similar LD decay that distinguished the three main cattle types. At short distances between markers (< 10 kb), taurine breeds showed higher LD (r² = 0.45) than their indicine (r² = 0.25) and composite (r² = 0.32) counterparts. This higher LD in taurine breeds was attributed to a smaller effective population size and a stronger bottleneck during breed formation. Using all SNPs on only the X chromosome, the three cattle types could still be distinguished. However for taurine breeds, the LD decay on the X chromosome was much faster and the background level much lower than for indicine breeds and composite populations. When using only SNPs that were polymorphic in all breeds, the analysis of the X chromosome mimicked that of the autosomes.

Conclusions

The pattern of LD mirrored some aspects of the history of breed populations and showed a sharp decay with increasing physical distance between markers. We conclude that the availability of the HD chip can be used to detect association signals that remained hidden when using lower density genotyping platforms, since LD dropped below 0.2 at distances of 50 kb.     

DNA sequence of SSCP haplotypes at the bovine growth hormone (bGH) gene

Previous studies using SSCP and PCR-RFLP methodologies uncovered nine polymorphic sites within the bGH gene, defining eight intragenic haplotypes falling into two main groups. In the present study we report the DNA sequence of these eight haplotypes. A total of 1494 bp were sequenced uncovering a total of 12 sequence variants. Haplotypes within groups differed among themselves at one or two sites, compared across groups, haplotypes of the two groups differed consistently at six sites, each of which was monomorphic within the respective groups. This comes to 4 differentiating sites per kb, suggesting that the two haplotype groups began to diverge about 400,000 years ago. This corresponds approximately to the estimated time of divergence of the Bos taurus and Bos indicus lineages, raising the possibility, supported by other evidence, that the two haplotype classes represent taurine and indicine haplotypes, respectively. Nucleotide sequence divergence of taurine and indicine genomes of this magnitude has far reaching implications with respect to QTL mapping and marker assisted selection in breeds derived from taurine x indicine crosses.     

Bos indicus introgression into (peri‐)alpine cattle breeds – evidence from the analysis of bovine whey protein variants

The major bovine whey proteins, α‐lactalbumin (α‐LA) and β‐lactoglobulin (β‐LG), exhibit breed‐specific genetic variation. The aim of this study was to identify possible new protein variants and determine the distribution of variants across a variety of 18 taurine and indicine cattle breeds applying a DNA‐based sequencing approach. To this end, the open reading frames of the respective genes (LALBA and LGB) were sequenced in 476 animals. Within the LALBA gene, a previously unknown synonymous and a previously undesignated non‐synonymous nucleotide exchange were identified. Furthermore, two known α‐LA variants (A and B) and four known β‐LG variants (A, B, C and W) were determined. The occurrence of typical indicine variants in some taurine cattle breeds, such as Suisse Eringer, German Hinterwälder and Hungarian Grey Steppe, further supports the hypothesis of ancient Bos indicus introgression into (peri‐)alpine cattle breeds.     

Low levels of taurine introgression in the current Brazilian Nelore and Gir indicine cattle populations

Background

Nelore and Gir are the two most important indicine cattle breeds for production of beef and milk in Brazil. Historical records state that these breeds were introduced in Brazil from the Indian subcontinent, crossed to local taurine cattle in order to quickly increase the population size, and then backcrossed to the original breeds to recover indicine adaptive and productive traits. Previous investigations based on sparse DNA markers detected taurine admixture in these breeds. High-density genome-wide analyses can provide high-resolution information on the genetic composition of current Nelore and Gir populations, estimate more precisely the levels and nature of taurine introgression, and shed light on their history and the strategies that were used to expand these breeds.

Results

We used the high-density Illumina BovineHD BeadChip with more than 777 K single nucleotide polymorphisms (SNPs) that were reduced to 697 115 after quality control filtering to investigate the structure of Nelore and Gir populations and seven other worldwide populations for comparison. Multidimensional scaling and model-based ancestry estimation clearly separated the indicine, European taurine and African taurine ancestries. The average level of taurine introgression in the autosomal genome of Nelore and Gir breeds was less than 1% but was 9% for the Brahman breed. Analyses based on the mitochondrial SNPs present in the Illumina BovineHD BeadChip did not clearly differentiate taurine and indicine haplotype groupings.

Conclusions

The low level of taurine ancestry observed for both Nelore and Gir breeds confirms the historical records of crossbreeding and supports a strong directional selection against taurine haplotypes via backcrossing. Random sampling in production herds across the country and subsequent genotyping would be useful for a more complete view of the admixture levels in the commercial Nelore and Gir populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-015-0109-5) contains supplementary material, which is available to authorized users.     

Paternal origins of Chinese cattle

To determine the genetic diversity and paternal origin of Chinese cattle, 302 males from 16 Chinese native cattle breeds as well as 30 Holstein males and four Burma males as controls were analysed using four Y‐SNPs and two Y‐STRs. In Chinese bulls, the taurine Y1 and Y2 haplogroups and indicine Y3 haplogroup were detected in seven, 172 and 123 individuals respectively, and these frequencies varied among the Chinese cattle breeds examined. Y2 dominates in northern China (91.4%), and Y3 dominates in southern China (90.8%). Central China is an admixture zone, although Y2 predominates overall (72.0%). The geographical distributions of the Y2 and Y3 haplogroup frequencies revealed a pattern of male indicine introgression from south to north China. The three Y haplogroups were further classified into one Y1 haplotype, five Y2 haplotypes and one Y3 haplotype in Chinese native bulls. Due to the interplay between taurine and indicine types, Chinese cattle represent an extensive reservoir of genetic diversity. The Y haplotype distribution of Chinese cattle exhibited a clear geographical structure, which is consistent with mtDNA, historical and geographical information.     

Y-SNPs Haplotype Diversity in Four Chinese Cattle Breeds

To investigate the genetic diversity of Chinese cattle, 96 male samples of 4 Chinese native cattle breeds were investigated using 5 single nucleotide polymorphisms specific to the bovine Y chromosome. Two previously described haplotypes (taurine Y2 and indicine Y3) were detected in 74 and 22 animals, respectively. The haplotype frequencies varied amongst the four native breeds. The taurine Y2 haplotype dominated in the Qinchuan, Dabieshan, and Yunba breeds. However, the indicine Y3 haplotype occurred in high frequency in the Enshi breed. Among the four native breeds, Yunba had the highest haplotype diversity (0.4330 ± 0.0750), followed by Qinchuan (0.2899 ± 0.1028) and Enshi (0.2222 ± 0.1662), Dabieshan was the least differentiated (0.1079 ± 0.0680). Compared with some foreign cattle breeds, the low level of haplotype diversity was detected in our breeds (0.2633 ± 0.1030).     

Y chromosome diversity and paternal origin of Chinese cattle

To determine the Y chromosome genetic diversity and paternal origin of Chinese cattle, 369 bulls from 17 Chinese native cattle breeds and 30 bulls from Holstein and four bulls from Burma were analyzed using a recently discovered USP9Y marker that could distinguish between taurine and indicine cattle more efficiently. In total, the taurine Y1, Y2 haplogroup and indicine Y3 haplogroup were detected in 7 (1.9 %), 193 (52.3 %) and 169 (45.8 %) individuals of 17 Chinese native breeds, respectively, although these frequencies varied amongst the Chinese native cattle breeds examined. Y2 dominates in northern China (91.4 %), while Y3 dominates in southern China (81.2 %). Central China is an admixture zone with Y2 predominating overall (72.0 %). Our results demonstrate that Chinese cattle have two paternal origins, one from B. taurus (Y2) and the other from B. indicus (Y3). The Y1 haplogroup may originate from the imported beef cattle breeds in western countries. The geographical distributions of the Y2 and Y3 haplogroup frequencies reveal a pattern of male indicine introgression from south to north China, and male taurine introgression from north to south China.     

Differences in the rumen methanogen populations of lactating Jersey and Holstein dairy cows under the same diet regimen

In the dairy cattle industry, Holstein and Jersey are the breeds most commonly used for production. They differ in performance by various traits, such as body size, milk production, and milk composition. With increased concerns about the impact of agriculture on climate change, potential differences in other traits, such as methane emission, also need to be characterized further. Since methane is produced in the rumen by methanogenic archaea, we investigated whether the population structure of methanogen communities would differ between Holsteins and Jerseys. Breed-specific rumen methanogen 16S rRNA gene clone libraries were constructed from pooled PCR products obtained from lactating Holstein and Jersey cows, generating 180 and 185 clones, respectively. The combined 365 sequences were assigned to 55 species-level operational taxonomic units (OTUs). Twenty OTUs, representing 85% of the combined library sequences, were common to both breeds, while 23 OTUs (36 sequences) were found only in the Holstein library and 12 OTUs (18 sequences) were found only in the Jersey library, highlighting increased diversity in the Holstein library. Other differences included the observation that sequences with species-like sequence identity to Methanobrevibacter millerae were represented more highly in the Jersey breed, while Methanosphaera-related sequences and novel uncultured methanogen clones were more frequent in the Holstein library. In contrast, OTU sequences with species-level sequence identity to Methanobrevibacter ruminantium were represented similarly in both libraries. Since the sampled animals were from a single herd consisting of two breeds which were fed the same diet and maintained under the same environmental conditions, the differences we observed may be due to differences in host breed genetics.     

Novel single nucleotide polymorphisms and haplotypes within the bovine<Emphasis Type="Italic"> CXCR2</Emphasis> gene

The objectives of this study were to identify single nucleotide polymorphisms (SNPs) and resulting haplotypes in the bovine CXCR2 gene. A 311-bp segment of the bovine CXCR2 gene was amplified and sequenced. Five SNPs at positions 612, 684, 777, 858, and 861 were expressed in both Holstein and Jersey dairy cattle. Four SNPs resulted in synonymous substitutions, while a non-synonymous switch at position 777 (GC) resulted in a glutamine to histidine substitution at amino acid residue 245. Strong linkage disequilibrium was exhibited for both breeds among all five loci (P<0.001). Both allele and genotype frequencies differed significantly with respect to breed at four of the five loci (P<0.001). The five polymorphisms generated ten distinct haplotypes. Six haplotypes were common between the two breeds, while Holsteins and Jerseys each uniquely expressed two haplotypes. Of the six common haplotypes, two represented 83% of the Jersey population; whereas four of these haplotypes represented 95% of the Holstein population.