Extensive MHC class II DRB3 diversity in African and European cattle

Genetic deversity at the highly polymorphic BoLA-DRB3 locus was investigated by DNA sequence analyses of 18 African cattle from two breeds representing the two subspecies of cattle, Bos primigenius indicus and Bos primigenius taurus. Yhe polymorphism was compared with that found in a sample ofd 32 European cattle from four breeds, all classified as B. p. taurus. Particularly extensive genetic diversity was found among African cattle, in which as many as 18 alleles were recognized in this small random sample of animals from two breeds. The observed similarity in allele frequency distribution between the two African populations, N'Dama and Zebu cattle, is consistent with the recent recognition of gene flow between B. p. indicus and B. P taurus cattle in Africa. A total of 30 DRB3 alleles were documented and as many as 26 of these were classified as major allelic types showing at least five amino acid substitutions compared with other major types. The observation of extensive genetic diversity at MHC loci in cattle, as well as in other farm animals, provides a compelling argument against matin-type preferences as a primary cause in maintaining major histocompatibility complex diversity, since the reproduction of these animals has been controlled by humans for many generations.The nucleotide sequence data reported in this paper have been submitted to the EMBL nucleotide sequence database and have been given the accession numbers X87641-X87670     

Methemoglobin reductase in three species of Bovidae

Isoelectric focusing of red cell hemolysates revealed several isozymes that stain for NADH-methemoglobin reductase. Evidence for two different genetic loci controlling the banding patterns was obtained. One locus controlled a single band present in all animals tested. The second locus controlled ten different banding patterns that could be accounted for by four codominant alleles. Band B occurred in Bison bison. Bands A and C occurred in Bos indicus and band D occurred in both Bos indicus and Bos taurus. Bands A, C, and D were not observed in Bison bison and bands A, B, and C were not observed in Bos taurus.     

Diversity analysis at MHC class II DQA locus in buffalo (Bubalus bubalis) indicates extensive duplication and trans-species evolution

Variation at MHC Class II-DQA locus in riverine and swamp buffaloes (Bubu) has been explored in this study. Through sequencing of buffalo DQA, 48 nucleotide variants identified from 17 individuals, reporting 42 novel alleles, including one pseudogene. Individual animal displayed two to seven variants, suggesting the presence of more than two Bubu-DQA loci, as an evidence of extensive duplication. dN values were found to be higher than dS values at peptide binding sites, separately for riverine and swamp buffaloes, indicating locus being under positive selection. Evolutionary analysis revealed numerous trans-species polymorphism with alleles from water buffalo assigned to at least three different loci (Bubu-DQA1, DQA2, DQA3). Alleles of both the sub-species intermixed within the cluster, showing convergent evolution of MHC alleles in bovines. The results thus suggest that both riverine and swamp buffaloes share con-current arrangement of DQA region, comparable to cattle in terms of copy number and population polymorphism.     

Complete mitogenome reveals genetic divergence and phylogenetic relationships among Indian cattle (Bos indicus) breeds

Indigenous cattle of India belong to the species, Bos indicus and they possess various adaptability and production traits. However, little is known about the genetic diversity and origin of these breeds. To investigate the status, we sequenced and analyzed the whole mitochondrial DNA (mtDNA) of seven Indian cattle breeds. In total, 49 single-nucleotide variants (SNVs) were identified among the seven breeds analyzed. We observed a common synonymous SNV in the COII gene (m.7583G?>?A) of all the breeds studied. The phylogenetic analysis and genetic distance estimation showed the close genetic relationship among the Indian cattle breeds, whereas distinct genetic differences were observed between Bos indicus and Bos taurus cattle. Our results indicate a common ancestor for European Zwergzebu breed and South Indian cattle. The estimated divergence time demonstrated that the Bos indicus and Bos taurus cattle lineages diverged 0.92 million years ago. Our study also demonstrates that ancestors of present zebu breeds originated in South and North India separately ～30,000 to 20,000 years ago. In conclusion, the identified genetic variants and results of the phylogenetic analysis may provide baseline information to develop appropriate strategies for management and conservation of Indian cattle breeds.     

The relative influence of temperature and humidity on cutaneous function in Bos indicus and Bos taurus females

The relative importance of dry- and wet-bulb temperatures on cutaneous function inBos indicus andBos taurus females under humid tropical climatic conditions was evaluated. The parameters investigated were sweating rate and skin temperature, while the species utilised were zebu White Fulani (Bos indicus) and German Brown and German Black and White (Bos taurus).The sweating rate, irrespective of breed, differed with the site of sampling and was more influenced by dry-bulb (59%) than by wet-bulb temperature (41%). Skin temperature responded more to wet-bulb temperature in White Fulani and German Black and White cattle, but not in German Brown cattle.It is concluded that the response of the animals, with respect to sweating, was similar but that the efficiency of sweating, judged by the lowering of skin temperature, was higher inBos indicus than inBos taurus. This, in part, may explain the higher degree of comfort demonstrated byBos indicus under tropical conditions.     

Characterization of Bison bison major histocompatibility complex class IIa haplotypes

American bison (Bison bison) and domestic cattle (Bos taurus and Bos indicus) evolved from a common ancestor 1–1.4 million years ago. Nevertheless, they show dramatic differences in their susceptibility to infectious diseases, including malignant catarrhal fever (MCF). Although bison are highly susceptible to ovine herpesvirus-2 (OvHV-2) associated MCF, about 20% of healthy domesticated and wild bison are positive for OvHV-2 antibody. We are interested in testing the hypothesis that, within the bison population, the polymorphism of major histocompatibility complex (MHC) class II genes influences resistance to MCF. However, since little was known about the MHC class II genes of bison, it was necessary to first characterize class II haplotypes present in Bi. bison (Bibi). Thus, the MHC class II haplotypes carried by 14 bison were characterized by the PCR-based cloning and sequencing of their DRB3, DQA, and DQB alleles. Twelve MHC class II haplotypes were identified in the 14 bison. These haplotypes comprised six previously reported and six new Bibi-DRB3 alleles, along with 11 Bibi-DQA and 10 Bibi-DQB alleles. For each bison class II allele, it was possible to identify closely related cattle sequences. The closest bison and bovine DQA, DQB, and DRB3 alleles, on average, differed by only 1.3, 3.5, and 5.8 amino acids, respectively. Furthermore, bison MHC haplotypes with both nonduplicated and duplicated DQ genes were identified; these haplotypes appear to have originated from the same ancestral haplotypes as orthologous cattle haplotypes. This study was supported by USDA-Agricultural Research Service grant CWU-5348-32000-018-00D. While working on this project, Dr. Bharat Bhushan was supported by a fellowship from the World-Bank-sponsored National Agricultural Technology Project of the Indian Council of Agricultural Research, Indian Ministry of Agriculture, New Delhi, India     

Sequence comparison of mitochondrial tRNA genes and origin of light strand replication in Bos taurus and Nellore (Bos indicus) breeds

Although the complete bovine mitochondrial DNA molecule has been previously sequencedand sequence comparisons of the mitochondrial displacement loop have been performed, detailed sequence information is limited on coding regions of mitochondrial DNA within and among breeds of Bos taurus and Bos indicus. This study analysed polymorphism of the mitochondrial DNA transfer RNA genes for trypto-phan, alanine, asparagine, cysteine, tyrosine and the origin of light strand replication among Ayrshire, Canadian, Belgium Blue, Brown Swiss, Hereford, Jersey, Limousine, Piedmon-taise, Red Angus, Simmental (Bos taurus) and a Nellore (Bos indicus). Nucleotide sequence analysis of a 420-bp fragment of mitochondrial DNA comprising the five transfer RNA genes showed 100% homology among single individuals of the Bos taurus breeds. The Nellore breed showed guanine to adenine substitutions in the DHU arm of asparagine tRNA and in the origin of light-strand replication. This equates to a 0.5% sequence difference between the Nellore andBos taurus breeds and may reflect an independent evolutionary origin of the species.     

A polymorphism in randomly amplified DNA that differentiates the Y chromosomes of Bos indicus and Bos taurus

A small number of west African Bos taurus cattle breeds, including the N'Dama, constitute a valuable genetic resource by virtue of their ability to remain productive under trypanosomiasis challenge. However, introgression of Bos indicus genes into the trypanotolerant breeds, particularly by introduction of zebu bulls, is a threat to this resource. This work describes the characterization and cloning of a bovine randomly amplified polymorphic DNA (RAPD) that is generated in polymerase chain reaction (PCR) with the 10 base primer ILO1065 from Bos indicus male templates, but not from B. taurus male templates or female templates of either type. Male-specific sequences with homology to the RAPD also occur in B. taurus breeds. This suggests that the polymorphism may be due to base substitution(s) in an ILO1065 priming site, or insertion/deletion events either affecting priming sites or occurring between sites on the cattle Y chromosome. We have shown that cattle, whether of B. indicus or B. taurus phenotype, which possess a typically B. indicus metaphase Y chromosome on the basis of QFQ banding, have a B. indicus ILO1065-generated genotype. The ILO1065-primed RAPD can be used in a simple dot blot assay as a probe of RAPD-PCR products, to provide a convenient, reliable and effective means of detecting introgression of zebu genes in B. taurus cattle populations.     

Polymorphism of mitochondrial DNA (mtDNA) in cattle and buffaloes

Mitochondrial DNA (mtDNA) from two breeds of cattle, viz., [Hariana (Bos indicus), Holstein (Bos taurus)] and Indian water buffalo (Bubalis bubalus), was analyzed using 13 restriction endonucleases which recognized an average of about 40 six-base sites. Polymorphism among cattle was detected with six of these enzymes. The two Holstein differed at six sites, whereas the Hariana breed (Bos indicus) did not show any site polymorphism. Surprisingly, the Hariana type differed by only one site from one of the Holstein types. The total size of buffalo mtDNA was estimated to be 16.4 kb. Polymorphism within the Murrah buffalo breed was observed with respect to aBglI site. Scarcely any of the restriction fragments of buffalo mtDNA matched those of cattle mtDNA.     

Detection of quantitative trait loci in Bos indicus and Bos taurus cattle using genome-wide association studies

Background

The apparent effect of a single nucleotide polymorphism (SNP) on phenotype depends on the linkage disequilibrium (LD) between the SNP and a quantitative trait locus (QTL). However, the phase of LD between a SNP and a QTL may differ between Bos indicus and Bos taurus because they diverged at least one hundred thousand years ago. Here, we test the hypothesis that the apparent effect of a SNP on a quantitative trait depends on whether the SNP allele is inherited from a Bos taurus or Bos indicus ancestor.

Methods

Phenotype data on one or more traits and SNP genotype data for 10 181 cattle from Bos taurus, Bos indicus and composite breeds were used. All animals had genotypes for 729 068 SNPs (real or imputed). Chromosome segments were classified as originating from B. indicus or B. taurus on the basis of the haplotype of SNP alleles they contained. Consequently, SNP alleles were classified according to their sub-species origin. Three models were used for the association study: (1) conventional GWAS (genome-wide association study), fitting a single SNP effect regardless of subspecies origin, (2) interaction GWAS, fitting an interaction between SNP and subspecies-origin, and (3) best variable GWAS, fitting the most significant combination of SNP and sub-species origin.

Results

Fitting an interaction between SNP and subspecies origin resulted in more significant SNPs (i.e. more power) than a conventional GWAS. Thus, the effect of a SNP depends on the subspecies that the allele originates from. Also, most QTL segregated in only one subspecies, suggesting that many mutations that affect the traits studied occurred after divergence of the subspecies or the mutation became fixed or was lost in one of the subspecies.

Conclusions

The results imply that GWAS and genomic selection could gain power by distinguishing SNP alleles based on their subspecies origin, and that only few QTL segregate in both B. indicus and B. taurus cattle. Thus, the QTL that segregate in current populations likely resulted from mutations that occurred in one of the subspecies and can have both positive and negative effects on the traits. There was no evidence that selection has increased the frequency of alleles that increase body weight.     

Cattle breed alters reproduction inAmblyomma americanum (Acari: Ixodidae)

FemaleAmblyomma americanum [L.] that feed onBos taurus cattle weigh significantly (P=0.05) more (x=505.4 mg) than females that feed onB. indicus x B. taurus cattle (x=450.7 mg) but not females that feed onB. indicus cattle (x=489.7 mg). Females that feed onB. taurus cattle lay more eggs (x=3492) and produce more larvae (x=3243) than females that feed onB. indicus cattle (x=3103 eggs; 2925 larvae) orB. indicus x B. taurus cattle (x=2961 eggs; 2759 larvae). The number of eggs produced per unit engorged female weight is not significantly different in ticks that feed onB. taurus andB. indicus x B. taurus cattle.     

mtDNA Diversity and genetic lineages of eighteen cattle breeds from <Emphasis Type="Italic">Bos taurus </Emphasis>and<Emphasis Type="Italic"> Bos indicus</Emphasis> in China

In order to clarify the origin and genetic diversity of indigenous cattle breeds in China, we carried out phylogenetic analysis of representatives of those breeds by employing mitochondrial gene polymorphism. Complete cyt b gene sequences, 1140 bp in length, were determined for a total of 136 individuals from 18 different breeds and these sequences were clustered into two distinct genetic lineages: taurine (Bos taurus) and zebu (Bos indicus). In analysis of the cyt b gene diversity, Chinese cattle showed higher nucleotide (0.00923) and haplotype diversity (0.848) than the reports from other studies, and the animals from the taurine lineage indicated higher nucleotide diversity (0.00330) and haplotype diversity (0.746) than the ones from the zebu lineage (0.00136; 0.661). The zebu mtDNA dominated in the southern breeds (63.3–100%), while the taurine dominated in the northern breeds (81.8–100%). Six cattle breeds from the central area of China exhibited intermediate frequencies of zebu mtDNA (25–71.4%). This polymorphism revealed a declining south-to-north gradient of female zebu introgression and a geographical hybrid zone of Bos taurus and Bos indicus in China.     

Population structure and ancestry of Qinchuan cattle

The aim of this study was to estimate population structure and ancestry of Qinchuan cattle by genotyping 27 individuals using the GeneSeek HD 77k BeadChip, and another 1355 cattle representing breeds distributed worldwide, which had been genotyped using the Illumina Bovine 50k BeadChip. Qinchuan cattle were characterized by a dominant Bos taurus ancestry, accompanied by a considerable proportion of Bos indicus ancestry based on principal components analysis and supervised admixture analysis. A small proportion of Bos javanicus ancestry was detected as well. A similar admixture pattern in both Qinchuan and Turkish cattle breeds reflects their similar degrees of zebu introgression. Our study presents a relatively clear view of the population structure and ancestry of Qinchuan cattle, serving to benefit our understanding of this breed and leading to better targeted conservation approaches moving forward.     

Population growth rates forAmblyomma americanum (Acari: Ixodidae) onBos indicus,B. taurus andB. indicus x B. taurus cattle

Densities ofAmblyomma americanum (L.) onBos indicus, B. taurus andB. indicus x B. taurus cattle are compared over a 3-year period, and the growth rate (rate of increase or decrease) of parasitic tick populations on each cattle genotype is estimated.Average log₁₀ densities of parasiticA. americanum larvae are significantly (P=0.05) lower onB. indicus cattle than onB. taurus andB. indicus x B. taurus cattle. Average log densities of nymphal and adult ticks onB. taurus cattle are significantly higher than onB. indicus cattle but neither cattle genotype differs in this regard fromB. indicus x B. taurus cattle.Estimated annual tick population growth rates (log₁₀) for parasiticA. americanum are positive onB. taurus cattle (+0.84 larvae, +0.09 nymphs, +0.22 adults calf^–1 year^–1), but are negative onB. indicus (–0.18 nymphs, –0.14 adults calf^–1 year^–1) andB. indicus x B. taurus cattle (–0.45 larvae, –0.24 nymphs, –0.14 adults calf^–1 year^–1). Populations of parasitic larvae were not detected onB. indicus cattle.     

Tick-susceptible Bos taurus cattle display an increased cellular response at the site of larval Rhipicephalus (Boophilus) microplus attachment, compared with tick-resistant Bos indicus cattle

Cattle demonstrate divergent and heritable phenotypes of resistance and susceptibility to infestation with the cattle tick Rhipicephalus (Boophilus) microplus. Bos indicus cattle are generally more resistant to tick infestation than Bos taurus breeds although large variations in resistance can occur within subspecies and within breed. Increased tick resistance has been previously associated with an intense hypersensitivity response in B. taurus breeds; however, the mechanism by which highly resistant B. indicus cattle acquire and sustain high levels of tick resistance remains to be elucidated. Using the commercially available Affymetrix microarray gene expression platform, together with histological examination of the larval attachment site, this study aimed to describe those processes responsible for high levels of tick resistance in Brahman (B. indicus) cattle that differ from those in low-resistance Holstein-Friesian (B. taurus) cattle. We found that genes involved in inflammatory processes and immune responsiveness to infestation by ticks, although up-regulated in tick-infested Holstein-Friesian cattle, were not up-regulated in Brahman cattle. In contrast, genes encoding constituents of the extracellular matrix were up-regulated in Brahmans. Furthermore, the susceptible Holstein-Friesian animals displayed a much greater cellular inflammatory response at the site of larval R. microplus attachment compared with the tick-resistant Brahman cattle.     

Characterization of the 3' untranslated region of bovine parathyroid hormone gene in N'Dama and Boran cattle

We describe a polymorphism in the bovine gene PTHG which can be readily typed by PCR assay. The polymorphism is codominantly inherited and the allele frequencies appear characteristic of Bos indicus and B. taurus cattle.     

Assessment of cattle genetic introgression into domestic yak populations using mitochondrial and microsatellite DNA markers

Hybridization between yak Poephagus grunniens and taurine Bos taurus or indicine B. indicus cattle has been widely practiced throughout the yak geographical range, and gene flow is expected to have occurred between these species. To assess the impact of cattle admixture on domestic yak, we examined 1076 domestic yak from 29 populations collected in China, Bhutan, Nepal, India, Pakistan, Kyrgyzstan, Mongolia and Russia using mitochondrial DNA and 17 autosomal microsatellite loci. A cattle diagnostic marker‐based analysis reveals cattle‐specific mtDNA and/or autosomal microsatellite allele introgression in 127 yak individuals from 22 populations. The mean level of cattle admixture across the populations, calculated using allelic information at 17 autosomal microsatellite loci, remains relatively low (mY_cattle = 2.66 ± 0.53% and Q_cattle = 0.69 ± 2.58%), although it varies a lot across populations as well as among individuals within population. Although the level of cattle admixture shows a clear geographical structure, with higher levels of admixture in the Qinghai‐Tibetan Plateau and Mongolian and Russian regions, and lower levels in the Himalayan and Pamir Plateau region, our results indicate that the level of cattle admixture is not significantly correlated with the altitude across geographical regions as well as within geographical region. Although yak‐cattle hybridization is primarily driven to produce F₁ hybrids, our results show that the subsequent gene flow between yak and cattle took place and has affected contemporary genetic make‐up of domestic yak. To protect yak genetic integrity, hybridization between yak and cattle should be tightly controlled.     

When and how did Bos indicus introgress into Mongolian cattle?

The Mongolian cattle are one of the most widespread breeds with strictly Bos taurus morphological features in northern China. In our current study, we presented a diversity of mitochondrial DNA (mtDNA) D-loop region and Y chromosome SNP markers in 25 male and 8 female samples of Mongolian cattle from the Xinjiang Uygur autonomous region in Western China, and detected 21 B. taurus and four Bos indicus (zebu) mtDNA haplotypes. Among four B. indicus mtDNA haplotypes, two haplotypes belonged to I1 haplogroup and the remaining two haplotypes belonged to I2 haplogroup. In contrast, all 25 male Mongolian cattle samples revealed B. taurus Y chromosome haplotype and no B. indicus haplotypes were found. Historical and archeological records indicate that B. taurus was introduced to Xinjiang during the second millennium BC and B. indicus appeared in this region by the second century AD. The two types of cattle coexisted for many centuries in Xinjiang, as depicted in clay and wooden figurines unearthed in the Astana cemetery in Turfan (3rd–8th century AD). Multiple lines of evidence suggest that the earliest B. indicus introgression in the Mongolian cattle may have occurred during the 2nd–7th centuries AD through the Silk Road around the Xinjiang region. This conclusion differs from the previous hypothesis that zebu introgression to Mongolian cattle happened during the Mongol Empire era in the 13th century.