Next Article in Journal
Associations between Serum Kallistatin Levels and Markers of Glucose Homeostasis, Inflammation, and Lipoprotein Metabolism in Patients with Type 2 Diabetes and Nondiabetic Obesity
Previous Article in Journal
Modified mRNA-Mediated CCN5 Gene Transfer Ameliorates Cardiac Dysfunction and Fibrosis without Adverse Structural Remodeling
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comprehensive Insights into the Remarkable Function and Regulatory Mechanism of FluG during Asexual Development in Beauveria bassiana

Institute of Plant Protection and Microbiology, Zhejiang Academy of Agricultural Sciences, Hangzhou 310021, China
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2024, 25(11), 6261; https://doi.org/10.3390/ijms25116261
Submission received: 29 April 2024 / Revised: 24 May 2024 / Accepted: 31 May 2024 / Published: 6 June 2024
(This article belongs to the Section Molecular Microbiology)

Abstract

:
Asexual development is the main propagation and transmission mode of Beauveria bassiana and the basis of its pathogenicity. The regulation mechanism of conidiation and the key gene resources for utilization are key links to improving the conidia yield and quality of Beauveria bassiana. Their clarification may promote the industrialization of fungal pesticides. Here, we compared the regulation of morphology, resistance to external stress, virulence, and nutrient utilization capacity between the upstream developmental regulatory gene fluG and the key genes brlA, abaA, and wetA in the central growth and development pathway. The results showed that the ΔbrlA and ΔabaA mutants completely lost the capacity to conidiate and that the ΔwetA mutant had seriously reduced conidiation capacity. Although the deletion of fluG did not reduce the conidiation ability as much as deletions of brlA, abaA, and wetA, it significantly reduced the fungal response to external stress, virulence, and nutrient utilization, while the deletion of the three other genes had little effect. Via transcriptome analysis and screening the yeast nuclear system library, we found that the differentially expressed genes in the ΔfluG mutants were concentrated in the signaling pathways of ABC transporters, propionate metabolism, tryptophan metabolism, DNA replication, mismatch repair, and fatty acid metabolism. FluG directly acted on 40 proteins that were involved in various signaling pathways such as metabolism, oxidative stress, and cell homeostasis. The analysis indicated that the regulatory function of fluG was mainly involved in DNA replication, cell homeostasis, fungal growth and metabolism, and the response to external stress. Our results revealed the biological function of fluG in asexual development and the responses to several environmental stresses as well as its influence on the asexual development regulatory network in B. bassiana.

1. Introduction

Among entomopathogenic fungi, Beauveria bassiana is a main source of the broadest-spectrum fungal pesticide that usually undergoes an asexual cycle [1]. Its conidia are formulated as active ingredients of fungal pesticides, produced with high quality to tolerate outdoor stresses, and are of special importance for biological control. A host infection by B. bassiana is initiated by conidia attaching to the insect cuticle, where conidial germination leads to hyphal extension for penetration through the cuticular layers via the expression of a variety of hydrolytic enzymes and other factors [2,3,4,5]. After entry into the host hemocoel, the hyphal development turns into a process for forming plenty of blastospores that propagate rapidly until the host is mummified to death [5,6,7]. The blastospores bind to the integument and grow as hyphae that penetrate the host cuticle again to form aerial conidia to begin a new infection cycle [8,9,10]. There are two processes of dimorphic transition during infection: conidial germination leads to hyphal penetration after conidial adhesion, and then the penetrating hyphae from the host hemocoel turn into blastospores. Although the mechanism of the dimorphic switch is critical for the virulence and asexual cycle of B. bassiana, much remains obscure.
The asexual cycle of B. bassiana comprises two distinct phases, hyphal growth and conidiation, which are precisely timed and genetically programmed in response to stimulation by internal and external signals [11,12,13,14]. Conidiation is genetically controlled by a central developmental pathway (CDP) consisting of the developmental regulators BrlA, AbaA, and WetA, which mediate the development of conidiophores and the formation of conidia in Aspergillus [15,16]. BrlA is a key activator that initiates conidiation, followed by the sequential activation of AbaA and WetA, respectively, in the middle and late phases [17,18,19,20]. PdbrlA, PdabaA, and PdwetA of Penicillium digitatum were also confirmed to control distinct stages of conidiogenesis [21]. B. bassiana needs to survive independently in the environment before infecting fungi, so the utilization of environmental nutrients, response to external stress, and conidiation will directly affect the vitality of fungi and their ability to infect insects [1,2,3,4]. Therefore, it is also very important whether asexual development-related genes are involved in responding to external environmental stimuli. In B. bassiana, conidiation is completely abolished in the absence of brlA or abaA [13] and is nearly abolished in the absence of wetA [11], and virulence is greatly attenuated in mutants with these deletions. Asexual development is governed by BrlA, AbaA, and WetA in the central pathway of B. bassiana. However, the response to environmental stress has not been detected after the deletion of brlA or abaA in B. bassiana.
The upstream developmental activation (UDA) pathway containing the signal transducer genes fluG and flbA–flbE could initiate conidiophore development by activating the CDP genes in Aspergillus [15,16,22]. The FluG protein was first discovered in Aspergillus niger, and has since been identified in other Aspergillus and Penicillium fungi [22,23,24,25]. In Aspergillus niger [19], FluG mainly activates a series of activators in the UDA pathway to positively regulate the CDP genes (brlA, abaA and wetA) in a non-linear transcriptional regulatory cascade mode, thereby terminating the nutritional growth and initiating the asexual sporulation process. The deficiency of FluG can inhibit the sporulation and petiole, while causing defects in conidiation structure leading to the fluffy phenotype [25,26]. In Aspergillus fumigatus, the fluG knockout strain showed similar growth morphology and normal sporulation compared to the wild strain on solid culture medium, but sporulation was absent in liquid culture, indicating that FluG only plays a special regulatory role in the sporulation process of A. fumigatus under liquid conditions [27,28]. The absence of fluG also did not lead to a decrease in fluffy colony morphology and sporulation in Penicillium digitatum [21]. The above research results indicate that the regulatory mechanisms of FluG in the asexual sporulation process of filamentous fungi are different and have species specificity.
In particular, fluG is the key genetic switch that terminates the growth of hyphae and initiates conidiation by activating flbAflbE in different ways after receiving the stimulation signal [29,30,31,32,33,34,35,36,37]. Meanwhile, fluG could be activated by inhibiting sfgA, which is a key suppressor participating in hyphal growth [24,25]. Transcription of the flb and CDP genes was still activated in the absence of fluG, and the changes were correlated with the yield changes of the conidia, implying that fluG is a core, but not unique, player in the UDA pathway [14]. As previously reported, fluG disruption in B. bassiana not only affects the conidial yield but also attenuates cell integrity and virulence and increases cellular sensitivities to different stresses. In view of the important role of fluG in the growth and development of B. bassiana, we comprehensively compared the differences between the upstream developmental regulatory gene fluG and the key genes brlA, abaA, and wetA in the central growth and development pathway in morphology, conidiation, nutrient utilization, the external stress response, and virulence and summarized the particular function of fluG located upstream of growth and development. Despite an in-depth and comprehensive analysis of the conidiation and environmental stresses of B. bassiana, the mechanism and signaling pathway of fluG remain poorly understood. This study sought to characterize the different roles of the fluG and CDP genes in morphology, conidiation, nutrient utilization, the external stress response, and virulence, including the signaling pathways they participate in and the downstream proteins they directly activate, in order to deepen our understanding of the mechanism of the upstream developmental regulatory gene fluG in growth and development, resistance to environmental stress, and virulence.

2. Results

2.1. Comparison of fluG Disruption with brlA, abaA, and wetA Disruption of Phenotypic and Asexual Cycle

To explore the unique role of fluG (BBA_04942) in the growth and development of B. bassiana, disruption mutants and corresponding complementation mutants of fluG and the CDP-related genes brlA (BBA_07544), abaA (BBA_00300), and wetA (BBA_06126) were constructed, as described above. All mutants were verified through PCR, RT-PCR, and qRT-PCR with paired primers (Figure S1 in the Supplementary Materials).
Compared with the control (wild and complemented) strains, the ΔfluG mutant grew slightly slower but without an apparently fluffy phenotype under the optimal regime (Figure 1A). Although the colonies of the ΔbrlA, ΔabaA, and ΔwetA mutants had sizes similar to the control strains, they were more cottony and thicker, especially the colonies of ΔbrlA and ΔabaA with greater mycelial density. Meanwhile, ΔwetA not only exhibited an obvious fluffy phenotype, but also had grooves on the back and heavy pigment accumulation.
The quantification of conidiation over 7 days on SDAY plates demonstrated the different degrees of reduction in conidiation yield in each deletion mutant (Figure 1B). Compared with the control strains, aerial conidiation was totally abolished in the ΔbrlA and ΔabaA mutants during the 7 days, while the ΔfluG and ΔwetA mutants suffered conidial yield losses of 73.9% and 97.1% on day 3, 90.0–97.0% and 98.3–99.7% on days 4–6, and 84.7% and 97.9% on day 7, respectively. The blastospore production level of each strain was monitored daily during an 8-day incubation in SDB culture. The ΔbrlA and ΔabaA mutants completely lost the ability to produce blastospores. The ΔfluG mutant showed 78.5–86.4% reductions on days 3–7, and the reduction diminished to 65.7% on day 8. The ΔwetA mutant showed only a 15.1% reduction on day 3. Subsequently, the reduction increased to 55.9–65.7% on days 4–8.
Aerial conidiation and blastospore production were restored by each gene complementation. The ΔbrlA and ΔabaA mutants completely lost the ability to produce aerial conidia on a plate or blastospores in a bath. The ΔwetA and ΔfluG mutants suffered severe defects in conidiation, and the ΔwetA mutant showed slightly more defects than the ΔfluG mutant.

2.2. Different Effects of FluG, BrlA, AbaA, and WetA on Multiple Stress Responses and Virulence

The deletion of the fluG, brlA, abaA, and wetA genes had different effects on the sensitivity of B. bassiana to external stress during 7 days of colony growth on CZA plates supplemented with chemical stressors (Figure 2). Compared with the WT strain, the EC50 values of the ΔfluG mutant for oxidative stress induced by menadione or H2O2 were reduced by 23.9% or 40.4%, respectively (Figure 2A,B). The ΔwetA mutant was also significantly more sensitive to menadione or H2O2, and its EC50 values were reduced by 13.6% or 7.6%, respectively (Tukey’s HSD, p = 0). The deletion of brlA and abaA did not cause a loss of fungal resistance to oxidative stress. The tolerance to the cell wall interference stress of Congo red or SDS decreased significantly by 40.6% or 26.7% in the ΔfluG mutant and by 24.8% or 32.1% in the ΔwetA mutant, respectively (Figure 2C,D). The ΔbrlA and ΔabaA mutants also did not have decreased resistance to cell wall interference stress. The resistance to NaCl hypertonic stress or carbendazim decreased significantly, by 26.0% or 68.0% in the ΔfluG mutant and by 8.2% or 12.1% in the ΔwetA mutant, while the resistance of the ΔbrlA and ΔabaA mutants did not change significantly (Figure 2E,F).
Since the ΔbrlA and ΔabaA mutants produced neither aerial conidia nor submerged blastospores, a blastospore-removed hyphal suspension of each strain was applied for normal infection through cuticular penetration by immersing (Figure 2G). In contrast to the WT strain, the semi-lethal times of the ΔbrlA and ΔabaA mutants increased by 1.86 times and 2.17 times, while those of the ΔfluG and ΔwetA mutants increased by 69.6% and 48.4%, respectively.

2.3. Different Effects of FluG, BrlA, AbaA, and WetA on Nutrient Utilization

Compared with the WT strain, all deletion mutants showed different degrees of growth defects in the CZA or CZA-derived media. After 7 days of standard incubation on the CZA, CZA-C, CZA-N, or CZA-C-N media, the colony sizes of the ΔfluG strain were significantly reduced to 75.5%, 74.5%, 77.2%, or 75.5% of those of the WT strain, respectively, while the deletion of brlA and wetA had little effect on growth (Figure 3A). When cultured on the CZA, CZA-C, or CZA-N media, the colony areas of the ΔabaA mutant increased significantly to 1.3, 1.3, or 1.4 times those of the WT strain, respectively. When the sole nitrogen source of CZA was replaced with 0.3% NH4Cl, NaNO2, or NH4NO3, among the mutant strains the ΔfluG mutant had the most significant damage due to the utilization of the substituted nitrogen source, and the colony sizes were reduced to 87.8%, 74.5%, or 92.4% of those of the WT strain (Figure 3B). When replacing the sole carbon source with 3% olive oil, maltose, trehalose, glycerol, glucose, fructose, mannitol, sorbitol, lactose, acetate, or ethanol, the ΔfluG mutant still had the most significant damage due to the utilization of the replaced carbon source among the mutated strains, and the colony area decreased by 13.9% to 41.7% (Figure 3C). These data implied that compared with brlA, abaA, and wetA, significant growth defects were observed on the CZA and CZA-derived media after the deletion of fluG.

2.4. Regulatory Roles of fluG in Global Gene Expression

The regulatory roles of the upstream developmental regulatory gene fluG were examined by analyzing the transcriptomes of the ΔfluG and WT strains. Three replicates were derived from 3-day SDAY cultures in which conidiation was rapidly developing. Compared with the WT strain, the ΔfluG mutant had 1281 upregulated (~12.4% of the genome) and 1257 downregulated (~12.2% of the genome) genes (Figure 4A). The gene ontology (GO) analysis (Figure 4B) revealed that these differentially expressed genes (DEGs) were enriched for three GO terms that were involved in biological processes, cellular components, and molecular function. There were 20 GO terms for biological processes, and among them the DEGs were mainly focused on cellular processes, metabolic processes, localization, and responses to stimuli. For eight enriched cellular component terms, the DEGs mainly participated in organelles, membranes, and extracellular regions. For nine GO terms related to molecular function, the DEGs were mainly enriched for binding, catalytic activity, and transporter activity.
In the ΔfluG mutant, the DEGs for the top 20 significantly enriched terms in the KEGG revealed that the upregulated genes were focused on ABC transporters, propanoate metabolism, tryptophan metabolism, and galactose metabolism, while the downregulated genes were concentrated in DNA replication, mismatch repair, and fatty acid biosynthesis (Figure 4C). The gene expression information of the ABC transporter, DNA replication, and fatty acid metabolism pathway on which the DEGs were significantly enriched and five other key signaling pathways were analyzed (Figure 4D). The results in Figure 5 revealed that the transcript levels of 22 DEGs involved in DNA replication and 8 DEGs involved in fatty acid biosynthesis were all downregulated, while the transcript levels of 17 DEGs affecting ABC transporters were all upregulated. All four DEGs involved in oxidative phosphorylation were downregulated, and all five DEGs involved in autophagy were upregulated. For the ubiquinone term, the transcript levels of 3-dehydroshikimate dehydratase (BBA_01589) and 4-hydroxyphenylpyruvate dioxygenase (BBA_08551) were upregulated to 5.92 and 4.23 times the normal levels, while 4-hydroxybenzoate polyprenyltransferase (BBA_09001) and flavin prenyltransferase (BBA_04668) were downregulated to 0.031 and 0.155 times the normal levels, which affected the normal development and metabolism of the fungi [38,39,40,41]. Four DEGs were upregulated and two were downregulated in the peroxisome pathway, which affected the ability of the fungi to resist oxidative stress. The age–rage signaling pathway is very important for maintaining cell carbohydrate and protein homeostasis. In this pathway, the expression level of NADPH (BBA_07926) was downregulated to 0.325 times the normal level, and two 1-phosphatidylinositol 4,5-diphosphate phosphodiesterases (BBA_03011 and BBA_06798) were upregulated by two-fold [42,43].
In order to verify the effectiveness of the transcriptome sequencing, six upregulated genes and six downregulated genes of B. bassiana were randomly selected for qRT-PCR analysis (Figure S2). The results showed that the relative expression trend of the 12 genes was consistent with the results of the transcriptome sequencing, which verified the validity of the results of the transcriptome analysis.

2.5. Screening and Analysis of Interacting Proteins

To investigate FluG-interacting proteins, we performed a yeast two-hybrid (Y2H) screening assay. In total, 53 positive clones were screened (Figure S3). Through library screening, positive clone sequence analysis, and comparison online, 40 proteins were directly affected by FluG in the 53 positive clones. The relevant information shown in Table 1 revealed that these proteins participated in various signaling pathways such as metabolism, oxidative stress, and cellular homeostasis. Through GO enrichment analysis, the results showed that the 40 proteins were mainly concentrated in metabolic and cellular processes in the biological process term; in cells, cell components, and organelle components in the cellular component term; and in binding and catalytic activities in the molecular function term (Figure 6A). Among them, 26 proteins were significantly enriched in six signaling pathways, according to the KEGG enrichment analysis. These pathways were translation, amino acid metabolism, energy metabolism, carbohydrate metabolism, the metabolism of cofactors and vitamins, and the biosynthesis of other secondary metabolites (Figure 6B). The library screening was verified by yeast two-hybridization between five randomly selected interacting proteins and FluG (Figure S4).

3. Discussion

In the upstream developmental activation pathway of aspergilli, fluG is required for the commencement of conidiation. It collaborates with different flb genes [15,16,44,45], then three sequentially active genes (brlA, abaA, and wetA) in the central developmental pathway are activated to mediate the development of conidiophores and conidia [17,20]. In B. bassiana, disruption of brlA and abaA, as described before, led to abolished conidiation and blastospore production and greatly attenuated virulence but had no negative impact on hyphal growth in various media [13]. Deletion of wetA and fluG led to reductions in conidiation of varying degrees, attenuated conidial virulence, and several defects in response to nutritional and abiotic stresses [11,14]. Since the responses of brlA and abaA deletion mutants to environmental stress had not been detected and the functions of genes located in the central developmental pathway and upstream developmental activation pathway had not been compared together, we compared the differences in morphology, conidiation, nutrient utilization, stress resistance, and virulence among mutants of these four genes (fluG, brlA, abaA, and wetA) and conducted in-depth research on the downstream regulatory mechanisms of fluG, as discussed below.
In the comparison of morphology, conidiation, and blastospore production, it was found that the genes in the CDP had stronger control over conidiation and blastospore production, and the mutants had a more cottony and thicker phenotype in B. bassiana. Disruption of fluG had a limited impact on conidiation and blastospore production and had little effect on radial growth. These results indicate that brlA and abaA in the CDP served as master regulators of asexual development, while the ability of wetA, located downstream, to regulate asexual development was stronger than that of fluG but weaker than that of brlA and abaA. Both brlA and abaA have highly conserved roles in the regulation of asexual development in other filamentous fungi [20,26,46,47,48]. WetA is also a crucial regulator of conidiation capacity in B. bassiana but exerts negative feedback control over conidiation in a way very different from Aspergillus fumigatus and Fusarium graminearum [11,20,49]. fluG has a strong ability to control conidiation, but it is always inferior to the regulatory genes located in the CDP. This indicates that the fluG-mediated regulatory pathway is not the only way that the UDA pathway regulates asexual development.
Apart from the severe defects in conidiation, the conidia of the ΔwetA and ΔfluG mutants led to different degrees of cell wall damage, slower germination, rapid viability loss, attenuated virulence, and reduced stress tolerance in B. bassiana [11,14]. Due to the conidiation defects of the ΔbrlA and ΔabaA mutants, experiments on conidial viability could not be carried out. In our study, we compared the phenotypic differences of the ΔbrlA, ΔabaA, ΔwetA, and ΔfluG mutants in nutrient utilization, environmental stress, and virulence of hyphae to find the differences between the UDA and CDP genes. The results confirm that brlA and abaA do not participate in nutrient utilization or responses to environmental stress, but wetA and fluG play important roles, especially fluG. The results indicate that compared with the genes in the CDP, fluG is more involved in sensing external stimuli. The fungal virulence was largely attenuated via hyphal infection in the ΔbrlA and ΔabaA mutants, and it was already confirmed that the difficulty of hyphal infection of an insect through cuticular penetration could be attributed to the markedly reduced transcripts of multiple genes critical for host adhesion and cuticle degradation [14]. In addition, the ΔfluG mutant was much more virulent than the ΔwetA mutant, more sensitive to environmental stress, and less capable of utilizing environmental carbon and nitrogen sources for growth. Different defects in nutrient utilization, responses to environmental stress, and hyphal virulence between the four genes indicate that fluG plays an important role in the fungus sensing stimuli from the host insect and the environment.
Here, a simple genetic model diagram for upstream and central regulators was drawn based on previous research on conidiation in B. bassiana (Figure 7). brlA, abaA and wetA, as key regulators in CDP, directly affect conidiation [11,13]. In the flb genes, flbA showed a much greater role than flbC in fungal conidiation, blastospore production and insect–pathogenic lifecycle. flbB, flbD and flbE were not influential on the expression of brlA or abaA irrespective of limited or little contribution to conidiation [50]. Therefore, the interaction relationship between flb genes and the mode of action on regulators in CDP are not yet clear, and further research is still needed. In a previous study, the relative transcriptional levels of all flb and CDP genes were consistently active in the ΔfluG mutant during continuous cultivation and were correlated with the yield changes of conidia and blastospores in B. bassiana, but were very different from those of other filamentous fungi [14,26,37,51,52]. This implied that fluG did not orchestrate the flb genes for the activation of brlA to facilitate conidiation.
Hence, transcriptome analysis and yeast library screening were used to further verify the signaling pathways and interacting proteins directly affected by FluG. In the absence of fluG, DNA replication, fatty acid biosynthesis, and oxidative phosphorylation were adversely affected, while ABC transporters and autophagy were promoted. In addition, some genes regulating normal development and metabolism, the response to oxidative stress, and cell carbohydrate and protein homeostasis were affected to varying degrees. Moreover, 40 proteins directly affected by FluG also mainly participated in metabolism, oxidative stress, and cellular homeostasis, which was consistent with the effects on their transcription levels. Although two-hybrid systems offer numerous advantages for the identification of novel protein–protein interactions and have been widely applied in the research of protein function, there are still some uncertainties [53,54]. Therefore, they provide us with reliable information on proteins interacting with FluG, but still require extensive research to determine the interaction patterns between FluG and downstream proteins.
In addition, we also found that the flb and CDP genes did not exist in the FluG-interacting proteins. All phenotypic and transcriptional changes in the ΔfluG mutants and proteins that interact with FluG imply that fluG did not directly interact with the flb genes or CDP genes to initiate conidiation. We speculated that there may be some factors regulating the flb genes or CDP genes in the proteins interacting with FluG, or there may be some other fluG-like gene(s) acting as a regulator(s) in the UDA pathway.

4. Materials and Methods

4.1. Microbial Cultivation

The wild-type ARSEF 2860 strain of B. bassiana and its mutants were grown in SDAY (4% glucose, 1% peptone, 1% yeast extract, and 1.5% agar) at 25 °C with a 12:12 h light/dark cycle for hyphal growth and conidiation. The stress responses were assayed in CZA (3% sucrose, 0.3% NaNO3, 0.1% K2HPO3, 0.05% KCl, 0.05% MgSO4, 0.001% FeSO4, and 1.5% agar) as a control and in CZA with different stresses. Escherichia coli DH5α and Top 10 (Invitrogen, Shanghai, China) were cultured in LB medium plus 100 μg mL−1 kanamycin or 50 μg mL−1 ampicillin, depending upon the resistance marker, at 37 °C for plasmid cloning and propagation. Agrobacterium tumefaciens AGL-1 for fungal transformation was cultivated at 28 °C in YEB medium [11].

4.2. Generation and Identification of fluG, brlA, abaA, and wetA Mutants

The genes fluG, brlA, abaA, and wetA (tag loci: BBA_04942, BBA_07544, BBA_00300, and BBA_06126) were deleted from the WT strain via homologous recombination of their 5′ (up) and 3′ (down) coding fragments separated by the bar marker in the vector p0380-5′x-bar-3′x (x = fluG, brlA, abaA, or wetA). For targeted gene complementation, a full-length coding sequence with flanking regions was inserted into the vector p0380-sur-x. All the 5′ and 3′ fragments and full-length sequences were cloned from the WT strain with paired primers up-F/R, down-F/R and fl-F/R, and digested with appropriate restriction enzyme sites (Table S1 in the Supplementary Materials). All the knockout plasmids and complement plasmids for fluG and wetA were integrated into corresponding strains via Agrobacterium-mediated transformation [11,14]. Since the ΔbrlA and ΔabaA mutants were completely unable to produce conidia, their resultant complement plasmids were ectopically integrated into protoplasts of each deletion mutant via polyethylene glycol-mediated transformation [13]. Putative deletion or complementary mutants were screened in terms of Bar resistance to phosphinothricin (200 μg mL−1) or Sur resistance to chlorimuron ethyl (10 μg mL−1), and then identified via PCR accompanied by sequencing with paired primers id-F/R.
The temporal transcript patterns of the genes in the wild-type, deletion, and complementary mutants were assessed during 3 days of growth at 25 °C with a 12:12 h light/dark cycle on cellophane overlaid on SDAY plates. Hyphal suspensions of each strain (100 μL aliquots) were spread on each plate to initiate the cultures. Total RNA was extracted separately using TRIzolTM Plus Reagent (Takara, Kusatsu, Japan) and treated with DNase I (Takara, Kusatsu, Japan) following the manufacturer’s instructions. Every 5 μg RNA sample was reverse-transcribed with a PrimeScriptTM RT reagent kit (Takara, Kusatsu, Japan). The temporal transcript patterns and relative transcript levels of the genes (18S rRNA was used as an internal standard) were assessed via triplicate RT-PCR and qRT-PCR assays with the primers listed in Table S1. Positive deletion mutants were analyzed in parallel with the WT strain and rescued mutants as control strains in the following experiments with three independent replicates.

4.3. Assays for Radial Growth, Conidiation, Blastospore Production, Hyphal Stress Responses, and Virulence

The conidia yield on SDAY were quantified as described previously [11]. Briefly, 100 μL aliquots of hyphal suspensions for each strain were spread on SDAY plates and incubated at 25 °C with a 12:12 h light/dark cycle for 7 days. From the third day of incubation onward, the conidia on the 5 mm diameter disks were washed daily into 1 mL of 0.02% Tween-80. After supersonic vibration, three samples were measured with a hemocytometer to assess the conidial concentration converted to the number of conidia cm−2 colony.
To examine the blastospore production, 100 μL aliquots of hyphal suspensions for each strain were inoculated in SDB broth and cultured for 7 days at 25 °C at 150 rpm. From the second day, the blastospore concentration in the broth was quantified via microscopy, and the blastospore yield was indicated as the number of cells per milliliter of broth.
Hyphal suspensions of each strain (100 μL aliquots) were spread on SDAY plates and incubated at 25 °C with a 12:12 h light/dark cycle. Hyphal mass plugs (5 mm diameter) were bored from cultures cultivated for 3 days and attached centrally to plates of SDAY; CZA alone (as control); and amended CZA supplemented with menadione (1~4 mM) or H2O2 (15~60 mM) for oxidative stress, NaCl (0.5~2 M) for osmotic stress, and Congo red (0.5~2 mg mL−1) for cell wall stress. After 10 days at 25 °C with a 12:12 h light/dark cycle, the diameters of all colonies were measured as indices of their radial growth rates using two measurements taken perpendicular to each other across their centers. Typical colonies of each strain were photographed. The ratio of the colony size under a given stress compared to that in the control condition was defined as the survival index (Is). For each of the tested strains, Is = 1/[1 + exp(a + bx)], where x is the concentration (C) of each stressful chemical. When Is = 0.5, the fitted equations gave the solutions (−a/b) for the effective concentrations for the stressful chemicals to suppress colony growth by 50% (EC50) [11,55].
The virulence of each strain was detected via hyphae as described previously [13]. Since the ΔbrlA and ΔabaA strains produced neither conidia nor blastospores, the fresh hyphae without blastospores in 3-day-old SDB cultures of each strain were suspended in 0.02% Tween-80 and standardized to a concentration of 10 mg mL−1. Three replicates of 40 larvae were separately immersed for 15 s in 40 mL aliquots of each strain. The same volume of 0.02% Tween-80 was used as a control. All treated samples were maintained at 25 °C with a 12:12 h light/dark cycle, and the survival records were monitored every 24 h until the records no longer changed for two consecutive days. LT50 (no. of days) was generated as an index of the fungal virulence via probit analysis of each time–mortality trend.

4.4. Assays for Relative Growth Areas on Different Media

Hyphal mass plugs (5 mm diameter) were bored from the cultures cultivated for 3 days on SDAY plates at 25 °C that were described above and attached centrally to the plates of CZA (3% sucrose, 0.3% NaNO3, 0.1% K2HPO4, 0.05% KCl, 0.05% MgSO4, and 0.001% FeSO4, plus 1.5% agar) and CZA-derived media with altered carbon/nitrogen sources. The CZA-derived media were prepared by removing the 3% sucrose, 0.3% NaNO3, or both from the CZA; replacing the sole carbon source with 3% olive oil, maltose, trehalose, glycerol, glucose, fructose, mannitol, sorbitol, lactose, acetate, or ethanol; and replacing the sole nitrogen source with 0.3% NH4Cl, NaNO2, or NH4NO3. After 7 days of incubation at 25 °C with a 12:12 h light/dark cycle, the diameter of each colony was cross-measured to calculate the growth area. The relative growth area was determined as the ratio of the mutant strain’s growth area to the wild strain’s growth area.

4.5. RNA Extraction, cDNA Library Construction, and RNA-Seq

The total RNA of the wild strain and the ΔfluG mutant was extracted from the cultures on SDAY plates after 3 days according to method 4.2. The total high-quality RNA isolated from three independent biological replicates of the samples was used to construct independent cDNA libraries.
The RNA-Seq of two samples, namely the wild strain (control) group and the ΔfluG group, each with three biological replicates, was performed on an Illumina Hiseq 2000 platform, from which 150 bp paired-end reads were generated. Raw sequences were deposited in the NCBI Short Read Archive (SRA) database (http://www.ncbi.nlm.nih.gov/Traces/sra/, accessed on 21 February 2024). The accession number of the RNA-Seq data was PRJNA1078435. Raw reads in the FASTQ format were first filtered by removing any reads containing adapter sequences and low-quality reads. At the same time, the Q20, Q30, GC content, and sequence duplication level of the clean data were each calculated. The cleaned reads were mapped onto the B. bassiana reference genome [56] using the Bowtie2 V2.5.2 software tool [57].

4.6. Transcriptome Annotation, Expression Profiling, Data Analysis, and Data Validation

Gene function was annotated based on homology searches in the NCBI non-redundant protein (Nr), NCBI nucleotide (Nt), Swiss-Prot protein, euKaryotic Orthologous Groups (KOG), Kyoto Encyclopedia of Genes and Genomes (KEGG), Gene Ontology (GO), and Protein family (Pfam) databases.
The expression levels of genes were measured using counts of reads normalized based on their respective lengths in the Cufflinks 2.0.2 package using its default settings (http://cole-trap-nell-lab.github.io/cufflinks/, accessed on 2 November 2023) for normalization (genometric), followed by a distribution analysis in terms of Fragments Per Kilobase of exon per Million mapped reads (FPKM) units. The differentially expressed genes (DEGs) seq package was used to identify DEGs between the control and gene mutant groups; it provided statistical routines for determining differential expression levels in the digital gene expression data by applying a model based on the negative binomial distribution. The resulting p-values were adjusted using Benjamini and Hochberg’s approach to control the false discovery rate (FDR). Genes with an adjusted p-value < 0.05, as detected by DEGSeq, were designated as differentially expressed. The fold-change in a given gene’s expression between samples was calculated as log2 (treatment FPKM value/control FPKM value).
The DEG-related signaling pathways were analyzed in B. bassiana using the KEGG [58,59,60]. Related maps were also obtained from the KEGG.
To confirm the results of the transcriptome comparisons, qRT-PCR was performed on 16 randomly selected genes in B. bassiana [61]. The samples were consistent with those used by the transcriptome. The primers for qRT-PCR are listed in Supplementary Table S2.

4.7. Nuclear cDNA Library Construction and Yeast Two-Hybrid Library Screening

For nuclear library construction, RNA was isolated from the wild strain after inoculation on SDAY for 3 days. A nuclear library of 1.12 × 107 mL−1 clones was constructed with the vector pGADT7 (Takara Bio, Japan) according to the CloneMiner instructions by OEBiotech (Shanghai, China). For library screening, an open reading frame with fluG as the bait construct was generated with the vector pGBKT7 (Takara Bio, Japan). Yeast two-hybrid (Y2H) screening was performed using a Matchmaker Two-Hybrid system (Clontech). Briefly, the bait and cDNA library plasmids were co-transformed into Y2HGold competent yeast cells. After transformation, the co-transformants were plated onto SD/-Trp/-Leu/X-α-Gal/AbA (DDO/X/A) agar plates. Single blue colonies were selected and plated onto higher-stringency SD/-Trp/-Leu/-His/-Ade/X-α-Gal/AbA (QDO/X/A) plates to test reporter gene expression. Positive colonies were thereafter re-seeded in SD/-Trp/-Leu/-His/-Ade liquid media. Prey plasmids were extracted from putatively positive clones using an Easy Yeast Plasmid Isolation Kit (Clontech, Mountain View, CA, USA) and sequenced. After sequencing with the primer T7, the B. bassiana genes corresponding to the inserts in these clones were identified via Blast searches and then analyzed using the KEGG database.
To further confirm the interactions, the full-length cDNA sequences of related genes were amplified and inserted into the pGADT7 vector, which was co-transformed with pGBKT7 bait plasmids into yeast strain Y2HGold and co-cultured on DDO/X/A and QDO/X/A plates to test for interactions. The same positive and negative controls used in the Y2H screening were included in these experiments.

4.8. Statistical Analyses

All samples from three repeated assays were quantitative indices for phenotypic changes among the tested strains in a one-way analysis of variance. Significant differences among the wild, gene disruption, and complementation mutant strains were determined with Tukey’s honest significance test (Tukey’s HSD).

5. Conclusions

After comparing the regulation of morphology, resistance to external stress, virulence, and nutrient utilization capacity between the upstream developmental regulatory gene fluG and the key genes brlA, abaA, and wetA in the central growth and development pathway in Beauveria bassiana, we found that fluG serves as an important regulator in the UDA pathway and has less ability to control asexual development than genes in the CDP, but can strongly sense stimuli from host insects and the environment. FluG mainly affects fungal metabolism, oxidative stress, and cellular homeostasis, but does not directly interact with the flb genes or CDP genes to initiate conidiation. These results led to speculation that warrants further research on the proteins interacting with FluG and other fluG-like genes existing in the UDA pathway of asexual development.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/ijms25116261/s1.

Author Contributions

Conceptualization, J.C., J.Z. and F.L.; methodology, F.L.; validation, F.L., H.Z. and K.Y.; formal analysis, F.L.; investigation, J.Z. and F.L.; data curation, H.Z. and K.Y.; writing—original draft preparation, F.L.; writing—review and editing, J.C., J.Z. and F.L.; visualization, F.L.; supervision, J.C. and J.Z.; project administration, J.C., J.Z. and F.L.; funding acquisition, J.Z. and F.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by the National Natural Science Foundation of China (32001963), the Zhejiang Provincial Natural Science Foundation of China (LY23C140002), the China Agriculture Research System of MOF and MARA (CARS-24-G-07), and the “Pioneer” and “Leading Goose” R&D Program of Zhejiang (No. 2023C02030).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data generated or analyzed during this study are included in the published paper and associated Supplementary Materials. Raw sequences were deposited in the NCBI Short Read Archive (SRA) database (http://www.ncbi.nlm.nih.gov/Traces/sra/, accessed on 21 February 2024). The accession number of RNA-seq data was PRJNA1078435.

Acknowledgments

Thanks to Shanghai OE Biotech Co., Ltd. (Shanghai, China) for help with RNA-seq data sequencing and analysis, nuclear cDNA library construction, and Y2H library screening.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Wang, C.S.; Feng, M.G. Advances in fundamental and applied studies in China of fungal biocontrol agents for use against arthropod pests. Biol. Control 2014, 68, 129–135. [Google Scholar] [CrossRef]
  2. Holder, D.J.; Keyhani, N.O. Adhesion of the entomopathogenic fungus Beauveria (Cordyceps) bassiana to substrata. Appl. Environ. Microbiol. 2005, 71, 5260–5266. [Google Scholar] [CrossRef]
  3. Ye, S.D.; Ying, S.H.; Chen, C.; Feng, M.G. New solid-state fermentation chamber for bulk production of aerial conidia of fungal biocontrol agents on rice. Biotechnol. Lett. 2006, 28, 799–804. [Google Scholar] [CrossRef] [PubMed]
  4. Holder, D.J.; Kirkland, B.H.; Lewis, M.W.; Keyhani, N.O. Surface characteristics of the entomopathogenic fungus Beauveria (Cordyceps) bassiana. Microbiology 2007, 153, 3448–3457. [Google Scholar] [CrossRef] [PubMed]
  5. Pedrini, N.; Ortiz-Urquiza, A.; Huarte-Bonnet, C.; Zhang, S.; Keyhani, N.O. Targeting of insect epicuticular lipids by the entomopathogenic fungus Beauveria bassiana: Hydrocarbon oxidation within the context of a host-pathogen interaction. Front. Microbiol. 2013, 4, 24. [Google Scholar] [CrossRef] [PubMed]
  6. Fang, W.G.; Leng, B.; Xiao, Y.H.; Jin, K.; Ma, J.C.; Fan, Y.H.; Feng, J.; Yang, X.Y.; Zhang, Y.J.; Pei, Y. Cloning of Beauveria bassiana chitinase gene Bbchit1 and its application to improve fungal strain virulence. Appl. Environ. Microbiol. 2005, 71, 363–370. [Google Scholar] [CrossRef] [PubMed]
  7. Zhang, Y.J.; Feng, M.G.; Fan, Y.H.; Luo, Z.B.; Yang, X.Y.; Wu, D.; Pei, Y. A cuticle-degrading protease (CDEP-1) of Beauveria bassiana enhances virulence. Biocontrol. Sci. Technol. 2008, 18, 551–563. [Google Scholar] [CrossRef]
  8. Wanchoo, A.; Lewis, M.W.; Keyhani, N.O. Lectin mapping reveals stage-specific display of surface carbohydrates in in vitro and haemolymph derived cells of the entomopathogenic fungus Beauveria bassiana. Microbiology 2009, 155, 3121–3133. [Google Scholar] [CrossRef] [PubMed]
  9. Wang, J.; Ying, S.H.; Hu, Y.; Feng, M.G. Mas5, a homologue of bacterial DnaJ, is indispensable for the host infection and environmental adaptation of a filamentous fungal insect pathogen. Environ. Microbiol. 2016, 18, 1037–1047. [Google Scholar] [CrossRef] [PubMed]
  10. He, P.H.; Dong, W.X.; Chu, X.L.; Feng, M.G.; Ying, S.H. The cellular proteome is affected by a gelsolin (BbGEL1) during morphological transitions in aerobic surface versus liquid growth in the entomopathogenic fungus Beauveria bassiana. Environ. Microbiol. 2016, 18, 4153–4169. [Google Scholar] [CrossRef]
  11. 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]
  12. 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]
  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. [Google Scholar] [CrossRef] [PubMed]
  14. Guo, C.T.; Peng, H.; Tong, S.M.; Ying, S.H.; Feng, M.G. Distinctive role of fluG in the adaptation of Beauveria bassiana to insect-pathogenic lifecycle and environmental stresses. Environ. Microbiol. 2021, 23, 5184–5199. [Google Scholar] [CrossRef]
  15. Etxebeste, O.; Garzia, A.; Espeso, E.A.; Ugalde, U. Aspergillus nidulans asexual development: Making the most of cellular modules. Trends. Microbiol. 2010, 18, 569–576. [Google Scholar] [CrossRef] [PubMed]
  16. Park, H.S.; Yu, J.H. Genetic control of asexual sporulation in filamentous fungi. Curr. Opin. Microbiol. 2012, 15, 669–677. [Google Scholar] [CrossRef] [PubMed]
  17. Adams, T.H.; Boylan, M.T.; Timberlake, W.E. brlA is necessary and sufficient to direct conidiophore development in Aspergillus nidulans. Cell 1988, 54, 353–362. [Google Scholar] [CrossRef] [PubMed]
  18. Chang, Y.C.; Timberlake, W.E. Identification of Aspergillus brlA response elements (BREs) by genetic selection in yeast. Genetics 1993, 133, 29–38. [Google Scholar] [CrossRef] [PubMed]
  19. Sewall, T.C.; Mims, C.W.; Timberlake, W.E. abaA controls phialide differentiation in Aspergillus nidulans. Plant. Cell. 1990, 2, 731–739. [Google Scholar] [CrossRef]
  20. Tao, L.; Yu, J.H. AbaA and WetA govern distinct stages of Aspergillus fumigatus development. Microbiology 2011, 157, 313–326. [Google Scholar] [CrossRef]
  21. Wang, M.S.; Sun, X.P.; Zhu, C.Y.; Xu, Q.; Ruan, R.X.; Yu, D.L.; Li, H.Y. PdbrlA, PdabaA and PdwetA control distinct stages of conidiogenesis in Penicillium digitatum. Res. Microbiol. 2015, 166, 56–65. [Google Scholar] [CrossRef] [PubMed]
  22. Lee, B.N.; Adams, T.H. FluG and flbA function interdependently to initiate conidiophore development in Aspergillus nidulans through brlA beta activation. EMBO J. 1996, 15, 299–309. [Google Scholar] [CrossRef]
  23. Adams, T.H.; Wieser, J.K.; Yu, J.H. Asexual sporulation in Aspergillus nidulans. Microbio. Mol. Biol. Rev. 1998, 62, 35–54. [Google Scholar] [CrossRef] [PubMed]
  24. Seo, J.A.; Guan, Y.J.; Yu, J.H. Suppressor mutations bypass the requirement of fluG for asexual sporulation and sterigmatocystin production in Aspergillus nidulans. Genetics 2003, 165, 1083–1093. [Google Scholar] [CrossRef] [PubMed]
  25. Seo, J.A.; Guan, Y.J.; Yu, J.H. FluG-dependent asexual development in Aspergillus nidulans occurs via derepression. Genetics 2006, 172, 1535–1544. [Google Scholar] [CrossRef] [PubMed]
  26. Wang, F.F.; Krijgsheld, P.; Hulsman, M.; de Bekker, C.; Muller, W.H.; Reinders, M.; de Vries, R.P.; Wösten, H.A.B. FluG affects secretion in colonies of Aspergillus Niger. Antonie. Leeuwenhoek 2015, 107, 225–240. [Google Scholar] [CrossRef] [PubMed]
  27. Park, H.S.; Yu, J.H. Developmental regulators in Aspergillus fumigatus. J. Microbiol. 2016, 54, 223–231. [Google Scholar] [CrossRef] [PubMed]
  28. Yu, J.H. Regulation of development in Aspergillus nidulans and Aspergillus fumigutes. Mycobiology 2010, 38, 229–237. [Google Scholar] [CrossRef]
  29. Etxebeste, O.; Ni, M.; Garzia, A.; Kwon, N.J.; Fischer, R.; Yu, J.H.; Espeso, E.A.; Ugalde, U. Basic-zipper-type transcription factor FlbB controls asexual development in Aspergillus nidulans. Eukaryot. Cell 2008, 7, 38–48. [Google Scholar] [CrossRef]
  30. Etxebeste, O.; Herrero-García, E.; Araújo-Bazan, L.; Rodríguez-Urra, A.B.; Garzia, A.; Ugalde, U.; Espeso, E.A. The bZIP-type transcription factor FlbB regulates distinct morphogenetic stages of colony formation in Aspergillus nidulans. Mol. Microbiol. 2009, 73, 775–789. [Google Scholar] [CrossRef]
  31. Xiao, P.; Shin, K.S.; Wang, T.; Yu, J.H. Aspergillus fumigatus flbB encodes two basic leucine zipper domain (bZIP) proteins required for proper asexual development and gliotoxin production. Eukaryot. Cell 2010, 9, 1711–1723. [Google Scholar] [CrossRef]
  32. Garzia, A.; Etxebeste, O.; Herrero-García, E.; Ugalde, U.; Espeso, E.A. The concerted action of bZip and cMyb transcription factors FlbB and FlbD induces brlA expression and asexual development in Aspergillus nidulans. Mol. Microbiol. 2010, 75, 1314–1324. [Google Scholar] [CrossRef]
  33. Kwon, N.J.; Garzia, A.; Espeso, E.A.; Ugalde, U.; Yu, J.H. FlbC is a putative nuclear C2H2 transcription factor regulating development in Aspergillus nidulans. Mol. Microbiol. 2010, 77, 1203–1219. [Google Scholar] [CrossRef] [PubMed]
  34. Arratia-Quijada, J.; Sanchez, O.; Scazzocchio, C.; Aguirre, J. FlbD, a Myb transcription factor of Aspergillus nidulans, is uniquely involved in both asexual and sexual differentiation. Eukaryot. Cell 2012, 11, 1132–1142. [Google Scholar] [CrossRef]
  35. Garzia, A.; Etxebeste, O.; Herrero-Garcia, E.; Fischer, R.; Espeso, E.A.; Ugalde, U. Aspergillus nidulans FlbE is an upstream developmental activator of conidiation functionally associated with the putative transcription factor FlbB. Mol. Microbiol. 2009, 71, 172–184. [Google Scholar] [CrossRef]
  36. Kwon, N.J.; Shin, K.S.; Yu, J.H. Characterization of the developmental regulator FlbE in Aspergillus fumigatus and Aspergillus nidulans. Fungal Genet. Biol. 2010, 47, 981–993. [Google Scholar] [CrossRef]
  37. Chang, P.K.; Scharfenstein, L.L.; Mack, B.; Ehrlich, K.C. Deletion of the Aspergillus flavus orthologue of A. nidulans fluG reduces conidiation and promotes production of sclerotia but does not abolish aflatoxin biosynthesis. Appl. Environ. Microbiol. 2012, 78, 7557–7563. [Google Scholar] [CrossRef]
  38. Wei, K.Y.; Long, L.K.; Lin, Q.Y.; Ding, S.J. Functional characterization of a new 3-dehydroshikimate dehydratase from Eupenicillium parvum and its potential for protocatechuic acid production. Biosci. Biotechnol. Biochem. 2022, 86, 1024–1030. [Google Scholar] [CrossRef]
  39. Moran, G.R. 4-hydroxyphenylpyruvate dioxygenase. Arch. Biochem. Biophys. 2005, 433, 117–218. [Google Scholar] [CrossRef]
  40. Liu, X.F.; Xia, Y.J.; Zhang, Y.; Liang, L.H.; Xiong, Z.Q.; Wang, G.Q.; Song, X.; Ai, L.Z. Enhancement of antroquinonol production via the overexpression of 4-hydroxybenzoate polyprenyltransferase biosynthesis-related genes in Antrodia cinnamomea. Phytochemistry 2021, 184, 112677. [Google Scholar] [CrossRef]
  41. Bloor, S.; Michurin, I.; Titchiner, G.R.; Leys, D. Prenylated flavins: Structures and mechanisms. FEBS J. 2023, 290, 2232–2245. [Google Scholar] [CrossRef]
  42. Anthony, W.S. NADPH oxidases as electrochemical generators to produce ion fluxes and turgor in fungi, plants and humans. Open Biol. 2016, 6, 160028. [Google Scholar] [CrossRef] [PubMed]
  43. Deng, L.; Sugiura, R.; Ohta, K.; Tada, K.; Suzuki, M.; Hirata, M.; Nakamura, S.I.; Shuntoh, H.; Kuno, T. Phosphatidylinositol-4-phosphate 5-kinase regulates fission yeast cell integrity through a phospholipase C-mediated protein kinase C-independent pathway. J. Biol. Chem. 2005, 280, 27561–27568. [Google Scholar] [CrossRef] [PubMed]
  44. Ogawa, M.; Tokuoka, M.; Jin, F.J.; Takahashi, T.; Koyama, Y. Genetic analysis of conidiation regulatory pathways in koji-mold Aspergillus oryzae. Fungal Genet. Biol. 2010, 47, 10–18. [Google Scholar] [CrossRef] [PubMed]
  45. 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]
  46. Borneman, A.R.; Hynes, M.J.; Andrianopoulos, A. The abaA homologue of Penicillium marneffei participates in two developmental programmes: Conidiation and dimorphic growth. Mol. Microbiol. 2000, 38, 1034–1047. [Google Scholar] [CrossRef] [PubMed]
  47. Mead, M.E.; Borowsky, A.T.; Joehnk, B.; Steenwyk, J.L.; Shen, X.X.; Sil, A.; Rokas, A. Recurrent loss of abaA, a master regulator of asexual development in filamentous fungi, correlates with changes in genomic and morphological traits. Genome Biol. Evol. 2020, 12, 1119–1130. [Google Scholar] [CrossRef] [PubMed]
  48. Zhang, J.G.; Xu, S.Y.; Ying, S.H.; Feng, M.G. Roles of BrlA and AbaA in mediating asexual and insect pathogenic lifecycles of Metarhizium robertsii. J. Fungi 2022, 8, 1110. [Google Scholar] [CrossRef]
  49. Son, H.; Kim, M.G.; Min, K.; Lim, J.Y.; Choi, G.J.; Kim, J.C.; Chae, S.K.; Lee, Y.W. WetA is required for conidiogenesis and conidium maturation in the ascomycete fungus Fusarium graminearum. Eukaryot. Cell 2014, 13, 87–98. [Google Scholar] [CrossRef]
  50. 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, 23, 334. [Google Scholar] [CrossRef]
  51. Adams, T.H.; Hide, W.A.; Yager, L.N.; Lee, B.N. Isolation of a gene required for programmed initiation of development by Aspergillus nidulans. Mol. Cell. Biol. 1992, 12, 3827–3833. [Google Scholar] [CrossRef] [PubMed]
  52. Iradi-Serrano, M.; Tola-Garcia, L.; Cortese, M.S.; Ugalde, U. The early asexual development regulator fluG codes for a putative bifunctional enzyme. Front. Microbiol. 2019, 10, 778. [Google Scholar] [CrossRef]
  53. Serebriiskii, I.G.; Golemis, E.A. Two-hybrid system and false positives. Approaches to detection and elimination. Methods Mol Biol. 2001, 177, 123–134. [Google Scholar] [CrossRef] [PubMed]
  54. Serebriiskii, I.G.; Estojak, J.; Berman, M.; Golemis, E.A. Approaches to detecting false positives in yeast two-hybrid systems. Biotechniques 2000, 28, 328–330, 332–336. [Google Scholar] [CrossRef] [PubMed]
  55. Li, F.; Shi, H.Q.; Ying, S.H.; Feng, M.G. Distinct contributions of one Fe- and two Cu/Zn- cofactored superoxide dismutases to antioxidation, UV tolerance and virulence of Beauveria bassiana. Fungal Genet. Biol. 2015, 81, 160–171. [Google Scholar] [CrossRef] [PubMed]
  56. 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]
  57. Langmead, B.; Salzberg, S.L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 2012, 9, 357–359. [Google Scholar] [CrossRef] [PubMed]
  58. Kanehisa, M.; Goto, S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000, 28, 27–30. [Google Scholar] [CrossRef] [PubMed]
  59. Kanehisa, M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019, 28, 1947–1951. [Google Scholar] [CrossRef]
  60. Kanehisa, M.; Furumichi, M.; Sato, Y.; Ishiguro-Watanabe, M.; Tanabe, M. KEGG: Integrating viruses and cellular organisms. Nucleic Acids Res. 2021, 49, D545–D551. [Google Scholar] [CrossRef]
  61. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCt method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Effects of fluG, brlA, abaA, and wetA on vegetative growth, aerial conidiation, and blastospore production of B. bassiana. (A) Top (row 1), bottom (row 2), and side (row 3) views of fungal colonies initiated with hyphal mass plugs (5 mm diameter) and cultivated for 10 days on SDAY at 25 °C. (B) Conidial yields from SDAY cultures initiated with 100 μL of hyphal suspension per plate and grown for 7 days at 25 °C with a 12:12 h light/dark cycle. (C) Blastospore yields in the submerged SDB cultures over 7 days of incubation at 25 °C and 150 rpm.
Figure 1. Effects of fluG, brlA, abaA, and wetA on vegetative growth, aerial conidiation, and blastospore production of B. bassiana. (A) Top (row 1), bottom (row 2), and side (row 3) views of fungal colonies initiated with hyphal mass plugs (5 mm diameter) and cultivated for 10 days on SDAY at 25 °C. (B) Conidial yields from SDAY cultures initiated with 100 μL of hyphal suspension per plate and grown for 7 days at 25 °C with a 12:12 h light/dark cycle. (C) Blastospore yields in the submerged SDB cultures over 7 days of incubation at 25 °C and 150 rpm.
Ijms 25 06261 g001
Figure 2. Contributions of fluG, brlA, abaA, and wetA to multiple stress responses and virulence of B. bassiana during growth. (AF): EC50 values of chemical stressors required to suppress radial growth by 50% after 7 days of cultivation on CZA plates at 25 °C with 12:12 h light/dark cycle. (G) LT50 (no. of days) for hyphal virulence to G. mellonella larvae inoculated via topical application (immersed). Note: different letters on the bars denote significant differences in each group (Tukey’s HSD, p < 0.05). Error bars: SDs from three replicates.
Figure 2. Contributions of fluG, brlA, abaA, and wetA to multiple stress responses and virulence of B. bassiana during growth. (AF): EC50 values of chemical stressors required to suppress radial growth by 50% after 7 days of cultivation on CZA plates at 25 °C with 12:12 h light/dark cycle. (G) LT50 (no. of days) for hyphal virulence to G. mellonella larvae inoculated via topical application (immersed). Note: different letters on the bars denote significant differences in each group (Tukey’s HSD, p < 0.05). Error bars: SDs from three replicates.
Ijms 25 06261 g002
Figure 3. Effects of fluG, brlA, abaA, and wetA on nutrient utilization of B. bassiana. (AD) Relative colony sizes of B. bassiana strains initiated with hyphal mass plugs with 5 mm diameters, measured as the ratio of the mutant strain’s growth area to the wild strain’s growth area after 7 days of cultivation on various media at 25 °C. Note: different letters on the bars denote significant differences in each group (Tukey’s HSD, p < 0.05). Error bars: SDs from three replicates.
Figure 3. Effects of fluG, brlA, abaA, and wetA on nutrient utilization of B. bassiana. (AD) Relative colony sizes of B. bassiana strains initiated with hyphal mass plugs with 5 mm diameters, measured as the ratio of the mutant strain’s growth area to the wild strain’s growth area after 7 days of cultivation on various media at 25 °C. Note: different letters on the bars denote significant differences in each group (Tukey’s HSD, p < 0.05). Error bars: SDs from three replicates.
Ijms 25 06261 g003
Figure 4. Overview of RNA-Seq data. (A) Histogram of differentially expressed genes (DEGs) from the ΔfluG mutant compared with the WT. Red and blue colored bars represent significantly upregulated and downregulated genes, respectively. (B) Gene ontology (GO) classification of DEGs from the ΔfluG mutant compared with the WT. The enriched GO terms are along the vertical axis, and the horizontal axis indicates the percentage of DEGs in a given term. (C,D) Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of up- and downregulated DEGs from the ΔfluG mutant compared with the WT. The ordinate represents the pathway name, the abscissa represents the enrichment factor, and the point size represents the number of DEGs in that pathway, while the point colors denote the differing Q-value ranges.
Figure 4. Overview of RNA-Seq data. (A) Histogram of differentially expressed genes (DEGs) from the ΔfluG mutant compared with the WT. Red and blue colored bars represent significantly upregulated and downregulated genes, respectively. (B) Gene ontology (GO) classification of DEGs from the ΔfluG mutant compared with the WT. The enriched GO terms are along the vertical axis, and the horizontal axis indicates the percentage of DEGs in a given term. (C,D) Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of up- and downregulated DEGs from the ΔfluG mutant compared with the WT. The ordinate represents the pathway name, the abscissa represents the enrichment factor, and the point size represents the number of DEGs in that pathway, while the point colors denote the differing Q-value ranges.
Ijms 25 06261 g004
Figure 5. Heat map analysis of DEGs in key signaling pathways. (A,B) The relative transcript levels of DEGs in eight important signaling pathways were analyzed and are shown using a heat map.
Figure 5. Heat map analysis of DEGs in key signaling pathways. (A,B) The relative transcript levels of DEGs in eight important signaling pathways were analyzed and are shown using a heat map.
Ijms 25 06261 g005
Figure 6. GO classification and KEGG pathway enrichment analysis of proteins that FluG interacted with. (A,B): GO classification and KEGG pathway enrichment analysis of proteins that were selected by yeast two-hybrid library screening.
Figure 6. GO classification and KEGG pathway enrichment analysis of proteins that FluG interacted with. (A,B): GO classification and KEGG pathway enrichment analysis of proteins that were selected by yeast two-hybrid library screening.
Ijms 25 06261 g006
Figure 7. Genetic model diagram for upstream and central regulators in B. bassiana.
Figure 7. Genetic model diagram for upstream and central regulators in B. bassiana.
Ijms 25 06261 g007
Table 1. Sequence analysis of screening of positive clones using fluG sequence as bait.
Table 1. Sequence analysis of screening of positive clones using fluG sequence as bait.
Clone No.Gene_idGene_NameClone No.Gene_idGene_Name
1BBA_09035phenol 2-monooxygenase34BBA_02100hypothetical protein
2BBA_05531cytochrome b-c1 complex subunit 235BBA_06761zinc finger and BTB domain-containing protein 7A
3BBA_10085SWI/SNF and RSC complexes subunit ssr336BBA_01481ubiquitin carrier protein
4BBA_03721alkaline phosphatase37BBA_06393malate dehydrogenase
5BBA_01606isocitrate/isopropylmalate dehydrogenase38BBA_04259carbon catabolite derepressing protein kinase
6BBA_08294biogenesis of lysosome-related organelles complex 1 subunit KXD139BBA_09482polynucleotide kinase 3 phosphatase
7–11BBA_08577ribosomal protein S740BBA_00422Actin-like protein ARP9
12–17BBA_081823-hydroxybenzoate 6-monooxygenase41BBA_01918hypothetical protein
18BBA_04319transcriptional regulatory protein pro142BBA_10102cyclic-amp-dependent transcription factor atf-2
19–20BBA_03167glycosyl hydrolase 53 domain-containing protein43BBA_06969prolyl-tRNA synthetase
21BBA_07782hypothetical protein44BBA_03016peptidase family protein
22BBA_03527tRNA intron endonuclease45BBA_01844tat pathway signal sequence
23–25BBA_00726Glutamine synthetase46BBA_00241Ssu72-like protein
26–27BBA_02314COP9 signalosome complex subunit 547BBA_08019threonyl-tRNA synthetase
28BBA_00241RNA polymerase II subunit A C-terminal domain phosphatase Ssu72-like protein48BBA_07743serine peptidase
29BBA_02591C6 zinc finger domain protein49BBA_06377RasGEF domain-containing protein
30BBA_08096putative Zn(II)2Cys6 transcription factor50BBA_05170acetamidase/formamidase family protein
31BBA_02631PIR protein repeat protein51BBA_09338catalase/peroxidase HPI/catalase P
32BBA_08791alkaline serine protease AorO52BBA_02558ferulic acid esterase (FaeA)
33BBA_01621nitrogen permease regulator 253BBA_00422Actin-like protein ARP9
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

Li, F.; Zhang, J.; Zhong, H.; Yu, K.; Chen, J. Comprehensive Insights into the Remarkable Function and Regulatory Mechanism of FluG during Asexual Development in Beauveria bassiana. Int. J. Mol. Sci. 2024, 25, 6261. https://doi.org/10.3390/ijms25116261

AMA Style

Li F, Zhang J, Zhong H, Yu K, Chen J. Comprehensive Insights into the Remarkable Function and Regulatory Mechanism of FluG during Asexual Development in Beauveria bassiana. International Journal of Molecular Sciences. 2024; 25(11):6261. https://doi.org/10.3390/ijms25116261

Chicago/Turabian Style

Li, Fang, Juefeng Zhang, Haiying Zhong, Kaili Yu, and Jianming Chen. 2024. "Comprehensive Insights into the Remarkable Function and Regulatory Mechanism of FluG during Asexual Development in Beauveria bassiana" International Journal of Molecular Sciences 25, no. 11: 6261. https://doi.org/10.3390/ijms25116261

APA Style

Li, F., Zhang, J., Zhong, H., Yu, K., & Chen, J. (2024). Comprehensive Insights into the Remarkable Function and Regulatory Mechanism of FluG during Asexual Development in Beauveria bassiana. International Journal of Molecular Sciences, 25(11), 6261. https://doi.org/10.3390/ijms25116261

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