Genome-wide selection of superior reference genes for expression studies in Ganoderma lucidum

Selection and validation of reference genes for normalization of gene expression in Floccularia luteovirens

Floccularia luteovirens is one of the rare edible fungi with high nutritional value found on the Qinghai-Tibet Plateau. However, research at the molecular level on this species is currently constrained due to the lack of reliable reference genes for this species. Thirteen potential reference genes (ACT, GAPDH, EF-Tu, SAMDC, UBI, CLN1, β-TUB, γ-TUB, GTP, H3, UBC, UBC-E2, and GTPBP1) were chosen for the present study, and their expression under various abiotic conditions was investigated. Stability of gene expression was tested using GeNorm, NormFinder, BestKeeper, Delta-Ct, and RefFinder. The results showed that the most suitable reference genes for salt treatment were ACT and EF-Tu. Under drought stress, γ-TUB and UBC-E2 would be suitable for normalization. Under oxidative stress, the reference genes H3 and GAPDH worked well. Under heat stress, the reference genes EF-Tu and γ-TUB were suggested. Under extreme pH stress, UBC-E2 and H3 were appropriate reference genes. Under cadmium stress, the reference genes ACT and UBC-E2 functioned well. In different tissues, H3 and GTPBP1 were appropriate reference genes. The optimal internal reference genes when analyzing all samples were H3 and SAMDC. The expression level of HSP90 was studied to further validate the applicability of the genes identified in this study.

Identification and Validation of Reference Genes for RT-qPCR Analysis in Switchgrass under Heavy Metal Stresses

Switchgrass (Panicum Virgatum L.) has been recognized as the new energy plant, which makes it ideal for the development of phytoremediation on heavy metal contamination in soils with great potential. This study aimed to screen the best internal reference genes for the real-time quantitative PCR (RT-qPCR) in leaves and roots of switchgrass for investigating its response to various heavy metals, such as cadmium (Cd), lead (Pb), mercury (Hg), chromium (Cr), and arsenic (As). The stability of fourteen candidate reference genes was evaluated by BestKeeper, GeNorm, NormFinder, and RefFinder software. Our results identified U2AF as the best reference gene in Cd, Hg, Cr, and As treated leaves as well as in Hg, Pb, As, and Cr stressed root tissues. In Pb treated leaf tissues, 18S rRNA was demonstrated to be the best reference gene. CYP5 was determined to be the optimal reference gene in Cd treated root tissues. The least stable reference gene was identified to be CYP2 in all tested samples except for root tissues stressed by Pb. To further validate the initial screening results, we used the different sets of combinatory internal reference genes to analyze the expression of two metal transport associated genes (PvZIP4 and PvPDB8) in young leaves and roots of switchgrass. Our results demonstrated that the relative expression of the target genes consistently changed during the treatment when CYP5/UBQ1, U2AF/ACT12, eEF1a/U2AF, or 18S rRNA/ACT12 were combined as the internal reference genes. However, the time-dependent change pattern of the target genes was significantly altered when CYP2 was used as the internal reference gene. Therefore, the selection of the internal reference genes appropriate for specific experimental conditions is critical to ensure the accuracy and reliability of RT-qPCR. Our findings established a solid foundation to further study the gene regulatory network of switchgrass in response to heavy metal stress.

CRISPR-Cas9 assisted functional gene editing in the mushroom Ganoderma lucidum

The genetic manipulation of basidiomycete mushrooms is notoriously difficult and immature, and there is a lack of research reports on clustered regularly interspaced short palindromic repeat (CRISPR) based gene editing of functional genes in mushrooms. In this work, Ganoderma lucidum, a famous traditional medicinal basidiomycete mushroom, which produces a type of unique triterpenoid-anti-tumor ganoderic acids (GAs), was used, and a CRISPR/CRISPR-associated protein-9 nuclease (Cas9) editing system for functional genes of GA biosynthesis was constructed in the mushroom. As proof of concept, the effect of different gRNA constructs with endogenous u6 promoter and self-cleaving ribozyme HDV on ura3 disruption efficiency was investigated at first. The established system was applied to edit a cytochrome P450 monooxygenase (CYP450) gene cyp5150l8, which is responsible for a three-step biotransformation of lanosterol at C-26 to ganoderic acid 3-hydroxy-lanosta-8, 24-dien-26 oic acid. As a result, precisely edited cyp5150l8 disruptants were obtained after sequencing confirmation. The fermentation products of the wild type (WT) and cyp5150l8 disruptant were analyzed, and a significant decrease in the titer of four identified GAs was found in the mutant compared to WT. Another CYP gene involved in the biosynthesis of squalene-type triterpenoid 2, 3; 22, 23-squalene dioxide, cyp505d13, was also disrupted using the established CRISPR-Cas9 based gene editing platform of G. lucidum. The work will be helpful to strain molecular breeding and biotechnological applications of G. lucidum and other basidiomycete mushrooms.

Validation of reference genes for quantitative gene expression analysis in Auricularia cornea

Auricularia cornea Ehrenb., previously named A. polytricha (Mont.) Sacc, has become one of the most widely cultivated mushrooms in China. Considerable research has been conducted on its cultivation, pathogen identification, proteomics, and more. However, to the best of our knowledge, no studies have been performed on reference-gene validation in this species. Formerly, reference genes were selected for their expression levels only relied upon from others species, owing to the fact that the gene stability in this species is unknown. In this study, nine candidate genes, including tubulin alpha-1A chain (TUBA1A), β-tubulin (Btu), phosphoglucomutase (Pgm), actin 1 (Act1), protein phosphatase 2A regulatory subunit (PP2A), polyubiquitin (UBQ), glyceraldehyde-3-phosphate dehydrogenase (Gapdh), 18S ribosomal protein (18S) and 28S ribosomal protein (28S), were evaluated among different strains and developmental stages. Four algorithms (i.e., geNorm, NormFinder, BestKeeper and RefFinder) were used to analyze candidate genes. The results revealed that UBQ was the most stable reference gene, while 18S was the least stable. Despite these results, the candidate genes were largely inadequate and only two were considered suitable. Based on candidate gene stability, PP2A and UBQ were identified as a set of usable interior control genes for future analyses in this species. This is the first systematic study conducted for selecting reference genes in A. cornea, and lays the foundation for identifying genes and quantifying gene expression in this species.

Validation of suitable reference genes for quantitative gene expression analysis in Tripterygium wilfordii

Validation of suitable reference genes is critical in quantitative real-time polymerase chain reaction (qRT-PCR) analysis. Suitable and reliable reference genes for the normalization of gene expression data are characterized by high gene expression stability across tissues and different experimental conditions. This study evaluated the gene expression stability of ten reference genes commonly used in Arabidopsis thaliana for their suitability in qRT-PCR analysis in Tripterygium wilfordii Hook.f. The orthologous sequences of these ten candidate genes were identified from T. wilfordii transcriptomic data (Project No. SRX472292). Five algorithms including GeNorm, NormFinder, BestKeeper, ΔCt, and RefFinder were used to assess the gene expression stability of these putative reference genes in different plant tissues and different stress conditions. The results identified ACTINT7 and TBP as the most suitable reference genes across all samples. The gene expressions of TwHMGR (3-hydroxy-3-methylglutaryl coenzyme A reductase, KU246037.1) and of TwDXR (1-deoxy-D-xylulose-5-phosphate reductoisomerase, KJ174341.1) were investigated to validate the suitability of the reference genes. The validation analysis confirmed the suitability of ACTINT7 and TBP as the best reference genes for elucidating secondary metabolite biosynthesis pathway in T. wilfordii. In summary, this study identified the most suitable and reliable reference genes for future qRT-PCR- based studies in T. wilfordii.

Selection and Evaluation of Appropriate Reference Genes for RT-qPCR Normalization of Volvariella volvacea Gene Expression under Different Conditions

Volvariella volvacea (V. volvacea), commonly referred to as Chinese (paddy straw) mushroom, is a basidiomycete with a protein-rich volva and pileus. Selecting appropriate reference genes is a crucial step in the normalization of quantitative real-time PCR data. Therefore, 12 candidate reference genes were selected from the V. volvacea transcriptome based on previous studies and then BestKeeper, geNorm, and NormFinder were used to identify reference genes stably expressed during different developmental stages and conditions. Of the 12 candidate reference genes, SPRY domain protein (SPRYp), alpha-tubulin (TUBα), cyclophilin (CYP), L-asparaginase (L-asp), and MSF1-domain-containing protein (MSF1) were the most stably expressed under different experimental conditions, while 18S ribosomal RNA (18S), 28S ribosomal RNA (28S), and beta-actin (ACTB) were the least stably expressed. This investigation not only revealed potential factors influencing the suitability of reference genes, but also identified optimal reference genes from a pool of candidate genes under a wide range of conditions.

Reference Gene Selection for Quantitative Real-Time PCR of Mycelia from Lentinula edodes under High-Temperature Stress

Housekeeping genes are important for measuring the transcription expression of functional genes; 10 traditional reference genes, TUB, TUA, GADPH, EF1, 18S, GTP, ACT, UBI, UBC, and H2A, were tested for their adequacy in Lentinula edodes (L. edodes). Using specific primers, mRNA levels of these candidate housekeeping genes were evaluated in mycelia of L. edodes, which were treated with high-temperature stress at 37°C for 0, 4, 8, 12, 18, and 24 hours. After treatment, expression stability of candidate genes was evaluated using three statistical software programs: geNorm, NormFinder, and BestKeeper. According to geNorm, TUB had the lowest M values in L. edodes strains 18 and 18N44. Using NormFinder, the best candidate reference gene in strain 18 was TUB (0.030), and the best candidate reference gene in strain 18N44 was UBI (0.047). In BestKeeper analysis, the standard deviation (SD) values of UBC, TUA, H2A, EF1, ACT, 18S, and GTP in strain 18 and those of GADPH and GTP in strain 18N44 were greater than 1; thus, these genes were disqualified as reference genes. Taken together, only UBI and TUB were found to be desirable reference genes by BestKeeper software. Based on the results of three software analyses, TUB was the most stable gene under all conditions and was verified as an appropriate reference gene for quantitative real-time polymerase chain reaction in L. edodes mycelia under high-temperature stress.

Identification of genes related to chlamydospore formation in Clonostachys rosea 67-1

Chlamydospores are specific structures that are of great significance to the commercialization of fungal biopesticides. To explore the genes associated with chlamydospore formation, a biocontrol fungus Clonostachys rosea 67‐1 that is capable of producing resistant spores under particular conditions was investigated by transcriptome sequencing and analysis. A total of 549,661,174 clean reads were obtained, and a series of differentially expressed genes potentially involved in fungal chlamydospore formation were identified. At 36 hr, 67 and 117 genes were up‐ and downregulated in C. rosea during chlamydospore production, compared with the control for conidiation, and 53 and 24 genes were up‐ and downregulated at 72 hr. GO classification suggested that the differentially expressed genes were related to cellular component, biological process, and molecular function categories. A total of 188 metabolism pathways were linked to chlamydospore production by KEGG analysis. Sixteen differentially expressed genes were verified by reverse transcription quantitative PCR, and the expression profiles were consistent with the transcriptome data. To the best of our knowledge, it is the first report on the genes associated with chlamydospore formation in C. rosea. The results provide insight into the molecular mechanisms underlying C. rosea sporulation, which will assist the development of fungal biocontrol agents.

De novo transcriptome sequencing and comprehensive analysis of the heat stress response genes in the basidiomycetes fungus Ganoderma lucidum

Ganoderma lucidum is a valuable basidiomycete with numerous pharmacological compounds, which is widely consumed throughout China. We previously found that the polysaccharide content of Ganoderma lucidum fruiting bodies could be significantly improved by 45.63% with treatment of 42 °C heat stress (HS) for 2 h. To further investigate genes involved in HS response and explore the mechanisms of HS regulating the carbohydrate metabolism in Ganoderma lucidum, high-throughput RNA-Seq was conducted to analyse the difference between control and heat-treated mycelia at transcriptome level. We sequenced six cDNA libraries with three from control group (mycelia cultivated at 28 °C) and three from heat-treated group (mycelia subjected to 42 °C for 2 h). A total of 99,899 transcripts were generated using Trinity method and 59,136 unigenes were annotated by seven public databases. Among them, 2790 genes were identified to be differential expressed genes (DEGs) under HS condition, which included 1991 up-regulated and 799 down-regulated. 176 DEGs were then manually classified into five main responsive-related categories according to their putative functions and possible metabolic pathways. These groups include stress resistance-related factors; protein assembly, transportation and degradation; signal transduction; carbohydrate metabolism and energy provision-related process; other related functions, suggesting that a series of metabolic pathways in Ganoderma lucidum are activated by HS and the response mechanism involves a complex molecular network which needs further study. Remarkably, 48 DEGs were found to regulate carbohydrate metabolism, both in carbohydrate hydrolysis for energy provision and polysaccharide synthesis. In summary, this comprehensive transcriptome analysis will provide enlarged resource for further investigation into the molecular mechanisms of basidiomycete under HS condition.

Functional analysis of an APSES transcription factor (GlSwi6) involved in fungal growth, fruiting body development and ganoderic-acid biosynthesis in Ganoderma lucidum

The APSES transcription factors have been identified as key regulators of fungal development and other biological processes in fungi. In the present study, the function of Ganoderma lucidum GlSwi6, a homolog of Saccharomyces cerevisiae Swi6, was characterized. RNAi was used to examine the function of GlSwi6 in G. lucidum. Silencing GlSwi6 resulted in multiple developmental defects, including reduced fungal growth and increased hyphal branching, and the GlSwi6-silenced strains did not exhibit primordium or fruiting body formation. In addition, the H₂O₂ and ganoderic-acid (GA) levels of the GlSwi6-silenced strains decreased approximately 50% and 25%, respectively, compared with those of the WT strain. Furthermore, the addition of H₂O₂ led to the recovery of the GA levels of GlSwi6-silenced strains, implying that GlSwi6 might regulate GA biosynthesis by regulating the intracellular ROS levels. Taken together, these results indicate that GlSwi6 is involved in fungal growth, development and GA biosynthesis in G. lucidum.