COPD 差异表达基因的生物信息学分析及在LUSC样本中的表达

为确定慢性阻塞性肺病(COPD)的分子标记物及COPD与肺鳞状细胞癌(LUSC)共存的差异表达基因,探寻COPD合并肺癌的预测因子,发现新的治疗靶点。本研究采用生物信息学方法,从GEO数据库中筛选3套基因芯片数据集,挖掘COPD患者小气道上皮细胞(SAEC)的差异表达基因(DEG)以及潜在的生物标记物,并通过基因本体(GO)、京都基因与基因组百科全书(KEGG)富集分析预测DEGs的功能及参与的代谢途径。继而对DEGs构建PPI网络,使用Cytoscape软件筛选子模块和Hub基因,并将Hub基因通过TCGA数据库分析其在LUSC中的差异表达情况及差异基因间的相关性。结果共获得52个上调基因和24个下调基因,代谢通路主要集中在细胞色素P450对外源物质的代谢、化学致癌、花生四烯酸代谢及甲状腺激素合成四条途径上,通过Cytoscape软件从PPI网络中筛选得到2个功能模块和10个Hub基因,进一步验证发现其中5个基因在TCGA数据库中的LUSC样本中同样差异表达。由此推测SPP1、ALDH3A1、SPRR3、KRT6A和SPRR1B 可能为COPD 分子标记物及COPD与LUSC共存的DEGs,从而为研究COPD和LUSC的发病机制及二者潜在关系奠定良好的基础。     

单核细胞与未成熟树突状细胞粘附相关差异表达基因筛选及相互作用分析

本研究目的是分析单核细胞(monocytes,mon)与未成熟树突状细胞(immature DCs,im DCs)两个不同发育阶段间差异表达基因的功能及其编码蛋白的相互作用,筛选出粘附相关的关键基因,并进行生物信息学分析。我们从NCBI(美国国立生物技术信息中心)公共数据平台GEO(Gene Expression Omnibus)下载基因芯片数据GSE15076,使用perl语言将探针id转换成gene symbol。利用R Bioconducto软件(RGui 3.3.1)解析原始数据并筛选差异基因。进一步利用STRING online工具筛选核心差异基因(可信度≥0.4),Cytoscape软件构造蛋白质相互作用网络图。对STRING online工具筛选核心差异基因进行京都基因与基因组百科全书(kyotoencyclopedia of genes and genomes,KEGG)分析。im DCs相对于mon表达的差异基因,过滤条件为Log FC2和adj.p.val0.01。通过分析,我们共找到1 223个差异基因,其中567个上调基因,656个下调基因,通过进一步筛选,在网络中的上调基因有399个,下调基因458个,主要涉及m TOR信号通路、细胞代谢、细胞粘附等功能。其中上调基因中CDH1,CD274等差异基因与细胞粘附密切相关,而下调基因中SELL、IL-6、TNF等差异基因与细胞粘附改变密切相关。因此,在单核细胞发育分化到未成熟树突细胞阶段,粘附相关基因的表达变化,对进一步深入理解im DCs独特的生物物理学特性和免疫学功能来说具有重要意义。为进一步研究DCs细胞粘附功能的分子机制提供指导。     

细胞凋亡和代谢基因在转移倾向结肠癌中的研究

目的:应用基因微阵列技术初步筛选与不同转移倾向结肠癌相关的细胞凋亡和代谢相关基因,研究转移相关基因功能.方法:取结肠癌肝转移和无转移结肠癌组织,采用人全基因组表达谱芯片获得两组织的基因表达谱,分析比较两者之间细胞凋亡和代谢基因的差异表达情况;利用基因数据库检索结肠癌相关基因,分析基因功能.结果:应用含有16450个克隆(其中3869个未知)的cDNA微阵列分析发现,细胞凋亡或肿瘤相关基因中,2倍以上(Ratio值小于0.5或大于2.0)差异基因共216个,上调基因85个,下调基因129个.表达差异5倍以上(Ratio值小于0.2或大于5.0)共32个,上调基因10个,下调基因22个.在细胞代谢相关基因中,2倍以上(Ratio值小于0.5或大于2.O)差异基因共205个,上调基因86个,下调基因119个.表达差异5倍以上(Ratio值小于0.2或大于5.0)共15个,上调基因10个,下调基因5个.利用基因数据库检索分析发现5个基因与结肠癌转移关系密切.结论:结肠癌的发生和转移是多基因参与的,本实验应用基因微阵列技术发现细胞凋亡和代谢相关基因中发现5个基因与结肠癌转移关系密切.     

黄秋葵花和果荚转录组测序及类黄酮代谢差异表达分析

采用Illumina Hi seq 2500转录组高通量测序,构建黄秋葵花和果荚转录组文库,并利用测序评估、差异基因功能注释等生物信息学方法进行分析。研究结果表明：(1)分别获得黄秋葵花和果荚有效数据7.12 Gb和8.14 Gb,碱基百分比(Q30)均达到91.0%以上。(2)获得差异表达基因(DEGs)1 336个,其中上调基因319个,下调基因1 017个。(3)获得功能注释的基因有1 131个,GO将455注释基因归成41个功能小类,主要涉及代谢过程、催化活性、单一生物过程和细胞过程等过程;KOG注释了472个DEGs,涉及23个功能分类,其中与次生代谢直接相关过程O和Q类别获得111个注释结果;KEGG将372个DEGs注释到80条代谢通路中,获得F3H、F3′5′H、DFR、ANR、ANS、GT、LAR共10个关键差异基因,其中F3H、DFR在黄秋葵花中表现上调效应,F3′5′H、ANR、LAR在黄秋葵果荚中表现显著上调效应,ANS、GT则分别在花和果荚中均有上调或下调效应。(4)实时定量PCR分析显示,其中7个差异表达基因得到的相对表达量与转录组表达谱分析趋势完全一致。(5)类黄酮代谢途径分析表明,黄秋葵花通过F3H、DFR、ANS、GT途径将NAR催化为生成花青素苷;果荚则将NAR通过F3′5′H将催化为DHM,后在FLS催化下生成黄酮醇类物质等;部分NAR在F3H、 DFR催化下,生成无色飞燕草素苷元,再分别在ANS、LAR作用下,进入花青素苷元和原花青素合成途径。该研究结果丰富了黄秋葵转录组信息,为黄秋葵类黄酮物质纯化和功能利用提供参考依据。     

基于GEO和TCGA数据库分析促癌基因INHBA和抑癌基因CLCA4、CA4在结直肠癌中表达

目的:通过分析GEO数据库结直肠癌相关芯片集,寻找差异基因,并在TCGA数据库和GEO数据库进行验证,为结直肠癌的早期诊断寻找标志物。方法:分析GEO数据库结直肠癌相关芯片集GSE21510、GSE25071、GSE32323。分别分析差异基因,采用文恩图软件查找共同差异基因。进一步在TCGA数据库查找差异基因在结直肠癌中的表达及生存曲线。最后通过GEO数据库GSE24514验证差异基因的表达。结果:GSE21510,包含104例样本,共筛选出251个差异基因,其中上调基因146个,下调基因105个。GSE25071,包含50例样本,共筛选出669个差异基因,其中上调基因312个,下调基因357个。GSE32323,包含10例样本,共筛选出353个差异基因,其中上调基因115个,下调基因238个。在样本中上调基因为促癌基因,下调基因为抑癌基因。经文恩图分析,3个基因集交集共有15个基因,其中上调基因3个,下调基因12个。在TCGA数据库中查找差异基因的表达量和生存曲线,生存曲线选择结肠癌数据集,选取279个样本进行分析。根据差异基因的表达和生存曲线,最终确定促癌基因INHBA和抑癌基因CLCA4、CA4为结直肠癌的标志物。最后在GSE24514芯片集验证差异基因的表达。结论:通过GEO和TCGA数据库筛选及验证,发现在结直肠癌组织中INHBA基因明显上调,CLCA4、CA4基因明显下调。最终确定促癌基因INHBA和抑癌基因CLCA4、CA4可作为结直肠癌早期诊断的标志物。     

茶树种子发育过程的转录组分析

该研究以‘铁观音’茶树品种的种子为试验材料,采用转录组测序技术分析种子发育的3个时期(幼果期、膨大期、成熟期)的表达差异,探究茶树种子油脂代谢的分子机制。结果表明:(1)经转录组测序、组装后共获得30 940 581个clean reads,经数据合并拼接最终得到36 951条非冗余Unigene序列,其中28 476个Unigene可得到功能注释;在转录本中能够被注释到GO分类的Unigene有11 201条(30.3%),KEGG分析发现共有17 172个基因参与了127个代谢通路。(2)经KEGG通路筛选出14条与脂肪酸代谢相关的通路,且随着茶籽的发育,大部分脂肪酸调控途径相关基因呈下调趋势,其中上调基因数最多的有α-亚麻酸代谢途径和脂肪酸降解途径(有17个基因表达量上调),下调基因数最多的是甘油磷脂代谢途径(有58个基因表达量下调);在茶籽发育幼果期α-亚麻酸代谢途径中表达量上调的基因数超过表达量下调的基因数。(3)研究发现茶籽脂肪酸合成相关的基因涉及14个脂类调控途径,共409条差异基因;随着茶树种子发育到成熟期,上调的差异表达基因数量在减少,下调的差异表达基因数量增加,其中α-亚麻酸途径中的基因PLA2G16、DAD1、pldA、FabF、FabI表达量上调显著,随后表达量下调。(4)qRT-PCR检测结果表明,7个茶树FAD和1个ACP差异表达基因的水平与转录组测序结果基本一致;随着茶籽的发育,基因CsFAD7和Δ6-CsFAD从幼果期、果实膨大期至果实成熟期都为差异下调表达,CsFAD2、CsFAD6和Δ7-CsFAD为差异上调表达,CsFAD8、Δ8-CsFAD和CsACP在幼果期至果实膨大期差异上调表达,在果实膨大期至果实成熟期差异下调表达。     

基于转录组测序的花烟草低温胁迫响应转录因子挖掘

为了鉴定参与花烟草低温胁迫的转录因子,对花烟草幼苗进行了4℃低温处理,并在处理后12 h采集其幼苗样本,提取总RNA后,采用高通量测序技术,进行了转录组测序。在对差异表达基因进行GO及KEGG分析的基础上,对参与其中的转录因子进行了挖掘,并采用qRT-PCR的方法,对转录组测序的结果进行了验证。结果表明,低温处理后,花烟草基因表达量变化在2倍以上的基因有8388个(P<0.01),其中,上调表达4229个,下调表达4159个。这些差异表达基因的功能归类于生物过程、细胞组分及分子功能3大类69个GO条目,并显著富集在40条KEGG代谢通路中。同时,在低温胁迫下,花烟草有118个转录因子的表达发生了显著改变,其中,上调表达82个,下调表达36个。这些转录因子属于28个家族,其中,数量最多的为NAC家族19个,其次为ERF家族16个、MYB家族15个、WRKY家族15个。本研究的结果为花烟草低温响应分子机制的研究提供了借鉴。     

单核细胞与未成熟树突细胞的细胞骨架相关的核心差异基因表达的生物信息学分析

为进一步深入理解不同分化阶段树突状细胞(dendritic cells,DCs)的生物物理学特性和免疫学功能的差异。从CEO数据库获得CD14~+单核细胞(monocytes,mon)及其诱导而成的未成熟DCs(immature DCs,im DCs)的m RNA表达数据(芯片编号:GSE15076),利用R语言软件(RGui 3.2.3)解析原始数据并筛选差异基因。进一步利用STRING online工具筛选核心差异基因(可信度≥0.4),Cytoscape软件构造蛋白质相互作用网络图。对STRING online工具筛选核心差异基因进行京都基因与基因组百科全书(kyoto encyclopedia of genes and genomes,KEGG)分析。实时荧光定量PCR技术检测mon及im DCs部分核心差异基因CCNB1、CDK5、IL-6和TNF在基因水平上的表达情况。im DCs相对于mon表达的差异的基因共3 371个(im DCs与mon对比,表达变化倍数≥2,p值0.05),其中上调基因为2 396个,下调基因为975个,通过进一步筛选后获得上调基因655个,下调基因110个,主要涉及m TOR信号通路、细胞代谢、细胞粘附等功能。其中上调基因中CDK5、CCNB1等差异基因与细胞骨架改变密切相关,而下调基因中IL-6、TNF、CXCR4等差异基因与细胞骨架改变密切相关,这对进一步深入理解im DCs独特的生物物理学特性和免疫学功能来说具有重要意义。     

烟草脉斑驳病毒侵染本氏烟的转录组研究

为了研究烟草脉斑驳病毒（TVMV）侵染本氏烟的抗病分子机制,本文建立了TVMV转录组数据库,挖掘抗病相关基因以及代谢和信号通路,采用Illumina Novaseq 6000高通量测序平台对侵染TVMV的本氏烟进行转录组测序,并进行生物信息学分析;利用R package:edge R等软件分析了有关抗病的差异表达基因,对差异表达基因进行基因本体论（GO）及京都基因和基因组百科全书（KEGG）基于途径的通路富集分析。基于TVMV侵染本氏烟转录组的测序结果中共获得4 593个差异表达基因,其中3 564个基因上调表达,1 029个基因下调表达。共筛选出10个抗病相关的差异表达显著基因,6个上调,4个下调。GO数据库中注释到生物学过程、细胞组分及分子功能等三大类共50个功能组。其中,在第三大类分子功能中,蛋白质结合功能以及ATP结合功能所涉及的差异表达基因最显著,分别为501个和453个差异基因。KEGG共富集20条通路,注释到核糖体途径的差异基因富集程度最高,基因最显著,差异表达基因数为307个。蛋白质结合、ATP结合、核糖体、真核生物核糖体的生物反应、DNA复制（DNA replicat...     

基于高通量测序的都匀地区福鼎大白种茶树根茎叶分析

为探究茶树中茶多酚等产物代谢途径的相关基因,该研究以贵州都匀地区福鼎大白种茶树的根茎叶为对象,利用高通量测序技术构建茶的转录组数据库并筛选其根茎叶差异表达基因。结果表明:共获得70.88 Gb Clean Data,各样品Clean Data均达到6.33 Gb,Q30碱基百分比在93.22%以上。将Clean Reads与中国种茶树参考基因组进行序列比对,比对效率从87.83%到91.14%。基于比对结果,进行可变剪接预测分析和基因结构优化分析,发掘新基因13 531个,其中10 244个得到功能注释。利用FPKM进行基因表达量分析,根据基因在不同样品中表达量识别差异表达基因。叶与茎的差异基因有5 595个,其中2 769个在茎中上调,2 826个下调,叶与根有9 650个差异基因,5 056个上调,4 594个下调,茎与根中有5 644个差异基因,2 938个上调,2 706个下调,并通过GO和KEGG分析,将差异基因进行功能注释和富集分析。上述结果为揭示都匀地区福鼎大白种茶参与类黄酮、茶氨酸和咖啡碱等代谢途径相关的基因提供了参考,为选育优良品种等提供了理论依据。     

利用RNA-seq技术分析淹水胁迫下转BnERF拟南芥差异表达基因

为探究淹水胁迫下BnERF调节的耐淹防御相关途径,应用RNA-seq技术,对淹水6小时后的拟南芥(Arabidopsis thaliana)野生型(WT)和转BnERF株系(E33)幼苗进行基因表达分析。结果表明,淹水3天后,E33表现出较强的耐淹性,地上部生长状况和根系发育均明显强于野生型。E33幼苗未淹水处理时相对于野生型单独上调的基因有9个,4个为膜结合蛋白,其中2个参与MAPK级联途径,其它5个参与氧化胁迫及水分调节途径;与未淹水野生型相比,无论是未淹水处理还是淹水6小时后的E33幼苗中缺氧响应、抗氧化防护及细胞、器官发育相关基因的表达量均上调。另外,淹水6小时后E33的差异基因并未完全覆盖淹水6小时后野生型的差异基因;E33幼苗中缺氧响应、氧化胁迫响应、能量的产生与转变、乙醇代谢途径中的基因以及乙烯响应因子基因的表达量都明显高于野生型。上述结果表明,BnERF直接或间接调节植物的淹水胁迫相关生理代谢途径,参与淹水胁迫的防御过程。     

高油酸花生发芽期低温胁迫转录组及差异表达基因分析

高油酸花生(ArachishypogaeaL.)油具有优异的营养成分和热氧化稳定性,有利于人体健康和工业生产。但是高油酸花生在发芽期间对温度比较敏感,限制了其在低温地区的引种。为了进一步了解高油酸花生在发芽期响应低温胁迫的分子机制,本研究选用田间试验中耐寒表现不同的4种高油酸花生品种,分析其在发芽期低温胁迫下的全基因组水平调控。通过转录组高通量测序共获得139429条Unigene,其中两组耐寒与不耐寒高油酸花生品种在低温胁迫下共产生差异表达基因(differentially expressed gene, DEG) 3520个,且耐寒花生中上调表达的DEG数目大于下调表达的数目。GO分析表明,耐寒高油酸花生中有关细胞膜代谢与完整性以及细胞外周蛋白差异表达基因的数量较多;KEGG通路分析表明,植物病原相互作用和植物激素信号转导通路在抗寒中起着重要作用。进一步筛选出4个低温诱导蛋白基因——时钟调控蛋白基因(TIC)、AGO4蛋白基因(AGO4)、组蛋白–赖氨酸N甲基类转移酶ATX3基因(ATX3)、FERONIA类受体蛋白激酶基因(FER)和3个转录因子基因——bHLH、3R-1MYB和EREB,采用qRT-PCR检测这些基因在低温胁迫下的表达量。结果显示,在低温处理3h后TIC、ATX3和AGO4以及转录因子基因bHLH、MYB和EREB的表达量明显升高,FER在胁迫12h后也有明显升高,表明这些基因在花生萌发期响应低温胁迫。本研究为深入了解高油酸花生发芽期低温胁迫转录调控机制和筛选花生抗寒基因提供了数据资源。     

小麦叶锈菌休眠与萌发夏孢子的差异表达

【目的】小麦叶锈菌引起的小麦叶锈病是小麦的重要病害之一,孢子萌发和侵入气孔是保证叶锈菌正常侵染寄主的重要前提。本研究旨在研究小麦叶锈菌夏孢子萌发的差异表达特性,为揭示萌发的机制及其与致病的关系提供依据。【方法】利用小麦叶锈菌致病类型THTT的休眠夏孢子和萌发夏孢子进行RNA-seq分析。【结果】在萌发夏孢子中筛选出相比休眠夏孢子上调表达的基因为4400个,下调基因5325个。GO富集分析发现,上调差异基因主要涉及细胞进程、有机物代谢、信号传导、单个有机体过程、催化酶活性等;下调基因主要涉及单组织过程、单个有机体细胞过程、有机物代谢过程、催化活性调节等过程。利用KEGG分析发现,差异基因共参与109条通路,从中筛选出两条与孢子萌发相关的通路——丝裂原活化蛋白激酶(MAPK)信号通路和囊泡运输中的SNARE蛋白交流。10个基因的定量分析结果与基因数字表达谱分析结果一致。【结论】研究获得了叶锈菌孢子萌发及侵入气孔过程中的重要差异基因,研究结果为研究叶锈菌致病机制奠定基础。     

胃癌差异表达基因在染色体上的定位及其功能分析

利用标准化的Affymetrix公司生产的U133A基因芯片检测胃癌（T）与切缘正常胃黏膜（C）基因表达谱差异,并利用生物信息学方法对检测结果进行差异基因在染色体定位和功能分析。结果表明：胃癌与正常胃黏膜比较差异8倍以上共有270个基因,其中表达上调[信号比的对数值（SLR）≥3]有157个,表达下调（SLR≤-3）有113个。从表达差异的基因在染色体定位分析,发现除4个基因未知其定位外,其余所有差异表达基因散在分布和各条染色体上,但以1号染色体为最多,有26个（占9．8％）,其次是11和19号染色体上分别有24个（各占9．1％）。而差异表达的基因发生在染色体短臂（q）上有173个（占65％）。从表达差异的基因功能分类看,属于酶和酶调控子基因最多（67个,24．8占％）,其次是信号传导基因（43个,占15．9％）,第3类是核酸结合基因（17个,占6．3％）,第4类是转运子基因（15个,占5．5％）,第5类是蛋白结合基因（12个,占4．4％）,还有功能未知的基因有50个,占18．5％。以上5大类共占基因总数56．9％。胃癌差异表达基因散在分布在各条染色体上,但以1、11、19号染色体差异表达基因居多。这5大类（酶和酶调控子、信号传导、核酸结合、转运子、蛋白结合）相关基因异常是今后研究胃癌的重要基因。     

金针菇单、双核菌丝差异表达基因分析

对金针菇Flammulina velutipes单核菌丝W23的菌丝体以及与L11质配后的双核菌丝H1123菌丝体进行了转录组测序,以本实验室已获得的W23基因组为参考基因组研究两样本间差异基因,并对这些差异基因进行了GO功能和Pathway显著性富集分析。差异基因分析显示,两个样本中共有显著性差异表达的基因3 504个,其中在双核菌丝中上调、下调的基因数分别为2 151和1 353个。研究发现差异表达基因含有很多的转录因子基因、蛋白激酶以及WD40 repeat-like蛋白。Gene Ontology（GO）功能分析结果表明,extracellular region和membrane-enclosed lumen条目下的差异基因全部为上调表达,而envelope下的差异基因全部为下调表达,以利于双核菌丝分裂时锁状联合的形成而便于核的迁移。Pathway 功能富集分析结果表明,脂肪酸、氨基酸以及大部分糖类合成相关基因具有比较活跃的上调表达。说明双核菌丝主要进行营养物质的富集,为下一步在合适条件下分化成原基,进入生殖生长阶段储备物质基础。     

用基因芯片筛选前列腺癌LNCaP和C4-2细胞系中差异表达的基因

目的：筛选与前列腺癌进展相关的功能基因。方法：利用Affymetrix公司的人类金基因组U133A芯片,筛选在前列腺癌细胞系LNCaP和C4-2中差异表达的基因,并用RT-PCR方法验证部分差异基因的表达。结果：芯片筛选结果表明,与LNCaP细胞系相比,C4-2细胞系中217个基因转录本表达上调,101个基因转录本表达下调;对其中45个基因进行了RT-PCR验证,证明35个基因转录本的表达结果与芯片筛选一致。结论：筛选出了在LNCaP和C4-2细胞系中显著差异表达的基因35个,为进一步研究前列腺癌进展的机制奠定了基础。     

转录组和代谢组分析瘤菌根菌对铁皮石斛的促生机制

为筛选出促进铁皮石斛(Dendrobium officinale)生长的差异表达基因和差异代谢物,对瘤菌根菌与无菌盆栽铁皮石斛苗共生后形成的侧根根系进行转录组、代谢组和双组学联合分析。结果表明,转录组分析共找到262条差异表达基因富集到了35条通路中,其中内质网蛋白质加工通路途径的差异基因最多,其次为氨基糖和核苷酸糖代谢通路。代谢组分析共检测出194个差异代谢物富集到33个KEGG通路中,其中代谢途径的差异代谢物最多有133个,其次为不同环境的微生物代谢途径的差异代谢物有70个。通过联合分析,有9个差异基因的差异表达导致丝氨酸、谷氨酸、D-甘露糖和激素等代谢物的积累量发生变化,这可能是瘤菌根菌促进铁皮石斛生长的重要原因。因此,推测瘤菌根菌促进铁皮石斛生长与氨基酸、糖、植物激素的积累及相关基因的表达变化有关。     

基于高通量转录组测序的草鱼雌雄性腺差异表达基因分析

草鱼雌性个体生长优势明显高于雄性,挖掘精巢和卵巢差异表达的功能基因对草鱼生殖调控具有重要的价值及应用前景。采用Illumina HiSep 2500转录组测序技术分别对12月龄的3尾雌性草鱼和3尾雄性草鱼的性腺组织的mRNA进行测序分析,精巢和卵巢共得到8 363个差异表达基因。与卵巢相比,精巢上调基因6 446个(77%),下调基因1 917个(23%),差异明显。将上述差异基因在Nr、GO和KEGG数据库进行比对,获得48种Go功能注释分类和313条KEGG代谢通路。分析发现,Go功能分类中"绑定"功能富集差异基因的数量最多,占总数的15.8%;KEGG代谢通路中"信号转导"通路富集差异基因的数量最多,占总数的10.2%。结合转录本数据,挑选与草鱼性腺发育相关的关键差异基因Dmrt1、amh、foxl2、cyp19a1a进行分析,研究显示,与卵巢相比,Dmrt1和amh基因在草鱼精巢中极显著高表达(p0.01);与精巢相比,cyp19a1a和foxl2基因在草鱼卵巢中显著高表达(p0.05),与转录组数据分析结论一致,说明Dmrt1和amh基因与早期草鱼精巢发育调控相关,cyp19a1a和foxl2基因与早期草鱼卵巢发育调控相关。本实验为进一步了解草鱼性别基因的调控机理提供数据依据。     

6天龄与10天龄大鼠的睾丸组织的基因差异表达

通过基因芯片技术,利用Roche-NimbleGen公司制作的大鼠12×135K全基因组表达谱芯片,对日龄为6d和10d的大鼠睾丸组织进行全基因组表达差异分析。结果显示:具有2倍以上的差异表达基因有4298个,其中表达上调的基因共1878个,表达下调的基因共2420个。这些差异表达的基因中有3154个基因具有基因本体注释,参与了154个生物学通路。进一步分析表明具有8倍以上差异表达的基因有13个,这些基因参与了生物学过程、细胞组分和分子功能等基因本体分类,进一步选择3个差异表达的基因,LOC686076、Cxcl6和Trib3,做了实时定量RT-PCR检测。其结果趋势与芯片数据一致。因此,我们初步认为精原干细胞的发生与增殖在大鼠早期的发育过程中已经有大量的基因参与,是一个多基因协调表达的过程。