Next Article in Journal
Emerging Roles of Exocyst Complex in Fungi: A Review
Previous Article in Journal
Metarhizium caribense sp. nov., a Novel Species of Entomopathogenic Metarhizium Fungi Associated with Weevils Impairing Coffee, Sugar Cane and Sweet Potato Cultivation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Rad6 Regulates Conidiation by Affecting the Biotin Metabolism in Beauveria bassiana

Fujian Key Laboratory of Marine Enzyme Engineering, Fuzhou University, Fuzhou 350108, China
*
Authors to whom correspondence should be addressed.
J. Fungi 2024, 10(9), 613; https://doi.org/10.3390/jof10090613
Submission received: 15 July 2024 / Revised: 24 August 2024 / Accepted: 25 August 2024 / Published: 28 August 2024

Abstract

:
Rad6 is a canonical ubiquitin-conjugating enzyme known for its role in regulating chromosome-related cellular processes in yeast and has been proven to have multiple functions in Beauveria bassiana, including insect-pathogenic lifestyle, UV damage repair, and conidiation. However, previous studies have only reported the key role of Rad6 in regulating conidial production in a nutrient-rich medium, without any deep mechanism analyses. In this study, we found that the disruption of Rad6 leads to a profound reduction in conidial production, irrespective of whether the fungus is cultivated in nutrient-rich or nutrient-poor environments. The absence of rad6 exerts a suppressive effect on the transcription of essential genes in the central developmental pathway, namely, brlA, abaA, and wetA, resulting in a direct downregulation of conidiation capacity. Additionally, mutant strains exhibited a more pronounced decline in both conidial generation and hyphal development when cultured in nutrient-rich conditions. This observation correlates with the downregulation of the central developmental pathway (CDP) downstream gene vosA and the upregulation of flaA in nutrient-rich cultures. Moreover, single-transcriptomics analyses indicated that irregularities in biotin metabolism, DNA repair, and tryptophan metabolism are the underlying factors contributing to the reduced conidial production. Comprehensive dual transcriptomics analyses pinpointed abnormal biotin metabolism as the primary cause of conidial production decline. Subsequently, we successfully restored conidial production in the Rad6 mutant strain through the supplementation of biotin, further confirming the transcriptomic evidence. Altogether, our findings underscore the pivotal role of Rad6 in influencing biotin metabolism, subsequently impacting the expression of CDP genes and ultimately shaping the asexual life cycle of B. bassiana.

1. Introduction

The ubiquitin post-modification system comprises ubiquitin-activating (E1), ubiquitin-conjugating (E2), and ubiquitin-ligating (E3) enzymes, which transfer ubiquitin molecules to substrates [1]. Ubiquitin-marked substrates typically undergo selective degradation, influencing their nature, function, and fate [1,2]. E2 proteins act as a bridge between E1 and E3 enzymes, regulating ubiquitin chain initiation, formation, and elongation. They determine the direction of ubiquitination, ranging from monoubiquitination to multiubiquitination, thereby dictating the fate of the modified proteins [3].
Biotin, also known as vitamin B7, vitamin H, and coenzyme R, is a water-soluble vitamin that interacts directly with a wide range of enzymes and is one of the essential cofactors involved in the central pathways of metabolism in prokaryotic and eukaryotic cells [4]. As one of the vital vitamins for all living organisms, biotin plays a crucial role in many important biological processes, such as biotin-dependent carboxylases, fatty acid biosynthesis, gluconeogenesis and amino acid metabolism, cell signaling, gene expression, and chromatin structure [4,5,6]. The increase or decrease in biotin affects the life activities of organisms, e.g., biotin is vital for the growth, virulence, and spore production of Alternaria alternata [7]. Biotin can enhance the osmotic tolerance of Torulopsis mogii by promoting the accumulation of alginate, increasing the K+/Na+ ratio, promoting ATPase activity, and increasing the synthesis of long-chain highly saturated fatty acids [8]. Biotin deficiency affects the expression of genes that synthesize biotin in yeast as well as precursor transport genes and the transcription factor Vhr1p gene [9]. Biotin influences Candida proliferative capacity and virulence by regulating the expression of the biotin transcription factor Vhr1 and the biotin transporter protein Vht1 and, thereby, affecting the proliferative capacity and virulence of Candida [10].
In filamentous fungi, the central developmental pathway (CDP) plays a pivotal role in influencing fungal asexual development by regulating the expression of downstream conidiation-specific genes through cascade regulation. BrlA, AbaA, and WetA are three key regulators in the CDP [11,12]. In B. bassiana, the disruption of brlA and abaA not only results in the elimination of conidia and blastospore production, but also inhibits its dimorphic transition [13]. WetA and its downstream velvet protein, VosA, are associated with conidial maturation [14]. In addition, flb genes, including flbA-E, are essential for CDP activation by mediating the transcription of brlA and abaA [15]. These findings indicate that the CDP has a complex regulatory network and is indispensable in the asexual development of B. bassiana.
The rad6 gene encodes a ubiquitin-conjugating enzyme in the model fungus Saccharomyces cerevisiae [16]. The disruption of Rad6 leads to multiple deficiencies in S. cerevisiae, affecting DNA repair, selective proteolysis, checkpoints, meiosis, and sporulation [17,18,19]. Entomopathogenic fungi are employed globally as biocontrol agents, yet few researchers have explored the role of the ubiquitin-conjugating enzyme Rad6 in their life cycles. Due to the role of Rad6 in yeast and the multiple functions of biotin in fungi, this raises the following question: would biotin have any relationship with Rad6 in regulating biological potential in B. bassiana? We hypothesized that the orthologous Rad6 in B. bassiana may play a role with biotin in conidiation similar to that in yeast. Transcriptomics is a widely utilized method in the study of filamentous fungi. Hou et al. applied transcriptomic investigation to reveal the survival mechanism for B. bassiana under linoleic acid stress [20]. Wang et al. uncovered the thousands of differentially expressed genes involved in a central developmental pathway and secondary metabolism under light or dark incubation conditions by transcriptomic analysis [21]. This study aimed to test this hypothesis by comparing the phenotypes and transcriptomic patterns of Rad6 mutants and wild types (WTs). As shown below, Rad6 mediates the conidiation capability of B. bassiana, possibly through the regulation of biotin metabolism.

2. Materials and Methods

2.1. Bioinformatic Analysis of Fungal Rad6 Homologs

The standard S. cerevisiae Rad6 (NP_011457) served as the bait in search of the genomes of B. bassiana [22], as well as certain human and plant pathogenic fungi, and entomopathogenic fungi at http://blast.ncbi.nlm.nih.gov/Blast.cgi/ (accessed on 8 October 2019). A neighbor-joining method in MEGA7 was used to reveal the phylogenetic relationship between B. bassiana Ras6 and its orthologs [23].

2.2. Generation of rad6 Mutants

The WT strain of B. bassiana ARSEF 2860 (Bb2860) used in this study was obtained from the U.S. Department of Agriculture Entomopathogenic Fungus Collection (Ithaca, NY, USA). The rad6 disruption strain, Δrad6, was generated by homologous recombination in the wild-type B. bassiana strain ARSEF 2860. This process involved separating its coding/flanking fragments with a bar marker, a resistant gene against phosphinothricin. The complementary strain (Δrad6::rad6) was designed by ectopic integration of a cassette containing the entire open reading frame (ORF) of the rad6 gene and a sur marker, a resistant gene against chlorimuron ethyl. Briefly, the upstream (1541 bp) and downstream (1553 bp) fragments of the rad6 ORF were amplified and ligated to the restriction enzyme sites EcoRI/BamHI and XbaI/SpeI of the p0380-GFP-bar vector using the ClonExpress II One Step Cloning Kit (Vazyme Biotech, Nanjing, China). For the complementary strain, the promoter, coding sequence, and downstream regions were amplified and inserted into p0380-sur-gateway, resulting in p0380-sur-rad6. The disruption plasmid p0380-GFP-rad6-bar and complementation plasmid p0380-sur-rad6 were transformed into the WT and Δrad6, respectively, by using Agrobacterium-mediated transformation method, following the procedures detailed in our previous studies [1,24,25]. The putative disruption transformants were initially screened on CDA plates (200 μg/mL phosphinothricin or 15 μg/mL chlorimuron ethyl), followed by a secondary LSCM screening to identify non-fluorescent strains. Positive transformants were subjected to PCR verification. The primers used for generating the mutants are shown in Supplemental File S1, Table S1. The results of the verification of the mutants are described in our previous publication [26].

2.3. Cultivation Conditions

The strains were cultured in Sabouraud dextrose agar (SDAY) medium (4% glucose, 1% peptone, 1% yeast extract, and 1.5% agar) at 25 °C on a 12:12 h light/dark cycle for 8 days, and then the concentration of conidia was adjusted to 107 conidia mL−1 for subsequent experiments. SDAY and the nutrient-poor medium, Czapek-Dox agar (CDA) (3% sucrose, 0.3% NaNO3, 0.1% K2HPO4, 0.05% KCl, 0.05% MgSO4, and 0.001% FeSO4, plus 1.5% agar), were used for phenotypic assay and the fungal culture for transcriptomic analysis.

2.4. Hyphal Biomass and Conidial Production

For hyphal biomass and conidial production experiments, 100 μL aliquots of a 107 conidia mL−1 suspension were evenly spread onto the cellophane-covered SDAY and CDA plates, followed incubation at 25 °C with a 12/12 light/dark cycle. For the hyphal biomass test, samples from each strain (WT, Δrad6, and Δrad6::rad6) were collected and dried to a constant weight for hyphal biomass analysis. For conidiation evaluation, three 5 mm fungal plugs were collected from each plate and then subjected to ultrasonic vibration in 1 mL of 0.02% Tween 80 liquid to quantify conidiation capability. Samples were collected on days 4, 5, 6, and 7, with three independent repeats.

2.5. Transcriptome Analysis and Verification

As described above, cultures of WT and Δrad6 strains were collected after incubating at 25 °C with a 12/12 light/dark cycle for four days. For transcriptomic analysis, three independent samples from each strain were submitted to Vazyme Biotech Co. (Nanjing, China). Briefly, the mRNA was extracted from each culture using magnetic oligo(dT) beads and fragmenting into little segments. The segments and random hexamer primers were used as templates to synthesize the first-strand cDNA, which were further used for the second-strand cDNA synthesis. The double-strand cDNA was end-repaired by adding singular adenines to the ends. After being mixed with proper adaptors, the cDNA library was sequenced on an Illumina HiSeq™ platform in Majorbio (Shanghai, China). To identify differentially expressed genes (DEGs) between WTs and rad6 mutant strains, the log2 R(Δrad6/WT) was calculated. DEGs were defined as downregulated (log2 R(Δrad6/WT) ≤ −1) or upregulated (log2 R(Δrad6/WT) ≥ 1) at a significance level of q (corrected p) < 0.05. All DEGs were further analyzed for GO and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment at http://vip.majorbio.com/ (accessed on 1 March 2024). The transcription level of key genes involved in the central development pathway was obtained from the date of SDAY and CDA transcriptomes (Supplemental File S2, Tables S4 and S5. The real-time quantitative (qRT-PCR) analyses were performed using SYBR Green assay for transcriptome data verification. The WT and the rad6 disruption strain were incubated under the same conditions described in Section 2.2, followed by RNA extraction. The one-step gDNA removal and DNA synthesis kit (TransGen Biotech, Beijing, China) was used for the reverse transcription and synthesis of the cDNA. Five downregulated and five upregulated genes from the transcriptome data were selected for qRT-PCR verification, and the primers were designed using the online tool Primer3 (https://bioinfo.ut.ee/primer3-0.4.0/) (accessed on 14 September 2023), followed by the primer specificity checked with NCBI Primer-Blast (https://www.ncbi.nlm.nih.gov/tools/primer-blast/, accessed on 14 September 2023) (Supplemental File S1, Table S2).

2.6. Phenotype Validation Experiments

For conidial production, 100 μL of conidial suspension at a concentration of 107 conidia mL−1 was uniformly inoculated into SDAY and CDAY Petri plates containing varying concentrations of biotin. Biotin was filtered and sterilized through a 0.22 μm sterile filter (Membrane: PES, Syringe filter, Tianjin, China) and then added to SDAY and CDA culture media that had been sterilized at high temperatures, with final concentrations of 1.5 µM, 3 µM, and 6 µM, respectively. Three culture samples (5 mm plugs) were collected from each dish, and ultrasonic vibration was applied to calculate conidia yield.

2.7. Statistical Analysis

Data collected from three independent experimental replicates were subjected to one-factor (strain) analysis of variance, followed by Tukey’s honestly significant difference (HSD) test, to assess the variance in means between the Δrad6 strain and its control strains.

3. Results

3.1. Bioinformatical Analysis of Rad6 Protein in B. bassiana

The study conducted a bioinformatic analysis of the Rad6 protein in B. bassiana. Specifically, in the B. bassiana, BBA_01469 was identified as a homologous gene to S. cerevisiae Rad6 (NP_011457), with an identity of 77.5% and an E value of 2 × 10−92, named Rad6 herein. rad6 is a gene with three introns and encodes a 151-amino acid protein (molecular mass: 17.2 kD; isoelectric point: 5.44). This protein includes a Ubiquitin-conjugating enzyme E2 domain (UBQ-conjugat_E2; residues 7–145) and two E3 interaction residues (residues 64–65; 98–99), separated by an active site cysteine at position 88 (Figure 1a). Notably, the typical domains and residues align precisely with Rad6 orthologs in S. cerevisiae, Aspergillus nidullans, and Aspergillus flavus. However, no nuclear localization signal was found, except for residues 9–19 in S. cerevisiae. In a comparative analysis of all fungal Rad6 homologous protein sequences, the phylogenetic tree revealed two branches: filamentous fungi and yeast. The identity of the filamentous fungal branch protein sequences exceeded 93%, while the yeast branch ranged from 77% to 84%. The filamentous fungal branch could be roughly subdivided into two smaller branches: pathogenic fungi (with a higher similarity ranging from 96% to 100%) and non-pathogenic Aspergillus and Penicillium (with a lower similarity ranging from 93% to 95%) (Figure 1b).

3.2. Vital Role of Rad6 Protein in Conidiation and Hyphal Growth

The Δrad6 strains exhibited significant defects in both conidial production and hyphal growth, irrespective of whether they were cultured in nutrient-rich or nutrient-poor media. When cultured on SDAY plates, the conidial production of Δrad6 was significantly dropped by approximately 94% after a 7-day culture (Figure 2a). The dry hyphal mass of the Δrad6 was also significantly lower than that of the control strains (Figure 2b). Similarly, on nutrient-poor CDA plates, the Δrad6 strains produced only 33.65% of the conidia that the WT strains did after four days of incubation and reached 33.78% after seven days (Figure 2c). The dry hyphal biomass of Rad6 disruption mutant ranged from 40.93% to 54.22% of that in WT strains (Figure 2d). Generally, the phenotypic defects caused by disruption of rad6 could be rescued by its complementation. The dry hyphal biomass and the conidial production of Δrad6::rad6 cultured on SDAY or CDA were restored to a level close to that of WT.

3.3. Transcriptomic View of Rad6-Mediated Regulation Pathway Differences

Transcriptomic analysis was carried out to reveal the gene expression patten differences between Δrad6 and WT, cultured on both the nutrient-rich medium SDAY and nutrient-poor medium CDA. The volcano plot exhibits a strikingly similar pattern in the disruption and WT strains incubated in both SDAY and CDA media. Specifically, 1976 downregulated and 2109 upregulated differentially expressed genes (DEGs) were identified in Δrad6 vs. WT strains incubated in SDAY, while 1437 and 1794 DEGs were found in CDA cultures (Figure 3a,b). KEEG enrichment analysis indicated that the disruption strains exhibited alterations in the DNA replication, tryptophan, and biotin metabolism pathways in both transcriptomes (Figure 3c,d).
Further investigation revealed the existence of 2006 total DEGs, comprising 939 upregulated and 813 downregulated DEGs, common to both transcriptomes (Figure 4a–c). Enrichment analysis of these common DEGs identified DNA replication and the biotin metabolism pathway as the top two significantly influenced pathways (Figure 4d). Interestingly, the transcriptional levels of genes flbC, brlA, abaA, and wetA were significantly downregulated in both nutrient-rich/poor transcriptomes (Supplemental File S1, Table S3). In B. bassiana, the function of flbC, brlA, abaA, and wetA was characterized as critical genes in the central developmental pathway (CDP), which affects the fungal asexual development [13,14,15]. These results show that rad6 may be involved in fungal conidiation by regulating the transcription of various genes in B. bassiana.

3.4. Conidiation Increases with Biotin Addition

To investigate the potential biotin metabolism disorder, conidia from Δrad6 strains and the control strains were inoculated onto SDAY and CDA plates containing various biotin concentrations and incubated for eight days under optimal conditions. On the SDAY plates, the conidial production of Δrad6 was 6.13 × 107 conidia per cm2, only 9.52% of that of WT. However, it showed a significant increase (p-value < 0.05) to 39.24 × 107 conidia per cm2 after the addition of 6 µM biotin (Figure 5a). The relative conidiation inhibition (RCI) rates due to the addition of biotin were significantly lower in Δrad6 than that in WT (Figure 5b). On CDA plates, the conidial production of Δrad6 was 0.89 × 107 conidia per cm2, 33.78% of that of the WT strains. However, it showed a significant increase (p-value < 0.05) to 1.03 × 107 conidia per cm2 after the addition of 6 µM biotin (Figure 5c). The relative conidiation inhibition (RCI) rates due to the addition of biotin were also significantly lower in Δrad6 strains than that in WT (Figure 5d).

3.5. Transcription Pattern of CDP Genes in Different Media with and without Exogenous Biotin

In the transcriptomic validation experiment, all ten genes that exhibited downregulation in the SDAY and CDA transcriptomics data were consistently and significantly downregulated in qRT-PCR. Similarly, all ten genes that showed upregulation in the transcriptomics data displayed significant upregulation in the qRT-PCR results (Supplemental File S1, Figure S1). A noteworthy observation emerged when assessing the expression levels of key genes involved in the central developmental pathway. Compared to strains cultivated on SDAY without the addition of biotin, those exposed to biotin exhibited a substantial increase in the expression of key genes associated with this central development pathway (Figure 6a,b). This finding strongly suggests that biotin enhances the expression of pivotal genes within the central development pathway, consequently bolstering conidial production. In the case of CDA, the transcription levels of brlA, abaA, wetA, and vosA were 51.76%, 62.27%, 54.60%, and 220.12% of that in WT (Figure 6c). When exogenous biotin was added, their transcription levels were 47.39%, 119.35%, 170.82%, and 286.60% of that in WT (Figure 6d). This observation suggests that biotin augments conidial production by stimulating the expression of abaA, wetA, and vosA. These results consistently support our prior analysis.

4. Discussion

Post-modification ubiquitination profoundly influences various aspects of cellular events in model organisms such as S. cerevisiae, Arabidopsis thaliana, and Homo sapiens [27,28], thereby exerting a significant influence on the fate of life across kingdoms. In this study, we investigated the function of Rad6 in B. bassiana, a widely used filamentous entomopathogenic fungus for biocontrol. Our experimental data demonstrated that the rad6 gene affects the conidial production capacity of B. bassiana under both nutrient-rich and nutrient-poor conditions. Interestingly, we observed that, in addition to being regulated by the central developmental pathway, the conidiation deficiency in Δrad6 can be partially alleviated through the exogenous addition of biotin, as discussed below.
The deletion of rad6 leads to a dramatic loss in conidial production in both SDAY and CDA incubated strains, which is consistent with the previous studies in S. cerevisiae and M. oryzae [17,29]. The deletion of Rad6 reduced about 67% of the conidiation potential but only 45% of the hyphal mass in poor CDA plates. Similar but slightly different, in rich SDAY media, the absence of Rad6 caused a 94% conidia yield decrease but 55% hyphal mass loss. The higher percentage of conidial production rate loss and lower hyphal mass rate loss of Δrad6 compared to that in WT indicated that the absence of Rad6 plays a more substantial role in inhibiting conidial production than hyphal growth.
In filamentous fungi Aspergillus and Penicillium, the strain conidial production is usually mediated by the central development pathway consisting of the activators of BrlA, AbaA, and WetA [30]. The deletion of BrlA and AbaA has no influence on hyphal growth but abolishes aerial conidiation in B. bassiana [13]. WetA contributes to nearly 98% of the conidiation capacities of B. bassiana [14]. The downregulated transcription level of CDP-related genes in Δrad6 grown on SDAY and CDA explained its dramatic conidiation decrease. VosA is a downstream gene of CDP, which can also contribute to 88% of the conidiation capacities in B. bassiana [14]. The significant downregulation of VosA in SDAY cultures and upregulation in CDA cultures may have caused the profound conidiation loss in the former. On the other hand, the more severe down-regulation of CDP genes in SDAY cultures could be another reason. Fluffy genes flbA–flbE are well-known activators of brlA in Aspergillus nidulans [30,31]. Among them, flbA significantly negatively regulated the hyphal biomass of B. bassiana [15]. Therefore, the up-/down-regulated transcriptional level of flbA in Δrad6 cultured on SDAY/CDA, respectively, may provide an important insight into the more severe hyphal mass defect observed in Δrad6 when cultured on SDAY compared to CDA.
The KEGG analysis based on the transcriptomic data of Δrad6 vs. WT revealed that biotin metabolism abnormality is one of the key reasons for conidial reproduction. Biotin, a water-soluble vitamin, is indispensable for the viability of a broad spectrum of organisms, from bacteria to eukaryotes [32]. In filamentous fungi, biotin enhances the colonization and proliferation capability of Candida albicans in the catheter lumen [33]; it is an essential growth factor for the pathogenic fungus Allescheria boydii [34]. Adding exogenous biotin to the minimal medium could not only increase the formation of wild-type conidias by nearly three-fold, but also alleviate the conidial formation defect in the bioB knockout strain of Alternaria alternata [8]. In this study, with the addition of 6 µM biotin, the conidia yield of Δrad6 significantly increased by 45.64% and 48.77% when cultured on SDAY and CDA media. These results verify the hypothesis that the exogenous addition of biotin benefits conidial production in B. bassiana, which is consistent with the results that biotin can restore the growth defects in Aspergillus oryzae [35]. Apart from being a covalently bound coenzyme for carboxylases, evidence has shown that biotin can be attached to histones [36] and plays a role in gene repression, heterochromatin structures, and repression of transposons [37,38]. We suggest that the conidial yield reduction is caused by downregulated CDP genes whose transcription level may be indirectly influenced by the biotin metabolism disorder. This hypothesis gains substantial support from the significant upregulation of CDP-related genes in Δrad6 cultured with biotin compared to that in Δrad6 cultured without biotin.
In summary, Rad6 profoundly influences the conidiation capacity of B. bassiana under nutrient-rich and nutrient-poor conditions. Transcriptomic analysis of cultures grown in these diverse nutritional environments indicates a pivotal role of biotin in modulating conidial production, a finding substantiated through biotin supplementation. Further exploration suggests that biotin may impact conidial generation via histone biotinylation, consequently affecting gene transcription within the central developmental pathway. While these discoveries underscore the significant involvement of Rad6 in fungal asexual development, the precise mechanisms by which Rad6 influences biotin metabolism through ubiquitination and how biotin suppresses the expression of CDP-related genes remain elusive, warranting further investigation in subsequent studies.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/jof10090613/s1. Supplementary Material File S1: Figure S1. RT verification of SDAY (a) and CPZ (b) boards. The primers utilized for dual identification are detailed in Table S2. Table S1. Paired primers designed for constructing deletion mutants and complementary mutants of rad6 in B. bassiana. Table S2. Paired primers used for assessing the reliability of transcriptome data in B. bassiana through qPCR analysis. Table S3. Transcript levels of Fluffy genes flbAflbE in SDAY and CDA transcriptomes. Supplementary Material Flie S2: Table S4. Ttranscriptome results of SDAY cultures (SDAY_Rad6/SDAY_WT). Table S5. Ttranscriptome results of CDA cultures (CDA_Rad6/CDA_WT).

Author Contributions

Y.G. (Yuhan Guo): Data curation; investigation; methodology; softwar; supervision; visualizatio; writing—original draf; writing—review and editing. H.H.: Data curation; validation; investigation; software; writing—review and editin. Y.G. (Yi Guan): Funding acquisition; data curatio; validatio; investigatio; methodolog; writing—review and editin. L.Z.: Funding acquisition; project administration; resource; software. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the National Natural Science Foundation of China (32072481) and the Fujian Provincial Natural Science Foundation General Program (2023J01262 and 2021J01605). This work was also supported by the Qishan Scholar Foundation from Fuzhou University (XRC-24067) and a visiting scholarship from the China Scholarship Council (CSC).

Institutional Review Board Statement

This article does not contain any studies with human participants or animals performed by any of the authors.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data generated or analyzed during this study are included in this published article (and its Supplementary Information Files). All engineered strains used in this study are available upon reasonable request. All data generated or analyzed during this study are deposited in figshare, DOI: https://doi.org/10.6084/m9.figshare.26156842. The transcriptome sequencing data were deposited in the NCBI database with the accession number PRJNA1037315.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wang, D.Y.; Mou, Y.N.; Du, X.; Guan, Y.; Feng, M.G. Ubr1-mediated ubiquitylation orchestrates asexual development, polar growth, and virulence-related cellular events in Beauveria bassiana. Appl. Microbiol. Biotechnol. 2021, 105, 2747–2758. [Google Scholar] [CrossRef]
  2. Ciechanover, A. The unravelling of the ubiquitin system. Nat. Rev. Mol. Cell Biol. 2015, 16, 322–324. [Google Scholar] [CrossRef]
  3. Ye, Y.H.; Rape, M. Building ubiquitin chains: E2 enzymes at work. Nat. Rev. Mol. Cell Biol. 2009, 10, 755–764. [Google Scholar] [CrossRef] [PubMed]
  4. Streit, W.R.; Entcheva, P. Biotin in microbes, the genes involved in its biosynthesis, its biochemical role and perspectives for biotechnological production. Appl. Microbiol. Biotechnol. 2003, 61, 21–31. [Google Scholar] [CrossRef]
  5. Zempleni, J. Uptake, localization, and noncarboxylase roles of biotin. Annu. Rev. Nutr. 2005, 25, 175–196. [Google Scholar] [CrossRef]
  6. Leon-Del-Rio, A. Biotin in metabolism, gene expression, and human disease. J. Inherit. Metab. Dis. 2019, 42, 647–654. [Google Scholar] [CrossRef] [PubMed]
  7. Wu, P.C.; Chen, C.W.; Choo, C.Y.L.; Chen, Y.K.; Yago, J.I.; Chung, K.R. Biotin biosynthesis affected by the NADPH oxidase and lipid metabolism is required for growth, sporulation and infectivity in the citrus fungal pathogen Alternaria alternata. Microbiol. Res. 2020, 241, 126566. [Google Scholar] [CrossRef]
  8. Chen, X.; Yan, M.; Xie, F.L.; Dai, J.; Li, D.S.; Wang, Z.; Li, Z.J.; Yao, J.; Li, X.; Zhao, X.J. Biotin enhances salt tolerance of Torulopsis mogii. Ann. Microbiol. 2015, 65, 393–398. [Google Scholar] [CrossRef]
  9. Pirner, H.M.; Stolz, J. Biotin sensing in Saccharomyces cerevisiae is mediated by a conserved DNA element and requires the activity of biotin-protein ligase. J. Biol. Chem. 2006, 281, 12381–12389. [Google Scholar] [CrossRef] [PubMed]
  10. Sprenger, M.; Hartung, T.S.; Allert, S.; Wisgott, S.; Niemiec, M.J.; Graf, K.; Jacobsen, I.D.; Kasper, L.; Hube, B. Fungal biotin homeostasis is essential for immune evasion after macrophage phagocytosis and virulence. Cell. Microbiol. 2020, 22, e13197. [Google Scholar] [CrossRef]
  11. Yu, J.H.; Mah, J.H.; Seo, J.A. Growth and developmental control in the model and pathogenic aspergilli. Eukaryot. Cell 2006, 5, 1577–1584. [Google Scholar] [CrossRef]
  12. Son, Y.E.; Yu, J.H.; Park, H.S. Regulators of the Asexual Life Cycle of Aspergillus nidulans. Cells 2023, 12, 1544. [Google Scholar] [CrossRef]
  13. Zhang, A.X.; Mouhoumed, A.Z.; Tong, S.M.; Ying, S.H.; Feng, M.G. BrlA and AbaA Govern Virulence-Required Dimorphic Switch, Conidiation, and Pathogenicity in a Fungal Insect Pathogen. mSystems 2019, 4, e00140-19. [Google Scholar] [CrossRef]
  14. Li, F.; Shi, H.Q.; Ying, S.H.; Feng, M.G. WetA and VosA are distinct regulators of conidiation capacity, conidial quality, and biological control potential of a fungal insect pathogen. Appl. Microbiol. Biotechnol. 2015, 99, 10069–10081. [Google Scholar] [CrossRef] [PubMed]
  15. Guo, C.T.; Luo, X.C.; Ying, S.H.; Feng, M.G. Differential Roles of Five Fluffy Genes (flbA-flbE) in the Lifecycle In Vitro and In Vivo of the Insect-Pathogenic Fungus Beauveria bassiana. J. Fungi 2022, 8, 334. [Google Scholar] [CrossRef]
  16. Jentsch, S.; McGrath, J.P.; Varshavsky, A. The yeast DNA repair gene RAD6 encodes a ubiquitin-conjugating enzyme. Nature 1987, 329, 131–134. [Google Scholar] [CrossRef]
  17. Prakash, L. The structure and function of RAD6 and RAD18 DNA repair genes of Saccharomyces cerevisiae. Genome 1989, 31, 597–600. [Google Scholar] [CrossRef] [PubMed]
  18. McDonough, M.; Sangan, P.; Gonda, D.K. Characterization of novel yeast RAD6 (UBC2) ubiquitin-conjugating enzyme mutants constructed by charge-to-alanine scanning mutagenesis. J. Bacteriol. 1995, 177, 580–585. [Google Scholar] [CrossRef]
  19. Game, J.C.; Chernikova, S.B. The role of RAD6 in recombinational repair, checkpoints and meiosis via histone modification. DNA Repair 2009, 8, 470–482. [Google Scholar] [CrossRef] [PubMed]
  20. Hou, J.; Zhang, H.; Ding, J.L.; Feng, M.G.; Ying, S.H. Transcriptomic investigation reveals a physiological mechanism for Beauveria bassiana to survive under linoleic acid stress. iScience 2023, 26, 106551. [Google Scholar] [CrossRef]
  21. Wang, D.Y.; Tong, S.M.; Guan, Y.; Ying, S.H.; Feng, M.G. The velvet protein VeA functions in asexual cycle, stress tolerance and transcriptional regulation of Beauveria bassiana. Fungal Genet. Biol. 2019, 127, 1–11. [Google Scholar] [CrossRef] [PubMed]
  22. Xiao, G.H.; Ying, S.H.; Zheng, P.; Wang, Z.L.; Zhang, S.W.; Xie, X.Q.; Shang, Y.F.; St Leger, R.J.; Zhao, G.P.; Wang, C.S.; et al. Genomic perspectives on the evolution of fungal entomopathogenicity in Beauveria bassiana. Sci. Rep. 2012, 2, 483. [Google Scholar] [CrossRef] [PubMed]
  23. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef]
  24. Fang, W.G.; Zhang, Y.J.; Yang, X.Y.; Zheng, X.L.; Duan, H.; Li, Y.; Pei, Y. Agrobacterium tumefaciens-mediated transformation of Beauveria bassiana using an herbicide resistance gene as a selection marker. J. Invertebr. Pathol. 2004, 85, 18–24. [Google Scholar] [CrossRef]
  25. Wang, D.Y.; Ren, K.; Tong, S.M.; Ying, S.H.; Feng, M.G. Pleiotropic effects of Ubi4, a polyubiquitin precursor required for ubiquitin accumulation, conidiation and pathogenicity of a fungal insect pathogen. Environ. Microbiol. 2020, 22, 2564–2580. [Google Scholar] [CrossRef]
  26. Guan, Y.; He, H.M.; Guo, Y.H.; Zhang, L.B. Essential roles of Rad6 in conidial property, stress tolerance, and pathogenicity of Beauveria bassiana. Virulence 2024, 15, 2362748. [Google Scholar] [CrossRef]
  27. Sridhar, V.V.; Kapoor, A.; Zhang, K.L.; Zhu, J.J.; Zhou, T.; Hasegawa, P.M.; Bressan, R.A.; Zhu, J.K. Control of DNA methylation and heterochromatic silencing by histone H2B deubiquitination. Nature 2007, 447, 735–738. [Google Scholar] [CrossRef]
  28. Giraldo, M.I.; Xia, H.J.; Aguilera-Aguirre, L.; Hage, A.; van Tol, S.; Shan, C.; Xie, X.P.; Sturdevant, G.L.; Robertson, S.J.; McNally, K.L.; et al. Envelope protein ubiquitination drives entry and pathogenesis of Zika virus. Nature 2020, 585, 414–419. [Google Scholar] [CrossRef]
  29. Shi, H.B.; Chen, G.Q.; Chen, Y.P.; Dong, B.; Lu, J.P.; Liu, X.H.; Lin, F.C. MoRad6-mediated ubiquitination pathways are essential for development and pathogenicity in Magnaporthe oryzae. Environ. Microbiol. 2016, 18, 4170–4187. [Google Scholar] [CrossRef] [PubMed]
  30. Park, H.S.; Yu, J.H. Genetic control of asexual sporulation in filamentous fungi. Curr. Opin. Microbiol. 2012, 15, 669–677. [Google Scholar] [CrossRef]
  31. Ojeda-López, M.; Chen, W.; Eagle, C.E.; Gutiérrez, G.; Jia, W.L.; Swilaiman, S.S.; Huang, Z.; Park, H.S.; Yu, J.H.; Cánovas, D.; et al. Evolution of asexual and sexual reproduction in the aspergilli. Stud. Mycol. 2018, 91, 37–59. [Google Scholar] [CrossRef]
  32. Rodionov, D.A.; Mironov, A.A.; Gelfand, M.S. Conservation of the biotin regulon and the BirA regulatory signal in Eubacteria and Archaea. Genome Res. 2002, 12, 1507–1516. [Google Scholar] [CrossRef]
  33. Ohara, H. Growth and Control of Catheter-related Bloodstream Infection Causing Bacteria in Nutrient Solutions. Yakugaku Zasshi J. Pharm. Soc. Jpn. 2022, 142, 303–315. [Google Scholar] [CrossRef]
  34. Villela, G.G.; Cury, A. Study on the vitamin nutrition of Allescheria boydii Shear. J. Bacteriol. 1950, 59, 1–11. [Google Scholar] [CrossRef] [PubMed]
  35. Tanabe, Y.; Maruyama, J.I.; Yamaoka, S.; Yahagi, D.; Matsuo, I.; Tsutsumi, N.; Kitamoto, K. Peroxisomes Are Involved in Biotin Biosynthesis in Aspergillus and Arabidopsis. J. Biol. Chem. 2011, 286, 30455–30461. [Google Scholar] [CrossRef]
  36. Zempleni, J.; Wijeratne, S.S.; Hassan, Y.I. Biotin. Biofactors 2009, 35, 36–46. [Google Scholar] [CrossRef] [PubMed]
  37. Camporeale, G.; Giordano, E.; Rendina, R.; Zempleni, J.; Eissenberg, J.C. Drosophila melanogaster holocarboxylase synthetase is a chromosomal protein required for normal histone biotinylation, gene transcription patterns, lifespan, and heat tolerance. J. Nutr. 2006, 136, 2735–2742. [Google Scholar] [CrossRef] [PubMed]
  38. Zempleni, J.; Chew, Y.C.; Bao, B.; Pestinger, V.; Wijeratne, S.S.K. Repression of Transposable Elements by Histone Biotinylation. J. Nutr. 2009, 139, 2389–2392. [Google Scholar] [CrossRef]
Figure 1. Structural analysis of B. bassiana Rad6 and its phylogenetic relationship with other fungal homologues. (a) Domain analysis of B. bassiana Rad6. (b) The phylogenetic relationship between B. basiana Rad6 and other fungal homologues. Bootstrap values obtained through 1000 replications are presented at the nodes for reference. The branch length in the scale is proportional to the genetic distance, as determined using the neighbor-joining method within the MEGA7 program, accessible at http://www.megasoftware.net/ (accessed on 8 October 2019) (note: the accession number of each protein in NCBI and its similarity to the corresponding protein of Beauveria bassiana Rad6 are shown in parentheses).
Figure 1. Structural analysis of B. bassiana Rad6 and its phylogenetic relationship with other fungal homologues. (a) Domain analysis of B. bassiana Rad6. (b) The phylogenetic relationship between B. basiana Rad6 and other fungal homologues. Bootstrap values obtained through 1000 replications are presented at the nodes for reference. The branch length in the scale is proportional to the genetic distance, as determined using the neighbor-joining method within the MEGA7 program, accessible at http://www.megasoftware.net/ (accessed on 8 October 2019) (note: the accession number of each protein in NCBI and its similarity to the corresponding protein of Beauveria bassiana Rad6 are shown in parentheses).
Jof 10 00613 g001
Figure 2. Impacts of Rad6 disruption on conidiation capability and hyphal growth of B. bassiana in different culture media. (a) Conidia yield of each strain cultured in nutrient-rich SDAY plates. (b) Dry hyphal biomass of strains collected from nutrient-rich SDAY plates. (c) Conidia yield of each strain cultured in nutrient-poor CDA plates. (d) Dry hyphal biomass of strains collected from nutrient-poor CDA plates. All experiments were initiated by spreading 100 μL of a 107 conidia mL−1 suspension. Asterisked bars in each group significantly differ from unmarked bars (Tukey’s HSD, p < 0.05). Error bar: SD from three replicates.
Figure 2. Impacts of Rad6 disruption on conidiation capability and hyphal growth of B. bassiana in different culture media. (a) Conidia yield of each strain cultured in nutrient-rich SDAY plates. (b) Dry hyphal biomass of strains collected from nutrient-rich SDAY plates. (c) Conidia yield of each strain cultured in nutrient-poor CDA plates. (d) Dry hyphal biomass of strains collected from nutrient-poor CDA plates. All experiments were initiated by spreading 100 μL of a 107 conidia mL−1 suspension. Asterisked bars in each group significantly differ from unmarked bars (Tukey’s HSD, p < 0.05). Error bar: SD from three replicates.
Jof 10 00613 g002
Figure 3. Transcriptome analysis of Δrad6 vs. WT collected from different culture media. (a,b) Distribution of log2 FC and p values for genes identified in the transcriptomes from SDAY and CDA cultures. (c,d) KEGG pathway enrichment analysis of DEGs from SDAY and CDA cultures. Differentially expressed genes (DEGs) are defined as log2 FC ≤ –1 or log2 FC ≥ 1 at the level of p < 0.05. The remaining genes are insignificantly affected (–1 ≤ log2 FC ≤ 1). (Note: DnR refers to downregulated, UpR refers to upregulated, NDR refers to non-significant differentially regulated. The pathways are ordered based on the significant results in the KEGG enrichment analysis. p-value < 0.001 labeled ***, p-value < 0.01 labeled **, and p-value < 0.05 labeled *).
Figure 3. Transcriptome analysis of Δrad6 vs. WT collected from different culture media. (a,b) Distribution of log2 FC and p values for genes identified in the transcriptomes from SDAY and CDA cultures. (c,d) KEGG pathway enrichment analysis of DEGs from SDAY and CDA cultures. Differentially expressed genes (DEGs) are defined as log2 FC ≤ –1 or log2 FC ≥ 1 at the level of p < 0.05. The remaining genes are insignificantly affected (–1 ≤ log2 FC ≤ 1). (Note: DnR refers to downregulated, UpR refers to upregulated, NDR refers to non-significant differentially regulated. The pathways are ordered based on the significant results in the KEGG enrichment analysis. p-value < 0.001 labeled ***, p-value < 0.01 labeled **, and p-value < 0.05 labeled *).
Jof 10 00613 g003
Figure 4. Analysis of genes co-occurring in the transcriptome of different media in B. bassiana. (ac) Overlapped in the total, upregulated, and downregulated DEGs from SDAY and CDA cultures in B. bassiana. (d) KEGG pathway enrichment analysis of common DEGs from SDAY and CDA cultures. The pathways are ordered based on the significant results in the KEGG enrichment analysis. p-value < 0.001 labeled ***, p-value < 0.01 labeled **, and p-value < 0.05 labeled *.
Figure 4. Analysis of genes co-occurring in the transcriptome of different media in B. bassiana. (ac) Overlapped in the total, upregulated, and downregulated DEGs from SDAY and CDA cultures in B. bassiana. (d) KEGG pathway enrichment analysis of common DEGs from SDAY and CDA cultures. The pathways are ordered based on the significant results in the KEGG enrichment analysis. p-value < 0.001 labeled ***, p-value < 0.01 labeled **, and p-value < 0.05 labeled *.
Jof 10 00613 g004
Figure 5. Phenotype verification of biotin addition. (a,c) Conidiation of each strain in SDAY and CDA plates amended with different concentrations of biotin. (b,d) Relative conidiation inhibition rate of biotin addition at different concentrations. Asterisked bars in each group significantly differ from unmarked bars (Tukey’s HSD, p < 0.05). Error bar: SD from three replicates.
Figure 5. Phenotype verification of biotin addition. (a,c) Conidiation of each strain in SDAY and CDA plates amended with different concentrations of biotin. (b,d) Relative conidiation inhibition rate of biotin addition at different concentrations. Asterisked bars in each group significantly differ from unmarked bars (Tukey’s HSD, p < 0.05). Error bar: SD from three replicates.
Jof 10 00613 g005
Figure 6. qRT-PCR verification with or without exogenous biotin. (a,c) Transcription levels of CDP-related genes in biotin-free SDAY and CDA plates. (b,d) Transcription levels of CDP-related genes in SDAY and CDA plates with biotin added. Asterisked bars in each group significantly differ from unmarked bars (Tukey’s HSD, p < 0.05). Error bar: SD from three replicates.
Figure 6. qRT-PCR verification with or without exogenous biotin. (a,c) Transcription levels of CDP-related genes in biotin-free SDAY and CDA plates. (b,d) Transcription levels of CDP-related genes in SDAY and CDA plates with biotin added. Asterisked bars in each group significantly differ from unmarked bars (Tukey’s HSD, p < 0.05). Error bar: SD from three replicates.
Jof 10 00613 g006
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Guo, Y.; He, H.; Guan, Y.; Zhang, L. Rad6 Regulates Conidiation by Affecting the Biotin Metabolism in Beauveria bassiana. J. Fungi 2024, 10, 613. https://doi.org/10.3390/jof10090613

AMA Style

Guo Y, He H, Guan Y, Zhang L. Rad6 Regulates Conidiation by Affecting the Biotin Metabolism in Beauveria bassiana. Journal of Fungi. 2024; 10(9):613. https://doi.org/10.3390/jof10090613

Chicago/Turabian Style

Guo, Yuhan, Haomin He, Yi Guan, and Longbin Zhang. 2024. "Rad6 Regulates Conidiation by Affecting the Biotin Metabolism in Beauveria bassiana" Journal of Fungi 10, no. 9: 613. https://doi.org/10.3390/jof10090613

APA Style

Guo, Y., He, H., Guan, Y., & Zhang, L. (2024). Rad6 Regulates Conidiation by Affecting the Biotin Metabolism in Beauveria bassiana. Journal of Fungi, 10(9), 613. https://doi.org/10.3390/jof10090613

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop