西藏牦牛mtDNA D-loop区的遗传多样性及其遗传分化

通过测定和分析西藏11个牦牛类群114个个体的mtDNA D-loop区全序列,对西藏牦牛的遗传多样性、类群间的亲缘关系及其遗传分化进行了研究。结果表明:①西藏牦牛mtDNA D-loop区全序列长度为890—896 bp,4种核苷酸T、C、A、G的平均比例分别为28.5%、25.3%、32.4%、13.8%,西藏牦牛mtDNA D-loop区富含碱基A+T,表现出一定的碱基偏好性。②共检测到130个变异位点,占分析总位点数的14.33%;其中单一多态位点85个,占多态位点总数的65.38%,简约信息位点45个,占多态位点总数的34.62%。序列变异中碱基缺失、插入和碱基替换等均有,其中碱基替换变异类型中转换114次,颠换12次,在转换变异类型中以A/G、T/C为主,占95.61%,在颠换变异类型中以A/T为主,占75%。③在114个个体中鉴定出90种单倍型,单倍型多样性为0.981±0.008,核苷酸多样性为0.01056±0.00701,均说明西藏牦牛具有丰富的单倍型类型。④90种单倍型分为2个聚类簇(Ⅰ、Ⅱ),聚类簇Ⅰ包含80种单倍型,占全部单倍型的88.89%,涵盖本研究中所有的西藏牦牛类群;聚类簇Ⅱ中有10种单倍型,占单倍型总数的11.11%,涉及的类群有工布江达、帕里、丁青、巴青、江达、类乌齐、桑桑、桑日、斯布,说明西藏牦牛可能有2个母系起源。⑤西藏牦牛类群间核苷酸分歧度(Dxy)在0.503%—1.416%之间,聚类分析和AMOVA分析显示西藏牦牛可分为两大类,康布牦牛、嘉黎牦牛为一类,其余的牦牛类群为另一类。     

家牦牛线粒体DNA(mtDNA)遗传多样性及其分类

通过分析包括我国10个家牦牛品种(类群)在内共296个样本的mtDNA控制(D-loop)区部分序列的遗传变异,对我国家牦牛的遗传多样性、遗传分化、聚类关系和分类进行了研究.所测序列经比对后,共检测到61个变异位点,定义了77种单倍型.分析显示青海环湖牦牛的单倍型多样性最高,达0.9848±0.0403,而巴州牦牛单倍型多样性最低,为0.8000±0.0825;核苷酸多样性方面,斯布牦牛存在最为丰富的核苷酸序列变异,核苷酸多样性值为0.022582±0.011767,而巴州牦牛仅为0.006856±0.002476,表明我国家牦牛品种(类群)遗传多样性水平存在较大差异.总体上,我国家牦牛单倍型多样性、核苷酸多样性分别为0.9251±0.0095和0.015265±0.007757,呈现出丰富的遗传多样性.聚类分析显示我国家牦牛存在两个聚类簇--斯布牦牛独立为一类;其余9个品种(类群)聚为一类,表明家牦牛品种(类群)间遗传距离与地理分布无明显相关.分子变异分析(AMOVA)显示九龙、嘉黎、斯布牦牛构成的组与其余7个家牦牛品种(类群)构成的组之间存在极显著的遗传分化(?CT=0.05285,P<0.01),且其品种间/组内遗传分化不显著(?SC=0.00648,P>0.05),支持依据遗传分化程度将我国家牦牛划分为两大类型.AMOVA支持的分组在品种(类群)组成上与蔡立的研究结果相符,首次为我国家牦牛划分为横断高山型和青藏高原型两种类型提供了源自分子遗传学的证据.     

新疆克孜勒苏柯尔克孜地区帕米尔牦牛mtDNA COⅠ 的遗传多样性分析

为了更好地保护和开发我国新疆的牦牛遗传资源,对新疆克孜勒苏柯尔克孜地区帕米尔牦牛的遗传多样性及群体遗传结构进行了分析。对3个帕米尔牦牛类群的线粒体COⅠ基因序列进行测定,分析其多态性,构建系统进化树。研究结果显示,帕米尔牦牛COⅠ基因序列全长为1 048 bp,共检测到13个变异位点,其中单态突变位点9个,简约信息位点4个。在帕米尔牦牛类群中共检测出10种单倍型,平均单倍型多样性（haplotype diversity,H_d）、平均核苷酸多样性（nucleotide diversity,P_i）分别为0.762和0.001 64。3个类群间平均遗传距离为0.004 30,且苏木塔什牦牛与布伦口牦牛的遗传距离最小（0.003 72）。结果表明,帕米尔牦牛具有丰富的遗传多样性,但各类群之间遗传分化不明显,且经人工或自然选择后群体发生了扩张和基因交流的情况。研究结果为新疆帕米尔牦牛的遗传改良及新疆优良畜禽遗传资源的保护和开发提供了科学依据。     

西藏牦牛遗传多样性的ISSR分析

选用10对ISSR引物对西藏地区17个牦牛类群(或品种)共850头牦牛进行了PCR扩增,利用Popgen 32软件对数据进行统计分析。结果表明:10对ISSR引物共扩增出134条条带,其中多态性条带132条,多态性比率为95%,扩增片段的大小范围在200~2000bp。17个牦牛类群的平均有效等位基因数为1.498,Nei遗传多样性指数平均值为0.294,平均Shannon多样性指数为0.449,各位点遗传多样性存在一定差别,说明17个西藏牦牛类群具有丰富的遗传多样性。聚类结果表明,西藏牦牛大致分为3类:错那牦牛、江达牦牛、工布江达牦牛、康布牦牛、巴青牦牛、类乌齐牦牛、桑桑牦牛、嘉离牦牛、斯布牦牛和隆子牦牛聚为1类;丁青牦牛、聂荣牦牛、帕里牦牛、日多牦牛、桑日牦牛和申扎牦牛聚为1类;仲巴牦牛单独聚为1类。聚类分析结果与牦牛的地理分布大体一致。     

西藏牦牛微卫星DNA的遗传多样性

为了解西藏牦牛品种或类群的遗传多样性和亲缘关系,文章选用8对微卫星标记引物,利用PCR和复合电泳银染技术,通过计算基因频率(P)、有效等位基因数(Ne)、群体杂合度(He)、多态信息含量(PIC)和遗传距离(D),对西藏11个牦牛类群共480个个体进行了遗传多样性和系统进化分析。结果表明:(1)8个微卫星标记在西藏牦牛类群中均表现出多态性,且均属高度多态位点,遗传多样性丰富;(2)8个微卫星标记的平均多态信息含量在西藏牦牛类群中均高于0.5,其中HEL13最高为0.8496,TGLA57最低为0.7349;在11个牦牛类群中,桑日牦牛的平均多态信息含量最高(0.7949),该群体内部存在较多的遗传变异;丁青牦牛最低(0.7505),则群体相对较纯;(3)在11个牦牛类群中,其杂合度大小分别为:桑日(0.8193)>江达(0.8190)>桑桑(0.8157)>巴青(0.8150)>康布(0.8123)>嘉黎(0.8087)>工布江达(0.8054)>斯布(0.8041)>类乌齐(0.8033)>帕里(0.8031)>丁青(0.7831),西藏东部牦牛的遗传多样性较西部的遗传多样性大,预示西藏东部可能是牦牛的发源地之一;(4)根据遗传距离,用UPGMA法构建聚类关系,表明西藏11个牦牛类群可以分为三大类,即嘉黎牦牛、帕里牦牛、桑桑牦牛、巴青牦牛、类乌齐牦牛、康布牦牛聚为一类,斯布牦牛、工布江达牦牛、桑日额牛、江达牦牛聚为一类,丁青牦牛单独成为一类。综上所述,西藏牦牛的遗传多样性较丰富,所选微卫星标记可用于西藏牦牛遗传多样性的评估。     

中国部分牦牛品种线粒体DNA遗传多态性研究

采用DNA测序技术首次测定了中国5个牦牛品种(类群)35个个体线粒体DNA控制区(D loop)全序列,结果表明:牦牛线粒体DNA控制区全序列长度为891～895bp,T、C、A、G4种核苷酸的含量分别为28.5%、25.3%、32.5%、13.7%。检测到55个变异位点,约占分析位点总数的6.16%,核甘酸变异类型有转换、颠换、插入/缺失4种。确定了24种单倍型,单倍型H4和H6为中国牦牛的主体单倍型,各单倍型在品种间分布不平衡,所测定的5个牦牛品种(类群)单倍型多样度平均为0.9697±0.0180,说明我国牦牛D loop单倍型类型丰富。牦牛品种内各序列平均核苷酸差异数为10.936,核苷酸多样度为1.231%;牦牛品种间核苷酸分歧度(Dxy)约为0.760%～2.155%,品种间双参数距离范围为0.000～0.029,表明我国牦牛遗传多态性丰富。对D环序列变异的分子方差分析和单倍型网络关系图的结果表明:我国牦牛品种间出现显著的遗传分化,牦牛单倍型网络关系图聚为2个聚类簇,表明我国牦牛有2个母系来源,或者有2个主要的驯化地点。     

西藏牦牛SRAP遗传多样性及分类进化研究

用4对SRAP分子标记引物对西藏11个牦牛类群和四川麦洼牦牛的DNA进行扩增,研究其遗传多样性和分类关系。结果表明,在335头牦牛中,共得到29个基因位点,其中有19个多态位点,多态率占65.52%。12个牦牛群体间的Nei’s遗传多样性和Shannon多样性指数分别为0.048 2和0.073 4,遗传相似系数在0.781 1-0.989 1。巴青牦牛和康布牦牛的遗传多样性指数较其他类群高,分别为0.095 5和0.090 0;桑日牦牛类群的遗传多样性指数最低,仅为0.008 5。这些结果表明,12个牦牛类群的SRAP遗传多样性较低。根据Nei’s遗传距离,利用UPGMAM构建聚类关系图结果显示,嘉黎牦牛、帕里牦牛、类乌齐牦牛、桑桑牦牛、康布牦牛、巴青牦牛、丁青牦牛、斯布牦牛和麦洼牦牛聚为一大类,然后依次才与桑日牦牛、工布江达牦牛和江达牦牛相聚在一起,显示在SRAP分子遗传标记所反映的牦牛基因组的遗传结构中,江达牦牛、工布江达牦牛和桑日牦牛与其他牦牛群体间的亲缘关系较远,牦牛的这种亲缘关系与其地理分布也不一致,说明这12个牦牛的起源、演化关系较复杂,有待于进一步研究分析。     

中国近海重要生态建群红藻真江蓠的群体遗传多样性

真江蓠(Gracilaria vermiculophylla)是中国近海潮间带生态系统结构组成和功能维持的重要支撑物种, 但有关其群体遗传结构和多样性分布模式的研究目前仍较缺乏。本研究利用线粒体cox1序列对我国近海19个真江蓠地理群体进行了系统发育和群体遗传分析。461个长度为641 bp的cox1序列片段共含有21个多态位点, 产生15个单倍型。基于cox1序列的系统进化分析、单倍型分析和主成分分析显示, 19个真江蓠群体分化为南北两个类群, 其中浙江嵊泗以北的13个群体形成北方类群, 福建厦门以南的6个群体形成南方类群。遗传距离和分子方差分析显示真江蓠南北各类群内的遗传分化较小, 南北类群间的遗传分化达到亚种水平。南北类群间的差异是我国近海真江蓠群体遗传变异的主要来源。     

淮河水系日本沼虾群体遗传结构和系统演化的线粒体COI序列分析

测定了淮河水系17个日本沼虾(Macrobrachium nipponense)野生群体共248个个体的线粒体细胞色素氧化酶亚基I(COI)部分序列,获得623 bp核苷酸片段,包括48个变异位点,定义了31个单倍型,共享单倍型有12个,整体单倍型多样性和平均核苷酸多样性均处于中间水平。AMOVA分析表明,17个群体间的遗传分化系数Fst=0.041 3(P0.05),群体间遗传分化较小。Kimura 2-paramter遗传距离在五河与焦岗湖、花家湖及瓦埠湖群体间最大,为0.014,在高邮和邵伯湖群体之间最小,为0.003。MP系统树与单倍型进化网络关系图具有较高的一致性,31个单倍型被分为3个进化枝,其中一个进化枝主要以下游群体为主,另外2个进化枝主要以中游群体为主。群体中性检验、错配分析表明,淮河日本沼虾近期曾经历过种群扩张。     

长江中游鳊群体的遗传多样性

采用线粒体DNA(mt DNA)细胞色素b基因(cytb)和控制区序列,分析长江中游宜都、沙市、燕窝、团风等4个产卵场鳊种群的遗传多样性和遗传结构。结果表明:长江中游鳊群体cytb序列共检出82个多态位点,86种单倍型,平均单倍型多样性指数(Hd)和核苷酸多样性指数(Pi)分别为0.930和0.00244;控制区序列共检出变异位点46个,单倍型76种,Hd指数和Pi指数分别为0.972和0.00505;构建的单倍型网络结构和分子方差分析(AMOVA)显示,4个群体的遗传变异绝大部分来自群体内部,群体间无显著遗传分化;群体间的分化指数(FST)、平均基因流(Nm)和平均K2-P遗传距离均表明,4个鳊地理群体间存在广泛的基因交流,未发生明显群体遗传分化;中性检验表明,鳊历史上发生了群体扩张,扩张时间在第四纪冰期后期。     

Phylogeographical analysis shows the need to protect the wild yaks' last refuge in Nepal

The wild yak Bos mutus was believed to be regionally extinct in Nepal for decades until our team documented two individuals from Upper Humla, north‐western Nepal, in 2014. The International Union for Conservation of Nature (IUCN) seeks further evidence for the conclusive confirmation of that sighting. We conducted line transects and opportunistic sign surveys in the potential wild yak habitats of Humla, Dolpa, and Mustang districts between 2015 and 2017 and collected genetic samples (present and historic) of wild and domestic yaks Bos grunniens. We also sighted another wild yak in Upper Humla in 2015. Phylogenetic and haplotype network analyses based on mitochondrial D‐loop sequences (~450 bp) revealed that wild yaks in Humla share the haplotype with wild yaks from the north‐western region of the Qinghai‐Tibetan Plateau in China. While hybridization with domestic yaks is a major long‐term threat, illegal hunting for meat and trophy put the very small populations of wild yaks in Nepal at risk. Our study indicates that the unprotected habitat of Upper Humla is the last refuge for wild yaks in Nepal. We recommend wild yak conservation efforts in the country to focus on Upper Humla by (i) assigning a formal status of protected area to the region, (ii) raising awareness in the local communities for wild yak conservation, and (iii) providing support for adaptation of herding practice and pastureland use to ensure the viability of the population.     

Comprehensive investigation of nucleotide diverdity in yaks

To understand the maternal genetic diversity of Tianzhu white yak better, we analyzed mtDNA D‐loop sequences of 209 Tianzhu white yaks, which are from the central region of Tianzhu white yak habitat. Accordingly, a total of 45 haplotypes were identified in Tianzhu white yaks in this study, and 18 of them were unique. The nucleotide diversity and haplotype diversity of population studied were 0.024 ± 0.003 and 0.946 ± 0.007 respectively, revealing that Tianzhu white yak possess a relatively high genetic diversity. The phylogenetic analysis, combining D‐loop sequences in this study with 533 previous published D‐loop sequences of 13 yak breeds, indicated that Tianzhu white yaks fell mainly into haplogroup A and that a small portion belonged to haplogroups B, C, D and E. Moreover, six haplotypes of 20 individuals identified in Tianzhu white yak were in the taurine haplogroup, indicating hybridization between Bos taurus and Tianzhu white yaks. In summary, this study supplies a comprehensive maternal genetic pattern for Tianzhu white yak and provides a basic reference for future breeding programs to conserve the purebred Tianzhu white yak.     

Comparative Study on the Genetic Diversity of GHR Gene in Tibetan Cattle and Holstein Cows

Due to the phenotype-based artificial selection in domestic cattle, the underlying functional genes may be indirectly selected and show decreasing diversity in theory. The growth hormone receptor (GHR) gene has been widely proposed to significantly associate with critical economic traits in cattle. In the present study, we comparatively studied the genetic diversity of GHR in Tibetan cattle (a traditional unselected breed, n = 93) and Chinese Holstein cow (the intensively selected breed, n = 94). The Tibetan yak (n = 38) was also included as an outgroup breed. A total of 21 variants were detected by sequencing 1279 bp genomic fragments encompassing the largest exon 9. Twelve haplotypes (H1～H12) constructed by 15 coding SNPs were presented as a star-like network profile, in which haplotype H2 was located at the central position and almost occupied by Tibetan yaks. Furthermore, H2 was also identical to the formerly reported sequence specific to African cattle. Only haplotype H5 was simultaneously shared by all three breeds. Tibetan cattle showed higher nucleotide diversity (0.00215 ± 0.00015) and haplotype diversity (0.678 ± 0.026) than Holstein cow. Conclusively, we found Tibetan cattle have retained relatively high genetic variation of GHR. The predominant presence of African cattle specific H2 in the outgroup yak breed would highlight its ancestral relationship, which may be used as one informative molecular marker in the phylogenetic studies.     

锁阳ITS序列遗传多样性分析

为阐明锁阳(Cynomorium songaricum)的遗传结构与遗传多样性, 以甘肃省河西走廊地区及青海省共18个居群的188个锁阳个体为研究对象, 利用现代分子生物学技术, 采用序列分析方法从核-质基因方面对遗传结构进行了分析。结果显示, 锁阳ITS序列总长度为687 bp, 含有7个变异位点, 定义9个单倍型, 整体单倍型多态性H_d=0.294 20, 核苷酸多样性π=0.000 49。在整个单倍型网络中介图中, 单倍型H1位于中心位置, 并在所有的居群中均有分布, 为核心的古老单倍型。分子方差分析结果显示, 锁阳种群变异主要来源于种群内。根据ITS序列得到的群体间遗传分化系数以及Mantel检验结果, 锁阳种群间的遗传距离与地理距离之间不存在相关性, 表明现存的锁阳居群是相对近期发生生境片段化的产物。中性检验结果表明, 锁阳拒绝中性进化, 群体历经扩张或者基因座位受到负选择作用, 其中性零假说不能被排除。该研究为锁阳的系统分类、资源鉴定以及保护措施的制定提供了分子证据。     

青藏高原濒危植物唐古特大黄的遗传多样性

唐古特大黄(Rheum tanguticum)是中国传统的中藏药材,近几年由于生境的严重破坏,已濒临灭绝,并被列入濒危植物名单。为了探索唐古特大黄物种濒危的原因并保护其野生资源,本研究采集了9个居群87个个体的唐古特大黄样本,基于该物种的叶绿体基因trn S-G序列对其进行了遗传多样性研究。结果表明,唐古特大黄物种具有较高的遗传多样性水平(Ht=0.694),其中95.97%的遗传分化来自于居群间(G_(ST)=0.960),4.03%的遗传分化来自于居群内(Hs=0.028)。AMOVA分析也显示唐古特大黄居群间基因流较小(N_m=0.01),存在较高的遗传分化(F_(ST)=0.9631)。唐古特大黄较高的遗传多样性水平可能与该物种较长的进化史和生活史有关,居群间较高的遗传分化可能与高山地区特殊的地理环境和人类活动有关。根据研究结果,建议对唐古特大黄所有野生居群进行就地保护,同时收集种质资源开展异地繁殖工作,以保护物种的遗传多样性,维持其进化潜力。     

藏羚羊mtDNA D-Loop区遗传多样性研究

该研究采用非损伤性DNA基因分型技术,对可可西里地区10个藏羚羊（Pantholops hodgsonii）个体的mtDNA非编码区部分片段（444~446bp）进行了序列分析,结果显示A、T%含量（61.8%）明显高于G、C%含量（38.2%）,共发现10种单倍型,包括48个多态位点,其中转换位点44个、颠换位点1个、插入位点1个、缺失位点2个。单倍型间平均遗传距离为0.031,单倍型多态性（h）为1.000,核苷酸多态性（π）为2.96%。说明藏羚羊线粒体控制区存在着丰富变异,最后从藏羚羊的生态习性及地理分布两方面对这一结果进行了分析探讨。      

Analysis of genetic diversity and population structure of Chinese yak breeds （Bos grunniens） using microsatellite markers

Nine Chinese yak breeds (Maiwa,Tianzhu White,Qinghai Plateau,Sibu,Zhongdian,Pall,Tibetan High Mountain,Jiulong,and Xin-jiang) and Gayal were analyzed by means of 16 microsatellite markers to determine the level of genetic variation within populations,genetic relationship between populations,and population structure for each breed.A total of 206 microsatellite alleles were observed.Mean F-statistics (0.056) for 9 yak breeds indicated that 94.4% of the genetic variation was observed within yak breeds and 5.6% of the genetic variation existed amongst breeds.The Neighbor-Joining phylogenetic free was constructed based on Nei's standard genetic dis-tances and two clusters were obtained.The Gayal separated from the yaks far away and formed one cluster and 9 yak breeds were grouped together.The analysis of population structure for 9 yak breeds and the Gayal showed that they resulted in four clusters; one clus-ter includes yaks from Tibet Autonomous Region and Qinghai Province,one cluster combines Zhongdian,Maiwa,and Tianzhu White,and Jiulong and Xinjiang come into the third cluster.Pali was mainly in the first cluster (90%),Jiulong was mainly in the second cluster (87.1%),Zhongdian was primarily in the third cluster (83%),and the other yak breeds were distributed in two to three clusters.The Gayal was positively left in the fourth cluster (99.3%).