Tracking the best reference genes for RT-qPCR data normalization in filamentous fungi

Act1 out of Action: Identifying Reliable Reference Genes in Trichoderma reesei for Gene Expression Analysis

Trichoderma reesei is a well-established industrial enzyme producer and has been the subject of extensive research for various applications. The basis of many research studies is the analysis of gene expression, specifically with RT-qPCR, which requires stable reference genes for normalization to yield reliable results. Yet the commonly used reference genes, act1 and sar1, were initially chosen based on reports from the literature rather than systematic validation, raising concerns about their stability. Thus, properly evaluated reference genes for T. reesei are lacking. In this study, five potentially new reference genes were identified by analyzing publicly available transcriptome datasets of the T. reesei strains QM6a and Rut-C30. Their expression stability was then evaluated under relevant cultivation conditions using RT-qPCR and analyzed with RefFinder. The two most stable candidate reference genes were further validated by normalizing the expression of the well-characterized gene cbh1 and comparing the results to those obtained using act1 and sar1. Additionally, act1 and sar1 were normalized against the new reference genes to assess the variability in their expression. All five new reference genes exhibited a more stable expression than act1 and sar1. Both in silico and RT-qPCR analysis ranked the so far uncharacterized gene, bzp1, as the most stable. Further, we found that act1 and sar1 have strain- and condition-dependent expression variability, suggesting that they are unsuitable as universal reference genes in T. reesei. Based on these results, we propose to use the combination of bzp1 and tpc1 for the normalization in RT-qPCR analysis instead of act1 and sar1.

Identification of putative fungal reference genes with stable expression from large RNA-seq datasets

RNA-sequencing (RNA-seq) is the dominant technology for genome-wide transcript quantification in various biological studies. The wide applications of RNA-seq have played an essential role in elucidating complex molecular mechanisms of fungal physiology, and have generated large volumes of related data that are valuable for further bioinformatic mining. In this study, we focus on identifying fungal reference genes from large available transcriptome datasets. In total, 44 candidate reference genes from Aspergillus niger were identified through strict statistical analysis of 332 transcriptomic samples. These candidates cover both newly identified genes and previously reported housekeeping genes and were enriched in several basic cellular pathways, such as genes encoding ubiquitin-conjugating enzyme, 26S proteasome regulatory subunits, vacuolar H+-ATPase subunits, mitochondrial import protein and Ras-related GTPase. Moreover, 26 of the newly identified reference genes with a single ortholog in four other fungi showed stable expression patterns across these fungi. Additionally, these new candidates showed more stable expression than the traditionally used reference genes in the tested datasets, such as gapdh, highlighting their potential to improve normalization of RT-qPCR and transcriptome data.

Neuropeptide Y at the crossroads of male reproductive functions in a seasonally breeding reptile, Hemidactylus flaviviridis

Neuropeptide Y (NPY) is known to be a key regulator of reproductive functions across vertebrates. Its role in the modulation of gonadotropin releasing hormone (GnRH) has been extensively explored in mammals and fishes. However, no such report exists in aves, amphibians, or reptiles. Hence, the present study aimed to develop an insight into reproductive phase-dependent expression and role of NPY in the diencephalon region of male wall lizards. Expression of ligand, npy, and its receptor, npyr varied across the annual reproductive cycle of Hemidactylus flaviviridis with the highest expression of both observed during recrudescence. Further, the diencephalon region of recrudescent wall lizards treated with NPY demonstrated an increased mRNA level of gnrh and its receptor, gnrhr. In addition, the current study also elucidates hormonal regulation of diencephalonic npy and npyr wherein neuropeptides like kisspeptin and substance P, the gonadotropin, FSH, as well as sex steroids, dihydrotestosterone (DHT) and 17β-estradiol (E2) inhibited npy and npyr expression in the diencephalon region of wall lizards. With regard to adipokines, leptin stimulated while nesfatin-1 inhibited diencephalonic npy and npyr expression. In conclusion, the current study is the first to present a comprehensive picture of reproductive phase-specific expression pattern, role, and hormonal regulation of neuropeptide Y in the diencephalon region of the male wall lizard, H. flaviviridis.

Screening and stability verification of reference genes in Botrytis cinerea ZX2 fermentation

Botrytis cinerea, an airborne plant pathogen, holds the potential to synthesize sesquiterpenes, which have been used for the industrial production of abscisic acid. Previously, through our genetic technology, we obtained strain ZX2, whose main product 1´,4´-trans-ABA-diol is physiologically active in plants. In this study, 50 L of fed fermentation was carried out with ZX2 strain to study the stability of expression of TUA, TUB, ATC, EF-1, GAPDH, UCE and GTP genes. Four kinds of software (GeNorm, NormFinder, BestKeeper and Delta Ct) were used to analyze the expression stability of candidate genes, and finally the best reference gene was screened by RefFinder. Based on the results, the ACT was the most stable gene. It was used to normalize the expression levels of two genes related to 1´,4´-trans-ABA-diol production (hmgr and bcaba3) when fed-batch fermentation. Guide the selection of appropriate internal reference genes during the fermentation process to accurately quantify the relative transcription levels of target genes in B.cinerea ZX2.

Selection and Validation of Reference Genes for Reverse-Transcription Quantitative PCR Analysis in Sclerotium rolfsii

Reference genes are important for the accuracy of gene expression profiles using reverse-transcription quantitative PCR (RT-qPCR). However, there are no available reference genes reported for Sclerotium rolfsii; it actually has a pretty diverse and wide host range. In this study, seven candidate reference genes (UBC, β-TUB, 28S, 18S, PGK, EF1α and GAPDH) were validated for their expression stability in S. rolfsii under conditions of different developmental stages, populations, fungicide treatments, photoperiods and pHs. Four algorithm programs (geNorm, Normfinder, Bestkeeper and ΔCt) were used to evaluate the gene expression stability, and RefFinder was used to integrate the ranking results of four programs. Two reference genes were recommended by RefFinder for RT-qPCR normalization in S. rolfsii. The suitable reference genes were GAPDH and UBC across developmental stages, PGK and UBC across populations, GAPDH and PGK across fungicide treatments, EF1α and PGK across photoperiods, β-TUB and EF1α across pHs and PGK and GAPDH across all samples. Four target genes (atrB, PacC, WC1 and CAT) were selected for the validation of the suitability of selected reference genes. However, using one or two reference genes in combination to normalize the expression of target genes showed no significant difference in S. rolfsii. In short, this study provided reliable reference genes for studying the expression and function of genes in S. rolfsii.

Identification and evaluation of suitable reference genes for RT-qPCR analyses in Trichoderma atroviride under varying light conditions

BACKGROUND

Trichoderma atroviride is a competitive soil-borne mycoparasitic fungus with extensive applications as a biocontrol agent in plant protection. Despite its importance and application potential, reference genes for RT-qPCR analysis in T. atroviride have not been evaluated. Light exerts profound effects on physiology, such as growth, conidiation, secondary metabolism, and stress response in T. atroviride, as well as in other fungi. In this study, we aimed to address this gap by identifying stable reference genes for RT-qPCR experiments in T. atroviride under different light conditions, thereby enhancing accurate and reliable gene expression analysis in this model mycoparasite. We measured and compared candidate reference genes using commonly applied statistical algorithms.

RESULTS

Under cyclic light-dark cultivation conditions, tbp and rho were identified as the most stably expressed genes, while act1, fis1, btl, and sar1 were found to be the least stable. Similar stability rankings were obtained for cultures grown under complete darkness, with tef1 and vma1 emerging as the most stable genes and act1, rho, fis1, and btl as the least stable genes. Combining the data from both cultivation conditions, gapdh and vma1 were identified as the most stable reference genes, while sar1 and fis1 were the least stable. The selection of different reference genes had a significant impact on the calculation of relative gene expression, as demonstrated by the expression patterns of target genes pks4 and lox1.

CONCLUSION

The data emphasize the importance of validating reference genes for different cultivation conditions in fungi to ensure accurate interpretation of gene expression data.

Lignin-Modifying Enzymes in Scedosporium Species

Scedosporium species are usually soil saprophytes but some members of the genus such as S. apiospermum and S. aurantiacum have been regularly reported as causing human respiratory infections, particularly in patients with cystic fibrosis (CF). Because of their low sensitivity to almost all available antifungal drugs, a better understanding of the pathogenic mechanisms of these fungi is mandatory. Likewise, identification of the origin of the contamination of patients with CF may be helpful to propose prophylactic measures. In this aim, environmental studies were conducted demonstrating that Scedosporium species are abundant in human-made environments and associated with nutrient-rich substrates. Although their natural habitat remains unknown, there is accumulated evidence to consider them as wood-decaying fungi. This study aimed to demonstrate the ability of these fungi to utilize lignocellulose compounds, especially lignin, as a carbon source. First, the lignolytic properties of Scedosporium species were confirmed by cultural methods, and biochemical assays suggested the involvement of peroxidases and oxidases as lignin-modifying enzymes. Scedosporium genomes were then screened using tBLASTn searches. Fifteen candidate genes were identified, including four peroxidase and seven oxidase genes, and some of them were shown, by real-time PCR experiments, to be overexpressed in lignin-containing medium, thus confirming their involvement in lignin degradation.

Adaptation to glucose starvation is associated with molecular reorganization of the circadian clock in Neurospora crassa

The circadian clock governs rhythmic cellular functions by driving the expression of a substantial fraction of the genome and thereby significantly contributes to the adaptation to changing environmental conditions. Using the circadian model organism Neurospora crassa, we show that molecular timekeeping is robust even under severe limitation of carbon sources, however, stoichiometry, phosphorylation and subcellular distribution of the key clock components display drastic alterations. Protein kinase A, protein phosphatase 2 A and glycogen synthase kinase are involved in the molecular reorganization of the clock. RNA-seq analysis reveals that the transcriptomic response of metabolism to starvation is highly dependent on the positive clock component WC-1. Moreover, our molecular and phenotypic data indicate that a functional clock facilitates recovery from starvation. We suggest that the molecular clock is a flexible network that allows the organism to maintain rhythmic physiology and preserve fitness even under long-term nutritional stress.

Characterization of Pseudogymnoascus destructans conidial adherence to extracellular matrix: Association with fungal secreted proteases and identification of candidate extracellular matrix binding proteins

Pseudogymnoascus destructans is the etiological agent of white-nose syndrome (WNS), a fungal skin infection of hibernating bats. Pathophysiology of the disease involves disruption of bat metabolism and hibernation patterns, which subsequently causes premature emergence and mortality. However, information on the mechanism(s) and virulence factors of P. destructans infection is minimally known. Typically, fungal adherence to host cells and extracellular matrix (ECM) is the critical first step of the infection. It allows pathogenic fungi to establish colonization and provides an entry for invasion in host tissues. In this study, we characterized P. destructans conidial adherence to laminin and fibronectin. We found that P. destructans conidia adhered to laminin and fibronectin in a dose-dependent, time-dependent and saturable manner. We also observed changes in the gene expression of secreted proteases, in response to ECM exposure. However, the interaction between fungal conidia and ECM was not specific, nor was it facilitated by enzymatic activity of secreted proteases. We therefore further investigated other P. destructans proteins that recognized ECM and found glyceraldehyde-3-phosphate dehydrogenase and elongation factor 1-alpha among the candidate proteins. Our results demonstrate that P. destructans may use conidial surface proteins to recognize laminin and fibronectin and facilitate conidial adhesion to ECM. In addition, other non-specific interactions may contribute to the conidial adherence to ECM. However, the ECM binding protein candidates identified in this study highlight additional potential fungal virulence factors worth investigating in the P. destructans mechanism of infection in future studies.

Identification of reference genes and their validation for gene expression analysis in phytopathogenic fungus Macrophomina phaseolina

Macrophomina phaseolina is a soil-borne pathogenic fungus that infects a wide range of crop species and causes severe yield losses. Although the genome of the fungus has been sequenced, the molecular basis of its virulence has not been determined. Identification of up-regulated genes during fungal infection is important to understand the mechanism involved in its virulence. To ensure reliable quantification, expression of target genes needs to be normalized on the basis of certain reference genes. However, in the case of M. phaseolina, reference genes or their expression analysis have not been reported in the literature. Therefore, the objective of this study was to evaluate 12 candidate reference genes for the expression analysis of M. phaseolina genes by applying three different fungal growth conditions: a) during root and stem infection of soybean, b) in culture media with and without soybean leaf infusion and c) by inoculating a cut-stem. Based on BestKeeper, geNorm and NormFinder algorithms, CYP1 was identified as the best recommended reference gene followed by EF1β for expression analysis of fungal gene during soybean root infection. Besides Mp08158, CYP1 gene was found suitable when M. phaseolina was grown in potato-dextrose broth with leaf infusion. In the case of cut-stem inoculation, Mp08158 and Mp11185 genes were found to be most stable. To validate the selected reference genes, expression analysis of two cutinase genes was performed. In general, the expression patterns were similar when the target genes were normalized against most or least stable gene. However, in some cases different expression pattern can be obtained when least stable gene is used for normalization. We believe that the reference genes identified and validated in this study will be useful for gene expression analysis during host infection with M. phaseolina.

Cytokinin Regulates Energy Utilization in Botrytis cinerea

The plant hormone cytokinin (CK) is an important developmental regulator. Previous work has demonstrated that CKs mediate plant immunity and disease resistance. Some phytopathogens have been reported to secrete CKs and may manipulate CK signaling to improve pathogenesis. In recent work, we demonstrated that CK directly inhibits the development and virulence of fungal phytopathogens by attenuating the cell cycle and reducing cytoskeleton organization. Here, focusing on Botrytis cinerea, we report that CK possesses a dual role in fungal biology, with role prioritization being based on sugar availability. In a sugar-rich environment, CK strongly inhibited B. cinerea growth and deregulated cytoskeleton organization. This effect diminished as sugar availability decreased. In its second role, we show using biochemical assays and transgenic redox-sensitive fungal lines that CK can promote glycolysis and energy consumption in B. cinerea, both in vitro and in planta. Glycolysis and increased oxidation mediated by CK were stronger in low sugar availability, indicating that sugar availability could indeed be one possible element determining the role of CK in the fungus. Transcriptomic data further support our findings, demonstrating significant upregulation to glycolysis, oxidative phosphorylation, and sucrose metabolism upon CK treatment. Thus, the effect of CK in fungal biology likely depends on energy status. In addition to the plant producing CK during its interaction with the pathogen for defense priming and pathogen inhibition, the pathogen may take advantage of this increased CK to boost its metabolism and energy production, in preparation for the necrotrophic phase of the infection.

IMPORTANCE The hormone cytokinin (CK) is a plant developmental regulator. Previous research has highlighted the involvement of CK in plant defense. Here, we report that CK has a dual role in plant-fungus interactions, inhibiting fungal growth while positively regulating B. cinerea energy utilization, causing an increase in glucose utilization and energy consumption. The effect of CK on B. cinerea was dependent on sugar availability, with CK primarily causing increases in glycolysis when sugar availability was low, and growth inhibition in a high-sugar environment. We propose that CK acts as a signal to the fungus that plant tissue is present, causing it to activate energy metabolism pathways to take advantage of the available food source, while at the same time, CK is employed by the plant to inhibit the attacking pathogen.

Antifungal effect of triclosan on Aspergillus fumigatus: quorum quenching role as a single agent and synergy with liposomal amphotericin-B

The purpose of this research was to determine Aspergillus fumigatus conidial viability and its biofilm formation upon treatment with triclosan and amphotericin-B loaded liposomes. A. fumigatus was treated with the antimicrobials, triclosan and liposomal amphotericin-B (L-AMB), in single and combined supplementation. To quantify the cells’ viability upon treatments, resazurin-based viability assay was performed. Confocal laser scanning microscopy was done by applying FUN-1 stain to screen the role of the agents on extracellular polymeric substances. Total A. fumigatus biomass upon treatments was estimated by using crystal violet-based assay. To study the agents’ effect on the conidial viability, flow cytometry analysis was performed. Expression levels of A. fumigatus genes encoding cell wall proteins, α-(1,3)-glucans and galactosaminogalactan were analysed by real-time polymerase chain reaction assay. A synergistic interaction occurred between triclosan and L-AMB when they were added sequentially (triclosan + L-AMB) at their sub-minimum inhibitory concentrations, the triclosan and L-AMB MICs were dropped to 0.6 and 0.2 mg/L, respectively, from 2 to 1 mg/L. Besides, L-AMB and triclosan contributed to the down-regulation of α-(1,3)-glucan and galactosaminogalactan in A. fumigatus conidia and resulted in less conidia aggregation and mycelia adhesion to the biotic/abiotic surfaces; A. fumigatus conidia-became hydrophilic upon treatment, as a result of rodlet layer being masked by a hydrophilic layer or modified by the ionic strength of the rodlet layer. In A. fumigatus, the potential mechanisms of action for L-AMB might be through killing the cells and for triclosan through interrupting the cells’ development as a consequence of quorum quenching.

ThhspA1 is involved in lacA transcriptional regulation of Trametes hirsuta AH28-2 exposed to o-toluidine

White rot fungi, especially Trametes spp., respond to a wide range of aromatic compounds and dramatically enhance laccase activity, while the activation mechanisms remain to be elucidated. Here, we show that an Hsp70 homolog named ThhspA1 regulates the transcription of laccase LacA in Trametes hirsuta AH28-2 when confronted with o-toluidine. ThhspA1 is pulled down by lacA promoter sequence from the nuclear mixture extracted from T. hirsuta AH28-2 induced by 2 mM o-toluidine. Silencing of ThhspA1 results in a sharp decrease in lacA transcripts and laccase activity in vivo. By comparison, ThhspA1 overexpression does not affect lacA transcription, and laccase activity shows slight enhancement or remains unchanged upon induction with o-toluidine. Electrophoretic mobility shift assays suggest a direct interaction between ThhspA1 and the lacA promoter region. Further investigation shows that the integrity of ThhspA1 is critical since its substrate binding domain (SBD) and nucleotide-binding domain (NBD) are both necessary for DNA binding, with a higher affinity of SBD than NBD based on fluorescence polarization assay. Our results demonstrate that ThhspA1 functions as an aromatic-stress-related DNA binding transcriptional factor required for LacA expression.

Cytokinin Inhibits Fungal Development and Virulence by Targeting the Cytoskeleton and Cellular Trafficking

Cytokinin (CK) is an important plant developmental regulator, having activities in many aspects of plant life and response to the environment. CKs are involved in diverse processes in the plant, including stem cell maintenance, vascular differentiation, growth and branching of roots and shoots, leaf senescence, nutrient balance, and stress tolerance. In some cases, phytopathogens secrete CKs. It has been suggested that to achieve pathogenesis in the host, CK-secreting biotrophs manipulate CK signaling to regulate the host cell cycle and nutrient allocation. CK is known to induce host plant resistance to several classes of phytopathogens from a few works, with induced host immunity via salicylic acid signaling suggested to be the prevalent mechanism for this host resistance. Here, we show that CK directly inhibits the growth, development, and virulence of fungal phytopathogens. Focusing on Botrytis cinerea (Bc), we demonstrate that various aspects of fungal development can be reversibly inhibited by CK. We also found that CK affects both budding and fission yeast in a similar manner. Investigating the mechanism by which CK influences fungal development, we conducted RNA next-generation sequencing (RNA-NGS) on mock- and CK-treated B. cinerea samples, finding that CK alters the cell cycle, cytoskeleton, and endocytosis. Cell biology experiments demonstrated that CK affects cytoskeleton components and cellular trafficking in Bc, lowering endocytic rates and endomembrane compartment sizes, likely leading to reduced growth rates and arrested developmental programs. Mutant analyses in yeast confirmed that the endocytic pathway is altered by CK. Our work uncovers a remarkably conserved role for a plant growth hormone in fungal biology, suggesting that pathogen-host interactions resulted in fascinating molecular adaptations on fundamental processes in eukaryotic biology.

Coprinopsis cinerea uses laccase Lcc9 as a defense strategy to eliminate oxidative stress during fungal-fungal interactions

Frequently, laccases are triggered during fungal cocultivation for overexpression. The function of these activated laccases during coculture has not been clarified. Previously, we reported that Gongronella sp. w5 (w5) (Mucoromycota, Mucoromycetes) specifically triggered the laccase Lcc9 overexpression in Coprinopsis cinerea (Basidiomycota, Agaricomycetes). To systematically analyze the function of the overexpressed laccase during fungal interaction, C. cinerea mycelia before and after the initial Lcc9 overexpression were chosen for transcriptome analysis. Results showed that accompanied by specific utilization of fructose as carbohydrate substrate, oxidative stress derived from antagonistic compounds secreted by w5 appears to be a signal critical for laccase production in C. cinerea. Reactive oxygen species (ROS) decrease in the C. cinerea wild-type strain followed the increase in laccase production and then, lcc9 transcription and laccase activity stopped. By comparison, increased H₂O₂ content and mycelial ROS levels were observed during the entire cocultivation in lcc9 silenced C. cinerea strains. Moreover, lcc9 silencing slowed down the C. cinerea mycelial growth, affected hyphal morphology, and decreased the asexual sporulation in coculture. Our results showed that intracellular ROS acted as signal molecules to stimulate defense responses by C. cinerea with the expression of oxidative stress response regulator Skn7 and various detoxification proteins. Lcc9 takes part as a defense strategy to eliminate oxidative stress during the interspecific interaction with w5. Importance: The overproduction of laccase during interspecific fungal interactions is notoriously known. However, the exact role of the up-regulated laccases remains underexplored. Based on comparative transcriptomic analysis of C. cinerea and gene silencing of laccase Lcc9, here we show that oxidative stress derived from antagonistic compounds secreted by Gongronella sp. w5 was a signal critical for laccase Lcc9 production in Coprinopsis cinerea. Intracellular ROS acted as signal molecules to stimulate defense responses by C. cinerea with the expression of oxidative stress response regulator Skn7 and various detoxification proteins. Ultimately, Lcc9 takes part as a defense strategy to eliminate oxidative stress and help cell growth and development during the interspecific interaction with Gongronella sp. w5. These findings deepened our understanding of fungal interactions in their natural population and communities.

Selection of suitable reference genes for quantitative real time PCR in different Tulasnella isolates and orchid-fungus symbiotic germination system

Under natural conditions, mycorrhizal symbiosis accompanies nearly the entire life cycle of orchids from seed germination through to flowering and fruiting. Tulasnella-like orchid mycorrhizal fungi are the most common mycorrhizal fungi found in association with orchid species. Presently suitable reference genes have not been systematically selected for the quantification of gene expression via Real-Time Quantitative Reverse Transcription PCR (RT-qPCR). We evaluated 12 candidate Tulasnella genes in nine different Tulasnella isolates and in the Dendrobium-fungal symbiotic germination associations followed by statistical analysis using the programs Bestkeeper, geNorm, and Normfinder to analyze the expression stability of the individual genes. The results showed that the EF2, UBC, and PP2A genes had the highest rankings with relatively stable expression levels across the different genotypes and during the symbiotic seed germination process by the three programs, and may be suitable for RT-qPCR normalization. Furthermore, the gene encoding C-5 Sterol desaturase (C5SD) was selected to verify the reliability of EF2, UBC, and PP2A expression during the Tulasnella-Dendrobium symbiotic seed germination process. This study is the first systematic exploration of optimal reference genes for gene expression studies during the colonization of orchid seeds by the mycorrhizal fungus Tulasnella.

Reference genes for gene expression analysis in the fungal pathogen Neonectria ditissima and their use demonstrating expression up-regulation of candidate virulence genes

European canker, caused by the necrotrophic fungal phytopathogen Neonectria ditissima, is one of the most damaging apple diseases worldwide. An understanding of the molecular basis of N. ditissima virulence is currently lacking. Identification of genes with an up-regulation of expression during infection, which are therefore probably involved in virulence, is a first step towards this understanding. Reverse transcription quantitative real-time PCR (RT-qPCR) can be used to identify these candidate virulence genes, but relies on the use of reference genes for relative gene expression data normalisation. However, no report that addresses selecting appropriate fungal reference genes for use in the N. ditissima-apple pathosystem has been published to date. In this study, eight N. ditissima genes were selected as candidate RT-qPCR reference genes for gene expression analysis. A subset of the primers (six) designed to amplify regions from these genes were specific for N. ditissima, failing to amplify PCR products with template from other fungal pathogens present in the apple orchard. The efficiency of amplification of these six primer sets was satisfactory, ranging from 81.8 to 107.53%. Analysis of expression stability when a highly pathogenic N. ditissima isolate was cultured under 10 regimes, using the statistical algorithms geNorm, NormFinder and BestKeeper, indicated that actin and myo-inositol-1-phosphate synthase (mips), or their combination, could be utilised as the most suitable reference genes for normalisation of N. ditissima gene expression. As a test case, these reference genes were used to study expression of three candidate virulence genes during a time course of infection. All three, which shared traits with fungal effector genes, had up-regulated expression in planta compared to in vitro with expression peaking between five and six weeks post inoculation (wpi). Thus, these three genes may well be involved in N. ditissima pathogenicity and are priority candidates for further functional characterization.

Enhancing Activities of Salt-Tolerant Proteases Secreted by Aspergillus oryzae Using Atmospheric and Room-Temperature Plasma Mutagenesis

Aspergillus oryzae 3.042 was mutagenized using atmospheric and room-temperature plasma (ARTP) technology to enhance its salt-tolerant proteases activity. Compared to the starting strain, mutant H8 subjected to 180 s of ARTP treatment exhibited excellent genetic stability (15 generations), growth rate, and significantly increased activities of neutral proteases, alkaline proteases, and aspartyl aminopeptidase during fermentation. Mutant H8 significantly enhanced the contents of 1-5 kDa peptides, aspartic acid, serine, threonine, and cysteine in soy sauce by 16.61, 7.69, 17.30, 8.61, and 45.00%, respectively, but it had no effects on the contents of the other 14 free amino acids (FAAs) due to its slightly enhanced acidic proteases activity. Analyses of transcriptional expressions of salt-tolerant alkaline protease gene (AP, gi: 217809) and aspartyl aminopeptidase gene (AAP, gi: 6165646) indicated that their expression levels were increased by approximately 30 and 27%, respectively. But no mutation was found in the sequences of AP and AAP expression cassettes, suggesting that the increased activities of proteases in mutant H8 should be partially attributed to the increased expression of proteases. ARTP technology showed great potential in enhancing the activities of salt-tolerant proteases from A. oryzae.

The CYPome of the model xenobiotic-biotransforming fungus Cunninghamella elegans

The fungus Cunninghamella elegans is recognised as a microbial model of mammalian drug metabolism owing to its ability to catabolise xenobiotic compounds in an analogous fashion to animals. Its ability to produce phase I (oxidative) metabolites of drugs is associated with cytochrome P450 (CYP) activity; however, almost nothing is known about these enzymes in the fungus. In this paper we report the in silico analysis of the genome sequence of C. elegans B9769, which contains 32 genes putatively coding for CYPs. Based on their predicted amino acid sequences these were classified as belonging to CYP509, 5203, 5208, 5313, 5210, 61 and 51 families. Reverse transcription-quantitative PCR revealed that the gene coding for CYP5313D1 was significantly upregulated when C. elegans DSM1908 was cultivated in sabouraud dextrose in contrast to its expression in cells grown in Roswell Park Memorial Institute medium. This corresponded to the fungus' xenobiotic biotransformation ability when grown in the two media. Heterologous expression of cyp5313D1 in Pichia pastoris resulted in a recombinant strain that biotransformed flurbiprofen to 4'-hydroxyflurbiprofen, the same metabolite generated by C. elegans cultures. This is the first report of a xenobiotic-biotransforming CYP from this biotechnologically important fungus.

Reference genes for qRT-PCR normalisation in different tissues, developmental stages, and stress conditions of Hypericum perforatum

Hypericum perforatum L. is a widely known medicinal herb used mostly as a remedy for depression because it contains high levels of naphthodianthrones, phloroglucinols, alkaloids, and some other secondary metabolites. Quantitative real-time PCR (qRT-PCR) is an optimized method for the efficient and reliable quantification of gene expression studies. In general, reference genes are used in qRT-PCR analysis because of their known or suspected housekeeping roles. However, their expression level cannot be assumed to remain stable under all possible experimental conditions. Thus, the identification of high quality reference genes is essential for the interpretation of qRT-PCR data. In this study, we investigated the expression of 14 candidate genes, including nine housekeeping genes (HKGs) (ACT2, ACT3, ACT7, CYP1, EF1-α, GAPDH, TUB-α, TUB-β, and UBC2) and five potential candidate genes (GSA, PKS1, PP2A, RPL13, and SAND). Three programs-GeNorm, NormFinder, and BestKeeper-were applied to evaluate the gene expression stability across four different plant tissues, four developmental stages and a set of abiotic stress and hormonal treatments. Integrating all of the algorithms and evaluations revealed that ACT2 and TUB-β were the most stable combination in different developmental stages samples and all of the experimental samples. ACT2, TUB-β, and EF1-α were identified as the three most applicable reference genes in different tissues and stress-treated samples. The majority of the conventional HKGs performed better than the potential reference genes. The obtained results will aid in improving the credibility of the standardization and quantification of transcription levels in future expression studies on H. perforatum.

Reference genes for accurate evaluation of expression levels in Trichophyton interdigitale grown under different carbon sources, pH levels and phosphate levels

Tinea pedis is a type of dermatophytosis caused by anthropophilic keratinolytic fungi such as Trichophyton interdigitale. Quantitative reverse transcription PCR (RT-qPCR) is a reliable and reproducible technique for measuring changes in target gene expression across various biological conditions. A crucial aspect of accurate normalization is the choice of appropriate internal controls. To identify reference genes for accurate evaluation of expression levels in T. interdigitale, the transcription levels of eight candidate reference genes (adp-rf, β-act, ef1-α, gapdh, psm1, sdha, rpl2 and ubc) and one target gene (Tri m4) were analysed by RT-qPCR after growing the dermatophyte under different environmental conditions. The results obtained from expression stability evaluations with NormFinder, geNorm, BestKeeper, and RefFinder software demonstrated that adp-rf and psm1 were the most stable internal control genes across all experimental conditions. The present study constitutes the first report of the identification and validation of reference genes for RT-qPCR normalization for T. interdigitale grown under different environmental conditions resembling the conditions encountered by fungi during invasion of skin.

Metabolic comparison of aerial and submerged mycelia formed in the liquid surface culture of Cordyceps militaris

An entomopathogenic fungus, Cordyceps sp. has been known to produce cordycepin which is a purine nucleoside antimetabolite and antibiotic with potential anticancer, antioxidant and anti-inflammatory activities. Interestingly, Cordyceps militaris produces significantly higher amount in a liquid surface culture than in a submerged culture. The liquid surface culture consists of mycelia growing into the air (aerial mycelia) and mycelia growing toward the bottom into the medium (submerged mycelia). In this study, to clarify roles of aerial and submerged mycelia of C. militaris in the cordycepin production the difference in metabolism between these mycelia was investigated. From transcriptomic analyses of the aerial and submerged mycelia at the culture of 5, 12 and 19 days, the metabolism of the submerged mycelia switched from the oxidative phosphorylation to the fermentation pathway. This activated the pentose phosphate pathway to provide building block materials for the nucleotide biosynthetic pathway. Under hypoxic conditions, the 5-aminolevulinic acid synthase (CCM_01504), delta-aminolevulinic acid dehydratase (CCM_00935), coproporphyrinogen III oxidase (CCM_07483) and cytochrome c oxidase 15 (CCM_05057) genes of heme biosynthesis were significantly upregulated. In addition, the liquid surface culture revealed that metabolite coproporhyrinogen III and glycine, the product and precursor of heme, were increased at 12th day and decreased at 19th day, respectively. These results indicate that the submerged mycelia induce the activation of iron acquisition, the ergosterol biosynthetic pathway, and the iron cluster genes of cordycepin biosynthesis in a hypoxic condition. Even though, the expression of the cluster genes of cordycepin biosynthesis was not significantly different in both types of mycelia.

Carbon sources and XlnR-dependent transcriptional landscape of CAZymes in the industrial fungus Talaromyces versatilis: when exception seems to be the rule

Background

Research on filamentous fungi emphasized the remarkable redundancy in genes encoding hydrolytic enzymes, the similarities but also the large differences in their expression, especially through the role of the XlnR/XYR1 transcriptional activator. The purpose of this study was to evaluate the specificities of the industrial fungus Talaromyces versatilis, getting clues into the role of XlnR and the importance of glucose repression at the transcriptional level, to provide further levers for cocktail production.

Results

By studying a set of 62 redundant genes representative of several categories of enzymes, our results underlined the huge plasticity of transcriptional responses when changing nutritional status. As a general trend, the more heterogeneous the substrate, the more efficient to trigger activation. Genetic modifications of xlnR led to significant reorganisation of transcriptional patterns. Just a minimal set of genes actually fitted in a simplistic model of regulation by a transcriptional activator, and this under specific substrates. On the contrary, the diversity of xlnR⁺ versus ΔxlnR responses illustrated the existence of complex and unpredicted patterns of co-regulated genes that were highly dependent on the culture condition, even between genes that encode members of a functional category of enzymes. They notably revealed a dual, substrate-dependant repressor-activator role of XlnR, with counter-intuitive transcripts regulations that targeted specific genes. About glucose, it appeared as a formal repressive sugar as we observed a massive repression of most genes upon glucose addition to the mycelium grown on wheat straw. However, we also noticed a positive role of this sugar on the basal expression of a few genes, (notably those encoding cellulases), showing again the strong dependence of these regulatory mechanisms upon promoter and nutritional contexts.

Conclusions

The diversity of transcriptional patterns appeared to be the rule, while common and stable behaviour, both within gene families and with fungal literature, the exception. The setup of a new biotechnological process to reach optimized, if not customized expression patterns of enzymes, hence appeared tricky just relying on published data that can lead, in the best scenario, to approximate trends. We instead encourage preliminary experimental assays, carried out in the context of interest to reassess gene responses, as a mandatory step before thinking in (genetic) strategies for the improvement of enzyme production in fungi.

Electronic supplementary material

The online version of this article (10.1186/s12934-019-1062-8) contains supplementary material, which is available to authorized users.

Reference genes for accurate normalization of gene expression in wood-decomposing fungi

Wood-decomposing fungi efficiently decompose plant lignocellulose, and there is increasing interest in characterizing and perhaps harnessing the fungal gene regulation strategies that enable wood decomposition. Proper interpretation of these fungal mechanisms relies on accurate quantification of gene expression, demanding reliable internal control genes (ICGs) as references. Commonly used ICGs such as actin, however, fluctuate among wood-decomposing fungi under defined conditions. In this study, by mining RNA-seq data in silico and validating ICGs in vitro using qRT-PCR, we targeted more reliable ICGs for studying transcriptional responses in wood-decomposing fungi, particularly responses to changing environments (e.g., carbon sources, decomposition stages) in various culture conditions. Using the model brown rot fungus Postia placenta in a first-pass study, our mining efforts yielded 15 constitutively-expressed genes robust in variable carbon sources (e.g., no carbon, glucose, cellobiose, aspen) and cultivation stages (e.g., 15 h, 72 h) in submerged cultures. Of these, we found 7 genes as most suitable ICGs. Expression stabilities of these newly selected ICGs were better than commonly used ICGs, analyzed by NormFinder algorithm and qRT-PCR. In a second-pass, multi-species study in solid wood, our RNA-seq mining efforts revealed hundreds of highly constitutively expressed genes among four wood-decomposing fungi with varying nutritional modes (brown rot, white rot), including a shared core set of ICGs numbering 11 genes. Together, the newly selected ICGs highlighted here will increase reliability when studying gene regulatory mechanisms of wood-decomposing fungi.

Enhancing cordycepin production in liquid static cultivation of Cordyceps militaris by adding vegetable oils as the secondary carbon source

This study evaluated different vegetable oils as the second carbon source in liquid static culture of Cordyceps militaris in terms of mycelial growth and cordycepin production. The maximum mycelial concentration and cordycepin production were observed under cottonseed oil and peanut oil induction, respectively. In the condition of adding 20 g/L of peanut oil at Day 0, the final concentration of cordycepin reached to the highest, about 5.29 g/L, which was about 3.17 times higher than that of the control. The qRT-PCR and enzyme activity analysis confirmed that addition of peanut oil up-regulated the expression of the genes encoding glucose-6-phosphate dehydrogenase and isocitrate lyase, as well as the genes in the cordycepin biosynthesis pathway, cns1 and cns2, during the cultivation in C. militaris.

Selection of suitable reference genes for qRT-PCR normalisation under different experimental conditions in Eucommia ulmoides Oliv

Normalisation of data, by choosing the appropriate reference genes, is fundamental for obtaining reliable results in quantitative real-time PCR (qPCR). This study evaluated the expression stability of 11 candidate reference genes with different varieties, developmental periods, tissues, and abiotic stresses by using four statistical algorithms: geNorm, NormFinder, BestKeeper, and RefFinder. The results indicated that ubiquitin-conjugating enzyme S (UBC) and ubiquitin-conjugating enzyme E2 (UBC E2) could be used as reference genes for different E. ulmoides varieties and tissues, UBC and histone H4 (HIS4) for different developmental periods, beta-tubulin (TUB) and UBC for cold treatment, ubiquitin extension protein (UBA80) and HIS4 for drought treatment, and ubiquitin-60S ribosomal protein L40 (UBA52) and UBC E2 for salinity treatment. UBC and UBC E2 for the group “Natural growth” and “Total”, UBA80 and UBC for the group “Abiotic stresses”. To validate the suitability of the selected reference genes in this study, mevalonate kinase (MK), phenylalanine ammonia-lyase (PAL), and 4-coumarate-CoA ligase (4CL) gene expression patterns were analysed. When the most unstable reference genes were used for normalisation, the expression patterns had significant biases compared with the optimum reference gene combinations. These results will be beneficial for more accurate quantification of gene expression levels in E. ulmoides.

Identification of suitable internal control genes for transcriptional studies in Eleusine coracana under different abiotic stress conditions

Finger millet [Eleusine coracana (L.) Gaertn] is an excellent food and forage crop of arid and semiarid areas in Africa and Asia. It is well adapted to drought, heat, high salinity, poor soil fertility and low pH with an efficient C₄ carbon fixation mechanism for high yield potential. To normalize the target gene expression data, the identification of suitable reference genes is essential. Ten candidate reference genes were selected and their expression stability was analyzed in various samples treated with different abiotic stress conditions. Five different statistical algorithms: geNorm, NormFinder, BestKeeper, ΔCt, and RefFinder were used to determine the stability of these genes. Our results revealed GAPDH, EEF1a, ACT and CYC as highly stable reference genes and PP2A and eIF4A as least stable reference genes across all the samples and suggesting that these genes could be used for accurate transcript normalization under abiotic stress. To the best of our knowledge, this is the first report on identification of suitable reference genes for accurate transcript normalization using qRT-PCR in finger millet under abiotic stress.

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.

Identification and application of exogenous dsRNA confers plant protection against Sclerotinia sclerotiorum and Botrytis cinerea

Sclerotinia sclerotiorum, the causal agent of white stem rot, is responsible for significant losses in crop yields around the globe. While our understanding of S. sclerotiorum infection is becoming clearer, genetic control of the pathogen has been elusive and effective control of pathogen colonization using traditional broad-spectrum agro-chemical protocols are less effective than desired. In the current study, we developed species-specific RNA interference-based control treatments capable of reducing fungal infection. Development of a target identification pipeline using global RNA sequencing data for selection and application of double stranded RNA (dsRNA) molecules identified single gene targets of the fungus. Using this approach, we demonstrate the utility of this technology through foliar applications of dsRNAs to the leaf surface that significantly decreased fungal infection and S. sclerotiorum disease symptoms. Select target gene homologs were also tested in the closely related species, Botrytis cinerea, reducing lesion size and providing compelling evidence of the adaptability and flexibility of this technology in protecting plants against devastating fungal pathogens.

Transcription of lignocellulose-decomposition associated genes, enzyme activities and production of ethanol upon bioconversion of waste substrate by Phlebia radiata

Previously identified twelve plant cell wall degradation-associated genes of the white rot fungus Phlebia radiata were studied by RT-qPCR in semi-aerobic solid-state cultures on lignocellulose waste material, and on glucose-containing reference medium. Wood-decay-involved enzyme activities and ethanol production were followed to elucidate both the degradative and fermentative processes. On the waste lignocellulose substrate, P. radiata carbohydrate-active enzyme (CAZy) genes encoding cellulolytic and hemicellulolytic activities were significantly upregulated whereas genes involved in lignin modification displayed a more complex response. Two lignin peroxidase genes were differentially expressed on waste lignocellulose compared to glucose medium, whereas three manganese peroxidase-encoding genes were less affected. On the contrary, highly significant difference was noticed for three cellulolytic genes (cbhI_1, eg1, bgl1) with higher expression levels on the lignocellulose substrate than on glucose. This indicates expression of the wood-attacking degradative enzyme system by the fungus also on the recycled, waste core board material. During the second week of cultivation, ethanol production increased on the core board to 0.24 g/L, and extracellular activities against cellulose, xylan, and lignin were detected. Sugar release from the solid lignocellulose resulted with concomitant accumulation of ethanol as fermentation product. Our findings confirm that the fungus activates its white rot decay system also on industrially processed lignocellulose adopted as growth substrate, and under semi-aerobic cultivation conditions. Thus, P. radiata is a good candidate for lignocellulose-based renewable biotechnology to make biofuels and biocompounds from materials with less value for recycling or manufacturing.

Genomic Organization and Expression of Iron Metabolism Genes in the Emerging Pathogenic Mold Scedosporium apiospermum

The ubiquitous mold Scedosporium apiospermum is increasingly recognized as an emerging pathogen, especially among patients with underlying disorders such as immunodeficiency or cystic fibrosis (CF). Indeed, it ranks the second among the filamentous fungi colonizing the respiratory tract of CF patients. However, our knowledge about virulence factors of this fungus is still limited. The role of iron-uptake systems may be critical for establishment of Scedosporium infections, notably in the iron-rich environment of the CF lung. Two main strategies are employed by fungi to efficiently acquire iron from their host or from their ecological niche: siderophore production and reductive iron assimilation (RIA) systems. The aim of this study was to assess the existence of orthologous genes involved in iron metabolism in the recently sequenced genome of S. apiospermum. At first, a tBLASTn analysis using A. fumigatus iron-related proteins as query revealed orthologs of almost all relevant loci in the S. apiospermum genome. Whereas the genes putatively involved in RIA were randomly distributed, siderophore biosynthesis and transport genes were organized in two clusters, each containing a non-ribosomal peptide synthetase (NRPS) whose orthologs in A. fumigatus have been described to catalyze hydroxamate siderophore synthesis. Nevertheless, comparative genomic analysis of siderophore-related clusters showed greater similarity between S. apiospermum and phylogenetically close molds than with Aspergillus species. The expression level of these genes was then evaluated by exposing conidia to iron starvation and iron excess. The expression of several orthologs of A. fumigatus genes involved in siderophore-based iron uptake or RIA was significantly induced during iron starvation, and conversely repressed in iron excess conditions. Altogether, these results indicate that S. apiospermum possesses the genetic information required for efficient and competitive iron uptake. They also suggest an important role of the siderophore production system in iron uptake by S. apiospermum.

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.

Proteasome storage granules protect proteasomes from autophagic degradation upon carbon starvation

26S proteasome abundance is tightly regulated at multiple levels, including the elimination of excess or inactive particles by autophagy. In yeast, this proteaphagy occurs upon nitrogen starvation but not carbon starvation, which instead stimulates the rapid sequestration of proteasomes into cytoplasmic puncta termed proteasome storage granules (PSGs). Here, we show that PSGs help protect proteasomes from autophagic degradation. Both the core protease and regulatory particle sub-complexes are sequestered separately into PSGs via pathways dependent on the accessory proteins Blm10 and Spg5, respectively. Modulating PSG formation, either by perturbing cellular energy status or pH, or by genetically eliminating factors required for granule assembly, not only influences the rate of proteasome degradation, but also impacts cell viability upon recovery from carbon starvation. PSG formation and concomitant protection against proteaphagy also occurs in Arabidopsis, suggesting that PSGs represent an evolutionarily conserved cache of proteasomes that can be rapidly re-mobilized based on energy availability.

Proteins perform many jobs within an organism, including providing structure and support, and protecting against infection. The levels of the many proteins in a cell need to be carefully controlled so that the correct amounts are present at the right place and time to perform these tasks. This control can be achieved by balancing the production of new proteins with the break down (or degradation) of proteins that are no longer required or become dysfunctional.

Most cells have two pathways for degrading proteins. One pathway breaks down individual proteins specifically marked for elimination; this causes them to be recognized by a structure called the proteasome, which chops proteins into smaller pieces. Larger protein assemblies – including the proteasome itself – are to big for the proteasome and thus need to be degraded by another pathway called autophagy. This process engulfs and delivers parts of a cell to a membrane-bound compartment called the vacuole, which ‘digests’ and recycles these larger constituents.

Proteasomes are degraded by autophagy when they are not working correctly and when nitrogen (a crucial nutrient) is in short supply. However, proteasomes are not degraded when cells lack carbon, even though this starvation is known to activate autophagy in the same way that an absence of nitrogen does. So how do proteasomes escape degradation when cells are starved for carbon?

Marshall and Vierstra now show that upon carbon starvation, proteasomes rapidly exit the cell nucleus and cluster together in the main part of the cell (termed the cytosol). These clusters are known as proteasome storage granules (PSGs). In fungi and plants, mutations or conditions inside the cell that make it difficult for PSGs to assemble cause proteasomes to instead be broken down in the vacuole when carbon availability is low. Clustering into PSGs therefore protects proteasomes from autophagy. This clustering appears advantageous to cells; yeast cells that could form PSGs were better able to start growing again when their nutrient supply improved.

Protein clustering (also known as aggregation) is an important strategy that cells use to survive stressful conditions. However, it can also be harmful when proteins aggregate inappropriately, such as occurs in Alzheimer’s disease. Researchers may be able to use PSG assembly as a convenient model to study the causes and consequences of protein aggregation; this knowledge could ultimately be applied to improve human health and crop productivity.

Identifying optimal reference genes for the normalization of microRNA expression in cucumber under viral stress

Cucumber green mottle mosaic virus (CGMMV) is an economically important pathogen and causes significant reduction of both yield and quality of cucumber (Cucumis sativus). Currently, there were no satisfied strategies for controlling the disease. A better understanding of microRNA (miRNA) expression related to the regulation of plant-virus interactions and virus resistance would be of great assistance when developing control strategies for CGMMV. However, accurate expression analysis is highly dependent on robust and reliable reference gene used as an internal control for normalization of miRNA expression. Most commonly used reference genes involved in CGMMV-infected cucumber are not universally expressed depending on tissue types and stages of plant development. It is therefore crucial to identify suitable reference genes in investigating the role of miRNA expression. In this study, seven reference genes, including Actin, Tubulin, EF-1α, 18S rRNA, Ubiquitin, GAPDH and Cyclophilin, were evaluated for the most accurate results in analyses using reverse transcription-quantitative polymerase chain reaction (RT-qPCR). Gene expression was assayed on cucumber leaves, stems and roots that were collected at different days post inoculation with CGMMV. The expression data were analyzed using algorithms including delta-Ct, geNorm, NormFinder, and BestKeeper as well as the comparative tool RefFinder. The reference genes were subsequently validated using miR159. The results showed that EF-1α and GAPDH were the most reliable reference genes for normalizing miRNA expression in leaf, root and stem samples, while Ubiquitin and EF-1α were the most suitable combination overall.

Stability evaluation of reference genes for gene expression analysis by RT-qPCR in soybean under different conditions

Real-time quantitative reverse transcription PCR is a sensitive and widely used technique to quantify gene expression. To achieve a reliable result, appropriate reference genes are highly required for normalization of transcripts in different samples. In this study, 9 previously published reference genes (60S, Fbox, ELF1A, ELF1B, ACT11, TUA5, UBC4, G6PD, CYP2) of soybean [Glycine max (L.) Merr.] were selected. The expression stability of the 9 genes was evaluated under conditions of biotic stress caused by infection with soybean mosaic virus, nitrogen stress, across different cultivars and developmental stages. ΔCt and geNorm algorithms were used to evaluate and rank the expression stability of the 9 reference genes. Results obtained from two algorithms showed high consistency. Moreover, results of pairwise variation showed that two reference genes were sufficient to normalize the expression levels of target genes under each experimental setting. For virus infection, ELF1A and ELF1B were the most stable reference genes for accurate normalization. For different developmental stages, Fbox and G6PD had the highest expression stability between two soybean cultivars (Tanlong No. 1 and Tanlong No. 2). ELF1B and ACT11 were identified as the most stably expressed reference genes both under nitrogen stress and among different cultivars. The results showed that none of the candidate reference genes were uniformly expressed at different conditions, and selecting appropriate reference genes was pivotal for gene expression studies with particular condition and tissue. The most stable combination of genes identified in this study will help to achieve more accurate and reliable results in a wide variety of samples in soybean.

Selection of reference genes from Shiraia bambusicola for RT-qPCR analysis under different culturing conditions

Stable reference genes are necessary to analyse quantitative real-time reverse transcription PCR (qRT-PCR) data and determine the reliability of the final results. For further studies of the valuable fungus Shiraia bambusicola, the identification of suitable reference genes has become increasingly urgent. In this study, three conventional reference genes and nine novel candidates were evaluated under different light conditions (all-dark, all-light and 12-h light/dark) and in different media (rice medium, PD medium, and Czapek–Dox medium). Three popular software programs (geNorm, NormFinder and BestKeeper) were used to analyse these genes, and the final ranking was determined using RefFinder. SbLAlv9, SbJsn1, SbSAS1 and SbVAC55 displayed the best stability among the genes, while SbFYVE and SbPKI showed the worst. These emerging genes exhibited significantly better properties than the three existing genes under almost all conditions. Furthermore, the most reliable reference genes were identified separately under different nutrient and light conditions, which would help accessible to make the most of the existing data. In summary, a group of novel housekeeping genes from S. bambusicola with more stable properties than before was explored, and these results could also provide a practical approach for other filamentous fungi.

Daily rhythms and enrichment patterns in the transcriptome of the behavior-manipulating parasite Ophiocordyceps kimflemingiae

Various parasite-host interactions that involve adaptive manipulation of host behavior display time-of-day synchronization of certain events. One example is the manipulated biting behavior observed in Carpenter ants infected with Ophiocordyceps unilateralis sensu lato. We hypothesized that biological clocks play an important role in this and other parasite-host interactions. In order to identify candidate molecular clock components, we used two general strategies: bioinformatics and transcriptional profiling. The bioinformatics approach was used to identify putative homologs of known clock genes. For transcriptional profiling, RNA-Seq was performed on 48 h time courses of Ophiocordyceps kimflemingiae (a recently named species of the O. unilateralis complex), whose genome has recently been sequenced. Fungal blastospores were entrained in liquid media under 24 h light-dark (LD) cycles and were harvested at 4 h intervals either under LD or continuous darkness. Of all O. kimflemingiae genes, 5.3% had rhythmic mRNAs under these conditions (JTK Cycle, ≤ 0.057 statistical cutoff). Our data further indicates that a significant number of transcription factors have a peaked activity during the light phase (day time). The expression levels of a significant number of secreted enzymes, proteases, toxins and small bioactive compounds peaked during the dark phase or subjective night. These findings support a model whereby this fungal parasite uses its biological clock for phase-specific activity. We further suggest that this may be a general mechanism involved in parasite-host interactions.

Insight into cordycepin biosynthesis of Cordyceps militaris: Comparison between a liquid surface culture and a submerged culture through transcriptomic analysis

Cordyceps militaris produces cordycepin, which is known to be a bioactive compound. Currently, cordycepin hyperproduction of C. militaris was carried out in a liquid surface culture because of its low productivity in a submerged culture, however the reason was not known. In this study, 4.92 g/L of cordycepin was produced at the 15^th day of C. militaris NBRC 103752 liquid surface culture, but only 1 mg/L was produced in the submerged culture. RNA-Seq was used to clarify the gene expression profiles of the cordycepin biosynthetic pathways of the submerged culture and the liquid surface culture. From this analysis, 1036 genes were shown to be upregulated and 557 genes were downregulated in the liquid surface culture compared with the submerged culture. Specifically, adenylosuccinate synthetase and phosphoribosylaminoimidazole-succinocarboxamide (SAICAR) synthase in purine nucleotide metabolism were significantly upregulated in the liquid surface culture. Thick mycelia formation in the liquid surface culture was found to induce the expression of hypoxia-related genes (GABA shunt, glutamate synthetase precursor, and succinate-semialdehyde dehydrogenase). Cytochrome P450 oxidoreductases containing heme were also found to be significantly enriched, suggesting that a hypoxic condition might be created in the liquid surface culture. These results suggest that hypoxic conditions are more suitable for cordycepin production in the liquid surface culture compared with the submerged culture. Our analysis paves the way for unraveling the cordycepin biosynthesis pathway and for improving cordycepin production in C. militaris.

Identification and validation of reference genes for qRT-PCR studies of the obligate aphid pathogenic fungus Pandora neoaphidis during different developmental stages

The selection of stable reference genes is a critical step for the accurate quantification of gene expression. To identify and validate the reference genes in Pandora neoaphidis–an obligate aphid pathogenic fungus—the expression of 13classical candidate reference genes were evaluated by quantitative real-time reverse transcriptase polymerase chain reaction(qPCR) at four developmental stages (conidia, conidia with germ tubes, short hyphae and elongated hyphae). Four statistical algorithms, including geNorm, NormFinder, BestKeeper and Delta Ct method were used to rank putative reference genes according to their expression stability and indicate the best reference gene or combination of reference genes for accurate normalization. The analysis of comprehensive ranking revealed that ACT1and 18Swas the most stably expressed genes throughout the developmental stages. To further validate the suitability of the reference genes identified in this study, the expression of cell division control protein 25 (CDC25) and Chitinase 1(CHI1) genes were used to further confirm the validated candidate reference genes. Our study presented the first systematic study of reference gene(s) selection for P. neoaphidis study and provided guidelines to obtain more accurate qPCR results for future developmental efforts.

Heavy Metals Induce Iron Deficiency Responses at Different Hierarchic and Regulatory Levels

In plants, the excess of several heavy metals mimics iron (Fe) deficiency-induced chlorosis, indicating a disturbance in Fe homeostasis. To examine the level at which heavy metals interfere with Fe deficiency responses, we carried out an in-depth characterization of Fe-related physiological, regulatory, and morphological responses in Arabidopsis (Arabidopsis thaliana) exposed to heavy metals. Enhanced zinc (Zn) uptake closely mimicked Fe deficiency by leading to low chlorophyll but high ferric-chelate reductase activity and coumarin release. These responses were not caused by Zn-inhibited Fe uptake via IRON-REGULATED TRANSPORTER (IRT1). Instead, Zn simulated the transcriptional response of typical Fe-regulated genes, indicating that Zn affects Fe homeostasis at the level of Fe sensing. Excess supplies of cobalt and nickel altered root traits in a different way from Fe deficiency, inducing only transient Fe deficiency responses, which were characterized by a lack of induction of the ethylene pathway. Cadmium showed a rather inconsistent influence on Fe deficiency responses at multiple levels. By contrast, manganese evoked weak Fe deficiency responses in wild-type plants but strongly exacerbated chlorosis in irt1 plants, indicating that manganese antagonized Fe mainly at the level of transport. These results show that the investigated heavy metals modulate Fe deficiency responses at different hierarchic and regulatory levels and that the interaction of metals with physiological and morphological Fe deficiency responses is uncoupled. Thus, this study not only emphasizes the importance of assessing heavy metal toxicities at multiple levels but also provides a new perspective on how Fe deficiency contributes to the toxic action of individual heavy metals.

Selection of reference genes for quantitative real-time RT-PCR assays in different morphological forms of dimorphic zygomycetous fungus Benjaminiella poitrasii

Benjaminiella poitrasii, a dimorphic non-pathogenic zygomycetous fungus, exhibits a morphological yeast (Y) to hypha (H) reversible transition in the vegetative phase, sporangiospores (S) in the asexual phase and zygospores (Z) in the sexual phase. To study the gene expression across these diverse morphological forms, suitable reference genes are required. In the present study, 13 genes viz. ACT, 18S rRNA, eEF1α, eEF-Tu,eIF-1A, Tub-α, Tub-b, Ubc, GAPDH, Try, WS-21, NADGDH and NADPGDH were evaluated for their potential as a reference, particularly for studying gene expression during the Y-H reversible transition and also for other asexual and sexual life stages of B. poitrasii. Analysis of RT-qPCR data using geNorm, normFinder and BestKeeper software revealed that genes such as Ubc, 18S rRNA and WS-21 were expressed at constant levels in each given subset of RNA samples from all the morphological phases of B. poitrasii. Therefore, these reference genes can be used to elucidate the role of morpho-genes in B. poitrasii. Further, use of the two most stably expressed genes (Ubc and WS-21) to normalize the expression of the ornithine decarboxylase gene (Bpodc) in different morphological forms of B. poitrasii, generated more reliable results, indicating that our selection of reference genes was appropriate.

Selection and Validation of Reference Genes for Accurate RT-qPCR Data Normalization in Coffea spp. under a Climate Changes Context of Interacting Elevated [CO2] and Temperature

World coffee production has faced increasing challenges associated with ongoing climatic changes. Several studies, which have been almost exclusively based on temperature increase, have predicted extensive reductions (higher than half by 2,050) of actual coffee cropped areas. However, recent studies showed that elevated [CO₂] can strongly mitigate the negative impacts of heat stress at the physiological and biochemical levels in coffee leaves. In addition, it has also been shown that coffee genotypes can successfully cope with temperatures above what has been traditionally accepted. Altogether, this information suggests that the real impact of climate changes on coffee growth and production could be significantly lower than previously estimated. Gene expression studies are an important tool to unravel crop acclimation ability, demanding the use of adequate reference genes. We have examined the transcript stability of 10 candidate reference genes to normalize RT-qPCR expression studies using a set of 24 cDNAs from leaves of three coffee genotypes (CL153, Icatu, and IPR108), grown under 380 or 700 μL CO₂ L^-1, and submitted to increasing temperatures from 25/20°C (day/night) to 42/34°C. Samples were analyzed according to genotype, [CO₂], temperature, multiple stress interaction ([CO₂], temperature) and total stress interaction (genotype, [CO₂], and temperature). The transcript stability of each gene was assessed through a multiple analytical approach combining the Coeficient of Variation method and three algorithms (geNorm, BestKeeper, NormFinder). The transcript stability varied according to the type of stress for most genes, but the consensus ranking obtained with RefFinder, classified MDH as the gene with the highest mRNA stability to a global use, followed by ACT and S15, whereas α-TUB and CYCL showed the least stable mRNA contents. Using the coffee expression profiles of the gene encoding the large-subunit of ribulose-1,5-bisphosphate carboxylase/oxygenase (RLS), results from the in silico aggregation and experimental validation of the best number of reference genes showed that two reference genes are adequate to normalize RT-qPCR data. Altogether, this work highlights the importance of an adequate selection of reference genes for each single or combined experimental condition and constitutes the basis to accurately study molecular responses of Coffea spp. in a context of climate changes and global warming.

Highly active promoters and native secretion signals for protein production during extremely low growth rates in Aspergillus niger

Background

The filamentous ascomycete Aspergillus niger is used in many industrial processes for the production of enzymes and organic acids by batch and fed-batch cultivation. An alternative technique is continuous cultivation, which promises improved yield and optimized pipeline efficiency.

Results

In this work, we have used perfusion (retentostat) cultivation to validate two promoters that are suitable for A. niger continuous cultivation of industrially relevant products. Firstly, promoters of genes encoding either an antifungal protein (Panafp) or putative hydrophobin (PhfbD) were confirmed as active throughout retentostat culture by assessing mRNA and protein levels using a luciferase (mluc) reporter system. This demonstrated the anafp promoter mediates a high but temporally variable expression profile, whereas the hfbD promoter mediates a semi-constant, moderate-to-high protein expression during retentostat culture. In order to assess whether these promoters were suitable to produce heterologous proteins during retentostat cultivation, the secreted antifungal protein (AFP) from Aspergillus giganteus, which has many potential biotechnological applications, was expressed in A. niger during retentostat cultivation. Additionally, this assay was used to concomitantly validate that native secretion signals encoded in anafp and hfbD genes can be harnessed for secretion of heterologous proteins. Afp mRNA and protein abundance were comparable to luciferase measurements throughout retentostat cultivation, validating the use of Panafp and PhfbD for perfusion cultivation. Finally, a gene encoding the highly commercially relevant thermal hysteresis protein (THP) was expressed in this system, which did not yield detectable protein.

Conclusion

Both hfbD and anafp promoters are suitable for production of useful products in A. niger during perfusion cultivation. These findings provide a platform for further optimisations for high production of heterologous proteins with industrial relevance.

Electronic supplementary material

The online version of this article (doi:10.1186/s12934-016-0543-2) contains supplementary material, which is available to authorized users.

Identification of novel and robust internal control genes from Volvariella volvacea that are suitable for RT-qPCR in filamentous fungi

The selection of appropriate internal control genes (ICGs) is a crucial step in the normalization of real-time quantitative PCR (RT-qPCR) data. Housekeeping genes are habitually selected for this purpose, despite accumulating evidence on their instability. We screened for novel, robust ICGs in the mushroom forming fungus Volvariella volvacea. Nine commonly used and five newly selected ICGs were evaluated for expression stability using RT-qPCR data in eight different stages of the life cycle of V. volvacea. Three different algorithms consistently determined that three novel ICGs (SPRYp, Ras and Vps26) exhibited the highest expression stability in V. volvacea. Subsequent analysis of ICGs in twenty-four expression profiles from nine filamentous fungi revealed that Ras was the most stable ICG amongst the Basidiomycetous samples, followed by SPRYp, Vps26 and ACTB. Vps26 was expressed most stably within the analyzed data of Ascomycetes, followed by HH3 and β-TUB. No ICG was universally stable for all fungal species, or for all experimental conditions within a species. Ultimately, the choice of an ICG will depend on a specific set of experiments. This study provides novel, robust ICGs for Basidiomycetes and Ascomycetes. Together with the presented guiding principles, this enables the efficient selection of suitable ICGs for RT-qPCR.

Algal carbohydrates affect polyketide synthesis of the lichen-forming fungus Cladonia rangiferina

Lichen secondary metabolites (polyketides) are produced by the fungal partner, but the role of algal carbohydrates in polyketide biosynthesis is not clear. This study examined whether the type and concentration of algal carbohydrate explained differences in polyketide production and gene transcription by a lichen fungus (Cladonia rangiferina). The carbohydrates identified from a free-living cyanobacterium (Spirulina platensis; glucose), a lichen-forming alga (Diplosphaera chodatii; sorbitol) and the lichen alga that associates with C. rangiferina (Asterochloris sp.; ribitol) were used in each of 1%, 5% and 10% concentrations to enrich malt yeast extract media for culturing the mycobiont. Polyketides were determined by high performance liquid chromatography (HPLC), and polyketide synthase (PKS) gene transcription was measured by quantitative PCR of the ketosynthase domain of four PKS genes. The lower concentrations of carbohydrates induced the PKS gene expression where ribitol up-regulated CrPKS1 and CrPKS16 gene transcription and sorbitol up-regulated CrPKS3 and CrPKS7 gene transcription. The HPLC results revealed that lower concentrations of carbon sources increased polyketide production for three carbohydrates. One polyketide from the natural lichen thallus (fumarprotocetraric acid) also was produced by the fungal culture in ribitol supplemented media only. This study provides a better understanding of the role of the type and concentration of the carbon source in fungal polyketide biosynthesis in the lichen Cladonia rangiferina.

Quantitative real-time PCR normalization for gene expression studies in the plant pathogenic fungi Lasiodiplodia theobromae

Lasiodiplodia theobromae is a highly virulent plant pathogen. It has been suggested that heat stress increases its virulence. The aim of this work was to evaluate, compare, and recommend normalization strategies for gene expression analysis of the fungus growing with grapevine wood under heat stress. Using RT-qPCR-derived data, reference gene stability was evaluated through geNorm, NormFinder and Bestkeeper applications. Based on the geometric mean using the ranking position obtained for each independent analysis, genes were ranked from least to most stable as follows: glyceraldehyde-3-phosphate dehydrogenase (GAPDH), actin (ACT), β-tubulin (TUB) and elongation factor-1α (EF1α). Using RNAseq-derived data based on the calculated tagwise dispersion these genes were ordered by increasing stability as follows: GAPDH, ACT, TUB, and EF1α. The correlation between RNAseq and RTqPCR results was used as criteria to identify the best RT-qPCR normalization approach. The gene TUB is recommended as the best option for normalization among the commonly used reference genes, but alternative fungal reference genes are also suggested.

Comparative Transcriptomic Approaches Exploring Contamination Stress Tolerance in Salix sp. Reveal the Importance for a Metaorganismal de Novo Assembly Approach for Nonmodel Plants

Metatranscriptomic study of nonmodel organisms requires strategies that retain the highly resolved genetic information generated from model organisms while allowing for identification of the unexpected. A real-world biological application of phytoremediation, the field growth of 10 Salix cultivars on polluted soils, was used as an exemplar nonmodel and multifaceted crop response well-disposed to the study of gene expression. Sequence reads were assembled de novo to create 10 independent transcriptomes, a global transcriptome, and were mapped against the Salix purpurea 94006 reference genome. Annotation of assembled contigs was performed without a priori assumption of the originating organism. Global transcriptome construction from 3.03 billion paired-end reads revealed 606,880 unique contigs annotated from 1588 species, often common in all 10 cultivars. Comparisons between transcriptomic and metatranscriptomic methodologies provide clear evidence that nonnative RNA can mistakenly map to reference genomes, especially to conserved regions of common housekeeping genes, such as actin, α/β-tubulin, and elongation factor 1-α. In Salix, Rubisco activase transcripts were down-regulated in contaminated trees across all 10 cultivars, whereas thiamine thizole synthase and CP12, a Calvin Cycle master regulator, were uniformly up-regulated. De novo assembly approaches, with unconstrained annotation, can improve data quality; care should be taken when exploring such plant genetics to reduce de facto data exclusion by mapping to a single reference genome alone. Salix gene expression patterns strongly suggest cultivar-wide alteration of specific photosynthetic apparatus and protection of the antenna complexes from oxidation damage in contaminated trees, providing an insight into common stress tolerance strategies in a real-world phytoremediation system.