基于转录组测序分析壮骨止痛方治疗去势大鼠骨质疏松症的基因表达差异

目的:基于转录组测序分析壮骨止痛方治疗去势大鼠骨质疏松模型的差异表达基因及相关通路.方法:将12只SD雌性大鼠随机分为实验组、模型组,每组各6只,采用国内外公认的雌性大鼠去卵巢3个月造成绝经后骨质疏松模型.造模后1周开始灌胃给药,实验组给予壮骨止痛胶囊(567 mg/kg),模型组给予同体积的蒸馏水(1 mL/150 ...     

基于转录组学的补骨脂酚致小鼠肝损伤种属差异

目的 探究补骨脂酚在Institute of Cancer Research（ICR）小鼠和昆明（KM）小鼠不同种属之间的肝毒性应答差异。方法 通过急性毒性和亚急性毒性动物实验证实补骨脂酚致小鼠肝毒性的客观表现,运用转录组学比较正常ICR小鼠和KM小鼠2种不同种属小鼠间的差异基因富集通路,在此基础上采用实时荧光定量聚合酶链式反应（Real-time PCR）验证相关通路上关键基因的mRNA表达证实补骨脂酚致肝损伤的种属差异。结果 补骨脂酚亚急性毒性实验,与正常组比较,ICR小鼠血清碱性磷酸酶（ALP）,5′-核苷酸酶（5′-NT）含量均有不同程度升高,差异无统计学意义,KM小鼠无明显变化;病理结果显示ICR小鼠以肝细胞肥大为主要病理特征,KM小鼠主要表现为肝细胞脂肪变性。补骨脂酚急性毒性实验,给药3 d,KM小鼠出现竖毛,精神萎靡,濒死,ICR小鼠无明显毒性表现。与正常组比较,KM小鼠（400 mg·kg^-1）血清丙氨酸氨基转移酶（ALT）,天门冬氨酸氨基转移酶（AST）显著升高（P<0.01）,肝脏总活性氧（SOD）酶活力显著降低（P<0.01）;ICR小鼠（400 mg·kg^-1）血清5′-NT,胆碱酯酶（CHE）含量显著升高（P<0.01）。ICR小鼠肝脑比增加20.34%,KM小鼠增加29.14%（P<0.01）。ICR小鼠肝脏病理主要表现为肝细胞肥大,KM小鼠以局灶性炎症为主,伴肝细胞肥大、肝细胞脂肪变性。转录组学京都基因和基因组数据库（KEGG）通路和Reactome通路富集分析,发现ICR小鼠和KM小鼠不同种属间差异基因表达主要涉及氧化磷酸化、胆汁分泌、胆汁酸及胆汁酸盐合成和代谢通路。与正常组比较,细胞色素P450家族成员7A1（CYP7A1）在各给药组均显著上调（P<0.01）,多耐药相关蛋白（MRP2）在KM小鼠各给药组均显著降低（P<0.01）,在ICR小鼠各给药组显著升高（P<0.01）;胆汁酸盐输出泵（BSEP）在ICR小鼠急性肝损伤（400 mg·kg^-1）组明显降低（P<0.05）;SHP1在KM小鼠急性肝损伤（400 mg·kg^-1）组高表达,法尼醇X受体（FXR）在ICR小鼠亚急性肝损伤（200 mg·kg^-1）组显著降低（P<0.01）;SOD1,过氧化氢酶（CAT）,核转录因子E₂相关因子2（NFR2）在KM小鼠急性肝损伤（400 mg·kg^-1）组显著降低（P<0.01）,CAT在KM小鼠亚急性肝损伤（200 mg·kg^-1）组显著降低（P<0.01）;谷胱甘肽巯基转移酶1（GSTA1）,谷胱甘肽过氧化物酶1（GPX1）在KM小鼠急性肝损伤（400 mg·kg^-1）组显著增强（P<0.01）,SOD1,CAT,NRF2,GSTA1在ICR小鼠亚急性肝损伤（200 mg·kg^-1）组显著升高（P<0.01）,CAT,NRF2在ICR小鼠急性肝损伤（400 mg·kg^-1）组显著升高（P<0.01）。结论 随着补骨脂酚给药剂量增加,KM小鼠氧化应激所致肝损伤逐渐加重;ICR小鼠表现出更强的抗氧化能力。通过比较2种小鼠对补骨脂酚的毒性应答反应,关于研究补骨脂酚所致小鼠肝脏肝毒性机制,ICR小鼠更适合研究胆汁分泌及胆汁酸代谢相关的机制,KM小鼠更易受氧化应激的影响而发生肝损伤。     

基于代谢组学和转录组学的不同生长年限下银杏萜类生物合成关键基因表达分析

目的 以不同发育阶段不同树龄的银杏叶（银杏幼树叶,多年生银杏树叶）为对象,分析次生代谢物变化规律,测定黄酮及萜内酯含量,分析萜类生物合成的通路和关键基因表达,为提高银杏萜内酯产量提供分子生药学数据。方法 采用液相色谱-质谱联用（LC-MS）技术,并进行高通量转录组学测序（RNA-seq）,从差异表达基因中寻找次生代谢物生物合成的关联酶基因。结果 代谢数据表明,不同树龄,不同发育阶段的银杏叶次生代谢物明显不同,黄酮类成分,老树叶片与幼树叶4～7月含量均呈下降趋势;同发育阶段,幼树叶含量高于老树叶片含量;萜内酯类成分,老树叶片4～7月含量略微呈上升趋势,幼树叶略微呈缓降趋势;同发育阶段,幼树叶含量高于老树叶片含量。转录组分析出萌发期,老树叶片与幼树叶片生长状态相似,5～7月老树叶与幼树叶基因表达发生较大变化。筛选出49条与萜类合成有关的候选基因,分析后发现甲羟戊酸途径（mevalonatepathway,MVA）中老树叶表达高于幼树叶;甲基赤藻糖醇磷酸途径（methylerythritol4-phosphate pathway,MEP）中幼树叶表达高于老树叶。结论 幼树叶比老树叶药用活性成...     

基于转录组测序的银杏类黄酮生物合成关键基因表达分析

目的对银杏进行转录组测序,分析银杏生长发育不同时期类黄酮生物合成关键基因的表达。方法以银杏幼年树和成年树不同时期叶片作为供试材料,利用Illumina HiSeq 2000进行转录组测序,对Unigene进行功能注释,分析银杏类黄酮生物合成关键基因的表达特征。结果转录组测序共获得43 073条Unigene,其中35 179条被注释,基因差异表达(DEG)筛选得到5 117个基因。通过对类黄酮合成相关KEGG途径进行分析,筛选出50条候选基因。分析50条候选基因的表达规律,发现类黄酮合成关键基因均在银杏幼叶中高表达,但在成年树与幼年树同时期叶片之间表达差异不显著。分析与类黄酮合成关系密切的13个基因发现,其中肉桂酸4-羟化酶(C4H)、查耳酮合成酶(CHS)、花青素合成酶(ANS)、花青素还原酶(ANR)和类黄酮O-甲基转移酶(FOMT)基因的表达量较高,类黄酮3′-羟化酶(F3′H)、类黄酮3′,5′-羟化酶(F3′5′H)和黄酮醇合成酶(FLS)基因的表达量则相对较低。结论通过转录组高通量测序,筛选分析了银杏类黄酮生物合成的关键基因及其表达特征,为提高银杏类黄酮产量提供了分子生药学理论基础。     

基于比较转录组的夏枯草组织差异表达分析

目的分析夏枯草Prunella vulgaris果穗、茎、叶片转录组,挖掘夏枯草次生代谢产物生物合成的功能及调节基因。方法利用Illumina高通量测序技术对夏枯草不同组织进行转录组测序,并从差异表达的基因中鉴定次生代谢物生物合成的相关酶基因。结果在夏枯草的3个不同组织的转录本中,共有8 270个Unigenes在至少2个样品间差异显著。对不同组织差异表达的基因进行KEGG富集分析,结果表明,不同组织中苯丙素类生物合成的基因表达均有较大的变化。在夏枯草差异基因中分别搜索三萜类和酚酸类生物合成途径的关键酶,共鉴定到31个三萜类生物合成相关的Unigene,16个酚酸类生物合成相关的Unigenes,113个P450相关的Unigenes。结论本研究为后续发掘夏枯草次生物质代谢合成途径相关功能基因提供依据,也为夏枯草次生代谢物的生物合成调控研究奠定基础。     

基于转录组测序的青天葵差异表达基因分析

目的在青天葵球茎和叶片两个组织的转录组测序数据的基础上,进行差异表达基因(differential expressed gene,DEG)的筛选和分析。方法采用Reads Per Kb per reads(RPKM)方法计算青天葵转录组测序获得的单基因(unigene)在球茎和叶片中表达量,根据球茎和叶片之间RPKM的log2倍数绝对值大于1且错误发现率小于0.01筛选DEG,将DEG归类至Nr、GO、KEGG等公共数据库获取其功能释义;并基于功能注释,挖掘参与青天葵激素合成和信号传导过程的DEG。结果筛选出符合条件的DEG共62 605条,包括了31 488条表达上调DEG和31 117条表达下调DEG。分别有51 652条和14 215条DEG在GO和KEGG数据库中获得注释,涵盖GO数据库的42个功能亚类和KEGG数据库的280条代谢途径。在挖掘到的16个与青天葵激素合成和信号传导过程的DEG中,7个为上调表达基因,9个为下调表达基因。结论通过转录组测序数据较全面地提供了青天葵球茎和叶片基因的差异表达情况,为青天葵的功能基因的克隆和功能分析等提供研究基础和理论依据。     

迷迭香精油诱导肝癌HepG2细胞凋亡后bcl-2和bax基因表达变化的研究

目的：探讨凋亡相关基因bcl-2和bax在迷迭香精油诱导肝癌HepG2细胞凋亡后表达的变化。方法：水蒸气蒸馏法提取迷迭备精油，GC—MS鉴定其成分。采用免疫组化法检测bcl-2蛋白和bax蛋白表达。结果：不同浓度的迷迭香精油处理细胞12、24、48h后bcl-2蛋白的表达降低，bax蛋白的表达升高，与对照组比较均有显著性差异，并且均现时间和剂量依赖性。结论：迷迭香精油诱导肝癌HepG2细胞凋亡与凋亡调控基因bax表达增强，bcl-2表达减少有关。     

猪苓菌核不同发育时期转录组学差异表达研究

药用真菌猪苓菌核的发育代谢过程一直备受关注。为了解猪苓菌核的发育历程,该研究对3个不同时期的猪苓菌核进行了转录组学分析,共得到88. 12 Gb测序数据,包含85 235条unigene。对差异基因进行深度挖掘,筛选了有关跨膜运输、防御、极性生长、形态发育、黑色素合成、细胞壁合成、药效成分麦角甾醇和猪苓多糖合成功能的DEGs,从而推测了菌核发育的分子机制。根据DEGs的通路富集注释结果和相关报道,进一步推测猪苓菌核的发育是由于非生物胁迫(温差和缺氧)和生物胁迫(共生菌蜜环菌和伴生菌的入侵)诱导所致,造成菌核的氧化应激状态,加速合成膜转运蛋白和麦角甾醇从而实现物质的跨膜运输和转换,并通过WD40蛋白实现自身的防御机制。小GTPase和细胞色素P450响应环境信号,进而调控细胞的极性生长和形态发生,期间表皮黑色素沉积,细胞壁加厚,猪苓多糖合成,最终使得猪苓菌核的发育成熟。猪苓菌核不同时期转录组的研究为今后解析其发育历程和相关代谢过程的分子机制奠定了坚实的基础。     

基于代谢组学及转录组学研究灰毡毛忍冬不同品种黄酮类成分差异积累的转录调控机制

为探究湘蕾型和野生型灰毡毛忍冬黄酮类成分差异积累的分子调控机制,该研究以灰毡毛忍冬2个品种花、茎和叶为材料,采用超高效液相色谱-质谱联用(UPLC-MS)技术和高通量转录组测序(RNA-seq)技术进行代谢组学和转录组学联合分析,筛选黄酮类差异代谢物、关键差异表达基因及相应转录因子,并随机挑选其中14条差异表达基因进行qRT-PCR验证。结果表明,代谢组学分析共筛选出黄酮类化合物生物合成相关途径差异代谢物17个,包括柚皮苷查尔酮、芹菜素、槲皮苷等;转录组学分析共筛选出黄酮类化合物生物合成相关途径差异表达基因19条,包括2条CHS、3条CHI等;调控网络分析和加权基因共表达网络分析(WGCNA)筛选出调节黄酮类成分积累的关键酶基因CHS1、FLS1和HCT;同时发现MYB12和LBD4这2个转录因子可能通过调控关键酶基因CHS1、FLS1和HCT的表达来调控黄酮类成分的生物合成;qRT-PCR结果显示,随机选取的14个差异表达基因的表达模式与RNA-seq结果基本一致。该研究初步解析了2个品种灰毡毛忍冬黄酮类成分差异积累的转录调控机制,为进一步阐明其中关键酶基因及相应转录因子对灰毡毛忍冬...     

乌头萜类生物合成代谢的转录组学分析

目的:对乌头的6个器官组织,主根、侧根、茎、叶、花、果实进行转录组测序分析,以研究参与其萜类次级代谢生物合成相关基因的表达谱,挖掘其功能基因。方法:采用Illumina Hi Seq 2500平台测序、组装得unigenes,与公共数据库NR,Swiss-Prot,GO,COG,KOG,KEGG比对、注释以获得差异基因,使用实时荧光定量聚合酶链式反应(Real-time PCR)验证参与萜类次级代谢生物合成的关键基因。结果:本研究共获得156 967 635条Clean Reads(28 254 174 300 bp),经序列合并拼接后获得103 337个unigenes,通过核苷酸和蛋白序列等方面的同源性分析,表明其中37.31%(38 554个unigenes)与其他生物的已知基因具有不同程度的同源性,通过功能分类研究和代谢途径分析(KEGG)获得参与萜类合成158个unigenes(编码5个关键酶)。结论:这些发现的基因将为乌头萜类化合物的生物合成途径及分子机制提供基础数据。     

款冬叶不同发育阶段的转录组学分析

目的比较不同生长时期款冬叶中苯丙素类成分生物合成相关基因的表达,推测苯丙素类成分的最佳合成积累时期,为款冬叶资源的合理利用提供科学依据。方法利用Illumina Hi Seq2500测序技术对不同时期的款冬叶进行转录组测序,获得转录组数据后进行基因功能注释等生物信息学分析,对苯丙素类生物合成相关基因在不同时期的表达进行比较。结果通过转录组测序技术获得46 793个Unigenes,平均长度为952.144 8 bp,其中有4 774个可以被NR、Swiss-Prot、eggNOG、GO、KEGG等公共数据库共同注释。对注释得到的Unigenes进行KEGG代谢通路分析,结果显示款冬叶转录组中有144条Unigenes参与萜类、黄酮类和苯丙素类生物合成,其中萜类65条,苯丙素类64条,黄酮类15条。进一步对参与苯丙素类生物合成相关基因进行动态比较,发现其关键酶PAL、4CL、HCT、CCoAOMT等在9月的表达量最高,即苯丙素类成分在9月的生物合成量可能最高。结论为苯丙素类成分的生物合成途径及调控分析奠定了基础,也为款冬叶资源的开发利用奠定了科学依据。     

基于转录组测序的越南安息香根、茎和叶基因表达分析

目的对越南安息香Styrax tonkinensis进行转录组测序,获得其根、茎和叶的转录组信息特征。方法以越南安息香的根、茎和叶作为研究对象,使用Illumina HiSeq TM 2000进行越南安息香根、茎和叶的转录组测序分析。结果转录组测序根、茎和叶共获得53 835 045条高质量序列(clean reads),Trinity de novo组装获得69 151条Unigenes,平均长度778.51nt。BLAST分析表明分别有41 412(59.89%)、31 189(45.10%)、25 539(36.93%)、16 749(24.22%)个Unigenes在Nr、Swiss-port、KOG、KEGG数据库中得到注释,可归入GO分类的细胞组分、生物过程和分子功能3大类46分支,涉及129条KEGG标准代谢通路,其中有31个次生代谢标准通路。蛋白编码框序列3 461个,涉及高等植物转录因子54个家族。使用MISA软件挖掘10 974个简单重复序列(SSRs),二碱基重复SSRs数量最为丰富,有6282个(57.24%),五碱基重复SSRs最少,占2.45%。结论利用高通...     

基于转录组测序分析干旱胁迫对党参不同组织基因表达的调控

目的 对不同干旱胁迫程度下党参Codonopsispilosula叶、茎、根进行转录组测序，探究干旱胁迫对党参差异基因和多糖合成相关酶基因的调控。方法 利用Illumina HiSeq 2500高通量测序技术对正常水分（85%～90%）和轻度（65%～70%）、中度（50%～55%）、重度干旱胁迫（35%～40%）下生长旺盛期党参叶、茎、根组织进行测序并建立c DNA数据库，经de novo拼接后得到Unigene，并进一步开展生物信息学分析。结果 共获得Unigene 60 377条，分别有51 477、29 896、33 479、35 912、47 378、18 649、31 833条被非冗余数据库（non-redundant protein sequence database,NR）、真核生物蛋白相邻类的聚簇（clusters of orthologous group for eukaryotic complete,KOG）、基因本体（gene ontology,GO）、Swiss-Prot、egg NOG、京都基因与基因组数据库（Kyoto encyclopedia of ge...     

基于代谢组学和转录组学探究草珊瑚叶和根中黄酮类成分差异积累的转录调控机制北大核心CSCD

为了探究草珊瑚叶和根中黄酮差异积累的分子调控机制,该研究以草珊瑚叶和根为材料,采用液相色谱-质谱联用(LC-MS)技术和高通量转录组测序(RNA-seq)技术进行转录组学和代谢组学联合分析,筛选黄酮类差异代谢物和关键差异代谢酶基因,并随机选取其中8个差异表达基因进行qRT-PCR验证。结果表明,通过代谢组学分析发现在草珊瑚叶和根中共获得37个黄酮相关差异代谢物,包括乔松素、根皮苷、柚皮素、山柰酚、无色矢车菊素、5-O-咖啡酰莽草酸等;通过转录组学分析发现草珊瑚叶和根中共获得36条黄酮相关差异基因,包括2条PAL、3条4CL、2条CHS、4条CHI、2条FLS、1条DFR、1条CYP73A、1条CYP75B1、3条PGT1、6条HCT、2条C3′H、1条CCOAOMT、1条ANR、1条LAR、2条3AT、1条BZ1、2条IFTM7、1条CYP81E9。同时,预测了C2H2、bHLH、bZIP等6个转录因子参与调控草珊瑚叶和根中黄酮合成差异积累。qRT-PCR结果显示,随机选取的8个参与黄酮合成的酶基因在草珊瑚叶和根中上下调表达趋势与转录组测序结果一致。该研究初步解析了草珊瑚叶和根中黄酮类成分差异积累的转录调控机制,为进一步阐明其中关键酶基因及相应转录因子对草珊瑚中黄酮类成分积累的调控作用奠定基础。     

基于转录组测序研究不同海拔生境神农香菊的基因表达差异＊

目的 比较神农香菊不同海拔生境的主要生态环境因子差异和基因表达谱差异，研究神农香菊适应不同海拔生境的分子机制。方法 利用WheatA获得海拔相关的主要生态因子如温度和降水等并进行差异统计分析，利用转录组测序获得黄龙堰与宋洛两个不同海拔生境的神农香菊样本的基因表达谱，并基于GO、KEGG等数据库进行基因注释、功能富集分析，获取样本间的差异基因及其关联的代谢通路。结果 光照和温度是影响神农香菊生长发育的最重要的生态因子。不同海拔种植点样本的差异表达基因主要富集到植物激素信号转导、萜类、苯丙素类、不饱和脂肪酸类等次生代谢生物合成途径。其中，参与萜类生物合成途径的相关基因在不同海拔种植样本间显著上调或下调表达。结论 神农香菊可能通过调控萜类关键表达基因如DXS、DXR、HMGCS等的表达适应生境变化。研究结果为深入研究神农香菊不同海拔适应机制及新品种培育提供数据支撑。     

鸡骨草与毛鸡骨草的比较转录组学分析

的比较鸡骨草和毛鸡骨草转录组信息,获得两者比较转录组学的差异。方法利用Illumina HiSeqTM 4000 高通量测序平台对鸡骨草和毛鸡骨草嫩叶样品进行转录组测序;将基因序列与蛋白库作blastx 比 对,通过预测编码蛋白框获得编码蛋白和蛋白序列;使用blastp 和OrthoMCL 软件联用法找出物种间的特有基 因和同源基因,开展特有基因KEGG 和GO 富集分析、单拷贝同源基因分析。结果鸡骨草和毛鸡骨草分别检 测出18 073、17189 个特有基因,同源基因有37881 个。特有基因KEGG 富集分析发现,鸡骨草的核糖体途 径富集基因数量最多且差异最显著,毛鸡骨草中ABC 转运蛋白和胞吞作用途径等显著性富集,植物-病原互作 在两个物种间都呈现显著性富集;GO 富集分析中发现,鸡骨草的胞吞作用调节、囊泡介导的转运调节、核糖 体等显著富集,毛鸡骨草中的微管结合复合物、细胞骨架、细胞膜和蛋白磷酸化作用呈现显著性富集。近缘物 种比对分析结果显示鸡骨草和毛鸡骨草的物种亲缘关系极为相近。此外,鸡骨草和毛鸡骨草分别发现了2 102、 1 794 个简单重复序列（SSR）,两者出现频率最高的串联重复单元都为AG/CT。结论研究揭示了鸡骨草与毛 鸡骨草在转录组学上的异同,两者具有极为接近的亲缘关系,而特有基因差异揭示了两者在生物代谢、信号传 导、发育和功能等方面又存在不同,研究结果将为两者在栽培、生物调控和品质等方面的后续研究奠定基础。     

天麻苯丙烷代谢途径的转录组学分析

目的:对天麻进行转录组学分析,从转录组水平阐释天麻次生代谢产物含量与相关酶活性的关系。方法:采用Illumina Hiseq 4000对新鲜天麻及采后敞放5 d天麻进行转录组测序分析,通过与基因数据库比对分析注释和发现差异表达基因。结果:新鲜天麻及采后敞放5 d天麻分别获得了18 772 600、19 353 329条Clean Reads,Trinity组装后得到63 455条Unigene,通过与NR、Swiss-Prot、GO、COG、KOG、KEGG数据库比对,最终获得18 913个有注释信息的Unigene。天麻采收放置5 d后有953个基因转录表达量上升,1 876个基因转录表达量下降,合计2 829个差异表达基因。634条差异表达基因被COG注释到25个功能类别中,612条差异性表达基因被KEGG注释归属于50条代谢通路。天麻转录组差异表达数据中共11个关键基因被KEGG数据库苯丙烷类合成途径注释,与新鲜天麻相比,放置5 d的天麻CM、PDT、PAL、4CL、CHS、F5H、F3'M基因表达量下降,C4H、CCR表达量上升,CAD表达量既有上升又有下降。结论:本研究获得了较完整的天麻苯丙烷类产物合成代谢通路图,并注释了大量相关基因序列。通过本研究为后续深入的基因功能研究及天麻药材品质形成的分子机理研究、天麻栽培过程的质量控制以及采收后的保鲜处理等奠定了基础。     

基于外周血单个核细胞转录组学的柴归颗粒抗抑郁作用机制研究<

目的 通过转录组测序（RNAsequencing,RNA-seq）技术分析柴归颗粒对慢性不可预知温和应激（chronic unpredictable mild stress,CUMS）模型大鼠外周血单个核细胞（peripheral blood mononuclear cells,PBMCs）的转录组学特征，并采用qRT-PCR法验证其作用机制。方法 将大鼠随机分为对照组、模型组、盐酸文拉法辛（35 mg/kg）组和柴归颗粒低、中、高剂量（4.2、8.3、16.6 g/kg）组，造模同时给药，连续28 d。通过抑郁症传统药效学指标（体质量、旷场实验、糖水偏爱实验及强迫游泳实验）评价柴归颗粒的药效。采用RNA-seq技术对大鼠PBMCs进行转录组学分析，筛选差异基因，对其进行京都基因与基因组百科全书（Kyoto encyclopedia of genes and genomes,KEGG）分析；通过qRT-PCR法对筛选出的关键表达基因进行实验验证。结果 行为学指标显示柴归颗粒能显著改善CUMS模型大鼠的抑郁行为。转录组学中KEGG通路富集分析表明抑郁发生时及柴归颗粒干预后，都显著富集到嘌...     

经典名方五子衍宗丸干预半去势雄性小鼠生精功能的转录组学分析

目的：筛选五子衍宗丸干预半去势雄性小鼠生精功能相关的转录组，探讨其在干预生精功能低下进展中的潜在作用机制。方法：Balb/c小鼠按体质量随机分为假手术组、模型组、丙酸睾丸素组、五子衍宗丸组，每组12只，模型组和给药组摘取右侧睾丸和附睾构建半去势生精功能低下模型，假手术组仅剪开毛皮和右侧阴囊并立即消毒缝合，饲养1周后，五子衍宗丸组每天灌胃五子衍宗丸混悬液1.56 g·kg-1，丙酸睾丸素组肌肉注射丙酸睾丸素0.2 mg·kg-1，每周2次，假手术组和模型组每天灌胃等体积生理盐水，连续给药3周。干预结束后，苏木素-伊红（HE）染色观察睾丸组织病理；酶联免疫吸附测定法（ELISA）检测血清睾酮（T）、促黄体生成素（LH）和卵泡刺激素（FSH）含量；小动物精子自动检测仪测定附睾精子数量和活力；利用转录组学芯片技术检测各组睾丸组织mRNA表达水平，筛选五子衍宗丸调控相关基因转录组，并从中任选3条mRNA进行实时荧光定量聚合酶链式反应（Real-time PCR）验证转录组数据。验证合格的转录组数据通过基因本体（GO）注释分析及京都基因和基因组百科全书（KEGG）信号通路富集分析，解析五子衍宗丸调...