RESUMO
SUMMARY: We introduce Eliater, a Python package for estimating the effect of perturbation of an upstream molecule on a downstream molecule in a biomolecular network. The estimation takes as input a biomolecular network, observational biomolecular data, and a perturbation of interest, and outputs an estimated quantitative effect of the perturbation. We showcase the functionalities of Eliater in a case study of E. coli transcriptional regulatory network. AVAILABILITY AND IMPLEMENTATION: The code, the documentation, and several case studies are available open source at https://github.com/y0-causal-inference/eliater.
RESUMO
Causal query estimation in biomolecular networks commonly selects a 'valid adjustment set', i.e. a subset of network variables that eliminates the bias of the estimator. A same query may have multiple valid adjustment sets, each with a different variance. When networks are partially observed, current methods use graph-based criteria to find an adjustment set that minimizes asymptotic variance. Unfortunately, many models that share the same graph topology, and therefore same functional dependencies, may differ in the processes that generate the observational data. In these cases, the topology-based criteria fail to distinguish the variances of the adjustment sets. This deficiency can lead to sub-optimal adjustment sets, and to miss-characterization of the effect of the intervention. We propose an approach for deriving 'optimal adjustment sets' that takes into account the nature of the data, bias and finite-sample variance of the estimator, and cost. It empirically learns the data generating processes from historical experimental data, and characterizes the properties of the estimators by simulation. We demonstrate the utility of the proposed approach in four biomolecular Case studies with different topologies and different data generation processes. The implementation and reproducible Case studies are at https://github.com/srtaheri/OptimalAdjustmentSet.
Assuntos
Biologia Computacional , Simulação por ComputadorRESUMO
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.
Assuntos
COVID-19 , Humanos , SARS-CoV-2 , Megacariócitos/metabolismo , NF-kappa B/metabolismo , Pulmão/metabolismoRESUMO
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.
RESUMO
Self-driving labs (SDLs) combine fully automated experiments with artificial intelligence (AI) that decides the next set of experiments. Taken to their ultimate expression, SDLs could usher a new paradigm of scientific research, where the world is probed, interpreted, and explained by machines for human benefit. While there are functioning SDLs in the fields of chemistry and materials science, we contend that synthetic biology provides a unique opportunity since the genome provides a single target for affecting the incredibly wide repertoire of biological cell behavior. However, the level of investment required for the creation of biological SDLs is only warranted if directed toward solving difficult and enabling biological questions. Here, we discuss challenges and opportunities in creating SDLs for synthetic biology.
Assuntos
Inteligência Artificial , Biologia Sintética , HumanosRESUMO
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.
RESUMO
MOTIVATION: Estimating causal queries, such as changes in protein abundance in response to a perturbation, is a fundamental task in the analysis of biomolecular pathways. The estimation requires experimental measurements on the pathway components. However, in practice many pathway components are left unobserved (latent) because they are either unknown, or difficult to measure. Latent variable models (LVMs) are well-suited for such estimation. Unfortunately, LVM-based estimation of causal queries can be inaccurate when parameters of the latent variables are not uniquely identified, or when the number of latent variables is misspecified. This has limited the use of LVMs for causal inference in biomolecular pathways. RESULTS: In this article, we propose a general and practical approach for LVM-based estimation of causal queries. We prove that, despite the challenges above, LVM-based estimators of causal queries are accurate if the queries are identifiable according to Pearl's do-calculus and describe an algorithm for its estimation. We illustrate the breadth and the practical utility of this approach for estimating causal queries in four synthetic and two experimental case studies, where structures of biomolecular pathways challenge the existing methods for causal query estimation. AVAILABILITY AND IMPLEMENTATION: The code and the data documenting all the case studies are available at https://github.com/srtaheri/LVMwithDoCalculus. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.
Assuntos
Algoritmos , Cálculos , Humanos , Modelos Teóricos , ProteínasRESUMO
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.
Assuntos
Histoplasma/genética , Histoplasma/metabolismo , Metabolismo dos Lipídeos , Ácidos Graxos/biossíntese , Proteínas Fúngicas/genética , Proteínas Fúngicas/metabolismo , Histoplasma/crescimento & desenvolvimento , Histoplasmose/microbiologia , Humanos , Lipidômica , Proteômica , Esfingolipídeos/biossínteseRESUMO
We need to effectively combine the knowledge from surging literature with complex datasets to propose mechanistic models of SARS-CoV-2 infection, improving data interpretation and predicting key targets of intervention. Here, we describe a large-scale community effort to build an open access, interoperable and computable repository of COVID-19 molecular mechanisms. The COVID-19 Disease Map (C19DMap) is a graphical, interactive representation of disease-relevant molecular mechanisms linking many knowledge sources. Notably, it is a computational resource for graph-based analyses and disease modelling. To this end, we established a framework of tools, platforms and guidelines necessary for a multifaceted community of biocurators, domain experts, bioinformaticians and computational biologists. The diagrams of the C19DMap, curated from the literature, are integrated with relevant interaction and text mining databases. We demonstrate the application of network analysis and modelling approaches by concrete examples to highlight new testable hypotheses. This framework helps to find signatures of SARS-CoV-2 predisposition, treatment response or prioritisation of drug candidates. Such an approach may help deal with new waves of COVID-19 or similar pandemics in the long-term perspective.
Assuntos
COVID-19/imunologia , Biologia Computacional/métodos , Bases de Dados Factuais , SARS-CoV-2/imunologia , Software , Antivirais/uso terapêutico , COVID-19/genética , COVID-19/virologia , Gráficos por Computador , Citocinas/genética , Citocinas/imunologia , Mineração de Dados/estatística & dados numéricos , Regulação da Expressão Gênica , Interações entre Hospedeiro e Microrganismos/genética , Interações entre Hospedeiro e Microrganismos/imunologia , Humanos , Imunidade Celular/efeitos dos fármacos , Imunidade Humoral/efeitos dos fármacos , Imunidade Inata/efeitos dos fármacos , Linfócitos/efeitos dos fármacos , Linfócitos/imunologia , Linfócitos/virologia , Redes e Vias Metabólicas/genética , Redes e Vias Metabólicas/imunologia , Células Mieloides/efeitos dos fármacos , Células Mieloides/imunologia , Células Mieloides/virologia , Mapeamento de Interação de Proteínas , SARS-CoV-2/efeitos dos fármacos , SARS-CoV-2/genética , SARS-CoV-2/patogenicidade , Transdução de Sinais , Fatores de Transcrição/genética , Fatores de Transcrição/imunologia , Proteínas Virais/genética , Proteínas Virais/imunologia , Tratamento Farmacológico da COVID-19RESUMO
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.
Assuntos
Yarrowia/metabolismo , Acetato-CoA Ligase/metabolismo , Acetilcoenzima A/metabolismo , Teorema de Bayes , Biocombustíveis/microbiologia , Anidrases Carbônicas/metabolismo , Homosserina Desidrogenase/metabolismo , Engenharia Metabólica/métodos , Pirofosfatases/metabolismoRESUMO
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.
RESUMO
BACKGROUND: Mitigation of climate change requires that new routes for the production of fuels and chemicals be as oil-independent as possible. The microbial conversion of lignocellulosic feedstocks into terpene-based biofuels and bioproducts represents one such route. This work builds upon previous demonstrations that the single-celled carotenogenic basidiomycete, Rhodosporidium toruloides, is a promising host for the production of terpenes from lignocellulosic hydrolysates. RESULTS: This study focuses on the optimization of production of the monoterpene 1,8-cineole and the sesquiterpene α-bisabolene in R. toruloides. The α-bisabolene titer attained in R. toruloides was found to be proportional to the copy number of the bisabolene synthase (BIS) expression cassette, which in turn influenced the expression level of several native mevalonate pathway genes. The addition of more copies of BIS under a stronger promoter resulted in production of α-bisabolene at 2.2 g/L from lignocellulosic hydrolysate in a 2-L fermenter. Production of 1,8-cineole was found to be limited by availability of the precursor geranylgeranyl pyrophosphate (GPP) and expression of an appropriate GPP synthase increased the monoterpene titer fourfold to 143 mg/L at bench scale. Targeted mevalonate pathway metabolite analysis suggested that 3-hydroxy-3-methyl-glutaryl-coenzyme A reductase (HMGR), mevalonate kinase (MK) and phosphomevalonate kinase (PMK) may be pathway bottlenecks are were therefore selected as targets for overexpression. Expression of HMGR, MK, and PMK orthologs and growth in an optimized lignocellulosic hydrolysate medium increased the 1,8-cineole titer an additional tenfold to 1.4 g/L. Expression of the same mevalonate pathway genes did not have as large an impact on α-bisabolene production, although the final titer was higher at 2.6 g/L. Furthermore, mevalonate pathway intermediates accumulated in the mevalonate-engineered strains, suggesting room for further improvement. CONCLUSIONS: This work brings R. toruloides closer to being able to make industrially relevant quantities of terpene from lignocellulosic biomass.
RESUMO
Counterfactual inference is a useful tool for comparing outcomes of interventions on complex systems. It requires us to represent the system in form of a structural causal model, complete with a causal diagram, probabilistic assumptions on exogenous variables, and functional assignments. Specifying such models can be extremely difficult in practice. The process requires substantial domain expertise, and does not scale easily to large systems, multiple systems, or novel system modifications. At the same time, many application domains, such as molecular biology, are rich in structured causal knowledge that is qualitative in nature. This article proposes a general approach for querying a causal biological knowledge graph, and converting the qualitative result into a quantitative structural causal model that can learn from data to answer the question. We demonstrate the feasibility, accuracy and versatility of this approach using two case studies in systems biology. The first demonstrates the appropriateness of the underlying assumptions and the accuracy of the results. The second demonstrates the versatility of the approach by querying a knowledge base for the molecular determinants of a severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-induced cytokine storm, and performing counterfactual inference to estimate the causal effect of medical countermeasures for severely ill patients.
RESUMO
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.
RESUMO
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.
RESUMO
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.
RESUMO
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.
Assuntos
Ritmo Circadiano , Proteínas Fúngicas/genética , Regulação Fúngica da Expressão Gênica , Redes e Vias Metabólicas , Neurospora crassa/genética , Saccharomyces cerevisiae/genética , Relógios Circadianos , Proteínas Fúngicas/metabolismo , Neurospora crassa/metabolismo , Proteômica , Saccharomyces cerevisiae/metabolismo , Transcrição GênicaRESUMO
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.
RESUMO
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.