Engineering transcriptional regulation of pentose metabolism in Rhodosporidium toruloides for improved conversion of xylose to bioproducts.

Efficient conversion of pentose sugars remains a significant barrier to the replacement of petroleum-derived chemicals with plant biomass-derived bioproducts. While the oleaginous yeast Rhodosporidium toruloides (also known as Rhodotorula toruloides) has a relatively robust native metabolism of pentose sugars compared to other wild yeasts, faster assimilation of those sugars will be required for industrial utilization of pentoses. To increase the rate of pentose assimilation in R. toruloides, we leveraged previously reported high-throughput fitness data to identify potential regulators of pentose catabolism. Two genes were selected for further investigation, a putative transcription factor (RTO4_12978, Pnt1) and a homolog of a glucose transceptor involved in carbon catabolite repression (RTO4_11990). Overexpression of Pnt1 increased the specific growth rate approximately twofold early in cultures on xylose and increased the maximum specific growth by 18% while decreasing accumulation of arabitol and xylitol in fast-growing cultures. Improved growth dynamics on xylose translated to a 120% increase in the overall rate of xylose conversion to fatty alcohols in batch culture. Proteomic analysis confirmed that Pnt1 is a major regulator of pentose catabolism in R. toruloides. Deletion of RTO4_11990 increased the growth rate on xylose, but did not relieve carbon catabolite repression in the presence of glucose. Carbon catabolite repression signaling networks remain poorly characterized in R. toruloides and likely comprise a different set of proteins than those mainly characterized in ascomycete fungi.

Non-canonical D-xylose and L-arabinose metabolism via D-arabitol in the oleaginous yeast Rhodosporidium toruloides.

R. toruloides is an oleaginous yeast, with diverse metabolic capacities and high tolerance for inhibitory compounds abundant in plant biomass hydrolysates. While R. toruloides grows on several pentose sugars and alcohols, further engineering of the native pathway is required for efficient conversion of biomass-derived sugars to higher value bioproducts. A previous high-throughput study inferred that R. toruloides possesses a non-canonical L-arabinose and D-xylose metabolism proceeding through D-arabitol and D-ribulose. In this study, we present a combination of genetic and metabolite data that refine and extend that model. Chiral separations definitively illustrate that D-arabitol is the enantiomer that accumulates under pentose metabolism. Deletion of putative D-arabitol-2-dehydrogenase (RTO4_9990) results in > 75% conversion of D-xylose to D-arabitol, and is growth-complemented on pentoses by heterologous xylulose kinase expression. Deletion of putative D-ribulose kinase (RTO4_14368) arrests all growth on any pentose tested. Analysis of several pentose dehydrogenase mutants elucidates a complex pathway with multiple enzymes mediating multiple different reactions in differing combinations, from which we also inferred a putative L-ribulose utilization pathway. Our results suggest that we have identified enzymes responsible for the majority of pathway flux, with additional unknown enzymes providing accessory activity at multiple steps. Further biochemical characterization of the enzymes described here will enable a more complete and quantitative understanding of R. toruloides pentose metabolism. These findings add to a growing understanding of the diversity and complexity of microbial pentose metabolism.

Multi-Omics Driven Metabolic Network Reconstruction and Analysis of Lignocellulosic Carbon Utilization in Rhodosporidium toruloides.

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.

Genomewide and Enzymatic Analysis Reveals Efficient d-Galacturonic Acid Metabolism in the Basidiomycete Yeast Rhodosporidium toruloides.

Biorefining of renewable feedstocks is one of the most promising routes to replace fossil-based products. Since many common fermentation hosts, such as Saccharomyces cerevisiae, are naturally unable to convert many component plant cell wall polysaccharides, the identification of organisms with broad catabolism capabilities represents an opportunity to expand the range of substrates used in fermentation biorefinery approaches. The red basidiomycete yeast Rhodosporidium toruloides is a promising and robust host for lipid- and terpene-derived chemicals. Previous studies demonstrated assimilation of a range of substrates, from C5/C6 sugars to aromatic molecules similar to lignin monomers. In the current study, we analyzed the potential of R. toruloides to assimilate d-galacturonic acid, a major sugar in many pectin-rich agricultural waste streams, including sugar beet pulp and citrus peels. d-Galacturonic acid is not a preferred substrate for many fungi, but its metabolism was found to be on par with those of d-glucose and d-xylose in R. toruloides A genomewide analysis by combined transcriptome sequencing (RNA-seq) and RB-TDNA-seq revealed those genes with high relevance for fitness on d-galacturonic acid. While R. toruloides was found to utilize the nonphosphorylative catabolic pathway known from ascomycetes, the maximal velocities of several enzymes exceeded those previously reported. In addition, an efficient downstream glycerol catabolism and a novel transcription factor were found to be important for d-galacturonic acid utilization. These results set the basis for use of R. toruloides as a potential host for pectin-rich waste conversions and demonstrate its suitability as a model for metabolic studies with basidiomycetes.IMPORTANCE The switch from the traditional fossil-based industry to a green and sustainable bioeconomy demands the complete utilization of renewable feedstocks. Many currently used bioconversion hosts are unable to utilize major components of plant biomass, warranting the identification of microorganisms with broader catabolic capacity and characterization of their unique biochemical pathways. d-Galacturonic acid is a plant component of bioconversion interest and is the major backbone sugar of pectin, a plant cell wall polysaccharide abundant in soft and young plant tissues. The red basidiomycete and oleaginous yeast Rhodosporidium toruloides has been previously shown to utilize a range of sugars and aromatic molecules. Using state-of-the-art functional genomic methods and physiological and biochemical assays, we elucidated the molecular basis underlying the efficient metabolism of d-galacturonic acid. This study identified an efficient pathway for uronic acid conversion to guide future engineering efforts and represents the first detailed metabolic analysis of pectin metabolism in a basidiomycete fungus.

Synthetic Gene Network with Positive Feedback Loop Amplifies Cellulase Gene Expression in Neurospora crassa.

Second-generation or lignocellulosic biofuels are a tangible source of renewable energy, which is critical to combat climate change by reducing the carbon footprint. Filamentous fungi secrete cellulose-degrading enzymes called cellulases, which are used for production of lignocellulosic biofuels. However, inefficient production of cellulases is a major obstacle for industrial-scale production of second-generation biofuels. We used computational simulations to design and implement synthetic positive feedback loops to increase gene expression of a key transcription factor, CLR-2, that activates a large number of cellulases in a filamentous fungus, Neurospora crassa. Overexpression of CLR-2 reveals previously unappreciated roles of CLR-2 in lignocellulosic gene network, which enabled simultaneous induction of approximately 50% of 78 lignocellulosic degradation-related genes in our engineered Neurospora strains. This engineering results in dramatically increased cellulase activity due to cooperative orchestration of multiple enzymes involved in the cellulose degradation pathway. Our work provides a proof of principle in utilizing mathematical modeling and synthetic biology to improve the efficiency of cellulase synthesis for second-generation biofuel production.

Functional genomics of lipid metabolism in the oleaginous yeast Rhodosporidium toruloides.

The basidiomycete yeast Rhodosporidium toruloides (also known as Rhodotorula toruloides) accumulates high concentrations of lipids and carotenoids from diverse carbon sources. It has great potential as a model for the cellular biology of lipid droplets and for sustainable chemical production. We developed a method for high-throughput genetics (RB-TDNAseq), using sequence-barcoded Agrobacterium tumefaciens T-DNA insertions. We identified 1,337 putative essential genes with low T-DNA insertion rates. We functionally profiled genes required for fatty acid catabolism and lipid accumulation, validating results with 35 targeted deletion strains. We identified a high-confidence set of 150 genes affecting lipid accumulation, including genes with predicted function in signaling cascades, gene expression, protein modification and vesicular trafficking, autophagy, amino acid synthesis and tRNA modification, and genes of unknown function. These results greatly advance our understanding of lipid metabolism in this oleaginous species and demonstrate a general approach for barcoded mutagenesis that should enable functional genomics in diverse fungi.

Network of nutrient-sensing pathways and a conserved kinase cascade integrate osmolarity and carbon sensing in Neurospora crassa.

Identifying nutrients available in the environment and utilizing them in the most efficient manner is a challenge common to all organisms. The model filamentous fungus Neurospora crassa is capable of utilizing a variety of carbohydrates, from simple sugars to the complex carbohydrates found in plant cell walls. The zinc binuclear cluster transcription factor CLR-1 is necessary for utilization of cellulose, a major, recalcitrant component of the plant cell wall; however, expression of clr-1 in the absence of an inducer is not sufficient to induce cellulase gene expression. We performed a screen for unidentified actors in the cellulose-response pathway and identified a gene encoding a hypothetical protein (clr-3) that is required for repression of CLR-1 activity in the absence of an inducer. Using clr-3 mutants, we implicated the hyperosmotic-response pathway in the tunable regulation of glycosyl hydrolase production in response to changes in osmolarity. The role of the hyperosmotic-response pathway in nutrient sensing may indicate that cells use osmolarity as a proxy for the presence of free sugar in their environment. These signaling pathways form a nutrient-sensing network that allows Ncrassa cells to tightly regulate gene expression in response to environmental conditions.

Network reconstruction and systems analysis of plant cell wall deconstruction by Neurospora crassa.

BACKGROUND: Plant biomass degradation by fungal-derived enzymes is rapidly expanding in economic importance as a clean and efficient source for biofuels. The ability to rationally engineer filamentous fungi would facilitate biotechnological applications for degradation of plant cell wall polysaccharides. However, incomplete knowledge of biomolecular networks responsible for plant cell wall deconstruction impedes experimental efforts in this direction. RESULTS: To expand this knowledge base, a detailed network of reactions important for deconstruction of plant cell wall polysaccharides into simple sugars was constructed for the filamentous fungus Neurospora crassa. To reconstruct this network, information was integrated from five heterogeneous data types: functional genomics, transcriptomics, proteomics, genetics, and biochemical characterizations. The combined information was encapsulated into a feature matrix and the evidence weighted to assign annotation confidence scores for each gene within the network. Comparative analyses of RNA-seq and ChIP-seq data shed light on the regulation of the plant cell wall degradation network, leading to a novel hypothesis for degradation of the hemicellulose mannan. The transcription factor CLR-2 was subsequently experimentally shown to play a key role in the mannan degradation pathway of N. crassa. CONCLUSIONS: Here we built a network that serves as a scaffold for integration of diverse experimental datasets. This approach led to the elucidation of regulatory design principles for plant cell wall deconstruction by filamentous fungi and a novel function for the transcription factor CLR-2. This expanding network will aid in efforts to rationally engineer industrially relevant hyper-production strains.

Direct target network of the Neurospora crassa plant cell wall deconstruction regulators CLR-1, CLR-2, and XLR-1.

UNLABELLED: Fungal deconstruction of the plant cell requires a complex orchestration of a wide array of intracellular and extracellular enzymes. In Neurospora crassa, CLR-1, CLR-2, and XLR-1 have been identified as key transcription factors regulating plant cell wall degradation in response to soluble sugars. The XLR-1 regulon was defined using a constitutively active mutant allele, resulting in hemicellulase gene expression and secretion under noninducing conditions. To define genes directly regulated by CLR-1, CLR-2, and XLR-1, we performed chromatin immunoprecipitation and next-generation sequencing (ChIPseq) on epitope-tagged constructs of these three transcription factors. When N. crassa is exposed to plant cell wall material, CLR-1, CLR-2, and XLR-1 individually bind to the promoters of the most strongly induced genes in their respective regulons. These include promoters of genes encoding cellulases for CLR-1 and CLR-2 (CLR-1/CLR-2) and promoters of genes encoding hemicellulases for XLR-1. CLR-1 bound to its regulon under noninducing conditions; however, this binding alone did not translate into gene expression and enzyme secretion. Motif analysis of the bound genes revealed conserved DNA binding motifs, with the CLR-2 motif matching that of its closest paralog in Saccharomyces cerevisiae, Gal4p. Coimmunoprecipitation studies showed that CLR-1 and CLR-2 act in a homocomplex but not as a CLR-1/CLR-2 heterocomplex. IMPORTANCE: Understanding fungal regulation of complex plant cell wall deconstruction pathways in response to multiple environmental signals via interconnected transcriptional circuits provides insight into fungus/plant interactions and eukaryotic nutrient sensing. Coordinated optimization of these regulatory networks is likely required for optimal microbial enzyme production.

The proteome and phosphoproteome of Neurospora crassa in response to cellulose, sucrose and carbon starvation.

Improving cellulolytic enzyme production by plant biomass degrading fungi holds great potential in reducing costs associated with production of next-generation biofuels generated from lignocellulose. How fungi sense cellulosic materials and respond by secreting enzymes has mainly been examined by assessing function of transcriptional regulators and via transcriptional profiling. Here, we obtained global proteomic and phosphoproteomic profiles of the plant biomass degrading filamentous fungus Neurospora crassa grown on different carbon sources, i.e. sucrose, no carbon, and cellulose, by performing isobaric tags for relative and absolute quantification (iTRAQ)-based LC-MS/MS analyses. A comparison between proteomes and transcriptomes under identical carbon conditions suggests that extensive post-transcriptional regulation occurs in N. crassa in response to exposure to cellulosic material. Several hundred amino acid residues with differential phosphorylation levels on crystalline cellulose (Avicel) or carbon-free medium vs sucrose medium were identified, including phosphorylation sites in a major transcriptional activator for cellulase genes, CLR1, as well as a cellobionic acid transporter, CBT1. Mutation of phosphorylation sites on CLR1 did not have a major effect on transactivation of cellulase production, while mutation of phosphorylation sites in CBT1 increased its transporting capacity. Our data provides rich information at both the protein and phosphorylation levels of the early cellular responses to carbon starvation and cellulosic induction and aids in a greater understanding of the underlying post-transcriptional regulatory mechanisms in filamentous fungi.

Analysis of a conserved cellulase transcriptional regulator reveals inducer-independent production of cellulolytic enzymes in Neurospora crassa.

Cellulose is recalcitrant to deconstruction to glucose for use in fermentation strategies for biofuels and chemicals derived from lignocellulose. In Neurospora crassa, the transcriptional regulator, CLR-2, is required for cellulolytic gene expression and cellulose deconstruction. To assess conservation and divergence of cellulase gene regulation between fungi from different ecological niches, we compared clr-2 function with its ortholog (clrB) in the distantly related species, Aspergillus nidulans. Transcriptional profiles induced by exposure to crystalline cellulose were similar in both species. Approximately 50% of the cellulose-responsive genes showed strict dependence on functional clr-2/clrB, with a subset of 28 genes encoding plant biomass degrading enzymes that were conserved between N. crassa and A. nidulans. Importantly, misexpression of clr-2 under noninducing conditions was sufficient to drive cellulase gene expression, secretion, and activity in N. crassa, to a level comparable to wild type exposed to Avicel. However, misexpression of clrB in A. nidulans was not sufficient to drive cellulase gene expression under noninducing conditions, although an increase in cellulase activity was observed under crystalline cellulose conditions. Manipulation of clr-2 orthologs among filamentous fungi may enable regulated cellulosic enzyme production in a wide array of culture conditions and host strains, potentially reducing costs associated with enzyme production for plant cell wall deconstruction. However, this functionality may require additional engineering in some species.

Induction of lignocellulose-degrading enzymes in Neurospora crassa by cellodextrins.

Neurospora crassa colonizes burnt grasslands in the wild and metabolizes both cellulose and hemicellulose from plant cell walls. When switched from a favored carbon source such as sucrose to cellulose, N. crassa dramatically upregulates expression and secretion of a wide variety of genes encoding lignocellulolytic enzymes. However, the means by which N. crassa and other filamentous fungi sense the presence of cellulose in the environment remains unclear. Here, we show that an N. crassa mutant carrying deletions of two genes encoding extracellular ß-glucosidase enzymes and one intracellular ß-glucosidase lacks ß-glucosidase activity, but efficiently induces cellulase gene expression in the presence of cellobiose, cellotriose, or cellotetraose as a sole carbon source. These data indicate that cellobiose, or a modified version of cellobiose, functions as an inducer of lignocellulolytic gene expression in N. crassa. Furthermore, the inclusion of a deletion of the catabolite repressor gene, cre-1, in the triple ß-glucosidase mutant resulted in a strain that produces higher concentrations of secreted active cellulases on cellobiose. Thus, the ability to induce cellulase gene expression using a common and soluble carbon source simplifies enzyme production and characterization, which could be applied to other cellulolytic filamentous fungi.

Conserved and essential transcription factors for cellulase gene expression in ascomycete fungi.

Rational engineering of filamentous fungi for improved cellulase production is hampered by our incomplete knowledge of transcriptional regulatory networks. We therefore used the model filamentous fungus Neurospora crassa to search for uncharacterized transcription factors associated with cellulose deconstruction. A screen of a N. crassa transcription factor deletion collection identified two uncharacterized zinc binuclear cluster transcription factors (clr-1 and clr-2) that were required for growth and enzymatic activity on cellulose, but were not required for growth or hemicellulase activity on xylan. Transcriptional profiling with next-generation sequencing methods refined our understanding of the N. crassa transcriptional response to cellulose and demonstrated that clr-1 and clr-2 were required for the bulk of that response, including induction of all major cellulase and some major hemicellulase genes. Functional CLR-1 was necessary for expression of clr-2 and efficient cellobiose utilization. Phylogenetic analyses showed that CLR-1 and CLR-2 are conserved in the genomes of most filamentous ascomycete fungi capable of degrading cellulose. In Aspergillus nidulans, a strain carrying a deletion of the clr-2 homolog (clrB) failed to induce cellulase gene expression and lacked cellulolytic activity on Avicel. Further manipulation of this control system in industrial production strains may significantly improve yields of cellulases for cellulosic biofuel production.

Association (micellization) and partitioning of aglycon triterpenoids.

Micellization and solution properties of the aglycon triterpenoids asiatic acid (AA) and madecassic acid (MA) were examined experimentally and in computational simulations. AA and MA belong to the large class of bioactive aglycon triterpenoids, for which limited physicochemical data are available. In this study, solubility, partition coefficient, critical micelle concentrations (CMC), and surface tensions of AA and MA were measured. Reverse phase HPLC data, supported by dye probe experiments and drop shape analysis, showed the CMC in phosphate buffered saline (PBS) to be 15+/-2 microM, and 86+/-9 microM for AA and MA, respectively. The surface tensions of AA and MA in PBS were 64.1 and 64.4 mN/m, respectively. Matrix-assisted laser desorption/ionization (MALDI) mass spectrometry indicated the aggregation numbers of AA and MA to be 5 to 7. Molecular dynamics simulations confirmed that molecular association could occur between 5 and 7 molecules in solution. The IC(50) of AA and MA on human small cell carcinoma and human glioblastoma cell lines was 25+/-5 microM and 66+/-13 microM, respectively. The IC(50) is within the range of calculated CMC of AA and MA in bioassay media, suggesting that the micellar aggregates may lead to their cytotoxicity.