猪毛菜不同部位转录组特征分析

本研究以猪毛菜的地上部分和地下部分为实验材料,通过转录组测序分析探讨猪毛菜不同组织部位的基因表达特征和代谢物质形成有关的基因调控。测序数据分析结果显示获得192 777条转录本序列。差异表达筛选得到14 424个差异基因,其中上调表达基因6 137个,下调表达基因8 287个;GO功能富集分析得到14 424个差异基因。差异基因与KEGG数据库进行比对,共有5 169个差异基因注释在138个通路上,其中参与次生代谢途径共11条通路呈显著性富集。进一步分析黄酮类化合物合成和异喹啉生物碱生物合成通路中的差异基因和其对应的相关酶。随机选取了黄酮类化合物合成通路上的8个基因,通过RT-qPCR分析验证了转录组数据的准确性。本研究利用高通量测序技术和生物信息分析获得猪毛菜不同部位的转录组信息特征,为猪毛菜的药用成分研究和基因功能鉴定提供了科学依据。     

高粱响应丝黑穗病侵染的转录组分析

高粱是中国重要的粮食作物,用途非常广泛。丝黑穗病在中国的高粱产区均有发生,已成为影响中国高粱生产最重要的病害。本研究采用高通量测序技术对被丝黑穗病病菌侵染的抗感病高粱进行转录组测序,筛选得到差异基因620个,其中229个上调表达,391个下调表达。差异基因中有330个获得GO数据库功能注释,主要涉及生物合成过程、催化活性、细胞过程等诸多生理生化过程;KEGG数据库富集分析发现,共有105个基因被注释到52个代谢通路中,其中类黄酮生物合成、苯丙素生物合成、二苯乙烯、二庚烷和姜辣素的生物合成、植物激素信号传导、苯丙氨酸代谢、植物病原互作、类胡萝卜素生物合成、油菜素内酯生物合成等通路显著富集。有26个差异基因被鉴定为转录因子,分布在10个转录因子家族中。RT-qPCR验证了高粱中差异基因的表达量变化趋势与转录组测序分析结果一致。本研究结果表明,在抗病高粱响应病原菌侵染的过程中,大量基因被诱导或抑制表达,与病原胁迫响应有关的部分转录因子显著上调表达,抗病相关的代谢和信号传导途径被激活。本研究为进一步挖掘高粱抗病基因提供参考。     

白芨块根生长发育转录组分析

本研究以1年生、2年生和3年生的白芨根为研究材料,通过BGISEQ-500高通量测序平台进行转录组测序,测序结果共获得76 552条Unigenes序列。与七大数据库相似性比对结果显示,得到注释的Unigenes共有52 219条,占所有Unigenes的68.12%。对3个不同年限的白芨根转录组数据进行比较分析,共发现37 825条差异基因,涉及了47种生物功能,可划分为分子功能、细胞组分和生物学过程三大类。通过KEGG富集分析发现共涉及134条信号通路,主要的信号通路为代谢通路和次生代谢产物合成通路。共有59条关键酶差异基因参与了甘露糖的合成途径,多条途径参与了葡萄糖的合成;12条根膨大相关的差异基因参与了白芨根的膨大生长调控。RT-qPCR分析验证了与根膨大和糖代谢相关的6个具有代表性的基因表达与转录组结果趋于一致。本研究结果为进一步深入研究白芨根生长发育和活性成分的生物合成的分子调控机制提供理论依据。     

转录组学鉴定高粱营养器官间功能差异的决定基因

本研究采用Illumina HiSeq 4000高通量测序技术,对高粱根、茎和叶的转录组进行测序以及生物信息学分析,以期鉴定出决定营养器官功能异质性的核心基因。结果表明,转录组测序获得165 413 338条原始测序序列(Raw reads),过滤后获得160195618条序列(Clean reads);共鉴定出33204个基因及230685个可变剪接事件;三种营养器官基因表达模式的相似度较低,共筛选出10 083个差异表达基因,其中在不同营养器官均有表达差异的决定基因1320个。对决定基因进行GO富集分析,得到2946个GO功能注释(p0.05),分布于光合作用、结构分子活性和细胞组成相关的GO项中;KEGG富集分析表明差异基因在光合作用、黄酮和黄酮醇生物合成以及叶酸参与的"一碳单位"代谢途径富集程度显著。本研究为深入探讨高粱营养器官间功能差异性机制及器官特异性启动子克隆等工作奠定了数据基础。     

粳糯谷胚乳灌浆期差异表达基因富集分析

小米直链淀粉与支链淀粉的含量是影响蒸煮和食味品质的重要性状。粳谷和糯谷在籽粒灌浆期淀粉合成的途径和成分不同。利用高通量测序技术Illumina Hiseq4000对糯谷‘公谷68’和粳谷‘赤谷4号’(灌浆期第1天和第5天)未成熟籽粒进行转录组测序。结果表明:(1)糯谷和粳谷GBSSⅠ酶活性均呈低-高-低的趋势,二者活性大小存在一定差异。糯性品种‘公谷68’灌浆期第5天与第1天筛选出665个上调差异表达基因,上调基因较下调基因多431个;粳性品种‘赤谷4号’灌浆期第5天与第1天上调基因较下调基因多97个。(2)糯谷A_2-VS-A_1比较组差异基因GO富集在生物过程和分子功能,如种子油体生物发生、17-β-酮类固醇还原酶活性等7个功能;粳谷B_2-VS-B_1差异基因GO富集在光系统Ⅰ中的光捕获、色素结合等8个功能。(3)糯谷A_2-VS-A_1差异基因KEGG主要富集于咖啡因代谢、亚油酸代谢、花青素的生物合成和黄曲霉毒素生物合成途径,粳谷B_2-VS-B_1差异基因KEGG主要富集在光合作用-天线蛋白、亚油酸代谢、咖啡因代谢、油菜素类固醇生物途径,粳糯两个比较组富集过程均包括咖啡因代谢和亚油酸代谢途径。(4)糯谷和粳谷两个比较组中分别筛选到3个(SSⅡ-3, PHO1, AS)和4个(PHO1-1, AS, AGP16, WAX)差异较大且与胚乳粳、糯性相关的基因。以Actin (Si001873)为内参基因,将上述7个差异基因进行RT-qPCR验证,与转录组结果相吻合,说明差异基因与胚乳粳糯性相关。     

基于转录组测序筛选黄牛低氧适应性相关差异基因

为探明藏黄牛对低氧适应的相关差异表达基因并分析差异表达基因功能和代谢通路,以藏黄牛和三江黄牛心脏组织为研究对象,提取总RNA,建立测序文库后利用Illumina HiSeqTM 4000进行测序,经生物信息学分析筛选出差异表达转录本,并对差异表达基因进行GO富集和KEGG通路分析,最后采用荧光定量PCR(RT-qPCR)验证测序数据的准确性。结果表明:藏黄牛相比于三江黄牛得到显著差异表达基因608个,其中上调基因467个,下调基因141个,推测表达量较高的B2M、COX7C、NFATC1、ECE1、CAST基因是与低氧适应性相关的基因;差异转录本GO功能富集分析可知,在分子功能、细胞组分、生物过程分别富集到327,253,2 191个条目,其中富集程度最高的条目分别是细胞发育过程、结合与转录因子活性、细胞前缘。KEGG分析结果表明,差异转录本显著富集在8条通路,筛选低氧适应性候选基因NFATC1是MP-PKG信号通路与T细胞受体信号通路成员,推测其通过cGMP-PKG通路与T细胞受体信号通路参与低氧适应性调控。采用高通量测序技术对藏黄牛和三江黄牛心脏组织进行转录组测序分析,获得大量差异表达基因,为研究动物低氧适应性机制奠定基础。     

基于转录组测序技术挖掘影响布莱凯特黑牛睾丸大小的功能基因

旨在通过对布莱凯特黑牛大小睾丸转录组数据进行分析,挖掘影响睾丸生殖功能的关键候选基因,探究睾丸大小对牛生殖功能维持的分子机制。采集12月龄布莱凯特黑牛睾丸,根据质量、长径和短径将其分为2组(大小睾丸组各3头)进行高通量转录组测序;以|log₂(Fold Change)|≥1且P-value≤0.01为阈值筛选差异基因,对差异基因进行GO和KEGG功能富集分析,筛选出影响睾丸生殖功能的关键基因;对其进行CDS区克隆测序,并运用qRT-PCR、Western Blot、免疫组化等方法进行进一步研究。结果显示,差异表达基因共353个,其中114个基因在小睾丸中表达上调,239个基因表达下调。GO功能注释表明,差异基因参与代谢、发育、生殖等过程;KEGG富集分析显示,差异基因主要富集在Wnt信号传导、PI3K-Akt信号通路等与睾丸生殖功能相关的途径。最终筛选出睾丸生殖功能基因CCN4进行进一步研究。qRT-PCR和Western Blot结果显示,CCN4基因和蛋白在大睾丸中表达量显著低于小睾丸。免疫组化结果显示,CCN4蛋白主要在精细胞外的区域(特别是支持细胞)中表达...     

甘蓝型油菜种子发育灌浆后期的转录组分析

油菜种子发育是产量和品质形成的关键发育阶段,包含了复杂的发育过程和调控网络,有效地解析种子发育的转录调控机制具有重要的意义。以甘蓝型春油菜品种青杂5号为研究材料,利用RNA-seq技术对种子发育的后期(30-DAF,40-DAF)2个发育时间进行转录组测序,筛选差异基因,并利用GO数据库和KEGG数据库注释差异基因功能和可能参与的调控途径。结果表明,从油菜种子灌浆后期的2个时间点的转录组中分别检测到70 850和65 193个表达基因,筛选得到2 654个差异表达基因,其中1 941个基因下调表达,713个基因上调表达,29个基因表达差异倍数|log2Ratio|≥10。GO基因功能分析显示,生物学途径中富集最显著的条目是染色质组装相关的等生物学过程,分子功能方面富集最显著的条目依次是蛋白质代谢、营养库活性等功能类别,而在细胞组件方面富集最显著的条目是染色体相关的等细胞组件。Pathway显著性富集分析显示注释基因最多的途径是次生代谢途径,其次是淀粉、蔗糖代谢途径、苯丙素生物合成途径中的、碳代谢途径和氨基酸生物合成途径。甘蓝型油菜种子发育后期的转录组分析表明,种子发育30-DAF时期次生代谢物、脂质代谢等表达活跃,40-DAF时期逐渐转变为蛋白质、氨基酸生物合成、光合碳代谢、碳代谢等表达活跃,提示油菜种子灌浆后期仍处于复杂的物质与能量代谢调控过程。     

普通小麦颖壳蜡质缺失突变体glossy1的转录组分析

为了明确突变体颖壳蜡质含量显著变化的分子机制,本研究对源自济麦22颖壳蜡质缺失突变体glossy1与野生型进行了转录组分析。结果表明,在glossy1突变体中,共筛选到12,230个差异表达基因,其中5811个基因在突变体中上调表达,6419个下调表达。GO(gene ontology)功能富集分析发现,差异基因主要富集在蜡质合成和转运途径,具体分布在酰基转移酶活性、脂质结合和水解酶活性等条目,由此推测这些途径与小麦穗部蜡质缺失性状是紧密相关的。我们还利用RT-qPCR检测了参与蜡质代谢途径部分基因的表达,结果与转录组结果是一致的。本研究为今后探究小麦蜡质代谢的分子机制和基因调控网络提供了数据支持,同时也为抗逆小麦育种奠定了理论基础。     

杜仲雌花芽2个发育时期转录组分析

杜仲(Eucommia ulmoides)的果实是提取杜仲胶和杜仲籽油的重要原材料,在工业、医疗、食品等领域开发前景广阔。为了揭示雌蕊原基发育相关基因的表达情况,本研究以杜仲果用良种‘华仲6号’(Huazhong No.6)为材料,利用lllumina Hiseq X-10高通量测序平台,分别对花序原基分化期和雌蕊原基分化期的花芽进行转录组测序,通过生物信息学对2个发育时期的雌花芽转录组进行比较分析,筛选出与杜仲雌蕊原基发育相关的差异基因。结果显示,转录组测序共获得68.11 Gb数据,各样品的Clean reads与杜仲基因组进行序列比对,比对效率为91.59%~92.05%,同时2个发育时期共筛选出485个差异表达基因,其中219个差异基因在雌蕊原基分化期表达上调,266个表达下调。差异基因GO富集结果表明,花的发育、光周期现象、激素生物合成以及蛋白结合转录因子活性等相关的生物途径被富集。KEGG富集结果表明淀粉和蔗糖代谢、碳代谢、昼夜节律和植物激素信号转导等代谢通路被富集。结果表明光周期途径是诱导杜仲成花的重要途径,且雌蕊原基发育受碳水化合物、植物激素和其他代谢物质调控。本研究为杜仲花器官发育调控基因的挖掘提供了基础数据,也为果用杜仲的分子育种提供了参考。     

基于转录组测序技术挖掘大豆蛋白质合成相关基因

通过转录组测序技术挖掘蛋白质合成相关基因,为解析蛋白质合成机制奠定基础。以遗传背景相近但蛋白含量差异较大的大豆高蛋白品系冀HJ117(蛋白质含量52. 99%)及其回交亲本冀豆12(蛋白质含量46. 48%)为研究材料,利用转录组测序技术和生物信息学分析,以期挖掘大豆籽粒蛋白质合成相关基因。通过对转录组数据中冀HJ117和冀豆12差异表达基因的分析筛选,共得到336个差异表达基因,其中冀HJ117较冀豆12有195个上调表达基因,141个下调表达基因。通过GO功能富集分析,发现在分子功能类型中注释的差异基因主要与催化和连接等功能相关; KEGG显著富集分析发现差异表达基因主要富集在蛋白质内质网合成途径中,该途径共筛选到34个差异表达基因,其中33个基因在冀HJ117中表达量较高。从该途径中筛选出9个候选基因,采用荧光定量PCR方法检测其表达量,表达趋势与RNA-Seq检测结果基本一致,证实了RNA-Seq数据的可靠性。RNA-Seq测序结果获得了与蛋白质合成相关基因的基本信息,蛋白质内质网合成途径可能是冀HJ117与冀豆12蛋白质含量产生差异的重要通路。     

洋葱转录组测序及黄酮类化合物合成相关基因的分析

洋葱是世界范围内种植的重要蔬菜作物,其鳞茎颜色是洋葱选育和食用的重要性状。本研究对紫皮洋葱和黄皮洋葱进行高通量转录组测序,挖掘黄酮类化合物合成的基因,探讨不同皮色形成的分子基础。本研究共获得5 809个差异表达基因,包括2 991个上调表达基因,2 818个下调表达基因。差异基因GO功能富集分析显示,差异基因主要参与细胞组分、生物过程和分子功能3大类。KEGG显著富集分析显示,差异表达基因主要富集在3条代谢途径上。黄酮类化合物合成相关的通路有类黄酮合成途径、花青素合成途径、异黄酮合成途径、黄酮和黄酮醇合成途径,并筛选出了F3’H、FLS、CHI和DFR等黄酮类合成相关的结构基因。本研究为进一步探讨洋葱皮色相关的代谢合成途径关键酶的功能及其调控机制提供了基础。     

干旱胁迫下枇杷叶片的转录组分析

分析干旱胁迫下枇杷叶片的转录组,挖掘功能基因并对其差异表达基因进行筛选和分析,为枇杷抗旱提供理论依据。利用新一代高通量测序技术测序,对测序结果进行de novo拼接、功能注释和ORF预测,将差异表达基因在COG、GO和KEGG数据库中进行比对注释。测序结果表明,获得转录本共88 530个,平均长度为740.64 bp,ORF41 748条。COG、GO和KEGG数据库将转录本分别划分为24,54个功能类别及291条代谢通路中。25 197个差异表达的基因在30条代谢通路中显著富集,与酶活性、激素合成代谢和信号转导等相关的差异基因积极响应枇杷干旱,其中双萜类、油菜素类固醇和类胡萝卜素3条生物合成途径中的差异基因呈现出较为一致的表达。采用实时荧光定量(qRT-RCR)对选取的差异基因进行验证,其中过氧化物酶、蛋白激酶byr2、丝氨酸苏氨酸蛋白激酶、脱落酸8'-羟化酶、吲哚-3-乙酰乙酸合酶相关基因上调表达,细胞色素P450 734A1基因下调表达。为枇杷提供了较为全面的基因信息和代谢途径数据,为干旱胁迫下枇杷分子调控机制的深入研究奠定了基础。     

不同硬度西瓜果皮的转录组测序及相关基因表达分析

为解析西瓜果皮硬度不同的相关分子机制,挖掘西瓜果皮硬度相关的关键基因提供理论基础,以生育期相似但果皮硬度差异显著的高硬度（901）4-1-1-M和低硬度BSH为试验材料,选取果皮硬度差异最大时期即西瓜授粉30 d后的果皮,利用Illumina HiSeq^TM测序平台进行转录组测序分析,采用实时荧光定量的PCR技术对测序结果进行验证。结果显示,转录组测序共筛选获得1 085个差异表达基因,其中,上调表达基因555个,下调表达基因530个。GO富集分析表明,1 085个差异表达基因显著富集在细胞组分、分子功能和生物过程中的细胞壁、细胞外围、外部封装结构、胞外区、四吡咯骨架、氧化还原酶活性、转移酶活性、果胶酯酶活性等功能类别。KEGG富集分析表明,Cla97C04G075800、Cla97C02G044950、Cla97C09G165820、Cla97C10G195660、Cla97C01G025380等19个差异基因被注释富集程度最高的苯丙烷代谢通路,该路径最终形成木质素、5羟基愈创木酚木质素、愈创木酚木质素和对羟基苯酚木质素代谢物,qRT-PCR和转录组测序数据相关...     

基于转录组挖掘不同性别贵州白山羊肌肉生长发育的关键基因及PI3K/ATK信号通路分析

胴体生长发育是影响山羊养殖效益的重要评价指标,旨在通过RNA-Seq技术挖掘不同性别贵州白山羊生长发育的关键候选基因,为贵州白山羊生长发育研究提供理论依据,以同等饲养水平条件下贵州白山羊为研究对象,测定2周龄不同性别贵州白山羊的屠宰性能指标,通过背最长肌转录组分析,筛选差异表达基因,分析相关基因的信号通路,最后通过RT-qPCR对筛选的基因进行验证。结果表明,测序得到的Raw reads进行过滤后,6个样本共得到78.99 Gb Clean Data,单个样本获得Clean reads在83 030 104～95 739 024条,各样本与参考基因组的比对效率在93.75%～94.79%,共获得25 089个表达基因,筛选差异表达基因共获得1 077个差异基因,其中194个为新发现转录本基因,发现的差异基因中563个在公羊背最长肌中表达显著上调,514个表达显著下调。对获得的1 077个差异基因进行GO功能富集分析,其中587个为显著富集(Q-value≥0.05),主要分为三大类35个亚类。KEGG信号通路分析发现,这些差异基因共参与243条信号通路,其中参与PI3K/AKT信号通路...     

不同生长年限桔梗根部差异表达基因挖掘分析

为深入了解不同生长年限桔梗根部基因表达差异,本研究以1年及3年生紫花桔梗根部为试验材料进行转录组分析。共获取41.84 G clean base,拼接后获得54 425个Unigene,平均长度为1 352 bp。筛选得到3 700个DEGs,其中2 514个在3年生样本中上调表达,1 186个下调表达。3 700个差异表达基因共富集到2 975个GO项,52项显著富集。52项显著富集中,分子功能28项,生物学过程21项,细胞组分3项。KEGG富集表明,上调DEGs显著富集到植物病原体互作、淀粉蔗糖代谢、苯丙烷类生物合成和植物激素信号传导途径,下调DEGs显著富集到光合作用-天线蛋白、光合作用、玉米素生物合成途径。同时对各个途径中显著差异表达基因注释到的差异相关蛋白、转录因子及关键酶基因进行挖掘。荧光定量PCR的表达与转录组测序结果的一致性说明转录组测序数据可靠。该研究通过分析不同生长年限桔梗根部转录水平上的差异,有助于理解不同生长年限桔梗产量、结实率及植株抗病性等方面存在的差异,为后期开展桔梗分子育种、生长发育及相关功能基因解析提供基础信息。     

烟草叶片响应镉胁迫的差异表达基因鉴定及分析

烟草具有超富集镉的能力,严重降低烟叶品质,影响其经济价值。为了阐释烟草响应镉胁迫的分子机制,本研究采集了镉浓度为0和500μmol L^-1培养条件下的烟草叶片进行转录组测序。共获得76.94 Gb有效数据(Clean data),Q30碱基百分比均达到95.43%以上;在镉胁迫的烟草叶片中,共筛选出7735个差异表达基因,其中4833个基因表达上调,2902个基因表达下调,并通过qRT-PCR分析验证了转录组数据的可靠性。对差异转录本进行GO和KEGG富集分析,GO注释表明差异基因涉及代谢过程、应激反应、细胞结构体、催化活性和转录调节活性等;KEGG富集分析表明上调差异基因主要富集在氨基酸的生物合成、碳代谢、氧化磷酸化和柠檬酸循环等通路,下调差异基因则主要富集在光合作用、次生代谢产物的生物合成、代谢途径和植物激素信号转导途径。进一步分析植物激素信号转导通路发现,共有8条植物激素途径以不同的表达方式参与烟草对镉胁迫的响应。激素喷施烟草的实验结果表明,叶片通过调控赤霉素、油菜素内酯和茉莉酸途径以应对镉胁迫;拟南芥激素信号缺失突变体验证实验表明赤霉素、油菜素内酯、茉莉酸...     

火力楠转录组测序与组织差异表达基因分析

火力楠具有重要经济价值与生态功能,但缺乏其基因组信息,限制了其分子生物学、基因功能与分子育种的相关研究。以火力楠根、茎和叶为材料,利用Illumina高通量测序技术进行转录组测序,并对得到的测序数据进行de novo组装,获得97 503条Unigenes,N50为1 010 bp、平均长度852 bp。与公共数据库进行比对,注释到NR、Swiss-Prot数据库的Unigenes分别为60 926、44 249条。将Unigenes与COG数据库比对,有42 344条Unigenes成功注释,根据功能大致分成25类;与GO数据库比对,有15 947条Unigenes获得注释,按功能可分为生物学过程、细胞组分和分子功能3大类55亚类,其中参与的生物学过程较多;以KEGG数据库参考,有26 267条Unigenes参与330条代谢途径分支,以代谢相关的途径较为集中;差异基因表达分析我们得出,根与茎的差异表达基因为2 826个,根与叶之间的差异表达基因为3 230个,这两组差异基因相对校多,除此之外茎与叶差异基因相对较少;差异基因的GO分类、GO富集以及pathway富集分析表明,主要在细胞组分、蛋白质活性、生物合成以及代谢过程等所占比例较高。这些研究结果极大地扩充了火力楠的基因资源,将有助于火力楠基因的发掘与利用、分子标记的开发及其种质资源遗传改良的研究等。     

外源6-BA调控孕穗期低温后小麦幼穗发育的转录组分析

小麦孕穗期对低温非常敏感,低温胁迫后外源喷施6-苄氨基腺嘌呤(6-BA),能够缓解低温胁迫对小麦造成的伤害,通过转录组测序技术分析6-BA提高小麦抗寒性的分子机制。选用低温敏感型品种皖麦52和低温迟钝型品种烟农19为试验材料,在孕穗期低温胁迫后喷施20 mg L–1的6-BA溶液,以喷施等量蒸馏水处理为对照。观察幼穗形态,并测定幼穗可溶性糖含量和淀粉含量。再通过转录组测序筛选并分析差异表达基因,探究差异表达基因的功能及可能参与的调控通路,采用qRT-PCR方法对测序结果进行验证。外源6-BA处理10d后,与对照相比小麦幼穗形态发育良好,可溶性糖、淀粉含量均升高。转录组结果表明,在皖麦52和烟农19中分别鉴定出22,770个和9866个差异基因,其中有661个基因在两个品种中均上调,从中筛选出ARF5、AGPL1、1-SST、SWEET15等基因,这些基因与调节植物激素水平、淀粉合成、糖代谢等有关。对筛选出的差异基因进行GO和KEGG富集分析。GO注释表明皖麦52和烟农19差异基因的功能都主要富集于细胞结构稳定性、代谢、催化活性等。KEGG富集分析表示信号转导、内源激素水平的调节、碳代谢...     

基于转录组测序筛选鸡皮肤毛囊性状相关基因和关键通路

旨在通过对崇仁麻鸡背部皮肤毛囊组织的转录组学分析，筛选关键的差异表达基因和信号通路，探讨鸡皮肤毛囊性状的分子调控机制。采用TRIzol法提取皮肤毛囊组织总RNA,通过Illumina平台进行转录组高通量测序(RNA-Seq)。最后采用实时荧光定量的方法对转录组测序进行验证。研究结果如下：共筛选出了3 856个差异表达基因，在公鸡(BGB vs YGB)比较组合中筛选出了971个差异表达基因，其中274个表达上调，697个表达下调。在母鸡(BMB vs YMB)比较组合中筛选出了3 529个差异表达基因，其中1 477个上调，2 052个下调。在公鸡、母鸡共同的差异表达基因中筛选出了3个与皮肤毛囊性状有关的高表达基因KRT75、KRT6A和KRT14。在GO功能显著富集中，BGB vs YGB比较组合富集到了516个GO term, BMB vs YMB比较组合富集到了1 020个GO term,基因主要集中参与细胞黏附等活动。显著富集了多条KEGG通路，在BGB vs YGB比较组合中富集了8条显著的KEGG通路，在BMB vs YMB比较组合中富集了20条显著的KEGG通路，筛选出了...