1. Introduction
The introduction should briefly place the study in a broad context and highlight why it is important. It should define the purpose of the work and its significance. The current state of the research field should be carefully reviewed and key publications cited. Please highlight controversial and diverging hypotheses when necessary. Finally, briefly mention the main aim of the work and highlight the principal conclusions. As far as possible, please keep the introduction comprehensible to scientists outside your particular field of research. References should be numbered in order of appearance and indicated by a numeral or numerals in square brackets—e.g., [
1] or [
2,
3], or [
4,
5,
6]. See the end of the document for further details on references.
Agaricus bisporusis one of the largest genera of macrofungi with a pleasant taste and diverse health benefits including of anticancer, antioxidant, anti-diabetic, antimicrobial, and anti-obesity activity [
1,
2]. Due to its high nutritional and medicinal value,
A. bisporus has been considered as one of the most important mushrooms [
3]. Furthermore, numerous studies have highlighted the potential of bioactive compounds derived from
A. bisporus for developing drugs against serious diseases, as demonstrated through bioactivity analyses[
4].
A. bisporus lectin (ABL) and mannose-binding protein (Abmb) have been investigated for their potential in medicinal or pharmaceutical application [
5]. China is the top producer of
A. bisporusand has been for the last two to there decades.
The visual appeal of mushrooms is a key determinant in impacting consumer purchasing decisions [
6]. It is difficult to conserve the quality of harvested
A. bisporus due to the absence of a protective cuticle layer and high respiration rate [
7,
8]. The quality degradation of mushrooms causes moisture and nutrient loss, browning, tissue softening, and cap opening [
9]. In recent years, various postharvest approaches have been implemented to prolong the preservation of
A. bisporus. For example, physical methods such as modified at mosphere packaging, irradiation and pulsed electric fields can extend the shelf life of
A. bisporus [
10,
11,
12]. The senescence of
A. bisporusis delayed by chemical methods such as antimicrobial agents or electrolytic water cleaning, ozone, and coating treatment [
13,
14,
15,
16]. Additionally, the storage temperature is the most crucial factor influencing the respiration, color, and ripening of mushrooms. Increasing the storage temperature accelerates the aging, browning, weight loss, and softening of the mushrooms [
17]. Therefore, for better results, most postharvest
A. bisporus preservation methods should be accompanied by low and constant temperature (above 0 °C) conditioning strategies [
18].
Generally, browning is one of the main causes of quality losses of
A. bisporus [
19]. Browning decreases the commercial value of mushrooms and ultimately leads to considerable economic losses. Substantial evidence indicates that enzymatic browning during storage is the main cause of browning in
A. bisporus during storage [
21]. Enzymatic browning of
A. bisporus is a complex processinvolving various phenolic compounds [
21], phenoloxidases [
22,
23] and multiple other enzymes[
24,
25].The expression of
Abppo genes in the
A. bisporus fruit body (pileus, stipe and gill) browning are inhibited by UV-C treatment [
26]. Moreover, a study integrating metabolomics and transcriptomics reported that the
A. bisporus browning can be affected by the dynamic content of organic acids and trehalose [
27]. Maintenance of cell membrane integrity is also an important factor in inhibiting postharvest browning of
A. bisporus. Lin et al. showed that phase change materials combined with melatonin treatment could regulate the expressions of antioxidant-associated genes and maintain the cell membrane integrity to effectively delay browning of
A. bisporus [
28]. Additionally, melatonin treatment can prominently inhibit electron leakage by significantly increasing the expression levels of
AbNdufB9 and
AbRIP1to protect the cell membrane [
29].
Significant research advances have been achieved in the preservation of
A. bisporus. However, understanding of postharvest mushroom physiology and metabolism remains largely lacking. In recent years, high-throughput techniques have become widely utilized for investigating the molecular-level biological processes and metabolic mechanisms in plants [
30,
31]. In this study, we utilized metabolomic and transcriptomic data to analyze the dynamics of major metabolites and the underlying regulatory pathways in
A. bisporus throughout various storage stages. Our research not only sheds light on the postharvest physiology of
A. bisporus but also contributes to enhancing post-harvest storage practices and the development of effective preservation technologies.
2. Materials and Methods
2.1. Mushroom Materials and Treatments
White mushrooms (A. bisporus, A15) were harvested from Longchen Modern Agriculture Technology Company, in Jiaxing City, Zhejiang Province, China. Mushrooms with consistent shape, size, and excellent quality were carefully selected and promptly transported to the laboratory, where they were pre-cooled at 4°C overnight. The samples were randomly divided into 3 groups, each group had 3 biological replicates, and each replicate contained six A. bisporus fruiting bodies. Subsequently, the mushrooms were put into boxes and stored at 23±1 °C (room temperature, RT), and the mushroom flesh (0.5 cm thick) was collected from three storage periods of 0 h (CK), 6 h (RT_6) and 24 h (RT_24). Samples were ground into powder with liquid nitrogen and stored frozen at 80°C until later use.
2.2. Measurements of Browning Index, Firmness and Weight Loss
The degree of browning, also referred to as the browning index, was accessed by calculating the ratio of the browned area to the cap surface and scored on a scale of one to five [
32]. Firmness was measured by the Stable Micro System (SMS) TAXT Plustexture analyser (Stable Micro System, United Kingdom). For each mushroom, the 2mm diameter cylindrical probe of the firmness mirror was pressed into the mushroom at a speed of 5mm/s, and the hardness was defined as the maximum force in the force versus time curve. Weight loss of whole mushrooms during storage was noted. For each weight measurement, the mass of 20 mushrooms was measured. Weight loss was calculated as reported by Nazir et al. [
33].
2.3. Determinations of Polyphenol Oxidase Activity, Total Phenolic Content, Malondialdehyde (MDA) and Electrolyte Leakage
The extraction and activity of polyphenol oxidase (PPO) were performed as described by Meng et al. [
34] with slight modifications. Two grams of frozen tissue was homogenized in 5 mL of phosphate buffer (0.2 M, pH=6.4), and centrifuged at 12000 rpm for 15 min at 4°C. Then, 400 µL of the supernatant was taken and reacted with 2 mL of phosphate buffer and 900 μL 50mM catechol (Sigma Aldrich, Shanghai, China). Changes in absorbance at 420 nm were recorded over 3 min. A change of 0.01 per minute at 420 nm is defined as one unit (U) of PPO activity. PPO activity was expressed as U per g fresh weight (g
−1 FW).
Total phenolic content was measured accordingto the method described by Liu
et al.[
35]. Half a gram of frozen tissue was homogenized with 5 mL ofice-cold 1% HCl–methanol solution and incubated at 4 °C for 20 min. The mixture was centrifuged at 10000 rpm for 10 min at 4 °C. The collected supernatant was diluted five-fold with 1% HCl–methanol solution. The absorbance was measured at 280 nm, and the phenolic compounds were expressed as OD
280 per gram of fresh weight (g
−1 FW).
Measurement of MDA content was carried out using the method described byPleșoianu and Nour [
36] with some modifications. Two grams of frozen tissue was homogenized with 5 mL of 10% trichloroacetic acid (TCA), and centrifuged at 12000 rpm for 15 minat 4 °C. One milliliter of the supernatant was reacted with 1 mL of 0.67% 2-thiobarbituric acid (TBA) in a water bath at 100 °C for 20 min, and then cooled immediately to room temperature.The mixture was centrifuged at 12000 rpm for 10 min. Absorbencies of the supernatant at 450 nm, 532 nm and 600 nm were measured. The MDA content was calculated using the following formula:
MDA content (μmol g-1)=6.45 × (A532 − A600) − 0.56 × A450.
The determination of electrolyte leakage was conducted according to the previous method proposed by Yan
et al.[
37]. Mushroom pulp pellets (1 cm in diameter and 0.3 cm in thickness) were placed in a triangular flask with 40 mL double distilled water for 10 min at ambient temperature. The initial conductivity (EC1) of the mixture was measured.The flask was boiled for 10 min and then cooled down to ambient temperatures before measuring the final conductivity (EC2). The conductivity of deionized water (EC0) was also measured. Electrolyte leakage was an indicator of the change in membrane permeability and was calculated as follows:
Electrolyte leakage (%) = 100% ×(EC1 − EC0)/(EC2 − EC0)
2.4. Ultrastructural Observation of Cells
The samples were cut into small pieces (1 mm × 1 mm × 3 mm) and immediately prefixed in 2.5 % (v/v) glutaraldehyde solution. After rinsing the stationary phase liquid with phosphate buffered saline (PBS), the samples were then post-fixed with 1% osmium tetroxide for 2 h. Subsequently, the samples were cleaned using PBS again and dehydrated in acetone solutionsat 30%, 50%, 70%, 85%, 90% and then twice at 100% acetone. Samples were then placed in SPI 812 epoxy resin embedding agent overnight. Ultrathin sections (70 nm) were cut using an ultramicrotome (EM UC7, Leica Biosystems, Solms, Germany) and observed under transmission electron microscopy (TEM; Hitachi HT7700, Hitachi Ltd., Tokyo, Japan).
2.5. Transcriptome Sequencing and Data Analysis
RNA-Seq sequencing libraries were generated by NEBNext UltraTM RNA Library Prep Kit for Illumina (New England BioLabs, USA) using high-quality RNA samples, and then sequenced on an Illumina platform to generate150 bp paired-end reads. RNA-Seq was performed using the method described previously[
31].The genome sequence of
A.bisporus (
https://www.ncbi.nlm.nih.gov/datasets/genome/GCF_000300575.1/) was used as the reference. Sequencing readswere filtered by fastpv0.19.3. HISAT2 v2.1.0 was used for mapping clean reads to the reference genome. Gene expression levels were calculated as FPKM(fragments per kilobase million). DESeq2 v1.22.1 was employed to identify the differentially expressed genes (DEGs) with a threshold of the log2 |fold change| ≥ 1 and the false discovery rate (FDR) < 0.05. Three biological replicates were assessed for each sample.
2.6. Widely Targeted Metabolomics Profiling and Analysis
The extraction, identification and quantification of metabolites in A. bisporus flesh were conducted by Wuhan Metware Biotechnology Co., ltd (Wuhan, China). Briefly, metabolite extracts were analyzed using a UPLC-ESI-MS/MS system (UPLC, ExionLC™ AD; MS, Applied Biosystems 4500 Q TRAP). The sample extracts were injected into a UPLC (Agilent SB-C18, 1.8 μm, 2.1 mmx 100 mm), the mobile phase A was pure water with 0.1% formic acid, and phase B was acetonitrile with 0.1 % formic acid. Sample measurements were performed using a gradient program that employed a 95% A, 5% B starting condition. A linear gradient was programmed to 5% A, 95% B over 9 minutes with a composition of 5% A, 95% B held for 1 minute followed by a composition of 95% A, 5.0% B reached in 1.1 minutes and held for 2.9 minutes. The flow rate was set to 0.35 mL/min; the column oven was set to 40 °C; and the injection volume was 4 μL. Alternatively, the effluent was connected to an ESI triple quadrupole linear ion trap (QTRAP)-MS.
The ESI source was operated with the following parameters: source temperature 550°C; ion spray voltage (IS) 5500 V (positive ion mode)/-4500 V (negative ion mode); ion source gas I (GSI), gas II(GSII), curtain gas (CUR) were set at 50, 60, and 25 psi, respectively; the collision-activated dissociation(CAD) was high.QQQ scans were obtained in MRM experiments with the collision gas (nitrogen) setting as the medium. DP (de-clustering potential) and CE (collisional energy) for individual MRM leaps were accomplished by further DP and CE optimization. A specific set of MRM transitions were monitored for each period based on the metabolites eluted during that period. Three biological replicates were performed.
2.7. RNA Extraction and Real-Time Quantitative PCR (qPCR)
Total RNA isolation and purification were conductedwith Trizol (Takara, Biotechnology, Dalian, China), following the manufacturer’s instructions. Reverse transcription was performed with a HiScript II 1st Strand cDNA Synthesis Kit (Vazyme Biotech, Nanjing, China) according to the manufacturer’s instructions.The relative expression level of each gene was determined by real-time quantitative PCR with HiScript II Q RT SuperMix (Vazyme Biotech, Nanjing, China) and the 2-ΔΔCt method. All the reactions were repeated three times, and the AbTubulin gene (AGABI2DRAFT_195658) was used as an internal control for normalization.
2.8. Statistical Analysis
All of the measurements were conducted in triplicate. All statistical analyses were performed by SPSS Statistics 19 (IBM, New York, USA) using one-way analysis of variance (ANOVA). P<0.05 was considered statistically significant for Duncan’s multiple comparison test. The data were fitted and plotted using GraphPad Prism 9.0.0 (GraphPad Software Inc., La Jolla, CA, USA).
4. Discussion
A. bisporus is the most cultivated mushroom worldwide [
38]. Due to high respiration and transpiration rates, it is difficult to conserve the quality of harvested button mushrooms, resulting in the loss of taste and high nutritional value [
39,
40].
A.bisporusis very perishable and has a short shelf life of 1-3 days at room temperature (RT) after harvest [
41]. Figuring out the physiological and molecular mechanisms behind the button mushroomsstorage at RTis quite helpful for breeding practice and postharvest management. The metabolome and transcriptome have been combined to comprehensively analyze the physiological and metabolic information of fruits and vegetables; e.g., jujube [
42], flowering cabbage [
43], radish [
30], and citrus [
44]. In this study, we conducted a comparative analysis of the metabolome and transcriptome of
A. bisporus flesh at various time points during ambient storage. We identified numerous differentially abundant metabolites (DAMs) and differentially expressed genes (DEGs). Our findings highlight that the deterioration of
A. bisporus quality is primarily linked to processes such as membrane lipid metabolism, energy metabolism, and oxidative reactions.
Browning index is the main quality indicator for evaluating the browning of edible mushrooms [
45]. In addition, firmness and weight loss are the main factors determining the storage quality of edible fungi, reflecting the degree of softening andquality loss during storage [
46,
47]. Similar to aprevious study [
48], button mushrooms stored under 23℃ began to slightly brown after 2 h (
Figure 1A,B), with a corresponding decline in firmness (
Figure 1C) and weight (
Figure 1D). Mushroom browning is a complex process that is regulated by a series of enzymes such as PAL and POD. It is widely believed that PPO oxidizes phenolic substrates into anthraquinone compounds, leading to the occurrence of browning [
49]. Here, the content of phenolic substrates showed an increasing trend in ambient storage unlike previous studies [
50,
51]. We speculate that the changing trend of phenolic substances is inconsistent due to different storage temperatures of mushrooms and different measurement times. The antioxidant activity of mushrooms is usually related to the total phenolic content. Ambient temperature may signal stimulation of mushrooms to produce more phenolic compounds as a defense response [
52]. Moreover, in the present study, the trend of PPO activity was consistent with previous studies in that there was a gradual increase during storage [
26,
51]. However, in the transcriptomic data, the transcript abundance of
Abppo1 and
Abppo3 were down-regulated after 6 h storage at RT temperature and did not change significantly after 24 h. Meanwhile, the transcript abundance of
Abppo2,
Abppo4,
Abppo5, and
Abppo6 in mushroom flesh remained unchanged at 6 h and 24 h of storage. These results indicated that PPO-encoding genes might be regulated by multiple factors [
3].
The cell membrane system plays an important role in the normal physiological metabolism of fruits and vegetables [
53]. MDA is the major product of membrane lipid peroxidation and its accumulation is regarded as an indicator ofcell membrane integrity, leading to increased membrane leakage and cellular senescence [
54,
55]. As previously reported [
51], MDA contents and electron leakage of mushroom increased constantly during the storage, accelerating senescence of
A. bisporus. The TEM results further demonstrated that the microstructure integrity of the button mushrooms was disrupted during ambient storage (
Figure 2). Additionally, in ourmetabolome data, content of phospholipids, sphingolipids and fatty acids significantlychanged, with 95% of the phospholipid content increasing significantly after 24 h of ambientstorage (
Figure 4).The phospholipids with elevated levels were lysophosphatidylcholine (LPC) and lysophosphatidylethanolamine (LPE), which are the final metabolic breakdown product of phosphatidylcholine (PC) and phosphatidylethanolamine (PE), the most abundant phospholipid components in cell membranes [
56,
57]. Taken together, the results suggest that the lipid metabolites pathway may be activated in
A. bisporus during ambient storage. Previous studies have provided evidence that postharvest treatments could reduce MDA content and benefit quality due to the maintenance of cell membrane stability in many different vegetables and fruits, such as pineapples, persimmons, and button mushrooms [
58,
59,
60,
61]. The activation of the lipid metabolite pathways in flesh during ambient storage might be a mechanism employed by
A. bisporusto cope with abiotic stress, such as temperature.
During postharvest storage, a large number of metabolites changed rapidly, causing browning in parallel with quality deterioration of the button mushroom. Based on the analysis of the metabolomics, we found a significant increase in the abundances of some major fatty acids, glycerophospholipids, and glycerolipids relative to those in CK. Additionally, we characterized differential regulation of key metabolic pathways, including glycolysis, linoleic acid metabolism, sphingolipid metabolism, and nucleotide metabolism (
Figure S1). Several strategies in microorganisms indicate that ATP from energy metabolic reactions such as glycolysis and linoleic acid metabolism is involved in the synthesis of fatty acids and phospholipids to repair cell membranes [
62,
63,
64]. We further attempted to establish a regulatory network related to lipid metabolism according to the transcriptomic data (
Figure 6). Within this network, transcript abundances of a series of key genes for lipid metabolism, such as PLD, PI3K, PSS, LACS, and IPUT, were significantly induced during storage This effect likely contributed to increasing the content of fatty acids and sphingolipids. Furthermore, there is evidence that in peach, AP2/ERF superfamily transcription factor ABR1plays a critical role in regulating lipid metabolism [
65,
66]. However, the AP2/ERF transcription factor family, is found only in plants [
67].It would be interestingto mine other transcription factors that can regulate lipid metabolism in response to ambient temperature in future studies.
In summary, the metabolic processes of A. bisporus related to postharvest quality during ambient storage were investigated using metabolomic and transcriptomic approaches. Changes in lipid metabolites, particularly the high accumulation of LPC and LPE, indicated that the cell membrane structure of A. bisporus was dynamically remodeled during ambient storage. The changes may have a great impact on the appearance and texture of A. bisporus. Future studies will focus on developing preservation techniques that inhibit cell membrane degradation. Together, our results provide a theoretical basis and reference for further post-harvest preservation management of A. bisporus.
Author Contributions
Conceptualization, G.Y. and N. L.; methodology, M. T. and Y.T.; validation, Y.F.; formal analysis, F.C.; investigation, Y.T. and Y.H.; resources, F.C. and Y.F.; data curation, M. T., Y.T. and F.C.; writing—original draft preparation, M. T., Y.T. and N. L.; writing—review and editing, N. L. and G.Y.; visualization, F.C.; supervision, G.Y. and N. L; project administration, N. L.; funding acquisition, G.Y. and N. L. All authors have read and agreed to the published version of the manuscript.
Figure 1.
Button mushroom change inappearance(A), browning index (B), firmness (C), and weight loss (D) during postharvest storage at 23 °C. Vertical bars represent standard deviation (n = 3). Different lower-case letters represent significant differences between samples under different storage points.
Figure 1.
Button mushroom change inappearance(A), browning index (B), firmness (C), and weight loss (D) during postharvest storage at 23 °C. Vertical bars represent standard deviation (n = 3). Different lower-case letters represent significant differences between samples under different storage points.
Figure 2.
Changes in the cellular ultrastructure of A. bisporus before storage at room temperature (A1, A2, A3) and after 24 h of storage (B1, B2, B3). CW, cell wall; CM, cell membrane; ER, endoplasmic reticulum; M, mitochondria; Magnification power: A, B1. ×2500;A2,B2. ×5000;A3,B3. ×10000. Images are representative of three replicates.
Figure 2.
Changes in the cellular ultrastructure of A. bisporus before storage at room temperature (A1, A2, A3) and after 24 h of storage (B1, B2, B3). CW, cell wall; CM, cell membrane; ER, endoplasmic reticulum; M, mitochondria; Magnification power: A, B1. ×2500;A2,B2. ×5000;A3,B3. ×10000. Images are representative of three replicates.
Figure 3.
Metabonomic analysis of the button mushrooms under different storage times at 23 ℃. (A) PCA and (B) Pearson’s correlation coefficients of different storage times. (C) Pie chart of the types and quantities of metabolites identified. (D) The blue and purple bars indicate the numbers of up- and down-regulated DAMs between each comparison group.
Figure 3.
Metabonomic analysis of the button mushrooms under different storage times at 23 ℃. (A) PCA and (B) Pearson’s correlation coefficients of different storage times. (C) Pie chart of the types and quantities of metabolites identified. (D) The blue and purple bars indicate the numbers of up- and down-regulated DAMs between each comparison group.
Figure 4.
Lipidomic analysis of button mushrooms flesh under different storage times at 23 ℃. (A) The contents of different lipid families and species. Heatmap of contents of phospholipids (B), fatty acids (C), and sphingolipids (D) differentially accumulated in the comparison of CK, RT_6 and RT_24.The red, yellow, and green denote the low, middle, and high content.PC, phosphatidylcholine; LPC, lysophosphatidylcholine; LPE, lysophosphatidylethanolamine.
Figure 4.
Lipidomic analysis of button mushrooms flesh under different storage times at 23 ℃. (A) The contents of different lipid families and species. Heatmap of contents of phospholipids (B), fatty acids (C), and sphingolipids (D) differentially accumulated in the comparison of CK, RT_6 and RT_24.The red, yellow, and green denote the low, middle, and high content.PC, phosphatidylcholine; LPC, lysophosphatidylcholine; LPE, lysophosphatidylethanolamine.
Figure 5.
Transcriptomic analysisofbutton mushrooms flesh under different storage times at 23 ℃. (A) The blue and purple bars indicate the numbers of up- and down-regulated DEGs between each comparison group, respectively. (B) Venn diagram of up- and down-regulated differentially expressed genes, respectively. (C) KEGG annotation and classification of differentially expressed genes.
Figure 5.
Transcriptomic analysisofbutton mushrooms flesh under different storage times at 23 ℃. (A) The blue and purple bars indicate the numbers of up- and down-regulated DEGs between each comparison group, respectively. (B) Venn diagram of up- and down-regulated differentially expressed genes, respectively. (C) KEGG annotation and classification of differentially expressed genes.
Figure 6.
Modulation of lipid metabolism pathway genes during storage at 23 ℃. Solid arrows indicate biosynthetic steps and dashed arrows indicate catabolic steps.The expression patterns are presented by heatmap on the basis of log2 FPKM. The color gradient from green to red corresponds to transcript levels from low to high. Abbreviations used are as follows: MCMT, malonyl-CoA: acyl carrier protein malonyltransferase; LACS, long-chain acyl-CoA synthetase; GPAT, glycerol-3-phosphate acyltransferase; LPA, lysophosphatidic acid; LPAT, lysophosphatidic acid acyltransferase; PLA, phospholipase A; PA, phosphatidic acid; PAP, PA phosphatase;CDS, CDP-diacylglycerol synthase; PSS, base-exchange-type phosphatidylserine synthase; PLD, phospholipase D; PI3K, PI 3-kinase; KSR, ketosphinganine reductase; SK, sphingosine kinase; SPL, sphingosine 1-phosphate lyase; LOH, lag1 longevity assurance homolog; SBH, sphingosine base hydroxylase; CDase, ceramidase; ∆8SLD, delta8 sphingolipid long-chain base desaturase; FAH, fatty acid alpha-hydroxylase; IPUT, inositol phosphoryl ceramide glucuronosyltransferase.
Figure 6.
Modulation of lipid metabolism pathway genes during storage at 23 ℃. Solid arrows indicate biosynthetic steps and dashed arrows indicate catabolic steps.The expression patterns are presented by heatmap on the basis of log2 FPKM. The color gradient from green to red corresponds to transcript levels from low to high. Abbreviations used are as follows: MCMT, malonyl-CoA: acyl carrier protein malonyltransferase; LACS, long-chain acyl-CoA synthetase; GPAT, glycerol-3-phosphate acyltransferase; LPA, lysophosphatidic acid; LPAT, lysophosphatidic acid acyltransferase; PLA, phospholipase A; PA, phosphatidic acid; PAP, PA phosphatase;CDS, CDP-diacylglycerol synthase; PSS, base-exchange-type phosphatidylserine synthase; PLD, phospholipase D; PI3K, PI 3-kinase; KSR, ketosphinganine reductase; SK, sphingosine kinase; SPL, sphingosine 1-phosphate lyase; LOH, lag1 longevity assurance homolog; SBH, sphingosine base hydroxylase; CDase, ceramidase; ∆8SLD, delta8 sphingolipid long-chain base desaturase; FAH, fatty acid alpha-hydroxylase; IPUT, inositol phosphoryl ceramide glucuronosyltransferase.
Figure 7.
Differentially expressed transcription factors (TFs) in button mushrooms. (A) The ratios of differentially expressed transcription factors from different classes. (B) Expression profiles of differentially expressed C2H2 and Zn(II)2Cys6 family members. Up-regulated (red) and down-regulated (green) genes are indicated.
Figure 7.
Differentially expressed transcription factors (TFs) in button mushrooms. (A) The ratios of differentially expressed transcription factors from different classes. (B) Expression profiles of differentially expressed C2H2 and Zn(II)2Cys6 family members. Up-regulated (red) and down-regulated (green) genes are indicated.
Figure 8.
This is a figure. Schemes follow the same formatting. Transcriptome (gene) and metabolome (meta)-combined KEGG enrichment analysis. Joint KEGG enrichment of histograms for RT_6 vs. CK (A) and RT_24 vs. RT_6 (B).
Figure 8.
This is a figure. Schemes follow the same formatting. Transcriptome (gene) and metabolome (meta)-combined KEGG enrichment analysis. Joint KEGG enrichment of histograms for RT_6 vs. CK (A) and RT_24 vs. RT_6 (B).