1株猪源动物联合乳杆菌S7全基因组测序及生物信息学分析

[目的]旨在对1株具有较强耐酸耐胆盐特性的猪源动物联合乳杆菌S7进行全基因组测序及生物信息学分析。[方法]结合三代PacBio RSⅡ和二代Illumina HiSeq 2000测序技术对动物联合乳杆菌S7进行了全基因组测序,在此基础上使用GO、EggNOG、KEGG、CAZy、VFDB和CARD等数据库注释功能基因,利用TYGS(Type Strain Genome Server)构建系统发育树,并与2株同源模式菌株进行了共线性分析。[结果]动物联合乳杆菌S7的基因组大小为2.03 Mb, GC含量为44.14%,共预测到1 993个编码基因,包含65个tRNA、19个rRNA、35个sRNA、29个持家基因、11个基因组岛、6个前噬菌体、4个CRISPR-Cas、8个插入序列和10个转座子;分别有1 521、1 567、1 185和60个基因在GO、EggNOG、KEGG和CAZy数据库中得到注释;在VFDB和CARD数据库中注释到181个毒力基因和110个耐药基因;另外有14个基因参与耐酸,2个基因参与耐胆盐,22个基因参与黏附和聚集;基因组中还包括与温度应激、氧化应激和细菌素合成等多种与益生功能相关的基因。系统发育树和共线性分析发现,该菌株与模式菌株Ligilactobacillus animalis P38的进化关系最近。[结论]动物联合乳杆菌S7基因组中具有与耐酸、耐胆盐、黏附和聚集、温度应激、氧化应激和细菌素等相关的功能基因,可为菌株在饲料添加剂中的科学应用提供理论依据。

葡萄牙栖盐田菌LLJ914的全基因组测序和比较基因组分析

栖盐田菌属细菌因具有较强渗透压耐受能力而备受关注。然而,目前对其适应高盐或其他海洋极端环境的机制尚未探明。为揭示分离自冲绳海槽热液区的葡萄牙栖盐田菌(Salinicola lusitanus)LLJ914与其同属菌株之间的遗传特征及代谢潜力的差异,探究其在极端海洋生境下的适应机制。通过全基因组测序获得葡萄牙栖盐田菌LLJ914全基因组序列,选取同属其他菌株的基因组进行比较基因组学分析,探究栖盐田菌属各菌株及不同环境来源的葡萄牙栖盐田菌代谢潜力的差异。菌株LLJ914基因组大小为4781556 bp,GC含量为64.0%,共编码4229个蛋白,69个tRNA,12个rRNA。通过构建系统发育树,并进行平均核苷酸和平均氨基酸一致性分析,发现该菌株与马齿苋(Halimione portulacoides)内共生的葡萄牙栖盐田菌CR50^(T)亲缘关系最近。通过功能基因注释分析,发现葡萄牙栖盐田菌具有抵抗各种重金属的相关基因和利用甲基膦酸产生甲烷的phn基因簇;相比于植物内共生菌CR50^(T),菌株LLJ914基因组中包含更多与氨基酸和碳水化合物的运输代谢、能量生产和转换以及转录相关的功能基因。此外,菌株LLJ914基因组中还包含特有的与重金属抵抗、免疫防御及有氧呼吸相关的基因,这可能与其适应复杂极端的深海热液环境有关。本研究揭示了葡萄牙栖盐田菌LLJ914适应深海热液环境的遗传特征及代谢潜力,为更好地认识热液微生物的生态功能提供了参考。

苗草专用型大麦品种鉴定及全基因组关联分析

苗草工厂为实现草食动物的饲草周年供应提供了新的解决方案。针对苗草生产的品种需求,本研究对124份中国大麦育成品种(系)与种质开展了苗草转化率鉴定和生物量调控位点发掘。结果显示,大麦种子萌动后,水培条件下苗草生物量呈指数增长趋势并在7 d后进入平台期。在植物工厂条件下,鉴定出冬青16、扎青6号等10个高苗草转化率品种且分析发现苗草生物量与籽粒千粒重呈一定的负相关。进一步通过全基因组关联分析,共鉴定到12个苗草生物量相关QTL位点,从中预测出8个调控苗草生物量的候选基因。本研究不仅为大麦苗草工厂化生产筛选出高转化率品种,同时也为苗草专用型大麦品种的遗传改良奠定了基础。

蓝莓转录因子VcICE全基因组鉴定及其在花果发育进程的表达分析

为探究ICE(inducer of CBF expression)转录因子在高丛蓝莓花果发育进程中可能的调控机制,通过对VcICE基因家族进行全基因组鉴定及生物信息学分析,采用实时荧光定量PCR(qRT-PCR)技术分析VcICE基因亚家族成员在蓝莓花芽膨大和果实发育进程中的相对表达水平.结果表明:在高丛蓝莓‘Draper’基因组中鉴定出8个VcICE成员,系统进化分析将其划分为2个亚家族,各亚家族成员均包含S-rich,bHLH和ACT等保守结构域;VcICE基因启动子区域富含光响应、植物激素、生长发育及非生物胁迫相关顺式作用元件;VcICE亚家族基因在高丛蓝莓‘O’Neal’和‘Bluerain’花芽膨大与果实发育进程中的相对表达差异显著.研究结果以期为深入了解VcICE转录因子在蓝莓花果发育进程中的功能提供理论参考.

大豆叶型性状全基因组关联分析与候选基因鉴定

大豆叶片的形状和垂直分布影响群体冠层结构和光合效率并最终影响大豆的产量。植物叶片大小、形态因着生位置不同而产生差异的现象称为异形叶,虽然异形叶现象在被子植物中广泛存在,但目前关于大豆叶片发育过程中异形叶的调控的研究还很有限。本研究通过对283份大豆种质资源的叶长、叶宽、叶形指数和异形叶指数等叶型相关的性状在江苏南京进行连续2年的考察,利用全基因组关联分析检测到181个叶型性状相关位点,其中能够在2个环境或多个性状中重复检测到的位点18个。利用检测到的与叶型相关的SNP位点,结合基因的表达谱数据、拟南芥中同源基因的功能,鉴定与大豆叶片发育和异形叶形成相关的候选基因。其中在20号染色体的相关位点Chr20:36152820上游发现已知的大豆叶形调控基因Ln(Glyma.20G116200)。此外,在19号染色体的相关位点Chr19:45155943附近鉴定到2个候选基因Glyma.19G192700、Glyma.19G194100,分别被注释为Growth-regulating factor 4(GRF4)和LITTLE ZIPPER 3(ZPR3)基因的同源基因,为阐明大豆异形叶等叶型性状遗传的分子机制奠定了基础。

大豆地方种质资源鲜籽粒可溶性糖含量的全基因组关联分析

【目的】可溶性糖含量是鲜食大豆重要的品质性状之一,研究大豆鲜籽粒可溶性糖含量的遗传变异,深入解析可溶性糖含量的遗传机制,为鲜食大豆种质创新及品质育种提供依据。【方法】利用来自东北大豆生态区、北方大豆生态区、黄淮海大豆生态区和南方大豆生态区的133份大豆地方种质,在2021年连江春季、福清春季和秋季3个环境下对鲜籽粒可溶性糖含量进行表型测定,结合82187个高质量SNP标记,基于混合线性模型MLM(Q+K)对可溶性糖含量进行全基因组关联分析,挖掘可溶性糖含量显著相关的SNP位点,并以显著SNP位点为中心,两端各扩展119.07 kb连锁不平衡衰减距离为候选区间,根据候选区间内基因的注释和组织表达信息预测候选基因。【结果】3个环境下,鲜籽粒可溶性糖含量的变异范围为3.37—33.84 mg·g^(-1),遗传变异系数为24.59%—32.69%,可溶性糖含量遗传率为68.14%。通过全基因组关联分析,连江春季、福清春季和秋季3个环境下分别检测到与鲜籽粒可溶性糖含量显著关联的SNP有6、8和22个,表型变异解释率为12.43%—29.27%,以表型变异解释率较高的9个显著SNP位点所在的候选区间进行搜索,共获得86个基因,结合基因注释和组织表达信息,进一步筛选到9个候选基因,主要涉及转录因子、糖蛋白家族和糖类合成转运等生物学过程。其中,Glyma.01g016500、Glyma.13g042100、Glyma.16g131800和Glyma.16g155300在大豆种子及荚中表达水平较高,可作为大豆鲜籽粒可溶性糖的最具潜力候选基因。【结论】通过全基因组关联分析,检测到36个与鲜籽粒可溶性糖含量显著关联的SNP,进一步筛选出9个候选基因可能参与大豆鲜籽粒可溶性糖含量的调控,其中,Glyma.01g016500、Glyma.13g042100、Glyma.16g131800和Glyma.16g155300可作为调控大豆鲜籽粒可溶性糖含量的关键目标基因。

基于全基因组重测序数据的新疆褐牛基因组结构变异检测及群体结构分析

旨在鉴定新疆褐牛基因组结构变异(structure variation,SVs)特征,并在此基础上探索新疆褐牛培育过程中受到影响的结构变异。本研究基于全基因组重测序数据,利用Manta、Delly、Lumpy三款软件对新疆褐牛及其父母本(瑞士褐牛、哈萨克牛)以及中国西门塔尔牛基因组结构变异进行检测,之后利用主成分分析、构建系统发育树和遗传分化指数(F_(ST))将新疆褐牛基因组与其余3个品种进行比较分析。结果显示,4个牛品种共检测出54969个SVs,缺失型变异占比最高(56.59%),插入型变异占比最少(0.03%)。从数量上看,这些SVs在染色体上的分布呈现出随染色体号变大而逐渐减少的趋势。不同品种间存在共有变异和品种特有变异,其中地方品种哈萨克牛拥有的特有SVs数量最多。主成分分析和系统发育树结果发现,新疆褐牛同哈萨克牛和瑞士褐牛具有较近的亲缘关系。经选择信号分析、基因注释和富集分析,挖掘出了一批与产奶性状(PI 4K2A、ELOVL3、ECHS1、SCD、TCF 7L2、PNLIPRP2、BTRC、PLCE 1)、生长发育(BMP6、TLL2、MAPK9、ROR 2)、适应性(PRKC B)以及免疫(GSTO2、GSTO 1)相关的关键基因。本研究从结构变异上解析新疆褐牛的特征,并揭示一些新疆褐牛培育过程中高度分化的基因,为推动新疆褐牛的遗传改良提供了基础资料。

全基因组测序分析生畜肉中沙门氏菌血清型、耐药性与毒力因子

了解市售生畜肉中沙门氏菌(Salmonella)血清型、序列型(sequence typing,ST)分布、耐药性和毒力因子特征,探讨ST与耐药性和毒力的关系。收集2020—2023年从生畜肉中分离的68株沙门氏菌,用血清凝集法鉴定其血清型,微量肉汤稀释法进行药物敏感检测;全基因组测序分析多位点序列分型、耐药基因与毒力因子,并结合系统发育树分析其相关性。结果表明:68株沙门氏菌可分为14种血清型,以鼠伤寒沙门氏菌和肠炎沙门氏菌为主,ST以ST19、ST34和ST11为主;多重耐药菌株占比69.12%(47/68),对氨苄西林(76.47%)、四环素(55.88%)和萘啶酸(52.94%)耐药率较高;68株沙门氏菌均携带菌毛吸附相关因子、镁离子转运相关因子、非菌毛类吸附相关因子、调控相关因子和分泌系统相关因子;生畜肉中污染的沙门氏菌血清型种类多样,耐药性严重且携带多种毒力因子,需要加强畜类养殖和生产销售的管理。

重瓣榆叶梅全叶绿体基因组遗传特征分析

【目的】对重瓣榆叶梅Prunus triloba‘Multiplex’全叶绿体基因组序列进行测序,探究其系统发育位置并分析其叶绿体基因组成特点。【方法】以重瓣榆叶梅叶片为材料,采用2×CTAB法提取叶绿体DNA,利用Illumina NovaSeq平台进行叶绿体基因组的测序,组装、注释并分析其叶绿体基因组遗传特征。联合美国国家生物技术信息中心(NCBI)数据库数据,基于全叶绿体基因组序列构建了重瓣榆叶梅系统进化关系。【结果】重瓣榆叶梅叶绿体基因组全长为157827 bp,NCBI登录号MT937181,其结构为经典的四分体结构,由1个大单拷贝区域(LSC),1个小单拷贝区域(SSC)及反向重复区域(IRa/IRb)构成,其序列长度分别为86032、19023、26386 bp。GC和AT的总占比分别为36.80%和63.20%。重瓣榆叶梅的完整叶绿体基因组序列注释到132个基因,包括tRNA基因、编码蛋白基因、rRNA基因,分别为37、87、8个。重瓣榆叶梅叶绿体基因组共编码26678个密码子和236个符合条件的SSR位点。SSR位点中A/T碱基占优势,碱基偏好性十分明显。【结论】系统进化树分析表明,重瓣榆叶梅和榆叶梅P.triloba聚合成一分支结构,与同属植物长柄扁桃P.pedunculata亲缘关系较近。

基于全基因组测序的耐多黏菌素和替加环素肺炎克雷伯菌特征研究

目的分析1株多黏菌素和替加环素均耐药的肺炎克雷伯菌基因组特征。方法肺炎克雷伯菌KP2016分离自我院临床痰液标本。采用微量肉汤稀释法测定KP2016对亚胺培南、美罗培南、多黏菌素和替加环素的最低抑菌浓度。对菌株进行第二代Illumina和第三代Oxford Nanopore全基因组测序。通过Unicycler对第二代和第三代序列进行混合拼接。通过ABRicate v0.8.13调用Resfinder数据库分析耐药基因。通过ISFinder分析移动元件。利用CGE(Center for Genomic Epidemiology)网站分析菌株的ST型和质粒复制子类型。利用Proksee对质粒结构进行可视化。结果KP2016菌株对多黏菌素和替加环素均耐药,MIC分别为512和16 mg/L。MLST分析显示KP2016属于ST656,携带多黏菌素耐药相关mcr-1和mcr-8基因和替加环素耐药相关tmexCD1-toprJ1基因簇。KP2016携带1个染色体和5个环形质粒,质粒大小3,991~275,345bp,GC含量44.35%~52.20%。mcr-1基因位于约44kb的质粒p1上,上下游环境为ISKpn26-mcr-1-pap2-ISApl1;mcr-8基因位于IncR/IncN型质粒p2上,遗传结构为ISEcl1-mcr-8-orf-ISKpn26;tmexCD1-toprJ1基因簇结构为IS26-tnfxB1-tmexC1-tmexcD1-toprJ1。此外,菌株还携带多黏菌素耐药相关的染色体介导的CrrB突变(L204V、V237I)。结论肺炎克雷伯菌KP2016多黏菌素耐药由mcr-1、mcr-8和crrB突变介导,替加环素耐药由tmexCD1-toprJ1基因簇介导。应加强合理监测,防止其在医疗机构中进一步传播。

成熟期水稻种子脱水速率全基因组关联分析

[目的]籽粒脱水速率直接影响种子的安全收获和快速干燥。选育成熟期籽粒含水量低、脱水速率快的品种,可保障种子质量,降低生产成本。[方法]采用来自82个不同国家和地区的165份水稻核心种质作为试验材料,将成熟期种子脱水速率表型与基因型相结合进行GWAS分析,挖掘调控种子脱水的关键基因,为培育和创制脱水快速品种奠定基础。[结果]1)对165份水稻核心种质群体脱水速率性状进行描述性统计分析,结果显示2年脱水速率性状均呈连续性偏正态分布,2年快速脱水的脱水速率性状在年际间具有显著相关性。不同亚群之间快速脱水与慢速脱水速率正相关。2)GWAS关联分析共获得与脱水速率显著关联的SNP 170个,QTL 36个。通过LD分析定位到6个与脱水速率密切相关的QTL,分别为qGDR2.3、qGDR4.1、qGDR4.2、qGDR6.1、qGDR6.4、qGDR10.1。[结论]这些QTL区间内主要候选基因OsPIP1;1、Os TIFY9、Osb ZIP48、Os ATG8b、OsDREB1C、Os SCP46与水分转运活性、信号转导、转录调控、抗氧化防御密切相关,推测它们与成熟期水稻种子脱水速率相关,可作为候选基因。

大豆籽粒Ve含量的全基因组关联分析

维生素E(Ve)是大豆油中一种天然抗氧化剂,是评价大豆油营养价值的重要指标。本研究利用含有264份的大豆自然群体在2021年和2022年测定了籽粒中α-、γ-和δ-生育酚含量,并进行全基因组关联分析(Genome-wide association study,GWAS)。本研究共检测到199个与大豆Ve含量显著关联的SNP位点,其中9个可在2个环境或者2个性状被重复检测到,分别位于3号、7号、11号、12号、13号、15号、17号和18号染色体上。其中位于7号染色体上的显著关联信号是控制α-生育酚含量的主效位点,可在2年环境中被检测到,表型变异解释率为9.83%。对该位点候选基因进行筛选,获得一个编码myb转录因子的基因Glyma.07G054000,可能是这个位点的效应基因。另外,在12号染色体上得到2个编码γ-生育酚甲基转移酶的基因Glyma.12G014200和Glyma.12G014300,有可能是影响Ve含量的重要基因。本研究结果有助于解析大豆籽粒Ve含量的遗传基础及其调控机制,为大豆品质遗传改良奠定了基础。

基于一步法全基因组关联分析解析蛋黄比率遗传结构

【目的】利用全基因组关联分析(genome-wide association study,GWAS)解析60周龄蛋鸡蛋黄比率遗传结构,旨在挖掘地方鸡种优异基因资源,进而为培育适应市场需求的特色蛋鸡品种奠定基础。【方法】数据采集自江苏省家禽科学研究所蛋鸡资源群体。依据F_(2)设计构建试验群体,分别以东乡绿壳蛋鸡、白来航鸡为亲本,经正反交获得F_(1)、F_(2)代。表型数据组包括开产日龄、胫长、体重、蛋重、哈氏单位以及蛋黄比率。收集F_(2)代产蛋鸡血液样本,以600K基因芯片对1534只试验鸡进行基因型分型。对SNPs数据进行质量控制、基因型数据填充。表型数据经去除离群值等清洗步骤后,以单因素方差分析检验影响性状表型值的因素,确定列入固定效应的因子。以系谱数据构建遗传关系矩阵,以多性状动物模型分析资源群体方差组分及遗传参数;作为对比,以系谱数据加基因型分型数据构建杂交遗传关系矩阵,用单性状动物模型计算遗传力,用双性状动物模型计算遗传相关系数。利用BLUPF90软件获得SNP效应值,以postGSf90获得SNP对应的P值以及权重,从而获得蛋黄比率关联的基因组显著水平SNP。针对蛋黄比率性状,计算染色体遗传力和单倍型片段长度。【结果】东乡绿壳蛋鸡蛋黄比率高于白来航鸡蛋黄比率,品种-批次效应对表型值影响显著,将其列入统计模型固定效应。应用杂交遗传关系矩阵分析,蛋黄比率、蛋重遗传力分别为0.348±0.041、0.500±0.038,蛋黄比率与蛋重、哈氏单位遗传相关系数分别为-0.463±0.075、-0.165±0.111;应用系谱遗传关系矩阵分析,蛋黄比率、蛋重遗传力分别为0.387±0.052、0.465±0.052,蛋黄比率与蛋重、哈氏单位遗传相关系数分别为-0.405±0.091、-0.166±0.121。蛋黄比率关联分析的膨胀系数为1.007,表明蛋鸡资源群体未检测到层化结构。鸡22号染色体存在与蛋黄比率显著关联的遗传座位,在显著性SNP之下有8个支持位点超过基因组潜在水平线。显著性以及潜在水平SNP可解释2.18%的表型方差,单倍型片段长度达到71 kb。此外,鸡2号染色体也检测到基因组潜在水平SNP。染色体遗传力分析表明,蛋黄比率受微效多基因控制,染色体遗传力与该染色体容纳的基因数呈正比例。【结论】应用加权一步法GWAS解析了蛋黄比率遗传结构,筛选得到ADAM9基因与60周龄蛋黄比率密切相关;相比于系谱遗传关系矩阵,杂交遗传关系矩阵减少孟德尔抽样误差,获得的遗传参数更准确。

大豆抗烟粉虱性状全基因组关联分析

近年来,烟粉虱在我国大豆主产区大面积发生,严重影响大豆产量。为筛选大豆抗烟粉虱候选基因,利用Illumina Soy6K SNP芯片对192份大豆自然群体进行基因分型,基于混合线性模型对大豆抗烟粉虱性状进行全基因组关联分析。共检测到3个与大豆烟粉虱抗性显著关联的SNP位点:Gm07_36488859、Gm11_37409028和Gm18_2054425。根据等位基因效应分析、基因注释和表达情况,最终获得大豆抗烟粉虱候选基因Glyma07g31470和Gly⁃ma11g35670。

花生含油量全基因组选择及近红外光谱筛选的育种技术探究

花生含油量对单位面积产油量至关重要。该性状受多个微效基因控制,但可用的紧密连锁标记十分有限,传统的分子标记辅助选择育种准确性不高。全基因组选择作为一种新的育种方法,可实现对数量性状的早期预测;近红外光谱分析可对作物品质性状(如含油量等)进行无损检测。通过两者优势互补,建立花生含油量全基因组选择和近红外光谱筛选联合的育种技术,探讨影响花生含油量全基因组选择预测准确性的因素,为花生分子育种奠定理论基础。本研究以216个重组自交系为材料构建训练群体;分别以139、464和505株F_(2)、F_(3)和F_(4)为材料构建育种群体;利用自主开发的“PeanutGBTS40K”液相芯片进行基因分型,开展含油量全基因组选择育种模型分析;通过联合全基因组选择和近红外光谱筛选技术,开展花生含油量性状的育种应用,并评价其育种效果。结果显示,对训练群体进行基因分型后,总共获得30,355个高质量SNPs,并用于11个全基因组预测的模型选择分析。含油量预测准确性最高的模型为rrBLUP,其次是randomforest和svmrbf。以重组自交系为预测群体,F_(2)、F_(3)和F_(4)各世代含油量的预测准确性分别为0.116、0.128和0.119;以重组自交系叠加上一轮的育种群体为预测群体,各世代含油量的预测准确性分别为0.116、0.131和0.160。全基因组选择联合近红外筛选要比单独的全基因组选择对各世代的含油量选择效果提高1.8%、2.7%和3.4%;与单独的近红外筛选相比,差异不显著(0.10%、0.06%和0.07%);而近红外筛选与全基因组选择相比,含油量可显著提高1.7%、2.6%和3.3%。通过联合全基因组选择和近红外光谱筛选育种,F_(3)和F_(4)分别比F_(2)的含油量提高1.2%和1.0%。F_(4)总共获得16个入选改良株系,有10个株系含油量≥55.0%,其中2个株系(SF_(4)_201和SF_(4)_379)的理论产量分别比对照增产7.0%和11.1%。本研究通过建立花生含油量性状的全基因组选择-近红外光谱筛选联合育种技术,可有效实现花生含油量性状的遗传改良。

高粱CIPK家族基因的全基因组鉴定及非生物胁迫下的表达特征

钙调磷酸酶B样蛋白互作蛋白激酶(CIPK)是一种重要的Ca^(2+)信号传感器,在植物应答逆境非生物胁迫过程中发挥着重要作用。为了探究高粱中CIPK家族基因的功能,本研究从高粱基因组中鉴定了31个SbCIPK基因,这些基因不均匀地分布在高粱的9条染色体上,编码蛋白质的氨基酸数量为403~519个,等电点为6.07~9.38,相对分子质量为46 357.31~58 316.97。基因结构分析结果表明,SbCIPK家族基因分为内含子缺失型和内含子富集型2类。进化树分析结果表明,SbCIPK家族蛋白质成员分为8个亚族。基于转录组数据的表达模式分析结果表明,SbCIPK基因广泛参与对盐胁迫、干旱胁迫等非生物胁迫的响应。本研究结果可以为高粱CIPK家族基因的功能研究奠定基础。

全基因组重测序解析重庆武隆中华蜜蜂遗传变异及种群分化特征

【目的】解析重庆武隆中华蜜蜂的遗传变异及种群分化特征,为该地区传粉昆虫的遗传多样性研究提供理论参考。【方法】对重庆武隆中华蜜蜂进行全基因组重测序,并结合我国其他地区57群中华蜜蜂的重测序原始数据进行种群遗传分化分析。【结果】测序共获得武隆中华蜜蜂有效数据(clean data)33.53 Gb,发现高质量单核苷酸多态性(SNP)位点共3886321个,小片段的插入与缺失(small indel)位点1392028个。在武隆中华蜜蜂样品中筛选出具有相同变异位点的非同义SNP突变基因589个,编码区发生small indel的基因214个。这些基因主要富集于赖氨酸降解、轴突导向、细胞外基质受体互作和丝裂原活化蛋白激酶信号通路。最终获得16个具有相同SNP和small indel的突变基因,其中包括嗅觉受体基因2个、细胞色素P450基因1个。武隆中华蜜蜂与我国其他地理型中华蜜蜂可能存在显著的遗传分化,与华中中蜂和北方中蜂(陕西安康)种群的亲缘关系较近。【结论】推测氨基酸代谢、神经系统发育、嗅觉进化、抗逆性能改善以及对温度偏好性的转变与武隆中华蜜蜂适应重庆本地自然环境有关。随意引种可能导致优质中华蜜蜂品种混杂,生产性能和抗逆性能下降。因此,需要保护本土优质中华蜜蜂种质资源。

娄门鸭胸肌性状相关指标的全基因组关联分析

为寻找与鸭胸肌性状相关的候选基因,鉴定鸭胸肌性状相关分子标记,试验以199只娄门鸭为试验材料,采集血液提取DNA,用基因组重测序方法进行基因分型,利用PLINK 1.90对分型结果进行质控后,采用rMVP软件中的FarmCPU模型对SNPs(Single nucleo⁃tide polymorphisims)与胸宽、胸深、胸肌重和胸肌率4个性状指标做GWAS分析,定位出显著位点。结果显示:与胸宽相关的SNPs位点有4个,位于1、11、18号染色体上;与胸肌重相关的SNPs位点有2个,位于2、3号染色体上;与胸肌率相关的SNPs位点有2个,位于2号染色体上。通过基因功能注释,这些位点对应的胸肌性状候选基因有8个,分别是LOC119715526、RLBP1、ABHD2、KYAT1、SPOUT1、LOC101800569、EVA1A和NOL4。研究表明,在不同染色体上共鉴定到8个与胸肌性状显著相关的SNPs位点,这些位点共涉及8个候选基因。

基于全基因组重测序评估郏县红牛保种群遗传多样性与群体结构

【目的】了解郏县红牛保种群的群体结构和遗传多样性,从而更好地保护和利用郏县红牛遗传资源。【方法】基于全基因组重测序技术,对30头郏县红牛的全基因组单核苷酸多态性(single nucleotide polymorphism,SNP)进行检测,并通过分析群体遗传多样性、群体结构、连续纯合片段(runs of homozygosity,ROH)分布特征和亲缘关系,对郏县红牛保种群的保种效果进行综合评估。【结果】在30头郏县红牛个体中共检测到41215797个高质量SNPs位点;郏县红牛群体遗传多样性丰富,最小等位基因频率为0.13±0.13,多态性标记比例为0.58±0.33,期望杂合度为0.192±0.152,观测杂合度为0.191±0.158,核苷酸多样性为0.003±0.001。郏县红牛群体的有效群体含量呈逐代下降趋势,在1000代前有效群体含量为2716头,在20代前为143头;郏县红牛个体间的遗传距离介于0.80~0.89之间,平均遗传距离为0.83±0.02。聚类分析显示,30头郏县红牛共分为7个家系,15头公牛可分为5个家系;30头郏县红牛个体中共检测到5947个ROH,总长度为2.8 Gb,长度为0~0.5 Mb的ROH占比最多(73.08%),2~4 Mb的ROH占比最少(0.40%),基于ROH的平均近交系数为0.19±0.07,15头公牛的平均近交系数为0.15±0.05。【结论】郏县红牛保种群遗传多样性丰富,群体结构没有出现明显分层,整体上保持了较高水平的近交,出现了一定的近交风险,需加强选种选配工作,以确保郏县红牛遗传资源的可持续发展。

弯曲菌宽谱烈性噬菌体生物学特性、全基因组测序分析及初步应用

为探索控制弯曲菌污染的新途径,以弯曲菌NCTC12662为宿主菌,从养殖场鸡粪便中分离纯化得到弯曲菌噬菌体,并对其进行电镜形态特征观察、生物学特性测定与基因组学分析,同时采用棋盘法,测试噬菌体与抗生素的联合应用效果。结果显示,从山东省胶东地区养鸡场粪便中分离到1株噬菌体P-61。透射电镜下显示该噬菌体颗粒具有1个直径约为102.6 nm的二十面体头部和1个长度约为112.9 nm的长尾部,属于有尾噬菌体目,长尾噬菌体科;裂解谱为48.48%,裂解效价为5.5×10^(11)pfu/mL,感染复数为0.01;一步生长曲线显示潜伏期为20 min,裂解期约为130 min,裂解量为240 pfu/cell;最适宜生长温度为42℃,pH耐受性较强,在pH 5~9范围内效价较稳定。噬菌体基因组全长为130425 bp,G+C含量为26.13%,共有176个开放阅读框、3个tRNA,软件预测基因组中不含毒力基因和抗生素耐药基因。利用棋盘法测得,当噬菌体为10^(6) pfu/mL时,环丙沙星最低浓度降至1/8 MIC,联合用药组具有明显的联合杀菌效果。综上可知,得到的噬菌体P-61具有效价高、裂解谱宽及对环境适应能力强等生物学特性,且基因组信息清晰,能与抗生素联合应用发挥协同作用,具有消毒、净化等临床应用可能性。本研究为噬菌体制剂的开发利用提供了生物资源。