Search

Aspergillus nidulans is an important model organism for eukaryotic biology and the reference for the section Nidulantes in comparative studies. In this study, we de novo sequenced the genomes of 25 species of this section. Whole-genome phylogeny of 34 Aspergillus species and Penicillium chrysogenum clarifies the position of clades inside section Nidulantes. Comparative genomics reveals a high genetic diversity between species with 684 up to 2433 unique protein families. Furthermore, we categorized 2118 secondary metabolite gene clusters (SMGC) into 603 families across Aspergilli, with at least 40 % of the families shared between Nidulantes species. Genetic dereplication of SMGC and subsequent synteny analysis provides evidence for horizontal gene transfer of a SMGC. Proteins that have been investigated in A. nidulans as well as its SMGC families are generally present in the section Nidulantes, supporting its role as model organism. The set of genes encoding plant biomass-related CAZymes is highly conserved in section Nidulantes, while there is remarkable diversity of organization of MAT-loci both within and between the different clades. This study provides a deeper understanding of the genomic conservation and diversity of this section and supports the position of A. nidulans as a reference species for cell biology.

Abstract 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.

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.

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.

ABSTRACT 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.

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.

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.

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.

Aspergillus fungi classified within the section Flavi include harmful and beneficial species. Here, Kjærbølling et al. analyse the genomes of 23 Flavi species, showing high genetic diversity and potential for synthesis of over 13,700 CAZymes and 1600 secondary metabolites.

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.

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.

An oleaginous yeast Rhodosporidium toruloides is a promising host for converting lignocellulosic biomass to bioproducts and biofuels. In this work, we performed multi-omics analysis of lignocellulosic carbon utilization in R. toruloides and reconstructed the genome-scale metabolic network of R. toruloides. High-quality metabolic network models for model organisms and orthologous protein mapping were used to build a draft metabolic network reconstruction. The reconstruction was manually curated to build a metabolic model using functional annotation and multi-omics data including transcriptomics, proteomics, metabolomics, and RB-TDNA sequencing. The multi-omics data and metabolic model were used to investigate R. toruloides metabolism including lipid accumulation and lignocellulosic carbon utilization. The developed metabolic model was validated against high-throughput growth phenotyping and gene fitness data, and further refined to resolve the inconsistencies between prediction and data. We believe that this is the most complete and accurate metabolic network model available for R. toruloides to date.

The widespread use of Next-Generation Sequencing has opened a new era in the study of biological systems by significantly increasing the catalog of fungal genomes sequences and identifying gene clusters for known secondary metabolites as well as novel cryptic ones. However, most of these clusters still need to be examined in detail to completely understand the pathway steps and the regulation of the biosynthesis of metabolites. Genome sequencing approach led to the identification of the biosynthetic genes cluster of ochratoxin A (OTA) in a number of producing fungal species. Ochratoxin A is a potent pentaketide nephrotoxin produced by Aspergillus and Penicillium species and found as widely contaminant in food, beverages and feed. The increasing availability of several new genome sequences of OTA producer species in JGI Mycocosm and/or GenBank databanks led us to analyze and update the gene cluster structure in 19 Aspergillus and 2 Penicillium OTA producing species, resulting in a well conserved organization of OTA core genes among the species. Furthermore, our comparative genome analyses evidenced the presence of an additional gene, previously undescribed, located between the polyketide and non-ribosomal synthase genes in the cluster of all the species analyzed. The presence of a SnoaL cyclase domain in the sequence of this gene supports its putative role in the polyketide cyclization reaction during the initial steps of the OTA biosynthesis pathway. The phylogenetic analysis showed a clustering of OTA SnoaL domains in accordance with the phylogeny of OTA producing species at species and section levels. The characterization of this new OTA gene, its putative role and its expression evidence in three important representative producing species, are reported here for the first time.

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.

Abstract Background Plant biomass is the most abundant carbon source for many fungal species. In the biobased industry fungi, are used to produce lignocellulolytic enzymes to degrade agricultural waste biomass. Here we evaluated if it would be possible to create an Aspergillus nidulans strain that releases, but does not metabolize hexoses from plant biomass. For this purpose, metabolic mutants were generated that were impaired in glycolysis, by using hexokinase (hxkA) and glucokinase (glkA) negative strains. To prevent repression of enzyme production due to the hexose accumulation, strains were generated that combined these mutations with a deletion in creA, the repressor involved in regulating preferential use of different carbon catabolic pathways. Results Phenotypic analysis revealed reduced growth for the hxkA1 glkA4 mutant on wheat bran. However, hexoses did not accumulate during growth of the mutants on wheat bran, suggesting that glucose metabolism is re-routed towards alternative carbon catabolic pathways. The creAΔ4 mutation in combination with preventing initial phosphorylation in glycolysis resulted in better growth than the hxkA/glkA mutant and an increased expression of pentose catabolic and pentose phosphate pathway genes. This indicates that the reduced ability to use hexoses as carbon sources created a shift towards the pentose fraction of wheat bran as a major carbon source to support growth. Conclusion Blocking the direct entry of hexoses to glycolysis activates alternative metabolic conversion of these sugars in A. nidulans during growth on plant biomass, but also upregulates conversion of other sugars, such as pentoses.

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.

Abstract Background The hydrolysis of biomass to simple sugars used for the production of biofuels in biorefineries requires the action of cellulolytic enzyme mixtures. During the last 50 years, the ascomycete Trichoderma reesei, the main source of industrial cellulase and hemicellulase cocktails, has been subjected to several rounds of classical mutagenesis with the aim to obtain higher production levels. During these random genetic events, strains unable to produce cellulases were generated. Here, whole genome sequencing and transcriptomic analyses of the cellulase-negative strain QM9978 were used for the identification of mutations underlying this cellulase-negative phenotype. Results Sequence comparison of the cellulase-negative strain QM9978 to the reference strain QM6a identified a total of 43 mutations, of which 33 were located either close to or in coding regions. From those, we identified 23 single-nucleotide variants, nine InDels, and one translocation. The translocation occurred between chromosomes V and VII, is located upstream of the putative transcription factor vib1, and abolishes its expression in QM9978 as detected during the transcriptomic analyses. Ectopic expression of vib1 under the control of its native promoter as well as overexpression of vib1 under the control of a strong constitutive promoter restored cellulase expression in QM9978, thus confirming that the translocation event is the reason for the cellulase-negative phenotype. Gene deletion of vib1 in the moderate producer strain QM9414 and in the high producer strain Rut-C30 reduced cellulase expression in both cases. Overexpression of vib1 in QM9414 and Rut-C30 had no effect on cellulase production, most likely because vib1 is already expressed at an optimal level under normal conditions. Conclusion We were able to establish a link between a chromosomal translocation in QM9978 and the cellulase-negative phenotype of the strain. We identified the transcription factor vib1 as a key regulator of cellulases in T. reesei whose expression is absent in QM9978. We propose that in T. reesei, as in Neurospora crassa, vib1 is involved in cellulase induction, although the exact mechanism remains to be elucidated. The data presented here show an example of a combined genome sequencing and transcriptomic approach to explain a specific trait, in this case the QM9978 cellulase-negative phenotype, and how it helps to better understand the mechanisms during cellulase gene regulation. When focusing on mutations on the single base-pair level, changes on the chromosome level can be easily overlooked and through this work we provide an example that stresses the importance of the big picture of the genomic landscape during analysis of sequencing data.

Filamentous fungi have been important as model organisms since the beginning of modern biological inquiry and have benefitted from open data since the earliest genetic maps were shared. From early origins in simple Mendelian genetics of mating types, parasexual genetics of colony colour, and the foundational demonstration of the segregation of a nutritional requirement, the contribution of research systems utilising filamentous fungi has spanned the biochemical genetics era, through the molecular genetics era, and now are at the very foundation of diverse omics approaches to research and development. Fungal model organisms have come from most major taxonomic groups although Ascomycete filamentous fungi have seen the most major sustained effort. In addition to the published material about filamentous fungi, shared molecular tools have found application in every area of fungal biology. Similarly, shared data has contributed to the success of model systems. The scale of data supporting research with filamentous fungi has grown by 10 to 12 orders of magnitude. From genetic to molecular maps, expression databases, and finally genome resources, the open and collaborative nature of the research communities has assured that the rising tide of data has lifted all of the research systems together.

ABSTRACT The yeast Yarrowia lipolytica undergoes a morphological transition from yeast-to-hyphal growth in response to environmental conditions. A forward genetic screen was used to identify mutants that reliably remain in the yeast phase, which were then assessed by whole-genome sequencing. All the smooth mutants identified, so named because of their colony morphology, exhibit independent loss of DNA at a repetitive locus made up of interspersed ribosomal DNA and short 10- to 40-mer telomere-like repeats. The loss of repetitive DNA is associated with downregulation of genes with stress response elements (5′-CCCCT-3′) and upregulation of genes with cell cycle box (5′-ACGCG-3′) motifs in their promoter region. The stress response element is bound by the transcription factor Msn2p in Saccharomyces cerevisiae. We confirmed that the Y. lipolytica msn2 (Ylmsn2) ortholog is required for hyphal growth and found that overexpression of Ylmsn2 enables hyphal growth in smooth strains. The cell cycle box is bound by the Mbp1p/Swi6p complex in S. cerevisiae to regulate G1-to-S phase progression. We found that overexpression of either the Ylmbp1 or Ylswi6 homologs decreased hyphal growth and that deletion of either Ylmbp1 or Ylswi6 promotes hyphal growth in smooth strains. A second forward genetic screen for reversion to hyphal growth was performed with the smooth-33 mutant to identify additional genetic factors regulating hyphal growth in Y. lipolytica. Thirteen of the mutants sequenced from this screen had coding mutations in five kinases, including the histidine kinases Ylchk1 and Ylnik1 and kinases of the high-osmolarity glycerol response (HOG) mitogen-activated protein (MAP) kinase cascade Ylssk2, Ylpbs2, and Ylhog1. Together, these results demonstrate that Y. lipolytica transitions to hyphal growth in response to stress through multiple signaling pathways. IMPORTANCE Many yeasts undergo a morphological transition from yeast-to-hyphal growth in response to environmental conditions. We used forward and reverse genetic techniques to identify genes regulating this transition in Yarrowia lipolytica. We confirmed that the transcription factor Ylmsn2 is required for the transition to hyphal growth and found that signaling by the histidine kinases Ylchk1 and Ylnik1 as well as the MAP kinases of the HOG pathway (Ylssk2, Ylpbs2, and Ylhog1) regulates the transition to hyphal growth. These results suggest that Y. lipolytica transitions to hyphal growth in response to stress through multiple kinase pathways. Intriguingly, we found that a repetitive portion of the genome containing telomere-like and rDNA repeats may be involved in the transition to hyphal growth, suggesting a link between this region and the general stress response.

The asexual filamentous fungus Aspergillus niger is an important industrial cell factory for citric acid production. In this study, we genetically characterized a UV-generated A. niger mutant that was originally isolated as a nonacidifying mutant, which is a desirable trait for industrial enzyme production. Physiological analysis showed that this mutant did not secrete large amounts of citric acid and oxalic acid, thus explaining the nonacidifying phenotype. As traditional complementation approaches to characterize the mutant genotype were unsuccessful, we used bulk segregant analysis in combination with high-throughput genome sequencing to identify the mutation responsible for the nonacidifying phenotype. Since A. niger has no sexual cycle, parasexual genetics was used to generate haploid segregants derived from diploids by loss of whole chromosomes. We found that the nonacidifying phenotype was caused by a point mutation in the laeA gene. LaeA encodes a putative methyltransferase-domain protein, which we show here to be required for citric acid production in an A. niger lab strain (N402) and in other citric acid production strains. The unexpected link between LaeA and citric acid production could provide new insights into the transcriptional control mechanisms related to citric acid production in A. niger. Interestingly, the secondary metabolite profile of a ΔlaeA strain differed from the wild-type strain, showing both decreased and increased metabolite levels, indicating that LaeA is also involved in regulating the production of secondary metabolites. Finally, we show that our systems genetics approach is a powerful tool to identify trait mutations.

ABSTRACT The yeast Yarrowia lipolytica is a potent accumulator of lipids, and lipogenesis in this organism can be influenced by a variety of factors, such as genetics and environmental conditions. Using a multifactorial study, we elucidated the effects of both genetic and environmental factors on regulation of lipogenesis in Y. lipolytica and identified how two opposite regulatory states both result in lipid accumulation. This study involved comparison of a strain overexpressing diacylglycerol acyltransferase (DGA1) with a control strain grown under either nitrogen or carbon limitation conditions. A strong correlation was observed between the responses on the transcript and protein levels. Combination of DGA1 overexpression with nitrogen limitation resulted in a high level of lipid accumulation accompanied by downregulation of several amino acid biosynthetic pathways, including that of leucine in particular, and these changes were further correlated with a decrease in metabolic fluxes. This downregulation was supported by the measured decrease in the level of 2-isopropylmalate, an intermediate of leucine biosynthesis. Combining the multi-omics data with putative transcription factor binding motifs uncovered a contradictory role for TORC1 in controlling lipid accumulation, likely mediated through 2-isopropylmalate and a Leu3-like transcription factor. IMPORTANCE The ubiquitous metabolism of lipids involves refined regulation, and an enriched understanding of this regulation would have wide implications. Various factors can influence lipid metabolism, including the environment and genetics. We demonstrated, using a multi-omics and multifactorial experimental setup, that multiple factors affect lipid accumulation in the yeast Yarrowia lipolytica. Using integrative analysis, we identified novel interactions between nutrient restriction and genetic factors involving regulators that are highly conserved among eukaryotes. Given that lipid metabolism is involved in many diseases but is also vital to the development of microbial cell factories that can provide us with sustainable fuels and oleochemicals, we envision that our report introduces foundational work to further unravel the regulation of lipid accumulation in eukaryal cells.

ABSTRACT Fungi accumulate lipids in a manner dependent on the quantity and quality of the nitrogen source on which they are growing. In the oleaginous yeast Yarrowia lipolytica, growth on a complex source of nitrogen enables rapid growth and limited accumulation of neutral lipids, while growth on a simple nitrogen source promotes lipid accumulation in large lipid droplets. Here we examined the roles of nitrogen catabolite repression and its regulation by GATA zinc finger transcription factors on lipid metabolism in Y. lipolytica. Deletion of the GATA transcription factor genes gzf3 and gzf2 resulted in nitrogen source-specific growth defects and greater accumulation of lipids when the cells were growing on a simple nitrogen source. Deletion of gzf1, which is most similar to activators of genes repressed by nitrogen catabolite repression in filamentous ascomycetes, did not affect growth on the nitrogen sources tested. We examined gene expression of wild-type and GATA transcription factor mutants on simple and complex nitrogen sources and found that expression of enzymes involved in malate metabolism, beta-oxidation, and ammonia utilization are strongly upregulated on a simple nitrogen source. Deletion of gzf3 results in overexpression of genes with GATAA sites in their promoters, suggesting that it acts as a repressor, while gzf2 is required for expression of ammonia utilization genes but does not grossly affect the transcription level of genes predicted to be controlled by nitrogen catabolite repression. Both GATA transcription factor mutants exhibit decreased expression of genes controlled by carbon catabolite repression via the repressor mig1, including genes for beta-oxidation, highlighting the complex interplay between regulation of carbon, nitrogen, and lipid metabolism. IMPORTANCE Nitrogen source is commonly used to control lipid production in industrial fungi. Here we identified regulators of nitrogen catabolite repression in the oleaginous yeast Y. lipolytica to determine how the nitrogen source regulates lipid metabolism. We show that disruption of both activators and repressors of nitrogen catabolite repression leads to increased lipid accumulation via activation of carbon catabolite repression through an as yet uncharacterized method.

Cochliobolus heterostrophus race T, causal agent of southern corn leaf blight, requires T-toxin (a family of C35 to C49 polyketides) for high virulence on T-cytoplasm maize. Production of T-toxin is controlled by two unlinked loci, Tox1A and Tox1B, carried on 1.2 Mb of DNA not found in race O, a mildly virulent form of the fungus that does not produce T-toxin, or in any other Cochliobolus spp. or closely related fungus. PKS1, a polyketide synthase (PKS)-encoding gene at Tox1A, and DEC1, a decarboxylase-encoding gene at Tox1B, are necessary for T-toxin production. Although there is evidence that additional genes are required for Ttoxin production, efforts to clone them have been frustrated because the genes are located in highly repeated, A+T-rich DNA. To overcome this difficulty, ligation specificity-based expression analysis display (LEAD), a comparative amplified fragment length polymorphism/gel fractionation/capillary sequencing procedure, was applied to cDNAs from a near-isogenic pair of race T (Tox1+) and race O (Tox1-) strains. This led to discovery of PKS2, a second PKS-encoding gene that maps at Tox1A and is required for both Ttoxin biosynthesis and high virulence to maize. Thus, the carbon chain of each T-toxin family member likely is assembled by action of two PKSs, which produce two polyketides, one of which may act as the starter unit for biosynthesis of the mature T-toxin molecule.