The FlbC transcription factor contributes to the generation of strain heterogeneity in Fumagillin mycotoxin production in Aspergillus fumigatus

The FlbC transcription factor contributes to the generation of strain heterogeneity in Fumagillin mycotoxin production in Aspergillus fumigatus | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Article The FlbC transcription factor contributes to the generation of strain heterogeneity in Fumagillin mycotoxin production in Aspergillus fumigatus Yutaka Tanaka, Ikumi Asagiri, Rei Komiya, Hiroyuki Yamazaki, Akiho Yagi, and 2 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-4435109/v1 This work is licensed under a CC BY 4.0 License Status: Posted Version 1 posted You are reading this latest preprint version Abstract Aspergillus fumigatus shows inter-strain heterogeneity in the repertoire of secondary metabolites such as mycotoxin fumagillin but the mechanism by which this heterogeneity arises in those production is still not understood. Here, we investigates the effect of the flbC gene on fumagillin production in A. fumigatus by introducing flbC deletions into laboratory strains, A1151, A1159, and A1280 from different backgrounds and examines the effect of the flbC gene on inter-strain heterogeneity. We found that, although all three laboratory strains were capable of producing fumagillin, there was heterogeneity in the effects of flbC gene deletion on fumagillin production. This heterogeneity may be dependent on differences in the expression levels of the fma gene family involved in fumagillin production and different levels of transcriptional activation by transcription factors FapR and LaeA. The flbC gene expression level peaks at the same time as peak mycelial growth but shows differential expression level, suggesting that the differences among strains in the range of expression levels are dependent on differences in the upstream expression of flbC . Thus, our findings show that the different interactions between flbC and factors regulating the expression of fumagillin gene cluster are the probable cause for heterogeneity in fumagillin production in the organism. Aspergillus fumigatus flbC fumagillin secondary metabolism heterogeneity Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Figure 7 Introduction Aspergillus fumigatus is an opportunistic pathogen, which infects human beings. Although A. fumigatus infection in healthy individuals has not established, persons with a history of chemotherapy, immunosuppressive drug use, or congenital immunodeficiency syndrome can have lethal lesions 1 . The human-pathogenicity characteristic of A. fumigatus is determined by the interaction of several virulence factors in the organism. The mycelia produce secreted proteins 2 , polysaccharides 3 , and a wide variety of secondary metabolites 4 . In recent years, inter-species, inter-strain, and even individual spore heterogeneity in these phenotypes and genotypes has been noted 5 – 10 . Heterogeneity has also been observed between clinical and environmental isolates 11 . In basic research, several strains with different backgrounds are currently used as reference strains; however, notable differences have been reported between these strains in terms of pathogenicity 12 , genome structure 12 , 13 , secondary-metabolite production capacity 14 , single nucleotide polymorphisms (SNPs) 15 , chromatin modifications 16 , and transposable element distribution 16 . Vigorous research is now underway to understand how this heterogeneity arises between species or strains, mainly through large-scale data. Af293 is the most commonly used A. fumigatus strain for phenotyping in basic research because it was the first strain to have its entire genome sequenced and is currently the only strain that has been allocated to the chromosome 17 . In addition, CEA17 18 , D141 19 , and ATCC46645 20 strains have been established and used as laboratory strains. Heterogeneity among strains has made it more complicated to understand the physiology of the species by comparison between each study. Although the secondary metabolites produced by Aspergillus spp. are not explicitly associated with pathogenicity in humans, a correlation between the repertoire of secondary metabolites that can be produced and pathogenicity has been reported 21 . Inter-species and inter-strain heterogeneity is also observed in the repertoire of secondary metabolites 21 . Many A. fumigatus strains that secrete fumagillin, pseurotin, fumitremorgin, and gliotoxin are pathogenic, and closely related species of A. fumigatus , A. oerlinghausenensis and A. fischeri can produce fumagillin and pseurotin as secondary metabolites; however, these organisms are not pathogenic to humans. Comparisons between strains with different repertoires of secondary metabolite secretion have shown that strains with a more severe clinical course tend to secrete a greater variety of metabolites. Although no single secondary metabolite is known to determine the human pathogenicity of A. fumigatus , the relationship between strain heterogeneity in secondary-metabolite production repertoire and pathogenicity remains unclear. The fma gene family 22 , 23 , required for fumagillin biosynthesis in A. fumigatus , forms a cluster in the subtelomeric region of chromosome 8 24 . The gene cluster has eight Fma enzymes that catalyze sequential reactions to biosynthesize the end product and the transcription factor FapR/FumR 24,25 , which positively regulates the expression of the fma gene family upstream. Because fumagillin precursors are secreted by strains lacking each gene of the fma gene family, each fma gene is essential for the biosynthesis of fumagillin 23 . Loss of the gene encoding the histone methyltransferase LaeA affects the production of fumagillin and other secondary metabolites 26 , 27 . Inside the fumagillin gene cluster, there is an intertwined group of genes involved in pseurotin biosynthesis, which is synchronously regulated by FapR and LaeA 24 . In addition, a gene cluster for the biosynthesis of fumitremorgin is present in the neighborhood 24 , suggesting that fumagillin-pseurotin-fumitremorgin is synchronously affected by the regulation of chromatin structure. Chromosome 8 has a lower frequency of trimethylation modifications of the ninth lysine residue of histone H3 than other genomic regions and has many LINE transposon sequences 16 , suggesting that the area around the fumagillin gene cluster may be a region of relatively high genomic instability. Consistent with this, approximately 60% of A. fumigatus isolates have lost a part of or almost all the fumagillin gene cluster in the subtelomeric region of chromosome 8 15 . Thus, the role of the fma gene cluster and its upstream transcription factors in the production of fumagillin, a typical secondary metabolite, is becoming clearer; however, the mechanism by which this heterogeneity arises in the fumagillin production is not yet known. Secondary metabolite production is often regulated synchronously between several secondary-metabolite-producing systems rather than independently of each other 28 . Lind et al. reported that the BrlA transcription factor, which is central to the central developmental pathway, affects the expression of almost all secondary metabolite products 29 . The central developmental pathway is based on sequential transcriptional activation of a series of transcription factors, namely BrlA, AbaA, and WetA, which are important for the transcription of genes involved in morphogenesis 30 . During the asexual life cycle of A. fumigatus , dormant spores first swell and germinate to form filamentous mycelia, followed by the asexual reproductive spores. These series of shifts are synchronized with the continuous activation of the central developmental pathway. The comparison of the mRNA expression profiles of laeA- and brlA- gene deficient strains and the secondary-metabolic production profiles have shown that despite the high correlation between the two profiles and the fumagillin-fumitremorgin-pseurotin gene cluster, there is negative correlation between the expression patterns of the these gene clusters 29 . The central developmental pathway regulated by BrlA is further regulated upstream by several Flb transcription factors, which are together termed the upstream developmental activator (UDA) pathway. In Aspergilus nidulans , FlbB, FlbC, and FlbE are transcription factors that constitute the UDA pathway 31 . FlbB and FlbE cooperate, while FlbC alone activates the transcription of brlA 31 . In A. fumigatus strains deficient in flbB and flbE genes, reduced brlA gene expression, delayed mycelial growth, and reduced fumagillin and gliotoxin production have been observed 32 . However, the effects of the flbC gene on A. fumigatus growth, brlA gene expression, and secondary metabolite production have not yet been investigated in detail. In addition, in the aforementioned comparisons of brlA , laeA , and flbB/E gene-deficient strains, experiments were conducted using a single reference strain for each strain and inter-strain heterogeneity was not considered. Therefore, this study aimed to investigate the effect of the flbC gene on fumagillin production in A. fumigatus by introducing flbC deletions into several laboratory strains from different backgrounds and to consider the effect of the flbC gene on inter-strain heterogeneity. In particular, changes in the expression pattern of the brlA gene over time, fumagillin production capacity, expression of the fma gene cluster, and expression of the upstream transcription factors laeA and fapR were compared. Results Expression patterns of the flbC gene in the three strains A1151, A1159, and A1280 as laboratory strains with suppressed non-homologous repair for genetic modification were deposited in FGSC. The strains CEA17, D141, and ATCC46645 were the parental strains of A1151, A1159, and A1280, respectively 33 . The flbC gene in strain Af293 ( Afu2g13770 ) was searched in the draft genomes of strains CEA17 and found AFUB029400-T . Using this sequence as a reference, primers were designed to amplify the entire CDS length of the flbC gene. Amplification was attempted using the genomic DNA of the three strains as a template, and amplification products were obtained in all groups. The sequences in the CDS of the amplified products were checked for point mutations and were found to be completely homologous in all three strains, suggesting that there is probably no difference in the amino acid sequence of the FlbC protein produced in the three strains. Next, to compare the expression patterns of the flbC gene in the three strains, RNA was extracted from spores and mycelia obtained at different incubation times, and the relative amounts of flbC gene transcripts were quantified over time. The flbC gene transcript was not amplified immediately after the start of incubation, when almost all spores formed buds, but was detected after approximately 8 h and reached peak expression in all three strains after 16 h (Fig. 1 ). Thereafter, a gradual downward trend was observed in A1151 and A1280, followed by a slight increase in A1159 at 72 h. These trends are consistent with previous research, in which the flbC transcript levels in A. nidulans during vegetative growth peaked at approximately 18–24 h and decreased thereafter 34 . However, some heterogeneity appeared among the three strains, with expression ratios of approximately 3-fold in A1151, 4-fold in A1159, and 9-fold in A1280 at 16 h. In summary, the flbC gene appears to reach peak expression level approximately 16 h after the start of culture, that is, at the same time as peak mycelial growth. Comparison of three representative laboratory strains suggested that the CDS internal sequence of the flbC gene and the timing of the maximum expression level were generally the same, but the range of expression changes was slightly different. These results suggest that the differences among strains in the range of expression changes are dependent on differences in upstream expression of flbC . Effect of flbC expression heterogeneity on strain phenotype To investigate the effect of heterogeneity in the flbC expression pattern on the phenotype of the three laboratory strains, we produced flbC gene disruption mutant strains (Δ flbC ) in each background. Because the internal sequence of flbC was exactly the same, we planned highly efficient gene disruption using the CRISPR-Cas9 method with the PAM sequence common to the three strains. The Δ flbCs were generated by inducing double-strand breaks via Cas9 (CRISPR#1, 2) binding, followed by insertion of pyrithiamine resistance gene markers through homologous recombination, using each of the three strains as parent (Fig. 2 a, b). Δ flbCs were slightly reduced at radial growth circle than the parental strains (WT) on MM at 37°C (Fig. 2 c). The three WT strains grew differentially under these conditions, with the A1280 strain growing slightly weaker on day 3 or earlier than the other two. The radial growth of the flbC gene-deficient strain was almost identical to this trend. Notably, when grown with the cell wall effector Congo red, A1280 showed weaker growth among the WT strains; however, growth was restored by flbC gene deletion (Fig. 2 d). This was in contrast to A1151 and A1159, whose growth was suppressed by Congo red treatment. These findings suggest that the differences in the growth phenotypes observed in the three laboratory strains may be due to the function of the flbC gene. Effect of flbC expression patterns on brlA expression BrlA, which plays a central- and entry point- role in the central developmental pathway, is regulated by FlbC transcription factor 30 , 34 . To investigate how different expression patterns of the flbC gene affect the brlA expression, we measured brlA transcript over time in WT and Δ flbC of the three laboratory strains. WT of the A1151 and A1159-background reached peak mRNA expression at 16–24 h, with an approximate 3- to 4-fold increase than the expression at 4 h, whereas this increase was suppressed in the Δ flbC (Fig. 3 , upper and middle panel). In contrast, the WT of A1280-background showed little increase in brlA gene expression until 48 h after the start of incubation, and Δ flbC in the A1280-background also showed little change in brlA gene expression (Fig. 3 , bottom panel). In summary, brlA gene expression level in WT of A1151 and A1159 peaked at 16–24 h, whereas this peak was removed by flbC gene deletion. In contrast, A1280 showed a flat expression up to 48 h in brlA gene expression in both WT and Δ flbC rather than a peak-like pattern. These results suggest that the brlA gene expression pattern is different among the three laboratory strains and that there is heterogeneity in the regulatory role of the flbC gene on brlA gene expression. Effect of flbC gene on the heterogeneity of secondary metabolite production To investigate the effect of the flbC gene on the heterogeneity of secondary metabolite production, we measured the fumagillin-production capacity of these strains. Remarkably, given that previous studies have shown the loss of fumagillin-producing capacity in several A. fumigatus strains 15 , all three WT strains produced fumagillin (Fig. 4 ). The effect of flbC gene deletion on fumagillin production was tripartite: in A1151, fumagillin production was below the detection limit due to flbC gene deletion; in A1159, fumagillin production was reduced but not eliminated due to flbC gene deletion; and A1280 produced sufficient amount of fumagillin in both WT and Δ flbC , with little difference in production level (Fig. 4 ). These results indicate that although all three laboratory strains are capable of producing fumagillin, there is heterogeneity in the effects of the flbC gene deletion. To investigate whether the heterogeneity in producing fumagillin is dependent on differences in the expression levels of the genes responsible for fumagillin biosynthesis, we measured the transcripts belonging to fma gene family involved in the fumagillin gene cluster ( fma-TC / fmaA , fma-P450 / fmaG , fma-C6H / fmaF , fma-AT / fmaC , and fma-ABM / fmaE ) and those not involved in fumagillin biosynthesis but located within the gene cluster ( Afu8g00430 , Afu8g00440 / psoF and Afu8g00500 ). In A1151-background, the expression levels of many fumagillin biosynthesis-related genes, including fma-TC / fmaA , fma-P450 / fmaG , fma-C6H / fmaF , and fma-AT / fmaC , were significantly reduced in Δ flbC than those in WT (Fig. 5 a). Some genes not involved in fumagillin biosynthesis were also downregulated (Fig. 5 b). A comparison of the A1159 background strains showed the same trend, except for fma-AT / fmaC (Fig. 5 a). In contrast, comparisons of A1280 background strains showed significant differences in the reduced expression of fma-TC / fmaA and fma-C6H / fmaF , but the extent of the reduction was weaker than that in A1151 and A1159, and differences in the expression of other fma genes was not significant (Fig. 5 a). In summary, in the Δ flbC of A1151 and A1159 backgrounds, the expression of the fma gene family was generally and coordinately suppressed, consistent with the reduced fumagillin production, while in A1280, where flbC gene deletion did not result in reduced fumagillin production, the degree of suppression of fma gene family expression was weaker than in the other two strain. Thus, the heterogeneity in the three laboratory strains between the WT and Δ flbC with regard to fumagillin production may be dependent on differences in the expression levels of the fma gene family and probably correlates with the degree of hetero-chromatinization of the fumagillin gene cluster. Expression of the upstream transcription factors laeA and fapR The production of secondary metabolites requires the coordinated action of global and cluster-specific transcription factors. In the case of fumagillin, the transcriptional activation of biosynthetic genes is achieved by the histone methyltransferase LaeA, which acts as a global transcription factor, and FapR, which acts as a cluster-specific transcription factor. The fapR gene is located within the fumagillin gene cluster, and production of the FapR transcription factor is initiated by the release of heterochromatin from the cluster by LaeA. To determine whether there are differences in the upstream gene of the fma gene family, we measured the expression levels of laeA and fapR genes in the A1151 and A1280 backgrounds between WTs and Δ flbC s, which had contrasting difference in the fumagillin production. In the A1151-background strains, laeA gene expression peaked after 16–24 h and then decreased in the WT, whereas the peak was delayed to 48–72 h in the Δ flbC (Fig. 6 a, left panel). In contrast, in A1280, the expression of the laeA gene peaked at 16 h in the WT as same as that in A1151, while the Δ flbC showed a similar pattern to that of the WT (Fig. 6 a, right panel). The fapR gene showed reduced gene expression at 16–24 h in the Δ flbC than that in the WT of the A1151 background (Fig. 6 b, left panel), while there was no significant difference in expression levels at 24 h in the A1280 background (Fig. 6 b, right panel). In summary, in the A1151 background strain, the flbC gene deletion resulted in delayed expression of the laeA gene and suppressed the expression of fapR , but not in the A1280 background. These results indicate that the heterogeneity in fumagillin production between strains A1151 and A1280 may be due to the different levels of transcriptional activation of FapR and LaeA. More specifically, this suggests a dependence on different direct or indirect interactions between flbC and genes regulating the expression of the fumagillin gene cluster, such as laeA and fapR . Effect of heterogeneities on virulence The analysis of the three laboratory strains shows that the deletion of the flbC gene may create new heterogeneity or reduce that which was present originally. To estimate the effects of these heterogeneities on virulence, a pathogenicity test was conducted on silkworms. We found that the deletion of the flbC gene probably does not contribute to silkworm-pathogenicity of the strains because there was no significant difference in the survival among WT and Δ flbC of any background (Fig. 7 ). Deletion of the flbC gene reduced the growth of the laboratory strains on MM and showed differences in the production of the mycotoxin fumagillin; however, these differences were not observed in the pathogenicity test. Discussion This study shows that the flbC gene in A. fumigatus may increase or decrease the heterogeneity of hyphal growth, the expression patterns of brlA transcripts and secondary metabolites. FlbC is a C 2 H 2 transcription factor that is widely conserved in Ascomycetes, but not in Basidiomycetes and Zygomycetes 34 . The flbC gene has been found in A. nidulans 34 , A. flavus 35 , Trichoderma guizhouense 36 , and Beauveria bassiana 37 . In A. oryzae , it regulates the production of amylolytic enzymes, such as glucoamylase and α-amylase, during individual media growth 38 . Fusarium graminearum deficient of flbC gene showed no phenotype 39 . FlbC regulates brlA gene expression, along with the concerted action of FlbB, FlbD, and FlbE 30 . Thus, the extent to which any flb family molecule contributes to the asexual life cycle varies among species. These Flbs also have differnt localizations in fungal cells, with FlbB and FlbE tending to accumulate at the tips of filamentous fungi, while FlbC is present in all nuclei, including metulae and phialides 34 . This is not surprising, given that the plasticity and velocity of filamentous growth, as well as the cell wall and cell surface structures, differ between species 40 . Genomic heterogeneity between clinical and environmental isolates is known, and transcriptome heterogeneity exists within the same strain, even within individual spores 5 , 9 – 11 . The reason for the various levels of heterogeneity arising in A. fumigatus will become increasingly clear in future studies. As we showed, in the A1151 and A1159 backgrounds, loss of the flbC gene was associated with reduced fumagillin production. These strains showed reduced expression of most members of the fma gene family and, in correlation, delayed the expression of laeA gene and suppressed the expression of fapR gene in A1151 cells. The effect of laeA and fapR gene expression on fumagillin production may be explained by the following facts presented by Dhingra et al. 25 . The fapR/fumR gene-deficient strain completely loses its ability to produce fumagillin, and the expression of the fma gene family is also almost completely suppressed. The expression levels of fapR/fumR and fmaB are suppressed by more than 90% in the laeA gene-deficient strain, and fmaB gene expression is also suppressed in the fapR gene-deficient strain. In strains deficient in the VeA transcription factor, which forms a heterotrimer with the LaeA protein, fapR gene expression is repressed, and fumagillin production is reduced by more than 80%. Thus the expression of the fma gene family; global transcription factors such as LaeA and VeA, which are upregulated; and FapR, a fumagillin cluster-specific transcription factor, play important roles in fumagillin production. However, the mechanism by which flbC gene deletion causes delayed expression of the laeA gene remains unclear. Several studies have suggested a link between central developmental pathways and secondary metabolic production. Lind et al. reported altered production patterns of many secondary metabolic products in brlA and laeA gene deletion strains 29 . Kim et al. have also shown that deletion of the flbB and flbE genes results in repression of the laeA gene and reduced production of gliotoxin and fumagillin 32 . However, they did not discuss the mechanisms of laeA repression, but showed an altered expression pattern of the brlA gene due to flb gene deletion 41 . Moon et al. showed by multi-omics analysis that the flbC gene may be a direct target of LaeA and VeA transcription factors 42 . Rocha et al. revealed the regulation of flbC gene expression by the MAPK transcription factor RlmA 43 . A feedback regulatory mechanism may be present between FlbC and global transcription factors such as LaeA and RlmA. As we have shown, the Congo red sensitivity of the A1280 WT strain was restored by the loss of flbC gene, which suggests a link between the cell wall integrity pathway and flbC . Loss of the brlA gene, which is thought to be located downstream of FlbC, results in increased production of fumagillin 29 . This seems to contradict our finding that flbC gene deletion reduced fumagillin production in A1151 and A1159. This may be due to perturbations in LaeA or FapR expression. However, we may have to consider the compensatory effects induced by gene loss (or overexpression), as is the case when certain important genes are involved in the basic properties and survival of mycelia. Considering that BrlA is an important transcription factor that plays a critical role in central developmental pathways, and that it may be important for brlA expression to be oscillatory, unidirectionally increasing or decreasing its expression may have drastic effects 44 , 45 . This, in turn, is also true for flbC , and is a refutation and claim of the present study. The fact that fumagillin was reduced in flbC -deficient strains may not necessarily indicate that the FlbC transcription factor is “directly” involved in fumagillin production. Notably, among the three standard laboratory strains, heterogeneity was observed not only in the gene-deficient strains but also in the WT strain. A1151, A1159, and A1280 were clinical isolates, but each isolate was isolated from a person with a different disease condition 33 . Comparisons of the three strains have already been made at various levels, including whole genome, transcriptome, metabolome, leukemia model mouse infectivity, and macrophage phagocytosis capacity 12 , 13 , 17 , 46 . The fumagillin gene cluster is located in the subtelomeric region of chromosome 8, where transposable elements and chromatin modifications are frequently found in Af293 and CEA17 backgrounds 16 , 29 . Perrin et al. focused on the genomic instability of secondary-metabolite gene clusters and suggested that the expression profiles of secondary metabolites may differ among strains, showing that the structure of the subtelomeric region varies among strains, even within the same species 27 . Thereafter, Lind et al. organized these mechanisms and divided them into different types 29 . The present study just demonstrates this hypothetical model. However, the heterogeneity in fumagillin production identified in this study is probably not due to genomic instability in the subtelomeric region of chromosome 8 because WT of A1151, A1159, and A1280 expressed fma genes and produced fumagillin. According to Lind et al., the fumagillin gene cluster is relatively stable because it has fewer SNPs with non-synonymous substitutions than the other genomic regions 29 . FlbC of the three strains behaves slightly differently with regard to the regulation of the fumagillin gene cluster; however, understanding this difference will probably require more than a comparison of genomic DNA between the strains. Nevertheless, with the advent of the CRISPR-Cas9, as we have used, recombinants of various strains from one species can be obtained efficiently and quickly. Although we obtained mutants by homologous recombination of drug resistance markers, “co-genome editing 47 ,” which does not use homologous recombination of markers, is now being developed. This enables the genome editing of a large number of clinical and environmental isolates and the use of these isolates as parental strains. Comparison of more strains than that in this study, which compared only three strains, is now possible. Thus, our findings show that the different interactions between flbC and factors (LaeA and FapR) regulating the expression of fumagillin gene cluster are the probable cause for heterogeneity in fumagillin production among the strains. Moreover, our study serves as a basis for future studies on co-genome editing of large number of clinical and environmental isolates and the use of these isolates as parental strains. Methods Strains, media, and culture conditions The A. fumigatus strains used in this study are listed in Supplementary Table S1 . A. fumigatus strains A1151, A1159, and A1280 were obtained from the Fungal Genetics Stock Center (FGSC). The strains were grown on minimal medium [MM; 1% (w/v) glucose, 0.6% (w/v) NaNO 3 , 0.052% (w/v) KCl, 0.052% (w/v) MgSO 4 ·7 H 2 O, 0.152% (w/v) KH 2 PO 4 , and Hunter's trace elements, pH 6.5]. Plasmids were amplified in E. coli DH5α. Colony growth rates were measured as previously described 48 . Briefly, conidia from each strain were point-inoculated onto the center of MM agar plates with or without 400 µg/mL CongoRed (Fujifilm, Japan). Colony diameters were measured after 3, 4, and 5 d of incubation at 37°C. The growth rates of all the individual strains were measured 5 times. Generation of flbC mutants Each sgRNA target sequence was explored using the web service CHOPCHOP ( https://chopchop.cbu.uib.no/ ) and two synthesized DNA oligonucleotides, primer No. 3 and 4 (Supplementary Table S2 ). These sequences contained the T7 promoter, an sgRNA target sequence, and a sequence overlapping the Cas9 scaffold. The synthesis and purification of sgRNAs were performed using the Guide-it sgRNA In Vitro Transcription Kit (TAKARA, Japan). After the synthesized sgRNAs (sgRNA-1 and − 2) were quantified using a NanoDrop spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA), they were used for ribonucleoprotein (RNP) formation with Cas9. Repair templates were constructed based on the method described by Arai et al., with slight modifications 49 . The pyrithiamine resistance cassette ( ptrA marker) was obtained from the pPTRII plasmid (TAKARA). The KOD One PCR Master Mix (TOYOBO, Japan) and primer No. 1 and 2 (Supplementary Table S2 ) were used for PCR amplification to construct the repair template. The PCR products were purified using a FastGene Gel/PCR Extraction kit (NIPPON Genetics Co., Ltd., Japan) and quantified using a NanoDrop spectrophotometer. A. fumigatus protoplasts were generated, and fungal transformation was performed as described by Kadooka et al., with slight modifications 50 . Briefly, conidia were incubated in YG mediem (2% glucose, 0.5% yeast extract, and Hunter’s trace element) for 16 h at 37°C. Following incubation, the cell walls of germlings were digested with 32 mg/mL Vino-Taste Pro (Novozymes, Denmark) for 4 h 30 min at 30°C. Subsequently, 10 µg Guide-it Recombinant Cas9 (TAKARA) and 10 pmol each of sgRNA-1 and − 2 were mixed and incubated for 5 min, generating RNPs as described by Arai et al.. 49 Protoplasts were transformed with 2–3 µg of repair templates and RNPs and plated onto MM agar supplemented with 0.6 M KCl. Following 24 h incubation at room temperature, agar plates were overlaid with MM top agar containing 0.1 µg/mL pyrithiamine hydrobromide (Sigma-Aldrich, USA). The flbC mutants were isolated and confirmed using colony PCR with primer No. 5 and 6 (Supplementary Table S2 ). At least three deletion strains were isolated in each case. Quantitative real-time PCR analysis Quantification of mRNA expression was performed using THUNDERBIRD SYBR qPCR mix (TOYOBO) and StepOnePlus Real-Time PCR system (Thermo Fisher Scientific). The strains were grown at 37°C in liquid MM at each time point. Fungal RNA was isolated using the Sepasol-RNA I Super G (Nacalai Tesque, Japan) and reverse-transcribed into cDNA using 1 µg total RNA and ReverTra Ace qPCR RT Kit (TOYOBO). RT-PCR was performed using primers No. 7–32 (Supplementary Table S2 ), and gene expression levels were calculated using ΔΔCt method. HPLC analysis Each strain, maintained on a PDA plate, was inoculated into a 100-mL Erlenmeyer flask containing 50 mL of the production medium (3.0% sucrose, 3.0% soluble starch, 1.0% malt extract, 0.30% Ebios, 0.50% KH 2 PO 4 , and 0.050% MgSO 4 ·7H 2 O; adjusted to pH 6.0 before sterilization). The production culture was carried out at 25°C for 7 days under agitation (150 rpm). Each culture broth was treated with 50 mL of acetone after 7 days and filtered. The filtrate, after the evaporation of acetone, was extracted with EtOAc (50 mL), and the extract was concentrated in vacuo . The production of secondary metabolites was monitored using a Chromaster HPLC system (Hitachi High Technologies Co., Ltd., Tokyo, Japan) [column, CAPCELL PAK C18 MG-II (Osaka Soda. Co., Ltd., Osaka, Japan), i.d. 4.6 mm × 250 mm; flow rate, 0.8 mL/min; mobile phase, linear gradient from 20 to 100% CH 3 CN containing 0.05% TFA for 30 min and 100% CH 3 CN containing 0.05% TFA for 5 min; detection, UV 365 nm]. Each extract was prepared to 10 mg/mL CH 3 OH solution, and 5.0 µL of each sample was injected for the HPLC analysis. Silkworm larva infection model Silkworm ( Bombyx mori ) larvae were maintained in an incubator at 27°C in the dark at a controlled temperature until conidial infection. A. fumigatus conidia were resuspended at a concentration of 1 × 10 8 conidia/mL in phosphate buffered saline (PBS), and then silkworm 5th-instar larvae (2.0 g, n = 10) were infected with 1 × 10 6 conidia (50 µL/larva) by injecting the hemolymph using a disposable 1 mL syringe with a 27G needle (TERUMO, Tokyo, Japan). The uninfected controls received PBS. Larvae were raised in the dark at 37°C for 4 d. Generalized Wilcoxon tests were performed using GraphPad Prism 8 software to compare survival curves, and differences were considered statistically significant at p < 0.05. Fertilized silkworm eggs of Bombyx mori (Hu·Yo × Tukuba·Ne) were purchased from Ehime Sansyu (Ehime, Japan) and fed an artificial diet (Silk Mate 2S; Nihon Nosan Kogyo, Kanagawa, Japan, and Silkmate; Katakura Industries, Tokyo, Japan). Declarations Acknowledgments This work was supported in part by Grant-in-Aid for Scientific Research (C) (22K06600 to Y.T., , 21K06631 and 24K09865 to H.Y.) from the Japan Society for the Promotion of Science (JSPS KAKENHI). We express our thanks to Ms. M. Iwamoto of Tohoku Medical and Pharmaceutical University for her technical assistance. We would like to thank Editage (www.editage.jp) for English language editing. Author contributions Y.T. conceived and designed the experiments. Y.T., I.A, R.K., H.Y. and A.Y. conducted the experiments. R.U. and M.S. helped with the design and coordination of the study. Y.T., H.Y. and A. Y. wrote the paper. All authors reviewed the manuscript. Competing interests The authors declare no competing interests. Data availability statement All data generated or analyzed during this study are included in Supplementary Table S3. A full-length version of the figure obtained by agarose gel electrophoresis is shown in Supplementary FigS1, and no contrast adjustment was made when used in the figures. References Fisher, M. C. & Denning, D. W. The WHO fungal priority pathogens list as a game-changer. Nat. Rev. Microbiol. 21 , 211-212, doi:10.1038/s41579-023-00861-x (2023). Wartenberg, D. et al. Secretome analysis of Aspergillus fumigatus reveals Asp-hemolysin as a major secreted protein. Int. J. Med. Microbiol. 301 , 602-611, doi:10.1016/j.ijmm.2011.04.016 (2011). Briard, B. et al. Galactosaminogalactan activates the inflammasome to provide host protection. Nature 588 , 688-692, doi:10.1038/s41586-020-2996-z (2020). Manfiolli, A. O. et al. Mitogen-Activated Protein Kinase Cross-Talk Interaction Modulates the Production of Melanins in Aspergillus fumigatus. mBio 10 , doi:10.1128/mBio.00215-19 (2019). Danion, F. et al. Aspergillus fumigatus, One Uninucleate Species with Disparate Offspring. J Fungi (Basel) 7 , doi:10.3390/jof7010030 (2021). Hewitt, S. K., Foster, D. S., Dyer, P. S. & Avery, S. V. Phenotypic heterogeneity in fungi: Importance and methodology. Fungal Biol. Rev. 30 , 176-184, doi:10.1016/j.fbr.2016.09.002 (2016). Ijadpanahsaravi, M., Teertstra, W. R. & Wosten, H. A. B. Inter- and intra-species heterogeneity in germination of Aspergillus conidia. Antonie Van Leeuwenhoek 115 , 1151-1164, doi:10.1007/s10482-022-01762-4 (2022). Keller, N. P. Heterogeneity Confounds Establishment of "a" Model Microbial Strain. mBio 8 , doi:10.1128/mBio.00135-17 (2017). Wang, F. et al. Transcription in fungal conidia before dormancy produces phenotypically variable conidia that maximize survival in different environments. Nat Microbiol 6 , 1066-1081, doi:10.1038/s41564-021-00922-y (2021). Dos Santos, R. A. C. et al. Genomic and Phenotypic Heterogeneity of Clinical Isolates of the Human Pathogens Aspergillus fumigatus, Aspergillus lentulus, and Aspergillus fumigatiaffinis. Front Genet 11 , 459, doi:10.3389/fgene.2020.00459 (2020). Horta, M. A. C. et al. Examination of Genome-Wide Ortholog Variation in Clinical and Environmental Isolates of the Fungal Pathogen Aspergillus fumigatus. mBio 13 , e0151922, doi:10.1128/mbio.01519-22 (2022). Kowalski, C. H. et al. Heterogeneity among Isolates Reveals that Fitness in Low Oxygen Correlates with Aspergillus fumigatus Virulence. mBio 7 , doi:10.1128/mBio.01515-16 (2016). Garcia-Rubio, R., Monzon, S., Alcazar-Fuoli, L., Cuesta, I. & Mellado, E. Genome-Wide Comparative Analysis of Aspergillus fumigatus Strains: The Reference Genome as a Matter of Concern. Genes (Basel) 9 , doi:10.3390/genes9070363 (2018). Lind, A. L. et al. Drivers of genetic diversity in secondary metabolic gene clusters within a fungal species. PLoS Biol. 15 , e2003583, doi:10.1371/journal.pbio.2003583 (2017). Barber, A. E. et al. Aspergillus fumigatus pan-genome analysis identifies genetic variants associated with human infection. Nat Microbiol 6 , 1526-1536, doi:10.1038/s41564-021-00993-x (2021). Colabardini, A. C. et al. Chromatin profiling reveals heterogeneity in clinical isolates of the human pathogen Aspergillus fumigatus. PLoS Genet. 18 , e1010001, doi:10.1371/journal.pgen.1010001 (2022). Fedorova, N. D. et al. Genomic islands in the pathogenic filamentous fungus Aspergillus fumigatus. PLoS Genet. 4 , e1000046, doi:10.1371/journal.pgen.1000046 (2008). Kress, M. R. V. Z. et al. The akuBKU80 Mutant Deficient for Nonhomologous End Joining Is a Powerful Tool for Analyzing Pathogenicity in Aspergillus fumigatus. EUKARYOTIC CELL 5 , 207-211, doi:10.1128/EC.5.1.207 (2006). Wagener, J. et al. The putative alpha-1,2-mannosyltransferase AfMnt1 of the opportunistic fungal pathogen Aspergillus fumigatus is required for cell wall stability and full virulence. Eukaryot Cell 7 , 1661-1673, doi:10.1128/EC.00221-08 (2008). Hartmann, T. et al. Validation of a self-excising marker in the human pathogen Aspergillus fumigatus by employing the beta-rec/six site-specific recombination system. Appl. Environ. Microbiol. 76 , 6313-6317, doi:10.1128/AEM.00882-10 (2010). Steenwyk, J. L. et al. Variation Among Biosynthetic Gene Clusters, Secondary Metabolite Profiles, and Cards of Virulence Across Aspergillus Species. Genetics 216 , 481-497, doi:10.1534/genetics.120.303549 (2020). Lin, H. C. et al. The fumagillin biosynthetic gene cluster in Aspergillus fumigatus encodes a cryptic terpene cyclase involved in the formation of beta-trans-bergamotene. J. Am. Chem. Soc. 135 , 4616-4619, doi:10.1021/ja312503y (2013). Lin, H. C. et al. Generation of complexity in fungal terpene biosynthesis: discovery of a multifunctional cytochrome P450 in the fumagillin pathway. J. Am. Chem. Soc. 136 , 4426-4436, doi:10.1021/ja500881e (2014). Wiemann, P. et al. Prototype of an intertwined secondary-metabolite supercluster. Proc. Natl. Acad. Sci. U. S. A. 110 , 17065-17070, doi:10.1073/pnas.1313258110 (2013). Dhingra, S. et al. The fumagillin gene cluster, an example of hundreds of genes under veA control in Aspergillus fumigatus. PLoS One 8 , e77147, doi:10.1371/journal.pone.0077147 (2013). Keller, N. P., Turner, G. & Bennett, J. W. Fungal secondary metabolism - from biochemistry to genomics. Nat. Rev. Microbiol. 3 , 937-947, doi:10.1038/nrmicro1286 (2005). Perrin, R. M. et al. Transcriptional regulation of chemical diversity in Aspergillus fumigatus by LaeA. PLoS Pathog. 3 , e50, doi:10.1371/journal.ppat.0030050 (2007). Brakhage, A. A. Regulation of fungal secondary metabolism. Nat. Rev. Microbiol. 11 , 21-32, doi:10.1038/nrmicro2916 (2013). Lind, A. L., Lim, F. Y., Soukup, A. A., Keller, N. P. & Rokas, A. An LaeA- and BrlA-Dependent Cellular Network Governs Tissue-Specific Secondary Metabolism in the Human Pathogen Aspergillus fumigatus. mSphere 3 , doi:10.1128/mSphere.00050-18 (2018). Park, H. S. & Yu, J. H. Developmental regulators in Aspergillus fumigatus. J. Microbiol. 54 , 223-231, doi:10.1007/s12275-016-5619-5 (2016). Son, Y. E., Yu, J. H. & Park, H. S. Regulators of the Asexual Life Cycle of Aspergillus nidulans. Cells 12 , doi:10.3390/cells12111544 (2023). Kim, S. S., Kim, Y. H. & Shin, K. S. The developmental regulators, FlbB and FlbE, are involved in the virulence of Aspergillus fumigatus. J. Microbiol. Biotechnol. 23 , 766-770, doi:10.4014/jmb.1301.01002 (2013). Bertuzzi, M. et al. On the lineage of Aspergillus fumigatus isolates in common laboratory use. Med. Mycol. 59 , 7-13, doi:10.1093/mmy/myaa075 (2021). 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. 77 , 1203-1219, doi:10.1111/j.1365-2958.2010.07282.x (2010). Cho, H. J. et al. Regulation of Conidiogenesis in Aspergillus flavus. Cells 11 , doi:10.3390/cells11182796 (2022). Li, Y. et al. Comprehensive analysis of the regulatory network of blue-light-regulated conidiation and hydrophobin production in Trichoderma guizhouense. Environ. Microbiol. 23 , 6241-6256, doi:10.1111/1462-2920.15748 (2021). 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. 23 , 5184-5199, doi:10.1111/1462-2920.15500 (2021). Tanaka, M. et al. The C2H2-type transcription factor, FlbC, is involved in the transcriptional regulation of Aspergillus oryzae glucoamylase and protease genes specifically expressed in solid-state culture. Appl. Microbiol. Biotechnol. 100 , 5859-5868, doi:10.1007/s00253-016-7419-6 (2016). Son, H., Kim, M. G., Chae, S. K. & Lee, Y. W. FgFlbD regulates hyphal differentiation required for sexual and asexual reproduction in the ascomycete fungus Fusarium graminearum. J. Microbiol. 52 , 930-939, doi:10.1007/s12275-014-4384-6 (2014). Fukuda, S. et al. Trade-off between Plasticity and Velocity in Mycelial Growth. mBio 12 , doi:10.1128/mBio.03196-20 (2021). 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 9 , 1711-1723, doi:10.1128/EC.00198-10 (2010). Moon, H. et al. Unraveling the Gene Regulatory Networks of the Global Regulators VeA and LaeA in Aspergillus nidulans. Microbiol Spectr 11 , e0016623, doi:10.1128/spectrum.00166-23 (2023). Rocha, M. C. et al. The Cell Wall Integrity Pathway Contributes to the Early Stages of Aspergillus fumigatus Asexual Development. Appl. Environ. Microbiol. 86 , doi:10.1128/AEM.02347-19 (2020). Stewart, J. I. P. et al. Reducing Aspergillus fumigatus Virulence through Targeted Dysregulation of the Conidiation Pathway. mBio 11 , doi:10.1128/mBio.03202-19 (2020). Mah, J. H. & Yu, J. H. Upstream and downstream regulation of asexual development in Aspergillus fumigatus. Eukaryot Cell 5 , 1585-1595, doi:10.1128/EC.00192-06 (2006). Ries, L. N. A. et al. Nutritional Heterogeneity Among Aspergillus fumigatus Strains Has Consequences for Virulence in a Strain- and Host-Dependent Manner. Front. Microbiol. 10 , 854, doi:10.3389/fmicb.2019.00854 (2019). Todokoro, T. et al. Identification of a novel pyrithiamine resistance marker gene thiI for genome co-editing in Aspergillus oryzae. J. Biosci. Bioeng. 130 , 227-232, doi:10.1016/j.jbiosc.2020.04.013 (2020). Onoue, T. et al. Identification of Two Mannosyltransferases Contributing to Biosynthesis of the Fungal-type Galactomannan alpha-Core-Mannan Structure in Aspergillus fumigatus. Sci. Rep. 8 , 16918, doi:10.1038/s41598-018-35059-2 (2018). Arai, T. et al. Hmg1 mutations in Aspergillus fumigatus and their contribution to triazole susceptibility. Med. Mycol. 59 , 980-984, doi:10.1093/mmy/myab026 (2021). Kadooka, C. et al. Identification of an alpha-(1-->6)-Mannosyltransferase Contributing To Biosynthesis of the Fungal-Type Galactomannan alpha-Core-Mannan Structure in Aspergillus fumigatus. mSphere 7 , e0048422, doi:10.1128/msphere.00484-22 (2022). Additional Declarations No competing interests reported. Supplementary Files SupplementaryFigS1.pdf SupplementaryTableS1S2StrainlistandPrimerlist.pdf SupplementaryTableS3Dataavailability.xlsx Cite Share Download PDF Status: Posted Version 1 posted You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-4435109","acceptedTermsAndConditions":true,"allowDirectSubmit":true,"archivedVersions":[],"articleType":"Article","associatedPublications":[],"authors":[{"id":309576180,"identity":"5e585631-c558-4518-a04a-ce2734389f4a","order_by":0,"name":"Yutaka Tanaka","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAABA0lEQVRIiWNgGAWjYPCCBAZ+BgY2BoYDCCFm3KqZIVok20jWYnAMTQtOYHAj/+Dngpo0eeP7zccefDhzmEF3do8Bw48aBnZznFqSmaVnHMsx3HaMLd1wxo3DDGZ3zhgw9hxjYLZswKmFQZqHrYJx2zEeM2meD4frt93IMWDgbWBgNsDhSpAtv3n+VdhvboNoYTADamH8i18LmzRvW07iBjaQlhsQLcz4bJE889jMmrcvLXnGsbQ0yRln0oF+OVZwWOaYBE6/8B1PfHyb51uybX/z4WMSH45ZM5jdbt748E2NTTKuEFPAtF0CHD0SyQY4tMhj2i4BoexwaRkFo2AUjIIRBwCdoVtuaefSVgAAAABJRU5ErkJggg==","orcid":"","institution":"Tohoku Medical and Pharmaceutical University","correspondingAuthor":true,"prefix":"","firstName":"Yutaka","middleName":"","lastName":"Tanaka","suffix":""},{"id":309576181,"identity":"dc7f96b1-a2e0-4a65-be60-6174a109597b","order_by":1,"name":"Ikumi Asagiri","email":"","orcid":"","institution":"Tohoku Medical and Pharmaceutical University","correspondingAuthor":false,"prefix":"","firstName":"Ikumi","middleName":"","lastName":"Asagiri","suffix":""},{"id":309576182,"identity":"0a232731-e2e0-4682-8275-fb165c5e4852","order_by":2,"name":"Rei Komiya","email":"","orcid":"","institution":"Tohoku Medical and Pharmaceutical University","correspondingAuthor":false,"prefix":"","firstName":"Rei","middleName":"","lastName":"Komiya","suffix":""},{"id":309576183,"identity":"57f29e4d-a75f-432a-bc00-b928a9bfae16","order_by":3,"name":"Hiroyuki Yamazaki","email":"","orcid":"","institution":"Tohoku Medical and Pharmaceutical University","correspondingAuthor":false,"prefix":"","firstName":"Hiroyuki","middleName":"","lastName":"Yamazaki","suffix":""},{"id":309576184,"identity":"610c04bc-e482-4b6c-9f58-f4a63471b5d1","order_by":4,"name":"Akiho Yagi","email":"","orcid":"","institution":"Tohoku Medical and Pharmaceutical University","correspondingAuthor":false,"prefix":"","firstName":"Akiho","middleName":"","lastName":"Yagi","suffix":""},{"id":309576185,"identity":"b3c516bd-9557-4ebe-9a91-7a96ee41e996","order_by":5,"name":"Ryuji Uchida","email":"","orcid":"","institution":"Tohoku Medical and Pharmaceutical University","correspondingAuthor":false,"prefix":"","firstName":"Ryuji","middleName":"","lastName":"Uchida","suffix":""},{"id":309576186,"identity":"71ee29b3-ec96-44a1-9a65-edbcaafcd627","order_by":6,"name":"Masato Sasaki","email":"","orcid":"","institution":"Tohoku Medical and Pharmaceutical University","correspondingAuthor":false,"prefix":"","firstName":"Masato","middleName":"","lastName":"Sasaki","suffix":""}],"badges":[],"createdAt":"2024-05-17 07:54:58","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-4435109/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-4435109/v1","draftVersion":[],"editorialEvents":[],"editorialNote":"","failedWorkflow":false,"files":[{"id":57715048,"identity":"2704d00e-3ecf-44af-9d9c-c4a846d7ff28","added_by":"auto","created_at":"2024-06-04 16:52:05","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":17633,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eExpression level of the \u003c/strong\u003e\u003cem\u003e\u003cstrong\u003eflbC\u003c/strong\u003e\u003c/em\u003e\u003cstrong\u003e gene peaks almost in line, but the ratio varies between strains.\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe expression of \u003cem\u003eflbC\u003c/em\u003e during growth was determined using quantitative real-time PCR. The cultures were incubated in a liquid minimal medium (MM) until each time point. The expression ratios were normalized at 8 h, and the expression level of actin gene was the internal control, according to the ΔΔCt method. Data are the mean ± standard deviation of three independent experiments (n=3).\u003c/p\u003e","description":"","filename":"Fig1ExpressionleveloftheflbCgenepeaksalmostinlinebuttheratiovariesbetweenstrains..png","url":"https://assets-eu.researchsquare.com/files/rs-4435109/v1/98a76b2478f30f3f1247128b.png"},{"id":57715049,"identity":"b4abd207-9741-4186-9a38-8fecc7cea0c3","added_by":"auto","created_at":"2024-06-04 16:52:05","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":231072,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eMycelial growth of the Δ\u003c/strong\u003e\u003cem\u003e\u003cstrong\u003eflbC\u003c/strong\u003e\u003c/em\u003e\u003cstrong\u003e strains varies between strains.\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003e(a)\u003cstrong\u003e \u003c/strong\u003eSchematic representation of \u003cem\u003eflbC \u003c/em\u003egene deletion via homologous recombination and RNP-based CRISPR-Cas9 system. Black box indicates CDS of \u003cem\u003eflbC\u003c/em\u003e gene; white box indicates pyrithiamine resistant gene (\u003cem\u003eptrA\u003c/em\u003e), used as a selective marker; and x bp indicates the sequence length in the specified range. The scissors icon (marked as CRISPR#1, 2) indicates the PAM sequences used in this study. The arrows indicate the positions of the primers for colony PCR. (b) Confirmation of the correct recombination of \u003cem\u003eflbC \u003c/em\u003egene via colony PCR analysis. (Marker); DNA size marker. (c, d) Radial growth phenotypes of each WT and ∆\u003cem\u003eflbC\u003c/em\u003estrains grown on a solid (c) MM or (d) MM with 400 μM CongoRed for 5 d at 37 °C. Solid lines with black dots indicate WT strains, and dashed lines with white dots indicate Δ\u003cem\u003eflbC\u003c/em\u003e strains. The data are representative of five independent experiments.\u003c/p\u003e","description":"","filename":"Fig2MycelialgrowthoftheflbCstrainsvariesbetweenstrains..png","url":"https://assets-eu.researchsquare.com/files/rs-4435109/v1/2c4dbd8aa489884ccee957e8.png"},{"id":57715055,"identity":"12caf29e-d1d4-4064-af18-29468b969421","added_by":"auto","created_at":"2024-06-04 16:52:05","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":20205,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eIntra-species heterogeneity of the \u003c/strong\u003e\u003cem\u003e\u003cstrong\u003ebrlA\u003c/strong\u003e\u003c/em\u003e\u003cstrong\u003e gene expression observed in WT strains is abolished by \u003c/strong\u003e\u003cem\u003e\u003cstrong\u003eflbC\u003c/strong\u003e\u003c/em\u003e\u003cstrong\u003e gene deletion.\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eExpression of the \u003cem\u003ebrlA\u003c/em\u003e gene during growth was determined by using quantitative real-time PCR. The cultures were incubated in a liquid MM until each time point. The expression ratios were normalized at 4 h, and the expression level of actin gene was the internal control, according to the ΔΔCt method. Data are the mean ± standard deviation of three independent experiments and compared among groups at each time point using 2-way ANOVA with Tukey’s test for multiple comparisons (n =3): *, \u003cem\u003ep\u003c/em\u003e \u0026lt; 0.05; n.s., not significant. The data are representative of two independent experiments.\u003c/p\u003e","description":"","filename":"Fig3IntraspeciesheterogeneityofthebrlAgeneexpressionobservedinWTstrainsisabolishedbyflbCgenedeletion..png","url":"https://assets-eu.researchsquare.com/files/rs-4435109/v1/1e5574e2b4aa02de0dd33cd4.png"},{"id":57715345,"identity":"05bb6eb7-cecc-4fdb-aea6-51d48b601b74","added_by":"auto","created_at":"2024-06-04 17:00:06","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":65544,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eΔ\u003c/strong\u003e\u003cem\u003e\u003cstrong\u003eflbC\u003c/strong\u003e\u003c/em\u003e\u003cstrong\u003e strains show intra-species heterogeneity in fumagillin production capacity.\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe production of secondary metabolites including fumagillin produced by the Δ\u003cem\u003eflbC \u003c/em\u003estrains (a) and wild-type strains (b) and an authentic fumagillin (c) were detected and compared by HPLC analysis.\u003c/p\u003e","description":"","filename":"Fig4flbCstrainsshowintraspeciesheterogeneityinfumagillinproductioncapacity.1.png","url":"https://assets-eu.researchsquare.com/files/rs-4435109/v1/8c4e631ecdf34e36c77c0bbd.png"},{"id":57715053,"identity":"4f356600-2fda-4f99-af04-01dd3798b040","added_by":"auto","created_at":"2024-06-04 16:52:05","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":49589,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eΔ\u003c/strong\u003e\u003cem\u003e\u003cstrong\u003eflbC\u003c/strong\u003e\u003c/em\u003e\u003cstrong\u003e strains show intra-species heterogeneity in transcriptional activity of the genes forming the fumagillin-pseurotin gene cluster.\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eExpression of the genes forming the fumagillin gene cluster that (a) contribute and (b) do not contribute to fumagillin production at 24 h of growth was determined using quantitative real-time PCR. The cultures were incubated in a liquid MM until each time point. The expression ratios were normalized at 4 h, and the expression level of actin gene was the internal control, according to the ΔΔCt method. Data are the mean ± standard deviation of three independent experiments. The indicated \u003cem\u003ep-\u003c/em\u003evalues were determined using Student’s \u003cem\u003et\u003c/em\u003e-test (n = 3): *, \u003cem\u003ep\u003c/em\u003e\u0026lt; 0.05; **, \u003cem\u003ep\u003c/em\u003e \u0026lt; 0.01; n.s., not significant. The data are representative of two independent experiments.\u003c/p\u003e","description":"","filename":"Fig5flbCstrainsshowintraspeciesheterogeneityintranscriptionalactivityofthegenesformingthefumagillinpseurotingenecluster..png","url":"https://assets-eu.researchsquare.com/files/rs-4435109/v1/892f4fa50f6166e6be4aef5b.png"},{"id":57715056,"identity":"f82ed251-ee55-4cad-9aa6-d026083de99b","added_by":"auto","created_at":"2024-06-04 16:52:05","extension":"png","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":31114,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eΔ\u003c/strong\u003e\u003cem\u003e\u003cstrong\u003eflbC\u003c/strong\u003e\u003c/em\u003e\u003cstrong\u003e strains show intra-species heterogeneity in the timing and total expression of \u003c/strong\u003e\u003cem\u003e\u003cstrong\u003elaeA\u003c/strong\u003e\u003c/em\u003e\u003cstrong\u003e and \u003c/strong\u003e\u003cem\u003e\u003cstrong\u003efapR\u003c/strong\u003e\u003c/em\u003e\u003cstrong\u003e.\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe expression of (a) \u003cem\u003elaeA\u003c/em\u003e and (b) \u003cem\u003efapR\u003c/em\u003e genes during growth was determined using quantitative real-time PCR. The cultures were incubated in a liquid MM until each time point. The expression ratios were normalized at 4 h, and the expression level of actin gene was the internal control, according to the ΔΔCt method. Data are the mean ± standard deviation of three independent experiments and compared among groups at each time point using 2-way ANOVA with Tukey’s test for multiple comparisons (n=3 for A1151, n=6 for A1280) ; * \u003cem\u003ep\u003c/em\u003e \u0026lt; 0.05; n.s., not significant. The data are representative of two independent experiments.\u003c/p\u003e","description":"","filename":"Fig6flbCstrainsshowintraspeciesheterogeneityinthetimingandtotalexpressionoflaeAandfapR..png","url":"https://assets-eu.researchsquare.com/files/rs-4435109/v1/cde9fef3aa2a634fe6baa350.png"},{"id":57715058,"identity":"cb68a6e7-9ca9-4ddc-8ce1-307ae2da969b","added_by":"auto","created_at":"2024-06-04 16:52:05","extension":"png","order_by":7,"title":"Figure 7","display":"","copyAsset":false,"role":"figure","size":10889,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eThe \u003c/strong\u003e\u003cem\u003e\u003cstrong\u003eflbC\u003c/strong\u003e\u003c/em\u003e\u003cstrong\u003e gene in any strains is dispensable for infecting silkworms.\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe survival of silkworms (\u003cem\u003eB. mori\u003c/em\u003e) infected with \u003cem\u003eA. fumigatus\u003c/em\u003e is shown as Kaplan–Meier curves (\u003cem\u003en\u003c/em\u003e = 10 larvae/strain). Solid lines indicate groups infected with WT, and dashed lines indicate groups infected with Δ\u003cem\u003eflbC\u003c/em\u003e. The indicated \u003cem\u003ep-\u003c/em\u003evalues were determined using Wilcoxon test comparing each Δ\u003cem\u003eflbC\u003c/em\u003e with WT.\u003c/p\u003e","description":"","filename":"Fig7TheflbCgeneinanystrainsisdispensableforinfectingsilkworms..png","url":"https://assets-eu.researchsquare.com/files/rs-4435109/v1/b424dbe2ad28d102fb9db5be.png"},{"id":58916798,"identity":"c450e922-abdb-46ff-9328-380e65f32a6f","added_by":"auto","created_at":"2024-06-24 06:02:40","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":1329760,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-4435109/v1/20948c73-71ad-4d99-a8b4-2fa46b4786c0.pdf"},{"id":57715050,"identity":"cab13955-c0aa-45d9-9414-718d3305c6bf","added_by":"auto","created_at":"2024-06-04 16:52:05","extension":"pdf","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":1250034,"visible":true,"origin":"","legend":"","description":"","filename":"SupplementaryFigS1.pdf","url":"https://assets-eu.researchsquare.com/files/rs-4435109/v1/4e466cad5d352abbdf4f0c64.pdf"},{"id":57715947,"identity":"7624b044-970e-4d25-a7b7-9473bab71d1b","added_by":"auto","created_at":"2024-06-04 17:08:06","extension":"pdf","order_by":2,"title":"","display":"","copyAsset":false,"role":"supplement","size":115300,"visible":true,"origin":"","legend":"","description":"","filename":"SupplementaryTableS1S2StrainlistandPrimerlist.pdf","url":"https://assets-eu.researchsquare.com/files/rs-4435109/v1/5529a8fb606da6816c3c9ebc.pdf"},{"id":57715052,"identity":"364ba1c4-f6f7-48b6-854b-b80e63ee8d2e","added_by":"auto","created_at":"2024-06-04 16:52:05","extension":"xlsx","order_by":3,"title":"","display":"","copyAsset":false,"role":"supplement","size":74567,"visible":true,"origin":"","legend":"","description":"","filename":"SupplementaryTableS3Dataavailability.xlsx","url":"https://assets-eu.researchsquare.com/files/rs-4435109/v1/763eabc18865e25c9a687d82.xlsx"}],"financialInterests":"No competing interests reported.","formattedTitle":"The FlbC transcription factor contributes to the generation of strain heterogeneity in Fumagillin mycotoxin production in Aspergillus fumigatus","fulltext":[{"header":"Introduction","content":"\u003cp\u003e \u003cem\u003eAspergillus fumigatus\u003c/em\u003e is an opportunistic pathogen, which infects human beings. Although \u003cem\u003eA. fumigatus\u003c/em\u003e infection in healthy individuals has not established, persons with a history of chemotherapy, immunosuppressive drug use, or congenital immunodeficiency syndrome can have lethal lesions\u003csup\u003e\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e\u003c/sup\u003e. The human-pathogenicity characteristic of \u003cem\u003eA. fumigatus\u003c/em\u003e is determined by the interaction of several virulence factors in the organism. The mycelia produce secreted proteins\u003csup\u003e\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e\u003c/sup\u003e, polysaccharides\u003csup\u003e\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e\u003c/sup\u003e, and a wide variety of secondary metabolites\u003csup\u003e\u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e4\u003c/span\u003e\u003c/sup\u003e. In recent years, inter-species, inter-strain, and even individual spore heterogeneity in these phenotypes and genotypes has been noted\u003csup\u003e\u003cspan additionalcitationids=\"CR6 CR7 CR8 CR9\" citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e10\u003c/span\u003e\u003c/sup\u003e. Heterogeneity has also been observed between clinical and environmental isolates\u003csup\u003e\u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e11\u003c/span\u003e\u003c/sup\u003e. In basic research, several strains with different backgrounds are currently used as reference strains; however, notable differences have been reported between these strains in terms of pathogenicity\u003csup\u003e\u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e\u003c/sup\u003e, genome structure\u003csup\u003e\u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e,\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e\u003c/sup\u003e, secondary-metabolite production capacity\u003csup\u003e\u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e14\u003c/span\u003e\u003c/sup\u003e, single nucleotide polymorphisms (SNPs) \u003csup\u003e\u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e15\u003c/span\u003e\u003c/sup\u003e, chromatin modifications\u003csup\u003e\u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e\u003c/sup\u003e, and transposable element distribution\u003csup\u003e\u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e\u003c/sup\u003e. Vigorous research is now underway to understand how this heterogeneity arises between species or strains, mainly through large-scale data. Af293 is the most commonly used \u003cem\u003eA. fumigatus\u003c/em\u003e strain for phenotyping in basic research because it was the first strain to have its entire genome sequenced and is currently the only strain that has been allocated to the chromosome\u003csup\u003e\u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e\u003c/sup\u003e. In addition, CEA17\u003csup\u003e18\u003c/sup\u003e, D141\u003csup\u003e19\u003c/sup\u003e, and ATCC46645\u003csup\u003e20\u003c/sup\u003e strains have been established and used as laboratory strains. Heterogeneity among strains has made it more complicated to understand the physiology of the species by comparison between each study.\u003c/p\u003e \u003cp\u003eAlthough the secondary metabolites produced by \u003cem\u003eAspergillus\u003c/em\u003e spp. are not explicitly associated with pathogenicity in humans, a correlation between the repertoire of secondary metabolites that can be produced and pathogenicity has been reported\u003csup\u003e\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e\u003c/sup\u003e. Inter-species and inter-strain heterogeneity is also observed in the repertoire of secondary metabolites\u003csup\u003e\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e\u003c/sup\u003e. Many \u003cem\u003eA. fumigatus\u003c/em\u003e strains that secrete fumagillin, pseurotin, fumitremorgin, and gliotoxin are pathogenic, and closely related species of \u003cem\u003eA. fumigatus\u003c/em\u003e, \u003cem\u003eA. oerlinghausenensis\u003c/em\u003e and \u003cem\u003eA. fischeri\u003c/em\u003e can produce fumagillin and pseurotin as secondary metabolites; however, these organisms are not pathogenic to humans. Comparisons between strains with different repertoires of secondary metabolite secretion have shown that strains with a more severe clinical course tend to secrete a greater variety of metabolites. Although no single secondary metabolite is known to determine the human pathogenicity of \u003cem\u003eA. fumigatus\u003c/em\u003e, the relationship between strain heterogeneity in secondary-metabolite production repertoire and pathogenicity remains unclear.\u003c/p\u003e \u003cp\u003eThe \u003cem\u003efma\u003c/em\u003e gene family\u003csup\u003e\u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e,\u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e23\u003c/span\u003e\u003c/sup\u003e, required for fumagillin biosynthesis in \u003cem\u003eA. fumigatus\u003c/em\u003e, forms a cluster in the subtelomeric region of chromosome 8\u003csup\u003e24\u003c/sup\u003e. The gene cluster has eight Fma enzymes that catalyze sequential reactions to biosynthesize the end product and the transcription factor FapR/FumR\u003csup\u003e24,25\u003c/sup\u003e, which positively regulates the expression of the \u003cem\u003efma\u003c/em\u003e gene family upstream. Because fumagillin precursors are secreted by strains lacking each gene of the \u003cem\u003efma\u003c/em\u003e gene family, each \u003cem\u003efma\u003c/em\u003e gene is essential for the biosynthesis of fumagillin\u003csup\u003e\u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e23\u003c/span\u003e\u003c/sup\u003e. Loss of the gene encoding the histone methyltransferase LaeA affects the production of fumagillin and other secondary metabolites\u003csup\u003e\u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e26\u003c/span\u003e,\u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e27\u003c/span\u003e\u003c/sup\u003e. Inside the fumagillin gene cluster, there is an intertwined group of genes involved in pseurotin biosynthesis, which is synchronously regulated by FapR and LaeA\u003csup\u003e\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e\u003c/sup\u003e. In addition, a gene cluster for the biosynthesis of fumitremorgin is present in the neighborhood\u003csup\u003e\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e\u003c/sup\u003e, suggesting that fumagillin-pseurotin-fumitremorgin is synchronously affected by the regulation of chromatin structure. Chromosome 8 has a lower frequency of trimethylation modifications of the ninth lysine residue of histone H3 than other genomic regions and has many LINE transposon sequences\u003csup\u003e\u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e\u003c/sup\u003e, suggesting that the area around the fumagillin gene cluster may be a region of relatively high genomic instability. Consistent with this, approximately 60% of \u003cem\u003eA. fumigatus\u003c/em\u003e isolates have lost a part of or almost all the fumagillin gene cluster in the subtelomeric region of chromosome 8\u003csup\u003e15\u003c/sup\u003e. Thus, the role of the \u003cem\u003efma\u003c/em\u003e gene cluster and its upstream transcription factors in the production of fumagillin, a typical secondary metabolite, is becoming clearer; however, the mechanism by which this heterogeneity arises in the fumagillin production is not yet known.\u003c/p\u003e \u003cp\u003eSecondary metabolite production is often regulated synchronously between several secondary-metabolite-producing systems rather than independently of each other\u003csup\u003e\u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e\u003c/sup\u003e. Lind et al. reported that the BrlA transcription factor, which is central to the central developmental pathway, affects the expression of almost all secondary metabolite products\u003csup\u003e\u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e\u003c/sup\u003e. The central developmental pathway is based on sequential transcriptional activation of a series of transcription factors, namely BrlA, AbaA, and WetA, which are important for the transcription of genes involved in morphogenesis\u003csup\u003e\u003cspan citationid=\"CR30\" class=\"CitationRef\"\u003e30\u003c/span\u003e\u003c/sup\u003e. During the asexual life cycle of \u003cem\u003eA. fumigatus\u003c/em\u003e, dormant spores first swell and germinate to form filamentous mycelia, followed by the asexual reproductive spores. These series of shifts are synchronized with the continuous activation of the central developmental pathway. The comparison of the mRNA expression profiles of \u003cem\u003elaeA-\u003c/em\u003e and \u003cem\u003ebrlA-\u003c/em\u003egene deficient strains and the secondary-metabolic production profiles have shown that despite the high correlation between the two profiles and the fumagillin-fumitremorgin-pseurotin gene cluster, there is negative correlation between the expression patterns of the these gene clusters\u003csup\u003e\u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e\u003c/sup\u003e. The central developmental pathway regulated by BrlA is further regulated upstream by several Flb transcription factors, which are together termed the upstream developmental activator (UDA) pathway. In \u003cem\u003eAspergilus nidulans\u003c/em\u003e, FlbB, FlbC, and FlbE are transcription factors that constitute the UDA pathway\u003csup\u003e\u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e31\u003c/span\u003e\u003c/sup\u003e. FlbB and FlbE cooperate, while FlbC alone activates the transcription of \u003cem\u003ebrlA\u003c/em\u003e\u003csup\u003e\u003cem\u003e\u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e31\u003c/span\u003e\u003c/em\u003e\u003c/sup\u003e. In \u003cem\u003eA. fumigatus\u003c/em\u003e strains deficient in \u003cem\u003eflbB\u003c/em\u003e and \u003cem\u003eflbE\u003c/em\u003e genes, reduced \u003cem\u003ebrlA\u003c/em\u003e gene expression, delayed mycelial growth, and reduced fumagillin and gliotoxin production have been observed\u003csup\u003e\u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e32\u003c/span\u003e\u003c/sup\u003e. However, the effects of the \u003cem\u003eflbC\u003c/em\u003e gene on \u003cem\u003eA. fumigatus\u003c/em\u003e growth, \u003cem\u003ebrlA\u003c/em\u003e gene expression, and secondary metabolite production have not yet been investigated in detail. In addition, in the aforementioned comparisons of \u003cem\u003ebrlA\u003c/em\u003e, \u003cem\u003elaeA\u003c/em\u003e, and \u003cem\u003eflbB/E\u003c/em\u003e gene-deficient strains, experiments were conducted using a single reference strain for each strain and inter-strain heterogeneity was not considered.\u003c/p\u003e \u003cp\u003eTherefore, this study aimed to investigate the effect of the \u003cem\u003eflbC\u003c/em\u003e gene on fumagillin production in \u003cem\u003eA. fumigatus\u003c/em\u003e by introducing \u003cem\u003eflbC\u003c/em\u003e deletions into several laboratory strains from different backgrounds and to consider the effect of the \u003cem\u003eflbC\u003c/em\u003e gene on inter-strain heterogeneity. In particular, changes in the expression pattern of the \u003cem\u003ebrlA\u003c/em\u003e gene over time, fumagillin production capacity, expression of the \u003cem\u003efma\u003c/em\u003e gene cluster, and expression of the upstream transcription factors \u003cem\u003elaeA\u003c/em\u003e and \u003cem\u003efapR\u003c/em\u003e were compared.\u003c/p\u003e"},{"header":"Results","content":"\u003cp\u003e \u003cb\u003eExpression patterns of the\u003c/b\u003e \u003cb\u003eflbC\u003c/b\u003e \u003cb\u003egene in the three strains\u003c/b\u003e\u003c/p\u003e \u003cp\u003eA1151, A1159, and A1280 as laboratory strains with suppressed non-homologous repair for genetic modification were deposited in FGSC. The strains CEA17, D141, and ATCC46645 were the parental strains of A1151, A1159, and A1280, respectively\u003csup\u003e\u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e33\u003c/span\u003e\u003c/sup\u003e. The \u003cem\u003eflbC\u003c/em\u003e gene in strain Af293 (\u003cem\u003eAfu2g13770\u003c/em\u003e) was searched in the draft genomes of strains CEA17 and found \u003cem\u003eAFUB029400-T\u003c/em\u003e. Using this sequence as a reference, primers were designed to amplify the entire CDS length of the \u003cem\u003eflbC\u003c/em\u003e gene. Amplification was attempted using the genomic DNA of the three strains as a template, and amplification products were obtained in all groups. The sequences in the CDS of the amplified products were checked for point mutations and were found to be completely homologous in all three strains, suggesting that there is probably no difference in the amino acid sequence of the FlbC protein produced in the three strains. Next, to compare the expression patterns of the \u003cem\u003eflbC\u003c/em\u003e gene in the three strains, RNA was extracted from spores and mycelia obtained at different incubation times, and the relative amounts of \u003cem\u003eflbC\u003c/em\u003e gene transcripts were quantified over time. The \u003cem\u003eflbC\u003c/em\u003e gene transcript was not amplified immediately after the start of incubation, when almost all spores formed buds, but was detected after approximately 8 h and reached peak expression in all three strains after 16 h (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e). Thereafter, a gradual downward trend was observed in A1151 and A1280, followed by a slight increase in A1159 at 72 h. These trends are consistent with previous research, in which the \u003cem\u003eflbC\u003c/em\u003e transcript levels in \u003cem\u003eA. nidulans\u003c/em\u003e during vegetative growth peaked at approximately 18\u0026ndash;24 h and decreased thereafter\u003csup\u003e\u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e34\u003c/span\u003e\u003c/sup\u003e. However, some heterogeneity appeared among the three strains, with expression ratios of approximately 3-fold in A1151, 4-fold in A1159, and 9-fold in A1280 at 16 h. In summary, the \u003cem\u003eflbC\u003c/em\u003e gene appears to reach peak expression level approximately 16 h after the start of culture, that is, at the same time as peak mycelial growth. Comparison of three representative laboratory strains suggested that the CDS internal sequence of the \u003cem\u003eflbC\u003c/em\u003e gene and the timing of the maximum expression level were generally the same, but the range of expression changes was slightly different. These results suggest that the differences among strains in the range of expression changes are dependent on differences in upstream expression of \u003cem\u003eflbC\u003c/em\u003e.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003cb\u003eEffect of\u003c/b\u003e \u003cb\u003eflbC\u003c/b\u003e \u003cb\u003eexpression heterogeneity on strain phenotype\u003c/b\u003e\u003c/p\u003e \u003cp\u003eTo investigate the effect of heterogeneity in the \u003cem\u003eflbC\u003c/em\u003e expression pattern on the phenotype of the three laboratory strains, we produced \u003cem\u003eflbC\u003c/em\u003e gene disruption mutant strains (Δ\u003cem\u003eflbC\u003c/em\u003e) in each background. Because the internal sequence of \u003cem\u003eflbC\u003c/em\u003e was exactly the same, we planned highly efficient gene disruption using the CRISPR-Cas9 method with the PAM sequence common to the three strains. The Δ\u003cem\u003eflbCs\u003c/em\u003e were generated by inducing double-strand breaks via Cas9 (CRISPR#1, 2) binding, followed by insertion of pyrithiamine resistance gene markers through homologous recombination, using each of the three strains as parent (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ea, b). Δ\u003cem\u003eflbCs\u003c/em\u003e were slightly reduced at radial growth circle than the parental strains (WT) on MM at 37\u0026deg;C (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ec). The three WT strains grew differentially under these conditions, with the A1280 strain growing slightly weaker on day 3 or earlier than the other two. The radial growth of the \u003cem\u003eflbC\u003c/em\u003e gene-deficient strain was almost identical to this trend. Notably, when grown with the cell wall effector Congo red, A1280 showed weaker growth among the WT strains; however, growth was restored by \u003cem\u003eflbC\u003c/em\u003e gene deletion (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ed). This was in contrast to A1151 and A1159, whose growth was suppressed by Congo red treatment. These findings suggest that the differences in the growth phenotypes observed in the three laboratory strains may be due to the function of the \u003cem\u003eflbC\u003c/em\u003e gene.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003cb\u003eEffect of\u003c/b\u003e \u003cb\u003eflbC\u003c/b\u003e \u003cb\u003eexpression patterns on\u003c/b\u003e \u003cb\u003ebrlA\u003c/b\u003e \u003cb\u003eexpression\u003c/b\u003e\u003c/p\u003e \u003cp\u003eBrlA, which plays a central- and entry point- role in the central developmental pathway, is regulated by FlbC transcription factor\u003csup\u003e\u003cspan citationid=\"CR30\" class=\"CitationRef\"\u003e30\u003c/span\u003e,\u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e34\u003c/span\u003e\u003c/sup\u003e. To investigate how different expression patterns of the \u003cem\u003eflbC\u003c/em\u003e gene affect the \u003cem\u003ebrlA\u003c/em\u003e expression, we measured \u003cem\u003ebrlA\u003c/em\u003e transcript over time in WT and Δ\u003cem\u003eflbC\u003c/em\u003e of the three laboratory strains. WT of the A1151 and A1159-background reached peak mRNA expression at 16\u0026ndash;24 h, with an approximate 3- to 4-fold increase than the expression at 4 h, whereas this increase was suppressed in the Δ\u003cem\u003eflbC\u003c/em\u003e (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e, upper and middle panel). In contrast, the WT of A1280-background showed little increase in \u003cem\u003ebrlA\u003c/em\u003e gene expression until 48 h after the start of incubation, and Δ\u003cem\u003eflbC\u003c/em\u003e in the A1280-background also showed little change in \u003cem\u003ebrlA\u003c/em\u003e gene expression (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e, bottom panel). In summary, \u003cem\u003ebrlA\u003c/em\u003e gene expression level in WT of A1151 and A1159 peaked at 16\u0026ndash;24 h, whereas this peak was removed by \u003cem\u003eflbC\u003c/em\u003e gene deletion. In contrast, A1280 showed a flat expression up to 48 h in \u003cem\u003ebrlA\u003c/em\u003e gene expression in both WT and Δ\u003cem\u003eflbC\u003c/em\u003e rather than a peak-like pattern. These results suggest that the \u003cem\u003ebrlA\u003c/em\u003e gene expression pattern is different among the three laboratory strains and that there is heterogeneity in the regulatory role of the \u003cem\u003eflbC\u003c/em\u003e gene on \u003cem\u003ebrlA\u003c/em\u003e gene expression.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003cb\u003eEffect of\u003c/b\u003e \u003cb\u003eflbC\u003c/b\u003e \u003cb\u003egene on the heterogeneity of secondary metabolite production\u003c/b\u003e\u003c/p\u003e \u003cp\u003eTo investigate the effect of the \u003cem\u003eflbC\u003c/em\u003e gene on the heterogeneity of secondary metabolite production, we measured the fumagillin-production capacity of these strains. Remarkably, given that previous studies have shown the loss of fumagillin-producing capacity in several \u003cem\u003eA. fumigatus\u003c/em\u003e strains\u003csup\u003e\u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e15\u003c/span\u003e\u003c/sup\u003e, all three WT strains produced fumagillin (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e). The effect of \u003cem\u003eflbC\u003c/em\u003e gene deletion on fumagillin production was tripartite: in A1151, fumagillin production was below the detection limit due to \u003cem\u003eflbC\u003c/em\u003e gene deletion; in A1159, fumagillin production was reduced but not eliminated due to \u003cem\u003eflbC\u003c/em\u003e gene deletion; and A1280 produced sufficient amount of fumagillin in both WT and Δ\u003cem\u003eflbC\u003c/em\u003e, with little difference in production level (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e). These results indicate that although all three laboratory strains are capable of producing fumagillin, there is heterogeneity in the effects of the \u003cem\u003eflbC\u003c/em\u003e gene deletion. To investigate whether the heterogeneity in producing fumagillin is dependent on differences in the expression levels of the genes responsible for fumagillin biosynthesis, we measured the transcripts belonging to \u003cem\u003efma\u003c/em\u003e gene family involved in the fumagillin gene cluster (\u003cem\u003efma-TC / fmaA\u003c/em\u003e, \u003cem\u003efma-P450 / fmaG\u003c/em\u003e, \u003cem\u003efma-C6H / fmaF\u003c/em\u003e, \u003cem\u003efma-AT / fmaC\u003c/em\u003e, and \u003cem\u003efma-ABM / fmaE\u003c/em\u003e) and those not involved in fumagillin biosynthesis but located within the gene cluster (\u003cem\u003eAfu8g00430\u003c/em\u003e, \u003cem\u003eAfu8g00440 / psoF\u003c/em\u003e and \u003cem\u003eAfu8g00500\u003c/em\u003e). In A1151-background, the expression levels of many fumagillin biosynthesis-related genes, including \u003cem\u003efma-TC / fmaA\u003c/em\u003e, \u003cem\u003efma-P450 / fmaG\u003c/em\u003e, \u003cem\u003efma-C6H / fmaF\u003c/em\u003e, and \u003cem\u003efma-AT / fmaC\u003c/em\u003e, were significantly reduced in Δ\u003cem\u003eflbC\u003c/em\u003e than those in WT (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003ea). Some genes not involved in fumagillin biosynthesis were also downregulated (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003eb). A comparison of the A1159 background strains showed the same trend, except for \u003cem\u003efma-AT / fmaC\u003c/em\u003e (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003ea). In contrast, comparisons of A1280 background strains showed significant differences in the reduced expression of \u003cem\u003efma-TC / fmaA\u003c/em\u003e and \u003cem\u003efma-C6H / fmaF\u003c/em\u003e, but the extent of the reduction was weaker than that in A1151 and A1159, and differences in the expression of other \u003cem\u003efma\u003c/em\u003e genes was not significant (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003ea). In summary, in the Δ\u003cem\u003eflbC\u003c/em\u003e of A1151 and A1159 backgrounds, the expression of the \u003cem\u003efma\u003c/em\u003e gene family was generally and coordinately suppressed, consistent with the reduced fumagillin production, while in A1280, where \u003cem\u003eflbC\u003c/em\u003e gene deletion did not result in reduced fumagillin production, the degree of suppression of \u003cem\u003efma\u003c/em\u003e gene family expression was weaker than in the other two strain. Thus, the heterogeneity in the three laboratory strains between the WT and Δ\u003cem\u003eflbC\u003c/em\u003e with regard to fumagillin production may be dependent on differences in the expression levels of the \u003cem\u003efma\u003c/em\u003e gene family and probably correlates with the degree of hetero-chromatinization of the fumagillin gene cluster.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003cb\u003eExpression of the upstream transcription factors\u003c/b\u003e \u003cb\u003elaeA\u003c/b\u003e \u003cb\u003eand\u003c/b\u003e \u003cb\u003efapR\u003c/b\u003e\u003c/p\u003e \u003cp\u003eThe production of secondary metabolites requires the coordinated action of global and cluster-specific transcription factors. In the case of fumagillin, the transcriptional activation of biosynthetic genes is achieved by the histone methyltransferase LaeA, which acts as a global transcription factor, and FapR, which acts as a cluster-specific transcription factor. The \u003cem\u003efapR\u003c/em\u003e gene is located within the fumagillin gene cluster, and production of the FapR transcription factor is initiated by the release of heterochromatin from the cluster by LaeA. To determine whether there are differences in the upstream gene of the \u003cem\u003efma\u003c/em\u003e gene family, we measured the expression levels of \u003cem\u003elaeA\u003c/em\u003e and \u003cem\u003efapR\u003c/em\u003e genes in the A1151 and A1280 backgrounds between WTs and Δ\u003cem\u003eflbC\u003c/em\u003es, which had contrasting difference in the fumagillin production. In the A1151-background strains, \u003cem\u003elaeA\u003c/em\u003e gene expression peaked after 16\u0026ndash;24 h and then decreased in the WT, whereas the peak was delayed to 48\u0026ndash;72 h in the Δ\u003cem\u003eflbC\u003c/em\u003e (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003ea, left panel). In contrast, in A1280, the expression of the \u003cem\u003elaeA\u003c/em\u003e gene peaked at 16 h in the WT as same as that in A1151, while the Δ\u003cem\u003eflbC\u003c/em\u003e showed a similar pattern to that of the WT (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003ea, right panel). The \u003cem\u003efapR\u003c/em\u003e gene showed reduced gene expression at 16\u0026ndash;24 h in the Δ\u003cem\u003eflbC\u003c/em\u003e than that in the WT of the A1151 background (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003eb, left panel), while there was no significant difference in expression levels at 24 h in the A1280 background (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003eb, right panel). In summary, in the A1151 background strain, the \u003cem\u003eflbC\u003c/em\u003e gene deletion resulted in delayed expression of the \u003cem\u003elaeA\u003c/em\u003e gene and suppressed the expression of \u003cem\u003efapR\u003c/em\u003e, but not in the A1280 background. These results indicate that the heterogeneity in fumagillin production between strains A1151 and A1280 may be due to the different levels of transcriptional activation of FapR and LaeA. More specifically, this suggests a dependence on different direct or indirect interactions between \u003cem\u003eflbC\u003c/em\u003e and genes regulating the expression of the fumagillin gene cluster, such as \u003cem\u003elaeA\u003c/em\u003e and \u003cem\u003efapR\u003c/em\u003e.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cdiv id=\"Sec3\" class=\"Section2\"\u003e \u003ch2\u003eEffect of heterogeneities on virulence\u003c/h2\u003e \u003cp\u003eThe analysis of the three laboratory strains shows that the deletion of the \u003cem\u003eflbC\u003c/em\u003e gene may create new heterogeneity or reduce that which was present originally. To estimate the effects of these heterogeneities on virulence, a pathogenicity test was conducted on silkworms. We found that the deletion of the \u003cem\u003eflbC\u003c/em\u003e gene probably does not contribute to silkworm-pathogenicity of the strains because there was no significant difference in the survival among WT and Δ\u003cem\u003eflbC\u003c/em\u003e of any background (Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003e). Deletion of the \u003cem\u003eflbC\u003c/em\u003e gene reduced the growth of the laboratory strains on MM and showed differences in the production of the mycotoxin fumagillin; however, these differences were not observed in the pathogenicity test.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e"},{"header":"Discussion","content":"\u003cp\u003eThis study shows that the \u003cem\u003eflbC\u003c/em\u003e gene in \u003cem\u003eA. fumigatus\u003c/em\u003e may increase or decrease the heterogeneity of hyphal growth, the expression patterns of \u003cem\u003ebrlA\u003c/em\u003e transcripts and secondary metabolites. FlbC is a C\u003csub\u003e2\u003c/sub\u003eH\u003csub\u003e2\u003c/sub\u003e transcription factor that is widely conserved in Ascomycetes, but not in Basidiomycetes and Zygomycetes\u003csup\u003e\u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e34\u003c/span\u003e\u003c/sup\u003e. The \u003cem\u003eflbC\u003c/em\u003e gene has been found in \u003cem\u003eA. nidulans\u003c/em\u003e\u003csup\u003e\u003cem\u003e\u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e34\u003c/span\u003e\u003c/em\u003e\u003c/sup\u003e, \u003cem\u003eA. flavus\u003c/em\u003e\u003csup\u003e\u003cem\u003e\u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e35\u003c/span\u003e\u003c/em\u003e\u003c/sup\u003e, \u003cem\u003eTrichoderma guizhouense\u003c/em\u003e\u003csup\u003e\u003cem\u003e\u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e36\u003c/span\u003e\u003c/em\u003e\u003c/sup\u003e, and \u003cem\u003eBeauveria bassiana\u003c/em\u003e\u003csup\u003e\u003cem\u003e\u003cspan citationid=\"CR37\" class=\"CitationRef\"\u003e37\u003c/span\u003e\u003c/em\u003e\u003c/sup\u003e. In \u003cem\u003eA. oryzae\u003c/em\u003e, it regulates the production of amylolytic enzymes, such as glucoamylase and α-amylase, during individual media growth\u003csup\u003e\u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e38\u003c/span\u003e\u003c/sup\u003e. \u003cem\u003eFusarium graminearum\u003c/em\u003e deficient of \u003cem\u003eflbC\u003c/em\u003e gene showed no phenotype\u003csup\u003e\u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e39\u003c/span\u003e\u003c/sup\u003e. FlbC regulates \u003cem\u003ebrlA\u003c/em\u003e gene expression, along with the concerted action of FlbB, FlbD, and FlbE\u003csup\u003e\u003cspan citationid=\"CR30\" class=\"CitationRef\"\u003e30\u003c/span\u003e\u003c/sup\u003e. Thus, the extent to which any \u003cem\u003eflb\u003c/em\u003e family molecule contributes to the asexual life cycle varies among species. These Flbs also have differnt localizations in fungal cells, with FlbB and FlbE tending to accumulate at the tips of filamentous fungi, while FlbC is present in all nuclei, including metulae and phialides\u003csup\u003e\u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e34\u003c/span\u003e\u003c/sup\u003e. This is not surprising, given that the plasticity and velocity of filamentous growth, as well as the cell wall and cell surface structures, differ between species\u003csup\u003e\u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e40\u003c/span\u003e\u003c/sup\u003e. Genomic heterogeneity between clinical and environmental isolates is known, and transcriptome heterogeneity exists within the same strain, even within individual spores\u003csup\u003e\u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e,\u003cspan additionalcitationids=\"CR10\" citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e11\u003c/span\u003e\u003c/sup\u003e. The reason for the various levels of heterogeneity arising in \u003cem\u003eA. fumigatus\u003c/em\u003e will become increasingly clear in future studies.\u003c/p\u003e \u003cp\u003eAs we showed, in the A1151 and A1159 backgrounds, loss of the \u003cem\u003eflbC\u003c/em\u003e gene was associated with reduced fumagillin production. These strains showed reduced expression of most members of the \u003cem\u003efma\u003c/em\u003e gene family and, in correlation, delayed the expression of \u003cem\u003elaeA\u003c/em\u003e gene and suppressed the expression of \u003cem\u003efapR\u003c/em\u003e gene in A1151 cells. The effect of \u003cem\u003elaeA\u003c/em\u003e and \u003cem\u003efapR\u003c/em\u003e gene expression on fumagillin production may be explained by the following facts presented by Dhingra et al. \u003csup\u003e\u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e25\u003c/span\u003e\u003c/sup\u003e. The \u003cem\u003efapR/fumR\u003c/em\u003e gene-deficient strain completely loses its ability to produce fumagillin, and the expression of the \u003cem\u003efma\u003c/em\u003e gene family is also almost completely suppressed. The expression levels of \u003cem\u003efapR/fumR\u003c/em\u003e and \u003cem\u003efmaB\u003c/em\u003e are suppressed by more than 90% in the \u003cem\u003elaeA\u003c/em\u003e gene-deficient strain, and \u003cem\u003efmaB\u003c/em\u003e gene expression is also suppressed in the \u003cem\u003efapR\u003c/em\u003e gene-deficient strain. In strains deficient in the VeA transcription factor, which forms a heterotrimer with the LaeA protein, \u003cem\u003efapR\u003c/em\u003e gene expression is repressed, and fumagillin production is reduced by more than 80%. Thus the expression of the \u003cem\u003efma\u003c/em\u003e gene family; global transcription factors such as LaeA and VeA, which are upregulated; and FapR, a fumagillin cluster-specific transcription factor, play important roles in fumagillin production. However, the mechanism by which \u003cem\u003eflbC\u003c/em\u003e gene deletion causes delayed expression of the \u003cem\u003elaeA\u003c/em\u003e gene remains unclear. Several studies have suggested a link between central developmental pathways and secondary metabolic production. Lind et al. reported altered production patterns of many secondary metabolic products in \u003cem\u003ebrlA\u003c/em\u003e and \u003cem\u003elaeA\u003c/em\u003e gene deletion strains\u003csup\u003e\u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e\u003c/sup\u003e. Kim et al. have also shown that deletion of the \u003cem\u003eflbB\u003c/em\u003e and \u003cem\u003eflbE\u003c/em\u003e genes results in repression of the \u003cem\u003elaeA\u003c/em\u003e gene and reduced production of gliotoxin and fumagillin\u003csup\u003e\u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e32\u003c/span\u003e\u003c/sup\u003e. However, they did not discuss the mechanisms of \u003cem\u003elaeA\u003c/em\u003e repression, but showed an altered expression pattern of the \u003cem\u003ebrlA\u003c/em\u003e gene due to \u003cem\u003eflb\u003c/em\u003e gene deletion\u003csup\u003e\u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e41\u003c/span\u003e\u003c/sup\u003e. Moon et al. showed by multi-omics analysis that the \u003cem\u003eflbC\u003c/em\u003e gene may be a direct target of LaeA and VeA transcription factors\u003csup\u003e\u003cspan citationid=\"CR42\" class=\"CitationRef\"\u003e42\u003c/span\u003e\u003c/sup\u003e. Rocha et al. revealed the regulation of \u003cem\u003eflbC\u003c/em\u003e gene expression by the MAPK transcription factor RlmA\u003csup\u003e\u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e43\u003c/span\u003e\u003c/sup\u003e. A feedback regulatory mechanism may be present between FlbC and global transcription factors such as LaeA and RlmA. As we have shown, the Congo red sensitivity of the A1280 WT strain was restored by the loss of \u003cem\u003eflbC\u003c/em\u003e gene, which suggests a link between the cell wall integrity pathway and \u003cem\u003eflbC\u003c/em\u003e.\u003c/p\u003e \u003cp\u003eLoss of the \u003cem\u003ebrlA\u003c/em\u003e gene, which is thought to be located downstream of FlbC, results in increased production of fumagillin\u003csup\u003e\u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e\u003c/sup\u003e. This seems to contradict our finding that \u003cem\u003eflbC\u003c/em\u003e gene deletion reduced fumagillin production in A1151 and A1159. This may be due to perturbations in LaeA or FapR expression. However, we may have to consider the compensatory effects induced by gene loss (or overexpression), as is the case when certain important genes are involved in the basic properties and survival of mycelia. Considering that BrlA is an important transcription factor that plays a critical role in central developmental pathways, and that it may be important for \u003cem\u003ebrlA\u003c/em\u003e expression to be oscillatory, unidirectionally increasing or decreasing its expression may have drastic effects\u003csup\u003e\u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e44\u003c/span\u003e,\u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e45\u003c/span\u003e\u003c/sup\u003e. This, in turn, is also true for \u003cem\u003eflbC\u003c/em\u003e, and is a refutation and claim of the present study. The fact that fumagillin was reduced in \u003cem\u003eflbC\u003c/em\u003e-deficient strains may not necessarily indicate that the FlbC transcription factor is \u0026ldquo;directly\u0026rdquo; involved in fumagillin production.\u003c/p\u003e \u003cp\u003eNotably, among the three standard laboratory strains, heterogeneity was observed not only in the gene-deficient strains but also in the WT strain. A1151, A1159, and A1280 were clinical isolates, but each isolate was isolated from a person with a different disease condition\u003csup\u003e\u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e33\u003c/span\u003e\u003c/sup\u003e. Comparisons of the three strains have already been made at various levels, including whole genome, transcriptome, metabolome, leukemia model mouse infectivity, and macrophage phagocytosis capacity\u003csup\u003e\u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e,\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e,\u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e,\u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e46\u003c/span\u003e\u003c/sup\u003e. The fumagillin gene cluster is located in the subtelomeric region of chromosome 8, where transposable elements and chromatin modifications are frequently found in Af293 and CEA17 backgrounds\u003csup\u003e\u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e,\u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e\u003c/sup\u003e. Perrin et al. focused on the genomic instability of secondary-metabolite gene clusters and suggested that the expression profiles of secondary metabolites may differ among strains, showing that the structure of the subtelomeric region varies among strains, even within the same species\u003csup\u003e\u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e27\u003c/span\u003e\u003c/sup\u003e. Thereafter, Lind et al. organized these mechanisms and divided them into different types\u003csup\u003e\u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e\u003c/sup\u003e. The present study just demonstrates this hypothetical model. However, the heterogeneity in fumagillin production identified in this study is probably not due to genomic instability in the subtelomeric region of chromosome 8 because WT of A1151, A1159, and A1280 expressed \u003cem\u003efma\u003c/em\u003e genes and produced fumagillin. According to Lind et al., the fumagillin gene cluster is relatively stable because it has fewer SNPs with non-synonymous substitutions than the other genomic regions\u003csup\u003e\u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e\u003c/sup\u003e. FlbC of the three strains behaves slightly differently with regard to the regulation of the fumagillin gene cluster; however, understanding this difference will probably require more than a comparison of genomic DNA between the strains. Nevertheless, with the advent of the CRISPR-Cas9, as we have used, recombinants of various strains from one species can be obtained efficiently and quickly. Although we obtained mutants by homologous recombination of drug resistance markers, \u0026ldquo;co-genome editing\u003csup\u003e\u003cspan citationid=\"CR47\" class=\"CitationRef\"\u003e47\u003c/span\u003e\u003c/sup\u003e,\u0026rdquo; which does not use homologous recombination of markers, is now being developed. This enables the genome editing of a large number of clinical and environmental isolates and the use of these isolates as parental strains. Comparison of more strains than that in this study, which compared only three strains, is now possible.\u003c/p\u003e \u003cp\u003eThus, our findings show that the different interactions between \u003cem\u003eflbC\u003c/em\u003e and factors (LaeA and FapR) regulating the expression of fumagillin gene cluster are the probable cause for heterogeneity in fumagillin production among the strains. Moreover, our study serves as a basis for future studies on co-genome editing of large number of clinical and environmental isolates and the use of these isolates as parental strains.\u003c/p\u003e"},{"header":"Methods","content":"\u003cp\u003eStrains, media, and culture conditions\u003c/p\u003e \u003cp\u003eThe \u003cem\u003eA. fumigatus\u003c/em\u003e strains used in this study are listed in Supplementary Table \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003e. \u003cem\u003eA. fumigatus\u003c/em\u003e strains A1151, A1159, and A1280 were obtained from the Fungal Genetics Stock Center (FGSC). The strains were grown on minimal medium [MM; 1% (w/v) glucose, 0.6% (w/v) NaNO\u003csub\u003e3\u003c/sub\u003e, 0.052% (w/v) KCl, 0.052% (w/v) MgSO\u003csub\u003e4\u003c/sub\u003e\u0026middot;7 H\u003csub\u003e2\u003c/sub\u003eO, 0.152% (w/v) KH\u003csub\u003e2\u003c/sub\u003ePO\u003csub\u003e4\u003c/sub\u003e, and Hunter's trace elements, pH 6.5]. Plasmids were amplified in \u003cem\u003eE. coli\u003c/em\u003e DH5α. Colony growth rates were measured as previously described\u003csup\u003e\u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e48\u003c/span\u003e\u003c/sup\u003e. Briefly, conidia from each strain were point-inoculated onto the center of MM agar plates with or without 400 \u0026micro;g/mL CongoRed (Fujifilm, Japan). Colony diameters were measured after 3, 4, and 5 d of incubation at 37\u0026deg;C. The growth rates of all the individual strains were measured 5 times.\u003c/p\u003e \u003cp\u003eGeneration of \u003cem\u003eflbC\u003c/em\u003e mutants\u003c/p\u003e \u003cp\u003eEach sgRNA target sequence was explored using the web service CHOPCHOP (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://chopchop.cbu.uib.no/\u003c/span\u003e\u003cspan address=\"https://chopchop.cbu.uib.no/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e) and two synthesized DNA oligonucleotides, primer No. 3 and 4 (Supplementary Table \u003cspan refid=\"MOESM2\" class=\"InternalRef\"\u003eS2\u003c/span\u003e). These sequences contained the T7 promoter, an sgRNA target sequence, and a sequence overlapping the Cas9 scaffold. The synthesis and purification of sgRNAs were performed using the Guide-it sgRNA In Vitro Transcription Kit (TAKARA, Japan). After the synthesized sgRNAs (sgRNA-1 and \u0026minus;\u0026thinsp;2) were quantified using a NanoDrop spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA), they were used for ribonucleoprotein (RNP) formation with Cas9.\u003c/p\u003e \u003cp\u003eRepair templates were constructed based on the method described by Arai et al., with slight modifications\u003csup\u003e\u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e49\u003c/span\u003e\u003c/sup\u003e. The pyrithiamine resistance cassette (\u003cem\u003eptrA\u003c/em\u003e marker) was obtained from the pPTRII plasmid (TAKARA). The KOD One PCR Master Mix (TOYOBO, Japan) and primer No. 1 and 2 (Supplementary Table \u003cspan refid=\"MOESM2\" class=\"InternalRef\"\u003eS2\u003c/span\u003e) were used for PCR amplification to construct the repair template. The PCR products were purified using a FastGene Gel/PCR Extraction kit (NIPPON Genetics Co., Ltd., Japan) and quantified using a NanoDrop spectrophotometer.\u003c/p\u003e \u003cp\u003e \u003cem\u003eA. fumigatus\u003c/em\u003e protoplasts were generated, and fungal transformation was performed as described by Kadooka et al., with slight modifications\u003csup\u003e\u003cspan citationid=\"CR50\" class=\"CitationRef\"\u003e50\u003c/span\u003e\u003c/sup\u003e. Briefly, conidia were incubated in YG mediem (2% glucose, 0.5% yeast extract, and Hunter\u0026rsquo;s trace element) for 16 h at 37\u0026deg;C. Following incubation, the cell walls of germlings were digested with 32 mg/mL Vino-Taste Pro (Novozymes, Denmark) for 4 h 30 min at 30\u0026deg;C. Subsequently, 10 \u0026micro;g Guide-it Recombinant Cas9 (TAKARA) and 10 pmol each of sgRNA-1 and \u0026minus;\u0026thinsp;2 were mixed and incubated for 5 min, generating RNPs as described by Arai et al..\u003csup\u003e\u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e49\u003c/span\u003e\u003c/sup\u003e Protoplasts were transformed with 2\u0026ndash;3 \u0026micro;g of repair templates and RNPs and plated onto MM agar supplemented with 0.6 M KCl. Following 24 h incubation at room temperature, agar plates were overlaid with MM top agar containing 0.1 \u0026micro;g/mL pyrithiamine hydrobromide (Sigma-Aldrich, USA). The \u003cem\u003eflbC\u003c/em\u003e mutants were isolated and confirmed using colony PCR with primer No. 5 and 6 (Supplementary Table \u003cspan refid=\"MOESM2\" class=\"InternalRef\"\u003eS2\u003c/span\u003e). At least three deletion strains were isolated in each case.\u003c/p\u003e \u003cp\u003eQuantitative real-time PCR analysis\u003c/p\u003e \u003cp\u003eQuantification of mRNA expression was performed using THUNDERBIRD SYBR qPCR mix (TOYOBO) and StepOnePlus Real-Time PCR system (Thermo Fisher Scientific). The strains were grown at 37\u0026deg;C in liquid MM at each time point. Fungal RNA was isolated using the Sepasol-RNA I Super G (Nacalai Tesque, Japan) and reverse-transcribed into cDNA using 1 \u0026micro;g total RNA and ReverTra Ace qPCR RT Kit (TOYOBO). RT-PCR was performed using primers No. 7\u0026ndash;32 (Supplementary Table \u003cspan refid=\"MOESM2\" class=\"InternalRef\"\u003eS2\u003c/span\u003e), and gene expression levels were calculated using ΔΔCt method.\u003c/p\u003e \u003cp\u003eHPLC analysis\u003c/p\u003e \u003cp\u003eEach strain, maintained on a PDA plate, was inoculated into a 100-mL Erlenmeyer flask containing 50 mL of the production medium (3.0% sucrose, 3.0% soluble starch, 1.0% malt extract, 0.30% Ebios, 0.50% KH\u003csub\u003e2\u003c/sub\u003ePO\u003csub\u003e4\u003c/sub\u003e, and 0.050% MgSO\u003csub\u003e4\u003c/sub\u003e\u0026middot;7H\u003csub\u003e2\u003c/sub\u003eO; adjusted to pH 6.0 before sterilization). The production culture was carried out at 25\u0026deg;C for 7 days under agitation (150 rpm). Each culture broth was treated with 50 mL of acetone after 7 days and filtered. The filtrate, after the evaporation of acetone, was extracted with EtOAc (50 mL), and the extract was concentrated \u003cem\u003ein vacuo\u003c/em\u003e.\u003c/p\u003e \u003cp\u003eThe production of secondary metabolites was monitored using a Chromaster HPLC system (Hitachi High Technologies Co., Ltd., Tokyo, Japan) [column, CAPCELL PAK C18 MG-II (Osaka Soda. Co., Ltd., Osaka, Japan), i.d. 4.6 mm \u0026times; 250 mm; flow rate, 0.8 mL/min; mobile phase, linear gradient from 20 to 100% CH\u003csub\u003e3\u003c/sub\u003eCN containing 0.05% TFA for 30 min and 100% CH\u003csub\u003e3\u003c/sub\u003eCN containing 0.05% TFA for 5 min; detection, UV 365 nm]. Each extract was prepared to 10 mg/mL CH\u003csub\u003e3\u003c/sub\u003eOH solution, and 5.0 \u0026micro;L of each sample was injected for the HPLC analysis.\u003c/p\u003e \u003cp\u003eSilkworm larva infection model\u003c/p\u003e \u003cp\u003eSilkworm (\u003cem\u003eBombyx mori\u003c/em\u003e) larvae were maintained in an incubator at 27\u0026deg;C in the dark at a controlled temperature until conidial infection. \u003cem\u003eA. fumigatus\u003c/em\u003e conidia were resuspended at a concentration of 1 \u0026times; 10\u003csup\u003e8\u003c/sup\u003e conidia/mL in phosphate buffered saline (PBS), and then silkworm 5th-instar larvae (2.0 g, n\u0026thinsp;=\u0026thinsp;10) were infected with 1 \u0026times; 10\u003csup\u003e6\u003c/sup\u003e conidia (50 \u0026micro;L/larva) by injecting the hemolymph using a disposable 1 mL syringe with a 27G needle (TERUMO, Tokyo, Japan). The uninfected controls received PBS. Larvae were raised in the dark at 37\u0026deg;C for 4 d. Generalized Wilcoxon tests were performed using GraphPad Prism 8 software to compare survival curves, and differences were considered statistically significant at \u003cem\u003ep\u003c/em\u003e\u0026thinsp;\u0026lt;\u0026thinsp;0.05. Fertilized silkworm eggs of Bombyx mori (Hu\u0026middot;Yo \u0026times; Tukuba\u0026middot;Ne) were purchased from Ehime Sansyu (Ehime, Japan) and fed an artificial diet (Silk Mate 2S; Nihon Nosan Kogyo, Kanagawa, Japan, and Silkmate; Katakura Industries, Tokyo, Japan).\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eAcknowledgments\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003e\u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; This work was supported in part by Grant-in-Aid for Scientific Research (C) (22K06600 to Y.T., , 21K06631 and 24K09865 to H.Y.) from the Japan Society for the Promotion of Science (JSPS KAKENHI). We express our thanks to Ms. M. Iwamoto of Tohoku Medical and Pharmaceutical University for her technical assistance. We would like to thank Editage (www.editage.jp) for English language editing.\u003c/p\u003e\n\u003cp\u003e\u0026nbsp;\u003cstrong\u003eAuthor contributions\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eY.T. conceived and designed the experiments. Y.T., I.A, R.K., H.Y. and A.Y. conducted the experiments. R.U. and M.S. helped with the design and coordination of the study. Y.T., H.Y. and A. Y. wrote the paper. All authors reviewed the manuscript.\u003c/p\u003e\n\u003cp\u003e\u0026nbsp;\u003cstrong\u003eCompeting interests\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors declare no competing interests.\u003c/p\u003e\n\u003cp\u003e\u0026nbsp;\u003cstrong\u003eData availability statement\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eAll data generated or analyzed during this study are included in Supplementary Table S3. A full-length version of the figure obtained by agarose gel electrophoresis is shown in Supplementary FigS1, and no contrast adjustment was made when used in the figures.\u003c/p\u003e\n\u003cp\u003e\u0026nbsp;\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n\u003cli\u003eFisher, M. C. \u0026amp; Denning, D. W. The WHO fungal priority pathogens list as a game-changer. \u003cem\u003eNat. Rev. Microbiol.\u003c/em\u003e \u003cstrong\u003e21\u003c/strong\u003e, 211-212, doi:10.1038/s41579-023-00861-x (2023).\u003c/li\u003e\n\u003cli\u003eWartenberg, D.\u003cem\u003e et al.\u003c/em\u003e Secretome analysis of Aspergillus fumigatus reveals Asp-hemolysin as a major secreted protein. \u003cem\u003eInt. J. Med. Microbiol.\u003c/em\u003e \u003cstrong\u003e301\u003c/strong\u003e, 602-611, doi:10.1016/j.ijmm.2011.04.016 (2011).\u003c/li\u003e\n\u003cli\u003eBriard, B.\u003cem\u003e et al.\u003c/em\u003e Galactosaminogalactan activates the inflammasome to provide host protection. \u003cem\u003eNature\u003c/em\u003e \u003cstrong\u003e588\u003c/strong\u003e, 688-692, doi:10.1038/s41586-020-2996-z (2020).\u003c/li\u003e\n\u003cli\u003eManfiolli, A. O.\u003cem\u003e et al.\u003c/em\u003e Mitogen-Activated Protein Kinase Cross-Talk Interaction Modulates the Production of Melanins in Aspergillus fumigatus. \u003cem\u003emBio\u003c/em\u003e \u003cstrong\u003e10\u003c/strong\u003e, doi:10.1128/mBio.00215-19 (2019).\u003c/li\u003e\n\u003cli\u003eDanion, F.\u003cem\u003e et al.\u003c/em\u003e Aspergillus fumigatus, One Uninucleate Species with Disparate Offspring. \u003cem\u003eJ Fungi (Basel)\u003c/em\u003e \u003cstrong\u003e7\u003c/strong\u003e, doi:10.3390/jof7010030 (2021).\u003c/li\u003e\n\u003cli\u003eHewitt, S. K., Foster, D. S., Dyer, P. S. \u0026amp; Avery, S. V. Phenotypic heterogeneity in fungi: Importance and methodology. \u003cem\u003eFungal Biol. Rev.\u003c/em\u003e \u003cstrong\u003e30\u003c/strong\u003e, 176-184, doi:10.1016/j.fbr.2016.09.002 (2016).\u003c/li\u003e\n\u003cli\u003eIjadpanahsaravi, M., Teertstra, W. R. \u0026amp; Wosten, H. A. B. Inter- and intra-species heterogeneity in germination of Aspergillus conidia. \u003cem\u003eAntonie Van Leeuwenhoek\u003c/em\u003e \u003cstrong\u003e115\u003c/strong\u003e, 1151-1164, doi:10.1007/s10482-022-01762-4 (2022).\u003c/li\u003e\n\u003cli\u003eKeller, N. P. Heterogeneity Confounds Establishment of \u0026quot;a\u0026quot; Model Microbial Strain. \u003cem\u003emBio\u003c/em\u003e \u003cstrong\u003e8\u003c/strong\u003e, doi:10.1128/mBio.00135-17 (2017).\u003c/li\u003e\n\u003cli\u003eWang, F.\u003cem\u003e et al.\u003c/em\u003e Transcription in fungal conidia before dormancy produces phenotypically variable conidia that maximize survival in different environments. \u003cem\u003eNat Microbiol\u003c/em\u003e \u003cstrong\u003e6\u003c/strong\u003e, 1066-1081, doi:10.1038/s41564-021-00922-y (2021).\u003c/li\u003e\n\u003cli\u003eDos Santos, R. A. C.\u003cem\u003e et al.\u003c/em\u003e Genomic and Phenotypic Heterogeneity of Clinical Isolates of the Human Pathogens Aspergillus fumigatus, Aspergillus lentulus, and Aspergillus fumigatiaffinis. \u003cem\u003eFront Genet\u003c/em\u003e \u003cstrong\u003e11\u003c/strong\u003e, 459, doi:10.3389/fgene.2020.00459 (2020).\u003c/li\u003e\n\u003cli\u003eHorta, M. A. C.\u003cem\u003e et al.\u003c/em\u003e Examination of Genome-Wide Ortholog Variation in Clinical and Environmental Isolates of the Fungal Pathogen Aspergillus fumigatus. \u003cem\u003emBio\u003c/em\u003e \u003cstrong\u003e13\u003c/strong\u003e, e0151922, doi:10.1128/mbio.01519-22 (2022).\u003c/li\u003e\n\u003cli\u003eKowalski, C. H.\u003cem\u003e et al.\u003c/em\u003e Heterogeneity among Isolates Reveals that Fitness in Low Oxygen Correlates with Aspergillus fumigatus Virulence. \u003cem\u003emBio\u003c/em\u003e \u003cstrong\u003e7\u003c/strong\u003e, doi:10.1128/mBio.01515-16 (2016).\u003c/li\u003e\n\u003cli\u003eGarcia-Rubio, R., Monzon, S., Alcazar-Fuoli, L., Cuesta, I. \u0026amp; Mellado, E. Genome-Wide Comparative Analysis of Aspergillus fumigatus Strains: The Reference Genome as a Matter of Concern. \u003cem\u003eGenes (Basel)\u003c/em\u003e \u003cstrong\u003e9\u003c/strong\u003e, doi:10.3390/genes9070363 (2018).\u003c/li\u003e\n\u003cli\u003eLind, A. L.\u003cem\u003e et al.\u003c/em\u003e Drivers of genetic diversity in secondary metabolic gene clusters within a fungal species. \u003cem\u003ePLoS Biol.\u003c/em\u003e \u003cstrong\u003e15\u003c/strong\u003e, e2003583, doi:10.1371/journal.pbio.2003583 (2017).\u003c/li\u003e\n\u003cli\u003eBarber, A. E.\u003cem\u003e et al.\u003c/em\u003e Aspergillus fumigatus pan-genome analysis identifies genetic variants associated with human infection. \u003cem\u003eNat Microbiol\u003c/em\u003e \u003cstrong\u003e6\u003c/strong\u003e, 1526-1536, doi:10.1038/s41564-021-00993-x (2021).\u003c/li\u003e\n\u003cli\u003eColabardini, A. C.\u003cem\u003e et al.\u003c/em\u003e Chromatin profiling reveals heterogeneity in clinical isolates of the human pathogen Aspergillus fumigatus. \u003cem\u003ePLoS Genet.\u003c/em\u003e \u003cstrong\u003e18\u003c/strong\u003e, e1010001, doi:10.1371/journal.pgen.1010001 (2022).\u003c/li\u003e\n\u003cli\u003eFedorova, N. D.\u003cem\u003e et al.\u003c/em\u003e Genomic islands in the pathogenic filamentous fungus Aspergillus fumigatus. \u003cem\u003ePLoS Genet.\u003c/em\u003e \u003cstrong\u003e4\u003c/strong\u003e, e1000046, doi:10.1371/journal.pgen.1000046 (2008).\u003c/li\u003e\n\u003cli\u003eKress, M. R. V. Z.\u003cem\u003e et al.\u003c/em\u003e The akuBKU80 Mutant Deficient for Nonhomologous End Joining Is a Powerful Tool for Analyzing Pathogenicity in Aspergillus fumigatus. \u003cem\u003eEUKARYOTIC CELL\u003c/em\u003e \u003cstrong\u003e5\u003c/strong\u003e, 207-211, doi:10.1128/EC.5.1.207 (2006).\u003c/li\u003e\n\u003cli\u003eWagener, J.\u003cem\u003e et al.\u003c/em\u003e The putative alpha-1,2-mannosyltransferase AfMnt1 of the opportunistic fungal pathogen Aspergillus fumigatus is required for cell wall stability and full virulence. \u003cem\u003eEukaryot Cell\u003c/em\u003e \u003cstrong\u003e7\u003c/strong\u003e, 1661-1673, doi:10.1128/EC.00221-08 (2008).\u003c/li\u003e\n\u003cli\u003eHartmann, T.\u003cem\u003e et al.\u003c/em\u003e Validation of a self-excising marker in the human pathogen Aspergillus fumigatus by employing the beta-rec/six site-specific recombination system. \u003cem\u003eAppl. Environ. Microbiol.\u003c/em\u003e \u003cstrong\u003e76\u003c/strong\u003e, 6313-6317, doi:10.1128/AEM.00882-10 (2010).\u003c/li\u003e\n\u003cli\u003eSteenwyk, J. L.\u003cem\u003e et al.\u003c/em\u003e Variation Among Biosynthetic Gene Clusters, Secondary Metabolite Profiles, and Cards of Virulence Across Aspergillus Species. \u003cem\u003eGenetics\u003c/em\u003e \u003cstrong\u003e216\u003c/strong\u003e, 481-497, doi:10.1534/genetics.120.303549 (2020).\u003c/li\u003e\n\u003cli\u003eLin, H. C.\u003cem\u003e et al.\u003c/em\u003e The fumagillin biosynthetic gene cluster in Aspergillus fumigatus encodes a cryptic terpene cyclase involved in the formation of beta-trans-bergamotene. \u003cem\u003eJ. Am. Chem. Soc.\u003c/em\u003e \u003cstrong\u003e135\u003c/strong\u003e, 4616-4619, doi:10.1021/ja312503y (2013).\u003c/li\u003e\n\u003cli\u003eLin, H. C.\u003cem\u003e et al.\u003c/em\u003e Generation of complexity in fungal terpene biosynthesis: discovery of a multifunctional cytochrome P450 in the fumagillin pathway. \u003cem\u003eJ. Am. Chem. Soc.\u003c/em\u003e \u003cstrong\u003e136\u003c/strong\u003e, 4426-4436, doi:10.1021/ja500881e (2014).\u003c/li\u003e\n\u003cli\u003eWiemann, P.\u003cem\u003e et al.\u003c/em\u003e Prototype of an intertwined secondary-metabolite supercluster. \u003cem\u003eProc. Natl. Acad. Sci. U. S. A.\u003c/em\u003e \u003cstrong\u003e110\u003c/strong\u003e, 17065-17070, doi:10.1073/pnas.1313258110 (2013).\u003c/li\u003e\n\u003cli\u003eDhingra, S.\u003cem\u003e et al.\u003c/em\u003e The fumagillin gene cluster, an example of hundreds of genes under veA control in Aspergillus fumigatus. \u003cem\u003ePLoS One\u003c/em\u003e \u003cstrong\u003e8\u003c/strong\u003e, e77147, doi:10.1371/journal.pone.0077147 (2013).\u003c/li\u003e\n\u003cli\u003eKeller, N. P., Turner, G. \u0026amp; Bennett, J. W. Fungal secondary metabolism - from biochemistry to genomics. \u003cem\u003eNat. Rev. Microbiol.\u003c/em\u003e \u003cstrong\u003e3\u003c/strong\u003e, 937-947, doi:10.1038/nrmicro1286 (2005).\u003c/li\u003e\n\u003cli\u003ePerrin, R. M.\u003cem\u003e et al.\u003c/em\u003e Transcriptional regulation of chemical diversity in Aspergillus fumigatus by LaeA. \u003cem\u003ePLoS Pathog.\u003c/em\u003e \u003cstrong\u003e3\u003c/strong\u003e, e50, doi:10.1371/journal.ppat.0030050 (2007).\u003c/li\u003e\n\u003cli\u003eBrakhage, A. A. Regulation of fungal secondary metabolism. \u003cem\u003eNat. Rev. Microbiol.\u003c/em\u003e \u003cstrong\u003e11\u003c/strong\u003e, 21-32, doi:10.1038/nrmicro2916 (2013).\u003c/li\u003e\n\u003cli\u003eLind, A. L., Lim, F. Y., Soukup, A. A., Keller, N. P. \u0026amp; Rokas, A. An LaeA- and BrlA-Dependent Cellular Network Governs Tissue-Specific Secondary Metabolism in the Human Pathogen Aspergillus fumigatus. \u003cem\u003emSphere\u003c/em\u003e \u003cstrong\u003e3\u003c/strong\u003e, doi:10.1128/mSphere.00050-18 (2018).\u003c/li\u003e\n\u003cli\u003ePark, H. S. \u0026amp; Yu, J. H. Developmental regulators in Aspergillus fumigatus. \u003cem\u003eJ. Microbiol.\u003c/em\u003e \u003cstrong\u003e54\u003c/strong\u003e, 223-231, doi:10.1007/s12275-016-5619-5 (2016).\u003c/li\u003e\n\u003cli\u003eSon, Y. E., Yu, J. H. \u0026amp; Park, H. S. Regulators of the Asexual Life Cycle of Aspergillus nidulans. \u003cem\u003eCells\u003c/em\u003e \u003cstrong\u003e12\u003c/strong\u003e, doi:10.3390/cells12111544 (2023).\u003c/li\u003e\n\u003cli\u003eKim, S. S., Kim, Y. H. \u0026amp; Shin, K. S. The developmental regulators, FlbB and FlbE, are involved in the virulence of Aspergillus fumigatus. \u003cem\u003eJ. Microbiol. Biotechnol.\u003c/em\u003e \u003cstrong\u003e23\u003c/strong\u003e, 766-770, doi:10.4014/jmb.1301.01002 (2013).\u003c/li\u003e\n\u003cli\u003eBertuzzi, M.\u003cem\u003e et al.\u003c/em\u003e On the lineage of Aspergillus fumigatus isolates in common laboratory use. \u003cem\u003eMed. Mycol.\u003c/em\u003e \u003cstrong\u003e59\u003c/strong\u003e, 7-13, doi:10.1093/mmy/myaa075 (2021).\u003c/li\u003e\n\u003cli\u003eKwon, N. J., Garzia, A., Espeso, E. A., Ugalde, U. \u0026amp; Yu, J. H. FlbC is a putative nuclear C2H2 transcription factor regulating development in Aspergillus nidulans. \u003cem\u003eMol. Microbiol.\u003c/em\u003e \u003cstrong\u003e77\u003c/strong\u003e, 1203-1219, doi:10.1111/j.1365-2958.2010.07282.x (2010).\u003c/li\u003e\n\u003cli\u003eCho, H. J.\u003cem\u003e et al.\u003c/em\u003e Regulation of Conidiogenesis in Aspergillus flavus. \u003cem\u003eCells\u003c/em\u003e \u003cstrong\u003e11\u003c/strong\u003e, doi:10.3390/cells11182796 (2022).\u003c/li\u003e\n\u003cli\u003eLi, Y.\u003cem\u003e et al.\u003c/em\u003e Comprehensive analysis of the regulatory network of blue-light-regulated conidiation and hydrophobin production in Trichoderma guizhouense. \u003cem\u003eEnviron. Microbiol.\u003c/em\u003e \u003cstrong\u003e23\u003c/strong\u003e, 6241-6256, doi:10.1111/1462-2920.15748 (2021).\u003c/li\u003e\n\u003cli\u003eGuo, C. T., Peng, H., Tong, S. M., Ying, S. H. \u0026amp; Feng, M. G. Distinctive role of fluG in the adaptation of Beauveria bassiana to insect-pathogenic lifecycle and environmental stresses. \u003cem\u003eEnviron. Microbiol.\u003c/em\u003e \u003cstrong\u003e23\u003c/strong\u003e, 5184-5199, doi:10.1111/1462-2920.15500 (2021).\u003c/li\u003e\n\u003cli\u003eTanaka, M.\u003cem\u003e et al.\u003c/em\u003e The C2H2-type transcription factor, FlbC, is involved in the transcriptional regulation of Aspergillus oryzae glucoamylase and protease genes specifically expressed in solid-state culture. \u003cem\u003eAppl. Microbiol. Biotechnol.\u003c/em\u003e \u003cstrong\u003e100\u003c/strong\u003e, 5859-5868, doi:10.1007/s00253-016-7419-6 (2016).\u003c/li\u003e\n\u003cli\u003eSon, H., Kim, M. G., Chae, S. K. \u0026amp; Lee, Y. W. FgFlbD regulates hyphal differentiation required for sexual and asexual reproduction in the ascomycete fungus Fusarium graminearum. \u003cem\u003eJ. Microbiol.\u003c/em\u003e \u003cstrong\u003e52\u003c/strong\u003e, 930-939, doi:10.1007/s12275-014-4384-6 (2014).\u003c/li\u003e\n\u003cli\u003eFukuda, S.\u003cem\u003e et al.\u003c/em\u003e Trade-off between Plasticity and Velocity in Mycelial Growth. \u003cem\u003emBio\u003c/em\u003e \u003cstrong\u003e12\u003c/strong\u003e, doi:10.1128/mBio.03196-20 (2021).\u003c/li\u003e\n\u003cli\u003eXiao, P., Shin, K. S., Wang, T. \u0026amp; Yu, J. H. Aspergillus fumigatus flbB encodes two basic leucine zipper domain (bZIP) proteins required for proper asexual development and gliotoxin production. \u003cem\u003eEukaryot Cell\u003c/em\u003e \u003cstrong\u003e9\u003c/strong\u003e, 1711-1723, doi:10.1128/EC.00198-10 (2010).\u003c/li\u003e\n\u003cli\u003eMoon, H.\u003cem\u003e et al.\u003c/em\u003e Unraveling the Gene Regulatory Networks of the Global Regulators VeA and LaeA in Aspergillus nidulans. \u003cem\u003eMicrobiol Spectr\u003c/em\u003e \u003cstrong\u003e11\u003c/strong\u003e, e0016623, doi:10.1128/spectrum.00166-23 (2023).\u003c/li\u003e\n\u003cli\u003eRocha, M. C.\u003cem\u003e et al.\u003c/em\u003e The Cell Wall Integrity Pathway Contributes to the Early Stages of Aspergillus fumigatus Asexual Development. \u003cem\u003eAppl. Environ. Microbiol.\u003c/em\u003e \u003cstrong\u003e86\u003c/strong\u003e, doi:10.1128/AEM.02347-19 (2020).\u003c/li\u003e\n\u003cli\u003eStewart, J. I. P.\u003cem\u003e et al.\u003c/em\u003e Reducing Aspergillus fumigatus Virulence through Targeted Dysregulation of the Conidiation Pathway. \u003cem\u003emBio\u003c/em\u003e \u003cstrong\u003e11\u003c/strong\u003e, doi:10.1128/mBio.03202-19 (2020).\u003c/li\u003e\n\u003cli\u003eMah, J. H. \u0026amp; Yu, J. H. Upstream and downstream regulation of asexual development in Aspergillus fumigatus. \u003cem\u003eEukaryot Cell\u003c/em\u003e \u003cstrong\u003e5\u003c/strong\u003e, 1585-1595, doi:10.1128/EC.00192-06 (2006).\u003c/li\u003e\n\u003cli\u003eRies, L. N. A.\u003cem\u003e et al.\u003c/em\u003e Nutritional Heterogeneity Among Aspergillus fumigatus Strains Has Consequences for Virulence in a Strain- and Host-Dependent Manner. \u003cem\u003eFront. Microbiol.\u003c/em\u003e \u003cstrong\u003e10\u003c/strong\u003e, 854, doi:10.3389/fmicb.2019.00854 (2019).\u003c/li\u003e\n\u003cli\u003eTodokoro, T.\u003cem\u003e et al.\u003c/em\u003e Identification of a novel pyrithiamine resistance marker gene thiI for genome co-editing in Aspergillus oryzae. \u003cem\u003eJ. Biosci. Bioeng.\u003c/em\u003e \u003cstrong\u003e130\u003c/strong\u003e, 227-232, doi:10.1016/j.jbiosc.2020.04.013 (2020).\u003c/li\u003e\n\u003cli\u003eOnoue, T.\u003cem\u003e et al.\u003c/em\u003e Identification of Two Mannosyltransferases Contributing to Biosynthesis of the Fungal-type Galactomannan alpha-Core-Mannan Structure in Aspergillus fumigatus. \u003cem\u003eSci. Rep.\u003c/em\u003e \u003cstrong\u003e8\u003c/strong\u003e, 16918, doi:10.1038/s41598-018-35059-2 (2018).\u003c/li\u003e\n\u003cli\u003eArai, T.\u003cem\u003e et al.\u003c/em\u003e Hmg1 mutations in Aspergillus fumigatus and their contribution to triazole susceptibility. \u003cem\u003eMed. Mycol.\u003c/em\u003e \u003cstrong\u003e59\u003c/strong\u003e, 980-984, doi:10.1093/mmy/myab026 (2021).\u003c/li\u003e\n\u003cli\u003eKadooka, C.\u003cem\u003e et al.\u003c/em\u003e Identification of an alpha-(1--\u0026gt;6)-Mannosyltransferase Contributing To Biosynthesis of the Fungal-Type Galactomannan alpha-Core-Mannan Structure in Aspergillus fumigatus. \u003cem\u003emSphere\u003c/em\u003e \u003cstrong\u003e7\u003c/strong\u003e, e0048422, doi:10.1128/msphere.00484-22 (2022).\u003c/li\u003e\n\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":true,"highlight":"","institution":"","isAcceptedByJournal":false,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true},"keywords":"Aspergillus fumigatus, flbC, fumagillin, secondary metabolism, heterogeneity","lastPublishedDoi":"10.21203/rs.3.rs-4435109/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-4435109/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003e \u003cem\u003eAspergillus fumigatus\u003c/em\u003e shows inter-strain heterogeneity in the repertoire of secondary metabolites such as mycotoxin fumagillin but the mechanism by which this heterogeneity arises in those production is still not understood. Here, we investigates the effect of the \u003cem\u003eflbC\u003c/em\u003e gene on fumagillin production in \u003cem\u003eA. fumigatus\u003c/em\u003e by introducing \u003cem\u003eflbC\u003c/em\u003e deletions into laboratory strains, A1151, A1159, and A1280 from different backgrounds and examines the effect of the \u003cem\u003eflbC\u003c/em\u003e gene on inter-strain heterogeneity. We found that, although all three laboratory strains were capable of producing fumagillin, there was heterogeneity in the effects of \u003cem\u003eflbC\u003c/em\u003e gene deletion on fumagillin production. This heterogeneity may be dependent on differences in the expression levels of the \u003cem\u003efma\u003c/em\u003e gene family involved in fumagillin production and different levels of transcriptional activation by transcription factors FapR and LaeA. The \u003cem\u003eflbC\u003c/em\u003e gene expression level peaks at the same time as peak mycelial growth but shows differential expression level, suggesting that the differences among strains in the range of expression levels are dependent on differences in the upstream expression of \u003cem\u003eflbC\u003c/em\u003e. Thus, our findings show that the different interactions between \u003cem\u003eflbC\u003c/em\u003e and factors regulating the expression of fumagillin gene cluster are the probable cause for heterogeneity in fumagillin production in the organism.\u003c/p\u003e","manuscriptTitle":"The FlbC transcription factor contributes to the generation of strain heterogeneity in Fumagillin mycotoxin production in Aspergillus fumigatus","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2024-06-04 16:52:00","doi":"10.21203/rs.3.rs-4435109/v1","editorialEvents":[{"type":"communityComments","content":0}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"0c127182-c236-4b18-9b08-37f01d9626d7","owner":[],"postedDate":"June 4th, 2024","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"posted","subjectAreas":[],"tags":[],"updatedAt":"2024-06-24T05:54:32+00:00","versionOfRecord":[],"versionCreatedAt":"2024-06-04 16:52:00","video":"","vorDoi":"","vorDoiUrl":"","workflowStages":[]},"version":"v1","identity":"rs-4435109","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-4435109","identity":"rs-4435109","version":["v1"]},"buildId":"8U1c8b4HqxoKbykW_rLl7","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}