Introgression,admixture, and selection facilitate genetic adaptation to high-altitude environments in cattle

Domestication and subsequent selection of cattle to form breeds and biological types that can adapt to different environments partitioned ancestral genetic diversity into distinct modern lineages. Genome-wide selection particularly for adaptation to extreme environments left detectable signatures genome-wide. We used high-density genotype data for 42 cattle breeds and identified the influence of Bos grunniens and Bos javanicus on the formation of Chinese indicine breeds that led to their divergence from India-origin zebu. We also found evidence for introgression, admixture, and migration in most of the Chinese breeds. Selection signature analyses between high-altitude (≥1800 m) and low-altitude adapted breeds (<1500 m) revealed candidate genes (ACSS2, ALDOC, EPAS1, EGLN1, NUCB2) and pathways that are putatively involved in hypoxia adaptation. Immunohistochemical, real-time PCR and CRISPR/cas9 ACSS2-knockout analyses suggest that the up-regulation of ACSS2 expression in the liver promotes the metabolic adaptation of cells to hypoxia via the hypoxia-inducible factor pathway. High altitude adaptation involved the introgression of alleles from high-altitude adapted yaks into Chinese Bos taurus taurus prior to their formation into recognized breeds and followed by selection. In addition to selection, adaptation to high altitude environments has been facilitated by admixture and introgression with locally adapted cattle populations.     

Expression pattern of bta-mir-2898 miRNA and their correlation with heat shock proteins during summer heat stress among native vs crossbred cattle

Earlier studies identified the role of bta-mir-2898 in bovine. Our earlier study identified that, bta-mir-2898 can be over expressed in crossbred cattle during heat stress. Nevertheless the differential expression of bta-mir-2898 among native vs crossbred cattle during summer stress along with it's correlation with different heat shock proteins (HSPs) is not yet studied. In the present context, we studied the differential expression of bta-mir-2898 among Frieswal (Bos indicus x Bos taurus) and Sahiwal (Bos indicus) breeds of cattle during a range of environmental air temperatures and further investigated the correlation of bta-mir-2898 with different HSPs (HSP70, HSP90, HSP60. HSF, HSPB8 and HSP27). It was observed that, at peak air temperature the relative miRNA expression level (p < 0.05) of bta-mir-2898 was 3.4 ± 0.41 and 0.79 ± 0.22 among Frieswal and Sahiwal, respectively. We also observed significant levels (p < 0.05) of mRNA abundance of HSP70, HSP90, HSPB8 and HSP27 among the breeds. In all the cases Sahiwal found to exhibited higher level of HSPs in comparison to Frieswal. Studies revealed that the expression profile of bta-mir-2898 was negatively correlated with the expression of all the HSPs during thermal stress in post anti-mir2898 treated PBMC invitro cultured model originated from both Frieswal and Sahiwal cattle breeds. However, significantly (p < 0.05) higher negative correlations were observed between bta-mir-2898 and HSP70, HSP60 and HSPB8. Present findings highlighted the preliminary role of overexpressed bta-mir-2898 in cattle during thermal stress and its impact on different heat shock proteins.     

Antimicrobial Susceptibility Profiles of <Emphasis Type="Italic">Salmonella</Emphasis><Emphasis Type="Italic">enterica</Emphasis> Serotypes Recovered from Pens of Commercial Feedlot Cattle Using Different Types of Composite Samples

Salmonella enterica in cattle production systems may be associated with important human and animal disease issues. However, tremendous diversity exists among Salmonella recovered, and more information is needed about strains of greatest potential health concern, particularly those that are multidrug resistant (MDR). By characterizing Salmonella isolates from commercial feedlot pens, this study aimed to evaluate the strain diversity and prevalence of MDR Salmonella from different types of composite pen samples. Antimicrobial susceptibility profiles, serotype, and presence or absence of the integron-encoded intI1 gene were determined for 530 Salmonella isolates recovered using composite rope (n = 335), feces (n = 59), and water (n = 136) samples from 21 pens in 3 feedlots. The study investigated only pens with available isolates from multiple sample types. Most isolates (83.0%) of the 19 Salmonella serotypes identified were susceptible or intermediately susceptible to all the antimicrobials evaluated. Resistance to sulfisoxazole (14.9%), streptomycin (3.8%), and tetracycline (3.6%) were the most common. None of the isolates tested positive for a class 1 integron, and only 2.5% were resistant to multiple antimicrobials. All the MDR isolates, namely, serotypes Uganda (n = 9), Typhimurium (n = 2), and Give (n = 2), were resistant to at least five antimicrobials. Most MDR isolates (n = 11) were from two pens during 1 week within one feedlot. Overall, many Salmonella isolates collected within a pen were similar in terms of serotype and antimicrobial susceptibility regardless of sample type. However, MDR Salmonella and rare serotypes were not recovered frequently enough to suggest a general strategy for appropriate composite sampling of feedlot cattle populations for Salmonella detection and monitoring.     

中国部分地方黄牛品种线粒体D-loop区的遗传变异与进化分析

测定了13个黄牛品种125个个体的线粒体D-loop区段的全序列,包括12个中国地方黄牛品种的123个个体和德国黄牛2个个体,并进行了分析。结果显示,共检测到93个变异位点,57个单倍型,平均核苷酸差异(average number ofnucleotide differences,k)为22.708,核苷酸多样度(nucleotide diversity,π)为0.0251±0.00479,单倍型多样度(haplotypediversity,Hd)为0.888±0.026,表明我国黄牛品种遗传多样性非常丰富。构建的Neighbor-Joining进化树显示这13个品种主要分成两大类型:普通牛和瘤牛;新发现的特殊类型Ⅲ只有一个西藏阿沛甲咂牛的个体,它与牦牛D-loop序列最相近,证明西藏地区的黄牛与牦牛之间存在基因渗入现象。普通牛和瘤牛在日喀则驼峰牛中占的比例分别是64.3%和35.7%,在阿沛甲咂牛中占的比例分别是50.0%和50.0%,证明了西藏的黄牛也有瘤牛类型。云南牛品种的单倍型非常丰富证明了云南在中国黄牛起源上的重要地位;在27个中国黄牛品种中(本研究11个品种以及GenBank上的16个品种)找到了中国瘤牛的核心单倍型i1,并且对它进行了讨论。同时证明了西藏瘤牛独立于中国瘤牛核心类群的特殊性。     

Haematobia irritans parasitism of F1 yak × beef cattle (Bos grunniens × Bos taurus) hybrids

The horn fly Haematobia irritans (Diptera: Muscidae) is a blood obligate ectoparasite of bovids that causes annual losses to the U.S. beef cattle industry of over US$1.75 billion. Climate warming, the anthropogenic dispersion of bovids and the cross‐breeding of beef cattle with other bovid species may facilitate novel horn fly–host interactions. In particular, hybridizing yaks [Bos grunniens (Artiodactyla: Bovidae)] with beef cows (Bos taurus) for heterosis and carcass improvements may increase the exposure of yak × beef hybrids to horn flies. The present paper reports on the collection of digital images of commingled beef heifers (n = 12) and F1 yak × beef hybrid bovids (heifers, n = 7; steers, n = 5) near Laramie, Wyoming (～ 2200 m a.s.l.) in 2018. The total numbers of horn flies on beef heifers and F1 yak × beef heifers [mean ± standard error (SE): 88 ± 13 and 70 ± 17, respectively] did not differ significantly; however, F1 yak × beef steers had greater total horn fly abundance (mean ± SE: 159 ± 39) than female bovids. The present report of this experiment is the first such report in the literature and suggests that F1 yak × beef bovids are as susceptible as cattle to horn fly parasitism. Therefore, similar monitoring and treatment practices should be adopted by veterinarians, entomologists and producers.     

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.     

Genetic parameters of infectious bovine keratoconjunctivitis and its relationship with weight and parasite infestations in Australian tropical Bos taurus cattle

Background

Infectious bovine keratoconjunctivitis (IBK) or ‘pinkeye’ is an economically important ocular disease that significantly impacts animal performance. Genetic parameters for IBK infection and its genetic and phenotypic correlations with cattle tick counts, number of helminth (unspecified species) eggs per gram of faeces and growth traits in Australian tropically adapted Bos taurus cattle were estimated.

Methods

Animals were clinically examined for the presence of IBK infection before and after weaning when the calves were 3 to 6 months and 15 to 18 months old, respectively and were also recorded for tick counts, helminth eggs counts as an indicator of intestinal parasites and live weights at several ages including 18 months.

Results

Negative genetic correlations were estimated between IBK incidence and weight traits for animals in pre-weaning and post-weaning datasets. Genetic correlations among weight measurements were positive, with moderate to high values. Genetic correlations of IBK incidence with tick counts were positive for the pre-weaning and negative for the post-weaning datasets but negative with helminth eggs counts for the pre-weaning dataset and slightly positive for the post-weaning dataset. Genetic correlations between tick and helminth eggs counts were moderate and positive for both datasets. Phenotypic correlations of IBK incidence with helminth eggs per gram of faeces were moderate and positive for both datasets, but were close to zero for both datasets with tick counts.

Conclusions

Our results suggest that genetic selection against IBK incidence in tropical cattle is feasible and that calves genetically prone to acquire IBK infection could also be genetically prone to have a slower growth. The positive genetic correlations among weight traits and between tick and helminth eggs counts suggest that they are controlled by common genes (with pleiotropic effects). Genetic correlations between IBK incidence and tick and helminth egg counts were moderate and opposite between pre-weaning and post-weaning datasets, suggesting that the environmental and (or) maternal effects differ between these two growth phases. This preliminary study provides estimated genetic parameters for IBK incidence, which could be used to design selection and breeding programs for tropical adaptation in beef cattle.     

Copy number variation detection in Chinese indigenous cattle by whole genome sequencing

Copy number variation (CNV) refers to a kind of structural variation, having functional and evolutionary effects on phenotypes. Thus far, further elucidation of the CNVs in different Chinese indigenous cattle breeds by whole genome sequencing have yet not been done. In this study, a comprehensive genomic analysis was performed on 75 cattle individuals including six Chinese indigenous cattle breeds and two non-native specialized beef cattle breeds. Based on the 11,486 CNVRs discovered, population analysis was performed, showed that all the cattle breeds clustered in to three clades, consistent with their lineages Bos taurus, Bos taurus × Bos indicus and Bos indicus. Importantly, a set of CNVRs related genes were found to be associated with the traits of interest, which include meat production or quality (CAST, ACTC1, etc.), adaption (BLA-DQB, EGLN2, etc.) and coat color (KIT, MITF, etc.). These results provide valuable full genome variation resources for Chinese bovine genome research and would be helpful for cattle breeding and selection programs in the future.     

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.     

Phylogenetic relationships among the European and American bison and seven cattle breeds reconstructed using the BovineSNP50 Illumina Genotyping BeadChip

Here we present the first attempt to use the BovineSNP50 Illumina Genotyping BeadChip for genome-wide screening of European bison Bison bonasus bonasus (EB), two subspecies of American bison: the plains bison Bison bison bison (PB), the wood bison Bison bison athabascae (WB) and seven cattle Bos taurus breeds. Our aims were to (1) reconstruct their evolutionary relationships, (2) detect any genetic signature of past bottlenecks and to quantify the consequences of bottlenecks on the genetic distances amongst bison subspecies and cattle, and (3) detect loci under positive or stabilizing selection. A Bayesian clustering procedure (STRUCTURE) detected ten genetically distinct clusters, with separation among all seven cattle breeds and European and American bison, but no separation between plain and wood bison. A linkage disequilibrium based program (LDNE) was used to estimate the effective population size (N _e) for the cattle breeds; N _e was generally low, relative to the census size of the breeds (cattle breeds: mean N _e = 299.5, min N _e = 18.1, max N _e = 755.0). BOTTLENECK 1.2 detected signs of population bottlenecks in EB, PB and WB populations (sign test and standardized sign test: p = 0.0001). Evidence for loci under selection was found in cattle but not in bison. All extant wild populations of bison have shown to have survived severe bottlenecks, which has likely had large effects on genetic diversity within and differentiation among groups.     

Geographic distribution of haplotype diversity at the bovine casein locus

The genetic diversity of the casein locus in cattle was studied on the basis of haplotype analysis. Consideration of recently described genetic variants of the casein genes which to date have not been the subject of diversity studies, allowed the identification of new haplotypes. Genotyping of 30 cattle breeds from four continents revealed a geographically associated distribution of haplotypes, mainly defined by frequencies of alleles at CSN1S1 and CSN3. The genetic diversity within taurine breeds in Europe was found to decrease significantly from the south to the north and from the east to the west. Such geographic patterns of cattle genetic variation at the casein locus may be a result of the domestication process of modern cattle as well as geographically differentiated natural or artificial selection. The comparison of African Bos taurus and Bos indicus breeds allowed the identification of several Bos indicus specific haplotypes (CSN1S1*C-CSN2*A²-CSN3*A_I/CSN3*H) that are not found in pure taurine breeds. The occurrence of such haplotypes in southern European breeds also suggests that an introgression of indicine genes into taurine breeds could have contributed to the distribution of the genetic variation observed.     

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.     

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.     

Genetic diversity,introgression and relationships among West/Central African cattle breeds

Genetic diversity, introgression and relationships were studied in 521 individuals from 9 African Bos indicus and 3 Bos taurus cattle breeds in Cameroon and Nigeria using genotype information on 28 markers (16 microsatellite, 7 milk protein and 5 blood protein markers). The genotypes of 13 of the 16 microsatellite markers studied on three European (German Angus, German Simmental and German Yellow) and two Indian (Nelore and Ongole) breeds were used to assess the relationships between them and the African breeds. Diversity levels at microsatellite loci were higher in the zebu than in the taurine breeds and were generally similar for protein loci in the breeds in each group. Microsatellite allelic distribution displayed groups of alleles specific to the Indian zebu, African taurine and European taurine. The level of the Indian zebu genetic admixture proportions in the African zebus was higher than the African taurine and European taurine admixture proportions, and ranged from 58.1% to 74.0%. The African taurine breed, Muturu was free of Indian zebu genes while its counter Namchi was highly introgressed (30.2%). Phylogenic reconstruction and principal component analysis indicate close relationships among the zebu breeds in Cameroon and Nigeria and a large genetic divergence between the main cattle groups – African taurine, European taurine and Indian zebu, and a central position for the African zebus. The study presents the first comprehensive information on the hybrid composition of the individual cattle breeds of Cameroon and Nigeria and the genetic relationships existing among them and other breeds outside of Africa. Strong evidence supporting separate domestication events for the Bos species is also provided.     

Genetics and domestic cattle origins

Genetics has the potential to provide a novel layer of information pertaining to the origins and relationships of domestic cattle. While it is important not to overstate the power of archeological inference from genetic data, some previously widespread conjectures are inevitably contradicted with the addition of new information. Conjectures regarding domesticated cattle that fall into this category include a single domestication event with the development of Bos indicus breeds from earlier Bos taurus domesticates; the domestication of a third type of cattle in Africa having an intermediate morphology between the two taxa; and the special status of the Jersey breed as a European type with some exotic influences. In reality, a wideranging survey of the genetic variation of modern cattle reveals that they all derive from either zebu or taurine progenitors or are hybrids of the two. The quantitative divergence between Bos indicus and Bos taurus strongly supports a predomestic separation; that between African and European taurines also suggests genetic input from native aurochsen populations on each continent. Patterns of genetic variants assayed from paternally, maternally, and biparentally inherited genetic systems reveal that extensive hybridization of the two subspecies is part of the ancestry of Northern Indian, peripheral European, and almost all African cattle breeds. In Africa, which is the most extensive hybrid zone, the sexual asymmetry of the process of zebu introgression into native taurine breeds is strikingly evident. © 1998 Wiley-Liss, Inc.     

Genetic diversity in Iranian indigenous sheep vis-à-vis selected exogenous sheep breeds and wild mouflon

The heterogeneity of climate and different agro-ecological conditions in Iran have resulted in development of 27 indigenous sheep breeds. Wild Asiatic mouflon (Ovis orientalis) is believed to be the ancestor of Iranian sheep. Evaluation of genetic diversity and population structure within and among domestic breeds has important implications for animal breeding programs and genetic resources management. Based on 50K SNP genotype data, we studied the genetic diversity of five indigenous Iranian sheep breeds: Afshari (n = 37), Moghani (n = 34), Qezel (n = 35), Zel (n = 46) and Lori-Bakhtiari (n = 46), and Asiatic mouflon (n = 8) sampled from Iran. Furthermore, genetic diversity and the breed admixture of Iranian sheep were assessed on a larger geographic scale using a reference panel comprising: three indigenous Afghan breeds – Arabi (n = 15), Balouchi (n = 15) and Gadik (n = 15); three indigenous breeds from Turkey and Cyprus – Cyprus Fat Tail (n = 30), Karakas (n = 18) and Norduz (n = 20); and three commercial European breeds – Suffolk (n = 19), Comisana (n = 24) and Engadine Red Sheep (n = 24). The results revealed that the investigated breeds are divided into five genetically distinct clusters according to their geographic origin. Afshari was closest to the local mouflon population and showed signs of mouflon admixture. Qezel was identified as a hybrid sheep breed. Much evidence supported the Afghan breeds being identical. Inbreeding values, which were estimated based on ROHs, were highest for Suffolk (F_ROH = 0.0544) and lowest for Balouchi (F_ROH = 0.0078). In conclusion, analysis of selected breeds from neighboring countries along with Asiatic mouflon gave a deeper insight into the evolutionary history and origin of Iranian sheep with important implications for future breed management.     

Prevalence and Antimicrobial Resistance of Thermophilic Campylobacter spp. from Cattle Farms in Washington State

The prevalence of thermophilic Campylobacter spp. was investigated in cattle on Washington State farms. A total of 350 thermophilic Campylobacter isolates were isolated from 686 cattle sampled on 15 farms (eight dairies, two calf rearer farms, two feedlots, and three beef cow-calf ranches). Isolate species were identified with a combination of phenotypic tests, hipO colony blot hybridization, and multiplex lpxA PCR. Breakpoint resistance to four antimicrobials (ciprofloxacin, nalidixic acid, erythromycin, and doxycycline) was determined by agar dilution. Campylobacter jejuni was the most frequent species isolated (34.1%), followed by Campylobacter coli (7.7%) and other thermophilic campylobacters (1.5%). The most frequently detected resistance was to doxycycline (42.3% of 350 isolates). Isolates from calf rearer facilities were more frequently doxycycline resistant than isolates from other farm types. C. jejuni was most frequently susceptible to all four of the antimicrobial drugs studied (58.8% of 272 isolates). C. coli isolates were more frequently resistant than C. jejuni, including resistance to quinolone antimicrobials (89.3% of isolates obtained from calves on calf rearer farms) and to erythromycin (72.2% of isolates obtained from feedlot cattle). Multiple drug resistance was more frequent in C. coli (51.5%) than in C. jejuni (5.1%). The results of this study demonstrate that C. jejuni is widely distributed among Washington cattle farms, while C. coli is more narrowly distributed but significantly more resistant.     

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.