Genetic Diversity Analysis of Shaanxi Moschus berezovskii Based on mtDNA Cytb Gene and D-loop Region Sequence

【Objective】 This study was aimed to explore the genetic diversity of Shaanxi Moschus berezovskii population,and understand the genetic information of Moschus berezovskii.【Method】 The hair of Moschus berezovskii was collected to extract DNA,the mitochondrial DNA(mtDNA) cytochrome b(Cytb) gene and D-loop sequences of 43 Moschus berezovskii individuals were determined,and the base composition was counted.All sequences were integrated and compared using ClustalX 2.0 software to obtain nucleotide polymorphic sites (SNPs) in the population.The nucleotide diversity (Pi),number of haplotype (H),haplotype diversity (Hd) and average number of nucleotide differences (K) were calculated by DNASP 5.10 software.The genetic distance among different haplotypes of Cytb gene and D-loop sequences was calculated by Mega 7.0 software,and Neighbor-Joining (NJ) phylogenetic tree was constructed.【Result】 The AT content of Cytb gene and D-loop region were higher than GC content,indicating there was bias in base composition.There were 241 and 383 SNPs of Cytb gene and D-loop region,respectively.The nucleotide diversity of Cytb gene and D-loop region were 0.28343 and 0.07707,and the haplotype diversity was 0.983 and 0.975,respectively,indicating that the population genetic diversity was rich.The genetic distances of 35 haplotypes of Cytb gene ranged from 0.002 to 0.831,and 29 haplotypes of D-loop region ranged from 0.006 to 1.342.The phylogenetic tree showed that there were two mitochondrial lineages,indicating that there were two mitochondrial maternal origins.The evolutionary analysis of D-loop region also supported this conclusion.【Conclusion】 The nucleotide diversity and haplotype diversity of Moschus berezovskii population were high,and the genetic diversity was rich.At the same time further supported the view of Moschus berezovskii and Moshus moschiferus belonged to a branch of the view.     

基于网络药理学和分子对接探究归芪益母口服液治疗猪气虚血瘀证的作用机制

【目的】运用网络药理学和分子对接技术探讨归芪益母口服液治疗猪气虚血瘀证的有效活性成分及其作用机理。【方法】通过中药系统药理学数据库与分析平台(TCMSP)、中医药综合数据库(TCMID)获取归芪益母口服液主要活性成分及其靶点,通过GeneCards、OMIM和Malacards数据库获得母猪气虚血瘀证的作用靶点。使用Cytoscape 3.9.0软件构建药物-成分-靶点网络图,利用STRING数据库和Cytoscape 3.9.0软件构建靶点蛋白-蛋白互作(PPI)网络图,采用DAVID数据库对共同作用靶点进行GO功能和KEGG信号通路富集分析。通过Autodock Tools 1.5.7软件对其核心成分与核心靶点进行分子对接。【结果】归芪益母口服液筛选得到主要的有效活性成分为槲皮素、山奈酚和异鼠李素,其对应76个关键作用靶点,6 335个疾病靶点,41个交集靶点;PPI网络显示肿瘤坏死因子(TNF)、白介素-6(IL6)和IL10为关键靶点。GO功能富集分析共筛选出12条目,其中生物过程7条目,包括免疫应答、STAT蛋白酪氨酸磷酸化的正向调节、对糖皮质激素的反应等;细胞组分2条目,包括细胞外组分和细胞外空隙;分子功能3条目,包括细胞因子、生长因子活性和IL2受体结合。KEGG信号通路富集分析显示,归芪益母口服液主要通过细胞因子受体互作、JAK-STAT信号通路、PI3K-Akt信号通路发挥作用。分子对接结果显示,槲皮素、山奈酚、异鼠李素等核心成分与母猪气虚血瘀证关键靶点的结合能均＜0 kJ/mol。【结论】归芪益母口服液可能通过槲皮素、山奈酚、异鼠李素等主要活性成分作用于核糖体结合蛋白(RPN1)、胰核糖核酸酶(RNASE1)、血管生成蛋白抑制因子1(RNH1)等靶点,通过参与细胞因子受体互作、JAK-STAT信号通路、PI3K-Akt信号通路等治疗猪气虚血瘀证。     

阿坝藏族羌族自治州若尔盖地区藏猪mtDNA D-Loop的遗传多样性分析

旨在通过获得和对比若尔盖地区藏猪mtDNA D-Loop高变区的部分序列,为该地区藏猪遗传资源的保护与利用提供参考。本试验收集了若尔盖地区9个乡镇共80头藏猪耳组织,以甘南州藏猪mtDNA D-Loop基因序列为模板设计引物,采用PCR扩增测序技术获得了80条核苷酸序列,并采用生物信息学的数据处理方法进行分析。结果表明,若尔盖地区藏猪mtDNA D-Loop高变区（435 bp）的A+T含量（56.46%）明显高于G+C含量（43.54%）,存在碱基偏倚性现象;在80条长度为435 bp的序列中,共检测到14个变异位点,鉴定了17个单倍型,单倍型多样度（Hd）、核苷酸多样度（Pi）和平均核苷酸差异数（k）分别为0.881、0.004 66和2.028,其中Hd、Pi和k在降扎乡藏猪群体中最高,Pi和k在占哇乡藏猪中最低;共享单倍型8个,特有单倍型9个,且若尔盖地区不同藏猪群体间特有单倍型数差异较大,其中,降扎乡藏猪的特有单倍型数量最多,占单倍型总数的17.65%（3/17）;Hap_1和Hap_15单倍型是4个乡镇（降扎乡、益哇乡、热尔乡和冻列乡）群体的共享单倍型,表明这4个乡镇藏猪群体存在两个共同的母系祖先单倍型。若尔盖地区藏猪的平均遗传距离为0.004 66,其中降扎乡和冻列乡藏猪间的遗传距离最大为0.006 90,包座乡和益哇乡藏猪间的遗传距离最小为0.002 16;构建的NJ分子系统进化树将若尔盖地区藏猪分为2支,而加入中国地方家猪、野猪和引种猪后,若尔盖地区藏猪在进化树中比较分散,说明该地区藏猪母源血统遗传复杂,彼此之间基因交流多。本研究进一步证实了若尔盖地区藏猪比西藏林芝、山南、日喀则、甘孜州、阿坝州藏猪的遗传多样性程度高,受到人工选择强度低,应强化对若尔盖地区藏猪遗传资源的保护与利用。     

猪TRIM3基因克隆、生物信息学及组织表达分析

【目的】克隆大白猪三基序结合蛋白3（tripartite motif-containing 3,TRIM3）基因,并对其进行生物信息学和组织表达分析。【方法】采用PCR技术扩增并克隆大白猪TRIM3基因CDS全长序列,连接pMD18-T载体并转化大肠杆菌DH5α感受态细胞,通过蓝白斑筛选阳性克隆,菌液PCR鉴定后测序,与不同物种TRIM3基因序列比对并构建系统进化树;应用多种在线软件对其编码蛋白进行生物信息学分析,并利用实时荧光定量PCR方法检测TRIM3基因在大白猪不同组织中的相对表达量。【结果】大白猪TRIM3基因CDS序列全长2 235 bp,编码744个氨基酸。相似性和遗传进化分析结果显示,大白猪与野猪的相似性最高,达99.7%,与鸭的相似性最低,为75.1%;大白猪TRIM3基因与野猪先聚为一类,与牛和山羊亲缘关系较近。生物信息学分析显示,大白猪TRIM3蛋白分子质量为80.58 ku,理论等电点（pI）为8.32,不稳定系数为40.85,为亲水性蛋白,但不是分泌蛋白,无糖基化位点,预测其有60个磷酸化位点,主要存在于细胞质内;在TRIM3蛋白二级结构中以无规则卷曲为主,占41.67%,三级结构模型预测结果与二级结构一致。组织表达分析表明,大白猪TRIM3基因在心脏、肝脏、脾脏、肺脏、肾脏、肌肉、气管、结肠中均有分布,肺脏中表达量最多且显著高于其他组织（P＜0.05）。【结论】本研究成功克隆大白猪TRIM3基因CDS全长序列,并进行了生物信息学和组织表达分析,为进一步研究大白猪TRIM3蛋白的免疫学功能提供理论依据,对探究大白猪TRIM3基因参与先天性免疫和抗病毒感染分子机制具有重要意义。     

水牛TRAIL基因克隆及生物信息学分析

【目的】对水牛肿瘤坏死因子(tumor necrosis factor,TNF)相关凋亡诱导配体(TNF-related apotosis-inducing ligand,TRAIL)基因CDS序列进行克隆及序列分析,并对其编码的蛋白进行生物信息学分析,为后期TRAIL蛋白调控水牛卵巢卵泡发育、颗粒细胞增殖及凋亡的研究奠定基础。【方法】利用RT-PCR方法克隆水牛TRAIL基因CDS序列,对所获序列进行核苷酸序列、氨基酸序列相似性比对,构建系统进化树,并通过生物信息学软件分析TRAIL基因编码蛋白的结构和功能。【结果】试验成功克隆水牛TRAIL基因CDS序列,长864 bp,编码287个氨基酸;水牛TRAIL基因与牦牛、普通牛、山羊、绵羊、野猪、马、人、黑猩猩和家鼠的核苷酸序列相似性分别为99.2%、99.3%、95.9%、96.3%、84.7%、84.8%、81.3%、81.3%和70.0%。系统进化树结果表明,水牛与牦牛、普通牛的亲缘关系最近,与家鼠亲缘关系最远。氨基酸序列比对结果表明,在不同物种间,其跨膜结构域和TNF结构域序列保守性较高。TRAIL蛋白属于亲水性蛋白,存在1个跨膜结构域,140―285位氨基酸处为TNF区,具有29个磷酸化位点,无信号肽和糖基化位点,主要定位于细胞质中。TRAIL蛋白二级结构主要以无规则卷曲为主,约占51.57%,其次为延伸链(24.39%)和α-螺旋(24.04%)。TRAIL蛋白三级结构与二级结构一致,且与模型蛋白人TRAIL蛋白的相似性为75.53%。【结论】本试验克隆得到水牛TRAIL基因CDS区序列,大小为864 bp,编码287个氨基酸,水牛与牦牛、普通牛亲缘关系最近,TRAIL蛋白跨膜结构域和TNF结构域在不同物种间序列保守性较高,这可能与其功能有关。     

利用微卫星遗传标记评估5个西藏牦牛群体遗传多样性及群体结构

【目的】 评估当前西藏主要牦牛群体遗传多样性，梳理不同地区群体间遗传结构，明确西藏5个牦牛群体（阿里牦牛、斯布牦牛、娘亚牦牛、类乌齐牦牛和帕里牦牛）的保种情况和种群间系统发育关系。【方法】 利用13个微卫星标记（SSR）对5个牦牛群体共计195个个体进行基因分型，并对各群体的等位基因数量、基因杂合度、多态信息含量及群体间遗传距离等遗传参数进行评估。【结果】 阿里牦牛群体等位基因数最多（6.43），类乌齐牦牛等位基因数最少（5.00）；观测杂合度范围为0.5311（娘亚牦牛）~0.5995（类乌齐牦牛）。各群体内偏离Hardy-Weinberg平衡的位点数为5（类乌齐牦牛）~9（阿里牦牛）个；群体内近交系数最高为0.172（阿里牦牛），且4个群体（阿里牦牛、娘亚牦牛、斯布牦牛和帕里牦牛）存在显著近交风险（P<0.05）。从遗传结构来看，所有群体间均为显著遗传分歧（P<0.05），STRUCTURE分析结果显示，5个牦牛群体划分为3个簇，其中阿里牦牛较其他牦牛群体具有更为丰富的遗传背景。系统发育树结果表明，不同群体间系统发育关系相对独立，且与种群栖息地分布不一致。主坐标分析（PCoA）结果显示，不同群体间部分个体存在较近亲缘关系，表明不同群体间存在遗传物质交流。【结论】 5个西藏牦牛群体具有丰富的遗传多样性水平，且种群系统发育关系相对独立，但多数群体存在群体事件风险。本研究不仅有助于明确西藏牦牛地方种群遗传多样性水平，同时可为今后的保种策略提升提供理论参考。     

蒙古马胎儿成纤维细胞永生化细胞体系建立的研究

【目的】试验旨在建立永生化蒙古马胎儿成纤维细胞(equine fetal fibroblasts, EFFs)系,为马繁殖与发育生物学研究提供参考。【方法】采用1 mg/mL胶原酶Ⅳ消化分离原代EFFs,利用表达人端粒酶逆转录亚基(hTERT)和潮霉素B抗性的pLV-hTERT-IRES-hygro慢病毒感染EFFs,经潮霉素B筛选得到阳性细胞hTERT-EFFs,以第3代EFFs为对照,采用实时荧光定量PCR检测hTERT mRNA的表达水平,CCK-8法绘制生长曲线检测细胞增殖能力,流式细胞仪检测细胞凋亡情况,β-半乳糖苷酶染色检测细胞衰老情况。使用DNAMAN软件比对分析人、马、小鼠、牛、绵羊、猪多物种TERT核苷酸、蛋白和端粒酶RNA组分(TERC)序列相似性。【结果】原代EFFs在实验室培养体系下可传至17代,但第17代细胞发生了衰老。hTERT基因可成功转入EFFs,并至少稳定表达20代。然而hTERT-EFFs增殖能力显著低于EFFs(P<0.05),hTERT-EFFs凋亡早期细胞率显著高于EFFs(P<0.05)。细胞衰老检测结果表明,hTERT-EFFs...     

关中奶山羊FGF2基因克隆、生物信息学分析及组织表达研究

【目的】对关中奶山羊成纤维细胞生长因子2（fibroblast growth factor 2，FGF2）基因进行克隆及生物信息学分析，并检测其在山羊各组织中的表达差异，为后续探究该基因的功能奠定基础。【方法】以关中奶山羊乳腺cDNA为模板，采用RT-PCR技术扩增并克隆关中奶山羊FGF2基因完整CDS区序列，进行相似性比对、系统发育树构建及生物信息学分析；运用实时荧光定量PCR方法检测FGF2基因在关中奶山羊心脏、肝脏、脾脏、肾脏、背最长肌、肺脏、乳腺和卵巢组织中的表达情况。【结果】关中奶山羊FGF2基因编码区长468 bp，可编码155个氨基酸，分子质量为17.26 ku，理论等电点为9.58，其中含量最高的是甘氨酸，占10.3%。相似性比对结果表明，关中奶山羊FGF2氨基酸序列与人、牛、马、兔、绵羊、虎鲸和野猪的相似性分别为98.7%、99.4%、100.0%、98.7%、99.4%、99.4%和62.6%。系统进化树显示，FGF2基因在不同物种间具有高度保守性，与马的亲缘关系最近。关中奶山羊FGF2蛋白不稳定指数为38.39，属于稳定亲水性蛋白质，不含跨膜结构和信号肽；FGF2蛋白二级结构含有α-螺旋、β-转角、延伸链、无规则卷曲，分别占比10.32%、14.19%、30.97%、44.52%。蛋白质互作预测分析显示，FGF2蛋白与FGFR、KDR、FGF1、FGFR4、MET和FGFR1等与乳腺生长发育相关的蛋白存在相互作用。实时荧光定量PCR检测到关中奶山羊FGF2基因在心脏、肝脏、脾脏、肾脏、背最长肌、肺脏、乳腺和卵巢中均有表达，在乳腺中表达水平最高，其次为卵巢，在背最长肌中的表达水平最低。【结论】关中奶山羊FGF2基因CDS区序列长468 bp，编码155个氨基酸，为亲水蛋白，二级结构以延伸链和无规则卷曲为主。FGF2基因在关中奶山羊泌乳高峰期的不同组织中均有表达。本试验结果为进一步探究FGF2基因在关中奶山羊乳腺发育中的作用及其具体功能提供了参考。     

藏羊MSTN基因克隆、生物信息学及组织表达分析

【目的】 对藏羊肌肉生长抑制素（myostatin，MSTN）基因进行克隆和生物信息学分析，检测其在藏羊不同组织中的表达，为探究MSTN基因在藏羊中的生物学功能提供参考。【方法】 以藏羊背最长肌组织cDNA为模板，克隆藏羊MSTN基因完整CDS区序列并测序，用SeqMan程序对测序结果进行拼接，并用BLAST在线程序对组装后的序列进行分析鉴定。用生物信息学软件进行相似性比对、系统进化树构建及生物信息学分析，用实时荧光定量PCR检测MSTN基因在藏羊不同组织中的表达量。【结果】 藏羊MSTN基因CDS全长为1 128 bp，编码375个氨基酸。藏羊MSTN基因氨基酸序列与绵羊、牦牛、牛、猪、恒河猴、人、黑猩猩、犬、鸡及斑马的相似性依次为100.0%、93.4%、93.4%、95.2%、94.4%、94.2%、94.4%、93.1%、87.8%和87.5%；系统进化树分析结果表明，藏羊与绵羊的亲缘关系最近，与斑马和鸡的亲缘关系最远。藏羊MSTN蛋白属于亲水性分泌蛋白，且具有不稳定性，不含跨膜结构，含1个信号肽，存在31个潜在的磷酸化位点、2个N-糖基化修饰位点，主要分布在线粒体和细胞质中；MSTN蛋白二级结构以无规卷曲为主，其次是α-螺旋、延伸链和β-转角；三级结构预测结果与二级结构一致。实时荧光定量PCR结果显示，MSTN基因在藏羊不同组织中均有表达，其中在臂三头肌和半腱肌的表达量显著高于肺脏、下丘脑、心脏、肝脏、十二指肠、瘤胃和肾脏（P<0.05）。【结论】 成功克隆了藏羊MSTN基因；该基因在藏羊肌肉中的表达量高于内脏组织。该结果为进一步研究MSTN基因对藏羊肌肉生长发育的调控机制奠定了基础。     

基于线粒体ATP6及Cytb基因多态性研究藏山羊高原适应性

为了探究藏山羊ATP6和Cytb是否具有潜在的高海拔适应性突变，本研究利用藏山羊（n=157，来自西藏地区的8个群体，海拔高度>3 500 m）与低海拔山羊（n=104，来自欧亚大陆50个群体，海拔高度<1 000 m）两组261条ATP6和Cytb序列，对核苷酸组成、遗传多样性、单倍型网络以及基因突变对蛋白结构域功能影响进行比较分析。结果，在ATP6和Cytb中分别鉴定出33和67个单核苷酸多态性位点（SNPs），其中错义突变均为6个。值得注意的是，ATP6 m.8102A>G（Ile→Met；P=0.000 6）和Cytb m.14794A>G（Thr→Ala；P=0.007 7）是低海拔山羊群体特有的错义突变，由此形成的单倍体H27和Ha36与高原适应性呈显著负相关（H27，P=0.007 0；Ha36，P=0.012 3）。进一步分析推测，这两个突变阻碍了质子跨膜的正常功能，并影响了质子或电子在OXPHOS中的转移，从而导致了低海拔山羊和藏山羊之间呼吸作用效率的差异。本研究通过对线粒体ATP6和Cytb基因编码区多态性到功能的相关分析，初步推测了藏山羊高海拔低氧适应性机制，为高原适应性相关研究提供了一定的见解。     

绵羊MyoG基因内含子Ⅱ多态性及其与生长性状的关联分析

【目的】 分析肌细胞生成素(myogenin,MyoG)基因内含子Ⅱ的多态性及其对绵羊生长性状的影响,筛选对绵羊生长发育有显著影响的分子标记,以期为绵羊的分子育种提供参考。【方法】 选取大尾寒羊、小尾寒羊、豫西脂尾羊、湖羊、杜泊羊5种绵羊为研究对象,采用PCR-RFLP技术对MyoG基因内含子Ⅱ(Eco72 Ⅰ)进行基因分型,利用PopGene32软件计算绵羊MyoG基因内含子Ⅱ的群体遗传多样性,利用SPSS 17.0软件对不同基因型与绵羊生长性状进行关联分析。【结果】 小尾寒羊、杜泊羊、湖羊群体中MyoG基因内含子Ⅱ均存在3种基因型:AA(368/540 bp)、AB(908/368/540 bp)和BB(908 bp);大尾寒羊和豫西脂尾羊群体中仅检测到2种基因型:AB(908/368/540 bp)和BB(908 bp)。大尾寒羊、小尾寒羊、豫西脂尾羊、湖羊、杜泊羊的AB基因型频率分别为0.845、0.633、0.917、0.706和0.811,BB基因型频率分别为0.155、0.033、0.083、0.176和0.054,AA基因型频率分别为0、0.333、0、0.118和0.135。小尾寒羊的杂合度最低(0.455),杜泊羊的杂合度最高(0.497),表明杜泊羊的遗传多样性要高于其他群体;5个群体多态信息含量(PIC)均为中度多态(0.25<PIC<0.5)。关联分析发现,小尾寒羊MyoG基因内含子Ⅱ BB基因型个体胸围、胸宽、头长、颈长均显著低于AA、AB基因型(P<0.05)。【结论】 MyoG基因内含子Ⅱ Eco72 Ⅰ位点可作为影响绵羊生长性状的分子标记,结果可为今后绵羊生长性状的分子标记辅助选择提供参考。     

基于RAD-Seq技术的大围山微型鸡遗传进化分析

【目的】 在分子水平上探讨大围山微型鸡的遗传进化。【方法】 以大围山微型鸡与其他18个地方鸡品种及2个国外引进品种为研究对象,每个品种选取公鸡10只,母鸡20只,采血,提取全基因组DNA,进行简化基因组测序(RAD-Seq),鉴定21个品种基因组SNP,计算大围山微型鸡遗传统计量,分析遗传多样性和遗传结构、筛选受选择基因并进行功能富集分析。【结果】 在大围山微型鸡群体中鉴定出331 892个SNPs。大围山微型鸡的平均杂合度(Ho)为0.219、核苷酸多样度(Pi)为0.245,近交系数(Fis)为0.107,与其他18个地方鸡种相比遗传多样性呈中度多态。聚类分析发现,大围山微型鸡与瓢鸡、茶花鸡和藏鸡聚为一类,亲缘关系较近,且与瓢鸡的群体分化指数(Fst)最低(0.0929),亲缘关系最近,与河南斗鸡的Fst最高(0.2179),亲缘关系最远。共筛选出200个受选择基因。GO分析结果显示,这些受选择基因主要富集在运动、刺激应答、信号传导等生物学过程;KEGG分析结果表明,这些受选择基因主要富集在代谢、肌动蛋白细胞骨架的调节、MAPK等信号通路。【结论】 大围山微型鸡遗传多样性呈中度多态,与瓢鸡亲缘关系最近,本研究筛选出了200个受到选择的基因,这些基因在代谢、信号传导、颅骨发育等多个方面发挥作用。     

Helminths of wild boar in the isolated population close to the northern border of its habitat area

One hundred wild boars (Sus scrofa) from a geographically isolated population on the island Saaremaa of western Estonia were examined for visceral helminths. Seven helminth species, Metastrongylus pudendotectus, M. salmi, M. elongatus, Ascaris suum, Trichuris suis, Dicrocoelium dendriticum and Taenia hydatigena larva, were found. The predominant helminths discovered were lung nematodes (prevalence 82%, mean intensity 96.2 per animal). A significant negative correlation was observed between the weight of wild boars and the number of lungworms and helminth species. The number of helminth species found in the wild boar population on the island was lower compared to that of the adjacent mainland.     

浙江6个地方猪种的mtDNA控制区多态性分析

本研究对淳安花猪等6个浙江地方猪种的116条mtDNA序列进行了多态性分析,参照猪群为长白猪和皖南花猪。结果显示,浙江地方猪种581 bp mtDNA控制区序列中共存在14个变异位点,群体的核苷酸多样度Pi值为0.00403±0.00036;单倍型多样度Hd值为0.835±0.018。地方猪的单倍型从Hap_5至Hap_19共15种,而长白猪为Hap_1至Hap_4。淳安花猪和皖南花猪的单倍型频率均以Hap_8为主(0.680和0.833);嘉兴黑猪以Hap_5(0.625)为主;嵊县花猪和碧湖猪均以Hap_6为主(0.500和0.375);岔路黑猪以Hap_5、Hap_8居多(0.429和0.429);金华猪有11种单倍型,频率最高的是Hap_5(0.241)。NJ系统发生树中,浙江地方猪种大致可分为5支(编号A~E),其中,淳安花猪主要出现在E支,个体数占淳安花猪样本数的68.0%(17/25);C支全为金华猪,个体数占金华猪样本数的29.3%(17/58);D支金华猪个体数占金华猪样本数的24.1%(14/58)。本研究结果表明,浙江地方猪种的mtDNA多态性较为丰富,且与长白猪在基因序列上存在明显差别。地方猪种的单倍型组成和频率与其地理分布及品种的形成特点有密切关系。     

基于线粒体控制区的鸡不同杂交组合遗传多样性研究

旨在探讨鸡不同杂交组合线粒体控制区（mtDNA D-loop区）的遗传多样性和单倍型特性。选取固始鸡和隐性白羽鸡及其正、反交F1代、藏鸡以及F2代等6个群体共387个个体的mtDNA D-loop区进行测序,分析其遗传规律和单倍型特性,并与不同红色原鸡亚种进行聚类,分析其母系起源。结果显示,6个群体D-loop区全序列大小为1 231 bp,共检测到28个多态位点和1个C碱基缺失,共构成19种单倍型,分为A、B、C和E 4个单倍型群,其中,固始鸡和反交F1代主要为A、C单倍型,固始鸡A、C单倍型比例分别为53.42%和46.58%,反交F1代A、C单倍型比例分别为50.75%和49.25%;隐性白羽鸡、正交F1代和F2代优势单倍型均为E单倍型,占比分别为48.89%、48.84%和50.00%。6个鸡群体单倍型多样度（Hd）在0.496~0.729之间,核苷酸多样度（Pi）在0.003 40~0.005 41之间,Hd值和Pi值最大的均为正交F1代,其次为隐性白羽鸡和F2代,固始鸡和反交F1代群体遗传多样性接近。聚类分析显示,A、B单倍型群与滇南亚种交叉聚为一枝;E单倍型群与印度亚种交叉聚为一枝;C单倍型群与印度亚种、指名亚种、印尼亚种以及滇南亚种聚为一枝。结果提示,mtDNA D-loop区遵循严格的母系遗传,后代的遗传多样性和单倍型比例与其母本基本一致;我国家鸡群体具有多个红色原鸡母系起源,且主要起源于原鸡滇南亚种。     

中畜草原白羽肉鸭与樱桃谷鸭的遗传差异分析

【目的】 通过基因组重测序技术对中畜草原白羽肉鸭与樱桃谷鸭的遗传差异进行分析，追溯两个品种鸭在不同人工选择下的基因组变异机制，以此阐明优异性状形成的遗传基础。【方法】 选择樱桃谷鸭商品代和中畜草原白羽肉鸭商品代各16只进行基因组重测序，过滤掉高缺失率与最小等位基因频率较低的位点，获得高质量SNPs用于后续分析；对两品种的基因型数据进行主成分分析（PCA）确定其遗传分化情况；采用群体遗传分化指数（Fst）和群体核酸多样性比值（Pi）两种分析方法综合筛选中畜草原白羽肉鸭和樱桃谷鸭的受选择信号。【结果】 主成分分析结果显示，中畜草原白羽肉鸭和樱桃谷鸭分化显著。以10 kb窗口5 kb步长分别计算Fst与Pi分析值，取前1%作为阈值（Fst>0.177，Pi>0.885），在两种分析方法的信号重叠区域共筛选到410 kb候选区域，共注释到21个候选基因。对候选基因进行GO与KEGG富集分析发现，12个基因显著富集到细胞组分和分子功能两大类中（P<0.05），2个基因显著富集到代谢相关通路（P<0.05）。这些基因中，与脂质代谢、氨基酸代谢、免疫调控相关的基因包括PDE3A、PRKAR2B、SEMA5A、SHANK2、STXBP6与LOC101803508（GOLGB1）受到了强选择。【结论】 虽然樱桃谷鸭与中畜草原白羽肉鸭均源自北京鸭，但经过不同强度、不同方向持续的人工选育，2个品种明显分化，且中畜草原白羽肉鸭具有更高的遗传多态性。筛选到了2个品种间一系列遗传分化候选基因，并重点挖掘了参与白羽肉鸭风味肉品质调控的相关基因，为后续研究不同白羽肉鸭品种特征、筛选品种特异性分子标记提供了参考。     

Investigation of Co-infection of PPV7 and PCV2 and Genetic Variation Analysis for PPV7 Cap Gene in Fujian and Guangdong

【Objective】 To investigate the co-infection situation of Porcine parvovirus (PPV7) and Porcine circovirus type 2 (PCV2)in Fujian and Guangdong, and to understand the molecular genetic characteristics of PPV7 Cap gene.【Method】 The blood sample of 432 infected pigs from 69 pig farms in Fujian and Guangdong were collected to detect PPV7 and PCV2 by PCR.The PPV7 Cap gene of positive samples was cloned and sequenced.The DNAStar software was used to analyze the nucleotide and amino acid sequences of PPV7 Cap gene, and Mega 7.0 software was used to draw the genetic evolution tree.【Result】 The results showed that the positive rate of PPV7 was 21.99% (96/432), the positive rate of farms was 53.62%(37/69), and the positive rate of PCV2 was 54.17% (234/432), and the co-infection rate of PCV2 and PPV7 was 13.43% (58/432).The 17 PPV7 Cap gene sequences were amplified using PCR.Nucleotide homology analysis revealed that the homology of the 17 PPV7 Cap gene sequences was 85.6%-100%, and the homology with reference strains was 85.8%-99.0%.Amino acid sequence comparison analysis revealed that the amino acid homology of the 17 PPV7 Cap protein sequences was 87.6%-100%, and the amino acid homology with reference strains was 82.6%-98.7%.Phylogenetic analysis of Cap gene showed that PPV7 could be divided into five main evolutionary branches of PPV7a-PPV7e, among which 9 isolates belonged to PPV7a subtype, 3 isolates belonged to PPV7b subtype, 4 isolates belong to PPV7c subtype, and only 1 isolates isolates belonged to PPV7e subtype.【Conclusion】 This study indicated that PPV7 was widely prevalent in Fujian and Guangdong regions, and had a high co-infection rate with PCV2, which might be the pathogenic factor of Porcine circovirus associated disease(PCVAD).The genetic diversity of PPV7 isolates was abudant in both regions, and PPV7a was the dominant strain at present.The findings of this study provided theoretical basis and data reference for PPV7 prevention and control and vaccine research.     

中国荷斯坦牛PIK3CB基因多态性及其与繁殖和产奶性状的关联分析

【目的】 探究磷脂酰肌醇-4,5-二磷酸3-激酶催化亚单位β（PIK3CB）基因多态性及其与中国荷斯坦牛繁殖和产奶性状的关系。【方法】 通过混池测序对中国荷斯坦牛PIK3CB基因进行单核苷酸多态性（SNP）位点筛选,采用竞争性等位基因特异性PCR （KASP）技术在1 160头健康泌乳中国荷斯坦牛中进行SNP分型并进行群体遗传学分析,采用线性模型进行SNP与11个繁殖和产奶性状基于单位点和单倍型组合的关联分析。【结果】 在PIK3CB基因中共检测到了17个SNPs,筛选出7个SNPs用于后续分析。关联分析发现,7个SNPs与多个目标性状存在显著或极显著的关联（P<0.05;P<0.01）;位于外显子区域的g.130433743 A>G位点AA基因型个体和位于可变剪接区域的g.130448069 G>A位点GG基因型个体,其经产牛首末次配种间隔、产奶量、乳蛋白量和乳脂量最低,体细胞评分最高,上述基因型个体具有较短的首末次配种间隔,而产奶性能相对较差;g.130387717 G>A位点AA基因型个体,其初配日龄和青年牛首末次配种间隔最低,产奶量、乳蛋白量和乳脂量最高,该基因型个体的繁殖和产奶性能均较好,上述3个SNPs位点可作为中国荷斯坦牛繁殖和产奶性状的候选位点重点关注。单倍型分析发现,PIK3CB基因的g.130387717 G>A、g.130430832 A>－、g.130433743 A>G、g.130433982 C>T、g.130446073 C>T和g.130448069 G>A 6个SNPs紧密连锁形成一个单倍型块,且与多个目标性状存在显著或极显著关联（P<0.05;P<0.01）,其中H2H3和H2H4单倍型组合个体的繁殖和产奶性能较好,为优势单倍型组合。【结论】 中国荷斯坦牛PIK3CB基因存在丰富的遗传变异,其多态性与繁殖和产奶性状存在关联,g.130433743 A>G、g.130448069 G>A和g.130387717 G>A位点可作为潜在分子标记,为中国荷斯坦牛的平衡育种提供理论依据。     

ZBED6基因敲除猪肾脏转录组分析

【目的】 探究锌指BED结构域结合蛋白6(zinc finger BED domain-containing protein 6,ZBED6)基因敲除对猪肾脏组织基因转录表达的影响,并解析ZBED6基因调控猪肾脏代谢的靶基因及其通路。【方法】 对8月龄ZBED6基因敲除巴马香猪(ZBED6 KO)和同月龄野生型巴马香猪(WT)的肾脏组织(n=3)进行总RNA提取,利用实时荧光定量PCR检测胰岛素样生长因子2(insulin-like growth factor 2,IGF2)和差异基因的表达量;以Illumina Hiseq高通量测序技术对ZBED6 KO和WT猪肾脏组织mRNA进行转录组测序(RNA-Seq),用猪Sus scrofa 11.1作为参考基因组序列,通过生物信息学软件TopHat、Cufflink、Cuffmerge和Cuffdiff对测序数据进行分析,筛选ZBED6 KO和WT猪肾脏组织中的差异表达基因,对差异表达基因进行层次聚类和KEGG通路富集分析,并通过实时荧光定量PCR验证RNA-Seq结果中差异表达基因的可靠性。【结果】 实时荧光定量PCR结果显示,ZBED6 KO猪肾脏IGF2基因mRNA表达量显著高于WT猪(P<0.05)。测序共获得78 G数据量,每个样本比对率在87.3%以上,测序质量良好;ZBED6 KO和WT猪肾脏组织之间共检测到25 213个基因,以调整后P<0.05和log₂|FoldChange|>2为筛选标准,得到299个差异表达基因,其中上调的基因103个,下调的基因199个。热图和主成分分析(PCA)结果显示,ZBED6 KO和WT组猪肾脏基因组内表达模式相似,组间表达模式不同;KEGG通路富集分析显示,差异基因主要参与视黄醇及疾病代谢相关通路。ZBED6基因敲除后,其中有9个差异表达基因(CYP2C42、AOX1、ENSSSCG00000036274、RDH16、CYP2A19、ENSSSCG00000022724、CYP26B1、CYP1A1、RDH5)被富集到视黄醇代谢通路中,可能参与调控猪肾脏代谢平衡。实时荧光定量PCR检测发现,7个差异表达基因(CYP2C42、AOX1、RDH16、CYP2A19、CYP26B1、CYP1A1、RDH5)的表达模式与RNA-Seq分析结果一致,证实RNA-Seq结果的可靠性。【结论】 本试验利用RNA-Seq技术分析了ZBED6基因敲除对巴马香猪肾脏代谢的影响,其通过调控肾脏多个下游基因表达,从而影响其代谢相关信号通路,试验结果为阐明ZBED6基因功能提供了材料。     

梅花鹿MSRB3基因多态性及其与茸重性状的关联分析

【目的】 本研究旨在探索梅花鹿甲硫氨酸亚砜还原酶B3（methionine sulfoxide reductase B3，MSRB3）基因多态性及其与茸重性状的关联性。【方法】 选取高、低产茸量的梅花鹿各8只，应用DNA直接测序法检测基因变异位点。采用MassARRAY^® SNP分型技术对314头24月龄梅花鹿MSRB3基因进行基因分型，并结合梅花鹿茸重性状数据进行了关联分析和单倍型分析。【结果】 在梅花鹿MSRB3基因中共发现6个SNPs位点，其中2个SNPs位点位于外显子区域，且突变均未引起氨基酸改变，属于同义突变，其余4个SNPs位点均存在于内含子区域。分型结果显示4个样本未分型成功，其余310个样本进行后续分析。各位点观测杂合度和期望杂合度基本一致，g.44455582 T>C、g.44455759 C>T、g.44414424 T>C、g.44350306 T>C及g.44340836 G>A等5个位点的杂合度较高，均属于中度多态（0.25<PIC<0.5），g.44340734 G>C位点杂合度较低，属于低度多态（PIC<0.25）。卡方检验结果表明，g.44455759 C>T位点偏离Hardy-Weinberg平衡状态（P<0.05），其他5个位点均处于Hardy-Weinberg平衡状态（P>0.05）。对6个SNPs位点的不同基因型与梅花鹿茸重的关联分析结果表明，g.44455582 T>C位点的TT基因型个体茸重极显著高于CC和TC基因型（P<0.01）。单倍型分析结果显示，g.44340734 G>C和g.44350306 T>C、g.44455582 T>C和g.44455759 C>T位点存在强连锁，各产生3种单倍型，在Block 2区块中单倍型TC茸重显著高于单倍型CT （P<0.05）。【结论】 MSRB3基因与梅花鹿茸重性状密切相关，g.44455582 T>C位点和单倍型TC可作为筛选高产梅花鹿的分子标记。