芋全长转录组测序分析及淀粉生物合成相关基因挖掘

【目的】获取芋(Colocasia esculenta)全长转录本序列和结构信息,并挖掘淀粉生物合成相关基因和转录本的结构信息,为后续阐明芋淀粉生物合成的分子机制及深入利用芋基因资源提供科学依据。【方法】以荔浦芋1号为试材,采集其3月龄植株的不同组织(叶片、叶柄、球茎、匍匐茎和根)开展混合样三代全长转录组测序,对获得的转录本进行功能注释,通过生物信息学方法分析转录本可变剪接(AS)、可变多聚腺苷酸化(APA)、长链非编码RNA(lncRNA)、转录因子(TF)和简单重复序列(SSR)等的结构信息,通过京都基因和基因组百科全书(KEGG)注释获得芋淀粉生物合成相关的基因和转录本,并对其进行AS和APA分析。【结果】通过对芋叶片、叶柄、球茎、匍匐茎和根混合样开展三代全长转录组测序,共获得全长转录本序列38 043条;通过与参考基因组进行比对,鉴定出新转录本31 058条,其中28 785条获得功能注释;通过对转录本结构分析,共鉴定出9360个AS事件、6436个基因存在poly(A)位点、25 081个SSR位点、304个lncRNA、1911个融合转录本和1608个TF。通过KEGG注释...     

杜鹃花叶片转录组测序数据组装及功能注释

采用2代 Illumina Hi-Seq 测序技术对2年生杜鹃花白凤4号(Rhododendron pulchurum cv. BaiFeng 4)扦插苗叶片的转录组进行测序,获得了213 723 424条 Clean reads数据,经过质控和De novo。组装共获得平均长度为930 nt的53 568个Unigenes。与Nr、Swiss-Prot、KOG、KEGG四大数据库进行BLAST信息比对(E-value≤10^-5),共获得28 877个注释基因。与Nr数据库序列同源性比较发现,杜鹃花与葡萄(Vitis vinifera)具有较高的同源性,而与其他物种的同源性较低。杜鹃花转录组的Unigenes在KOG数据库比对上30 508个注释,分为25类。在GO数据库比对上17 218个Unigenes,其中与抗逆有关的Unigenes有11 928个。在KEGG数据库的132条代谢通路中富集了6 475个杜鹃花Unigenes,其中注释到植物激素生物合成代谢途径和信号转导途径的有384个。同时,有1 062个Unigenes被注释为转录因子,有8 738条SSR标记被挖掘。     

冬瓜高通量转录组测序及分析

【目的】获得冬瓜转录组序列、遗传变异等信息,从中挖掘冬瓜基因数据及SSR分子标记,为冬瓜后续研究提供数据支撑。【方法】以冬瓜嫩叶为材料,利用Illumina HiSeq~(TM)2000技术对冬瓜进行转录组测序,构建数据库从中获得干净序列。经De novo拼接组装后,将获得的单基因簇(Unigene)数据在非冗余蛋白数据库(nonredundant protein database,Nr)、蛋白质序列数据库(Swiss Prot protein database,Swiss Prot)、基因本体论数据库(gene ontology,GO)、蛋白质真核同源数据库(eukaryotic orthologous groups,KOG)、东京基因与基金组百科全书(Kyoto encyclopedia of genes and genomes,KEGG)、蛋白质家族域数据库(protein families database,Pfam)6个公共数据库中进行比对,最终得到冬瓜单基因簇注释信息。利用MISA软件对转录组单基因簇进行搜索,获得单基因簇中的SSR位点。【结果】从冬瓜嫩叶中得到62 021 032条高品质序列,组装后获得40 611条单基因簇,平均长度955 bp。将所有单基因簇在Nr和Swiss Prot数据库中进行比对,结果分别比对到27 474及19 573条单基因簇;在GO数据库中,所注释到的10 659条单基因簇分别匹配到生物功能、分子功能和细胞组分3个本体的47个功能组中;与KOG数据库进行注释比对,根据其功能将注释到的单基因簇划分为25类;KEGG数据库比对注释到10 799条冬瓜的单基因簇,可分为5个大类、19个亚类、125条代谢途径;在Pfam数据库中比对到17 990条单基因簇,分属于369个类群。SSR位点搜索发现,有5 086条单基因簇包含SSR序列,获得5 474个SSR位点。【结论】利用高通量测序获得大量冬瓜转录组信息,有助于从分子水平对冬瓜进行深入研究。     

杜仲雌雄株转录组测序数据组装及基因功能注释

以10年以上树龄的杜仲雌株当年新发枝条上的幼果、嫩芽、叶片和树皮和雄株新发枝条上嫩芽、叶片和树皮为材料,采用Illumina Hi SeqTM2000高通量测序技术进行转录组测序,获得雌株51,574,000条、雄株52,430,502条Clean Reads数据,分别包含总长度为4,641,660,000nt和4,718,745,180nt核苷酸序列数据信息;经拼接组装,获得雌株基因信息长达69,461,730nt的423,339个Contig片段,获得雄株基因信息长达94,814,201nt的542,383个Contig片段;经进一步拼接,分别获得平均长度为288nt的雌株159,434个Unigene片段和平均长度为231nt的雄株257,288个Unigene片段,共有48,761个表达序列标签(EST)。以BLAST(E-value1.0E-5)将Unigene对NR、NT、KEGG和COG数据库进行比对,获得CDS序列35,541条,再通过ESTscan分析获得CDS片段13,220条,共获得48,761条CDS片段。与NR数据库比对发现杜仲雌、雄株转录组Unigene与葡萄相似序列最多(33.8%),其次是蓖麻(11.4%)和杨树(11.2%),与拟南芥的相似序列仅2.3%;根据Unigene与COG数据库比对结果,可将有COG功能的7,571条Unigene分为24类,而根据GO数据库注释,杜仲转录组有GO功能注释的23,314条Unigene可分为生物过程、细胞组分和分子功能3大类55分支。与KEGG数据库比对,杜仲雌、雄株转录组17,468条Ungenes分属128类代谢通路,其中有2,399条属于次生物质代谢途径,314条参与萜类化合物生物合成途径。     

玉木耳转录组测序及褐变相关基因的挖掘

为了从分子水平研究玉木耳褐变的机理,对褐变和未褐变的玉木耳子实体进行了转录组测序分析,共筛选到5915个差异基因,其中有924个Unigene上调表达,有4991个Unigene下调表达。通过基因本体(Gene Ontology, GO)富集分析发现,差异基因显著富集于细胞及代谢过程、单一有机体、细胞和细胞组分、细胞器和大分子复合体、结合作用和催化活性。KEGG代谢通路富集分析表明,富集程度比较显著的通路包括翻译通路、信号转导途径、内分泌系统、折叠排序及退化、运输和分解代谢通路。挑选与组织褐变相关的差异基因进行qRT-PCR表达模式验证,获得的结果与表达谱分析结果一致。     

玉木耳转录组测序及褐变相关基因的挖掘

为了从分子水平研究玉木耳褐变的机理,对褐变和未褐变的玉木耳子实体进行了转录组测序分析,共筛选到5915个差异基因,其中有924个Unigene上调表达,有4991个Unigene下调表达。通过基因本体(Gene Ontology, GO)富集分析发现,差异基因显著富集于细胞及代谢过程、单一有机体、细胞和细胞组分、细胞器和大分子复合体、结合作用和催化活性。KEGG代谢通路富集分析表明,富集程度比较显著的通路包括翻译通路、信号转导途径、内分泌系统、折叠排序及退化、运输和分解代谢通路。挑选与组织褐变相关的差异基因进行qRT-PCR表达模式验证,获得的结果与表达谱分析结果一致。     

藤茶高通量转录组分析及黄酮类化合物合成相关基因挖掘

[目的]分析藤茶高通量转录组序列,从中挖掘出黄酮类化合物合成相关基因,为进一步揭示藤茶黄酮类化合物生物合成调控机制提供理论参考.[方法]分别采集藤茶的幼叶和成熟叶,提取其总RNA构建cDNA文库,采用Il-lumina HiSeqTM 4000高通量测序平台对藤茶叶片进行转录组测序,经过滤处理后运用Trinity组装,将获得的Unigene与Nr、Nt、Pfam、Swiss-Prot、GO、KO和KOG 7个数据库进行比对注释,并预测Unigenes的编码区序列(CDS);基于KEGG信号通路富集分析,发掘藤茶黄酮类化合物合成相关基因.[结果]藤茶叶片转录组测序获得82126236条原始测序序列(Raw reads),过滤处理后得到80156972条高质量序列(Clean reads),进一步组装拼接得到92472条Unige-nes,平均长度为1208 bp,N50长度为1780 bp,其中,至少在1个数据库注释的Unigenes有84217条,占Unigenes总数的91.07%,有8944条Unigenes在7个数据库均被注释,占Unigenes总数的9.67%.在GO数据库成功注释的41116条Unige-nes可分为生物学过程、细胞组分和分子功能三大类,共56个小类;在KOG数据库注释的14553条Unigenes可分成25类,其中,一般功能预测注释成功的Unigenes最多(1946条);其次是翻译后修饰、蛋白质翻转、分子伴侣(1776条),参与次生代谢物质的生物合成、转运和降解的Unigenes较少,仅有319条;KEGG信号通路富集分析发现,共有15262条Unigenes注释到128条KEGG信号通路,以注释为代谢的Unigenes最多,为8694条,其中筛选获得有98个黄酮类化合物合成相关基因,分别编码苯丙烷代谢通路的3种关键酶和类黄酮代谢通路的14种关键酶.藤茶叶片转录组Unigenes与Swiss-Prot和Nr数据库比对,获得52582条CDS序列,ESTScan 3.0.3预测获得35535条CDS序列.[结论]藤茶在细胞过程、代谢过程、单有机体过程、细胞和细胞部分、结合和催化活性能力分布的基因较丰富,在一般功能、翻译、翻译后修饰、蛋白质翻转及分子伴侣的基因表达量较高,具有较强的碳水化合物代谢能力.多种关键酶基因参与藤茶黄酮类化合物的生物合成,推测其生物合成途径存在多条分支,调控机制也较复杂.     

葡萄牙牡蛎外套膜转录组测序及壳色基因挖掘

采用Illumina高通量测序平台对黑壳与金壳葡萄牙牡蛎(Crassostrea angulata)外套膜组织开展转录组测序研究,其中黑壳组总计产出5 857 718 092 nt数据,筛选后的干净阅读子58 004 560个;金壳组产出6 727 447 326 nt数据,干净阅读子66 614 854个。将发掘的新基因与COG、GO、KEGG、Swiss-Prot、NR数据库进行序列比对,最终得到各数据库注释的新基因数量分别为147、397、203、567、1 347个。通过对差异表达基因的序列进行生物学功能注释和聚类分析以及利用已知贝壳基质蛋白序列进行搜索,鉴定出14条可能与壳色表达有关的贝壳基质蛋白基因,其中5条下调基因可能与金壳基因表达相关。     

中国濒危植物金钱松转录组测序 及生物信息学分析

金钱松是中国特有的孑遗单种属裸子植物,现存的自然种群数量很少,多被引种栽培,也是著名的庭院观赏树种。迄今为止,其遗传背景和基因组信息并不清楚,对于金钱松的保护及其遗传结构研究迫切需要基因组资源。采用Illumina HiSeqTM2500高通量测序平台对金钱松叶片进行转录组测序,经de novo组装共获得70 761条Unigene,平均长度为699 bp,N50的长度为1 300 bp,Q20和Q30序列分别占96.59%和91.29%。通过对7个不同的蛋白质和功能域数据库进行比对和功能注释,有43 674条Unigene(61.72%)注释成功。在GO数据库中,有28 355条Unigene按功能被划分成3大类56个小类,以执行生物过程的类区所占比例最多。通过KEGG pathway分析,有14 623条Unigene注释成功,发现了显著性富集的32条代谢通路,以代谢相关的基因最多。在KOG数据库中,有15 931条Unigene被分配到26个基因功能大类中,其中以参与一般功能、转录、翻译、修饰及蛋白运输的基因最为丰富。此外,利用MISA软件对转录组序列进行EST-SSR位点搜索与分析,共检测到2 260条Unigene含有2 462个EST-SSR位点,分布频率为3.48%,其中有180条序列含有一个以上EST-SSR位点,83条序列含有复合EST-SSR位点,以三核苷酸重复基元类型最为丰富,占42.53%(1 047个EST-SSR),重复次数主要以5～8次为主。这些重要的转录组序列为进一步了解金钱松生物学过程的分子机制提供了有价值的信息,并为未来的功能基因组分析、分子标记开发和群体遗传学分析提供了丰富的资源。     

UV-C处理穿心莲转录组分析及穿心莲内酯合成相关基因挖掘

【目的】对紫外线C（UV-C）处理的穿心莲叶片进行转录组测序分析,挖掘穿心莲内酯合成相关基因,为深入探究穿心莲内酯合成相关基因的调控机制和生物学功能提供理论参考。【方法】选用生长60 d的穿心莲幼苗为材料,经UV-C照射2 h后,利用Illumina HiSeq二代测序平台进行转录组测序,并结合生物信息学软件和实时荧光定量PCR （qRT-PCR）对所得序列进行功能注释和相对表达水平验证,挖掘响应UV-C的穿心莲内酯合成相关基因。【结果】共注释到21545个穿心莲基因,有2137个基因呈显著差异表达模式,其中,1147个基因显著上调表达,990个基因显著下调表达。GO功能注释结果显示,UV-C处理造成了氧化胁迫,线粒体氧化磷酸化电子传递链的NADH脱氢酶基因上调表达。KEGG代谢通路富集结果显示,二萜类物质穿心莲内酯的前体（E,E,E）-香叶基香叶基二磷酸（GGPP）合成途径[甲羟戊酸（MVA）途径和2-C-甲基-D-赤藓糖醇-4-磷酸（MEP）途径]差异表达基因显著富集。转录组测序结果和实时荧光定量PCR（qRT-PCR）检测结果均显示,MVA途径的4种合成酶基因HMGS（CXN00016849）、HMGR（CXN00000058）、MVK（CXN00002900）和MVD（CXN00004398）及GGPP合成关键酶基因FPS的表达量显著上调;MEP途径中,仅DXS（CXN00013302）的表达量显著下调。蛋白互作预测结果显示,2个途径中的差异表达基因编码蛋白之间存在互作,HMGS与CYP71家族成员之间,以及FPS1与UGT71、UGT73和UGT76之间的互作网络丰富。【结论】 UV-C照射调控穿心莲重要生命进程,其中,次生代谢、氧化磷酸化、光合作用等途径基因转录本的相对表达量呈UV-C特异性响应模式。合成穿心莲内酯前体GGPP的细胞质MVA途径基因显著上调表达,推测UV-C诱导穿心莲内酯合成的主要场所为细胞质。穿心莲MVA途径基因及CYP71、UGT71、UGT73和UGT76可作为穿心莲内酯合成途径研究的候选分子靶标。     

番石榴果实品质评价及黄酮类化合物合成相关基因挖掘

【目的】综合评价番石榴Psidium guajava不同品种间的果实品质差异并挖掘黄酮类化合物合成的关键基因。【方法】对番石榴6个品种的11个果实品质指标进行测定,并结合主成分分析方法综合评价其品质差异,运用转录组测序技术比较各品种间的差异表达基因(Differentially expressed gene, DEG),通过GO和KEGG富集分析,挖掘黄酮类化合物合成的DEG,利用实时荧光定量PCR (Quantitative real-time PCR,qRT-PCR)研究DEG在不同品种间的特异性表达。【结果】6种番石榴试材中‘金斗香’和‘胭脂红’品质最优,得分较高,‘水晶’和‘西瓜红’较低,‘珍珠’和‘红宝石’居中;‘金斗香’和‘胭脂红’的类黄酮质量分数较高,分别为9.76和10.05 mg/g,是‘水晶’(5.74 mg/g)的1.5倍以上,显著高于其他品种(P>0.05)。转录组测序分析显示,‘金斗香’和‘胭脂红’的DEG聚为一类,其余4种的DEG聚为一类。黄酮类化合物的生物合成途径中CHS、FLS、CYP73A、CYP98A3、DFR、E2.1.1.104、E1.14...     

西瓜果皮蜡粉相关基因定位研究

试验选用果皮无蜡粉西瓜品系‘W1-1’为母本,果皮有蜡粉西瓜品系‘1061’为父本,配制杂交组合获得F2代群体,分析果皮蜡粉性状遗传规律发现,西瓜果皮蜡粉性状受一对显性基因控制。利用群体分离分析(BSA)法筛选得到西瓜果皮有蜡粉基因位于第8号染色体,以双亲材料基因组重测序为基础,在定位区域开发30对引物,其中19对引物具有多态性。利用19对引物标记F2代群体,构建分子遗传图谱,该图谱全长230.29 cM,定位到一个距离为2.78 cM控制西瓜果皮蜡粉的位点,两个与该位点紧密连锁的CAPS标记,分别为w8-6149和w8-0794。利用100份西瓜自然群体材料分析2个与西瓜果皮有蜡粉紧密连锁标记有效性,分子数据与田间数据吻合率分别为79.0%和72.0%。     

An optimized protocol using Steedman's wax for high-sensitivity RNA in situ hybridization in shoot apical meristems and flower buds of cucumber

In situ mRNA hybridization (ISH) is a powerful tool for examining the spatiotemporal expression of genes in shoot apical meristems and flower buds of cucumber. The most common ISH protocol uses paraffin wax; however, embedding tissue in paraffin wax can take a long time and might result in RNA degradation and decreased signals. Here, we developed an optimized protocol to simplify the process and improve RNA sensitivity. We combined embedding tissue in low melting-point Steedman's wax with processing tissue sections in solution, as in the whole-mount ISH method in the optimized protocol. Using the optimized protocol, we examined the expression patterns of the CLAVATA3 (CLV3) and WUSCHEL (WUS) genes in shoot apical meristems and floral meristems of Cucumis sativus (cucumber) and Arabidopsis thaliana (Arabidopsis). The optimized protocol saved 4–5 days of experimental period compared with the standard ISH protocol using paraffin wax. Moreover, the optimized protocol achieved high signal sensitivity. The optimized protocol was successful for both cucumber and Arabidopsis, which indicates it might have general applicability to most plants.     

春大豆花芽cDNA文库的构建及质量分析

【目的】为了解大豆多荚、多粒相关基因的调控机制,构建了大豆花芽cDNA文库,根据要求初步鉴定了所构建文库的质量.【方法】以大豆吉农18突变体的幼嫩花芽为材料提取总RNA,采用SMART技术合成双链cDNA.经蛋白酶K的消化及SfiⅠ酶切后,将所得cDNA克隆到λTriplEx质粒载体中,成功构建了大豆花芽全长cDNA文库.【结果和结论】将初始文库经扩增后保存,检测扩增文库滴度为2.13×108pfu/mL,重组率接近95.3%,菌落PCR鉴定插入片段主要分布在0.5~2.0 kb.插入片段平均大小在1.0 kb左右.表明本研究所构建的文库既满足了目的基因的分离筛选,又可保证全长cDNA文库的获得,该文库的构建为进一步开展相关基因的克隆及分子生物学研究奠定基础.     

金边鲤皮肤黑色素合成及转运相关基因表达分析

[目的]研究影响金边鲤皮肤色素沉着的关键基因及信号通路,揭示其皮肤颜色差异形成的原因,为金边鲤良种培育提供候选基因信息.[方法]基于突变型金边鲤黄色皮肤(Y组)和黑色皮肤(B组)及野生型金边鲤黑色皮肤(W组)的转录组测序数据,筛选分析获得部分共有差异表达基因,利用实时荧光定量PCR检测差异表达基因和黑色素合成相关基因的表达情况,并分析各基因间的相关性.[结果]依据转录组数据筛选获得ACT3、MYH2、TPM1、TPM2和TPM3共5个来源于心肌肾上腺素信号通路及与细胞动力相关的基因;除ACT3基因外,其余4个基因(MYH2、TPM1、TPM2和TPM3)在Y组皮肤中的相对表达量均显著低于B组和W组(P<0.05,下同),与转录组分析的差异基因表达趋势一致.MC-1R、TYR、TYRP1和TYRP2等4个黑色素合成相关基因在金边鲤皮肤中的表达情况为:B组皮肤中MC-1R基因的相对表达量显著高于Y组和W组;TYR、TYRP1和TYRP2基因在3组皮肤中的相对表达量排序均为W组>B组>Y组,其中W组的相对表达量显著高于B组和Y组.相关性分析结果显示,MC-1R和TYR基因与ACT3、MYH2、TPM1、TPM2和TPM3基因呈显著或极显著(P<0.01)正相关,而TYRP1和TYRP2基因与这5个差异表达基因均无显著相关性(P>0.05).[结论]ACT3、MYH2、TPM1、TPM2和TPM3基因及其所在的心肌肾上腺素信号通路可能影响黑色素合成与运动,且与黑色素合成相关基因存在相关性,可作为金边鲤体色变异研究中的标记基因.     

凡纳滨对虾转录组测序分析及肌肉生长发育相关基因的筛选

【目的】筛选出凡纳滨对虾生长发育的功能基因及代谢调控网络,揭示其生长发育的分子机制,为后续开展凡纳滨对虾分子生物学研究提供宝贵的基因数据来源。【方法】以快速生长群体和慢速生长群体的凡纳滨对虾肌肉组织为研究材料,通过Illumina HiSeqTM2500平台对构建的cDNA文库进行高通量测序分析,以StringTie进行拼接组装后利用DESeq2筛选差异表达基因,并基于KOG、GO、nr、COG、Swiss-Pro、KEGG和Pfam等数据库进行差异表达基因功能注释分析。【结果】经拼接组装共获得53458844条Clean reads,各样品Clean reads的Q30均在93.00%以上;Illumina测序获得的Clean reads与参考基因组的比对效率在87.75%~87.80%,说明转录组测序数据真实可靠。采用SnpEff进行SNP/InDel变异注释分析,结果显示,SNP位点中以A>G、G>A、C>T和T>C等4种类型的数量较多（14950~21562个）,InDel位点则以SYNONYMOUS_CODING的数量最多（33630个）。使用StringTie对Mapped reads（比对到参考基因组的Reads）进行拼接,共发掘到4607个新基因,分别输入COG、GO、KEGG、KOG、Pfam、Swiss-Prot、eggNOG和nr数据库中进行序列比对,最终发现共有1098个新基因被注释,以被nr数据库注释的新基因数量最多（1077个）,而被COG数据库注释的新基因数量最少（416个）。基于Q-value<0.05且Fold Change>2的筛选条件,共获得1408个差异表达基因（661个为显著上调表达基因,747个为显著下调表达基因）;1408个差异表达基因被注释到53个GO功能条目中,其中,22条被注释到生物学过程（Biological process）,16条被注释到细胞组分（Cellular component）,15条被注释到分子功能（Molecular function）;KEGG信号通路富集分析发现3条重要的信号通路,分别是溶酶体通路（Lysosome）、氨基糖和核苷酸糖新陈代谢（Amino sugar and nucleotide sugar metabolism）及鞘脂类代谢（Sphingolipid metabolism）。在溶酶体通路中,CTSL基因、Nramp基因、MyoG基因和Myf5基因在凡纳滨对虾快速生长群体肌肉组织中呈上调表达,而Trypsin基因呈下调表达。【结论】通过转录组测序分析从快速生长群体和慢速生长群体的凡纳滨对虾肌肉组织中筛选出1408个差异表达基因（661个为显著上调表达基因,747个为显著下调表达基因）,主要富集在溶酶体、氨基糖和核苷酸糖新陈代谢及鞘脂类代谢等通路上,在对虾肌肉生长发育过程中发挥重要作用。     

牛优势卵泡颗粒细胞与膜细胞中生殖激素合成相关基因的筛选与分析

【目的】探究牛优势卵泡颗粒细胞(GCs)和膜细胞(TCs)中基因的表达差异,筛选并研究与生殖激素合成相关的基因。【方法】在GEO数据库获得GSE83524芯片数据,通过在线软件GEO2R筛选牛优势卵泡中GCs和TCs之间的差异表达基因。使用DAVID软件对获得的差异表达基因进行GO功能富集分析;使用KOBAS软件进行KEGG信号通路分析,筛选出GCs和TCs之间参与生殖激素合成相关信号通路;使用String结合Cytoscape软件构建蛋白与蛋白之间的相互作用网络(PPI)。以牛的优势卵泡为材料,采用Real-time PCR对筛选出的差异表达基因进行验证。【结果】经在线软件GEO2R分析,共获得247个差异表达基因,其中43个表达上调,204个表达下调;247个差异表达基因的GO功能分析结果共分为3大类41组,其中参与生物学过程(BP)的有22组,参与细胞组分(CC)的有12组,参与分子功能(MF)的有7组;KEGG信号通路分析共发现17条通路,其中腺苷酸环化酶3基因(ADCY3)、细胞色素P450 17A1基因(CYP17A1)、腺苷酸环化酶8基因(ADCY8)、鸟嘌呤核苷酸结合蛋白基因(GNAS)、骨形态发生蛋白6基因(BMP6)5个差异表达基因参与了卵巢类固醇生成通路。经PPI网络互作分析,筛选出了前10位的Hub基因;Real-time PCR结果显示,CYP17A1、ADCY8在TCs中的表达量极显著高于GCs (P0.01);ADCY3、BMP6在TCs中的表达量显著高于GCs(P0.05),GNAS在GCs中的表达量极显著高于TCs(P0.01),这5个基因在GCs和TCs中的表达趋势与GEO芯片数据结果一致。【结论】ADCY3、CYP17A1、ADCY8、GNAS、BMP6在GCs与TCs之间表达量存在显著差异,其可能在牛优势卵泡GCs和TCs中对生殖激素的合成及分泌产生了重要作用。