牦牛Myf-5基因的克隆测序及其系统进化研究

本研究对牦牛(九龙牦牛)的生肌决定因子5(Myf-5)基因进行了T-A克隆测序和分析,并与多个物种的相应基因编码区核苷酸序列、氨基酸序列进行了比对分析,构建了物种间的系统进化树.结果 表明:①牦牛的Myf-5基因大小为3313 bp,由3个外显子和2个内含子组成,与普通牛等9个物种比较,在基因大小上有较大的差异,但外显子和内含子的组成一致.②牦牛、大额牛、瘤牛、水牛、普通牛、人、恒河猴、黑猩猩、狗、家鼠、软体贝壳、鸡、斑马鱼等物种间Myf-5基因编码区的核苷酸序列同源性较高,其中,牦牛、大额牛、瘤牛、水牛、普通牛间的同源性最高,达98.4%以上,说明Myf-5基因编码区核苷酸序列在动物物种间具有较高的保守性.③牦牛、大额牛、瘤牛、水牛、普通牛、黑猩猩、恒河猴、人、狗、家鼠、软体贝壳、鸡、斑马鱼等物种间Myf-5基因编码蛋白的氨基酸序列具有较高的同源性,保守性强,这一结果与编码区核苷酸序列的比对结果基本一致.④根据核苷酸序列,用NJ法构建的牦牛、大额牛、瘤牛、水牛、普通牛、人、恒河猴、黑猩猩、狗、家鼠、软体贝壳、鸡、斑马鱼等13个物种的分子系统进化树显示:牦牛、普通牛、瘤牛、水牛和大额牛在较近的亲缘关系下聚为一大类,人、黑猩猩和恒河猴聚为另一大类,然后这两类再和其他物种相聚.这一分类结果与各物种的动物学分类结果和血液蛋白、mtDNA水平上的聚类结果基本一致,支持牦牛、普通牛和瘤牛3个物种间不应该是属间或亚属间关系,而应是同一属下的不同种,将牦牛、普通牛和瘤牛划分在同一个属--牛属(Bos),而将水牛划分在另一个单独的属的观点.同时也显示该基因序列适合用于动物学分类.     

雷琼牛GH基因第五外显子遗传多样性及其分化

对16头雷琼牛GH基因第5外显子序列进行分析,发现了1个变异位点,定义了2种单倍型。引用巴州牦牛2个个体GH基因同源区序列并结合GenBank中牛属普通牛、其它瘤牛和牦牛3个种群与水牛1个远缘种GH基因同源区序列,分别采用邻接（NJ）法和最大简约（MP）法构建分子系统发育树,得到基本一致的拓扑结构,结果显示GH基因的分化早于雷琼牛（瘤牛）、其它瘤牛、普通牛、牦牛和水牛的分化,瘤牛物种内存在多型,同时证实了GH基因第5外显子区有着较高的突变率。     

水牛、大额牛和牦牛抑制素α亚基编码基因外显子1多态性分析

为了揭示牛科物种INHA基因的遗传特征,该文采用PCR产物直接测序法对水牛、大额牛和牦牛INHA基因外显子1及其侧翼序列进行多态性检测,并结合已发表的包括牛科物种在内的一些哺乳动物数据进行了比较分析。结果表明,在水牛INHA基因外显子1中存在c.73C>A替换,为同义替换,河流型和沼泽型水牛编码产物一致;在大额牛的INHA基因外显子1中发现c.62C>T、c.187G>A替换,分别引起INHA中氨基酸发生p.P21L、p.V63M改变,两者均为相同性质氨基酸的替换;在牦牛中发现c.62C>T、c.129A>G替换,前者也引起编码氨基酸发生p.P21L替换,后者为同义替换。在INHA基因5’侧翼区所测出的序列中,水牛、大额牛和牦牛等物种内均未发现SNP位点,但在种间发现存在c.-6T>G的替换,大额牛、牦牛和普通牛均为c.-6G,而水牛为c.-6T。在INHA基因内含子中,水牛的第31～36位核苷酸处发现有6个碱基的缺失,即c.262+31₂62+36delTCTGAC;该位点在河流型水牛中野生型（+/+）占主体,而在沼泽型水牛中则缺失型（-/-）占主体。在大额牛、牦牛和普通牛等其它牛科物种的内含子中均未发现该缺失,但与水牛相比,大额牛、牦牛和普通牛内含子中发现缺失c.262+78₂62+79delTG。序列比对显示,INHA基因外显子1序列中c.43A和c.67G为水牛中所特有,而c.173A和c.255G为大额牛、牦牛和普通牛所共有,c.24C、c.47G、c.174T和c.206T为山羊所特有。大额牛、牦牛和普通牛间INHA基因外显子1序列差异较小,而山羊和水牛与它们间的差异相对较大。     

牦牛与其他物种ZFX/ZFY基因片段间的进化关系

利用PCR扩增、克隆和序列分析法对牦牛ZFX/ZFY基因第11外显子部分片段进行了研究,并同来自于NCBI GenBank中人、猩猩、普通牛等9个物种的ZFX/ZFY基因核苷酸及其氨基酸序列进行了进化分析.结果表明,牦牛ZFX、ZFY基因间核苷酸序列同源性为94.1%,显示同一物种同源基因ZFX/ZFY间存在变异;比较的10个物种间ZFX基因核苷酸序列同源性为87.7%、ZFY基因为81.7%,相应ZFX、ZFY氨基酸同源性分别为96.6%、91.0%,ZFY基因的变异性大于ZFX基因,显示X染色体与Y染色体可能是独立进化.     

以Cytb基因部分序列分析大额牛系统地位及其种群现状

对12头大额牛Cytb基因部分序列447bp进行分析,共发现33个变异位点,定义了4种单倍型;结合GenBank中牛属普通牛、瘤牛、牦牛和印度褐牛4个近缘种的Cytb基因同源区序列进行分析,以水牛作为外群,分别采用邻接(NJ)法和最大简约(MP)法构建分子系统发育树,得到基本一致的拓扑结构,结果表明大额牛是独立于普通牛、瘤牛和印度褐牛之外的Bos属的一个独立的种。同时还指出目前我国有相当比例的大额牛受到外来物种的入侵,我国的大额牛保护工作正面临非常严峻的考验。     

务川黑牛mtDNA遗传变异及起源进化分析

目的:测定务川黑牛mtDNA D-loop区全序列,以了解其遗传背景。方法:采用PCR直接测序方法测定务川黑牛mtD-NA D-loop区全序列。结果:32个个体D-loop区的全序列中,四种碱基A、T、C和G含量分别为:33.1%、28.4%、25.0%和13.5%,A+T平均含量61.5%。经过序列比对,共检测到57个变异位点,占所测mtDNA序列总长的6.26%,其中转换54个,颠换3个;产生17种单倍型,其中6种是普通牛单倍型,11种是瘤牛单倍型;平均核苷酸差异为22.623,单倍型多样度0.901±0.039,核苷酸多样度0.0249±0.00147。结论:务川黑牛这一品种遗传多样性非常丰富,据此构建的NJ进化树显示务川黑牛含有普通牛和瘤牛的血统。     

基于CO Ⅰ基因的龟鳖目动物鉴定与分类进化分析

为探索COⅠ基因条形码技术在龟鳖目物种分类与鉴定中的适用性,本研究应用通用引物PCR扩增获得龟鳖目动物9科26属45种共205个样品的634 bp COⅠ基因序列片段进行研究分析。结果显示,COⅠ基因在龟鳖目动物中存在少量的碱基插入和缺失,共计12种20个位点,物种占比34.28%;平均碱基含量(A+T)为55.42%,显著高于(G+C)的44.58%。通过Kimura's 2-parameter计算45个物种的种内遗传距离为0～0.023 98,平均为0.003 2;种间遗传距离为0.022 3～0.334 9,平均为0.1 637,种间平均遗传距离约为种内平均遗传距离的51.16倍;有3个物种的种内遗传距离大于Hebert推荐区分物种的最小种间遗传距离0. 02,其中1个物种的种内遗传距离大于45个物种的种间最小遗传距离0.022 3。同科内不同属间及同目内不同科间的遗传距离分别为0.146 7、0.212 7,显示物种间遗传距离随着分类阶元的上升而逐步扩大。系统进化树聚类结果显示,不同个体依据亲缘关系程度不同分别聚为独立的一支,没有出现物种交叉分布,能准确反映个体间的亲缘关系。结果表明利用COⅠ基因鉴定龟鳖目物种具有可行性,但科内属间及以上分类阶元进化节点自展率不高,无法通过COⅠ基因明确其进化关系,须结合其他DNA指标加以确认。     

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

中国荷斯坦牛和鲁西黄牛mtDNAD-loop序列多态性分析

为了从母系遗传角度深入阐明中国荷斯坦牛和鲁西黄牛的群体遗传多样性以及起源进化,本研究采用PCR测序法测定了9头中国荷斯坦牛和11头鲁西黄牛的线粒体DNA D-loop区的部分序列,并经剪切后进行生物信息学软件分析.结果表明,20个个体D-loop区共711 bp,共检测到19种单倍型和50个多态位点.核苷酸多样性(Pi)为0.021 33±0.004 54,单倍型多样性(Hd)为0.994±0.019,平均核苷酸差异数(k)为15.146 20.构建的NJ网络进化树共分为两大支系,其中部分鲁西黄牛与瘤牛聚为一支,而中国荷斯坦牛和部分鲁西黄牛与普通黄牛聚为一支.说明中国荷斯坦牛和鲁西黄牛群体遗传多样性均较高;鲁西黄牛同时含有瘤牛和普通黄牛的血统,而中国荷斯坦牛只含有普通黄牛的血统.     

利用微卫星标记分析贵州白水牛与普通水牛的遗传多样性与亲缘关系

利用11对微卫星引物对贵州白水牛和6个普通水牛群体的有效等位基因数、基因杂合度、多态信息含量和遗传距离进行了分析。结果表明,11个微卫星基因座在7个水牛群体中均存在多态性,可以用于水牛的遗传多样性评估。贵州白水牛和6个普通水牛群体的平均杂合度和平均多态信息含量分别为0·7450～0·7922和0·7021～0·7605。贵州白水牛和遵义普通水牛的亲缘关系最近,遗传距离为0·0910。由UPGMA聚类法得到的系统聚类图反映了贵州白水牛和6个普通水牛群体的亲缘关系远近程度,贵州白水牛没有单独聚为一类,而是与同分布区的遵义普通水牛首先聚类,然后依次与其余地区的普通水牛聚类。研究提出了贵州白水牛应是其产地的普通水牛群体中的一个突变群的见解,为开展贵州白水牛的遗传资源评估、保护与利用提供了分子水平的遗传学依据。     

野牦牛线粒体基因组序列测定及其系统进化

野牦牛属高寒地区的特有物种,是我国最珍贵的野生动物遗传资源之一,已被列为国家一级重点保护动物。对野牦牛mtDNA进行全序列测定和结构分析,并基于线粒体基因组序列对其系统发生进行了探讨。结果表明:(1)野牦牛线粒体基因组全序列的大小为16 322 bp,整个基因组由37个编码基因和D-loop区组成;22个tRNA基因序列长度为1 524 bp、2个RNA基因序列长度为2 528 bp、13个编码蛋白基因序列长度为11420 bp、D-loop区长度为892 bp。基因组中无间隔序列,基因间排列紧密,基因内无内含子。(2)野牦牛具有较丰富的遗传多样性。(3)分子系统发生关系显示牦牛为牛亚科中的一个独立属,即牦牛属(Poephagus),牦牛属包括家牦牛(Poephagus grunniens)和野牦牛(Poephagus mutus)2个种。野牦牛线粒体基因组全序列的获得和结构解析对研究牦牛的起源、演化和分类,以及野牦牛遗传资源的保护、开发和利用均具有重要的理论和实际意义。     

Complete sequence of the yak (Bos grunniens) mitochondrial genome and its evolutionary relationship with other ruminants

The yak (Bos grunniens) is the most important domesticated species in the Qinhai-Tibetan Plateau. In present study, the complete sequence of the yak mitochondrial genome was determined. Sequence analysis revealed that there are no differences with cattle in the yak mitochondrial genome organization. Interestingly, within the D-loop, the conserved sequence blocks are less conserved than surrounding regions. Neighbor-Joining (NJ) trees based on single genes, gene sets and concatenated genes of mitochondrial genome were constructed. The analysis identified the yak as a sister group of a cattle/zebu clade. Based on substitutions in 22 tRNA genes, 12S rRNA gene and 16S rRNA gene, the dating of divergence between yak and cattle/zebu, and yak and water buffalo, was proposed to have occurred 4.38-5.32 and 10.54-13.85 million years before present, respectively. This is consistent with the paleontologyical data. Yak and sheep/goat divergent dating predicts that their divergence occurred at 13.14-27.99 million years before the present day.     

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.     

Phylogeny of bovine species based on AFLP fingerprinting

The Bovini species comprise both domestic and wild cattle species. Published phylogenies of this tribe based on mitochondrial DNA contain anomalies, while nuclear sequences show only low variation. We have used amplified fragment length polymorphism (AFLP) fingerprinting in order to detect variation in loci distributed over the nuclear genome. Computer-assisted scoring of electrophoretic fingerprinting patterns yielded 361 markers, which provided sufficient redundancy to suppress stochastic effects of intraspecies polymorphisms and length homoplasies (comigration of non-homologous fragments). Tree reconstructions reveal three clusters: African buffalo with water buffalo, ox with zebu, and bison with wisent. Similarity values suggest a clustering of gaur and banteng, but bifurcating clustering algorithms did not assign consistent positions to these species and yak. We propose that because of shared polymorphisms and reticulations, tree topologies are only partially adequate to represent the phylogeny of the Bovini. Principal-coordinate analysis positions zebu between a gaur/banteng cluster and taurine cattle. This correlates with the region of origin of these species and suggests that genomic distances between the cattle species have been influenced by genetic exchange between neighbouring ancestral populations.     

Evolution and domestication of the Bovini species

Domestication of the Bovini species (taurine cattle, zebu, yak, river buffalo and swamp buffalo) since the early Holocene (ca. 10 000 BCE) has contributed significantly to the development of human civilization. In this study, we review recent literature on the origin and phylogeny, domestication and dispersal of the three major Bos species – taurine cattle, zebu and yak – and their genetic interactions. The global dispersion of taurine and zebu cattle was accompanied by population bottlenecks, which resulted in a marked phylogeographic differentiation of the mitochondrial and Y-chromosomal DNA. The high diversity of European breeds has been shaped through isolation-by-distance, different production objectives, breed formation and the expansion of popular breeds. The overlapping and broad ranges of taurine and zebu cattle led to hybridization with each other and with other bovine species. For instance, Chinese gayal carries zebu mitochondrial DNA; several Indonesian zebu descend from zebu bull × banteng cow crossings; Tibetan cattle and yak have exchanged gene variants; and about 5% of the American bison contain taurine mtDNA. Analysis at the genomic level indicates that introgression may have played a role in environmental adaptation.     

Analysis of conserved microsatellite sequences suggests closer relationship between water buffalo Bubalus bubalis and sheep Ovis aries.

The distribution and evolutionary pattern of the conserved microsatellite repeat sequences (CA)n, (TGG)6, and (GGAT)4 were studied to determine the divergence time and phylogenetic position of the water buffalo, Bubalus bubalis. The mean allelic frequencies of these repeat loci showed a high level of heterozygosity among the euartiodactyls (buffalo, cattle, sheep, and goat). Genetic distances calculated from the allelic frequencies of these microsatellites were used to position Bubalus bubalis in the phylogenetic tree. The tree topology revealed a closer proximity of the Bubalus bubalis to the Ovis aries (sheep) genome than to other domestic species. The estimated time of divergence of the water buffalo genome relative to cattle, goat, sheep, pig, rabbit, and horse was found to be 21, 0.5, 0.7, 94, 20.3, and 408 million years (Myr), respectively. Although water buffaloes share morphological and biochemical similarities with cattle, our study using the microsatellite sequences places the bubaline species in an entirely new phylogenetic position. Our results also suggest that with respect to these repeat loci, the water buffalo genome shares a common ancestry with sheep and goat after the divergence of subfamily Bovinae (Bos taurus) from the family Bovidae.     

Species authentication of commercial beef jerky based on PCR-RFLP analysis of the mitochondrial 12S rRNA gene

In this study,we determined species-specific variations by analyzing the mitochondrial 12S rRNA gene sequence variation(～440 bp) in 17 newly obtained sequences and 90 published cattle,yak,buffalo,goat,and pig sequences,which represent 62 breeds and 17 geographic regions.Based on the defined species-specific variations,two endonucleases,Alu I and Bfa I,were selected for species authentication using raw meat/tissue samples and the PCR-RFLP method.Goat and pig were identified using the Alu I enzyme,while cattle,yak,and buffalo were identified by digestion with Bfa I.Our approach had relatively high detection sensitivity of cattle DNA in mixed cattle and yak products,with the lowest detectable threshold equaling 20% of cattle DNA in a mixed cattle/yak sample.This method was successfully used to type commercial beef jerky products,which were produced by different companies utilizing various processing technologies.Our results show that several yak jerky products might be implicated in commercial fraud by using cattle meat instead of yak meat.     

Mitochondrial DNA sequence diversity and origin of Chinese domestic yak

In order to clarify the origin and genetic diversity of yak in China, we analysed mitochondrial DNA (mtDNA) control region sequences (approximately 891 bp) in 52 individuals from four domestic yak (Poephagus grunniens) breeds, as well as from a hybrid between yak and cattle (Pianniu). Twenty-five samples were further selected for partial (420 bp) cytochrome b sequencing based on control region sequence information. Two yak samples shared sequences with Chinese cattle (Bos taurus); the remaining yak mtDNAs converged into two major clades in the phylogenetic analysis. Genetic diversity varied substantially among the breeds, with the hybrid Pianniu yak demonstrating the highest diversity. Our results suggest that the Chinese yak was domesticated from two distinct matrilineal sources or from a heterogeneous pool containing both divergent lineages, with occasional gene introgression from cattle.