双孢蘑菇不同潮次采收期转录组差异表达分析

Differential expression analysis of the transcriptome at different cycles of harvesting periods of Agaricus bisporus

  • 摘要: 【目的】分析双孢蘑菇不同潮次采收期转录组差异表达,为筛选影响双孢蘑菇产量和品质的关键基因及双孢蘑菇的高效栽培提供理论参考。【方法】以双孢蘑菇A15为试验材料,利用转录组测序(RNA-Seq)技术对双孢蘑菇第1潮次采收期(CZ1)、第2潮次采收期(CZ2)、第3潮次采收期(CZ3)样品进行分析。将RNA-Seq数据序列比对参考基因组筛选差异表达基因,对差异表达基因进行GO功能注释和KEGG信号通路富集分析,并通过实时荧光定量PCR验证RNA-Seq结果的准确性。【结果】不同潮次样品共获得58.54Gb的Cleandata,Q20为98.41%~98.49%,Q30为95.01%~95.24%,有效碱基占比为49.00%~49.70%。以CZ1样品为对照组,在CZ2和CZ3样品中分别筛选到849和1297个差异表达基因,CZ1 vs CZ2组和CZ1 vs CZ3组中共有545个差异表达基因,CZ1 vs CZ2组、CZ1 vs CZ3组和CZ2vsCZ3组中共有90个差异表达基因。差异表达基因GO功能注释分析结果显示,CZ1vsCZ2组、CZ1vsCZ3组中差异表达基因被注释到的GO功能分为生物过程、分子功能和细胞组分三大功能类别。KEGG信号通路分析结果显示,CZ1vsCZ2组和CZ1vsCZ3组中代谢途径、抗生素的生物合成以及次生代谢产物的合成等信号通路被显著富集。从高表达量的差异表达基因中筛选出10个差异表达基因进行分析,CZ3样品中酪氨酸酶基因(TYR)表达量是CZ1样品的7.04倍,CZ1样品中过氧化氢酶基因(CAT3)表达量达到峰值,为CZ3样品的3.20倍。4个差异表达基因的表达趋势与RNA-Seq分析结果基本一致。【结论】从3个不同潮次采收期A15双孢蘑菇样品差异表达基因中筛选到10个高表达量基因,其中包括催化黑色素生成的关键限速酶基因TYR及过氧化氢酶的编码基因CAT3。不同潮次采收期双孢蘑菇的抗氧化能力不同,第1潮次采收期双孢蘑菇的抗氧化能力最强,而第3潮次采收期的抗氧化能力最弱。

     

    Abstract: 【Objective】The transcriptome differential expression of Agaricum bisporus at different cycles of harvesting periods was analyzed to provide theoretical reference for screening key genes affecting yield and quality of A. bisporus and efficient cultivation of A. bisporus. 【Method】Using A. bisporus A15 as experimental material, transcriptome sequencing(RNA-Seq)technology was used to analyze the samples of A. bisporus A15 at the first cycle of harvesting period(CZ1), the second cycle of harvesting period(CZ2)and the third cycle of harvesting period(CZ3). RNA-Seq data sequences were compared with reference genomes to screen for differentially expressed genes. GO functional annotation and KEGG signaling pathway enrichment analysis were performed for differentially expressed genes. The accuracy of RNASeq results was verified by real-time fluorescence quantitative PCR. 【Result】A total of 58.54 Gb Clean data was obtained from samples of different cycles,with Q20 ranging from 98.41% to 98.49%,Q30 ranging from 95.01% to 95.24%,and effective base proportion ranging from 49.00% to 49.70%. Taking CZ1 sample as the control group, 849 and 1297 differentially expressed genes were screened in CZ2 and CZ3 samples respectively. A total of 545 differentially expressed genes were detected in CZ1 vs CZ2 group and CZ1 vs CZ3 group. There were 90 differentially expressed genes in CZ1 vs CZ2 group, CZ1 vs CZ3 group and CZ2 vs CZ3 group. The GO function annotation analysis of differentially expressed genes showed that the GO function of differentially expressed genes in CZ1 vs CZ2 group and CZ1 vs CZ3 group could be divided into 3 categories:biological process, molecular function and cell component. KEGG signaling pathway analysis showed that metabolic pathways, antibiotic biosynthesis and secondary metabolites synthesis were significantly enriched in CZ1 vs CZ2 group and CZ1 vs CZ3 group. Ten differentially expressed genes with high expression were selected from the differentially expressed genes for analysis. The expression level of tyrosinase gene(TYR) in CZ3 sample was as 7.04 times as that of CZ1 sample, and the expression level of catalase gene(CAT3)in CZ1 reached the peak, and was as 3.20 times as that in CZ3. The expression trend of 4 differentially expressed genes was basically consistent with the results of RNA-Seq analysis. 【Conclusion】Ten genes with high expression were selected from the samples of A. bisporus A15 in 3 different cycles of harvesting periods, including TYR, a key rate-limiting enzyme gene that catalyzed melanin production, and CAT3, a catalase encoding gene. The antioxidant capacity of A. bisporus is different at different cycles of harvesting periods. The antioxidant capacity of A. bisporus at the first cycle of harvesting period is the strongest, while the antioxidant capacity at the third cycle of harvesting period is the weakest.

     

/

返回文章
返回