ORIGINAL RESEARCH article

Front. Bioeng. Biotechnol., 07 August 2023

Sec. Synthetic Biology

Volume 11 - 2023 | https://doi.org/10.3389/fbioe.2023.1222812

FungalBraid 2.0: expanding the synthetic biology toolbox for the biotechnological exploitation of filamentous fungi

  • 1. Food Biotechnology Department, Instituto de Agroquímica y Tecnología de Alimentos (IATA), Consejo Superior de Investigaciones Científicas (CSIC), Valencia, Spain

  • 2. Instituto de Biología Molecular y Celular de Plantas (IBMCP), Consejo Superior de Investigaciones Científicas (CSIC)-Universitat Politècnica de València (UPV), Valencia, Spain

  • 3. Preventive Medicine and Public Health, Food Science, Toxicology and Forensic Medicine Department. Faculty of Pharmacy. Universitat de València. Vicente Andrés Estellés s/n, Valencia, Spain

Abstract

Fungal synthetic biology is a rapidly expanding field that aims to optimize the biotechnological exploitation of fungi through the generation of standard, ready-to-use genetic elements, and universal syntax and rules for contributory use by the fungal research community. Recently, an increasing number of synthetic biology toolkits have been developed and applied to filamentous fungi, which highlights the relevance of these organisms in the biotechnology field. The FungalBraid (FB) modular cloning platform enables interchangeability of DNA parts with the GoldenBraid (GB) platform, which is designed for plants, and other systems that are compatible with the standard Golden Gate cloning and syntax, and uses binary pCAMBIA-derived vectors to allow Agrobacterium tumefaciens-mediated transformation of a wide range of fungal species. In this study, we have expanded the original FB catalog by adding 27 new DNA parts that were functionally validated in vivo. Among these are the resistance selection markers for the antibiotics phleomycin and terbinafine, as well as the uridine-auxotrophic marker pyr4. We also used a normalized luciferase reporter system to validate several promoters, such as PpkiA, P7760, Pef1α, and PafpB constitutive promoters, and PglaA, PamyB, and PxlnA inducible promoters. Additionally, the recently developed dCas9-regulated GB_SynP synthetic promoter collection for orthogonal CRISPR activation (CRISPRa) in plants has been adapted in fungi through the FB system. In general, the expansion of the FB catalog is of great interest to the scientific community since it increases the number of possible modular and interchangeable DNA assemblies, exponentially increasing the possibilities of studying, developing, and exploiting filamentous fungi.

1 Introduction

Filamentous fungi have acquired a great biotechnological relevance as biofactories for the sustainable production of organic acids, proteins, enzymes, and metabolites with applications in the agri-food, chemical, pharmaceutical, textile, paper, and biofuel industries (Meyer et al., 2016). Their ability to grow on many distinct and economic substrates and plant residues, and their high secretory capacity justify the biotechnological interest of these microorganisms, which have become essential contributors to the so-called circular bio-economy (Meyer et al., 2020). Enzymes produced by fungi currently make up more than half of the enzymes used in the industry (de Vries et al., 2020). Additionally, fungal genomes contain a large number of biosynthetic gene clusters encoding potentially useful biomolecules to be exploited (Robey et al., 2021), reflecting the relevance of filamentous fungi as cell factories. However, there are still aspects that need to be improved since the conditions and levels of production of different biomolecules are highly variable, and some of them are difficult to produce in a cost-efficient manner.

Synthetic biology (SynBio) is an ever-expanding scientific field that has revolutionized genetic and metabolic engineering. SynBio provides new tools for the generation of ready-to-use, standardized, modular genetic elements to obtain microbial strains with optimized properties either by the production of specific proteins or by fine-tuning the expression of specific metabolic pathway-related genes (Benner and Sismour, 2005). In this context, fungal SynBio is rapidly evolving. Our group has adapted the GoldenBraid (GB) modular cloning platform originally developed for plants (Sarrion-Perdigones et al., 2013) to filamentous fungi, a variant called FungalBraid (FB) (https://gbcloning.upv.es/fungal/). This modular cloning method is based on type IIS restriction enzymes and pCAMBIA-derived binary vectors for Agrobacterium tumefaciens-mediated transformation (ATMT), with the main advantages of the full reusability of its DNA parts and their interchangeability between plants and fungi as long as they are functionally compatible (Hernanz-Koers et al., 2018; Vazquez-Vilar et al., 2020). The domestication or incorporation of new DNA parts into the FB system is achieved by cloning them into level 0 pUPD2 vectors, and transcriptional units (TUs) are then formed by combining different level 0 parts in a multipartite assembly into level 1 pDGB3α vectors. The GB and FB systems allow the combination of different TUs contained in two compatible pDGB3α vectors in a bipartite assembly into level 2 pDGB3Ω vectors, which can then be combined in the same way back into pDGB3α vectors, allowing the indefinite expansion of the multigene construct and designs of increased complexity.

Since the development of the FB system, an increasing number of SynBio-based applications in fungi have been reported (Dahlmann et al., 2021; Mózsik et al., 2021; Mózsik et al., 2022), which highlights the need for a boost in the SynBio toolkit for these organisms. However, there is still a shortage of tools for orthogonal and fine-tuned expression of genes applied to filamentous fungi. In this sense, an increase in the repertoire of promoters is required. Promoters with different expression levels or which are inducible and/or cell-specific would increase the flexibility and the ability to optimize expression systems, especially for proteins which can be toxic. These promoters may come from different organisms or may be created using synthetic designs. Considering that constitutive and inducible promoters are commonly used among the scientific community, synthetic promoters have been less exploited. These promoters often comprised a core or minimal promoter and an upstream region in which cis-regulatory elements are incorporated (Martins-Santana et al., 2018). These cis-regulatory elements are typically obtained from the binding sites of transcriptional regulators which activate or inactivate gene expression. Although natural transcriptional regulators limit the freedom in the design of cis-regulatory elements, the use of CRISPR activation (CRISPRa) strategies allows the use of virtually any 20-base pair (bp) sequence as a cis-regulatory box in the development of synthetic promoters (Moreno-Giménez et al., 2022). In this regard, the collection of nuclease-deactivated Cas9 (dCas9)-regulated synthetic promoter GB_SynP that has been recently developed for plants (Moreno-Giménez et al., 2022) could easily be adapted to fungi, given the interchangeability of DNA parts between GB and FB systems. Additionally, GB/FB systems provide a standard measurement using a luciferase/renilla transient assay to estimate relative expression levels of promoters, including the synthetic promoters (Vazquez-Vilar et al., 2017; Gandía et al., 2022).

In this study, we have incorporated 27 new genetic parts into the FB system, which include native strong and inducible fungal promoters, synthetic promoters, terminators, and selection markers. All these components have been validated in vivo in two economically relevant fungi: the non-model postharvest pathogen of citrus Penicillium digitatum (Palou, 2014), and in the well-known fungus with the Generally Recognized as Safe (GRAS) status and a long record of industrial use Penicillium chrysogenum (Fierro et al., 2022). The strength of the constitutive promoters has been characterized and compared in a nanoluciferase-normalized luciferase-based reporter system; the induction levels of the inducible promoters have also been quantified, and the activation of the synthetic promoters has been studied using programmable transcriptional factors based on CRISPRa (Mózsik et al., 2021). Overall, the expansion of the FB toolkit will be of great interest to the scientific community to further aid the exploitation of fungal workhorses and accelerate the discovery and production of (novel) bioactive molecules for multiple biotechnological applications.

2 Materials and methods

2.1 Strains, media, and growth conditions

The fungal strains used in this study were P. digitatum CECT 20796 (isolate PHI26) (Marcet-Houben et al., 2012) and P. chrysogenum wild-type ATCC 10002 (Q176) (Hegedüs et al., 2011). Fungi were routinely cultured on potato dextrose agar (PDA, Difco-BD Diagnostics) plates for 7 days at 25°C. For transformation, the generated vectors were amplified in the bacterium Escherichia coli JM109 and grown in Luria–Bertani (LB) medium at 37°C with either 50 μg/mL chloramphenicol, 50 μg/mL kanamycin, 100 μg/mL spectinomycin, or 100 μg/mL ampicillin depending on the vector. A. tumefaciens AGL-1 was cultured in LB medium at 28°C with 20 μg/mL rifampicin and the corresponding antibiotic depending on the vector used.

For growth profiles, 5 µL of conidial suspension (5 × 104 conidia/mL) was deposited on the center of PDA plates, and colony morphology was assessed and compared daily by visual inspection.

2.2 Design, domestication, and DNA assemblies of genetic elements

All the genetic elements that have been incorporated on the FB platform are listed in Table 1. New DNA parts were domesticated according to GB rules and tools (https://gbcloning.upv.es) and ordered from an external company as synthetic genes (gBlocks™, IDT). In the case of the coding sequence (CDS) of pyr4, the gene from Trichoderma reesei was codon-optimized according to the optimal codon frequency of Penicillium genera prior to GB/FB domestication. Domesticated elements were ligated into the pUPD2 entry vector via the restriction–ligation protocol, as previously described (Hernanz-Koers et al., 2018; Vazquez-Vilar et al., 2020). Positive E. coli clones were confirmed by routine PCR amplifications and Sanger sequencing using external specific primers OJM524 and OJM525, which are designed for pUPD2 vectors (Hernanz-Koers et al., 2018) (Table 2). Multiple assemblies into pDGB3α vectors of the DNA parts contained in pUPD2 vectors were carried out to obtain different TUs, and binary assemblies were subsequently performed to combine different TUs into multigenetic constructs within these pDGB3α or pDGB3Ω vectors, as previously described (Hernanz-Koers et al., 2018; Vazquez-Vilar et al., 2020).

TABLE 1

Selection marker
Auxotrophy
CodeNamePlasmidDescriptionReference
FB271*Ppyr4pUPD2Promoter of the pyr4 gene from T. reeseiThis study
FB272*pyr4pUPD2Coding sequence of the pyr4 gene from T. reeseiThis study
FB273*Tpyr4pUPD2Terminator of the pyr4 gene from T. reeseiThis study
FB293*TU_pyr4pDGB3α2Assembly of the transcriptional unit for the auxotrophy marker pyr4 from T. reseeiThis study
FB359*5′ upstream Pdig pyrGpUPD25′ upstream region of the pyrG gene in P. digitatumThis study
FB361*3′ downstream Pdig pyrGpUPD23′ downstream region of the pyrG gene in P. digitatumThis study
FB372*FB359+FB361pDGB3α1Assembly for pyrG deletion in P. digitatumThis study
Resistance
CodeNamePlasmidDescriptionReference
FB413*blepUPD2Coding sequence for phleomycin resistanceThis study
FB414*ergApUPD2Coding sequence for terbinafine resistanceThis study
FB411*PpcbCpUPD2Promoter of isopenicillin N synthase from P. rubensThis study
FB416*TamdSpUPD2Terminator of acetamidase-encoding gene amdS from A. nidulansThis study
FB430*Ppcbc:ble:TamdSpDGB3α2TU for the expression of phleomycin resistanceThis study
FB431*PgpdA:ergA:TamdSpDGB3α2TU for the expression of terbinafine resistanceThis study
Constitutive/inducible promoters
CodeNamePlasmidDescriptionReference
FB007PgpdApUPD2Promoter of glyceraldehyde-3-phosphate dehydrogenase from A. nidulansHernanz-Koers et al. (2018)
FB291*PxlnApUPD2Promoter of the endo-1,4-beta-xylanase A gene from A. nidulans, xylose-inducibleThis study
FB389*PpkiApUPD2Promoter of the highly expressed pyruvate kinase gene from A. nigerThis study
FB404*PafpBpUPD2Promoter of the antifungal protein afpB gene from P. digitatumThis study
FB405*PglaApUPD2Promoter of the glucoamylase gene from A. niger, maltose/starch-inducibleThis study
FB406*PamyBpUPD2Promoter of the TAKA-amylase A gene from A. oryzae, maltose/starch-inducibleThis study
FB407*Pef1αpUPD2Promoter of the elongation factor 1-α gene from P. digitatum (PDIG_59570)This study
FB408*P07760pUPD2Promoter of the ubiquitin ligase gene from P. digitatum (PDIG_07760)This study
GB0096Luciferase (Luc)pUPDCoding sequence for the firefly luciferase proteinSarrion-Perdigones et al. (2013)
FB001PtrpCpUPD2Promoter of the multifunctional tryptophan biosynthesis protein coding gene trpC from A. nidulansHernanz-Koers et al. (2018)
FB002TtubpUPD2Terminator of the tubulin-encoding gene from N. crassaHernanz-Koers et al. (2018)
FB008TtrpCpUPD2Terminator of the trpC gene from A. nidulansHernanz-Koers et al. (2018)
FB009PtrpC:nptII:TtubpDGB3α2TU for the expression of geneticin resistanceHernanz-Koers et al. (2018)
FB312PgpdA:Nluc:TtubpDGB3α1RTU for the expression of nanoluciferase (Nluc) under PgpdAGandía et al. (2022)
FB323PgpdA:Nluc:Ttub::PtrpC:nptII:Ttub:: Ppaf:Luc:TtubpDGB3α1Module for the expression of geneticin resistance and nanoluciferase, and the expression of firefly luciferase under the PpafGandía et al. (2022)
FB367*PgpdA:Nluc:Ttub::PtrpC:nptII:TtubpDGB3α1Module for the expression of geneticin resistance and nanoluciferaseThis study
FB417PafpB:Luc:TtrpCpDGB3α2TU for the expression of luciferase under PafpBThis study
FB418PglaA:Luc:TtrpCpDGB3α2TU for the expression of luciferase under PglaAThis study
FB419PamyB:Luc:TtrpCpDGB3α2TU for the expression of luciferase under PamyBThis study
FB420Pef1α:Luc:TtrpCpDGB3α2TU for the expression of luciferase under Pef1αThis study
FB421P07760:Luc:TtrpCpDGB3α2TU for the expression of luciferase under P07760This study
FB423PgpdA:Luc:TtrpCpDGB3α2TU for the expression of luciferase under PgpdAThis study
FB424PxlnA:Luc:TtrpCpDGB3α2TU for the expression of luciferase under PxlnAThis study
FB426PpkiA:Luc:TtrpCpDGB3α2TU for the expression of luciferase under PpkiAThis study
FB432FB367+FB417pDGB3Ω1Module for the expression of geneticin resistance, Nluc under the PgpdA promoter, and luciferase under PafpBThis study
FB433FB367+FB418pDGB3Ω1Module for the expression of geneticin resistance, Nluc under the PgpdA promoter, and luciferase under PglaAThis study
FB434FB367+FB419pDGB3Ω1Module for the expression of geneticin resistance, Nluc under the PgpdA promoter, and luciferase under PamyBThis study
FB435FB367+FB420pDGB3Ω1Module for the expression of geneticin resistance, Nluc under PgpdA promoter, and luciferase under Pef1αThis study
FB436FB367+FB421pDGB3Ω1Module for the expression of geneticin resistance, Nluc under the PgpdA promoter, and luciferase under P07760This study
FB438FB367+FB423pDGB3Ω1Module for the expression of geneticin resistance, Nluc under the PgpdA promoter, and luciferase under PgpdAThis study
FB439FB367+FB424pDGB3Ω1Module for the expression of geneticin resistance, Nluc under the PgpdA promoter, and luciferase under PxlnAThis study
FB441FB367+FB426pDGB3Ω1Module for the expression of geneticin resistance, Nluc under the PgpdA promoter, and luciferase under PpkiAThis study
dCas9-activated synthetic promoters
CodeNamePlasmidDescriptionReference
GB2815RandomSequence R1pUPD2Random sequence R1 of 1240 bp for the A1 distal promoter positionMoreno-Giménez et al. (2022)
GB2878G1aG2b.1pUPD2A2 proximal promoter sequence containing the target sequence for gRNA1 flanked by random sequencesMoreno-Giménez et al. (2022)
GB2885G1ab.1pUPD2A2 proximal promoter sequence consisting of two times the target sequence for gRNA1 flanked by random sequencesMoreno-Giménez et al. (2022)
GB3276G1abc.3pUPD2A2 proximal promoter sequence containing three times the target sequence for gRNA1 flanked by random sequencesMoreno-Giménez et al. (2022)
GB3413mPAFpUPD2Minimal promoter of the paf gene from P. chrysogenum, containing 62 bp upstream the transcription start site and the 5′UTR regionMoreno-Giménez et al. (2022)
FB395*R1:G1aG2b.1:mPAF:Luc:TtrpcpDGB3α2TU for the expression of luciferase under a synthetic promoter containing the target sequence for gRNA1 (1xLuc)This study
FB396*R1:G1ab.1:mPAF:Luc:TtrpcpDGB3α2TU for the expression of luciferase under a synthetic promoter containing two times the target sequence for gRNA1 (2xLuc)This study
FB397*R1:G1abc.3:mPAF:Luc:TtrpcpDGB3α2TU for the expression of luciferase under a synthetic promoter containing three times the target sequence for gRNA1 (3xLuc)This study
FB398*FB367+FB395pDGB3Ω1Module for the expression of geneticin resistance, Nluc under the PgpdA promoter, and luciferase under a synthetic promoter containing the target sequence for gRNA1This study
FB399*FB367+FB396pDGB3Ω1Module for the expression of geneticin resistance, Nluc under the PgpdA promoter, and luciferase under a synthetic promoter containing two times the target sequence for gRNA1 (2xLuc)This study
FB400*FB367+FB397pDGB3Ω1Module for the expression of geneticin resistance, Nluc under the PgpdA promoter, and luciferase under a synthetic promoter containing three times the target sequence for gRNA1 (3xLuc)This study
FB403pAMA18.0_gRNA1pAMA18.0Expression plasmid for the dCas9 activation system and the (GB_SynP) gRNA1This study

FB parts reported in this study. DNA parts are grouped according to the purpose for which they were used.

*DNA parts from this study which are deposited in Addgene.

TABLE 2

IDUse*Sequence 5′-3′**Tm (°C)OriginPurposeReference
OJM371FATA​GAT​CTA​ACT​GAT​ATT​GAA​GGA​GCA52PtrpCMolecular characterizationThis study
OJM509FGCGCCGTCTCGCTCGGGAGTGGCGCATGCGGACAGACGG64PgpdAMolecular characterizationHernanz-Koers et al. (2018)
OJM522RGCG​CCG​TCT​CGC​TCA​AGC​GCA​TGT​CTC​AGA​CGGTCGATG62TtrpCMolecular characterizationHernanz-Koers et al. (2018)
OJM524FGCT​TTC​GCT​AAG​GAT​GAT​TTC​TGG70pUPD2Molecular characterizationHernanz-Koers et al. (2018)
OJM525RCAG​GGT​GGT​GAC​ACC​TTG​CC66pUPD2Molecular characterizationHernanz-Koers et al. (2018)
OJM555RTCA​TCA​TGC​AAC​ATG​CAT​GTA58TtubMolecular characterizationHernanz-Koers et al. (2018)
OJM655RCAT​CCA​TAC​TCC​ATC​CTT​CCC60pAMA18.0Molecular characterization and sequencingThis study
OJM656FCAT​TTT​TGT​CGT​CAT​GTG​CTG​G555′ Pdig pyrGMolecular characterizationThis study
OJM657RGAA​GGC​TGA​ACT​CAC​TGT​GG553′ Pdig pyrGMolecular characterizationThis study
OJM662FGCT​TTT​GCT​AAC​CAT​TTG​GGA​CAC52GB B6 codeCloning of FB372 into pDGB3α1This study
OJM663RAGC​GGT​GTC​CCA​AAT​GGT​TAG​CAA52GB C1 codeCloning of FB372 into pDGB3α1This study
OJM696FTTG​TCT​CAC​TCT​CTC​TTT​TCC51pyr4Molecular characterizationThis study
OJM697RATT​CCA​TGC​TTC​CAG​ATC​C51pyr4Molecular characterizationThis study
OJM698FATGGTCTCACCGACCAGTCCTGATGAGTCCGTGAGGACGAAACGAG60pAMA18.0Cloning of gRNA1 into pAMA18.0This study
OJM699RATGGTCTCTAAACTCTTCTCTCACCAACCAGTC GACGAGCTTACTCGTTTCGTCCTCACGGACTCA60pAMA18.0Cloning of gRNA1 into pAMA18.0This study
OJM705FTCC​TGG​AAG​TGC​GTT​GAT​CA51PxlnAMolecular characterizationThis study
OJM706FGGA​AGA​GAA​AAC​CTC​CGA​GTA​C54PpkiAMolecular characterizationThis study
OJM707FATG​AAT​TCC​ACC​GAA​TGC​AC53PafpBMolecular characterizationThis study
OJM708FTGC​CAT​TGG​CGG​AGG​GGT​CC53PglaAMolecular characterizationThis study
OJM709FTCA​ACT​GAT​TAA​AGG​TGC​CG53PamyBMolecular characterizationThis study
OJM710FGTG​AAA​AAA​CGG​ATG​GGG​AC53Pef1αMolecular characterizationThis study
OJM711FGAT​AAT​GGT​GAT​TCG​GCG​CG53P07760Molecular characterizationThis study
OJM715FGTA​TCT​GCA​TGT​TGC​ATC​GG53PpcbCMolecular characterizationThis study
OJM716RTAC​CGC​TCG​TAC​CAT​GGG​TT53TamdsMolecular characterizationThis study

Primers used in this study.

*F, forward; R, reverse.

**gRNA1 sequence is underlined.

The sequence of the single-guide RNA 1 (gRNA1) required to activate GB_SynP promoters was checked for the absence of off-target mutations in the P. digitatum genome using Geneious Prime software (https://www.geneious.com/), and was cloned into the pAMA18.0 vector, as described previously (Mózsik et al., 2021). Briefly, a primer pair was designed (OJM698 and OJM699, Table 2) which contained the target sequence of gRNA1, the hammerhead ribozyme, the inverted repetition of the 5′-end of the spacer sequence, and the recognition sites of BsaI. The resulting PCR product was purified (Wizard SV Gel and PCR Clean-Up System, Promega) and inserted into pAMA18.0 via the restriction–ligation reaction with BsaI and T4 ligase. Correct assemblies of the resulting pAMA18.0_gRNA1 vector were confirmed by Sanger sequencing.

2.3 Fungal transformation and mutant confirmation

Transformation of P. digitatum CECT 20796 (PHI26) and P. chrysogenum ATTC 10002 (Q176) with the corresponding FB binary vectors described in Table 1 was performed through ATMT, as previously described (Khang et al., 2007), with some modifications (Harries et al., 2015; Vazquez-Vilar et al., 2020).

In the case of P. digitatum uridine-auxotrophic ΔpyrG mutants, in which the pyrG gene (gene ID PDIG_38390) was deleted by homologous recombination without the insertion of any positive selection marker, mutants were selected on PDA supplemented with 1.22 g/L uridine (Sigma-Aldrich) and 1.25 g/L of 5-fluoroorotic acid (5-FOA, Formedium). For the ectopically complemented P. digitatum ΔpyrG:pyr4 strains, mutants were selected on PDA plates. P. digitatum and P. chrysogenum ectopic transformants containing ble TU (phleomycinR) were selected on PDA plates supplemented with 35 μg/mL and 25 μg/mL phleomycin (InvivoGen), respectively. Finally, P. chrysogenum ectopic transformants carrying ergA TU (terbinafineR) were selected on PDA plates supplemented with 0.5 μg/mL terbinafine hydrochloride (Sigma-Aldrich).

P. digitatum and P. chrysogenum transformants carrying different luciferase reporter constructs to test the constitutive, inducible, and GB_SynP promoters were selected on PDA plates containing 25 μg/mL geneticin (G418) (InvivoGen).

All transformants were molecularly confirmed by PCR reactions using NZYTaq II DNA polymerase (NZYTech) (Figures 13; Supplementary Figures S1–S3) from genomic DNA isolated with the NZY Tissue gDNA isolation kit (NZYTech) and primers purchased from IDT (Table 2).

FIGURE 1

FIGURE 3

For the validation of GB_SynP promoters, P. digitatum protoplasts from strains carrying luciferase reporter vectors for each of the three synthetic promoters tested (FB398, FB399, and FB400) (Table 1) were transformed with the self-replicative AMA1-based dCas9-containing plasmid pAMA18.0_gRNA1 (Mózsik et al., 2021), as previously described (Garrigues et al., 2022). Transformants were selected on PDA plates containing 0.95 M sucrose and 35 μg/mL phleomycin. Regarding the reusability of the system, the loss of the pAMA18.0_gRNA1 plasmid was confirmed after three consecutive streaks of the transformants in non-selective PDA plates, as previously described (Garrigues et al., 2022).

2.4 Luciferase/nanoluciferase assays

P. digitatum and P. chrysogenum strains carrying the luciferase reporter vectors for each of the tested promoters (FB433, FB434, FB435, FB436, FB438, FB439, and FB441) and P. digitatum strains with the luciferase reporter for the GB_SynP synthetic promoters (FB398, FB399, and FB400), carrying or not pAMA18.0_gRNA1, were grown in duplicate for 2 days in 100 mL flasks with 25 mL of liquid potato dextrose broth (PDB, Difco-BD Diagnostics) at 25°C, with shaking (150 rpm). For induction assays, transformants were grown in duplicate for 4 days in 100 mL flasks with 25 mL of either P. digitatum or P. chrysogenum minimal medium (PdMM or PcMM, respectively) (Sonderegger et al., 2016) using 2% D-glucose (PanReac), 2% maltose (Sigma-Aldrich) or 2% D-xylose (Sigma-Aldrich) as the sole carbon source. Grown mycelia were filtered, and a sample of 20 mg was collected and immediately frozen in liquid nitrogen. Luciferase and nanoluciferase measurements were performed using the Dual-Glo® Luciferase Assay System (Promega), as previously described (Gandía et al., 2022). Briefly, frozen samples were homogenized in 180 µL passive lysis buffer with a pestle and centrifuged (12,000 ×g, 10 min at 4°C). Then, 10 µL of the supernatant was transferred to a white 96-well plate (Thermo Fisher Scientific) and mixed with 40 µL of the luciferase reagent to measure luciferase luminescence in a CLARIOStar microplate reader (BMG LABTECH GmbH) with a measurement of 10 s and a delay of 2 s. Nanoluciferase luminescence was quantified thereafter by adding 40 µL of the Stop & Glo reagent and measured in the same way.

The luciferase/nanoluciferase ratio was determined for each sample, and normalized luminescence was calculated as the mean value of the ratios obtained from each duplicate. Statistical analyses were performed with GraphPad Prism 8.0.1 software. Differences between the strains were analyzed using one-way ANOVA, followed by post hoc multiple comparison Tukey’s test (p < 0.05). For induction experiments, we analyzed the differences in the growth of each strain in the presence of different carbon sources using Student’s t-test (p < 0.05).

2.5 Fruit infection assays

P. digitatum parental and mutant strains were inoculated on freshly harvested oranges (Citrus sinensis L. Osbeck cv Lane late), as previously described (González-Candelas et al., 2010). Briefly, three replicates of five orange fruits were inoculated with 5 μL of fungal conidial suspension (104 conidia/mL) at four equidistant wounds around the equator. Control mock inoculations were performed with 5 μL of sterile Milli-Q H2O. Once inoculated, fruits were maintained at 20°C and 90% relative humidity for up to 6 days. Each inoculated wound was scored daily for infection symptoms on consecutive days post-inoculation (dpi). We repeated the experiments twice. Differences in the percentage of infection for each strain compared to the control CECT 20796 were analyzed using Student’s t-test (p < 0.05) for each individual dpi.

3 Results

3.1 Selection markers for antibiotic resistance

The FB platform already contains some commonly used positive fungal selection markers based on antibiotic resistance, such as hph (hygromycinR, FB003) or nptII (geneticinR, FB009) (Hernanz-Koers et al., 2018). However, in the case of integrative approaches, multiple genetic modifications often depend on the availability of different antibiotic resistance genes for transformant selection, which can be a bottleneck for the exploitation of filamentous fungi. In this study, we expand the range of selection markers available in the FB platform by including two alternative antibiotic resistance-inducing genes, the ble resistance gene from the bacterial transposon Tn5 and the squalene epoxidase ergA gene from P. chrysogenum. The expression of the ble gene provides the selection of the antibiotic phleomycin (Austin et al., 1990), whereas the expression of the ergA gene provides resistance against the antibiotic terbinafine in a broad range of filamentous fungi (Austin et al., 1990; Sigl et al., 2010).

In order to include ble resistance in the FB platform, a functional TU was generated. For this, we assembled the ble coding sequence (FB413), together with the promoter of isopenicillin N synthase (PpcbC) from Penicillium rubens (FB411) (Polli et al., 2016) and the terminator from the acetamidase (Tamds) from Aspergillus nidulans (FB416) (Kelly and Hynes, 1985) into the pDGB3α2 vector to obtain FB430 (Table 2; Figure 1A) via restriction–ligation reactions. To functionally validate the resulting construct, we transformed P. chrysogenum and P. digitatum wild-type strains with the same FB430 via ATMT for the ectopic integration of the ble TU. P. chrysogenum transformants grown in the presence of 25 μg/mL phleomycin were selected and analyzed by PCR for the presence of the ble cassette (Figure 1B). The positive transformants, PCEM43053 and PCEM43062, showed growth on phleomycin-containing plates when compared to the parental ATCC 10002 (Figure 1C), further demonstrating the functionality of FB430. In parallel, FB430 was also validated in P. digitatum (Figures 1D–G). P. digitatum transformants grown in the presence of 35 μg/mL phleomycin were selected and confirmed by PCR (Figure 1D). The positive transformants PDZS43023, PDZS43041, and PDZS43051 were able to grow on phleomycin-containing plates (Figure 1E) and showed the same pathogenicity as the parental CECT 20796 in orange fruits (Figures 1F,G).

Similarly, to incorporate the terbinafine resistance-inducing gene in the FB platform, a functional TU for ergA was generated and validated in P. chrysogenum (Figure 2). We assembled the ergA coding sequence (FB414), together with the PgpdA promoter (FB007) and the Tamds terminator (FB416), into the pDGB3α2 vector to obtain the FB431 construct (Table 2; Figure 2A). P. chrysogenum transformants grown on 0.5 μg/mL terbinafine were chosen and confirmed by PCR (Figure 2B). The positive transformants PDZS43122, PDZS43131, and PDZS43142 could grow on phleomycin-containing PDA plates in contrast to the parental ATCC 10002 (Figure 2C), demonstrating the functionality of FB431.

FIGURE 2

Overall, both resistances were transformed ectopically to avoid any bias regarding the targeting of specific loci, and these experiments validated the use of FB430 and FB431 as standardized TUs for conferring positive selection in the transformation of different fungal species, expanding the antibiotic resistance selection markers currently available in the FB system.

3.2 Selection markers based on fungal auxotrophy

To date, no auxotrophic markers have been included in the FB platform despite the fact that they are sustainable alternatives for the use of antibiotics in transformant selection. The orotidine 5′-phosphate decarboxylase pyr4 gene from T. reesei is widely used as an auxotrophic selection marker that can be counter-selected using 5-FOA or fully supplemented using uridine (Díez et al., 1987; Derntl and Kiesenhofer, 2015). In this study, we set up several experiments to design, test, and validate pyr4 as a selection marker in pyr4/pyrG-deficient fungal strains. As a first step, uridine-auxotrophic P. digitatum ΔpyrG mutants were generated through ATMT using FB372 as the template for homologous recombination at the pyrG locus (Table 2; Supplementary Figure S4). Transformants were selected on PDA plates supplemented with 1.22 g/L uridine and 1.25 g/L 5-FOA, and were molecularly and phenotypically characterized (Supplementary Figure S4). Growth profiles showed that after pyrG deletion, P. digitatum mutants could no longer grow on PDA plates unless supplemented with uridine, confirming their auxotrophic condition. Additionally, these mutants could also grow in the presence of uridine and 5-FOA in contrast to the parental CECT 20796, further confirming pyrG deletion. Finally, infection assays on orange fruits revealed that P. digitatum ΔpyrG mutants showed highly reduced pathogenesis compared to the control (Supplementary Figure S4). Once the pyrG deletion mutants were obtained, a functional TU for the T. reesei pyr4 gene was generated. For this, we assembled the Penicillium codon-optimized and extensively domesticated pyr4 coding sequence (FB272), promoter (FB271), and terminator (FB276) into the pDGB3α2 vector to obtain FB293 (Table 2; Figure 3A) via restriction–ligation reactions. To functionally validate the resulting construct FB293, we transformed the P. digitatum ΔpyrG mutant PDSG37213 with FB293 via ATMT for the ectopic integration of pyr4 TU. Transformants grown on PDA plates were assessed by PCR (Figure 3B) and phenotypically analyzed to confirm pyrG:pyr4 complementation and, therefore, the absence of the auxotrophy. As shown in Figure 3C, complemented mutants PDSG29312, PDSG29321, and PDSG29333 were all able to grow on PDA plates without uridine, in contrast to the auxotrophic parental PDSG37213. Remarkably, P. digitatum complemented mutants that fully recovered their original pathogenicity (Figures 3D,E), which validates pyr4 as an auxotrophic selectable marker also for (phyto)pathogenic fungi, in which the deletion of pyrG orthologs has been demonstrated to decrease pathogenicity and virulence in the corresponding fungi (Zameitat et al., 2007; Higashimura et al., 2022).

3.3 Constitutive and inducible promoters

In order to further expand and characterize the promoter catalog available in the FB platform, a series of promoters from different fungal species were included in the collection (Table 1) and were functionally validated using a luciferase reporter system, as previously described (Gandía et al., 2022). This reporter consists of two TUs: the nanoluciferase coding sequence (Nluc, FB310) under the regulation of the PgpdA promoter that serves as an internal standard for normalization, and the luciferase sequence (Luc, GB0096) under the regulation of the promoter to be tested. Promoter strength is expressed as the ratio of the Luc signal divided by the Nluc internal standard. The promoters to be evaluated included the previously characterized strong pyruvate kinase gene promoter (PpkiA) from Aspergillus niger (FB389) (de Graaff et al., 1992). Novel promoters from P. digitatum included the antifungal protein AfpB gene promoter PafpB (FB404) (Garrigues et al., 2017), and two promoters with high expression levels were reported in a previous transcriptomic study: the elongation factor 1α gene promoter (Pef1α) (FB407) and the ubiquitin ligase PDIG_07760 gene promoter (P07760) (FB408) (Ropero-Pérez et al., 2023). The inducible promoters included in this study are the endo-1,4-β-xylanase A gene promoter (PxlnA) from A. nidulans (FB291, xylose-responsive) (Orejas et al., 1999), the glucoamylase gene promoter (PglaA) from A. niger, and the TAKA-amylase A gene promoter (PamyB) from Aspergillus oryzae (FB405 and FB406, respectively, both maltose/starch-responsive) (Fowler et al., 1990; Tsuchiya et al., 1992b). The widely used glyceraldehyde-3-phosphate dehydrogenase promoter PgpdA from A. nidulans, which was already available in the FB collection (FB007) (Hernanz-Koers et al., 2018), was also included in the analysis, as well as the luciferase reporter construct for the P. chrysogenum antifungal protein PAF promoter (Ppaf) from previously published data (FB323) (Gandía et al., 2022) to serve as references. To facilitate the cloning of new luciferase reporter constructs, the nanoluciferase reference gene and the nptII resistance gene were cloned into a pDGB3α1 vector (FB367) to be combined in a single reaction with the luciferase TUs cloned into pDGB3α2 vectors that included the promoters to be tested (Figure 4A). Due to the requirements of the FB binary assembly, an insulator sequence (GB3458) was also included at the 3’ end of the FB367 vector, which also helps prevent the interaction between Nluc and Luc TUs. The luciferase reporter constructs for each of the assayed promoters (FB432 to FB441) showed different normalized luciferase expression levels in P. digitatum after 2 days of growth in PDB (Figure 4B). The lowest expression levels were observed for the inducible promoters PxlnA, PamyB, and PglaA, in this order, due to the lack of inducers in this medium. Their expression levels, together with that driven by PafpB, were slightly above the basal signal observed in the control strain, but these were not statistically significant. The expression driven by PpkiA, PgpdA, and the new Pef1α were similar to those observed for Ppaf, while P07760 showed intermediate expression values between these expression and the inducible promoters.

FIGURE 4

The selected constructs carrying PxlnA, PpkiA, PglaA, PamyB, Pef1α, and P07760 were also transformed into P. chrysogenum, in which the luciferase expression level was about 10 times lower than the overall levels observed in P. digitatum (Figure 4C), except for Pef1α and PamyB, which showed a similar signal in both fungal chassis (0.27 and 0.1 for Pef1α in P. digitatum and P. chrysogenum, respectively, and 0.0007 and 0.0004 for PamyB). Unlike P. digitatum, the signal driven by PamyB was above the basal signal in P. chrysogenum and showed similar expression levels to that of the P7760 promoter. Relative expression levels among the other promoters were nevertheless maintained in both fungi, with PpkiA and Pef1α signals similar to that of the Ppaf reference promoter and with PxlnA-, PglaA-, and PamyB-driven signals similar to that of the control strains.

3.4 Induction of PglaA, PamyB, and PxlnA promoters

In order to further characterize the inducible promoters included in this study, we analyzed the induction of luciferase expression directed by PglaA, PamyB, and PxlnA promoters in P. digitatum (Figure 5A) and P. chrysogenum strains (Figure 5B) after 4 days of growth in minimal medium (PdMM for P. digitatum and PcMM for P. chrysogenum) using different inducers as the sole carbon source (2% maltose for PglaA and PamyB, and 2% xylose for PxlnA). When the fungi were grown in the presence of the inducer, the expression levels driven by all three promoters were significantly higher than those observed in the reference media with 2% glucose, which increase by 8x, 4x, and 10x for PglaA, PamyB, and PxlnA, respectively, in P. digitatum (Figure 5A); and by 2x, 9x, and 8x in P. chrysogenum (Figure 5B). Signals observed for the PxlnA promoter in the reference media (MM + glucose) were similar to the basal signal of the reference strains in both fungi, while PamyB and PglaA promoters showed higher basal expression in the same media, especially in P. chrysogenum. Remarkably, expression levels in the reference Ppaf promoter were found to increase significantly in P. chrysogenum when grown in the PcMM supplemented with maltose or xylose (0.22 a.u. on average) when compared to the medium supplemented with glucose (0.015 a.u. on average, 15 times lower). The same was observed, but to a lesser extent, in P. digitatum when Ppaf was expressed in PdMM glucose (average values of 0.6 a.u.) compared to its expression in PdMM xylose (0.15 a.u. average values, four times higher).

FIGURE 5

3.5 dCas9-activated synthetic promoters

Finally, we tested the recently developed GB_SynP (Moreno-Giménez et al., 2022) in our fungal chassis in combination with the pAMA18.0_gRNA1 plasmid, which delivers the CRISPRa system necessary to activate GB_SynP promoters in a non-integrative manner (Mózsik et al., 2021). To this end, we developed luciferase reporter constructs following the same procedure as for natural promoters (Figure 6A). In these constructs, luciferase expression was regulated by synthetic promoters consisting of an A1 distal promoter part formed by a random sequence (GB2815), and an A2 proximal promoter part including the target sequence for gRNA1 that was repeated once (GB2878), twice (GB2885), or three times (GB3276), and the minimal promoter mPAF (GB3423) derived from the native Ppaf from the fungus P. chrysogenum, which was previously found to drive a strong induction when exposed to the dCas9 system loaded with gRNA1 in plants (Moreno-Giménez et al., 2022). The resulting constructs containing one gRNA1 target (FB395, 1xLuc), two targets (FB396, 2xLuc), or three targets (FB397, 3xLuc) were stably transformed into P. digitatum via ATMT. Protoplasts obtained from these strains were re-transformed with the CRISPRa expression vector pAMA18.0_gRNA1 (Figure 6B). The expression levels of 1xLuc and 2xLuc constructs were not significantly higher than the basal signal of the control strain despite the presence of pAMA18.0_gRNA1, except for one of the 1xLuc re-transformants, which showed a low, but statistically significant increase in the luciferase expression when compared to the same strain in the absence of the CRISPRa system. A higher and significant increase of approximately 16 times on average in the luciferase signal was observed for the 3xLuc construct (Figure 6C) in all tested re-transformants. Although the expression driven by this promoter was approximately 10 times lower than that observed with the reference Ppaf promoter, it is comparable to that achieved with PglaA and PamyB, showing the functionality of these dCas9-activated synthetic promoters in fungi.

FIGURE 6

4 Discussion

The FB cloning platform allows for the open exchange of standardized, ready-to-use DNA parts in the fungal research community (Hernanz-Koers et al., 2018). Moreover, if the platform is functionally compatible and validated, it also allows for the exchange of parts between plants and fungi, as occurred with the fluorescent YFP protein or the hygromycin selection marker reported previously (Hernanz-Koers et al., 2018) or the synthetic promoters reported in this study. However, the number of validated genetic elements present in the FB platform was very limited to date, which hindered the biotechnological exploitation of filamentous fungi. In this study, we have expanded the available genetic elements in FB platforms by incorporating one auxotrophic selection marker (pyr4), two additional antibiotic resistance markers (ergA and ble), two strong promoters (PpkiA and Pef1α), two intermediate promoters (PafpB and P7760), three inducible promoters (PglaA, PamyB, and PxlnA), and three versions of the dCas9-regulated GB_SynP synthetic promoters. Even though the validation of these new parts has been performed in Penicillium species, the FB system has been demonstrated to mediate the expression of the same construct in different fungal genera. For instance, the FB027 construct used for the expression of YFP has been functionally validated in P. digitatum, Penicillium expansum, and A. niger (Hernanz-Koers et al., 2018; Vazquez-Vilar et al., 2020). Therefore, the FB system and, subsequently, the new FB parts described here, are expected to be of use in a wide range of fungal species of different genera.

Since the FB release, there have been an increasing number of SynBio-based genetic toolkits developed for filamentous fungi (Dahlmann et al., 2021; Mózsik et al., 2021; Mózsik et al., 2022). In this sense, FB, which derives from the GB cloning framework, shares most of the codes and type IIS restriction enzymes with these alternative SynBio collections (Weber et al., 2011), making it possible to combine code-compatible level 0 plasmids between these systems to assemble TUs into level 1 plasmids. However, these Golden Gate-based collections alternative to FB use plasmids derived from pAMA1 or pEHN8, which are introduced into fungal cells via protoplast transformation. In contrast, the FB collection is based on pCAMBIA-derived vectors and can be applied to a broad spectrum of fungal species that are compatible with ATMT (de Groot et al., 1998), which is considered to be a more advantageous transformation method than protoplasts as spores can be used directly for genetic transformation and transformation efficiencies are generally higher (Li et al., 2017 and references, therein). Moreover, unlike these other Golden Gate cloning systems, FB/GB systems allow for the indefinite expansion of multigenetic constructs via bipartite assemblies between pDG3α and pDGB3Ω vectors (Sarrion-Perdigones et al., 2011). To date, the GB system has permitted the assembly of up to 10 TUs (GB3243) (Selma et al., 2022) and inserted an assembly as large as 20 kb (GB4559–GB4585) (Moreno-Giménez et al., 2022), yet the transformation and propagation of larger constructs into E. coli might be hampered by the limitations to this host (Weber et al., 2011). In this regard, the adaptation of other ATMT-compatible vectors into FB/GB could be considered, such as the binary-BAC (BIBAC) vector reported by Hamilton (1997), which can carry >100 kb and has already been used to transform Fusarium, Aspergillus, or Ustilago species (Takken et al., 2004; Ali and Bakkeren, 2011).

Among the new genetic elements in the FB system, we included three commonly used fungal selection markers, two of them based on antibiotic resistance (ble and ergA, which confer resistance to phleomycin and terbinafine, respectively) and one based on fungal auxotrophy (pyr4), further expanding the possibilities for fungal transformation and mutant selection within the frame of the FB platform. Although antibiotic resistance markers are among the most widely used approaches for positive transformant selection, auxotrophic markers are more sustainable alternatives for the use of antibiotics, which can have undesired side effects on the fitness of the organism under study or cause unwanted spontaneous resistance. The orotidine 5′-phosphate decarboxylase-encoding gene pyr4 from T. reesei, which is an ortholog of the Aspergillus and Penicillium pyrA/pyrG gene, is widely applied as a strong auxotrophic selection marker that can be counter-selected using 5-FOA or fully supplemented using uracil or uridine (Díez et al., 1987; Derntl and Kiesenhofer, 2015). Interchangeability of these two orthologs has already been demonstrated between fungi from different phylogenetic classes, from Sordariomycetes to Ascomycetes and vice versa (Ballance and Turner, 1985; Díez et al., 1987; Gruber et al., 1990). Therefore, the FB-adapted pyr4 TU is expected to restore uridine/uracil auxotrophy in a broad range of fungal species. In the case of fungal (plant) pathogens, for which pyr disruption has been reported to reduce pathogenicity and virulence (Zameitat et al., 2007; Higashimura et al., 2022), and as also demonstrated here for P. digitatum ΔpyrG for the first time (Supplementary Figure S4), pyr4 complementation completely restored pathogenicity, thus demonstrating the suitability of this genetic element as an auxotrophic selectable marker, as well as for (phyto)pathogenic fungi.

The luciferase reporter system has allowed us to functionally validate and characterize seven distinct promoters that have been incorporated as standard DNA parts to the FB toolbox. Among these are novel promoter sequences for which their functionality had never been validated before (PafpB, P7760, and Pef1α). These promoters were selected either for their interesting behavior in a P. digitatum-based transcriptome analysis (Ropero-Pérez et al., 2023) or because they are well-known promoters, such as the strong PpkiA or the maltose-responsive PglaA and PamyB, all of which have been extensively used in Aspergillus species (Storms et al., 2005; Oliveira et al., 2008; Song et al., 2018). Among the inducible promoters included in this study, PglaA and PxlnA have already been implemented in alternative Golden Gate-compatible collections (Polli et al., 2016; Mózsik et al., 2021), yet the validation of either basal or induced states has not been described in Penicillium species. Herein, a wide expression range was found for all tested promoters, from the highly expressed PpkiA and Pef1α promoters, with levels similar to those of the well-known Ppaf and PgpdA, to lower- or almost no-expressed promoters, such as P07760 and PafpB. The expression levels of almost all these promoters were reduced in P. chrysogenum compared to P. digitatum, except for the newly characterized Pef1α promoter from P digitatum, which showed similar values in both fungal backgrounds. This likely reflects a greater orthogonality in this promoter, which may be of preferable use to ensure strong expression in other fungal chassis. Inducibility of PglaA, PamyB, and PxlnA was also validated using the luciferase reporter system in both Penicillium species, showing different expression ranges both in the presence and absence of the inducer. This allows for multiple options for the custom design of future experiments, allowing promoters with lower background expression to choose, such as PxlnA, when basal expression needs to be almost completely avoided or to prioritize activation over background expression with promoters such as PamyB. Interestingly, PamyB induction in P. chrysogenum was similar to that shown in the industrial workhorse A. oryzae using β-glucuronidase (GUS) as the reporter. The expression of this promoter was reported to increase 10 times in a maltose-containing medium compared to glucose (Ozeki et al., 1996). This expression was slightly higher (1250 U/mg) (Tada et al., 1991) than that reported for PglaA (903 mg U/mg) (Hata et al., 1992), which correlates with our results shown in Figure 5B. Both PglaA and PamyB promoters are commonly applied for the production of different proteins of interest, such as human tissue plasminogen (Wiebe et al., 2001), bovine chymosin (Ohno et al., 2011), or synthetic human lysozyme (Tsuchiya et al., 1992a), which further shows their relevance in the field of fungal biotechnology. On the other hand, the use of the PxlnA promoter is very limited to date, being its ortholog, PxylP from P. chrysogenum more extensively used (Yap et al., 2022). Herein, we demonstrate the possibility of implementing this promoter in the Penicillium genus, with more modest induction levels than those of PglaA and PamyB but with the lack of basal expression in the absence of the inducer.

Unexpectedly, the Ppaf expression was found to significantly increase in the presence of maltose in P. chrysogenum and xylose in both P. chrysogenum and P. digitatum. This would suggest that (i) maltose and xylose themselves or any of the maltose/xylose catabolic intermediates serve as inducers for Ppaf or (ii) Ppaf expression is partially repressed by glucose, which can be attributed to the presence of carbon catabolite repression CREA motifs in the Ppaf sequence, as previously described (Marx et al., 1995). This repression is nevertheless almost completely lost when Ppaf is expressed in a different fungal chassis, such as P. digitatum, suggesting different regulatory mechanisms between both fungal species despite their phylogenetic proximity.

The activation of GB_SynP promoters in P. digitatum was addressed using the luciferase reporter system and the CRISPRa system included in the pAMA18.0_gRNA1 vector (Mózsik et al., 2021). The non-integrative nature of this pAMA1-based plasmid makes it possible to revert promoter activation upon plasmid loss in the absence of selection pressure (Garrigues et al., 2022). Additionally, this CRISPRa system provides a method to easily assay expression variations within the same background strain, either by testing different activation domains or inducible systems, or by analyzing the induction level under different culture conditions. The activation of 1xLuc and 2xLuc constructs, however, was not achieved in P. digitatum using the pAMA18.0_gRNA1 vector as signals of all but one 1xLuc re-transformant were on the same range as the basal signal of the reference strain. The expression of 1xLuc was, however, not different from the expression observed in 3xLuc strains in the absence of the CRISPRa system, which could mean that this expression is within the range of basal expression of the synthetic promoters. On the other hand, in the case of the 3xLuc construct, we did observe an increase of more than one order of magnitude in the presence of the CRISPRa system compared to the non-activation control for all 3xLuc re-transformants tested. These results indicate that the activation of these promoters in fungi requires the presence of at least three repetitions of the gRNA target, which highly differ from what was observed in plants, where one repetition of the target sequence for gRNA1 was sufficient to drive a significant increase of synthetic promoter expression (Moreno-Giménez et al., 2022). The discrepancies in the GB_SynP behavior between plants and fungi could be attributed to the differences in CRISPRa systems used in each organism. Although pAMA18.0_gRNA1 used for activation in P. digitatum comprised the dCas9 protein fused to the VPR activation domain, the dCasEV2.1 complex used for activation in Nicotiana benthamiana plants includes a dCas9 protein fused to an EDLL activation domain and an extra MS2 protein fused to the VPR domain that is able to recognize and bind the modified gRNA scaffold (Selma et al., 2019; Moreno-Giménez et al., 2022). Although VPR showed a major contribution in the activation as the expression levels dropped significantly when MS2 was fused to other activation domains (Moreno-Giménez et al., 2022), in fungi, a second activation component might be required to reach higher activation levels. Another explanation for low expression levels in fungi may reside in gRNA1 used to trigger the activation of GB_SynP promoters, which was originally designed for plants. Although no off-targets were found for gRNA1 in the P. digitatum CECT 20796 genome, the efficiency of this gRNA may not be optimal for this chassis, and therefore, a gRNA designed specifically for fungi might enhance the activated expression levels. Additionally, expression levels in fungi could also be enhanced by creating new A2 proximal promoter parts with more than three repetitions of the gRNA1 target sequence. Further optimization of GB_SynP promoters for filamentous fungi, following these guidelines, will be explored in the near future to better characterize this tool and its potential for wide-range expression of customizable synthetic promoters in filamentous fungi.

5 Conclusion

FungalBraid 2.0 aims to accelerate the development of fungal SynBio by the inclusion of a new repertoire of 27 domesticated DNA parts, resulting in functionally validated resistance and auxotrophic markers, as well as strong, inducible, and synthetic promoters for their contributory use by the fungal research community. With the newly expanded FB toolkit, a greater number of modular DNA assemblies are possible, exponentially increasing the possibilities for the study, development, and exploitation of filamentous fungi as cell biofactories.

Statements

Data availability statement

The datasets presented in this study can be found in the online repository DIGITAL CSIC (https://digital.csic.es/handle/10261/308066). DNA parts generated in this work can be found at: https://gbcloning.upv.es/ and https://www.addgene.org/ accession number(s) can be found in the Supplementary Material.

Author contributions

EM-G: methodology, investigation, formal analysis, visualization, and writing—original draft. MG: investigation and writing—review and editing. ZS: investigation. PM: conceptualization, supervision, funding acquisition, and writing—review and editing. LY: supervision and writing—review and editing. DO: resources, supervision and writing—review and editing. JM: conceptualization, methodology, supervision, funding acquisition, and writing—review and editing. SG: conceptualization, investigation, formal analysis, visualization, supervision, writing—original draft, and writing—review and editing. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by PROMETEO/2018/066 from “Conselleria d’Educació” (Generalitat Valenciana, Comunitat Valenciana, Spain), grant PID 2021-125858OB-100, and the Severo Ochoa Excellence Program CEX 2021-001189-S funded by MCIN/AEI/10.13039/501100011033 and by “ERDF A way of making Europe.” EM-G was the recipient of a predoctoral grant FPU18/02019 funded by MCIN/AEI/10.13039/501100011033 and by “ESF Investing in your future.” SG holds a Juan de la Cierva Incorporación grant (IJC 2020-042749-I) funded by MCIN/AEI/10.13039/501100011033 and the European Union NextGenerationEU/PRTR.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fbioe.2023.1222812/full#supplementary-material

References

  • 1

    AliS.BakkerenG. (2011). Introduction of large DNA inserts into the barley pathogenic fungus, Ustilago hordei, via recombined binary BAC vectors and agrobacterium-mediated transformation. Curr. Genet.57, 6373. 10.1007/s00294-010-0324-0

  • 2

    AustinB.HallR. M.TylerB. M. (1990). Optimized vectors and selection for transformation of Neurospora crassa and Aspergillus nidulans to bleomycin and phleomycin resistance. Gene93, 157162. 10.1016/0378-1119(90)90152-h

  • 3

    BallanceD. J.TurnerG. (1985). Development of a high-frequency transforming vector for Aspergillus nidulans. Gene36, 321331. 10.1016/0378-1119(85)90187-8

  • 4

    BennerS. A.SismourA. M. (2005). Synthetic biology. Nat. Rev. Genet.6, 533543. 10.1038/nrg1637

  • 5

    DahlmannT. A.TerfehrD.BeckerK.TeichertI. (2021). Golden gate vectors for efficient gene fusion and gene deletion in diverse filamentous fungi. Curr. Genet.67, 317330. 10.1007/s00294-020-01143-2

  • 6

    De GraaffL.Van Den BroeckH.VisserJ. (1992). Isolation and characterization of the Aspergillus niger pyruvate kinase gene. Curr. Genet.22, 2127. 10.1007/bf00351737

  • 7

    De GrootM. J. A.BundockP.HooykaasP. J. J.BeijersbergenA. G. M. (1998). Agrobacterium tumefaciens-mediated transformation of filamentous fungi. Nat. Biotechnol.16, 839842. 10.1038/nbt0998-839

  • 8

    De VriesR. P.PatyshakuliyevaA.GarriguesS.Agarwal-JansS. (2020). “The current biotechnological status and potential of plant and algal biomass degrading/modifying enzymes from ascomycete fungi,” in Grand challenges in fungal biotechnology. Editor NevalainenH. (Cham: Springer International Publishing), 81120.

  • 9

    DerntlC.Kiesenhofer DanielP., Mach robert, L., and mach-aigner astrid, R. (2015). novel strategies for genomic manipulation of Trichoderma reesei with the purpose of strain engineering. Appl. Environ. Microbiol.81, 63146323. 10.1128/aem.01545-15

  • 10

    DíezB.AlvarezE.CantoralJ. M.BarredoJ. L.MartínJ. F. (1987). Selection and characterization of pyrG mutants of Penicillium chrysogenum lacking orotidine-5′-phosphate decarboxylase and complementation by the pyr4 gene of Neurospora crassa. Curr. Genet.12, 277282. 10.1007/bf00435290

  • 11

    FierroF.VacaI.CastilloN. I.García-RicoR. O.ChávezR. (2022). Penicillium chrysogenum, a vintage model with a cutting-edge profile in biotechnology. Microorganisms10, 573. 10.3390/microorganisms10030573

  • 12

    FowlerT.BerkaR. M.WardM. (1990). Regulation of the glaA gene of Aspergillus niger. Curr. Genet.18, 537545. 10.1007/bf00327025

  • 13

    GandíaM.Moreno-GiménezE.Giner-LlorcaM.GarriguesS.Ropero-PérezC.LocascioA.et al (2022). Development of a FungalBraid Penicillium expansum-based expression system for the production of antifungal proteins in fungal biofactories. Microb. Biotechnol.15, 630647. 10.1111/1751-7915.14006

  • 14

    GarriguesS.GandíaM.PopaC.BoricsA.MarxF.CocaM.et al (2017). Efficient production and characterization of the novel and highly active antifungal protein AfpB from Penicillium digitatum. Sci. Rep.7, 14663. 10.1038/s41598-017-15277-w

  • 15

    GarriguesS.ManzanaresP.MarcosJ. F. (2022). Application of recyclable CRISPR/Cas9 tools for targeted genome editing in the postharvest pathogenic fungi Penicillium digitatum and Penicillium expansum. Curr. Genet.68, 515529. 10.1007/s00294-022-01236-0

  • 16

    González-CandelasL.AlamarS.Sánchez-TorresP.ZacaríasL.MarcosJ. F. (2010). A transcriptomic approach highlights induction of secondary metabolism in citrus fruit in response to Penicillium digitatum infection. BMC Plant Biol.10, 194. 10.1186/1471-2229-10-194

  • 17

    GruberF.VisserJ.KubicekC. P.De GraaffL. H. (1990). The development of a heterologous transformation system for the cellulolytic fungus Trichoderma reesei based on a pyrG-negative mutant strain. Curr. Genet.18, 7176. 10.1007/bf00321118

  • 18

    HamiltonC. M. (1997). A binary-BAC system for plant transformation with high-molecular-weight DNA. Gene200, 107116. 10.1016/s0378-1119(97)00388-0

  • 19

    HarriesE.GandíaM.CarmonaL.MarcosJ. F. (2015). The Penicillium digitatum protein O-mannosyltransferase Pmt2 is required for cell wall integrity, conidiogenesis, virulence and sensitivity to the antifungal peptide PAF26. Mol. Plant Pathol.16, 748761. 10.1111/mpp.12232

  • 20

    HataY.KitamotoK.GomiK.KumagaiC.TamuraG. (1992). Functional elements of the promoter region of the Aspergillus oryzae glaA gene encoding glucoamylase. Curr. Genet.22, 8591. 10.1007/bf00351466

  • 21

    HegedüsN.SiglC.ZadraI.PócsiI.MarxF. (2011). The paf gene product modulates asexual development in Penicillium chrysogenum. J. Basic Microbiol.51, 253262. 10.1002/jobm.201000321

  • 22

    Hernanz-KoersM.GandíaM.GarriguesS.ManzanaresP.YenushL.OrzaezD.et al (2018). FungalBraid: A GoldenBraid-based modular cloning platform for the assembly and exchange of DNA elements tailored to fungal synthetic biology. Fungal Genet. Biol.116, 5161. 10.1016/j.fgb.2018.04.010

  • 23

    HigashimuraN.HamadaA.OharaT.SakuraiS.ItoH.BanbaS. (2022). The target site of the novel fungicide quinofumelin, Pyricularia oryzae class II dihydroorotate dehydrogenase. J. Pestic. Sci.47, 190196. 10.1584/jpestics.d22-027

  • 24

    KellyJ. M.HynesM. J. (1985). Transformation of Aspergillus niger by the amdS gene of Aspergillus nidulans. EMBO J.4, 475479. 10.1002/j.1460-2075.1985.tb03653.x

  • 25

    KhangC. H.ParkS.-Y.RhoH.-S.LeeY.-H.KangS. (2007). “Filamentous fungi (Magnaporthe grisea and Fusarium oxysporum),” in Agrobacterium protocols volume 2. Editor WangK. (Totowa, NJ: Humana Press), 403420.

  • 26

    LiD.TangY.LinJ.CaiW. (2017). Methods for genetic transformation of filamentous fungi. Microb. Cell Fact.16, 168. 10.1186/s12934-017-0785-7

  • 27

    Marcet-HoubenM.BallesterA.-R.De La FuenteB.HarriesE.MarcosJ. F.González-CandelasL.et al (2012). Genome sequence of the necrotrophic fungus Penicillium digitatum, the main postharvest pathogen of citrus. BMC Genom13, 646. 10.1186/1471-2164-13-646

  • 28

    Martins-SantanaL.NoraL. C.Sanches-MedeirosA.LovateG. L.CassianoM. H. A.Silva-RochaR. (2018). Systems and synthetic biology approaches to engineer fungi for fine chemical production. Front. Bioeng. Biotechnol.6, 117. 10.3389/fbioe.2018.00117

  • 29

    MarxF.HaasH.ReindlM.StöfflerG.LottspeichF.RedlB. (1995). Cloning, structural organization and regulation of expression of the Penicillium chrysogenum paf gene encoding an abundantly secreted protein with antifungal activity. Gene167, 167171. 10.1016/0378-1119(95)00701-6

  • 30

    MeyerV.AndersenM. R.BrakhageA. A.BrausG. H.CaddickM. X.CairnsT. C.et al (2016). Current challenges of research on filamentous fungi in relation to human welfare and a sustainable bio-economy: A white paper. Fungal Biol. Biotechnol.3, 6. 10.1186/s40694-016-0024-8

  • 31

    MeyerV.BasenkoE. Y.BenzJ. P.BrausG. H.CaddickM. X.CsukaiM.et al (2020). Growing a circular economy with fungal biotechnology: A white paper. Fungal Biol. Biotechnol.7, 5. 10.1186/s40694-020-00095-z

  • 32

    Moreno-GiménezE.SelmaS.CalvacheC.OrzáezD. (2022). GB_SynP: A modular dCas9-regulated synthetic promoter collection for fine-tuned recombinant gene expression in plants. ACS Synth. Biol.11, 30373048. 10.1021/acssynbio.2c00238

  • 33

    MózsikL.IacovelliR.BovenbergR. a. L.DriessenA. J. M. (2022). Transcriptional activation of biosynthetic gene clusters in filamentous fungi. Front. Bioeng. Biotechnol.10, 901037. 10.3389/fbioe.2022.901037

  • 34

    MózsikL.PohlC.MeyerV.BovenbergR. a. L.NygårdY.DriessenA. J. M. (2021). Modular synthetic biology toolkit for filamentous fungi. ACS Synth. Biol.10, 28502861. 10.1021/acssynbio.1c00260

  • 35

    OhnoA.MaruyamaJ.-I.NemotoT.AriokaM.KitamotoK. (2011). A carrier fusion significantly induces unfolded protein response in heterologous protein production by Aspergillus oryzae. Appl. Microbiol. Biotechnol.92, 11971206. 10.1007/s00253-011-3487-9

  • 36

    OliveiraJ. M.Van Der VeenD.De GraaffL. H.QinL. (2008). Efficient cloning system for construction of gene silencing vectors in Aspergillus niger. Appl. Microbiol. Biotechnol.80, 917924. 10.1007/s00253-008-1640-x

  • 37

    OrejasM.MaccabeA. P.Pérez GonzálezJ. A.KumarS.RamónD. (1999). Carbon catabolite repression of the Aspergillus nidulans xlnA gene. Mol. Microbiol.31, 177184. 10.1046/j.1365-2958.1999.01157.x

  • 38

    OzekiK.KandaA.HamachiM.NunokawaY. (1996). Construction of a promoter probe vector autonomously maintained in Aspergillus and characterization of promoter regions derived from A. niger and A. oryzae genomes. Biosci. Biotechnol. Biochem.60, 383389. 10.1271/bbb.60.383

  • 39

    PalouL. (2014). “Penicillium digitatum, Penicillium italicum (green mold, blue mold),” in Postharvest decay. Editor Bautista-BañosS. (San Diego: Academic Press), 45102.

  • 40

    PolliF.MeijrinkB.BovenbergR. A. L.DriessenA. J. M. (2016). New promoters for strain engineering of Penicillium chrysogenum. Fungal Genet. Biol.89, 6271. 10.1016/j.fgb.2015.12.003

  • 41

    RobeyM. T.CaesarL. K.DrottM. T.KellerN. P.KelleherN. L. (2021). An interpreted atlas of biosynthetic gene clusters from 1,000 fungal genomes. PNAS118, e2020230118. 10.1073/pnas.2020230118

  • 42

    Ropero-PérezC.BolósB.Giner-LlorcaM.LocascioA.GarriguesS.GandíaM.et al (2023). Transcriptomic profile of Penicillium digitatum reveals novel aspects of the mode of action of the antifungal protein AfpB. Microbiol. Spectr.11 (3), e0484622. 10.1128/spectrum.04846-22

  • 43

    Sarrion-PerdigonesA.-VilarM.PalacíJ.CastelijnsB.FormentJ.ZiarsoloP.et al (2013). GoldenBraid 2.0: A comprehensive DNA assembly framework for plant synthetic biology. Plant Physiol.162, 16181631. 10.1104/pp.113.217661

  • 44

    Sarrion-PerdigonesA.FalconiE. E.ZandalinasS. I.JuárezP.Fernández-Del-CarmenA.GranellA.et al (2011). GoldenBraid: An iterative cloning system for standardized assembly of reusable genetic modules. PLOS ONE6, e21622. 10.1371/journal.pone.0021622

  • 45

    SelmaS.Bernabé-OrtsJ. M.Vazquez-VilarM.Diego-MartinB.AjenjoM.Garcia-CarpinteroV.et al (2019). Strong gene activation in plants with genome-wide specificity using a new orthogonal CRISPR/Cas9-Based programmable transcriptional activator. Plant Biotechnol. J.17, 17031705. 10.1111/pbi.13138

  • 46

    SelmaS.SanmartínN.Espinosa‐RuizA.GianoglioS.Lopez‐GresaM. P.Vázquez‐VilarM.et al (2022). Custom‐made design of metabolite composition in N. benthamiana leaves using CRISPR activators. Plant Biotechnol. J.20, 15781590. 10.1111/pbi.13834

  • 47

    SiglC.HandlerM.SprengerG.KürnsteinerH.ZadraI. (2010). A novel homologous dominant selection marker for genetic transformation of Penicillium chrysogenum: Overexpression of squalene epoxidase-encoding ergA. J. Biotechnol.150, 307311. 10.1016/j.jbiotec.2010.09.941

  • 48

    SondereggerC.GalgóczyL.GarriguesS.FizilÁ.BoricsA.ManzanaresP.et al (2016). A Penicillium chrysogenum-based expression system for the production of small, cysteine-rich antifungal proteins for structural and functional analyses. Microb. Cell Fact.15, 192. 10.1186/s12934-016-0586-4

  • 49

    SongL.OuedraogoJ. P.KolbuszM.NguyenT. T. M.TsangA. (2018). Efficient genome editing using tRNA promoter-driven CRISPR/Cas9 gRNA in Aspergillus niger. PLOS ONE13, e0202868. 10.1371/journal.pone.0202868

  • 50

    StormsR.ZhengY.LiH.SillaotsS.Martinez-PerezA.TsangA. (2005). Plasmid vectors for protein production, gene expression and molecular manipulations in Aspergillus niger. Plasmid53, 191204. 10.1016/j.plasmid.2004.10.001

  • 51

    TadaS.GomiK.KitamotoK.TakahashiK.TamuraG.HaraS. (1991). Construction of a fusion gene ComprisingThe taka-amylase A promoter and the Escherichia coli β-glucuronidase gene and analysis of its expression in Aspergillus oryzae. Mol. Genet. Genom.229, 301306. 10.1007/bf00272170

  • 52

    TakkenF. L.Van WijkR.MichielseC. B.HoutermanP. M.RamA. F.CornelissenB. J. (2004). A one-step method to convert vectors into binary vectors suited for agrobacterium-mediated transformation. Curr. Genet.45, 242248. 10.1007/s00294-003-0481-5

  • 53

    TsuchiyaK.TadaS.GomiK.KitamotoK.KumagaiC.JigamiY.et al (1992a). High level expression of the synthetic human lysozyme gene in Aspergillus oryzae. Appl. Microbiol. Biotechnol.38, 109114. 10.1007/bf00169428

  • 54

    TsuchiyaK.TadaS.GomiK.KitamotoK.KumagaiC.TamuraG. (1992b). Deletion analysis of the taka-amylase A gene promoter using a homologous transformation system in Aspergillus oryzae. Biosci. Biotechnol. Biochem.56, 18491853. 10.1271/bbb.56.1849

  • 55

    Vazquez-VilarM.GandíaM.García-CarpinteroV.MarquésE.Sarrion-PerdigonesA.YenushL.et al (2020). Multigene engineering by GoldenBraid cloning: From plants to filamentous fungi and beyond. Curr. Protoc. Mol. Biol.130, e116. 10.1002/cpmb.116

  • 56

    Vazquez-VilarM.Quijano-RubioA.Fernandez-Del-CarmenA.Sarrion-PerdigonesA.Ochoa-FernandezR.ZiarsoloP.et al (2017). GB3.0: A platform for plant bio-design that connects functional DNA elements with associated biological data. Nucleic Acids Res.45, 21962209. 10.1093/nar/gkw1326

  • 57

    WeberE.EnglerC.GruetznerR.WernerS.MarillonnetS. (2011). A modular cloning system for standardized assembly of multigene constructs. PLOS ONE6, e16765. 10.1371/journal.pone.0016765

  • 58

    WiebeM. G.KarandikarA.RobsonG. D.TrinciA. P. J.CandiaJ.-L. F.TrappeS.et al (2001). Production of tissue plasminogen activator (t-PA) in Aspergillus niger. Biotechnol. Bioeng.76, 164174. 10.1002/bit.1156

  • 59

    YapA.GlarcherI.MisslingerM.HaasH. (2022). Characterization and engineering of the xylose-inducible xylP promoter for use in mold fungal species. Metab. Eng. Commun.15, e00214. 10.1016/j.mec.2022.e00214

  • 60

    ZameitatE.FreymarkG.DietzC. D.LöfflerM.BölkerM. (2007). Functional expression of human dihydroorotate dehydrogenase (DHODH) in pyr4 mutants of Ustilago maydis allows target validation of DHODH inhibitors in vivo. Appl. Environ. Microbiol.73, 33713379. 10.1128/aem.02569-06

Summary

Keywords

fungal synthetic biology, GoldenBraid, promoters, selection markers, luciferase-based reporter system, CRISPR activation, filamentous fungi

Citation

Moreno-Giménez E, Gandía M, Sáez Z, Manzanares P, Yenush L, Orzáez D, Marcos JF and Garrigues S (2023) FungalBraid 2.0: expanding the synthetic biology toolbox for the biotechnological exploitation of filamentous fungi. Front. Bioeng. Biotechnol. 11:1222812. doi: 10.3389/fbioe.2023.1222812

Received

15 May 2023

Accepted

12 July 2023

Published

07 August 2023

Volume

11 - 2023

Edited by

Yaojun Tong, Shanghai Jiao Tong University, China

Reviewed by

Mario Andrea Marchisio, Tianjin University, China

Ahmad Bazli Ramzi, National University of Malaysia, Malaysia

Updates

Copyright

*Correspondence: Jose F. Marcos, ; Sandra Garrigues,

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics