Association of bovine β‐casein protein variant I with milk production and milk protein composition

The aim of this study was to detect new polymorphisms in the bovine β‐casein (β‐CN) gene and to evaluate association of (new) β‐CN protein variants with milk production traits and milk protein composition. Screening of the β‐CN gene in genomic DNA from 72 Holstein Friesian (HF) bulls resulted in detection of 19 polymorphisms and revealed the presence of β‐CN protein variant I in the Dutch HF population. Studies of association of β‐CN protein variants with milk composition usually do not discriminate protein variant I from variant A2. Association of β‐CN protein variants with milk composition was studied in 1857 first‐lactation HF cows and showed that associations of protein variants A2 and I were quite different for several traits. β‐CN protein variant I was significantly associated with protein percentage and protein yield, and with α_s1‐casein (α_s1‐CN), α_s2‐casein (α_s2‐CN), κ‐casein (κ‐CN), α‐lactalbumin (α‐LA), β‐lactoglobulin (β‐LG), casein index and casein yield. Inferring β‐κ‐CN haplotypes showed that β‐CN protein variant I occurred only with κ‐CN variant B. Consequently, associations of β‐κ‐CN haplotype IB with protein percentage, κ‐CN, α‐LA, β‐LG and casein index are likely resulting from associations of κ‐CN protein variant B, while associations of β‐κ‐CN haplotype IB with α_s1‐CN and α_s2‐CN seem to be resulting from associations of β‐CN variant I.     

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.     

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

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

A genomic map of climate adaptation in Mediterranean cattle breeds

Domestic species such as cattle (Bos taurus taurus and B. t. indicus) represent attractive biological models to characterize the genetic basis of short‐term evolutionary response to climate pressure induced by their post‐domestication history. Here, using newly generated dense SNP genotyping data, we assessed the structuring of genetic diversity of 21 autochtonous cattle breeds from the whole Mediterranean basin and performed genome‐wide association analyses with covariables discriminating the different Mediterranean climate subtypes. This provided insights into both the demographic and adaptive histories of Mediterranean cattle. In particular, a detailed functional annotation of genes surrounding variants associated with climate variations highlighted several biological functions involved in Mediterranean climate adaptation such as thermotolerance, UV protection, pathogen resistance or metabolism with strong candidate genes identified (e.g., NDUFB3, FBN1, METTL3, LEF1, ANTXR2 and TCF7). Accordingly, our results suggest that main selective pressures affecting cattle in Mediterranean area may have been related to variation in heat and UV exposure, in food resources availability and in exposure to pathogens, such as anthrax bacteria (Bacillus anthracis). Furthermore, the observed contribution of the three main bovine ancestries (indicine, European and African taurine) in these different populations suggested that adaptation to local climate conditions may have either relied on standing genomic variation of taurine origin, or adaptive introgression from indicine origin, depending on the local breed origins. Taken together, our results highlight the genetic uniqueness of local Mediterranean cattle breeds and strongly support conservation of these populations.     

5α‐Androst‐16‐en‐3α‐ol β‐D‐Glucuronide,Precursor of 5α‐Androst‐16‐en‐3α‐ol in Human Sweat

5α‐Androst‐16‐en‐3α‐ol (α‐androstenol) is an important contributor to human axilla sweat odor. It is assumed that α‐andostenol is excreted from the apocrine glands via a H₂O‐soluble conjugate, and this precursor was formally characterized in this study for the first time in human sweat. The possible H₂O‐soluble precursors, sulfate and glucuronide derivatives, were synthesized as analytical standards, i.e., α‐androstenol, β‐androstenol sulfates, 5α‐androsta‐5,16‐dien‐3β‐ol (β‐androstadienol) sulfate, α‐androstenol β‐glucuronide, α‐androstenol α‐glucuronide, β‐androstadienol β‐glucuronide, and α‐androstenol β‐glucuronide furanose. The occurrence of α‐androstenol β‐glucuronide was established by ultra performance liquid chromatography (UPLC)/MS (heated electrospray ionization (HESI)) in negative‐ion mode in pooled human sweat, containing eccrine and apocrine secretions and collected from 25 female and 24 male underarms. Its concentration was of 79 ng/ml in female secretions and 241 ng/ml in male secretions. The release of α‐androstenol was observed after incubation of the sterile human sweat or α‐androstenol β‐glucuronide with a commercial glucuronidase enzyme, the urine‐isolated bacteria Streptococcus agalactiae, and the skin bacteria Staphylococcus warneri DSM 20316, Staphylococcus haemolyticus DSM 20263, and Propionibacterium acnes ATCC 6919, reported to have β‐glucuronidase activities. We demonstrated that if α‐ and β‐androstenols and androstadienol sulfates were present in human sweat, their concentrations would be too low to be considered as potential precursors of malodors; therefore, the H₂O‐soluble precursor of α‐androstenol in apocrine secretion should be a β‐glucuronide.     

Y chromosome diversity and paternal origin of Chinese cattle

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

Composition and Chemical Variability of Ivoirian Polyalthia oliveri Leaf Oil

The chemical composition of 45 essential oil samples isolated from the leaves of Polyalthia oliveri harvested in three Ivoirian forests was investigated by GC‐FID (retention indices measured on two columns of different polarities), and by ¹³C‐NMR, following a method developed in our laboratory. In total, 41 components were identified. The content of the main components varied drastically from sample to sample: (E)‐β‐caryophyllene (1.2 – 50.8%), α‐humulene (0.6 – 47.7%), isoguaiene (0 – 27.9%), alloaromadendrene (0 – 24.7%), germacrene B (0 – 18.3%), δ‐cadinene (0.4 – 19.3%), and β‐selinene (0.2 – 18.5%). The analysis of six oil samples selected in function of their chromatographic profiles is reported in detail. The 45 oil compositions were submitted to hierarchical cluster and principal components analysis, which allowed the distinction of three groups within the oil samples. The compositions of the oils from group I (15 samples) and II (12 samples) were dominated by (E)‐β‐caryophyllene and α‐humulene, respectively. Oil samples of group III (18 samples) needed to be partitioned into four subgroups III.1–III.4 whose compositions were dominated by alloaromadenrene, isoguaiene, germacrene B, and δ‐cadinene, respectively.     

Worldwide Patterns of Ancestry,Divergence, and Admixture in Domesticated Cattle

The domestication and development of cattle has considerably impacted human societies, but the histories of cattle breeds and populations have been poorly understood especially for African, Asian, and American breeds. Using genotypes from 43,043 autosomal single nucleotide polymorphism markers scored in 1,543 animals, we evaluate the population structure of 134 domesticated bovid breeds. Regardless of the analytical method or sample subset, the three major groups of Asian indicine, Eurasian taurine, and African taurine were consistently observed. Patterns of geographic dispersal resulting from co-migration with humans and exportation are recognizable in phylogenetic networks. All analytical methods reveal patterns of hybridization which occurred after divergence. Using 19 breeds, we map the cline of indicine introgression into Africa. We infer that African taurine possess a large portion of wild African auroch ancestry, causing their divergence from Eurasian taurine. We detect exportation patterns in Asia and identify a cline of Eurasian taurine/indicine hybridization in Asia. We also identify the influence of species other than Bos taurus taurus and B. t. indicus in the formation of Asian breeds. We detect the pronounced influence of Shorthorn cattle in the formation of European breeds. Iberian and Italian cattle possess introgression from African taurine. American Criollo cattle originate from Iberia, and not directly from Africa with African ancestry inherited via Iberian ancestors. Indicine introgression into American cattle occurred in the Americas, and not Europe. We argue that cattle migration, movement and trading followed by admixture have been important forces in shaping modern bovine genomic variation.     

Composition and Chemical Variability of Cleistopholis patens Trunk Bark Oil from Côte d'Ivoire

The chemical composition of trunk bark oil from Cleistopholis patens (Benth .) Engl . & Diels , growing wild in Côte d'Ivoire, has been investigated by GC (FID) in combination with retention indices, GC/MS and ¹³C‐NMR. Moreover, one oil sample has been subjected to CC and all the fractions analyzed by GC (RI) and ¹³C‐NMR. In total, 61 components have been identified, including various sesquiterpene esters scarcely found in essential oils. ¹³C‐NMR was particularly efficient for the identification of a component not eluted on GC and for the quantification of heat‐sensitive compounds. Then, 36 oil samples, isolated from trunk bark harvested in six Ivoirian forests have been analyzed. The content of the main components varied drastically from sample to sample: (E)‐β‐caryophyllene (0.4 – 69.1%), β‐pinene (0 – 57%), α‐phellandrene (0 – 33.2%), α‐pinene (0.1 – 30.6%), β‐elemol (0.1 – 29.9%), germacrene D (0 – 25.4%), juvenile hormone III (0 – 22.9%), germacrene B (0 – 20.6%) and sabinene (tr‐20.3%). Statistical analysis, hierarchical clustering and principal components analysis, carried out on the 36 compositions evidenced a fair chemical variability of the stem bark oil of this species. Indeed, three clusters have been distinguished: the composition of group I (ten samples) was dominated by β‐pinene and α‐pinene, group II (nine samples) was represented by α‐phellandrene and p‐cymene and group III (16 samples) by β‐elemol. A sample displayed an atypical composition dominated by (E)‐β‐caryophyllene.     

11α‐Ethoxy‐β‐boswellic Acid and Nizwanone,a New Boswellic Acid Derivative and a New Triterpene,Respectively, from Boswellia sacra

A new boswellic acid derivative, 11α‐ethoxy‐β‐boswellic acid (EBA; 1 ) and a new ursane‐type triterpene, named nizwanone ( 2 ), were isolated from Omani frankincense Boswellia sacra Flueck . together with two known compounds papyriogenin B and rigidenol. The structures of 1 and 2 were elucidated by detailed spectroscopic analysis using ¹H‐ and ¹³C‐NMR, ¹H,¹H‐COSY, HMQC, HMBC, and HR‐EI‐MS techniques. The relative configurations of 1 and 2 were assigned by comparative analysis of the NMR spectral data with those of known analogs together with NOESY experiments. Structures of known compounds were identified by comparison with the reported data.     

Polymorphisms of two Y chromosome microsatellites in Chinese cattle

Two Y chromosome specific microsatellites UMN2404 and UMN0103 were genotyped and assessed for polymorphisms in a total of 423 unrelated males from 25 indigenous Chinese cattle breeds. Consistently, both microsatellites displayed specific indicine and taurine alleles in each bull examined. The indicine and taurine alleles were detected in 248 males (58.6%), and 175 males (41.4%), respectively, although these frequencies varied amongst different breeds examined. The indicine alleles dominated in the southern group (92.4%), while the taurine alleles dominated in the northern group (95.5%). Hainan Island was possibly the site for the origin of Chinese zebu, and Tibetan cattle were probably independently domesticated from another strain of Bos primigenius. The geographical distribution of these frequencies reveals a pattern of male indicine introgression and a hybrid zone of indicine and taurine cattle in China. The declining south-to-north and east-to-west gradient of male indicine introgression in China could be explained by historical data, geographical segregation and temperature and weather conditions.     

Genome‐wide detection of signatures of selection in Korean Hanwoo cattle

The Korean Hanwoo cattle have been intensively selected for production traits, especially high intramuscular fat content. It is believed that ancient crossings between different breeds contributed to forming the Hanwoo, but little is known about the genomic differences and similarities between other cattle breeds and the Hanwoo. In this work, cattle breeds were grouped by origin into four types and used for comparisons: the Europeans (represented by six breeds), zebu (Nelore), African taurine (N'Dama) and Hanwoo. All animals had genotypes for around 680 000 SNPs after quality control of genotypes. Average heterozygosity was lower in Nelore and N'Dama (0.22 and 0.21 respectively) than in Europeans (0.26–0.31, with Shorthorn as outlier at 0.24) and Hanwoo (0.29). Pairwise F_ST analyses demonstrated that Hanwoo are more related to European cattle than to Nelore, with N'Dama in an intermediate position. This finding was corroborated by principal components and unsupervised hierarchical clustering. Using genome‐wide smoothed F_ST, 55 genomic regions potentially under positive selection in Hanwoo were identified. Among these, 29 were regions also detected in previous studies. Twenty‐four regions were exclusive to Hanwoo, and a number of other regions were shared with one or two of the other groups. These regions overlap a number of genes that are related to immune, reproduction and fatty acid metabolism pathways. Further analyses are needed to better characterize the ancestry of the Hanwoo cattle and to define the genes responsible to the identified selection peaks.     

Genetic diversity and genomic signatures of selection among cattle breeds from Siberia,eastern and northern Europe

Domestication in the near eastern region had a major impact on the gene pool of humpless taurine cattle (Bos taurus). As a result of subsequent natural and artificial selection, hundreds of different breeds have evolved, displaying a broad range of phenotypic traits. Here, 10 Eurasian B. taurus breeds from different biogeographic and production conditions, which exhibit different demographic histories and have been under artificial selection at various intensities, were investigated using the Illumina BovineSNP50 panel to understand their genetic diversity and population structure. In addition, we scanned genomes from eight breeds for signatures of diversifying selection. Our population structure analysis indicated six distinct breed groups, the most divergent being the Yakutian cattle from Siberia. Selection signals were shared (experimental P‐value < 0.01) with more than four breeds on chromosomes 6, 7, 13, 16 and 22. The strongest selection signals in the Yakutian cattle were found on chromosomes 7 and 21, where a miRNA gene and genes related to immune system processes are respectively located. In general, genomic regions indicating selection overlapped with known QTL associated with milk production (e.g. on chromosome 19), reproduction (e.g. on chromosome 24) and meat quality (e.g. on chromosome 7). The selection map created in this study shows that native cattle breeds and their genetic resources represent unique material for future breeding.     

New Alkaloids and α‐Glucosidase Inhibitory Flavonoids from Ficus hispida

Two new pyrrolidine alkaloids, ficushispimines A ( 1 ) and B ( 2 ), a new ω‐(dimethylamino)caprophenone alkaloid, ficushispimine C ( 3 ), and a new indolizidine alkaloid, ficushispidine ( 4 ), together with the known alkaloid 5 and 11 known isoprenylated flavonoids 6  –  16 , were isolated from the twigs of Ficus hispida. Their structures were elucidated by spectroscopic methods. Isoderrone ( 8 ), 3′‐(3‐methylbut‐2‐en‐1‐yl)biochanin A ( 11 ), myrsininone A ( 12 ), ficusin A ( 13 ), and 4′,5,7‐trihydroxy‐6‐[(1R*,6R*)‐3‐methyl‐6‐(1‐methylethenyl)cyclohex‐2‐en‐1‐yl]isoflavone ( 14 ) showed inhibitory effects on α‐glucosidase in vitro.     

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.     

miR‐29a modulates tumor necrosis factor‐α‐induced osteogenic inhibition by targeting Wnt antagonists

We previously found that miR‐29a was significantly downregulated in Ankylosing spondylitis (AS) patients, a chronic inflammatory disease associated with bone metabolic disorder, however, the underlying mechanism remains unclear. In this study, we demonstrated that miR‐29a regulates tumor necrosis factor‐α (TNF‐α) mediated bone loss mainly by targeting DKK1 and GSK3β, thus activating the Wnt/β‐catenin pathway. Our findings may provide new insight into the pathogenesis of the bone metabolism disorder in inflammation environment and provide promising therapeutic target.     

Chemotaxonomic Importance of the Essential‐Oil Composition in Two Subspecies of Teucrium stocksianum Boiss. from Iran

The hydrodistilled essential oils obtained from aerial flowering parts of Teucrium stocksianum ssp. stocksianum (TSS) and T. stocksianum ssp. gabrielae (TSG) from Iran were analyzed by capillary GC and GC/MS. The oil analysis of two subspecies led to the identification of 65 compounds that accounted for 93.3 and 95.1% of the total oil compositions, respectively. Sesquiterpenoids (52.9%) constituted the main compounds in the essential oil of TSS represented mainly by cis‐sesquisabinene hydrate (12.0%), followed by epi‐β‐bisabolol (6.6%), guaiol (5.4%), and β‐eudesmol (4.4%), whilst monoterpenoids (61.2%) were found to be the major components of the oil of TSG, represented by α‐pinene (23.0%), β‐pinene (13.0%), myrcene (6.3%), and sabinene (6.3%). The principal component in both subspecies, TSS and TSG, was α‐pinene (22.0 and 23.0%, resp.) and β‐pinene (6.5 and 13.0%, resp.). epi‐α‐Cadinol, myrcene, and sabinene, which were detected as principal compounds of TSG, were characterized in lower amounts (<1.5%) in the oil of TSS. Seven components were identified in the oil of TSS corresponding to 25.9% of total oil, which were totally absent in the oil of TSG, of which cis‐sesquisabinene hydrate (12.0%), guaiol (5.4%), and β‐eudesmol (4.4%) were in considerable amounts. Taxonomic position of the subspecies is discussed on the basis of phytochemical data.     

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.     

Paternal origins of Chinese cattle

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