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

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

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

微卫星标记对12个中外牛品种群体遗传结构的研究

选用联合国粮农组织(FAO)和国际动物遗传学会(ISAG)推荐的12对微卫星引物, 采用荧光标记–多重PCR技术, 检测了9个中国地方黄牛品种和3个外来牛品种的遗传多样性。利用等位基因频率计算出各群体的平均遗传杂合度(H)、多态信息含量(PIC)和群体间的DA及DS遗传距离。基于DA遗传距离, 用UPGMA法进行聚类分析, 结果12个中外牛品种被聚为4类: Ⅰ类属于南方黄牛品种, 包括恩施牛、黎平牛、昭通牛和川南山地牛; Ⅱ类属于中原黄牛品种, 包括郏县红牛、早胜牛和平陆山地牛; Ⅲ类属于北方黄牛, 包括延边牛和长白地方牛; Ⅳ类属于外来牛品种, 包括西门塔尔牛、夏洛来牛和德国黄牛。研究结果为中国地方牛品种的保护和利用提供了理论依据。     

采用微卫星标记分析13个中外牛品种的遗传变异和品种间的遗传关系

本研究应用联合国粮农组织(FAO)和国际动物遗传学会(ISAG)推荐的10对微卫星引物,结合荧光–多重PCR技术,检测了10个中国地方黄牛品种和3个外来牛品种的基因型。通过计算基因频率、多态信息含量和遗传杂合度,以Nei’s遗传距离和Nei’s标准遗传距离为基础,采用非加权组对算术平均聚类法构建了聚类图,分析了13个牛品种的群体内遗传变异和群体间遗传关系。并以聚类分析和群体结构分析为基础,将13个中外黄牛品种分为三类:Ⅰ类属于普通黄牛品种,包括延边牛、沿江牛、长白地方牛、蒙古牛、阿勒泰白头牛、哈萨克牛、复州牛和西藏牛;Ⅱ类属于含有瘤牛血统的黄牛品种,包括日喀则驼峰牛和阿沛甲咂牛;Ⅲ类属于外来牛品种,包括德国黄牛、西门塔尔牛和夏洛来牛。研究结果为加强我国地方黄牛品种种质特性研究以及地方牛品种资源的保护与利用提供了科学的依据。     

中国四品种黄牛性染色体多态性的研究

对4个地方黄牛品种（蒙古牛,秦川牛,岭南牛和西镇牛）共63头牛（37♂,26♀）染色体核型的比较研究表明,性染色体在品种间和品种内均存在多态性。Y染色体有中、亚中和近端着丝点染色体,其频率有一定的分布规律。X染色体均为一大的亚中着丝点染色体,但臂比和着丝点指数有品种差异。本文从性染色体多态性上揭示了中国黄牛起源的多元性,为中国黄牛品种的系统分类增添了又一新的依据和线索。     

14.从Cyt b基因全序列分析中国10个黄牛品种的系统进化关系

基于Cyt b基因全长序列,对我国10个地方黄牛品种的遗传多态性和系统发育关系进行了分析.中国黄牛Cyt b基因序列富含碱基A和T,A+T平均含量为56.7%,密码子第3位点核苷酸表现出较强的碱基组成偏倚,该位点上G的含量最低,为 3.5%,而A的含量最高,为44.6%,核苷酸组成偏差指数为0.45.在47个个体1 140 bp的序列中检测到30个变异位点,约占核苷酸总数的2.63%,密码子第3位点的突变率(76.7%)远高于第1位点(16.7%)和第2位点(6.6%),第2位点的非同义替代率(100%)远高于第1位点(40%)和第3位点(4.3%).定义了15种单倍型.将10个品种作为一个大的种群分析时,发现中国黄牛Cyt b基因具有较高的单倍型多样性(Hd=0.783±0046)和一定程度的核苷酸多样性(Pi=0.00828±0.00030),这与中国黄牛地理分布广、种群庞大和起源复杂有关系.基于Kimura 2 parameter距离采用UPGMA法构建的10个黄牛品种分子系统树中将47个个体分别与瘤牛和普通牛聚在一起,说明中国黄牛主要起源于瘤牛和普通牛.北方牛受普通牛的影响较大,南方牛受瘤牛的影响较大,但是对于亲缘关系复杂的中原地区黄牛而言,瘤牛和普通牛对它们的影响力大小因品种而异,因而它们分别和北方牛和南方牛之间的亲缘关系远近程度也存在差异.不同品种之间的亲缘关系差异可能与它们的地理分布差异具有一定的关系.     

中国地方黄牛的遗传多样性*

中国黄牛起源复杂,我国地方黄牛群体不同品种在毛色、形态外貌、细胞遗传学、血液蛋白座位分析均表现出多样性。计算我国黄牛群体6个毛色座位平均杂合度和6个血液蛋白座位平均杂合度分别为0.3144和0.4873,表明我国地方黄牛群体的遗传多样性非常丰富。计算我国黄牛群体的6个毛色座位和6个血液蛋白座位的基因分化系数分别为0.3404和0.095,表明我国黄牛群体毛色差异中有34.04％是由品种间的差异造成的。血液蛋白的多态性有9.5％是由品种间的差异造成的。我国黄牛群体的遗传多样性主要来自品种内的遗传多样性。保存我国黄牛品种资源多样性不仅要从整个中国黄牛群体上考虑,而且要针对每个品种（或类群）进行保种。     

10个牛品种线粒体12S rRNA基因多态性分析

以鲁西牛、渤海黑牛、大别山牛、南阳牛、晋南牛、秦川牛、中国西门塔尔牛、日本和牛、海福特牛和安格斯牛10个牛品种共计353头牛为研究对象,利用单链构象多态性（SSCP）分子标记方法对其线粒体DNA 12S rRNA基因进行了多态性分析,结果发现该基因出现分布频率不同的3种单倍型（Ⅰ、Ⅱ、Ⅲ）,其中单倍型Ⅰ出现的频率最低,仅在晋南牛1个个体中发现;单倍型Ⅱ在中国本地黄牛中都存在,在中国西门塔尔牛中的分布频率很低,但在3个国外品种海福特牛、安格斯牛和日本和牛中不存在;单倍型Ⅲ在这10个牛品种中均有分布,但以国外品种最丰富。计算平均多态信息含量,发现我国地方品种在0.232～0.423之间,基本属于中度多态,说明我国地方黄牛品种在12S rRNA上多态性比较丰富,国外品种较为贫乏;但是单倍型Ⅰ频率很低有丢失的趋势,需要对其进行保护。     

中国部分黄牛品种mtDNA遗传多态性研究

对我国8个黄牛品种22个个体的mtDNA D-loop区910bp全序列进行了分析。结果表明：8个黄牛品种D-loop区序列中,A T平均含量为61．65％;经比对,共检测到66个核苷酸多态位点,约占核苷酸总数的7．25％;D-loop全序列突变类型有5种,即转换、颠换、插入、缺失及转换与颠换共存,它们分别占核苷酸多态位点的81．82％、6．06％、7．57％、3．03％及1．52％。以欧洲牛mtDNA D-loop全序列为标准,8个黄牛群体D-loop的平均核苷酸变异率分3个层次：西镇牛、蒙古牛、黑白花牛及秦川牛的核苷酸变异率最低,分别为0．37％、0．44％、0．52％和0．66％;南阳牛与郏县红牛的核苷酸变异率居中,分别为1．91％和2．02％;晋南牛与岳阳牛的核苷酸变异率最高,分别为4．47％和4．73％。中国黄牛品种内D-loop区序列歧异度为0．55％～5．39％,品种间序列歧异度为1．21％～6．59％。在所测黄牛个体中,mtDNA D-loop序列由19种单倍型组成,单倍型比例为86．36%,说明中国黄牛mtDNA遗传多态性很丰富。由此构建了中国8个黄牛品种的NJ分子系统树,聚类分析表明：所测黄牛的mtDNA D-loop序列表现为3个单倍型组,从而揭示中国黄牛可能有3个母系起源,以普通牛起源和瘤牛起源为主。     

中国17个黄牛品种mtDNA变异特征与多态性分析

中国黄牛品种资源丰富,尚有28个地方固有品种.为了进一步深入了解这些宝贵遗传资源,本研究通过mtDNA变异特征与多态性分析揭示这些来自中国不同地域地方黄牛的母系起源与分子系统学特征.在17个品种84个体的mtDNA D-loop全序列中,一共检测到了102个核苷酸替代突变位.由此定义的53个单倍型被类聚为2个明显的单倍群：普通牛和瘤牛.mtDNA D-loop全序列变异的第一个特征是转换发生的频率远高于颠换;第二个特征是缺失与替代突变共存;第三个特征是缺失突变率比较高.所有D-loop全序列的核苷酸多样性和单倍型多样性分别为0.026 78±0.000 50和0.919±0.027.普通牛D-loop单倍型在北方牛种群中占有优势(80%～100%),而瘤牛单倍型在南方牛种群中占有优势(42.9%～100%),2种不同单倍型在中原牛种群中的分布也存在差异.2种不同单倍型在中国不同地域17个黄牛品种中的差异性分布揭示出了瘤牛mtDNA基因在中国黄牛中自南而北、由高到低的流动模式,这种基因流动模式的形成可能是由历史事件、地理隔离以及气候环境差异等造成的.     

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

Genetic diversity and origin of Chinese cattle revealed by mtDNA D-loop sequence variation

To determine the origin and genetic diversity of Chinese cattle, we analyzed the complete mtDNA D-loop sequences of 84 cattle from 14 breeds/populations from southwest and west China, together with the available cattle sequences in GenBank. Our results showed that the Chinese cattle samples converged into two main groups, which correspond to the two species Bos taurus and Bos indicus. Although a dominant lineage was clearly discerned in both B. taurus and B. indicus mtDNAs, network analysis of the lineages in each of the two species further revealed multiple clades that presented regional difference. The B. taurus samples in China could be grouped into clades T2, T3, and T4, whereas B. indicus harbored two clades I1 and I2. Age estimation of these discerned clades showed a time range of 14,100-44,500 years before present (YBP). The phylogenetic pattern of Chinese cattle was consistent with the recently described cattle matrilineal pool from northeast Asia, but suggested that B. indicus contributed more to the cattle from south and southwest China. The genetic diversity of Chinese cattle varied among the breeds studied.     

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.     

Origin and phylogeographical structure of Chinese cattle

Complete mitochondrial D-loop sequences of 231 samples were used to explore the origin and genetic diversity of Chinese cattle. Phylogenetical analysis of these sequences revealed both Bos taurus and Bos indicus mitochondrial types in Chinese cattle. Four of the previously identified mitochondrial DNA lineages (T1–T4) were identified in the Bos taurus type, including lineage T1, which was found for the first time in Chinese cattle. Two lineages (I1 and I2) were identified in the Bos indicus type. Our results support the suggestion that the Yunnan-Guizhou Plateau is the domestication site of Chinese zebu. We also found evidence that Tibetan cattle originated from taurine and zebu cattle. The distribution pattern of Chinese cattle breeds was closely related to the geographical and climatic background. It was possible to divide Chinese cattle in this study into two major groups: northern and southern cattle.     

凉地方牛群体母系遗传背景研究

黄牛（Bos taurus）的驯养作为农耕文化重要内容之一,在我国历史悠久,长久以来形成的役用型逐步被肉用型所替代。为了揭示平凉地方牛群体的遗传背景,分析其是否具有生产优质牛肉的遗传基础,本研究测定了88头平凉地方牛群体mtDNA D-loop HVS序列,并对包括平凉地方牛群体在内的我国23个地方牛群体单倍型分布及系统发生关系进行了分析。结果表明,mtDNA D-loop 高变区,在平凉地方牛群体共有52个单倍型,而23个地方牛群体共有95个单倍型,这些单倍型在系统发生树和中介网络关系中分布于两个分枝,即瘤牛型和普通牛进化枝。因此,我们认为平凉地方牛群体和我国其他牛种一样,存在普通牛和瘤牛两个母系起源的遗传背景。     

贵州黄牛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大母系起源,其影响较为均一。并探讨了用核苷酸多样度π值的大小来衡量黄牛群体遗传分化程度的可行性。      

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.     

Complete mitochondrial genomes of Bos taurus and Bos indicus provide new insights into intra-species variation, taxonomy and domestication

The taurine and zebuine cattle breeds comprise the majority of the world cattle population but their taxonomic status is still controversial. The two forms of cattle are currently classified as Bos taurus and Bos indicus species and are differentiated primarily by the presence or absence of a hump. However, these two species hybridize readily, producing fully fertile offspring. We have determined and analyzed complete B. taurus and B. indicus mitochondrial genome sequences to investigate the extent of sequence divergences and to study their taxonomic status by molecular dating. The sequences encompassed 16,338 and 16,339 nucleotides, respectively, and differed at 237 positions. Estimated divergence times indicated that the two cattle lineages separated 1.7-2.0 million years ago. Combined phylogenetic analyses of 18 new and 130 previously reported extant B. taurus and B. indicus control region sequences with data from 32 archaeological specimens of the extinct wild aurochs (Bos primigenius) identified four major maternal lineages. B. primigenius haplotypes were present in all but the B. indicus lineage, and one B. taurus sequence clustered with B. primigenius P haplotypes that were not previously linked with domestic cattle. The B. indicus cluster and a recently reported new B. primigenius haplotype that represents a new lineage were approximately equidistant from the B. taurus cluster. These data suggest domestications from several differentiated populations of B. primigenius and a subspecies status for taurine (B. primigenius taurus) and zebuine (B. primigenius indicus) cattle.     

African-derived mitochondria in South American native cattle breeds (Bos taurus): evidence of a new taurine mitochondrial lineage

This article reports the nucleotide diversity within the control region of 42 mitochondrial chromosomes belonging to five South American native cattle breeds (Bos taurus). Analysis of these data in conjunction with B. taurus and B. indicus sequences from Africa, Europe, the Near East, India, and Japan allowed the recognition of eight new mitochondrial haplotypes and their relative positions in a phylogenetic network. The structure of genetic variation among different hypothetical groupings was tested through the molecular variance decomposition, which was best explained by haplotype group components. Haplotypes surveyed were classified as European-related and African-related. Unexpectedly, two haplotypes within the African cluster were more divergent from the African consensus than the latter from the European consensus. A neighbor-joining tree shows the position of two haplotypes compared to European/African mitochondrial lineage splitting. This different and putatively ancestral mitochondrial lineage (AA) is supported by the calibration of sequence divergence based on the Bos-Bison separation. The European/African mitochondria divergence might be subsequent (67,100 years before present) to that between AA and Africans (84,700 years before present), also preceding domestication times. These genetic data could reflect the haplotype distribution of Iberian cattle five centuries ago.