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.     

Application of mathematical expectation (ME) strategy for detecting low frequency mutations: An example for evaluating 14-bp insertion/deletion (indel) within the bovine PRNP gene

The detection method based on the mathematical expectation (ME) strategy is fast and accuracy for low frequency mutation screening in large samples. Previous studies have found that the 14-bp insertion/deletion (indel) variants of the 3′ untranslated region (3′ UTR) within bovine PRNP gene have been characterized with low frequency (≤5%) in global breeds outside China, which has not been determined in Chinese cattle breeds yet. Therefore, this study aimed to identify the 14-bp indel within PRNP gene in 5 major Chinese indigenous cattle breeds and to evaluate its associations with phenotypic traits. It was the first time to use ME strategy to detect low frequency indel polymorphisms and found that minor allele frequency was 0.038 (Qinchuan), 0.033 (Xianan), 0.013 (Nanyang), 0.003 (Jiaxian), and zero (Ji'an), respectively. Compared to the traditional detection method by which the sample was screened one by one, the reaction time by using the ME method was decreased 62.5%, 64.9%, 77.6%, 88.9% and 66.4%, respectively. In addition, the 14-bp indel was significantly associated with the growth traits in 2 cattle breeds, with the body length of Qinchuan cattle as well as the body weight and waistline of Xianan cattle. Our results have uncovered that the method based on ME strategy is rapid, reliable, and cost-effective for detecting the low frequency mutation as well as our findings provide a potential valuable theoretical basis for the marker-assisted selection (MAS) in beef cattle.     

Polymorphism of the prion protein gene (<Emphasis Type="Italic">PRNP</Emphasis>) in two Chinese indigenous cattle breeds

Prion protein (PRNP) gene has been located at position q17 of chromosome 13 in cattle. The polymorphisms of PRNP gene might be associated with BSE susceptibility. In the present work, we investigated the polymorphisms of PRNP gene, including SNP in exon 3, 23-bp indel in promoter region, 12-bp indel in intron 1 in 2 Chinese indigenous cattle breeds of northeast China. Eighty-six animals from Yanbian (34) and Chinese Red Steppes (52) were genotyped at PRNP locus by analyzing genomic DNA. A total of 4 single nucleotide polymorphism (SNP) sites were revealed in the PRNP gene exon 3 of the 2 cattle breeds investigated. Three of these SNPs were non-synonymous mutations that resulted in the amino acid exchanges (K119N, S154N, and M177V), and one is silent nucleotide substitutions (A234G). The two amino acid mutations of S154N and M177V were detected only in Yanbian with a very low frequency (0.0147), and they appears to be absent in Chinese Red Steppes. The average gene heterozygosity (H _e), effective allele numbers (N _e), Shannon’s information index (I) and polymorphism information content (PIC) were 0.3088, 1.5013, 0.3814 and 0.2000 in Yanbian, respectively, being relatively higher than that of Chinese Red Steppes (0.2885, 1.4985, 0.3462 and 0.1873, respectively). In 23-bp indel and 12-bp indel loci, three different genotypes were identified in both Yanbian and Chinese Red Steppes breeds. Based 23- and 12-bp indels, four haplotypes was constructed in the 2 Chinese cattle breeds, of which the 23-bp (−)/12-bp (−) was main haplotypes accounting for more than 50% of the total in both Yanbian and Chinese Red Steppes breeds. These results might be useful in understanding the genetic characteristics of PRNP gene in Chinese indigenous cattle breeds.     

Haplotype analysis of TLR4 gene and its effects on milk somatic cell score in Chinese commercial cattle

Bovine mastitis is a very complex and common disease of dairy cattle and a major source of economic losses to the dairy industry worldwide. In this study, the bovine TLR4 was taken as a candidate gene for mastitis resistance. This study aimed to analyze the associations of single nucleotide polymorphisms (SNP) or haplotype and somatic cell score (SCS) in 404 Chinese commercial dairy cattle including Chinese Holstein, Sanhe cattle and Chinese Simmental breeds. The polymerase chain reaction and sequencing methods were used for detecting genotype and allele frequency distribution of the two SNPs (rs8193062, rs8193064), statistical results showed that T allele at rs8193062 and C allele at rs8193064 were the predominate alleles. Moreover, six SNPs, including two SNPs (rs8193062, rs8193064) and four SNPs (rs8193060, rs8193069, rs29017188, rs8193046) which were chosen according the polymorphism level for the same cattle populations in previous studies, were used for haplotype analysis, the results revealed that twenty-one haplotypes were found in the mentioned animals, of which, Hap1 (30.5 %) and Hap2 (30.4 %) were the most common haplotypes. Hap2, Hap4 and Hap12 might negatively effect on milk SCS, whereas Hap13 might positively effect on milk SCS. The results in this study might assist in marker assisted selection and provided some reference to be implemented in breeding programs to improve the mastitis resistance of dairy cattle.     

Variation in the XKR4 gene was significantly associated with subcutaneous rump fat thickness in indicine and composite cattle

Variation in the XK, Kell blood group complex subunit–related family, member 4 (XKR4) gene on BTA14 was associated with rump fat thickness in a recent genome‐wide association study. This region is also of interest because it is known to show evidence of a signature of population genetic selection. In this study, additional variation in this gene was genotyped in a sample of a total of 1283 animals of the Belmont Red (BEL) and Santa Gertrudis (SGT) breeds. The SNP rs41724387 was significantly (P < 0.001) associated with rump fat thickness and explained 5.9% of the genetic variance for the trait in this sample. Using the 4466 genotypes for the SNP rs42646708 from several data sets to estimate effects in seven breeds, this relatively large quantitative trait locus effect appears to be a result of the variation in indicine and taurine–indicine composite cattle. However, the only DNA variant found in Brahman cattle that altered the predicted amino acid sequence of XKR4 was not associated with rump fat thickness. This suggests that causative mutations lie outside the coding sequence of this gene.     

<Emphasis Type="Italic">IGF-1</Emphasis> and <Emphasis Type="Italic">IGF-1R</Emphasis> gene polymorphisms in East Anatolian Red and South Anatolian Red cattle breeds

The aim of this study was to determine the allele and genotype frequencies of the Insulin like growth factor 1 gene (IGF-1) and Insulin like growth factor 1 receptor gene (IGF-1R) polymorphisms in East Anatolian (EAR) and South Anatolian Red cattle (SAR). Polymorphisms in both genes are claimed to affect economic parameters like body weight and subcutaneous back fat. For the study, blood samples were collected from 50 SAR and 50 EAR cattle. In both breeds, high frequencies of allele B of IGF-1, which was supposed to have positive effect on carcass traits and allele A of IGF-1R, related to milk traits, were observed. Therefore, no relationship of the polymorphisms studied with economic traits was observed, as both breeds have respectively low carcass and milk parameters. However, distribution of allele frequencies of IGF-1 and IGF-1R in SAR and EAR cattle was similar with Zebu cattle that support Bos indicus introgression to Anatolian breeds.     

Genetic Variant of MYLK4 Gene and its Association with Growth Traits in Chinese Cattle

As a member of MYLK family, MYLK4 gene may play a vital role in muscle development. In this study, one novel single-nucleotide polymorphism (SNP) was identified the bovine MYLK4 by sequencing pooled DNA samples (pool-Seq) and forced polymerase chain reaction-restriction fragment length polymorphism (forced PCR-RFLP) methods. Overall, we reported one mutation (SNP1) in the intron 10 region within the bovine MYLK4 gene in 559 individuals representing five main cattle breeds from China (Nanyang, NY; Qinchuan; Jiaxian, JX; Pinan cattle; and Caidamu cattle, CDM). Genotype AA and allele A were predominant in the QC, PN, and XN populations. Association analysis with growth traits in the QC breed showed that the animals with genotype GG had significantly greater chest breadth and hip width (P?<?0.05). Meanwhile, the genotype GG was strongly associated with withers height and body length than those with genotype AA (P?<?0.01 or P?<?0.05) at 12 months in the NY breed. These statistical results exhibited that the MYLK4 gene might be a potential candidate gene to improve cattle’s growth traits, and the SNP could be used as molecular markers in early marker-assisted selection (MAS) in beef cattle breeding program.     

Genetic diversity of Chinese cattle revealed by Y‐SNP and Y‐STR markers

With its vast territory and complex natural environment, China boasts rich cattle genetic resources. To gain the further insight into the genetic diversity and paternal origins of Chinese cattle, we analyzed the polymorphism of Y‐SNPs (UTY19 and ZFY10) and Y‐STRs (INRA189 and BM861) in 34 Chinese cattle breeds/populations, including 606 males representative of 24 cattle breeds/populations collected in this study as well as previously published data for 302 bulls. Combined genotypic data identified 14 Y‐chromosome haplotypes that represented three haplogroups. Y2‐104‐158 and Y2‐102‐158 were the most common taurine haplotypes detected mainly in northern and central China, whereas the indicine haplotype Y3‐88‐156 predominates in southern China. Haplotypes Y2‐108‐158, Y2‐110‐158, Y2‐112‐158 and Y3‐92‐156 were private to Chinese cattle. The population structure revealed by multidimensional scaling analysis differentiated Tibetan cattle from the other three groups of cattle. Analysis of molecular variance showed that the majority of the genetic variation was explained by the genetic differences among groups. Overall, our study indicates that Chinese cattle retain high paternal diversity (H = 0.607 ± 0.016) and probably much of the original lineages that derived from the domestication center in the Near East without strong admixture from commercial cattle carrying Y1 haplotypes.     

Fine‐mapping the POLL locus in Brahman cattle yields the diagnostic marker CSAFG29

The POLL locus has been mapped to the centromeric region of bovine chromosome 1 (BTA1) in both taurine breeds and taurine–indicine crosses in an interval of approximately 1 Mb. It has not yet been mapped in pure‐bred zebu cattle. Despite several efforts, neither causative mutations in candidate genes nor a singular diagnostic DNA marker has been identified. In this study, we genotyped a total of 68 Brahman cattle and 20 Hereford cattle informative for the POLL locus for 33 DNA microsatellites, 16 of which we identified de novo from the bovine genome sequence, mapping the POLL locus to the region of the genes IFNAR2 and SYNJ1. The 303‐bp allele of the new microsatellite, CSAFG29, showed strong association with the POLL allele. We then genotyped 855 Brahman cattle for CSAFG29 and confirmed the association between the 303‐bp allele and POLL. To determine whether the same association was found in taurine breeds, we genotyped 334 animals of the Angus, Hereford and Limousin breeds and 376 animals of the Brangus, Droughtmaster and Santa Gertrudis composite taurine–zebu breeds. The association between the 303‐bp allele and POLL was confirmed in these breeds; however, an additional allele (305 bp) was also associated but not fully predictive of POLL. Across the data, CSAFG29 was in sufficient linkage disequilibrium to the POLL allele in Australian Brahman cattle that it could potentially be used as a diagnostic marker in that breed, but this may not be the case in other breeds. Further, we provide confirmatory evidence that the scur phenotype generally occurs in animals that are heterozygous for the POLL allele.     

Genome-wide assessment of genetic diversity and population structure insights into admixture and introgression in Chinese indigenous cattle

Background

China exhibits a great diversity of ecosystems and abundant cattle resources, with nearly 30 million cattle from 53 indigenous breeds reared in specific geographic regions. To explore the genetic diversity and population structure of Chinese indigenous cattle, a population genetic analysis at both the individual and population levels was conducted and the admixture analysis was performed. We genotyped 572 samples from 20 Chinese indigenous cattle breeds using GeneSeek Genomic Profiler Bovine LD (GGP-LD, 30?K) and downloaded the published data of 77 samples from 4 worldwide commercial breeds genotyped with Illumina BovineSNP50 Beadchip (SNP50, 50?K).

Results

In principal component analysis (PCA) and neighbour-joining (NJ) tree analysis, samples of the same breeds were grouped together, leading to clear separation from other breeds. And Chinese indigenous cattle were clustered into two groups of southern and northern breeds, originated from Asian Bos indicus lineage and Eurasian Bos taurus lineage, respectively. In STRUCTURE K?=?2, a clear transition occurred from the northern breeds to the southern breeds. Additionally, the northern breeds contained a smaller Eurasian taurine (62.5%) descent proportion than that reported previously (more than 90%). In STRUCTURE K?=?3, a distinct descent was detected in the southern Chinese breeds, which could reflect a long-term selection history of Chinese indigenous cattle. The results from TreeMix and f3 statistic provided the evidence of an admixture history between southern breeds and northern breeds.

Conclusions

Consistent with the observed geographical distributions, Chinese indigenous cattle were divided into two genetic clusters, northern indigenous cattle and southern indigenous cattle. Three improved breeds in the northern area also exhibited northern indigenous ancestry. We found that the breeds distributed in the northern China showed more southern lineage introgression than previously reported. Central-located populations appeared to the admixture between southern and northern lineages, and introgression events from European cattle were observed in Luxi Cattle, Qinchuan Cattle and Jinnan Cattle. The study revealed the population structures and levels of admixture pattern among Chinese indigenous cattle, shedding light on the origin and evolutionary history of these breeds.

     

Polymorphism of Bovine Prolactin Gene: Microsatellites,PCR–RFLP

In the samples of Russian Ayrshire and Gorbatov Red cattle breeds, distribution of frequencies of prolactin (PRL) gene alleles generated due to the presence of polymorphic RsaI site in exon 3 were studied. In the breeds, the frequencies of the Ballele of the PRLgene (with RsaI(+) site) detected by the PCR–RFLP method were 14.1 and 8.6%, respectively. In Black Pied, Ayrshire and Gorbatov Red cattle breeds, variation of the microsatellite dinucleotide repeat in the regulatory region of the gene PRLwas also studied. Gorbatov Red breed was monomorphic at the microsatellite locus with the only allele 164 bp in length. Two alleles (164 bp and 162 bp) were detected in the other breeds studied. The frequencies of 164-bp allele of the microsatellite locus were 93.7 and 90.0% in Black Pied and Ayrshire breeds, respectively. In Gorbatov Red breed of dairy type with good beef qualities and low milk-fat yield, lower level of heterozygosity for PRLgene was demonstrated compared to Ayrshire and Black Pied breeds that have high milk-fat yield. In three cattle breeds, higher mean estimate of polymorphism information content of PCR–RFLP in exon 3 (PIC = 0.21) was revealed compared with the same estimate (PIC = 0.09) for the microsatellite locus variability in the regulatory region of the PRLgene. Characteristics of allele Bdistribution of thePRLgene in the representatives of the Bovidae family are considered.     

A genome‐wide association study reveals a quantitative trait locus for age at first calving in delta/notch‐like EGF repeat containing on chromosome 2 in Japanese Black cattle

Age at first calving (AFC) is an important trait for achieving earlier reproductive performance in cattle. To identify quantitative trait loci for AFC in Japanese Black cattle, we conducted a genome‐wide association study using 866 animals with extreme AFC values selected from a larger group of 52 009 animals. We identified single nucleotide polymorphisms (SNPs) on bovine chromosome 2 that were associated with AFC. These SNPs were located within 112.8‐kbp intronic region of delta/notch‐like EGF repeat containing (DNER) and proved to be in a state of high linkage disequilibrium. The association was replicated in an independent sample of 2963 animals. In the replicated population, the frequency of the reduced AFC allele (Q) was 0.463, and the allele accounts for 8% of the total genetic variance. The effect of allele substitution on AFC was a decrease of 11.54 days. The results suggest that the Q allele could serve as a useful marker in Japanese Black cattle to select animals with superior AFC performance.     

Genetic Variants in ADD1 Gene and their Associations with Growth Traits in Cattle

The α-adducin (ADD1) is a subunit of adducin which is a cytoskeleton heterodimeric protein. Adducin participates in oocytes chromosome meiosis of mice, prompting adducin has an effect on embryonic development. Adducin gene mutation has significantly functional change. So the present study was to identify and characterize polymorphisms within the coding region of the bovine ADD1 gene among different cattle breeds. Here, 11 novel single nucleotide polymorphisms (SNPs 1–11) were identified by DNA sequencing and polymerase chain reaction-single stranded conformational polymorphism, there were one synonymous mutation in exon 1 (SNP1); four missense mutations in exons 4, 7, and 8 (SNPs 3–6); and six mutations in introns 4, 12, 13, and 14 (SNPs 2, 7–10). The statistical analyses indicated that the some SNPs are associated with the growth traits (body length, body height, chest circumference, and hucklebone width) in Chinese Jiaxian cattle population. Our results provide evidence that polymorphisms in the ADD1 gene are associated with growth traits, and may be used for marker-assisted selection in beef cattle breeding program.     

Validation of genetic polymorphisms on BTA14 associated with carcass trait in a commercial Hanwoo population

The objective of this study was to validate the association of significant SNPs identified from a previous genome‐wide association study with carcass weight (CWT) in a commercial Hanwoo population. We genotyped 13 SNPs located on BTA14 in 867 steers from Korea Hanwoo feedlot bulls. Of these 13 SNPs, five SNPs, namely rs29021868, rs110061498, rs109546980, rs42404006 and rs42303720, were found to be significantly associated (P < 0.001) with CWT. These five significant markers spanned the 24.3 to 29.4 Mb region of BTA14. The most significant marker (rs29021868) for CWT in this study had a 13.07 kg allele substitution effect and accounted for 2.4% of the additive genetic variance in the commercial Hanwoo population. The SNP marker rs109546980 was found to be significantly associated with both CWT (P < 0.001) and eye muscle area (P < 0.001) and could potentially be exploited for marker‐assisted selection in Hanwoo cattle. We also genotyped the ss319607402 variation, which maps to intron2 of PLAG1 gene and which is already reported to be associated with height, to identify any significant association with carcass weight; however, no such association was observed in this Hanwoo commercial population.     

Research on association between variants and haplotypes of TRPV1 and TRPA1 genes with growth traits in three cattle breeds

The objective of this study was to examine the association of transient receptor potential vanilloid 1 (TRPV1) and transient receptor potential ankyrin 1 (TRPA1) genes polymorphisms with growth traits in three Chinese cattle breeds (Jiaxian red cattle, Qinchuan cattle and Luxi cattle). Through experiments we identified three single nucleotide polymorphisms (SNPs) in these three cattle breeds TRPV1 and TRPA1 genes using PCR-SSCP, (forced) PCR-RFLP methods. Three of these polymorphisms are all synonymous mutation which includes (NW_003104493.1: 30327 C?>?T), (NW_003104493.1: 33394 A?>?G) and (NW_003104493.1: 38471?G?>?A) are in exons. The other three polymorphisms are located at 3'UTR. Furthermore, we evaluated the haplotype frequency and the statistical analyses indicated that these SNPs of TRPV1 and TRPA1 genes were associated with bovine body height, body length, waist angle width, hucklebone width, cross ministry height, chest qingwidth (p?<?0.05) and recommendation height, cannon circumference (p?<?0.01) of Qingchuan cattle; body length, waist angle width (p?<?0.05) of Jiaxian red cattle; body weight, Body length, cannon circumference, chest circumference (p?<?0.05) and body height (p?<?0.01) of Luxi cattle. Our result confirms the polymorphisms in the TRPV1 and TRPA1 genes are associated with growth traits that may be used for marker-assisted selection (MAS) in three beef breeds programs.     

PRLH and SOD1 gene variations associated with heat tolerance in Chinese cattle

With the proposed global climate change, heat tolerance is becoming increasingly important to the sustainability of livestock production systems. Results from previous studies showed that variants in the prolactin releasing hormone (PRLH) (AC_000160.1:g.11764610G>A) and superoxide dismutase 1 (SOD1) (AC_000158.1:g.3116044T>A) genes play an important role in heat tolerance in African indicine cattle. However, it is unknown whether or not the mutations are associated with heat tolerance in Chinese cattle. In this study, PCR and DNA sequencing were used to genotype two missense mutations in 725 individuals of 30 cattle breeds. Analysis results demonstrated that two classes of base substitution were detected at two loci: AC_000160.1:g.11764610G>A and AC_000158.1:g.3116044T>A or T>C respectively, with amino acid substitutions arginine to histidine and phenylalanine to isoleucine or leucine. The frequencies of the G and T alleles of the two loci gradually diminished from northern groups to southern groups of native Chinese cattle, whereas the frequencies of A and A or C alleles showed a contrary pattern, displaying a significant geographical difference across native Chinese cattle breeds. Additionally, analysis of these two loci in Chinese indigenous cattle revealed that two SNPs were significantly associated with mean annual temperature (T), relative humidity (RH) and temperature humidity index (THI) (P < 0.01), suggesting that cattle with A or C alleles were distributed in regions with higher T, RH and THI. Our results suggest that the two mutations of PRLH and SOD1 genes in Chinese cattle were associated with the heat tolerance.     

Genetic Variants in SDC3 Gene are Significantly Associated with Growth Traits in Two Chinese Beef Cattle Breeds

Identification of the genes and polymorphisms underlying quantitative traits, and understanding these genes and polymorphisms affect economic growth traits, are important for successful marker-assisted selection and more efficient management strategies in commercial cattle (Bos taurus) population. Syndecan-3 (SDC3), a member of the syndecan family of type I transmembrane heparan sulfate proteoglycans is a novel regulator of feeding behavior and body weight. The aim of this study is to examine the association of the SDC3 polymorphism with growth traits in Chinese Jiaxian and Qinchuan cattle breeds (). Four single nucleotide polymorphisms (SNPs: 1–4) were detected in 555 cows from three Chinese native cattle breeds by means of sequencing pooled DNA samples and polymerase chain reaction-single stranded conformational polymorphism (PCR-SSCP) methods. We found one SNP (g.28362A > G) in intron and three SNPs (g.30742T > G, g.30821C > T and 33418 A > G) in exons. The statistical analyses indicated that these SNPs of SDC3 gene were associated with bovine body height, body length, chest circumference, and circumference of cannon bone (P < 0.05). The mutant-type variant was superior for growth traits; the heterozygote was associated with higher growth traits compared to wild-type homozygote. Our result confirms the polymorphisms in the SDC3 gene are associated with growth traits that may be used for marker-assisted selection in beef cattle breeding programs.     

A genome‐wide association study reveals a quantitative trait locus for days open on chromosome 2 in Japanese Black cattle

Days open (DO), which is the interval from calving to conception, is an important trait related to reproductive performance in cattle. To identify quantitative trait loci for DO in Japanese Black cattle, we conducted a genome‐wide association study with 33 303 single nucleotide polymorphisms (SNPs) using 459 animals with extreme DO values selected from a larger group of 15 488 animals. We identified a SNP on bovine chromosome 2 (BTA2) that was associated with DO. After imputation using phased haplotype data inferred from 586 812 SNPs of 1041 Japanese Black cattle, six SNPs associated with DO were located in an 8.5‐kb region of high linkage disequilibrium on BTA2. These SNPs were located on the telomeric side at a distance of 177 kb from the parathyroid hormone 2 receptor (PTH2R) gene. The association was replicated in a sample of 1778 animals. In the replicated population, the frequency of the reduced‐DO allele (Q) was 0.63, and it accounted for 1.72% of the total genetic variance. The effect of a Q‐to‐q allele substitution on DO was a decrease of 3.74 days. The results suggest that the Q allele could serve as a marker in Japanese Black cattle to select animals with superior DO performance.     

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.