基于mRNA-miRNA相互作用网络研究阿尔茨海默病机理

本研究基于表达谱数据系统分析了阿尔茨海默病中的差异表达基因;进一步的GO和KEGG富集分析研究了差异表达基因参与的生物学功能和生物学过程;最后通过mRNA-miRNA相互作用网络,挖掘了阿尔茨海默病中的关键基因和调控机理。研究表明,990个基因在阿尔茨海默病中差异表达(p-value≤0.01,|log FC|≥2),其中332个基因(33.5%)上调,658个基因(66.5%)下调。功能富集(GO)表明,差异基因参与了Regulation of macrophage activation、Myelin sheath和structural constituent of cytoskeleton等功能。KEGG通路富集分析表明,差异基因参与了,Synaptic vesicle cycle、PI3K-Akt signaling pathway、Nicotine addiction、Dopaminergic synapse和Retrograde endocannabinoid signaling等重要的生物学通路。最后,mRNA-miRNA相互作用网络鉴定了在阿尔茨海默病中6个关键的基因,包括TP53INP1、MET、MBNL3、CEBPB、GNAS和SMARCA4,其中TP53INP1和MET在前人的研究中有报道与阿尔茨海默病密切相关,MBNL3、CEBPB、GNAS和SMARCA4是新发现的一些基因。这些基因可能参与了阿尔茨海默病的发生和发展,可以作为其调控位点和药物靶点。     

丙酸通路基因多态性与猪肉质及胴体性状的关联分析

为了挖掘新的猪肉品质及胴体性状的候选基因,揭示猪肉质及胴体性状的遗传机制,文章将丙酸代谢通路作为候选通路,将通路内基因与猪肉质及胴体性状进行关联分析。实验采用37头三元杂交商品猪作为研究对象,首次针对丙酸通路中7个基因的36个SNP位点利用SNaPshot方法进行基因分型,分别用最小二乘模型及MB-MDR模型与肉质及胴体性状进行关联分析。结果发现,基因PCCB、MUT、MCEE及ACSS2上的4个SNP位点分别与肌内脂肪含量、背膘厚等性状显著相关(P<0.05),ACSS2与猪脂肪含量显著相关;MCEE及MUT与猪的背膘厚显著相关;PCCB基因与脂重显著相关。通过MB-MDR方法检测到多个SNP位点具有互作效应,并与背膘厚、水分含量、脂肪含量显著相关(P<0.05)。另外,丙酸代谢通路中的基因间的互作效应对猪肉品质有显著影响。     

重症肌无力风险基因挖掘

目的:挖掘重症肌无力(Myasthenia gravis,MG)可能的风险基因。方法:通过人工挖掘在Pub Med数据库收集重症肌无力风险基因,通过Gene数据库获取重症肌无力风险基因编号,用以表示基因或者其相应的蛋白。应用基因功能分析软件DAVID(http://david.abcc.ncifcrf.gov/)对重症肌无力风险基因进行KEGG通路富集分析,挖掘重症肌无力风险通路,进而对任意两个通路进行关联分析。应用基因功能分析软件DAVID的Gene Ontology,对MG风险基因进行功能注释,以P0.01来判定注释是否有显著意义。结果:(1)本研究挖掘出97个重症肌无力的风险基因,KEGG基因富集分析共筛选出44条与重症肌无力显著相关的通路,主要包括多种自身免疫性疾病相关通路、信号转导相关通路、肿瘤相关通路、抗原的加工提呈通路等等。(2)以上44条风险通路两两通路间均具有相关性(P.0.01)。结论:本研究共挖掘出44条重症肌无力风险通路,8个重症肌无力风险基因,分别为:NF-kB、TNFR、MEK、AP-1、Raf、MEK1/2、MSK1、TAPBP。其中,MEK同时出现在多个风险通路中,考虑其风险性更高。     

IgA 肾病患者与膜性肾病患者外周血单核细胞mRNA分析

为了寻找诊断、鉴别IgA肾病(IgAN)和膜性肾病(MN)的血液特异性标记物,利用公共数据库中的IgAN和MN患者的外周血单核细胞(PBMCs)的转录组表达谱数据集识别特异性生物标记物,为诊断和鉴别提供简便、可靠的依据补充。从公共基因表达数据库(GEO)下载IgAN患者组(n=15)和MN患者组(n=8)芯片数据集,筛选前250个差异表达基因(DEGs)。通过分析筛选关键基因和途径,进行基因本体(GO)富集分析、京都基因与基因组百科全书(KEGG)通路分析和蛋白质与蛋白质相互作用关系(PPI)分析等进一步了解DEGs。通过分析共发现75个显著DEGs,其中73个上调基因,2个下调基因。GO富集分析的生物学过程(BP)主要包括蛋白质转运、内溶酶体到溶酶体转运、趋化因子介导的信号通路作用等。显著富集差异表达基因KEGG通路分析包括Endocytosis和Hepatitis B的相关信号通路。PPI筛选出EPS15、STAT4、CCL2、SUN2、SEC24C、SEC31A、GOLGB1、F2R,RAB12和PTK2B等关键基因。成功筛选出核心差异表达基因,为IgAN和MN的诊断和鉴别提供简便、可靠的依据补充,甚至提供治疗的新靶点。     

基于蛋白质互作知识的生物学通路扩充新方法

生物学通路被广泛应用于基因功能学研究, 但现有的生物学通路知识并不完善, 仍需进一步扩充。生物信息学预测为通路扩充提供了一种有效且经济的途径。文章提出了一种融合蛋白质-蛋白质互作知识以及Gene Ontology(GO)数据库信息进行基因通路预测的新方法。首先选取目标基因在蛋白质-蛋白质互作层面上的邻居所在的Kyoto Encyclopedia of Genes and Genomes(KEGG)通路为候选通路, 然后通过检验候选通路中的基因是否在与目标基因关联的GO节点富集来判断目标基因的通路归属。分别利用Human Protein Reference Database (HPRD)和Biological General Repository for Interaction Datasets(BioGRID)数据库中的蛋白质-蛋白质互作信息进行预测。结果表明, 在两套数据中, 随着互作邻居个数的增加, 预测的平均准确率(在所有目标基因注释的通路中被成功预测的比例)及相对准确率(在至少有一个注释通路被成功预测的基因集中, 所有注释通路均被预测正确的基因所占的比例)均呈现上升趋势。当互作邻居个数达到22时, 预测的平均准确率分别达到96.2%(HPRD)和96.3%(BioGRID), 而相对准确率分别为93.3%(HPRD)和84.1%(BioGRID)。进一步利用新版数据库对旧版数据库中被更新的89个基因进行验证, 至少有一个更新通路被预测正确的基因有50个, 其中43个基因的更新通路被完全正确预测, 相对准确率为86.0%。这些结果显示该方法是一种可靠且有效的通路扩充方法。     

Noggin基因沉默对BMP和Wnt信号通路表达的影响

目的检测分析Noggin基因沉默对毛囊发育中BMP和Wnt信号通路的影响。方法采用实时荧光定量PCR和western blot技术对Noggin基因沉默的MC3T3-E1稳转细胞系中BMP-2、BMP-4、BMPR-IA、BMP-6、BMP-7、LEF-1、β-catenin的表达情况进行检测分析。结果实时荧光定量PCR结果显示,BMP信号通路中的五个基因的表达都受到Noggin基因沉默的显著的影响,其中BMP-2(P0.001)、BMP-4(P0.01)、BMP-6(P0.001)、BMP-7(P0.001)表达量均升高;BMPR-IA(P0.01)表达量降低。同时Wnt信号通路中的两个基因LEF-1(P0.001)、β-catenin(P0.001)的表达也都显著降低。Western blot结果显示,两条信号通路中这几种蛋白的表达也都受到影响,其中显著升高的有BMP-2(P0.05)、BMP-4(P0.05)、BMP-6(P0.05)和BMP-7(P0.05);显著降低的是β-catenin(P0.05)、BMPR-IA(P0.01)和LEF-1(P0.001)。结论在体外Noggin基因对BMP信号通路可能存在反馈性抑制机制,而对Wnt信号通路存在反馈性激活机制,为下一步探究在体内Noggin基因对BMP和Wnt信号通路表达的作用提供一定的依据。     

球囊菌胁迫中华蜜蜂幼虫肠道过程中病原的转录组学研究

【目的】本研究利用RNA-seq技术对球囊菌胁迫的中华蜜蜂(中蜂)幼虫肠道进行深度测序,经趋势分析得到差异表达基因(DEGs)的显著表达模式,进而对胁迫过程中的球囊菌进行转录组学分析。【方法】利用Illumina HiSeq 2500平台对球囊菌胁迫的中蜂幼虫肠道进行深度测序,并利用相关软件进行了深入分析。最后,通过RT-qPCR对RNA-seq数据进行了验证。【结果】本研究共得到球囊菌的41133932条高质量clean reads。22865个DEGs共聚类为8个基因表达模式,其中,16769个DEGs聚类为2个显著上调趋势与2个显著下调趋势。GO富集分析结果显示,显著上调与显著下调趋势中的DEGs分别富集于40与37个GO term,基因富集数最多的为细胞进程(2486 unigenes)。KEGG代谢通路(pathway)富集分析结果显示,显著上调与显著下调趋势中的DEGs分别富集于119和112个pathway,基因富集数最多的分别是氨基酸生物合成(127 unigenes)与核糖体(98 unigenes)。进一步分析表明球囊菌在胁迫中蜂幼虫肠道的过程中通过提高物质合成促进其增殖,而宿主通过抑制球囊菌的蛋白合成抵御病原入侵。富集在MAPK信号通路的11个DEGs的表达水平随着胁迫时间的延长而逐渐下降,推测中蜂幼虫通过抑制该通路而阻遏球囊菌增殖。【结论】本研究不仅为揭示白垩病过程中的球囊菌-中蜂幼虫互作提供了重要信息,也为阐明不同抗性蜂种的球囊菌抗性差异奠定了基础。     

多个GEO芯片联合分析阿尔茨海默病内嗅皮层的关键基因

运用生物信息学方法筛选阿尔茨海默病(Alzheimer's disease,AD)在内嗅皮层部位的关键基因和通路,探讨AD发生发展的潜在机制.从GEO数据库下载AD的内嗅皮层芯片数据,利用R软件中的sva包和limma包进行批次矫正,并筛选差异基因;运用Cytoscape软件中的ClueGO插件进行GO(Gene On...     

基于蛋白互作网络识别非小细胞肺癌相关基因功能模块

本研究对非小细胞肺癌(non-small cell lung carcinoma,NSCLC)基因表达数据进行差异表达分析,并与蛋白质相互作用网络(PPIN)数据进行整合,进一步利用Heinz搜索算法识别NSCLC相关的基因功能模块,并对模块中的基因进行功能(GO term)和通路(KEGG)富集分析,旨在探究肺癌发病分子机制。蛋白互作网络分析得到一个包含96个基因和117个相互作用的功能模块,以及8个对NSCLC的发生和发展起到关键作用候选基因标志物。富集分析结果表明,这些基因主要富集于基因转录催化及染色质调控等生物学过程,并在基础转录因子、黏着连接、细胞周期、Wnt信号通路及HTLV-Ⅰ感染等生物学通路中发挥重要作用。本研究对非小细胞肺癌相关的基因和生物学通路进行预测,可用于肺癌的早期诊断和早期治疗,以降低肺癌死亡率。     

基于RNA-Seq转录组测序技术揭示肉兔生长发育调控的分子机制

为研究肉兔肌肉生长发育调控的分子机制,本试验以84日龄齐卡巨型白兔和齐兴肉兔为研究对象,屠宰后取背最长肌组织并提取总RNA,反转录建立c DNA文库,利用Illumina HiSeq 2500测序平台进行测序,并进行序列分析和注释。筛选与肉兔肌肉生长发育相关的信号通路及差异表达基因,最后进行荧光定量PCR验证。结果表明:齐卡巨型白兔和齐兴肉兔背最长肌中显著差异表达基因总数为833个,其中齐卡巨型白兔中显著上调基因325个,显著下调基因508个。KEGG功能注释发现差异基因显著富集到PI3K-AKT通路、癌症通路和黏着斑通路,最终筛选出4个可能与生长发育显著相关的差异表达基因,为进一步对肉兔进行遗传改良提供了理论基础。     

未成熟树突状细胞与成熟树突状细胞的细胞骨架调控相关基因表达的差异分析

分析未成熟树突状细胞(immature dendritic cells,im DCs)向成熟的树突状细胞(mature dendritic cells,m DCs)分化的过程中,细胞骨架的调控相关基因及信号通路的表达变化,为进一步理解不同分化阶段树突状细胞(dendritic cells,DCs)的生物物理学特性、迁移能力和免疫学相关功能的改变。利用CEO数据库获得经CD14+单核细胞(monocytos,monos)诱导而成的im DCs和m DCs的m RNA的表达数据(芯片编号:GSE15076),通过R语言软件对原始数据进行处理,筛选差异表达基因(Log|FC|≥2,p0.05),利用STRING online-工具对差异表达基因进一步筛选(可信度≥0.4);Cytoscape软件构建蛋白质相互作用网络图(protein-protein interaction network,PPI),通过Cluster ONE对筛选后的差异表达基因进行聚类分析,筛选功能相关性密切的差异表达基因,并利用DAVID在线分析进一步进行GO分析和KEGG信号通路分析。m DCs相对于im DCs差异表达的基因共3 351个,上调基因1 801个,下调基因1 550个,其中C-C趋化因子受体活性、G-蛋白偶联的嘌呤核苷酸受体活性和异源三聚体G蛋白复合物等与细胞骨架调控密切相关。主要涉及的信号通路包括趋化因子/趋化因子受体信号通路、Rap1信号通路、Jak-STAT信号通路和PI3K-Akt信号通路等,其中的相关基因与细胞骨架的调控关系密切。这对进一步深入理解DCs的生物物理学特性、迁移能力和免疫学功能有一定的帮助。     

Humanin--一种新的神经保护性肽

Humanin(HN)是一种新近发现的神经保护性多肽,能够特异且有效地阻断由家族性阿尔茨海默病(FAD)基因和β淀粉样肽(Aβ)引起的细胞死亡.其作用机制可能是HN通过特异地阻断某些细胞死亡的信号通路来起到对细胞的保护作用.这些研究结果为抗阿尔茨海默病(AD)药物的研发提供了非常有价值的信息.     

基因芯片筛选新疆维吾尔妇女宫颈癌差异表达基因

新疆南部维吾尔族聚居区是宫颈癌高发区. 本文旨在利用基因芯片技术筛选与维吾尔族妇女宫颈癌发生相关的基因. 首先,分别提取5例新疆维吾尔妇女宫颈癌和5例子宫肌瘤组织(对照)的mRNA,逆转录成cDNA,并用Cy3-dUTP标记子宫肌瘤组织的cDNA,用Cy5 dUTP标记宫颈癌组织的cDNA,制成芯片杂交探针.为筛选出宫颈癌组织中差异表达的基因,上述标记探针分别与含有20 000条人类基因的Affymetrix基因芯片进行杂交,杂交信号用GeneChip Scanner 3000扫描仪扫描,并用芯片图像分析软件(SAM software)分析扫描结果.筛选出的差异表达基因经GO(Gene Ontology)分析和KEGG（Kyoto Encyclopedia of Genes and Genomes）信号通路分析,确定其在宫颈癌中的作用.基因芯片筛选结果显示,在宫颈癌组织中发现2 758个差异表达基因,其中1 326个上调基因,1 432个下调基因.GO分析和KEGG信号通路分析表明,表达差异在两倍以上的基因涉及168个信号通路,包括细胞粘附分子、细胞周期以及MAPK和mTOR信号通路等.上述结果表明,基因芯片技术筛选出大量与宫颈癌发生相关的基因,其中表达差异显著的基因涉及细胞粘附分子、细胞周期和mTOR等信号通路.     

西伯利亚蝗卵发育过程的比较转录组分析及滞育关联基因筛选

【目的】通过筛选西伯利亚蝗Gomphocerus sibiricus卵滞育基因及代谢通路,初步明确西伯利亚蝗卵滞育分子机理。【方法】利用Illumina NovaSeq 6000测序平台对西伯利亚蝗早期发育(ES)、滞育(DS)和滞育后发育阶段(PS)的卵进行转录组测序;采用GO和KEGG富集分析并结合文献,预测滞育相关通路并聚类热图分析,筛选蝗卵滞育关联基因,并选取其中6个重要基因进行qRT-PCR验证。【结果】DSvs ES和PS vs DS两比较组分别富集到12 419和4 789个差异表达基因,且多上调表达;两比较组共表达差异基因为2 206个,主要与糖代谢、环境胁迫和生长发育相关。DS vs ES组富集最显著的GO条目为蛋白结合,PS vs DS组GO条目主要包含酶活性、细胞骨架构建及蛋白结合等过程。滞育关联基因主要集中在Wnt信号通路、胰岛素信号通路、细胞周期通路以及昆虫激素生物合成等通路。6个重要的滞育关联基因的表达量变化趋势与转录组数据结果基本一致。【结论】本研究初步明确了西伯利亚蝗卵滞育调控的重要代谢途径,并筛选出20个滞育关联基因,为后续深入研究该物种的适应机理奠定了基础。     

准噶尔小胸鳖甲短时低温胁迫响应的转录组分析

【目的】拟步甲科昆虫小胸鳖甲Microdera punctipennis是分布于中国新疆古尔班通古特沙漠的特有物种,具有很强的耐寒性,但其低温响应的分子机制尚不明确。本研究旨在利用转录组测序技术丰富小胸鳖甲的已知基因信息,分析在短时低温胁迫下,上调表达基因的主要类群及参与的主要代谢通路。【方法】利用Trinity软件对分别在4℃低温和25℃(对照)下处理3 h的小胸鳖甲成虫RNA-Seq数据进行从头组装,利用Blast2go软件进行基因注释。通过DESeq和GFOLD软件分析差异表达基因,并对上调表达基因进行GO(Gene Ontology)和KEGG(Kyoto Encyclopedia of Genes and Genomes)代谢途径富集分析。【结果】测序过滤后得到68 165 001个有效读序,51 712个平均长度为984 bp的非冗余基因,其中29 925个基因(约57.87%)有同源序列(E-value10-5),与拟步甲科的赤拟谷盗Tribolium castaneum相应基因相似性最高(核苷酸序列一致性为72.75%)。低温响应上调表达基因有514个,富集到35个GO类群,18个KEGG途径。其中胁迫响应类群、生物调节类群和免疫系统过程类群都占有重要比例,嘌呤代谢、硫胺素代谢、糖酵解/糖异生等代谢通路显著富集。在胁迫响应类群中,热激蛋白Hsp90基因、超氧化物歧化酶(SOD)基因和钙联接蛋白Calnexin基因等8个非生物胁迫相关基因显著上调表达。【结论】荒漠昆虫小胸鳖甲低温转录组揭示,在4℃冷驯化过程中,代谢过程、胁迫响应、生物调节和免疫系统等生物学过程相关基因表达显著上调,其中几个非生物胁迫响应基因提示小胸鳖甲可能从分子伴侣、细胞周期阻滞、线粒体稳定性、抗氧化胁迫等多方面应对低温胁迫。KEGG富集分析所揭示的小胸鳖甲在低温下的转录组代谢通路与其他物种有较大差异。研究结果为后续生物信息学分析及小胸鳖甲耐寒关键基因的发掘及耐寒机制的揭示提供了基础数据。     

沙葱萤叶甲卵滞育的转录组学分析

【目的】建立沙葱萤叶甲Galeruca daurica滞育卵转录组数据库,挖掘卵滞育相关的基因以及代谢和信号通路,在转录组水平探讨卵滞育的分子机制。【方法】采用Illumina NovaSeq6000高通量测序平台对沙葱萤叶甲滞育卵与解除滞育卵进行转录组测序,并进行生物信息学分析;利用DESeq软件分析沙葱萤叶甲滞育卵与解除滞育卵中的差异表达基因,对差异表达基因进行KEGG通路富集分析;利用qRT PCR技术对10个差异表达基因的表达模式进行验证。【结果】基于沙葱萤叶甲滞育卵与解除滞育卵转录组测序结果,共获得53 389个unigene,其中差异表达基因2 145个,24个差异表达基因与保幼激素信号及脂肪酸生物合成和降解相关。与解除滞育卵相比,滞育卵转录组中1 297个基因上调表达,富集于124条KEGG通路,其中核糖体通路显著富集;848个基因下调表达,富集于73条KEGG通路,其中MAPK信号通路和糖胺聚糖生物合成通路显著富集。qRT-PCR结果表明,随机选取的10个差异表达基因的表达趋势与RNA-Seq转录组测序结果完全一致。【结论】保幼激素,脂肪酸生物合成和降解,核糖体,MAPK信号及糖胺聚糖生物合成等通路可能在沙葱萤叶甲卵滞育调控中起着重要的作用。     

基于转录组的八角挥发油合成相关基因挖掘与分析

为更好地挖掘八角(Illicium verum)挥发油合成相关基因，该文对挥发油性状差异显著的优良无性系桂角69号及普通品种砧01号叶片进行了转录组测序及组装注释，并对差异表达基因进行了GO分类和KEGG通路分析。结果表明：(1)转录本经组装后获得84 182条序列，使用NR、NT、Swiss-Prot、KEGG、KOG、GO和Pfam数据库进行序列比对，共注释了59 161条序列，筛选出30 572个差异表达基因。与砧01号相比，桂角69号叶片中上调基因有15 025个，下调基因有15 547个。(2)GO分类结果显示共有20 287个差异基因被注释。KEGG分析结果表明，有21 600个差异基因被注释到133条KEGG通路上，其中挥发油合成相关的单萜生物合成通路、萜类骨架生物合成通路、苯丙素合成通路中的芳樟醇合酶、月桂烯合酶、香叶基香叶基焦磷酸合酶、肉桂酰辅酶A还原酶、咖啡酸3-O-甲基转移酶、肉桂醇脱氢酶等关键酶基因呈差异表达。(3)转录因子分析发现差异表达基因分布于31个转录因子家族，其中MYB家族序列数量最多。该文利用转录组测序技术分析八角优良无性系与普通品种叶片的差异基因及...     

急性髓系白血病mRNA与非编码RNA差异表达谱及CeRNA调控网络分析

利用GEO数据库(gene expression omnibus database)通过生物信息学分析方法探讨急性髓系白血病(acute myelogenous leukemia,AML)的发病机制。检索GEO数据库中AML相关芯片数据集GSE142698、GSE142699和GSE96535。利用GEO2R分析得到差异mRNAs、miRNAs以及差异lncRNAs。利用在线生物信息学分析工具DAVID对差异mRNAs进行GO富集分析和KEGG通路分析。利用miRWalk数据库预测AML相关miRNAs的靶向mRNAs,利用Spongescan数据库预测AML相关miRNAs的靶向lncRNAs,构建lncRNA-miRNA-mRNA竞争性内源RNA （competing endogenous RNA,ceRNA）调控网络。共筛选出29个显著差异mRNAs、70个显著差异miRNAs和20 005个显著差异lncRNAs。GO富集分析和KEGG通路分析显示,差异表达基因主要涉及蛋白磷酸化、细胞分裂、细胞增殖的负调控、基因表达的正向调节、周期蛋白依赖的丝氨酸/苏氨酸激酶活性的调节等生物过程以及细胞周期、细胞衰老、癌症通路、PI3K-Akt通路等信号通路。将miRWalk数据库预测的靶向mRNAs与差异mRNAs取交集,Spongescan数据库预测的靶向lncRNAs与差异lncRNAs取交集,分别确定了25个mRNAs、6个lncRNAs参与AML相关ceRNA调控网络的构建。结果表明,lncRNAs可能作为关键的ceRNA,通过调控miRNA和相关靶基因参与AML的发生与发展,研究结果为AML诊断和治疗的分子生物学研究提供了新的依据。     

根据在癌突变谱中的共突变基因对识别癌基因

癌是一种涉及多基因变异的遗传异质性疾病,涉及多种生物学功能通路中不同基因的遗传变异。因此,识别癌基因是一项富有挑战性的工作。提出通过寻找在癌样本中突变显著共发生的基因筛选候选癌基因的方法。应用该方法,通过分析蛋白激酶基因在癌组织中的突变谱数据,发现了167个显著共发生突变的基因对,包含85个基因。分析这167个基因对发现:(1)发生共突变的基因富集已知的癌基因;(2)共突变基因对倾向于共扰动与癌症相关的通路对。以上结果提示,在癌样本中显著共发生突变的基因倾向于候选癌基因;在癌发生过程中起重要作用的基因倾向于协同扰动不同的癌相关细胞生物学过程。     

MDCK细胞和MDCK-G1细胞感染H1N1流感病毒后病毒滴度和细胞转录组表达差异分析

目的通过分析细胞表达基因及其相关信号通路的差异等信息,探索MDCK和MDCK-G1细胞株在接种H1N1流感病毒后出现病毒滴度和细胞生长趋势差异的内在生物学特性等原因。方法通过Illumina测序平台对2株细胞进行转录组测序(RNA sequencing, RNA-seq)和测序结果生物信息学分析后,选择了17个差异基因进行实时定量聚合酶链反应(quantitative real-time polymerase chain reaction, qRT-PCR)验证。结果 MDCK和MDCK-G1细胞差异基因(differentially expressed genes, DEGs)|log_2(FoldChange)|0和padj≤0.05总数为2 786个。相比于MDCK细胞,MDCK-G1细胞有967个基因的表达显著上升,1 819个基因的表达显著下降。差异基因通过基因本体(Gene Ontology, GO)数据库富集功能注释和京都基因和基因组数据库(Kyoto Encyclopedia of Genes and Genomes, KEGG)进行通路分析发现,2株细胞在免疫调控和细胞生长周期调控上均存在差异,qRT-PCR验证的17个基因中有16个基因的表达趋势与转录组测序结果一致。结论 2株细胞的转录谱存在显著差异。