Genomic Analysis of Aspergillus Section Terrei Reveals a High Potential in Secondary Metabolite Production and Plant Biomass Degradation.

Aspergillus terreus has attracted interest due to its application in industrial biotechnology, particularly for the production of itaconic acid and bioactive secondary metabolites. As related species also seem to possess a prosperous secondary metabolism, they are of high interest for genome mining and exploitation. Here, we present draft genome sequences for six species from Aspergillus section Terrei and one species from Aspergillus section Nidulantes. Whole-genome phylogeny confirmed that section Terrei is monophyletic. Genome analyses identified between 70 and 108 key secondary metabolism genes in each of the genomes of section Terrei, the highest rate found in the genus Aspergillus so far. The respective enzymes fall into 167 distinct families with most of them corresponding to potentially unique compounds or compound families. Moreover, 53% of the families were only found in a single species, which supports the suitability of species from section Terrei for further genome mining. Intriguingly, this analysis, combined with heterologous gene expression and metabolite identification, suggested that species from section Terrei use a strategy for UV protection different to other species from the genus Aspergillus. Section Terrei contains a complete plant polysaccharide degrading potential and an even higher cellulolytic potential than other Aspergilli, possibly facilitating additional applications for these species in biotechnology.

Genome-scale model development and genomic sequencing of the oleaginous clade Lipomyces.

The Lipomyces clade contains oleaginous yeast species with advantageous metabolic features for biochemical and biofuel production. Limited knowledge about the metabolic networks of the species and limited tools for genetic engineering have led to a relatively small amount of research on the microbes. Here, a genome-scale metabolic model (GSM) of Lipomyces starkeyi NRRL Y-11557 was built using orthologous protein mappings to model yeast species. Phenotypic growth assays were used to validate the GSM (66% accuracy) and indicated that NRRL Y-11557 utilized diverse carbohydrates but had more limited catabolism of organic acids. The final GSM contained 2,193 reactions, 1,909 metabolites, and 996 genes and was thus named iLst996. The model contained 96 of the annotated carbohydrate-active enzymes. iLst996 predicted a flux distribution in line with oleaginous yeast measurements and was utilized to predict theoretical lipid yields. Twenty-five other yeasts in the Lipomyces clade were then genome sequenced and annotated. Sixteen of the Lipomyces species had orthologs for more than 97% of the iLst996 genes, demonstrating the usefulness of iLst996 as a broad GSM for Lipomyces metabolism. Pathways that diverged from iLst996 mainly revolved around alternate carbon metabolism, with ortholog groups excluding NRRL Y-11557 annotated to be involved in transport, glycerolipid, and starch metabolism, among others. Overall, this study provides a useful modeling tool and data for analyzing and understanding Lipomyces species metabolism and will assist further engineering efforts in Lipomyces.

Unveiling a classical mutant in the context of the GH3 ß-glucosidase family in Neurospora crassa.

Classical fungal mutant strains obtained by mutagenesis have helped to elucidate fundamental metabolic pathways in the past. In the filamentous fungus Neurospora crassa, the gluc-1 strain was isolated long ago and characterized by its low level of ß-glucosidase activity, which is essential for the degradation of cellulose, the most abundant biopolymer on Earth and the main polymeric component of the plant cell wall. Based on genomic resequencing, we hypothesized that the causative mutation resides in the ß-glucosidase gene gh3-3 (bgl6, NCU08755). In this work, growth patterns, enzymatic activities and sugar utilization rates were analyzed in several mutant and overexpression strains related to gluc-1 and gh3-3. In addition, different mutants affected in the degradation and transport of cellobiose were analyzed. While overexpression of gh3-3 led to the recovery of ß-glucosidase activity in the gluc-1 mutant, as well as normal utilization of cellobiose, the full gene deletion strain Δgh3-3 was found to behave differently than gluc-1 with lower secreted ß-glucosidase activity, indicating a dominant role of the amino acid substitution in the point mutated gh3-3 gene of gluc-1. Our results furthermore confirm that GH3-3 is the major extracellular ß-glucosidase in N. crassa and demonstrate that the two cellodextrin transporters CDT-1 and CDT-2 are essential for growth on cellobiose when the three main N. crassa ß-glucosidases are absent. Overall, these findings provide valuable insight into the mechanisms of cellulose utilization in filamentous fungi, being an essential step in the efficient production of biorefinable sugars from agricultural and forestry plant biomass.

T-Toxin Virulence Genes: Unconnected Dots in a Sea of Repeats.

In 1970, the Southern Corn Leaf Blight epidemic ravaged U.S. fields to great economic loss. The outbreak was caused by never-before-seen, supervirulent, Race T of the fungus Cochliobolus heterostrophus. The functional difference between Race T and O, the previously known, far less aggressive strain, is production of T-toxin, a host-selective polyketide. Supervirulence is associated with ~1 Mb of Race T-specific DNA; only a fraction encodes T-toxin biosynthetic genes (Tox1). Tox1 is genetically and physically complex, with unlinked loci (Tox1A, Tox1B) genetically inseparable from breakpoints of a Race O reciprocal translocation that generated hybrid Race T chromosomes. Previously, we identified 10 genes for T-toxin biosynthesis. Unfortunately, high-depth, short-read sequencing placed these genes on four small, unconnected scaffolds surrounded by repeated A+T rich sequence, concealing context. To sort out Tox1 topology and pinpoint the hypothetical Race O translocation breakpoints corresponding to Race T-specific insertions, we undertook PacBio long-read sequencing which revealed Tox1 gene arrangement and the breakpoints. Six Tox1A genes are arranged as three small islands in a Race T-specific sea (~634 kb) of repeats. Four Tox1B genes are linked, on a large loop of Race T-specific DNA (~210 kb). The race O breakpoints are short sequences of race O-specific DNA; corresponding positions in race T are large insertions of race T-specific, A+T rich DNA, often with similarity to transposable (predominantly Gypsy) elements. Nearby, are 'Voyager Starship' elements and DUF proteins. These elements may have facilitated Tox1 integration into progenitor Race O and promoted large scale recombination resulting in race T. IMPORTANCE In 1970 a corn disease epidemic ravaged fields in the United States to great economic loss. The outbreak was caused by a never-before seen, supervirulent strain of the fungal pathogen Cochliobolus heterostrophus. This was a plant disease epidemic, however, the current COVID-19 pandemic of humans is a stark reminder that novel, highly virulent, pathogens evolve with devastating consequences, no matter what the host-animal, plant, or other organism. Long read DNA sequencing technology allowed in depth structural comparisons between the sole, previously known, much less aggressive, version of the pathogen and the supervirulent version and revealed, in meticulous detail, the structure of the unique virulence-causing DNA. These data are foundational for future analysis of mechanisms of DNA acquisition from a foreign source.

Identifying the gluc-1 and gluc-2 mutations in Neurospora crassa by genome resequencing.

Genome resequencing is an efficient strategy for associating mutant phenotypes with physical genomic loci (Baker 2009). A pilot study of this approach demonstrated that the Neurospora crassa genetic map was critical in narrowing the possible candidate mutations in a strain to a small number in a limited, defined region of the genome (McCluskey et al. 2011). In this study, we utilize a resequencing strategy to identify the mutations underlying the gluc-1 and gluc-2 genes in N. crassa.

Tuning a high performing multiplexed-CRISPRi Pseudomonas putida strain to further enhance indigoidine production.

In this study, a 14-gene edited Pseudomonas putida KT2440 strain for heterologous indigoidine production was examined using three distinct omic datasets. Transcriptomic data indicated that CRISPR/dCpf1-interference (CRISPRi) mediated multiplex repression caused global gene expression changes, implying potential undesirable changes in metabolic flux. 13C-metabolic flux analysis (13C-MFA) revealed that the core P. putida flux network after CRISPRi repression was conserved, with moderate reduction of TCA cycle and pyruvate shunt activity along with glyoxylate shunt activation during glucose catabolism. Metabolomic results identified a change in intracellular TCA metabolites and extracellular metabolite secretion profiles (sugars and succinate overflow) in the engineered strains. These omic analyses guided further strain engineering, with a random mutagenesis screen first identifying an optimal ribosome binding site (RBS) for Cpf1 that enabled stronger product-substrate pairing (1.6-fold increase). Then, deletion strains were constructed with excision of the PHA operon (ΔphaAZC-IID) resulting in a 2.2-fold increase in indigoidine titer over the optimized Cpf1-RBS construct at the end of the growth phase (∼6 h). The maximum indigoidine titer (at 72 h) in the ΔphaAZC-IID strain had a 1.5-fold and 1.8-fold increase compared to the optimized Cpf1-RBS construct and the original strain, respectively. Overall, this study demonstrated that integration of omic data types is essential for understanding responses to complex metabolic engineering designs and directly quantified the effect of such modifications on central metabolism.

Itaconic acid production is regulated by LaeA in Aspergillus pseudoterreus.

The global regulator LaeA controls secondary metabolism in diverse Aspergillus species. Here we explored its role in regulation of itaconic acid production in Aspergillus pseudoterreus. To understand its role in regulating metabolism, we deleted and overexpressed laeA, and assessed the transcriptome, proteome, and secreted metabolome prior to and during initiation of phosphate limitation induced itaconic acid production. We found that secondary metabolite clusters, including the itaconic acid biosynthetic gene cluster, are regulated by laeA and that laeA is required for high yield production of itaconic acid. Overexpression of LaeA improves itaconic acid yield at the expense of biomass by increasing the expression of key biosynthetic pathway enzymes and attenuating the expression of genes involved in phosphate acquisition and scavenging. Increased yield was observed in optimized conditions as well as conditions containing excess nutrients that may be present in inexpensive sugar containing feedstocks such as excess phosphate or complex nutrient sources. This suggests that global regulators of metabolism may be useful targets for engineering metabolic flux that is robust to environmental heterogeneity.

Phenotype to genotype in Neurospora crassa: Association of the scumbo phenotype with mutations in the gene encoding ceramide C9-methyltransferase.

Using a legacy of genetic mutants of Neurospora crassa, paired with resequencing efforts through JGI, we have identified the gene responsible for the 'scumbo' mutant. This early morphological mutant was described as "Irregular flat, spreading growth with knobby protrusions and abnormal conidiation, but no free conidia. Mycelium usually appears yellowish rather than orange. Female fertile." (Perkins, Radford et al. 2000). Our further investigation has found new insights into the identity and associated functions of scumbo as a ceramide C9 methyltransferase, previously annotated as "similar to cyclopropane-fatty-acyl-phospholipidsynthase", encoded by the gene NCU07859. This enzyme performs a fungal-specific methyl modification of glycosyl-ceramides and has implications for membrane homeostasis and hyphal polarity in filamentous fungi.

Online Hydrophilic Interaction Chromatography (HILIC) Enhanced Top-Down Mass Spectrometry Characterization of the SARS-CoV-2 Spike Receptor-Binding Domain.

SARS-CoV-2 cellular infection is mediated by the heavily glycosylated spike protein. Recombinant versions of the spike protein and the receptor-binding domain (RBD) are necessary for seropositivity assays and can potentially serve as vaccines against viral infection. RBD plays key roles in the spike protein's structure and function, and thus, comprehensive characterization of recombinant RBD is critically important for biopharmaceutical applications. Liquid chromatography coupled to mass spectrometry has been widely used to characterize post-translational modifications in proteins, including glycosylation. Most studies of RBDs were performed at the proteolytic peptide (bottom-up proteomics) or released glycan level because of the technical challenges in resolving highly heterogeneous glycans at the intact protein level. Herein, we evaluated several online separation techniques: (1) C2 reverse-phase liquid chromatography (RPLC), (2) capillary zone electrophoresis (CZE), and (3) acrylamide-based monolithic hydrophilic interaction chromatography (HILIC) to separate intact recombinant RBDs with varying combinations of glycosylations (glycoforms) for top-down mass spectrometry (MS). Within the conditions we explored, the HILIC method was superior to RPLC and CZE at separating RBD glycoforms, which differ significantly in neutral glycan groups. In addition, our top-down analysis readily captured unexpected modifications (e.g., cysteinylation and N-terminal sequence variation) and low abundance, heavily glycosylated proteoforms that may be missed by using glycopeptide data alone. The HILIC top-down MS platform holds great potential in resolving heterogeneous glycoproteins for facile comparison of biosimilars in quality control applications.

Major involvement of two laccase genes in conidial pigment biosynthesis in Aspergillus oryzae.

Wild-type strains of Aspergillus oryzae develop yellow, yellow-green, green, or brown conidia. Previous reports suggested that the conidiation initiates with the biosynthesis of a yellow pigment YWA1 from acetyl-CoA by a polyketide synthase encoded by wA (AO090102000545). This is followed by the conversion to other pigment by a laccase encoded by yA (AO090011000755). Based on orthologous pathways in other Aspergilli, it is reasonable to hypothesize that in addition to yA, AO090102000546 encoding laccase and AO090005000332 encoding Ayg1-like hydrolase play a role in A. oryzae conidial pigment biosynthesis. However, the involvement of these two genes in conidial pigmentation remains unclear. In this study, we tested this hypothesis by assessing the conidial colors of both disruption and overexpression mutants to verify whether AO090102000546 and AO090005000332 were associated with the conidial pigmentation. Observation of single, double, and triple disruptants of these three genes suggested that conidial pigments were synthesized by two laccase genes, AO090011000755 and AO090102000546, whereas Ayg1-like hydrolase gene AO090005000332 was proven to have no obvious association with the synthesis. This was corroborated by observing the phenotype of each overexpression mutant. Interestingly, AO090005000332 overexpression mutant produced smoky yellow-green conidia, different from the wild-type strain. Thus, the AO090005000332-encoded protein is likely to maintain the enzymatic activity. However, the expression level was observed to be one-third of that of AO090102000546 and one-seventh of that of AO090011000755. Consequently, apparent lack of obvious contribution of AO090005000332 to conidial pigmentation could be attributed to its low expression level. Expression analysis indicated similar profiles in several wild-type strains. KEY POINTS: • Conidial pigment biosynthesis after YWA1 mainly involves two laccases in A. oryzae. • Ayg1-like hydrolase in A. oryzae is not obviously involved in conidial pigmentation. • Conidial color is deemed dependent on expression level of two laccases and hydrolase.

Evidence of the Involvement of a Cyclase Gene in the Biosynthesis of Ochratoxin A in Aspergillus carbonarius.

Ochratoxin A (OTA) is a well-known mycotoxin with wide distribution in food and feed. Fungal genome sequencing has great utility for identifying secondary metabolites gene clusters for known and novel compounds. A comparative analysis of the OTA-biosynthetic cluster in A. steynii, A. westerdijkiae, A. niger, A. carbonarius, and P. nordicum has revealed a high synteny in OTA cluster organization in five structural genes (otaA, otaB, ota, otaR1, and otaD). Moreover, a recent detailed comparative genome analysis of Aspergilli OTA producers led to the identification of a cyclase gene, otaY, located in the OTA cluster between the otaA and otaB genes, encoding for a predicted protein with high similarity to SnoaLs domain. These proteins have been shown to catalyze ring closure steps in the biosynthesis of polyketide antibiotics produced in Streptomyces. In the present study, we demonstrated an upregulation of the cyclase gene in A. carbonarius under OTA permissive conditions, consistent with the expression trends of the other OTA cluster genes and their role in OTA biosynthesis by complete gene deletion. Our results pointed out the involvement of a cyclase gene in OTA biosynthetic pathway for the first time. They represent a step forward in the understanding of the molecular basis of OTA biosynthesis in A. carbonarius.

Bayesian Inference for Integrating Yarrowia lipolytica Multiomics Datasets with Metabolic Modeling.

Optimizing the metabolism of microbial cell factories for yields and titers is a critical step for economically viable production of bioproducts and biofuels. In this process, tuning the expression of individual enzymes to obtain the desired pathway flux is a challenging step, in which data from separate multiomics techniques must be integrated with existing biological knowledge to determine where changes should be made. Following a design-build-test-learn strategy, building on recent advances in Bayesian metabolic control analysis, we identify key enzymes in the oleaginous yeast Yarrowia lipolytica that correlate with the production of itaconate by integrating a metabolic model with multiomics measurements. To this extent, we quantify the uncertainty for a variety of key parameters, known as flux control coefficients (FCCs), needed to improve the bioproduction of target metabolites and statistically obtain key correlations between the measured enzymes and boundary flux. Based on the top five significant FCCs and five correlated enzymes, our results show phosphoglycerate mutase, acetyl-CoA synthetase (ACSm), carbonic anhydrase (HCO3E), pyrophosphatase (PPAm), and homoserine dehydrogenase (HSDxi) enzymes in rate-limiting reactions that can lead to increased itaconic acid production.

The F-box protein gene exo-1 is a target for reverse engineering enzyme hypersecretion in filamentous fungi.

Carbohydrate active enzymes (CAZymes) are vital for the lignocellulose-based biorefinery. The development of hypersecreting fungal protein production hosts is therefore a major aim for both academia and industry. However, despite advances in our understanding of their regulation, the number of promising candidate genes for targeted strain engineering remains limited. Here, we resequenced the genome of the classical hypersecreting Neurospora crassa mutant exo-1 and identified the causative point of mutation to reside in the F-box protein-encoding gene, NCU09899. The corresponding deletion strain displayed amylase and invertase activities exceeding those of the carbon catabolite derepressed strain Δcre-1, while glucose repression was still mostly functional in Δexo-1 Surprisingly, RNA sequencing revealed that while plant cell wall degradation genes are broadly misexpressed in Δexo-1, only a small fraction of CAZyme genes and sugar transporters are up-regulated, indicating that EXO-1 affects specific regulatory factors. Aiming to elucidate the underlying mechanism of enzyme hypersecretion, we found the high secretion of amylases and invertase in Δexo-1 to be completely dependent on the transcriptional regulator COL-26. Furthermore, misregulation of COL-26, CRE-1, and cellular carbon and nitrogen metabolism was confirmed by proteomics. Finally, we successfully transferred the hypersecretion trait of the exo-1 disruption by reverse engineering into the industrially deployed fungus Myceliophthora thermophila using CRISPR-Cas9. Our identification of an important F-box protein demonstrates the strength of classical mutants combined with next-generation sequencing to uncover unanticipated candidates for engineering. These data contribute to a more complete understanding of CAZyme regulation and will facilitate targeted engineering of hypersecretion in further organisms of interest.

Strain Construction for Intracellular Metabolic Pathway Localization in Y. lipolytica.

Pathway localization by fluorophore or epitope tagging can be accomplished through a multi-staged DNA construct and confirmation process, to generate a series of successfully tagged protein targets. Prerequisite conditions for this process in Y. lipolytica are auxotrophic selection (leu2 or ura3), impaired non-homologous end joining by deletion or impairment of ku70, and plasmids or gene pieces for epitope-selection cassette construction. The general approach for gene tagging can work for C- or N-terminal tags. Gene overexpression from an episomal plasmid can be accomplished through transcript amplification and cloning. C-terminal tagging allows expression of a gene-GFP fusion to be regulated from the endogenous promoter. The epitope-selection cassette also includes a constitutive or highly expressed promoter driving the auxotrophic or other selectable marker gene such as one conferring antifungal or antibiotic resistance. Strains for pathway localization utilize overlap PCR, PEG-based transformation, and a fast DNA preparation for rapid colony screening. Successful transformants can be used for pathway localization and condition-specific response analysis.

Integration of Proteomics and Metabolomics Into the Design, Build, Test, Learn Cycle to Improve 3-Hydroxypropionic Acid Production in Aspergillus pseudoterreus.

Biological engineering of microorganisms to produce value-added chemicals is a promising route to sustainable manufacturing. However, overproduction of metabolic intermediates at high titer, rate, and yield from inexpensive substrates is challenging in non-model systems where limited information is available regarding metabolic flux and its control in production conditions. Integrated multi-omic analyses of engineered strains offers an in-depth look at metabolites and proteins directly involved in growth and production of target and non-target bioproducts. Here we applied multi-omic analyses to overproduction of the polymer precursor 3-hydroxypropionic acid (3HP) in the filamentous fungus Aspergillus pseudoterreus. A synthetic pathway consisting of aspartate decarboxylase, beta-alanine pyruvate transaminase, and 3HP dehydrogenase was designed and built for A. pseudoterreus. Strains with single- and multi-copy integration events were isolated and multi-omics analysis consisting of intracellular and extracellular metabolomics and targeted and global proteomics was used to interrogate the strains in shake-flask and bioreactor conditions. Production of a variety of co-products (organic acids and glycerol) and oxidative degradation of 3HP were identified as metabolic pathways competing with 3HP production. Intracellular accumulation of nitrogen as 2,4-diaminobutanoate was identified as an off-target nitrogen sink that may also limit flux through the engineered 3HP pathway. Elimination of the high-expression oxidative 3HP degradation pathway by deletion of a putative malonate semialdehyde dehydrogenase improved the yield of 3HP by 3.4 × after 10 days in shake-flask culture. This is the first report of 3HP production in a filamentous fungus amenable to industrial scale biomanufacturing of organic acids at high titer and low pH.

Re-routing of Sugar Catabolism Provides a Better Insight Into Fungal Flexibility in Using Plant Biomass-Derived Monomers as Substrates.

The filamentous ascomycete Aspergillus niger has received increasing interest as a cell factory, being able to efficiently degrade plant cell wall polysaccharides as well as having an extensive metabolism to convert the released monosaccharides into value added compounds. The pentoses D-xylose and L-arabinose are the most abundant monosaccharides in plant biomass after the hexose D-glucose, being major constituents of xylan, pectin and xyloglucan. In this study, the influence of selected pentose catabolic pathway (PCP) deletion strains on growth on plant biomass and re-routing of sugar catabolism was addressed to gain a better understanding of the flexibility of this fungus in using plant biomass-derived monomers. The transcriptome, metabolome and proteome response of three PCP mutant strains, ΔlarAΔxyrAΔxyrB, ΔladAΔxdhAΔsdhA and ΔxkiA, grown on wheat bran (WB) and sugar beet pulp (SBP), was evaluated. Our results showed that despite the absolute impact of these PCP mutations on pure pentose sugars, they are not as critical for growth of A. niger on more complex biomass substrates, such as WB and SBP. However, significant phenotypic variation was observed between the two biomass substrates, but also between the different PCP mutants. This shows that the high sugar heterogeneity of these substrates in combination with the high complexity and adaptability of the fungal sugar metabolism allow for activation of alternative strategies to support growth.

A comparative genomics study of 23 Aspergillus species from section Flavi.

Section Flavi encompasses both harmful and beneficial Aspergillus species, such as Aspergillus oryzae, used in food fermentation and enzyme production, and Aspergillus flavus, food spoiler and mycotoxin producer. Here, we sequence 19 genomes spanning section Flavi and compare 31 fungal genomes including 23 Flavi species. We reassess their phylogenetic relationships and show that the closest relative of A. oryzae is not A. flavus, but A. minisclerotigenes or A. aflatoxiformans and identify high genome diversity, especially in sub-telomeric regions. We predict abundant CAZymes (598 per species) and prolific secondary metabolite gene clusters (73 per species) in section Flavi. However, the observed phenotypes (growth characteristics, polysaccharide degradation) do not necessarily correlate with inferences made from the predicted CAZyme content. Our work, including genomic analyses, phenotypic assays, and identification of secondary metabolites, highlights the genetic and metabolic diversity within section Flavi.

High-Throughput Large-Scale Targeted Proteomics Assays for Quantifying Pathway Proteins in Pseudomonas putida KT2440.

Targeted proteomics is a mass spectrometry-based protein quantification technique with high sensitivity, accuracy, and reproducibility. As a key component in the multi-omics toolbox of systems biology, targeted liquid chromatography-selected reaction monitoring (LC-SRM) measurements are critical for enzyme and pathway identification and design in metabolic engineering. To fulfill the increasing need for analyzing large sample sets with faster turnaround time in systems biology, high-throughput LC-SRM is greatly needed. Even though nanoflow LC-SRM has better sensitivity, it lacks the speed offered by microflow LC-SRM. Recent advancements in mass spectrometry instrumentation significantly enhance the scan speed and sensitivity of LC-SRM, thereby creating opportunities for applying the high speed of microflow LC-SRM without losing peptide multiplexing power or sacrificing sensitivity. Here, we studied the performance of microflow LC-SRM relative to nanoflow LC-SRM by monitoring 339 peptides representing 132 enzymes in Pseudomonas putida KT2440 grown on various carbon sources. The results from the two LC-SRM platforms are highly correlated. In addition, the response curve study of 248 peptides demonstrates that microflow LC-SRM has comparable sensitivity for the majority of detected peptides and better mass spectrometry signal and chromatography stability than nanoflow LC-SRM.