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.

Simple and Effective Squash-PCR for Rapid Genotyping of Industrial Microalgae.

Microalgae are recognized for their versatility in providing renewable energy, biopharmaceuticals, and nutraceuticals, attributed to their sustainable, renewable, and cost-effective nature. Genetic engineering has proven highly effective in enhancing microalgae production. PCR-based genotyping is the primary method for screening genetically transformed microalgae cells. Recently, we developed a novel PCR method, namely Squash-PCR, and employed it for the molecular analysis of industrially important fungi and yeasts. In this study, we successfully implemented the Squash-PCR technique in 12 industrially significant algae species. This approach offers a quick and reliable means of obtaining DNA templates directly from squashed algal cells, eliminating the need for time-consuming and labor-intensive cultivation and genomic DNA extraction steps. Our results demonstrate the effectiveness of Squash-PCR in detecting and characterizing target genes of interest in 12 different algae species. Overall, this study establishes the Squash-PCR method as a valuable tool for molecular studies in algae, enabling researchers to rapidly screen and manipulate genetic traits in diverse algal species.

Advances in genome-scale metabolic models of industrially important fungi.

Many fungal species have been used industrially for production of biofuels and bioproducts. Developing strains with better performance in biomanufacturing contexts requires a systematic understanding of cellular metabolism. Genome-scale metabolic models (GEMs) offer a comprehensive view of interconnected pathways and a mathematical framework for downstream analysis. Recently, GEMs have been developed or updated for several industrially important fungi. Some of them incorporate enzyme constraints, enabling improved predictions of cell states and proteome allocation. Here, we provide an overview of these newly developed GEMs and computational methods that facilitate construction of enzyme-constrained GEMs and utilize flux predictions from GEMs. Furthermore, we highlight the pivotal roles of these GEMs in iterative design-build-test-learn cycles, ultimately advancing the field of fungal biomanufacturing.

Identification of a specific exporter that enables high production of aconitic acid in Aspergillus pseudoterreus.

Aconitic acid is an unsaturated tricarboxylic acid that is attractive for its potential use in manufacturing biodegradable and biocompatible polymers, plasticizers, and surfactants. Previously Aspergillus pseudoterreus was engineered as a platform to produce aconitic acid by deleting the cadA (cis-aconitic acid decarboxylase) gene in the itaconic acid biosynthetic pathway. In this study, the aconitic acid transporter gene (aexA) was identified using comparative global discovery proteomics analysis between the wild-type and cadA deletion strains. The protein AexA belongs to the Major Facilitator Superfamily (MFS). Deletion of aexA almost abolished aconitic acid secretion, while its overexpression led to a significant increase in aconitic acid production. Transportation of aconitic acid across the plasma membrane is a key limiting step in its production. In vitro, proteoliposome transport assay further validated AexA's function and substrate specificity. This research provides new approaches to efficiently pinpoint and characterize exporters of fungal organic acids and accelerate metabolic engineering to improve secretion capability and lower the cost of bioproduction.

Rapid and robust squashed spore/colony PCR of industrially important fungi.

BACKGROUND: Fungi have been utilized for centuries in medical, agricultural, and industrial applications. Development of systems biology techniques has enabled the design and metabolic engineering of these fungi to produce novel fuels, chemicals, and enzymes from renewable feedstocks. Many genetic tools have been developed for manipulating the genome and creating mutants rapidly. However, screening and confirmation of transformants remain an inefficient step within the design, build, test, and learn cycle in many industrial fungi because extracting fungal genomic DNA is laborious, time-consuming, and involves toxic chemicals. RESULTS: In this study we developed a rapid and robust technique called "Squash-PCR" to break open the spores and release fungal genomic DNA as a template for PCR. The efficacy of Squash-PCR was investigated in eleven different filamentous fungal strains. Clean PCR products with high yields were achieved in all tested fungi. Spore age and type of DNA polymerase did not affect the efficiency of Squash-PCR. However, spore concentration was found to be the crucial factor for Squash-PCR in Aspergillus niger, with the dilution of starting material often resulting in higher PCR product yield. We then further evaluated the applicability of the squashing procedure for nine different yeast strains. We found that Squash-PCR can be used to improve the quality and yield of colony PCR in comparison to direct colony PCR in the tested yeast strains. CONCLUSION: The developed technique will enhance the efficiency of screening transformants and accelerate genetic engineering in filamentous fungi and yeast.

Engineering Rhodosporidium toruloides for production of 3-hydroxypropionic acid from lignocellulosic hydrolysate.

Microbial production of valuable bioproducts is a promising route towards green and sustainable manufacturing. The oleaginous yeast, Rhodosporidium toruloides, has emerged as an attractive host for the production of biofuels and bioproducts from lignocellulosic hydrolysates. 3-hydroxypropionic acid (3HP) is an attractive platform molecule that can be used to produce a wide range of commodity chemicals. This study focuses on establishing and optimizing the production of 3HP in R. toruloides. As R. toruloides naturally has a high metabolic flux towards malonyl-CoA, we exploited this pathway to produce 3HP. Upon finding the yeast capable of catabolizing 3HP, we then implemented functional genomics and metabolomic analysis to identify the catabolic pathways. Deletion of a putative malonate semialdehyde dehydrogenase gene encoding an oxidative 3HP pathway was found to significantly reduce 3HP degradation. We further explored monocarboxylate transporters to promote 3HP transport and identified a novel 3HP transporter in Aspergillus pseudoterreus by RNA-seq and proteomics. Combining these engineering efforts with media optimization in a fed-batch fermentation resulted in 45.4 g/L 3HP production. This represents one of the highest 3HP titers reported in yeast from lignocellulosic feedstocks. This work establishes R. toruloides as a host for 3HP production from lignocellulosic hydrolysate at high titers, and paves the way for further strain and process optimization towards enabling industrial production of 3HP in the future.

PeakDecoder enables machine learning-based metabolite annotation and accurate profiling in multidimensional mass spectrometry measurements.

Multidimensional measurements using state-of-the-art separations and mass spectrometry provide advantages in untargeted metabolomics analyses for studying biological and environmental bio-chemical processes. However, the lack of rapid analytical methods and robust algorithms for these heterogeneous data has limited its application. Here, we develop and evaluate a sensitive and high-throughput analytical and computational workflow to enable accurate metabolite profiling. Our workflow combines liquid chromatography, ion mobility spectrometry and data-independent acquisition mass spectrometry with PeakDecoder, a machine learning-based algorithm that learns to distinguish true co-elution and co-mobility from raw data and calculates metabolite identification error rates. We apply PeakDecoder for metabolite profiling of various engineered strains of Aspergillus pseudoterreus, Aspergillus niger, Pseudomonas putida and Rhodosporidium toruloides. Results, validated manually and against selected reaction monitoring and gas-chromatography platforms, show that 2683 features could be confidently annotated and quantified across 116 microbial sample runs using a library built from 64 standards.

Metabolic engineering to improve production of 3-hydroxypropionic acid from corn-stover hydrolysate in Aspergillus species.

BACKGROUND: Fuels and chemicals derived from non-fossil sources are needed to lessen human impacts on the environment while providing a healthy and growing economy. 3-hydroxypropionic acid (3-HP) is an important chemical building block that can be used for many products. Biosynthesis of 3-HP is possible; however, low production is typically observed in those natural systems. Biosynthetic pathways have been designed to produce 3-HP from a variety of feedstocks in different microorganisms. RESULTS: In this study, the 3-HP ß-alanine pathway consisting of aspartate decarboxylase, ß-alanine-pyruvate aminotransferase, and 3-hydroxypropionate dehydrogenase from selected microorganisms were codon optimized for Aspergillus species and placed under the control of constitutive promoters. The pathway was introduced into Aspergillus pseudoterreus and subsequently into Aspergillus niger, and 3-HP production was assessed in both hosts. A. niger produced higher initial 3-HP yields and fewer co-product contaminants and was selected as a suitable host for further engineering. Proteomic and metabolomic analysis of both Aspergillus species during 3-HP production identified genetic targets for improvement of flux toward 3-HP including pyruvate carboxylase, aspartate aminotransferase, malonate semialdehyde dehydrogenase, succinate semialdehyde dehydrogenase, oxaloacetate hydrolase, and a 3-HP transporter. Overexpression of pyruvate carboxylase improved yield in shake-flasks from 0.09 to 0.12 C-mol 3-HP C-mol-1 glucose in the base strain expressing 12 copies of the ß-alanine pathway. Deletion or overexpression of individual target genes in the pyruvate carboxylase overexpression strain improved yield to 0.22 C-mol 3-HP C-mol-1 glucose after deletion of the major malonate semialdehyde dehydrogenase. Further incorporation of additional ß-alanine pathway genes and optimization of culture conditions (sugars, temperature, nitrogen, phosphate, trace elements) for 3-HP production from deacetylated and mechanically refined corn stover hydrolysate improved yield to 0.48 C-mol 3-HP C-mol-1 sugars and resulted in a final titer of 36.0 g/L 3-HP. CONCLUSIONS: The results of this study establish A. niger as a host for 3-HP production from a lignocellulosic feedstock in acidic conditions and demonstrates that 3-HP titer and yield can be improved by a broad metabolic engineering strategy involving identification and modification of genes participated in the synthesis of 3-HP and its precursors, degradation of intermediates, and transport of 3-HP across the plasma membrane.

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.

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.

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.

Genetically Engineered Oleaginous Yeast Lipomyces starkeyi for Sesquiterpene α-Zingiberene Production.

Oleaginous yeast, such as Lipomyces starkeyi, are logical organisms for production of higher energy density molecules like lipids and terpenes. We demonstrate that transgenic L. starkeyi strains expressing an α-zingiberene synthase gene from lemon basil or Hall's panicgrass can produce up to 17 mg/L α-zingiberene in yeast extract peptone dextrose (YPD) medium containing 4% glucose. The transgenic strain was further examined in 8% glucose media with C/N ratios of 20 or 100, and YPD. YPD medium resulted in 59 mg/L α-zingiberene accumulation. Overexpression of selected genes from the mevalonate pathway achieved 145% improvement in α-zingiberene synthesis. Optimization of the growth medium for α-zingiberene production led to 15% higher titer than YPD medium. The final transgenic strain produced 700 mg/L α-zingiberene in fed-batch bioreactor culture. This study opens a new synthetic route to produce α-zingiberene or other terpenoids in L. starkeyi and establishes this yeast as a platform for jet fuel biosynthesis.

Deletion analysis of the itaconic acid biosynthesis gene cluster components in Aspergillus pseudoterreus ATCC32359.

The filamentous fungus Aspergillus terreus has been successfully used for industrial production of itaconic acid (IA) for many years. The IA biosynthesis pathway has recently been characterized at a molecular genetic level as an IA gene cluster by a clone-based transcriptomic approach. The cluster consists of four genes, including genes for cis-aconitic acid decarboxylase (cadA), a predicted transcription factor (tf), a mitochondrial organic acid transporter (mttA) and an MFS (major facilitator superfamily) type transporter (mfsA). In this research, we performed expressed sequence tag (EST) analysis and systematic gene deletions to further investigate the role of those genes during IA biosynthesis in A. pseudoterreus ATCC32359. EST analysis showed a similar expression pattern among those four genes that were distinct from neighboring genes and further confirmed that they belong to the same biosynthesis cluster. Systematic gene deletion analysis demonstrated that tf, cadA, mttA and mfsA genes in the cluster are essential for IA production; deletion of any of them will either completely abolish the IA production or dramatically decrease the amount of IA produced. The tf gene plays a regulatory role in this cluster. Deletion of tf led to decreased expression levels of cadA, mttA and mfsA. More importantly, a significant amount of aconitic acid was detected in the cadA deletion strain but not in the other deletion strains. Therefore, by deleting only one gene, the cadA, we established a novel microbial host for the production of aconitic acid and other value-added chemicals from sugars in lignocellulosic biomass.

Transcriptomic analysis of the oleaginous yeast Lipomyces starkeyi during lipid accumulation on enzymatically treated corn stover hydrolysate.

BACKGROUND: Efficient and economically viable production of biofuels from lignocellulosic biomass is dependent on mechanical and chemical pretreatment and enzymatic hydrolysis of plant material. These processing steps yield simple sugars as well as plant-derived and process-added organic acids, sugar-derived dehydration products, aldehydes, phenolics and other compounds that inhibit the growth of many microorganisms. Lipomyces starkeyi is an oleaginous yeast capable of robust growth on a variety of sugars and lipid accumulation on pretreated lignocellulosic substrates making it attractive as an industrial producer of biofuels. Here, we examined gene expression during batch growth and lipid accumulation in a 20-L bioreactor with either a blend of pure glucose and xylose or pretreated corn stover (PCS) that had been enzymatically hydrolyzed as the carbon sources. RESULTS: We monitored sugar and ammonium utilization as well as biomass accumulation and found that growth of L. starkeyi is inhibited with PCS hydrolysate as the carbon source. Both acetic acid and furfural are present at concentrations toxic to L. starkeyi in PCS hydrolysate. We quantified gene expression at seven time-points for each carbon source during batch growth and found that gene expression is similar at physiologically equivalent points. Analysis of promoter regions revealed that gene expression during the transition to lipid accumulation is regulated by carbon and nitrogen catabolite repression, regardless of carbon source and is associated with decreased expression of the translation machinery and suppression of the cell cycle. We identified 73 differentially expressed genes during growth phase in the bioreactor that may be involved in detoxification of corn stover hydrolysate. CONCLUSIONS: Growth of L. starkeyi is inhibited by compounds present in PCS hydrolysate. Here, we monitored key metabolites to establish physiologically equivalent comparisons during a batch bioreactor run comparing PCS hydrolysate and purified sugars. L. starkeyi's response to PCS hydrolysate is primarily at the beginning of the run during growth phase when inhibitory compounds are presumably at their highest concentration and inducing the general detoxification response by L. starkeyi. Differentially expressed genes identified herein during growth phase will aid in the improvement of industrial strains capable of robust growth on substrates containing various growth inhibitory compounds.

Deletion of the KU70 homologue facilitates gene targeting in Lipomyces starkeyi strain NRRL Y-11558.

The objective of this study was to disrupt the non-homologous end-joining (NHEJ) pathway gene (Lsku70Δ) and evaluate the effects of selected gene deletions related to glycogen synthesis (LsGSY1) and lipid degradation (LsMFE1, LsPEX10, and LsTGL4) on lipid production in the oleaginous yeast Lipomyces starkeyi. Disruption of the NHEJ pathway to reduce the rate of non-homologous recombination is a common approach used to overcome low-efficiency targeted deletion or insertion in various organisms. Here, the homologue of the LsKU70 gene was identified and disrupted in L. starkeyi NRRL Y-11558. The LsGSY1, LsMFE1, LsPEX10, LsTGL4, and LsURA3 genes were then replaced with a resistance marker in the Lsku70Δ strain and several site-specific insertions were assessed for targeted over-expression of selected genes. The targeted disruption efficiency of five selected genes (LsGSY1, LsMFE1, LsPEX10, LsTGL4, and LsURA3) was increased from 0 to 10% in the parent to 50-100% of transformants screened in the Lsku70Δ strain with 0.8-1.4 kb homologous flanking sequences, while the efficiency of site-specific gene insertion with the ß-glucuronidase reporter gene was 100% in the locus near the 3'-end coding (LsKU70) and non-coding (LsGSY1, LsMFE1, and LsPEX10) regions. Disruption of LsKU70 in isolation and in conjunction with LsGSY1, LsMFE1, LsPEX10, or LsTGL4 did not affect lipid production in L. starkeyi. Furthermore, ß-glucuronidase reporter gene activity was similar in strains containing site-specific targeted insertions. Therefore, over-expression of genes related to lipid synthesis at targeted loci can be further examined for improvement of total lipid production in L. starkeyi.

Regulation of Yeast-to-Hyphae Transition in Yarrowia lipolytica.

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. lipolyticamsn2 (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.

Draft Genome Sequence of Aspergillus oryzae ATCC 12892.

The draft genome sequence of Aspergillus oryzae ATCC 12892 is presented here. A. oryzae produces 3-nitropropionic acid, which has been investigated with regard to understanding the biosynthesis of nitroorganic compounds.

Multimodal microfluidic platform for controlled culture and analysis of unicellular organisms.

Modern live-cell imaging approaches permit real-time visualization of biological processes, yet limitations exist for unicellular organism isolation, culturing, and long-term imaging that preclude fully understanding how cells sense and respond to environmental perturbations and the link between single-cell variability and whole-population dynamics. Here, we present a microfluidic platform that provides fine control over the local environment with the capacity to replace media components at any experimental time point, and provides both perfused and compartmentalized cultivation conditions depending on the valve configuration. The functionality and flexibility of the platform were validated using both bacteria and yeast having different sizes, motility, and growth media. The demonstrated ability to track the growth and dynamics of both motile and non-motile prokaryotic and eukaryotic organisms emphasizes the versatility of the devices, which should enable studies in bioenergy and environmental research.

Omics Analyses of Trichoderma reesei CBS999.97 and QM6a Indicate the Relevance of Female Fertility to Carbohydrate-Active Enzyme and Transporter Levels.

The filamentous fungus Trichoderma reesei is found predominantly in the tropics but also in more temperate regions, such as Europe, and is widely known as a producer of large amounts of plant cell wall-degrading enzymes. We sequenced the genome of the sexually competent isolate CBS999.97, which is phenotypically different from the female sterile strain QM6a but can cross sexually with QM6a. Transcriptome data for growth on cellulose showed that entire carbohydrate-active enzyme (CAZyme) families are consistently differentially regulated between these strains. We evaluated backcrossed strains of both mating types, which acquired female fertility from CBS999.97 but maintained a mostly QM6a genetic background, and we could thereby distinguish between the effects of strain background and female fertility or mating type. We found clear regulatory differences associated with female fertility and female sterility, including regulation of CAZyme and transporter genes. Analysis of carbon source utilization, transcriptomes, and secondary metabolites in these strains revealed that only a few changes in gene regulation are consistently correlated with different mating types. Different strain backgrounds (QM6a versus CBS999.97) resulted in the most significant alterations in the transcriptomes and in carbon source utilization, with decreased growth of CBS999.97 on several amino acids (for example proline or alanine), which further correlated with the downregulation of genes involved in the respective pathways. In combination, our findings support a role of fertility-associated processes in physiology and gene regulation and are of high relevance for the use of sexual crossing in combining the characteristics of two compatible strains or quantitative trait locus (QTL) analysis.IMPORTANCETrichoderma reesei is a filamentous fungus with a high potential for secretion of plant cell wall-degrading enzymes. We sequenced the genome of the fully fertile field isolate CBS999.97 and analyzed its gene regulation characteristics in comparison with the commonly used laboratory wild-type strain QM6a, which is not female fertile. Additionally, we also evaluated fully fertile strains with genotypes very close to that of QM6a in order to distinguish between strain-specific and fertility-specific characteristics. We found that QM6a and CBS999.97 clearly differ in their growth patterns on different carbon sources, CAZyme gene regulation, and secondary metabolism. Importantly, we found altered regulation of 90 genes associated with female fertility, including CAZyme genes and transporter genes, but only minor mating type-dependent differences. Hence, when using sexual crossing in research and for strain improvement, it is important to consider female fertile and female sterile strains for comparison with QM6a and to achieve optimal performance.

Regulation of Nitrogen Metabolism by GATA Zinc Finger Transcription Factors in Yarrowia lipolytica.

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.