浙江缢蛏群体糖原含量变化与gys基因多态性的关联分析

为了探讨缢蛏糖原含量与糖原合酶基因(sc-gys)的特性及其用于分子标记辅助育种的可行性,实验检测了浙江缢蛏群体糖原含量的周年变化,克隆获得糖原合酶基因的cDNA全长序列,分析其在不同组织、不同月份的表达差异,并进一步筛选出与糖原含量相关的单核苷酸多态性(single nucleotide polymorphism, SNP)位点。结果显示,养殖缢蛏在2—6月糖原含量较高,4月达到最大值。sc-gys基因的cDNA全长为2 402 bp,开放阅读框为2 136 bp,编码711个氨基酸。荧光定量PCR结果显示,sc-gys基因主要在外套膜和足中表达,其不同月份的表达水平与糖原含量周年变化规律一致,4月表达量显著高于其他月份,表明该基因与糖原合成关系密切。在"甬乐1号"品种和台州野生群体sc-gys基因的外显子区共筛选到2个与糖原含量相关的SNP位点,且两位点之间没有连锁不平衡,优势基因型组合个体的糖原含量比群体均值提高11.8%,这些位点和优势基因型组合可为缢蛏糖原含量的遗传改良提供候选标记。     

CYP2J3样基因在不同生长速率凡纳滨对虾中的差异表达

为了获得对虾生长差异的分子机制,通过比较不同生长速率凡纳滨对虾的转录组(速长组和慢长组),获得了一条编码P450的差异表达基因,经过cDNA全长验证后命名为LvCYP2J3-like,其开放阅读框为1 488 bp,编码495个aa,含有一个完整的P450功能域,功能分析表明,该基因可能参与了蜕皮激素的灭活。方差分析表明,LvCYP2J3-like在不同组织、发育阶段和蜕皮周期的表达量都存在显著性差异,并且速长群体的表达量显著性低于慢长群体。表明该基因参与了发育和蜕皮过程,其表达量与环境胁迫关系密切。在LvCYP2J3-like和其预测转录因子基因中分别发现了4和11个速长组和慢长组等位基因频率显著差异的单核苷酸多态性(SNP)标记,说明这些基因遗传差异与其活性和表达量相关,从而影响生长速率。研究结果为凡纳滨对虾遗传和育种研究提供了参考资料。     

大口黑鲈组织蛋白酶B基因SNP筛选及其与生长性状关联分析

在动物体内组织蛋白酶B(cathepsin B,CTSB)是以酶原形式存在于溶酶体中,能水解多种蛋白,是调控动物生长发育的重要候选功能基因。为开发与大口黑鲈(Micropterus salmoides)生长性状相关的单核苷酸多态性(single nucleotide polymorphism,SNP)标记,本研究对大口黑鲈EST-SNP库中CTSB基因上的SNP位点进行筛选。经PCR扩增和测序分析,从该基因中筛选到1个SNP位点(命名为C+598G),该位点属于错义突变位点,编码的甘氨酸突变成为精氨酸。采用Sna Pshot方法对C+598G位点在165尾大口黑鲈养殖样本中的基因型进行检测,并分析其在群体中的遗传多样性。结果显示,C+598G位点在养殖群体中存在CC、GC和GG3种基因型,其频率分别是40.6%、47.3%和12.1%,C和G等位基因频率分别是64.2%和35.8%;该位点在群体中的观测杂合度(Ho)和多态信息含量(PIC)分别为0.473和0.354,符合哈代-温伯格平衡定律。进一步分析不同基因型与大口黑鲈生长性状相关性,结果表明,3种基因型个体的体质量、全长、体高、尾柄长和尾柄高性状平均值从小到大依次为GG型、GC型和CC型,其中GG型个体在体质量上显著小于CC型个体(P<0.05)。本研究获得的C+598G标记与大口黑鲈生长性状显著相关,可作为大口黑鲈分子标记辅助育种研究中的生长标记。     

牙鲆(Paralichthys olivaceus)丝氨酸蛋白酶Ⅰ-1基因多态性

采用重测序方法,利用牙鲆(Paralichthys olivaceus)野生个体,扩增了牙鲆丝氨酸蛋白酶(Serine protease) I-1基因的编码区(CDS)和2605 bp的启动子序列。在牙鲆丝氨酸蛋白酶I-1基因(PoSP I-1)编码区中筛选出9个单核苷酸多态性位点(SNPs),同时在其启动子区筛选出28个SNP位点。对感染鳗弧菌的抗病牙鲆和易感牙鲆进行了部分SNP位点的分型,通过分析CDS中SNP位点的分布情况,发现SNP365A/G位点在抗病个体中的AG基因型的基因频率是60%,高于易感个体的40% (P=0.01)。qRT-PCR结果显示,易感个体在细菌感染后,PoSP I-1相对表达量相比对照组呈下降趋势。在抗病个体中,PoSP I-1相对表达量相比对照组呈上升趋势,且高于易感个体。结果表明,丝氨酸蛋白酶基因在牙鲆抗鳗弧菌感染中发挥重要作用,该基因的SNP365A/G位点是与牙鲆鳗弧菌感染相关的候选位点,该位点可作为牙鲆抗鳗弧菌选择育种的潜在标记。     

三角帆蚌HcTyr基因内壳色性状相关SNP筛选及图谱定位

为研究三角帆蚌HcTyr基因与内壳色性状的相关性,本实验根据已分离的三角帆蚌HcTyr基因进行引物设计和片段扩增,并用144只三角帆蚌对其多态性进行筛选和验证,卡方检验分析了SNP位点和三角帆蚌内壳色相关性,并对HcTyr基因进行了图谱定位。结果显示,在HcTyr基因上筛选出8个SNP位点,其中有7个SNP位点与三角帆蚌内壳色L、a及d E呈显著相关性,用这7个SNP位点做单倍型构建和分析,发现Ⅱ、Ⅲ及Ⅳ这3种单倍型在白色群体中出现的频率极显著高于在紫色群体中出现的频率,而Ⅴ和Ⅶ这2两种单倍型在紫色群体中出现的频率极显著高于白色群体。进一步在商业养殖群体中验证发现,Ⅳ和Ⅶ这2种单倍型可分别作为白色和紫色群体的优势单倍型。研究表明,三角帆蚌HcTyr基因可作为内壳色相关的候选基因,其部分SNP位点可用于三角帆蚌分子辅助育种。另外本实验还将HcTyr基因定位在三角帆蚌LG16连锁群上,这为进一步解析该基因调控珍珠颜色形成的分子机制奠定了基础。     

大口黑鲈肌球蛋白重链基因SNPs的筛选及与生长性状的关联

肌球蛋白是肌肉细胞的重要组成成分,影响肌纤维的组成和肌肉生长。为了研究肌球蛋白重链(MYH)基因多态性与大口黑鲈生长性状的相关性,本研究采用PCR技术扩增得到编码区序列全长为9759 bp的大口黑鲈MYH基因,该基因包含37个外显子和36个内含子,编码1940个氨基酸。采用直接测序法在MYH基因上筛选到8个单核苷酸多态性标记(SNP)位点(A-305G、G-558C、A-2784C、A-2816G、T-4765A、C-6206T、C-6811T和G-6935T),有4个位于外显子上,其中2个属于同义突变。用SNa Pshot方法对从同批繁殖、同塘养殖的大口黑鲈"优鲈1号"群体中随机选取的430尾个体中各位点的基因型进行检测。结果显示,内含子上的A-2784C和A-2816G位点完全连锁,所有位点在大口黑鲈"优鲈1号"群体中的平均有效等位基因数、平均观测杂合度和平均期望杂合度分别为1.635、0.406和0.373,仅C-6206T、C-6811T和T-4765A位点的基因型频率分布符合Hardy-Weinberg定律。采用一般线性模型分析各位点与大口黑鲈生长性状之间的相关性,研究发现,C-6811T位点CC基因型个体的体质量和全长显著大于TT基因型,CC基因型个体的体高和尾柄长显著大于CT和TT基因型,其余位点不同基因型个体间的生长性状均不存在显著差异。C-6811T位点与生长性状显著相关,可作为大口黑鲈分子标记辅助育种的候选标记。     

辽东湾斑海豹线粒体ND4、tRNAArg、ND4L和ND3基因序列分析

用PCR产物直接测序法,测定渤海辽东湾斑海豹15个个体线粒体的一段从ND3到ND4基因长度为625 bp的DNA序列。其中的1~118 bp是ND3基因的后部分,120~188 bp是tRNAArg基因完整序列,189~485 bp是ND4L基因完整序列,479~625 bp是ND4基因的前部分。在ND4L和ND4之间,有7个碱基是重叠编码的。15个序列间,只有ND4L基因中有一个位点发生C/T碱基转换。与GenBank中的斑海豹序列比较,tRNAArg基因序列完全一致,而ND4L基因序列则有两处发生转换。利用从GenBank中下载的22个海豹科动物的ND4L基因序列的系统发生分析表明,海豹科动物分为南半球和北半球两大类群,斑海豹与港海豹为北半球种类,亲缘关系最近。     

施氏鲟性别相关的微卫星标记筛选

用32个微卫星标记分析了30尾（15尾♀;15尾♂）性成熟（12龄）施氏鲟（Acipenser schrenckii）雌雄个体的遗传差异。结果表明：雌、雄群体的遗传多样性差异没有达到显著水平。在6个微卫星标记中有9个等位基因在雌、雄群体中分布比例差异达到50%以上,其中雌性群体HLJSX 194-239bp、HLJSX201-207bp、HLJSX209-210bp、HLJSX210-178bp和HLJSX215-221bp 5个等位基因频率显著高于雄性群体;雄性群体HLJSX194-184bp、HLJSX201-214bp、HLJSX210-233bp和HLJSX226-194bp 4个等位基因频率显著高于雌性群体,初步表明这些基因位点可能与性别决定位点存在一定连锁关系。本研究结果为施氏鲟早期性别鉴定积累了基础数据。     

三角帆蚌HcTyp-1基因珍珠层颜色性状相关SNP筛选及图谱定位

本实验室根据前期研究发现,HcTyp-1基因可能参与了三角帆蚌贝壳珍珠层颜色的形成。本实验根据HcTyp-1基因cDNA设计引物,比较筛选获得了8个SNP位点并分析了这些位点与内壳色的相关性,然后对HcTyp-1基因进行了图谱定位。通过研究这些多态性位点与贝壳珍珠层颜色性状相关性发现,C+3057T位点基因型仅在a参数上存在显著差异,G+2985T和T+3006C 2个位点基因型分别在L、b及a、dE上存在显著差异,A+2834C、C+2919T和G+2986T这3个SNPs位点基因型在L、a及dE上均存在显著差异。连锁不平衡结果显示,HcTyp-1基因上除C+2912T、C+2983A这2个位点,其他6个位点具有显著差异,位点之间都存在强连锁不平衡。单倍型分析结果显示,单倍型T2和T4在白色群体中出现的频率极显著高于在紫色群体中出现的频率,T6和T8两种单倍型在紫色群体中出现的频率极显著高于白色群体。因此HcTyp-1基因的部分SNPs可作为三角帆蚌不同内壳色贝壳辅助育种的候选分子标记位点。另外,本研究还将HcTyp-1基因定位在三角帆蚌LG16连锁群上,这为进一步解析该基因调控珍珠颜色形成的分子机制奠定了基础。     

大口黑鲈POU1F1启动子区域SNPs对生长的影响

采用基因组步移技术克隆了大口黑鲈垂体特异性转录因子POU1F1启动子序列,该序列全长1 629 bp,预测存在4个八聚体转录因子1(Oct-1)结合位点和TATA框等基本转录元件,1个Homeobox转录因子结合位点和2个cAMP反应元件结合蛋白(CREB)结合位点。采用直接测序法研究发现,大口黑鲈POU1F1启动子序列的-18和-183两个位点存在SNPs,该两个突变位点只检测到A、B两种单倍型。对养殖群体进行不同单倍型的基因频率和生长关联分析,结果表明,单倍型A的等位基因频率为0.246,单倍型B的等位基因频率为0.754;群体关联分析表明,基因型为AA型和AB型的个体在体质量、体长、全长、体高和体宽等方面明显高于BB型个体,但基因型为AB型和AA型的个体之间在体质量方面没有明显差异。该研究结果为下一步进行分子标记辅助选择提供了基础。     

仿刺参线粒体全基因组序列结构及比较研究

利用已构建的仿刺参cDNA文库得到的线粒体DNA(mtDNA)相关基因序列设计扩增引物,测定了大连仿刺参线粒体基因组全序列,并对其进行了基因构成和进化分析。仿刺参线粒体基因组序列长16 109bp,其基因构成与其他后口动物基本一致,包括37个基因(2个rRNA基因、22个tRNA基因和13个蛋白质编码基因)和3个主要的非编码区。在其37个基因中,ND6、tRNASer(AGN)、tRNAGln、tRNAAla、tRNAVal、TrnaAsp位于L链上,其余均位于H链上。在13个蛋白质编码基因中,除ND1的起始密码子为GTG外,其余均以ATG作为起始密码子;除Cytb以"T"作为终止密码子外,其他蛋白质基因均具有完全的终止密码子,且在已知的棘皮动物线粒体蛋白质基因中,部分基因的起始和终止密码子表现出一定的纲内特异性。比较分析了大连、青岛、威海仿刺参线粒体基因组,三者的基因组成和排列相同,碱基组成相近,蛋白质编码基因的起始和终止密码子完全一致,但存在核苷酸和氨基酸序列的差异。三者的控制区序列存在多个插入/缺失和SNP位点。根据COI、Cytb和ND4计算了三者之间的遗传距离为0.006～0.018,遗传距离分析和系统进化关系分析都显示青岛仿刺参和威海仿刺参的关系较大连仿刺参更近。     

基于线粒体COII和ND4基因序列的塔里木裂腹鱼群体遗传多样性分析

为探究塔里木裂腹鱼群体的遗传多样性和遗传分化现状,选用车尔臣河、克孜勒河和阿克苏河3个不同地理种群的塔里木裂腹鱼共计126尾样本,进行线粒体DNA COII和ND4基因序列的对比分析,探讨了2种标记对塔里木裂腹鱼遗传多样性分析结果的差异和关系。结果显示,塔里木裂腹鱼线粒体DNA COII和ND4基因序列的A+T含量均高于G+C含量,碱基组成具有偏倚性。基于线粒体DNA COII和ND4基因序列分析,126个样本中分别确定了6个和23个单倍型,其中线粒体DNA COII基因序列中3个群体存在共享单倍型现象,线粒体DNA ND4基因序列中未发现3个群体间存在共享单倍型。2种标记下群体的单倍型多样性和核苷酸多样性分别为0.738±0.019/0.904±0.014和0.02129±0.00298/0.04876±0.00149,呈现出高单倍型多态性和高核苷酸多态性的特征。2种标记的分子方差分析（AMOVA）均表明,在所有群体中遗传变异主要来源于群体间,且车尔臣河群体的塔里木裂腹鱼与其他2个群体的遗传分化均达到显著水平（P<0.05）。贝叶斯法（BI）构建的BI树与单倍型网络图结构一致,塔里木裂腹鱼形成了2个分支。COII和ND4基因序列的岐点分布图均呈现双峰型,表明塔里木裂腹鱼现在的分布是先前分化种群发生二次交流的结果。种群结构分析结果亦表明塔里木裂腹鱼已分化出2个明显的地理种群,故建议将车尔臣河群体作为塔里木裂腹鱼的亚种。     

多巴色素异构酶基因对瓯江彩鲤黑斑体色的影响

多巴色素异构酶（dopachrome tautomerase，Dct）是酪氨酸酶家族的重要成员之一，在黑色素合成过程中具有重要作用。为了解多巴色素异构酶基因（）的“粉花”体色类型为研究对象，在对该基因进行克隆和表达分析的基础上，应用CRISPR/Cas9基因编辑技术研究了该基因对瓯江彩鲤黑斑体色的影响。结果发现，瓯江彩鲤存在2个Dct1与在“粉花”体色瓯江彩鲤的黑斑皮肤中的表达量显著高于白色皮肤，而Dct1和基因突变对瓯江彩鲤黑斑体色存在重要影响，与在鲤的黑斑体色形成中发挥着更大的作用。     

基于黄河鲤新品系IGF2b基因SNPs获取体系的建立及在该品系选育中的应用

为了建立鲤IGF2b (insulin-like growth factor 2)基因全基因组信息SNPs (single nucleotide polymorphism)的获取体系,并验证该获取体系在黄河鲤新品系选育中运用的适用性,本研究首先对10个不同鲤品种的33个鲤个体重测序数据进行分析,整体上了解IGF2b基因的单核苷酸位点及其在基因组上的分布情况。然后在基因组DNA水平上分段获取IGF2b基因的SNPs信息,并以不同的鲤品种验证引物的适用性,最后在黄河鲤新品系中检测其应用。研究结果:获得8对聚合酶链式反应(polymerase chain reaction, PCR)引物进行分段克隆的结果,并通过直接测序和限制性片段长度多态性来获取相应的SNPs。经过直接测序,本文对黄河鲤、建鲤、黄河鲤和建鲤的正反交后代进行PCR检测发现条带单一,目地片段正确;结合黄河鲤新品系的体重数据,本文又以IGF2b~(3#)引物检测到一个与体重明显相关的SNPs,同时检测到一个与该基因表达量相关的另一个SNPs。本研究的方法能够在该基因全基因组范围内开展鲤分子育种中SNPs的检测,并验证发现黄河鲤新品系IGF2b~(3#)引物227号位点处,出现的突变位点使体重降低,且有一处SNPs与其表达量有关。这将为其他物种、其他功能基因的多态性分子标记研究提供思路。     

黄河鲤酪氨酸酶基因的克隆、表达模式及定位分析

为探究Tyrosinase(Tyr)基因在黄河鲤体色形成中的作用,利用生物显微镜对黄河鲤黑色鳞片、银色鳞片、鳍条色素细胞的组成进行观察,并利用RACE克隆黄河鲤Tyr基因全长cDNA序列,荧光定量PCR分析其在不同组织中的表达情况,Western blot印迹和免疫组织化学方法检测其蛋白的表达定位。结果显示,黄河鲤具有黑色素细胞、黄色素细胞和虹彩细胞3种色素细胞,3种色素细胞的数量、种类和分布在不同组织中存在较大差异。黄河鲤Tyr基因cDNA全长1717 bp(GenBank ID:No.KY305667),其中ORF长1605 bp,编码535个氨基酸,5′-UTR区41 bp,3′-UTR区71 bp。蛋白同源分析发现,黄河鲤Tyr与鲤科鱼类相似性最高,与哺乳类及鸟类等脊椎动物相似性最低。系统进化分析显示,黄河鲤Tyr与鲤、瓯江彩鲤、锦鲤和斑马鱼等鲤科鱼类的亲缘关系最近。实时定量PCR分析表明,Tyr基因在黄河鲤不同组织中均有表达,肌肉中的表达量最高,其次是黑色皮肤,另外在黄色皮肤、臀鳍和尾鳍等黑色素细胞存在的组织中也有较高表达。Western blot印迹结果显示,Tyr基因在黑色皮肤中表达最高,其次是臀鳍和尾鳍,银色皮肤中没有表达。免疫组织化学结果显示,Tyr基因在黑色皮肤的黑色素细胞和黏液细胞中有表达。研究表明,Tyr基因表达水平与黑色素细胞数量和分布具有一定相关性,参与黄河鲤体色的形成。     

三疣梭子蟹线粒体基因组单核苷酸多态性

为获得三疣梭子蟹线粒体全基因组(mt DNA)SNP突变位点,本研究针对其mt DNA全序列设计了179对引物,采用高分辨率熔解曲线(HRM)检测技术对SNP进行了筛查。在179对引物中,有89对引物具有特异性扩增结果;进一步研究表明,89对引物中共有22对引物的扩增区域含有SNP突变位点,共包含24个SNPs。统计结果显示,三疣梭子蟹线粒体基因组中的SNP分布频率为0.15/100 bp,其中转换突变比例为79.1%、颠换突变比例为16.6%;C/T(G/A)突变比率为79.1%、G/T(C/A)突变比率为8.3%、A/T突变与G/C突变的比率均各占4.16%。本研究中所获得的SNP位点分布于三疣梭子蟹线粒体基因组的11个区域,且分布不均。其中COX1基因区域的SNP数目最多,其次在D-LOOP、ND1基因区域分别为3和4个SNP位点;而t RNA、12S r RNA等其他区域尚未发现突变位点。本研究所建立的三疣梭子蟹线粒体全基因组SNP的快速筛查方法及发现的SNPs,为进一步开展三疣梭子蟹种质遗传资源多样性、增殖放流标记等方面的研究提供了分析工具。     

艾氏蛇鳗线粒体基因组全序列结构分析和系统发育关系探讨

本研究采用高通量测序技术获得了艾氏蛇鳗(Ophichthus evermanni)线粒体基因组全序列, 并对其结构和特征进行了分析。结果表明, 艾氏蛇鳗线粒体基因组全长 17759 bp, 包含了 13 个蛋白编码基因(PCGs)、22 个转运 RNA 基因(tRNA)、2 个核糖体 RNA 基因(rRNA)、2 个控制区(D-loop)和 1 个轻链复制起始区(O_L)。线粒体 DNA 全序列的碱基组成分别为 A (31.27%)、G (16.19%)、C (26.22%)和 T (26.32%), 其中 A+T 含量(57.59%)大于 G+C 含量 (42.41%), 呈现出明显的 A+T 偏好性。与大多数硬骨鱼类不同, 艾氏蛇鳗线粒体基因组中发生了基因重排现象, ND6 基因和 tRNA-Glu 移到了 tRNA-Thr 和 tRNA-Pro 之间, 且 ND6 基因上游还存在另一个高度同源的 D-loop 区。 tRNA-Gln (Q)、tRNA-Ala (A)、tRNA-Asn (N)、tRNA-Cys (C)、tRNA-Tyr (Y)、tRNA-Ser^UCA (S1)、tRNA-Glu (E)、 tRNA-Pro (P)和 ND6 9 个基因位于 L 链, 其余基因均位于 H 链。除 tRNA-Ser (AGC)外, 其余 21 个 tRNA 均为典型的三叶草二级结构。分别采用邻接法和最大似然法, 基于 12 个蛋白编码基因(ND6 除外)构建了蛇鳗科鱼类系统发育关系树。结果显示艾氏蛇鳗与短尾蛇鳗(O. brevicaudatus)和食蟹豆齿蛇鳗(Pisodonophis cancrivorus)的亲缘关系较近, 蛇鳗属是蛇鳗科鱼类中分化较晚的一个类群。研究结果丰富了蛇鳗科鱼类线粒体基因组数据库, 也为该类群鱼类的系统分类研究提供了参考资料。     

长牡蛎MITF基因表达及与壳色的关联

为了解长牡蛎MITF基因的表达及其与壳色的关联,验证了长牡蛎中的4个MITF基因,对长牡蛎MITF氨基酸序列进行了序列分析和多序列比对,分析了长牡蛎幼体发育各时期的转录组,采用荧光定量PCR的方法研究长牡蛎各组织及黑壳、白壳牡蛎特定部位mRNA表达情况。4个MITF基因中有3个基因可能为假基因,有表达的长牡蛎MITF基因共编码448个氨基酸,为亲水性不稳定蛋白,含有N端结构域(MITFTFEBC3N superfamily)和高度保守的功能性结构域HLH结构域。转录组分析发现MITF在长牡蛎个体发育的各个时期均有表达,在稚贝期达到最高。组织表达结果显示MITF在外套膜中的表达水平显著高于其他组织。MITF在黑壳牡蛎闭壳肌中的表达量显著低于白壳牡蛎,而在黑壳牡蛎外套膜中的表达量高于白壳牡蛎,但不显著。在黑壳、白壳牡蛎外套膜边缘的表达量都显著高于内侧。研究表明,长牡蛎MITF基因可能在牡蛎壳形成早期就参与了黑色素的生成调控和个体的生长发育,可调控牡蛎酪氨酸酶Tyr2基因,参与外套膜和贝壳中黑色素的形成。本研究为进一步研究牡蛎壳色形成机制奠定了基础。     

草鱼mst2在免疫应答中的作用机制

为了初步阐明草鱼哺乳动物STE20样蛋白激酶2基因(mammalian sterile 20-like kinase 2, mst2)在机体免疫中的作用机制,实验采用RNA-Seq技术对干扰mst2后经脂多糖(lipopolysaccharide,LPS)应激的草鱼肾细胞系(Ctenopharyngodon idella kidney cell lines,CIK)进行了转录组测序与验证分析。测序原始数据经De novo拼接与组装后共获得22 374个独立功能基因(unigenes),其中已知功能基因为21 199个,预测的新基因为1 175个。干扰mst2后经LPS应激的unigenes表达差异分析表明,对照组与实验组之间共存在38个差异基因(differentially expressed genes,DEGs),其中上调基因16个,下调基因22个。利用实时荧光定量PCR技术(quantitative real-time PCR, qRT-PCR)对38个DEGs的RNASeq结果进行验证,结果显示,qRT-PCR和RNA-Seq分析一致,说明RNA-Seq分析结果可靠。采用RNA干扰技术干扰mst2后经LPS处理,CIK细胞转录组中DEGs参与免疫代谢的途径主要有MAPK信号通路、内吞作用途径、自噬途径和细胞因子受体相互作用途径。凋亡相关基因检测结果显示,干扰mst2并经LPS处理后,促凋亡基因(fas、bad1、bad2、caspase-3、caspase-8和caspase-9)转录水平上调,抗凋亡基因(bcl2)转录水平下调,证明干扰mst2后经LPS处理会诱发细胞发生凋亡。综上所述,mst2可通过调控凋亡相关过程参与机体免疫反应。本研究结果初步阐明了草鱼mst2参与机体免疫应答的分子机理,可为草鱼细菌性疾病防控提供一定的基础理论参考。