Origins and genetic diversity of New World Creole cattle: inferences from mitochondrial and Y chromosome polymorphisms

The ancestry of New World cattle was investigated through the analysis of mitochondrial and Y chromosome variation in Creoles from Argentina, Brazil, Mexico, Paraguay and the United States of America. Breeds that influenced the Creoles, such as Iberian native, British and Zebu, were also studied. Creoles showed high mtDNA diversity (H = 0.984 ± 0.003) with a total of 78 haplotypes, and the European T3 matriline was the most common (72.1%). The African T1a haplogroup was detected (14.6%), as well as the ancestral African‐derived AA matriline (11.9%), which was absent in the Iberian breeds. Genetic proximity among Creoles, Iberian and Atlantic Islands breeds was inferred through their sharing of mtDNA haplotypes. Y‐haplotype diversity in Creoles was high (H = 0.779 ± 0.019), with several Y1, Y2 and Y3 haplotypes represented. Iberian patrilines in Creoles were more difficult to infer and were reflected by the presence of H3Y1 and H6Y2. Y‐haplotypes confirmed crossbreeding with British cattle, mainly of Hereford with Pampa Chaqueño and Texas Longhorn. Male‐mediated Bos indicus introgression into Creoles was found in all populations, except Argentino1 (herd book registered) and Pampa Chaqueño. The detection of the distinct H22Y3 patriline with the INRA189‐90 allele in Caracú suggests introduction of bulls directly from West Africa. Further studies of Spanish and African breeds are necessary to elucidate the origins of Creole cattle, and determine the exact source of their African lineages.     

我国部分黄牛品种线粒体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%的序列均发生了鸟嘌呤变异, 推测这些变异很可能是我国瘤牛固有的变异。     

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.     

Y-SNPs Haplotype Diversity in Four Chinese Cattle Breeds

To investigate the genetic diversity of Chinese cattle, 96 male samples of 4 Chinese native cattle breeds were investigated using 5 single nucleotide polymorphisms specific to the bovine Y chromosome. Two previously described haplotypes (taurine Y2 and indicine Y3) were detected in 74 and 22 animals, respectively. The haplotype frequencies varied amongst the four native breeds. The taurine Y2 haplotype dominated in the Qinchuan, Dabieshan, and Yunba breeds. However, the indicine Y3 haplotype occurred in high frequency in the Enshi breed. Among the four native breeds, Yunba had the highest haplotype diversity (0.4330 ± 0.0750), followed by Qinchuan (0.2899 ± 0.1028) and Enshi (0.2222 ± 0.1662), Dabieshan was the least differentiated (0.1079 ± 0.0680). Compared with some foreign cattle breeds, the low level of haplotype diversity was detected in our breeds (0.2633 ± 0.1030).     

贵州黄牛mtDNA D-loop 遗传多样性研究

对贵州4个地方黄牛品种共计82个个体的线粒体DNA D-loop区全序列910 bp进行分析,检测到31种单倍型,其核苷酸多态位点65个,约占所测核苷酸总长的7.14%,其中有62个转换,2个颠换,1个转换/颠换共存。贵州4个黄牛品种mtDNA D-loop区核苷酸多样度（π值）为2.16%～2.61%,单倍型多样度（H）为0.695～0.909,表明贵州黄牛mtDNA遗传多样性比较丰富。根据单倍型构建了贵州4个黄牛品种的NJ分子系统树。聚类表明,贵州黄牛有普通牛和瘤牛2大母系起源,其影响较为均一。并探讨了用核苷酸多样度π值的大小来衡量黄牛群体遗传分化程度的可行性。      

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.     

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.     

Dual origins of dairy cattle farming--evidence from a comprehensive survey of European Y-chromosomal variation

Background

Diversity patterns of livestock species are informative to the history of agriculture and indicate uniqueness of breeds as relevant for conservation. So far, most studies on cattle have focused on mitochondrial and autosomal DNA variation. Previous studies of Y-chromosomal variation, with limited breed panels, identified two Bos taurus (taurine) haplogroups (Y1 and Y2; both composed of several haplotypes) and one Bos indicus (indicine/zebu) haplogroup (Y3), as well as a strong phylogeographic structuring of paternal lineages.

Methodology and Principal Findings

Haplogroup data were collected for 2087 animals from 138 breeds. For 111 breeds, these were resolved further by genotyping microsatellites INRA189 (10 alleles) and BM861 (2 alleles). European cattle carry exclusively taurine haplotypes, with the zebu Y-chromosomes having appreciable frequencies in Southwest Asian populations. Y1 is predominant in northern and north-western Europe, but is also observed in several Iberian breeds, as well as in Southwest Asia. A single Y1 haplotype is predominant in north-central Europe and a single Y2 haplotype in central Europe. In contrast, we found both Y1 and Y2 haplotypes in Britain, the Nordic region and Russia, with the highest Y-chromosomal diversity seen in the Iberian Peninsula.

Conclusions

We propose that the homogeneous Y1 and Y2 regions reflect founder effects associated with the development and expansion of two groups of dairy cattle, the pied or red breeds from the North Sea and Baltic coasts and the spotted, yellow or brown breeds from Switzerland, respectively. The present Y1-Y2 contrast in central Europe coincides with historic, linguistic, religious and cultural boundaries.     

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

Male-mediated introgression of Bos indicus genes into Argentine and Bolivian Creole cattle breeds

The geographic distribution and frequency of Bos taurus and Bos indicus Y chromosome haplotypes amongst Argentine and Bolivian Creole cattle breeds were studied, using cytogenetic and molecular genetic techniques. A complete correspondence between Y chromosome morphology and the haplotype of the Y-linked microsatellite marker INRA 124 was found in all males examined. The taurine and indicine haplotypes were detected in 85.7 and 14.3% of the males studied, respectively, although these frequencies varied amongst the different breeds examined. The geographic distribution of this polymorphism suggests a pattern of zebu introgression in South America. The highest frequencies of the Zebu Y-chromosome are found in Brazilian populations (43-90%), in the eastern part of the continent, while it is absent in the southernmost breeds from Uruguay and Argentina. Bolivian breeds, at the centre of the continent, exhibit intermediate values (17-41%). This east/west and north/south gradient of male Zebu introgression could be explained by historical events and environmental factors.     

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

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

Mitochondrial DNA variation in cattle of South China: origin and introgression

Ten restriction endonucleases were used to investigate the mitochondrial DNA restriction fragment length polymorphism (mtDNA RFLP) of 11 native cattle breeds and one cultivated cattle breed in South China. Twenty-three restriction morphs were detected, which can be sorted into five haplotypes. A phylogenetic tree of the haplotypes was constructed by using the 'upgMa' method. Our study showed that haplotype I and II are identical to the zebu (Bos indicus) and taurine (Bos taurus) haplotypes, respectively. Zebu and taurine were the two major origins of cattle populations in South China, and the zebu probably had more influence on the native cattle population than taurine did. Haplotype III is identical to haplotype I of yak (Bos grunniens), which was only detected in the Diqing cattle breed. Haplotype IV was detected for the first time. This haplotype, found only in Dehong cattle, might be from an independent domestication event, probably from another Bos indicus population. Divergence of haplotypes I and IV occurred about 268,000-535,000 years ago, much earlier than the 10,000-year history of cattle husbandry. Our results also suggest a secondary introgession of mtDNA from yak to Diqing cattle.     

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.     

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.     

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.     

中国地方黄牛的Y染色体遗传多样性及其进化起源

中国黄牛的进化起源与遗传多样性一直是国内外动物遗传学家感兴趣的课题之一.本文主要从Y染色体的形态多样性和Y染色体特异性微卫星标记遗传多态性两个方面对中国地方黄牛的遗传多样性和进化起源进行了综述.中国地方黄牛Y染色体具有中着丝粒、亚中着丝粒和近端着丝粒3种类型,这说明中国地方黄牛起源于普通牛和瘤牛.利用Y染色体特异性微卫星标记对中国地方黄牛Y染色体单倍型分布特征及Y染色体基因流模式的分析表明,北方种群中普通牛单倍型频率最高,瘤牛单倍型在南方种群中占优势;在中国不同地域,瘤牛Y染色体单倍型频率呈现自南而北、自东而西逐渐降低的趋势,这再次证实了中国黄牛主要来源于普通牛和瘤牛,这可能是这两类牛群在长期的历史进化过程中,分别从东南方向和西北方向进入我国,并在中原地区汇合的结果.本文为中国地方黄牛品种资源保护和杂交育种工作提供了参考依据.     

Assessment of the nucleotide sequence variability in the bovine T-cell receptor alpha delta joining gene region

The sequence of 2,193 nucleotides from the bovine T-cell receptor alpha/delta joining gene region (TCRADJ) was determined and compared with the corresponding human and murine sequences. The identity was 75.3% for the comparison of the Bos taurus vs. the Homo sapiens sequence and 63.8% for the Bos taurus vs. the Mus musculus sequence. This comparison permitted the identification of the putatively functional elements within the bovine sequence. Direct sequencing of 2,110 nucleotides in nine animals revealed 12 variable sites. Estimates, based on direct sequencing in three Holstein Friesian animals, for the two measures of sequence variability, nucleotide polymorphism (u) and nucleotide diversity (p), were 0.00050 (60.00036) and 0.00077 (60.00056), respectively. The test statistic, Tajima's D, for the comparison of the two measures indicates that the difference between u and p is close to significance (P < 0.05), suggesting the possibility of selective forces acting on the studied genomic region. Allelic variation at 5 of the 12 variable sites was analysed in 359 animals (48 Anatolian Black, 56 Braunvieh, 115 Fleckvieh, 47 Holstein Friesian, 50 Simmental and 43 Pinzgauer) using the oligonucleotide ligation assay (OLA) in combination with the enzyme linked immunoabsorbant assay (ELISA). Nine unambiguous haplotypes could be derived based on animals with a maximum of one heterozygous site. Four to seven haplotypes were present in the different breeds. When taking into account the frequencies of the haplotypes in the different breeds, especially in Anatolian Black, an ancestral cattle population, we could establish the likely phylogenetic relationships of the haplotypes. Such haplotype trees are the basis for cladistic candidate gene analysis. Our study demonstrates that the systematic search of single nucleotide polymorphisms (SNPs) is useful for analysing all aspects of variability of a given genomic region.     

Mitochondrial DNA sequence variation in Portuguese native dog breeds: diversity and phylogenetic affinities

In an extensive survey of the genetic diversity in Portuguese dogs, we have examined an 887-bp fragment of the mitochondrial DNA (mtDNA) from 8 Portuguese, 1 Spanish, and 2 North African native dog breeds, including village dogs from Portugal and Tunisia. Forty-nine haplotypes were found in the 164 individuals analyzed, with private haplotypes being found in several breeds. For example, the Castro Laboreiro Watchdog, a rare breed from a small and isolated region in Portugal, was monomorphic for mtDNA and possessed a new haplotype, which may be provisionally considered a breed-specific marker. Phylogenetic analyses recapitulated 4 major clades identified in other studies, but new haplotypes, grouping within a clade that was previously thought as geographically restricted, were detected in Portugal and Morocco. Portuguese village dogs showed no genetic differentiation from nonnative dogs or from local breeds of the areas in which the village dogs were sampled. Although Iberian and North African dog breeds possessed breed-specific mtDNA haplotypes, no significant geographic structure could be detected among them. There is no evidence for introgression of North African haplotypes in Iberian dogs, contrary to previous results for other domestic animals.     

婆罗门牛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)为独立驯化、起源于印度次大陆的假说。     

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.