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.     

我国部分黄牛品种线粒体D-loop区遗传多样性与起源分化

为了解我国地方黄牛品种线粒体DNA的遗传变异情况, 文章测定了16个地方黄牛品种206个个体线粒体D-loop区的全序列, 共检测到101个变异位点; 99种单倍型, 其中73种是普通牛单倍型, 26种是瘤牛单倍型; 平均核苷酸差异为22.6920, 单倍型多样度为0.9320, 核苷酸多样度为0.0227, 表明我国黄牛品种遗传多样性非常丰富。构建的NJ进化树显示16个品种来源于两大母系: 普通牛和瘤牛; 构建的Network图表明73种普通牛单倍型可以分为3大单倍型群; 26种瘤牛单倍型分为5种单倍型群, 推测我国瘤牛在迁移过程中, 至少经历了4次群体扩张事件。通过分析比较地方黄牛品种与内罗门牛共有的 H3单倍型, 发现其中只有16%的序列与内罗门牛的H3单倍型非常相似, 其余84%的序列均发生了鸟嘌呤变异, 推测这些变异很可能是我国瘤牛固有的变异。     

Evidence for a geographical cline of casein haplotypes in Portuguese cattle breeds

Genetic variants of bovine milk proteins have been intensively used to characterize breeds and as markers for population/QTL studies throughout the world. However, a large number of cattle breeds including those found in Portugal, remain unstudied. In this work, we have analysed the genetic variation of six milk protein loci in 10 Portuguese cattle breeds by isoelectric focusing. High genetic diversities were generally found across breeds, with the exception of Mirandesa that showed a trend to fixation of the most common alleles in five loci, as well as of the rarer CSN3B allele. The casein haplotype BA2A was often the most frequent, followed by haplotypes BA2B and BA1A. Remarkably, CA2A was found to be the second most frequent haplotype in Southern breeds, supporting a geographical cline between Central-Northern European breeds and Bos indicus populations. Our data suggest that high genetic similarity among neighbouring Portuguese breeds is mainly caused by gene flow, and that the geographical distribution of particular casein haplotypes may indicate an influence of African cattle.     

Independent mitochondrial origin and historical genetic differentiation in North Eastern Asian cattle

In order to clarify the origin and genetic diversity of cattle in North Eastern Asia, this study examined mitochondrial displacement loop sequence variation and frequencies of Bos taurus and Bos indicus Y chromosome haplotypes in Japanese, Mongolian, and Korean native cattle. In mitochondrial analyses, 20% of Mongolian cattle carried B. indicus mitochondrial haplotypes, but Japanese and Korean cattle carried only B. taurus haplotypes. In contrast, all samples revealed B. taurus Y chromosome haplotypes. This may be due to the import of zebu and other cattle during the Mongol Empire era with subsequent crossing with native taurine cattle. B. taurus mtDNA sequences fall into several geographically distributed haplogroups and one of these, termed here T4, is described in each of the test samples, but has not been observed in Near Eastern, European or African cattle. This may have been locally domesticated from an East Eurasian strain of Bos primigenius.     

婆罗门牛mtDNA D-loop变异及其遗传背景

对10头原种婆罗门牛mtDNAD-loop全序列912 bp测序,婆罗门牛遗传多样性丰富,检测到的9种单倍型兼有瘤牛(B.indicus)与普通牛(B.taurus)的遗传背景,核苷酸变异率为6.25 %,单倍型多态度为0.978±0.054 ,核苷酸多态度为0.014 30±0.008 68。所有单倍型聚为明显的两大分支,婆罗门牛的大部分单倍型为普通牛单倍型类群,并占绝对优势(90 %) ,仅Brah-6与亚洲瘤牛聚在一起,属于亚洲瘤牛线粒体单倍型,表明婆罗门牛的确是集亚洲瘤牛、欧洲普通牛等优良特性于一身(易产犊、产肉性能好、耐热与体表寄生虫等)的瘤牛品种之一。育种学家引种瘤牛的目的是改善当地牛的生产力与适应性,现代普通牛表现出明显又普遍的瘤牛渐渗现象。对现代的瘤牛品种而言,除亚洲瘤牛品种外,普通牛对其他瘤牛品种育成的贡献同样高。支持瘤牛(B.indicus)为独立驯化、起源于印度次大陆的假说。     

Admixture and diversity in West African cattle populations

We present a population genetic analysis of microsatellite variation in 16 West African cattle populations. West Africa represents a unique juxtaposition of different climatic and ecological zones in a relatively small geographical area. While more humid coastal regions are inhabited by the tsetse fly, a vector which spreads trypanosomiasis among cattle, the disease is not transmitted in the drier areas outside this zone. This is the most thorough study of genetic diversity in cattle within this area, which contains genetically important trypanotolerant Bos taurus breeds. Genetic relationships among the many breeds are examined and levels of diversity are assessed. Admixture levels were determined using a variety of methods. Ancestry informative or population-associated alleles (PAAs) were selected using populations from India, the Near East and Europe. Multivariate analysis, the admix program and model-based Bayesian admixture analysis approaches were also employed. These analyses reveal the direct impact of ecological factors and the profound effect of admixture on the cattle of this region. They also highlight the importance of efforts to prevent further dilution of African taurine breeds by B. indicus cattle.     

中国黄牛mtDNA D-loop遗传多样性及起源

黄牛自古以来就是我国一个重要的畜种,其经济、文化价值很高。我国是世界上黄牛品种资源最丰富的国家之一。据《中国牛品种志》介绍,把一些地区同种异名的黄牛品种合并以后,尚有28个地方黄牛品种,按照其地理分布区域分为北方黄牛、中原黄牛和南方黄牛三大类型(邱怀,1986)。如果把中国地方黄牛品种分得更细,则有49个固有品种(常洪,1995)。关于中国黄牛的起源,历来有不同的观点。一般认为,中国黄牛是多元起源的,但究竟起源于哪几个牛种,观点不一(陈宏等,1993;于汝梁等,1993;Yu et al.,1999;陈幼春,1990)。主要的观点有:(1)中国黄牛主要起源于…     

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.     

Polymorphisms in Blood Proteins of Bos indicus and Bos taurus Cattle Breeds of Cameroon and Nigeria, and Description of New Albumin Variants

Polymorphisms in the five blood protein loci albumin (ALB), carbonic anhydrase (CA II), vitamin D binding protein (GC), haemoglobin (HBB), and transferrin (TF) were investigated in 520 individuals from 12 cattle populations (Bos indicus and Bos taurus) in Cameroon and Nigeria by isoelectric focusing with carrier ampholytes in ultrathin polyacrylamide gels (PAG-IEF) and by linear gradient polyacrylamide gel electrophoresis (PAGE). While all loci in nine populations were polymorphic with up to six alleles at the ALB and TF loci: the Namchi population showed monomorphism at the CA II locus and Muturu at the ALB, CA II, and HBB loci. There was a clear distinction between Bos indicus and Bos taurus breeds at the ALB locus with ALBB predominating in indicine and ALBA predominating in taurtine breeds. CA IIS, GCA, and HBBA were the most commonly occurring alleles in all populations. Two variants not described before were demonstrated by PAG-IEF at the ALB locus and named ALBJ and ALBK. Mean effective number of alleles as measure of intrabreed diversity was higher in zebu populations (2.040-2.288) as compared to taurine breeds (1.349-1.836). Significant deviations from Hardy-Weinberg equilibrium occurred in some populations at the HBB and TF loci. More haplotypes of ALB/GC occurred in the zebu than taurine breeds. ALBAGCA predominated in the taurine populations and ALBBGCA in the indicine populations. Influence of zebu genes on the Namchi and N'Dama taurine breeds was detected at the ALB, CA II, HBB, and TF loci, and estimated at 61.5% and 5.7%, respectively. The high resolution of PAG-IEF in screening for polymorphisms within diversity studies was demonstrated.     

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.     

Extensive intragenic recombination and patterns of linkage disequilibrium at the CSN3 locus in European rabbit

Kappa-casein (CSN3) plays an important role in stabilising the Ca-sensitive caseins in the micelle. The European rabbit (Oryctolagus cuniculus) CSN3 has previously been shown to possess two alleles (A and B), which differ deeply in their intronic regions (indels of 100 and 1550 nucleotides in introns 1 and 4, respectively). Furthermore, a correlation between several reproductive performance traits and the different alleles was described. However, all these data were exclusively collected in rabbit domestic breeds, preventing a deeper understanding of the extensive polymorphism observed in the CSN3 gene. Additionally, the techniques available for the typing of both indel polymorphisms were until now not suitable for large-scale studies. In this report, we describe a simple, PCR-based typing method to distinguish rabbit CSN3 alleles. We analyse both ancient wild rabbit populations from the Iberian Peninsula and France, and the more recently derived English wild rabbits and domestic stocks. A new allele (C) showing another major indel (250 bp) in intron 1 was found, but exclusively detected in Iberian wild rabbits. In addition, our survey revealed the occurrence of new haplotypes in wild populations, suggesting that intragenic recombination is important in creating genetic diversity at this locus. This easy and low cost single-step PCR-based method results in an improvement over previous described techniques, can be easily set up in a routine molecular laboratory and would probably be a valuable tool in the management of rabbit domestic breeds.     

Y‐specific microsatellites reveal an African subfamily in taurine (Bos taurus) cattle

Five cattle Y‐specific microsatellites, totalling six loci, were selected from a set of 44 markers and genotyped on 608 Bos taurus males belonging to 45 cattle populations from Europe and Africa. A total of 38 haplotypes were identified. Haplogroups (Y1 and Y2) previously defined using single nucleotide polymorphisms did not share haplotypes. Nine of the 27 Y2‐haplotypes were only present in African cattle. Network and correspondence analyses showed that this African‐specific subfamily clustered separately from the main Y2‐subfamily and the Y1 haplotypes. Within‐breed genetic variability was generally low, with most breeds (78%) showing haplotypes belonging to a single haplogroup. amova analysis showed that partitioning of genetic variation among breeds can be mainly explained by their geographical and haplogroup assignment. Between‐breed genetic variability summarized via Principal Component Analysis allowed the identification of three principal components explaining 94.2% of the available information. Projection of principal components on geographical maps illustrated that cattle populations located in mainland Europe, the three European Peninsulas and Mediterranean Africa presented similar genetic variation, whereas those breeds from Atlantic Europe and British Islands (mainly carrying Y1 haplotypes) and those from Sub‐Saharan Africa (belonging to Y2‐haplogroup) showed genetic variation of a different origin. Our study confirmed the existence of two large Y‐chromosome lineages (Y1 and Y2) in taurine cattle. However, Y‐specific microsatellites increased analytical resolution and allowed at least two different Y2‐haplotypic subfamilies to be distinguished, one of them restricted to the African continent.     

Cattle demographic history modelled from autosomal sequence variation

The phylogeography of cattle genetic variants has been extensively described and has informed the history of domestication. However, there remains a dearth of demographic models inferred from such data. Here, we describe sequence diversity at 37 000 bp sampled from 17 genes in cattle from Africa, Europe and India. Clearly distinct population histories are suggested between Bos indicus and Bos taurus, with the former displaying higher diversity statistics. We compare the unfolded site frequency spectra in each to those simulated using a diffusion approximation method and build a best-fitting model of past demography. This implies an earlier, possibly glaciation-induced population bottleneck in B. taurus ancestry with a later, possibly domestication-associated demographic constriction in B. indicus. Strikingly, the modelled indicine history also requires a majority secondary admixture from the South Asian aurochs, indicating a complex, more diffuse domestication process. This perhaps involved multiple domestications and/or introgression from wild oxen to domestic herds; the latter is plausible from archaeological evidence of contemporaneous wild and domestic remains across different regions of South Asia.     

中国部分地方黄牛品种线粒体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,并且对它进行了讨论。同时证明了西藏瘤牛独立于中国瘤牛核心类群的特殊性。     

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.     

Detection of copy number variations and their effects in Chinese bulls

Background

Copy number variations (CNVs) are a main source of genomic structural variations underlying animal evolution and production traits. Here, with one pure-blooded Angus bull as reference, we describe a genome-wide analysis of CNVs based on comparative genomic hybridization arrays in 29 Chinese domesticated bulls and examined their effects on gene expression and cattle growth traits.

Results

We identified 486 copy number variable regions (CNVRs), covering 2.45% of the bovine genome, in 24 taurine (Bos taurus), together with 161 ones in 2 yaks (Bos grunniens) and 163 ones in 3 buffaloes (Bubalus bubalis). Totally, we discovered 605 integrated CNVRs, with more “loss” events than both “gain” and “both” ones, and clearly clustered them into three cattle groups. Interestingly, we confirmed their uneven distributions across chromosomes, and the differences of mitochondrion DNA copy number (gain: taurine, loss: yak & buffalo). Furthermore, we confirmed approximately 41.8% (253/605) and 70.6% (427/605) CNVRs span cattle genes and quantitative trait loci (QTLs), respectively. Finally, we confirmed 6 CNVRs in 9 chosen ones by using quantitative PCR, and further demonstrated that CNVR22 had significantly negative effects on expression of PLA2G2D gene, and both CNVR22 and CNVR310 were associated with body measurements in Chinese cattle, suggesting their key effects on gene expression and cattle traits.

Conclusions

The results advanced our understanding of CNV as an important genomic structural variation in taurine, yak and buffalo. This study provides a highly valuable resource for Chinese cattle’s evolution and breeding researches.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-480) contains supplementary material, which is available to authorized users.     

Characterization and genetic analysis of bovine αs1-casein I variant

The aim of this study was to identify the molecular genetic origin underlying the I variant of α_s1-casein and to develop a DNA-based test for this polymorphism as a tool for genetic analyses independent of milk sample testing. All coding exons and flanking regions of the α _{s 1} -casein gene were sequenced in DNA samples from cattle of known α _{s 1} -casein genotypes ( BI , CI , II , CC ), determined by isoelectric focusing of milk samples. A nucleotide substitution (A>T) in exon 11 (g.19836A>T) leads to the exchange of Glu with Asp at amino acid position 84 of the mature protein (p.Glu84Asp) and perfectly co-segregated with the presence of the α_s1-casein I variant in the milk of the analysed animals. Genotyping of a total of 680 DNA samples from 31 Bos taurus and Bos indicus cattle breeds and from Bos grunniens , Bison bison and Bison bonasus by restriction fragment length polymorphism analysis revealed the occurrence of Asp at position 84 at low frequencies in Bos taurus and Bos indicus breeds and established its origin from the α _s1-casein C variant (p.Glu192Gly). Ten different intragenic haplotypes in the gene region from intron 8 to intron 12 were observed by sequencing, of which two occurred in Bison bison and one in Bison bonasus only. Using available casein gene complex information, an association of Asp at position 84 to β-casein A ² and κ-casein B was shown in the Bos indicus breed Banyo Gudali . Taken together, we can postulate that the α _s1-casein variant I is caused by a non-synonymous nucleotide substitution in exon 11 of the gene and that it originated within Bos indicus and spread to Bos taurus subsequently.     

The Origins of Domesticated Cattle

The origin of domestic cattle has perplexed archaeologists for more than a century. Researchers have proposed various theories, which offer alternative spatial and chronological models for the origin and spread of domesticated cattle. One point of discussion is whether domestic cattle had a single or multiple origins; however, most authorities considered that the first steps towards cattle domestication were taken in southwest Asia and that domesticated cattle entered Europe with pastoralists migrating from this region. Domesticated taurine cattle were thought to have entered Africa in successive waves from southwest Asia, while zebu cattle migrated into Africa at a later date from Arabia and the Indian subcontinent. Analysis of mitochondrial DNA (mtDNA) shows that taurine and zebu cattle divergence before the Holocene and were probably domesticated independently. Recent mtDNA sequence data shows that African and European taurine cattle were probably domesticated independently, but that there was a process of genetic introgression between taurine and zebu cattle in Africa. Ancient DNA studies over the last 10 years suggest that Northern European aurochsens apparently contributed little or nothing to domestic cattle while Southern European aurochsens apparently made a significant input. However, Middle Eastern aurochsen, unfortunately not typed yet, are expected to be to be very similar to European breeds as well, both because archeological data suggest that the major center of domestication for European Bos taurus breeds was the Fertile Crescent (9), and also because a mtDNA sequence from a Syrian specimen dated at 8,000–9,000 years ago shows a typical European haplotype found both in modern breeds and the Italian aurochsen. Evidence seems to suggest that small to moderate levels of local gene flow from wild Bos primigenius females in selected breeds were either accepted or may be reinforced by Neolithic breeders.     

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.