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.

Circulating SARS-CoV-2+ megakaryocytes are associated with severe viral infection in COVID-19.

Several independent lines of evidence suggest that megakaryocytes are dysfunctional in severe COVID-19. Herein, we characterized peripheral circulating megakaryocytes in a large cohort of inpatients with COVID-19 and correlated the subpopulation frequencies with clinical outcomes. Using peripheral blood, we show that megakaryocytes are increased in the systemic circulation in COVID-19, and we identify and validate S100A8/A9 as a defining marker of megakaryocyte dysfunction. We further reveal a subpopulation of S100A8/A9+ megakaryocytes that contain severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) protein and RNA. Using flow cytometry of peripheral blood and in vitro studies on SARS-CoV-2-infected primary human megakaryocytes, we demonstrate that megakaryocytes can transfer viral antigens to emerging platelets. Mechanistically, we show that SARS-CoV-2-containing megakaryocytes are nuclear factor κB (NF-κB)-activated, via p65 and p52; express the NF-κB-mediated cytokines interleukin-6 (IL-6) and IL-1ß; and display high surface expression of Toll-like receptor 2 (TLR2) and TLR4, canonical drivers of NF-κB. In a cohort of 218 inpatients with COVID-19, we correlate frequencies of megakaryocyte subpopulations with clinical outcomes and show that SARS-CoV-2-containing megakaryocytes are a strong risk factor for mortality and multiorgan injury, including respiratory failure, mechanical ventilation, acute kidney injury, thrombotic events, and intensive care unit admission. Furthermore, we show that SARS-CoV-2+ megakaryocytes are present in lung and brain autopsy tissues from deceased donors who had COVID-19. To our knowledge, this study offers the first evidence implicating SARS-CoV-2+ peripheral megakaryocytes in severe disease and suggests that circulating megakaryocytes warrant investigation in inflammatory disorders beyond COVID-19.

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.

A Histoplasma capsulatum Lipid Metabolic Map Identifies Antifungal Targets.

Lipids play a fundamental role in fungal cell biology, being essential cell membrane components and major targets of antifungal drugs. A deeper knowledge of lipid metabolism is key for developing new drugs and a better understanding of fungal pathogenesis. Here, we built a comprehensive map of the Histoplasma capsulatum lipid metabolic pathway by incorporating proteomic and lipidomic analyses. We performed genetic complementation and overexpression of H. capsulatum genes in Saccharomyces cerevisiae to validate reactions identified in the map and to determine enzymes responsible for catalyzing orphan reactions. The map led to the identification of both the fatty acid desaturation and the sphingolipid biosynthesis pathways as targets for drug development. We found that the sphingolipid biosynthesis inhibitor myriocin, the fatty acid desaturase inhibitor thiocarlide, and the fatty acid analog 10-thiastearic acid inhibit H. capsulatum growth in nanomolar to low-micromolar concentrations. These compounds also reduced the intracellular infection in an alveolar macrophage cell line. Overall, this lipid metabolic map revealed pathways that can be targeted for drug development. IMPORTANCE It is estimated that 150 people die per hour due to the insufficient therapeutic treatments to combat fungal infections. A major hurdle to developing antifungal therapies is the scarce knowledge on the fungal metabolic pathways and mechanisms of virulence. In this context, fungal lipid metabolism is an excellent candidate for developing drugs due to its essential roles in cellular scaffolds, energy storage, and signaling transductors. Here, we provide a detailed map of Histoplasma capsulatum lipid metabolism. The map revealed points of this fungus lipid metabolism that can be targeted for developing antifungal drugs.

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.

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

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

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.

Metaphenomic Responses of a Native Prairie Soil Microbiome to Moisture Perturbations.

Climate change is causing shifts in precipitation patterns in the central grasslands of the United States, with largely unknown consequences on the collective physiological responses of the soil microbial community, i.e., the metaphenome. Here, we used an untargeted omics approach to determine the soil microbial community's metaphenomic response to soil moisture and to define specific metabolic signatures of the response. Specifically, we aimed to develop the technical approaches and metabolic mapping framework necessary for future systematic ecological studies. We collected soil from three locations at the Konza Long-Term Ecological Research (LTER) field station in Kansas, and the soils were incubated for 15 days under dry or wet conditions and compared to field-moist controls. The microbiome response to wetting or drying was determined by 16S rRNA amplicon sequencing, metatranscriptomics, and metabolomics, and the resulting shifts in taxa, gene expression, and metabolites were assessed. Soil drying resulted in significant shifts in both the composition and function of the soil microbiome. In contrast, there were few changes following wetting. The combined metabolic and metatranscriptomic data were used to generate reaction networks to determine the metaphenomic response to soil moisture transitions. Site location was a strong determinant of the response of the soil microbiome to moisture perturbations. However, some specific metabolic pathways changed consistently across sites, including an increase in pathways and metabolites for production of sugars and other osmolytes as a response to drying. Using this approach, we demonstrate that despite the high complexity of the soil habitat, it is possible to generate insight into the effect of environmental change on the soil microbiome and its physiology and functions, thus laying the groundwork for future, targeted studies.IMPORTANCE Climate change is predicted to result in increased drought extent and intensity in the highly productive, former tallgrass prairie region of the continental United States. These soils store large reserves of carbon. The decrease in soil moisture due to drought has largely unknown consequences on soil carbon cycling and other key biogeochemical cycles carried out by soil microbiomes. In this study, we found that soil drying had a significant impact on the structure and function of soil microbial communities, including shifts in expression of specific metabolic pathways, such as those leading toward production of osmoprotectant compounds. This study demonstrates the application of an untargeted multi-omics approach to decipher details of the soil microbial community's metaphenotypic response to environmental perturbations and should be applicable to studies of other complex microbial systems as well.

Circadian Proteomic Analysis Uncovers Mechanisms of Post-Transcriptional Regulation in Metabolic Pathways.

Transcriptional and translational feedback loops in fungi and animals drive circadian rhythms in transcript levels that provide output from the clock, but post-transcriptional mechanisms also contribute. To determine the extent and underlying source of this regulation, we applied newly developed analytical tools to a long-duration, deeply sampled, circadian proteomics time course comprising half of the proteome. We found a quarter of expressed proteins are clock regulated, but >40% of these do not arise from clock-regulated transcripts, and our analysis predicts that these protein rhythms arise from oscillations in translational rates. Our data highlighted the impact of the clock on metabolic regulation, with central carbon metabolism reflecting both transcriptional and post-transcriptional control and opposing metabolic pathways showing peak activities at different times of day. The transcription factor CSP-1 plays a role in this metabolic regulation, contributing to the rhythmicity and phase of clock-regulated proteins.

Prediction of Metabolite Concentrations, Rate Constants and Post-Translational Regulation Using Maximum Entropy-Based Simulations with Application to Central Metabolism of Neurospora crassa.

We report the application of a recently proposed approach for modeling biological systems using a maximum entropy production rate principle in lieu of having in vivo rate constants. The method is applied in four steps: (1) a new ordinary differential equation (ODE) based optimization approach based on Marcelin's 1910 mass action equation is used to obtain the maximum entropy distribution; (2) the predicted metabolite concentrations are compared to those generally expected from experiments using a loss function from which post-translational regulation of enzymes is inferred; (3) the system is re-optimized with the inferred regulation from which rate constants are determined from the metabolite concentrations and reaction fluxes; and finally (4) a full ODE-based, mass action simulation with rate parameters and allosteric regulation is obtained. From the last step, the power characteristics and resistance of each reaction can be determined. The method is applied to the central metabolism of Neurospora crassa and the flow of material through the three competing pathways of upper glycolysis, the non-oxidative pentose phosphate pathway, and the oxidative pentose phosphate pathway are evaluated as a function of the NADP/NADPH ratio. It is predicted that regulation of phosphofructokinase (PFK) and flow through the pentose phosphate pathway are essential for preventing an extreme level of fructose 1,6-bisphophate accumulation. Such an extreme level of fructose 1,6-bisphophate would otherwise result in a glassy cytoplasm with limited diffusion, dramatically decreasing the entropy and energy production rate and, consequently, biological competitiveness.

Moleculo Long-Read Sequencing Facilitates Assembly and Genomic Binning from Complex Soil Metagenomes.

Soil metagenomics has been touted as the "grand challenge" for metagenomics, as the high microbial diversity and spatial heterogeneity of soils make them unamenable to current assembly platforms. Here, we aimed to improve soil metagenomic sequence assembly by applying the Moleculo synthetic long-read sequencing technology. In total, we obtained 267 Gbp of raw sequence data from a native prairie soil; these data included 109.7 Gbp of short-read data (~100 bp) from the Joint Genome Institute (JGI), an additional 87.7 Gbp of rapid-mode read data (~250 bp), plus 69.6 Gbp (>1.5 kbp) from Moleculo sequencing. The Moleculo data alone yielded over 5,600 reads of >10 kbp in length, and over 95% of the unassembled reads mapped to contigs of >1.5 kbp. Hybrid assembly of all data resulted in more than 10,000 contigs over 10 kbp in length. We mapped three replicate metatranscriptomes derived from the same parent soil to the Moleculo subassembly and found that 95% of the predicted genes, based on their assignments to Enzyme Commission (EC) numbers, were expressed. The Moleculo subassembly also enabled binning of >100 microbial genome bins. We obtained via direct binning the first complete genome, that of "Candidatus Pseudomonas sp. strain JKJ-1" from a native soil metagenome. By mapping metatranscriptome sequence reads back to the bins, we found that several bins corresponding to low-relative-abundance Acidobacteria were highly transcriptionally active, whereas bins corresponding to high-relative-abundance Verrucomicrobia were not. These results demonstrate that Moleculo sequencing provides a significant advance for resolving complex soil microbial communities. IMPORTANCE Soil microorganisms carry out key processes for life on our planet, including cycling of carbon and other nutrients and supporting growth of plants. However, there is poor molecular-level understanding of their functional roles in ecosystem stability and responses to environmental perturbations. This knowledge gap is largely due to the difficulty in culturing the majority of soil microbes. Thus, use of culture-independent approaches, such as metagenomics, promises the direct assessment of the functional potential of soil microbiomes. Soil is, however, a challenge for metagenomic assembly due to its high microbial diversity and variable evenness, resulting in low coverage and uneven sampling of microbial genomes. Despite increasingly large soil metagenome data volumes (>200 Gbp), the majority of the data do not assemble. Here, we used the cutting-edge approach of synthetic long-read sequencing technology (Moleculo) to assemble soil metagenome sequence data into long contigs and used the assemblies for binning of genomes. Author Video: An author video summary of this article is available.

Summary of the DREAM8 Parameter Estimation Challenge: Toward Parameter Identification for Whole-Cell Models.

Whole-cell models that explicitly represent all cellular components at the molecular level have the potential to predict phenotype from genotype. However, even for simple bacteria, whole-cell models will contain thousands of parameters, many of which are poorly characterized or unknown. New algorithms are needed to estimate these parameters and enable researchers to build increasingly comprehensive models. We organized the Dialogue for Reverse Engineering Assessments and Methods (DREAM) 8 Whole-Cell Parameter Estimation Challenge to develop new parameter estimation algorithms for whole-cell models. We asked participants to identify a subset of parameters of a whole-cell model given the model's structure and in silico "experimental" data. Here we describe the challenge, the best performing methods, and new insights into the identifiability of whole-cell models. We also describe several valuable lessons we learned toward improving future challenges. Going forward, we believe that collaborative efforts supported by inexpensive cloud computing have the potential to solve whole-cell model parameter estimation.

Reconstruction and validation of a genome-scale metabolic model for the filamentous fungus Neurospora crassa using FARM.

The filamentous fungus Neurospora crassa played a central role in the development of twentieth-century genetics, biochemistry and molecular biology, and continues to serve as a model organism for eukaryotic biology. Here, we have reconstructed a genome-scale model of its metabolism. This model consists of 836 metabolic genes, 257 pathways, 6 cellular compartments, and is supported by extensive manual curation of 491 literature citations. To aid our reconstruction, we developed three optimization-based algorithms, which together comprise Fast Automated Reconstruction of Metabolism (FARM). These algorithms are: LInear MEtabolite Dilution Flux Balance Analysis (limed-FBA), which predicts flux while linearly accounting for metabolite dilution; One-step functional Pruning (OnePrune), which removes blocked reactions with a single compact linear program; and Consistent Reproduction Of growth/no-growth Phenotype (CROP), which reconciles differences between in silico and experimental gene essentiality faster than previous approaches. Against an independent test set of more than 300 essential/non-essential genes that were not used to train the model, the model displays 93% sensitivity and specificity. We also used the model to simulate the biochemical genetics experiments originally performed on Neurospora by comprehensively predicting nutrient rescue of essential genes and synthetic lethal interactions, and we provide detailed pathway-based mechanistic explanations of our predictions. Our model provides a reliable computational framework for the integration and interpretation of ongoing experimental efforts in Neurospora, and we anticipate that our methods will substantially reduce the manual effort required to develop high-quality genome-scale metabolic models for other organisms.

TB database 2010: overview and update.

The Tuberculosis Database (TBDB) is an online database providing integrated access to genome sequence, expression data and literature curation for TB. TBDB currently houses genome assemblies for numerous strains of Mycobacterium tuberculosis (MTB) as well assemblies for over 20 strains related to MTB and useful for comparative analysis. TBDB stores pre- and post-publication gene-expression data from M. tuberculosis and its close relatives, including over 3000 MTB microarrays, 95 RT-PCR datasets, 2700 microarrays for human and mouse TB related experiments, and 260 arrays for Streptomyces coelicolor. To enable wide use of these data, TBDB provides a suite of tools for searching, browsing, analyzing, and downloading the data. We provide here an overview of TBDB focusing on recent data releases and enhancements. In particular, we describe the recent release of a Global Genetic Diversity dataset for TB, support for short-read re-sequencing data, new tools for exploring gene expression data in the context of gene regulation, and the integration of a metabolic network reconstruction and BioCyc with TBDB. By integrating a wide range of genomic data with tools for their use, TBDB is a unique platform for both basic science research in TB, as well as research into the discovery and development of TB drugs, vaccines and biomarkers.