Genomic clues of the evolutionary history of Bos indicus cattle

Together with their sister subspecies Bos taurus, zebu cattle (Bos indicus) have contributed to important socioeconomic changes that have shaped modern civilizations. Zebu cattle were domesticated in the Indus Valley 8000 years before present (YBP). From the domestication site, they expanded to Africa, East Asia, southwestern Asia and Europe between 4000 and 1300 YBP, intercrossing with B. taurus to form clinal variations of zebu ancestry across the landmass of Afro‐Eurasia. In the past 150 years, zebu cattle reached the Americas and Oceania, where they have contributed to the prosperity of emerging economies. The zebu genome is characterized by two mitochondrial haplogroups (I1 and I2), one Y chromosome haplogroup (Y3) and three major autosomal ancestral groups (Indian‐Pakistani, African and Chinese). Phenotypically, zebu animals are recognized by their hump, large ears and excess skin. They are rustic, resilient to parasites and capable of bearing the hot and humid climates of the tropics. Many resources are available to study the zebu genome, including commercial arrays of SNP, reference assemblies and publicly available genotypes and whole‐genome sequences. Nevertheless, many of these resources were initially developed to support research and subsidize industrial applications in B. taurus, and therefore they can produce bias in data analysis. The combination of genomics with precision agriculture holds great promise for the identification of genetic variants affecting economically important traits such as tick resistance and heat tolerance, which were naturally selected for millennia and played a major role in the evolution of B. indicus cattle.     

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.     

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.     

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.     

A novel polymorphism of GDF5 gene and its association with body measurement traits in <Emphasis Type="Italic">Bos taurus</Emphasis> and <Emphasis Type="Italic">Bos indicus</Emphasis> breeds

Body measurement traits, influenced by genes and environmental factors, play numerous important roles in the value assessment of productivity and economy. Growth differentiate factor 5 (GDF5), involved in the development and maintenance of bone and cartilage, is an important candidate gene for body measurement traits selection through marker-assisted selection (MAS). In this study, based on the PCR-RFLP technology, we discovered and evaluated the potential association of the single nucleotide polymorphism (SNP) (T586C in exon 1) of the bovine GDF5 gene with body measurement traits in 985 Bos taurus breed, 42 Bos indicus breed and 76 Bos indicus × Bos taurus individuals. As the SNP marker, there were the significant effects on the Body length (BL) in the Bos taurus (BT) and Bos indicus × Bos taurus (BMY) populations (P < 0.05). In BT population, animals with the genotype TT had lower mean values for BL and Hip width (HW) than these with the TC and CC genotype (P < 0.01). In BMY population, animals with the genotype TC had lower mean values for BL than these with the genotype CC (P < 0.05). These results suggest that the SNP of the GDF5 gene could be a very useful genetic marker for body measurement traits in the bovine reproduction and breeding.     

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.     

Comprehensive analysis of the mitochondrial DNA diversity in Chinese cattle

Complete mitochondrial DNA D‐loop sequences of 1105 individuals were used to assess the diversity of maternal lineages of cattle populations in China. In total, 250 taurine and 88 zebu haplotypes were identified. Five main haplogroups—T1a, T2, T3, T4 and T5—were identified in Bos taurus, whereas Bos indicus harbored two haplogroups—I1 and I2. Our results suggest that the distribution of T1a in Asia was concentrated mainly in the northeast region (northeast China, Korea and Japan); haplogroups T2, T3 and T4 were predominant in Chinese cattle; and T5 was sporadically detected in Mongolian and Pingwu cattle. In contrast to the widespread presence of I1, I2 was distributed only in southwestern China (Yunnan‐Guizhou Plateau and the Tibet Autonomous Region) and Xinjiang Uygur Autonomous Region. This is the first time that all five taurine haplogroups and two zebu haplogroups have been found in Mongolian cattle. In addition, eight individuals in Tibetan cattle carried the Bos grunniens mtDNA type. The high mtDNA diversity (H = 0.904 ± 0.008) and the weak genetic structure among the 57 Chinese cattle breeds/populations are consistent with their complex historical background, migration route and ecological environment.     

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.     

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.     

A comparison of Ixodes rubicundus (Acari: Ixodidae) infestations on Friesian and Bonsmara cattle in South Africa

This study was conducted to compare the relative resistance of crossbred Bos indicus X B. taurus Bonsmara and B. taurus Friesian cattle to Ixodes rubicundus (Karoo paralysis tick) infestations. During periods of peak abundance of the ticks, Friesian oxen harboured almost twice or more than twice as many ticks as either Bonsmara oxen or cows. During periods of low tick abundance tick burdens on both cattle breeds were closely similar. It is envisaged that cattle can play an important role in an integrated control strategy against the Karro paralysis tick.     

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.     

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.     

Evidence for pleiotropism and recent selection in the PLAG1 region in Australian Beef cattle

A putative functional mutation (rs109231213) near PLAG1 (BTA14) associated with stature was studied in beef cattle. Data from 8199 Bos taurus, Bos indicus and Tropical Composite cattle were used to test the associations between rs109231213 and various phenotypes. Further, 23 496 SNPs located on BTA14 were tested for association with these phenotypes, both independently and fitted together with rs109231213. The C allele of rs109231213 significantly increased hip height, weight, net food intake, age at puberty in males and females and decreased IGF‐I concentration in blood and fat depth. When rs109231213 was fitted as a fixed effect in the model, there was an overall reduction in associations between other SNPs and these traits but some SNPs remained associated (P < 10^?4). Frequency of the mutant C allele of rs109231213 differed among B. indicus (0.52), B. taurus (0.96) and Tropical Composite (0.68). Most chromosomes carrying the C allele had the same surrounding 10 SNP haplotype, probably because the C allele was introgressed into Brahman from B. taurus cattle. A region of reduced heterozygosity surrounds the C allele; this is small in B. taurus but 20 Mb long in Brahmans, indicating recent and strong selection for the mutant allele. Thus, the C allele appears to mark a mutation that has been selected almost to fixation in the B. taurus breeds studied here and introduced into Brahman cattle during grading up and selected to a frequency of 0.52 despite its negative effects on fertility.     

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.     

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     

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.     

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.     

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.     

Local immune response against larvae of Rhipicephalus (Boophilus) microplus in Bos taurus indicus and Bos taurus taurus cattle

Bos taurus indicus cattle are less susceptible to infestation with Rhipicephalus (Boophilus) microplus than Bos taurus taurus cattle but the immunological basis of this difference is not understood. We compared the dynamics of leukocyte infiltrations (T cell subsets, B cells, major histocompatibility complex (MHC) class II-expressing cells, granulocytes) in the skin near the mouthparts of larvae of R. microplus in B. t. indicus and B. t. taurus cattle. Previously naïve cattle were infested with 50,000 larvae (B. t. indicus) or 10,000 larvae (B. t. taurus) weekly for 6 weeks. One week after the last infestation all of the animals were infested with 20,000 larvae of R. microplus. Skin punch biopsies were taken from all animals on the day before the primary infestation and from sites of larval attachment on the day after the first, second, fourth and final infestations. Infiltrations with CD3⁺, CD4⁺, CD8⁺ and γδ T cells followed the same pattern in both breeds, showing relatively little change during the first four weekly infestations, followed by substantial increases at 7 weeks post-primary infestation. There was a tendency for more of all cell types except granulocytes to be observed in the skin of B. t. indicus cattle but the differences between the two breeds were consistently significant only for γδ T cells. Granulocyte infiltrations increased more rapidly from the day after infestation and were higher in B. t. taurus cattle than in B. t. indicus. Granulocytes and MHC class II-expressing cells infiltrated the areas closest to the mouthparts of larvae. A large volume of granulocyte antigens was seen in the gut of attached, feeding larvae.