Skip to main content
  • ASM Journals
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Latest Articles
    • Special Issues
    • COVID-19 Special Collection
    • Editor's Picks
    • Special Series: Sponsored Minireviews and Video Abstracts
    • Archive
  • Topics
    • Applied and Environmental Science
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Novel Systems Biology Techniques
    • Early-Career Systems Microbiology Perspectives
  • For Authors
    • Getting Started
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics
  • About the Journal
    • About mSystems
    • Editor in Chief
    • Board of Editors
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • ASM Journals
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
mSystems
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Latest Articles
    • Special Issues
    • COVID-19 Special Collection
    • Editor's Picks
    • Special Series: Sponsored Minireviews and Video Abstracts
    • Archive
  • Topics
    • Applied and Environmental Science
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Novel Systems Biology Techniques
    • Early-Career Systems Microbiology Perspectives
  • For Authors
    • Getting Started
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics
  • About the Journal
    • About mSystems
    • Editor in Chief
    • Board of Editors
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
Research Article | Applied and Environmental Science

Experimentally Validated Reconstruction and Analysis of a Genome-Scale Metabolic Model of an Anaerobic Neocallimastigomycota Fungus

St. Elmo Wilken, Jonathan M. Monk, Patrick A. Leggieri, Christopher E. Lawson, Thomas S. Lankiewicz, Susanna Seppälä, Chris G. Daum, Jerry Jenkins, Anna M. Lipzen, Stephen J. Mondo, Kerrie W. Barry, Igor V. Grigoriev, John K. Henske, Michael K. Theodorou, Bernhard O. Palsson, Linda R. Petzold, Michelle A. O’Malley
Stephen R. Lindemann, Editor
St. Elmo Wilken
aDepartment of Chemical Engineering, University of California Santa Barbara, Santa Barbara, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jonathan M. Monk
bDepartment of Bioengineering, University of California San Diego, San Diego, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Patrick A. Leggieri
aDepartment of Chemical Engineering, University of California Santa Barbara, Santa Barbara, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Christopher E. Lawson
cJoint BioEnergy Institute, Lawrence Berkeley National Laboratory, Emeryville, California, USA
dBiological Systems and Engineering Division, Lawrence Berkeley National Laboratory, Berkeley, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Thomas S. Lankiewicz
cJoint BioEnergy Institute, Lawrence Berkeley National Laboratory, Emeryville, California, USA
eDepartment of Evolution Ecology and Marine Biology, University of California Santa Barbara, Santa Barbara, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Susanna Seppälä
aDepartment of Chemical Engineering, University of California Santa Barbara, Santa Barbara, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Chris G. Daum
fUS Department of Energy Joint Genome Institute, Lawrence Berkeley National Laboratory, Berkeley, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jerry Jenkins
fUS Department of Energy Joint Genome Institute, Lawrence Berkeley National Laboratory, Berkeley, California, USA
gHudsonAlpha Institute of Biotechnology, Huntsville, Alabama, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Anna M. Lipzen
fUS Department of Energy Joint Genome Institute, Lawrence Berkeley National Laboratory, Berkeley, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Stephen J. Mondo
fUS Department of Energy Joint Genome Institute, Lawrence Berkeley National Laboratory, Berkeley, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Kerrie W. Barry
fUS Department of Energy Joint Genome Institute, Lawrence Berkeley National Laboratory, Berkeley, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Igor V. Grigoriev
fUS Department of Energy Joint Genome Institute, Lawrence Berkeley National Laboratory, Berkeley, California, USA
hDepartment of Plant and Microbial Biology, University of California Berkeley, Berkeley, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
John K. Henske
aDepartment of Chemical Engineering, University of California Santa Barbara, Santa Barbara, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Michael K. Theodorou
iAgriculture Centre for Sustainable Energy Systems, Animal Production, Department of Agriculture and the Environment, Harper Adams University, Newport, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Bernhard O. Palsson
bDepartment of Bioengineering, University of California San Diego, San Diego, California, USA
cJoint BioEnergy Institute, Lawrence Berkeley National Laboratory, Emeryville, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Bernhard O. Palsson
Linda R. Petzold
jDepartment of Computer Science, University of California Santa Barbara, Santa Barbara, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Michelle A. O’Malley
aDepartment of Chemical Engineering, University of California Santa Barbara, Santa Barbara, California, USA
cJoint BioEnergy Institute, Lawrence Berkeley National Laboratory, Emeryville, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Michelle A. O’Malley
Stephen R. Lindemann
Purdue University
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Snorre Sulheim
SINTEF Industry
Roles: ad hoc peer reviewer
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Snorre Sulheim
Rajib Saha
University of Nebraska-Lincoln
Roles: ad hoc peer reviewer
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Rajib Saha
DOI: 10.1128/mSystems.00002-21
  • Article
  • Figures & Data
  • Info & Metrics
  • Review History
  • PDF
Loading

ABSTRACT

Anaerobic gut fungi in the phylum Neocallimastigomycota typically inhabit the digestive tracts of large mammalian herbivores, where they play an integral role in the decomposition of raw lignocellulose into its constitutive sugar monomers. However, quantitative tools to study their physiology are lacking, partially due to their complex and unresolved metabolism that includes the largely uncharacterized fungal hydrogenosome. Modern omics approaches combined with metabolic modeling can be used to establish an understanding of gut fungal metabolism and develop targeted engineering strategies to harness their degradation capabilities for lignocellulosic bioprocessing. Here, we introduce a high-quality genome of the anaerobic fungus Neocallimastix lanati from which we constructed the first genome-scale metabolic model of an anaerobic fungus. Relative to its size (200 Mbp, sequenced at 62× depth), it is the least fragmented publicly available gut fungal genome to date. Of the 1,788 lignocellulolytic enzymes annotated in the genome, 585 are associated with the fungal cellulosome, underscoring the powerful lignocellulolytic potential of N. lanati. The genome-scale metabolic model captures the primary metabolism of N. lanati and accurately predicts experimentally validated substrate utilization requirements. Additionally, metabolic flux predictions are verified by 13C metabolic flux analysis, demonstrating that the model faithfully describes the underlying fungal metabolism. Furthermore, the model clarifies key aspects of the hydrogenosomal metabolism and can be used as a platform to quantitatively study these biotechnologically important yet poorly understood early-branching fungi.

IMPORTANCE Recent genomic analyses have revealed that anaerobic gut fungi possess both the largest number and highest diversity of lignocellulolytic enzymes of all sequenced fungi, explaining their ability to decompose lignocellulosic substrates, e.g., agricultural waste, into fermentable sugars. Despite their potential, the development of engineering methods for these organisms has been slow due to their complex life cycle, understudied metabolism, and challenging anaerobic culture requirements. Currently, there is no framework that can be used to combine multi-omic data sets to understand their physiology. Here, we introduce a high-quality PacBio-sequenced genome of the anaerobic gut fungus Neocallimastix lanati. Beyond identifying a trove of lignocellulolytic enzymes, we use this genome to construct the first genome-scale metabolic model of an anaerobic gut fungus. The model is experimentally validated and sheds light on unresolved metabolic features common to gut fungi. Model-guided analysis will pave the way for deepening our understanding of anaerobic gut fungi and provides a systematic framework to guide strain engineering efforts of these organisms for biotechnological use.

INTRODUCTION

Anaerobic gut fungi in the early-branching phylum Neocallimastigomycota are found in the digestive tracts of large mammalian herbivores, where they play an integral role in the lignocellulolytic microbiome of their host (1). These fungi have an unusual lifestyle involving both a vegetative state and motile zoospores, and contain mitochondrion-like hydrogenosomes that are similar to organelles observed in other anaerobic eukaryotes (2). Recent transcriptomic and genomic analyses revealed that these fungi harbor an incredible diversity of carbohydrate-active enzymes (CAZymes) that are tailored to excel at decomposing lignocellulosic plant biomass (3–7). Moreover, these fungi appear to organize their biomass-degrading enzymes into unique eukaryotic cellulosomes, which are large and dynamic extracellular complexes of enzymes that likely contribute to the efficiency of biomass degradation by substrate channeling and modularity (4). Given their ability to metabolize raw lignocellulose, anaerobic gut fungi are an appealing biotechnological platform to drive the conversion of lignocellulose into hydrolyzed sugars and, ultimately, into renewable chemicals via fermentation (8–10).

While the lignocellulolytic capabilities of anaerobic gut fungi motivate their biotechnological interest, they are temperature sensitive, anaerobic, relatively slow growing, and hindered by requirements for specialized media. Moreover, their genomes are extremely AT and repeat rich, which makes sequencing and genetic engineering challenging (7, 11, 12). No robust genetic engineering tools have been developed for this class of fungi, hampering classic molecular biology techniques that can be used to investigate, understand, and engineer their metabolism. Despite these challenges, experimental and omics data sets have emerged to elucidate some aspects of their metabolism (7, 13–16). However, a framework to combine these data in a systematic manner to understand and engineer anaerobic gut fungal metabolism for industrial applications is not available. In particular, there is no clear consensus on the pathways operating in their hydrogenosome, a mitochondrion-like organelle involved in their energy and hydrogen metabolism. Current hypotheses suggest that either an energetically unfavorable pathway involving pyruvate formate lyase is used to produce H2 or a pathway involving pyruvate ferredoxin oxidoreductase that is not supported by extracellular metabolite measurements is used (14, 16).

Genome-scale metabolic models (GEMs) can be used to address these shortcomings, as they are well suited to act as knowledge base platforms for integrating multi-omic data sets and have been successfully used to drive the engineering of both prokaryotes and eukaryotes (17–19). Moreover, by experimentally testing the predictions of a GEM, it is possible to systematically refine the understanding of the metabolism of an organism. This ability to probe the metabolism in silico is particularly appealing in the context of nonmodel microbes, such as the anaerobic gut fungi, where direct metabolic manipulation is challenging.

Here, we introduce a high-quality PacBio-sequenced genome (200 Mbp, 62× sequencing depth) of the anaerobic gut fungus Neocallimastix lanati. Comparative genomic analyses revealed that N. lanati is metabolically similar to the other sequenced Neocallimastigomycota, including the presence of many CAZymes (∼1,788 CAZymes, with 585 associated with the fungal cellulosome), suggesting that insights gained from understanding its genome may be generalizable to other species in the clade. Therefore, we used the genome of N. lanati to construct the first genome-scale metabolic model of an anaerobic gut fungus. This fungus is well suited to act as a model system to investigate the metabolism of anaerobic gut fungi because, unlike many of them, it grows relatively well in completely defined (M2) medium (13, 20), it is relatively fast growing among gut fungal strains (μ = 0.045 ± 0.003 h−1 in M2 medium) (21–23), and it can be cryopreserved.

The 3-compartment (extracellular, cytosolic, and hydrogenosomal compartments) model introduced here, named iNlan20, is composed of 1,018 genes, 1,023 reactions, and 816 metabolites and models the primary metabolism of N. lanati. The model is stoichiometrically consistent as well as mass and charge balanced. Experimental, genomic, transcriptomic, and metabolic flux analysis data were used to build and validate the model, which recapitulates extracellular metabolite production rates and accurately models the observed fungal growth rate. Furthermore, the model refines and expands on previous hypotheses regarding the metabolism of the gut fungal hydrogenosome. Both the model and experimental data suggest that pyruvate formate lyase (PFL) is significantly more active than pyruvate ferredoxin oxidoreductase (PFO) in the hydrogenosome but that hydrogen formation can only occur via the latter pathway. Going forward, this fungus and its associated model can be used to guide efforts to further refine aspects of gut fungal metabolism that remain unclear and direct metabolic engineering strategies. Indeed, model-based analysis could be invaluable in designing stable consortia between anaerobic gut fungi and other industrially utilized organisms—something that has not yet been fully realized (21, 24, 25).

RESULTS

The genome of N. lanati is rich in carbohydrate-active enzymes and metabolically similar to other anaerobic gut fungi.Given the large repeat-rich genomes inherent to anaerobic gut fungi (12), PacBio sequencing was used to obtain a high-quality genome of the isolate N. lanati; this isolate was sourced from a fecal pellet of a sheep obtained from an enclosure at the Santa Barbara Zoo (the Index Fungorum identification number is IF557810). While the genome of this fungus is large (Table 1), it is the second least fragmented of all 5 published gut fungal genomes (see Table S1 in the supplemental material). The N. lanati genome encodes a rich array of carbohydrate-active enzymes (CAZymes) in numbers similar to those reported from other gut fungal genomes (Table S1) (3, 4, 7). In total, 1,788 CAZymes were identified in the genome, of which, 1,253 were expressed in the transcriptome. Like other anaerobic gut fungi, N. lanati deploys both complexed (cellulosomes) and uncomplexed CAZymes through its rhizoidal network (Table 1; Fig. 1), both of which contribute to the decomposition of lignocellulose. Figure S1 shows the breakdown of CAZyme domains identified in the genome of N. lanati. This, in combination with its relatively high growth rate on defined M2 media, suggests that N. lanati is a good model anaerobic gut fungus.

FIG 1
  • Open in new tab
  • Download powerpoint
FIG 1

The morphology of Neocallimastix lanati aids in the decomposition of unpretreated lignocellulose by disrupting the lignocellulosic plant biomass to increase the surface area available for enzymatic attack. A micrograph of a mature N. lanati sporangium growing on corn stover in M2 medium after 3 days of growth at 39°C. The filamentous rhizoidal network is used to increase the surface area for its lignocellulolytic enzymes that decompose the lignocellulosic corn stover into its fermentable sugar constituents.

View this table:
  • View inline
  • View popup
  • Download powerpoint
TABLE 1

Summary of the features of the genome of N. lanatia

TABLE S1

Annotation data for each publicly available anaerobic gut fungal genome shows that their genomes are large and fragmented and that many of the putative genes have not been assigned a function. Data taken from the MycoCosm (excludes bidirectional blast assigned annotations) (48). Here, annotated genes are defined as any gene that has at least one protein family (PF) domain assigned to it based on the annotation data supplied by the JGI database. Download Table S1, DOCX file, 0.01 MB.
Copyright © 2021 Wilken et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

FIG S1

The CAZyme repertoire (number of annotated domains shown in brackets across the genome) of N. lanati include numerous glycoside hydrolases (GHs), dockerins (DOCs), and carbohydrate binding modules (CBMs). These domains facilitate the hydrolysis of cellulose- and hemicellulose-based substrates, suggesting that N. lanati has an expansive cellulolytic capability. PL, polysaccharide lyase; EXPN, expansin; GT, glycosyltransferase; CE, carbohydrate esterase. Annotation data for N. lanati can be found at https://mycocosm.jgi.doe.gov/Neolan1/Neolan1.info.html. Download FIG S1, TIF file, 0.2 MB.
Copyright © 2021 Wilken et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Despite advances in sequencing and annotation, a large number of putative gut fungal genes remain completely unannotated (∼48% of the 27,677 predicted genes of N. lanati) (Table S1), which is consistent with previous genomic annotations in this clade. These unannotated genes contributed to gaps found in the draft reconstructed metabolism of N. lanati. A comparative genomic analysis within the primary metabolism across all high-quality publicly available gut fungal genomes (Anaeromyces robustus, Neocallimastix californiae, Pecoramyces ruminatium, and Piromyces finnis) revealed that the gut fungi are metabolically similar. Of the 1,023 unique EC numbers identified across these 5 genomes, fewer than ∼3% are unique to each isolate (Fig. 2). This suggests that gut fungi share a similar primary metabolism. Thus, some metabolic gaps in the draft reconstruction could be filled by searching for genes in N. lanati that are homologous to those encoded, and annotated, in the genomes of the other gut fungal isolates. In this way, key enzymes in the biosynthesis pathways of arginine, asparagine, biotin, riboflavin, lipids, and fatty acids were identified and included in the metabolic reconstruction of N. lanati. In total, 35 gaps in the primary metabolic pathways were identified and annotated in this manner, as noted in the confidence score and homologous gene annotation fields in the model.

FIG 2
  • Open in new tab
  • Download powerpoint
FIG 2

Anaerobic gut fungi have very similar genetic metabolic potentials, suggesting that metabolic gaps can be filled by looking for homologous genes found in the other sequenced isolates. Each Venn diagram was generated by inspecting the intersection of the annotated EC numbers contained in the genome of each fungus for each metabolic module. Overlapping regions imply that those isolates share the EC assignments contained in each of the metabolic modules. The EC numbers contained in each module are based the KEGG database (60) (see supplemental data 4 in the iNlan20 GitHub repository available at https://github.com/stelmo/iNlan20 for the list of modules encompassing each Venn diagram), while the EC assignments for each fungus are based on the JGI and bidirectional annotation data as described in Materials and Methods.

The curated metabolic model of N. lanati captures the carbon, amino acid, vitamin, fatty acid, nucleotide, and lipid metabolism.Based on the draft metabolic reconstruction of N. lanati, a manually curated genome-scale metabolic model of N. lanati was built (iNlan20) according to an established protocol for generating high-quality reconstructions (26). The model contains 1,023 reactions, 816 metabolites, and 1,018 genes distributed across 3 compartments (hydrogenosome, cytoplasm, and extracellular space). Where possible, experimental data were used to curate the model. Materials and Methods details specifics on the curation process, as well as experiments used to construct the biomass objective function of the model. Briefly, Table 2 shows the experimentally measured macromolecular components of N. lanati that were used to construct the biomass objective function for the genome-scale metabolic model. Further simplifying assumptions were made to construct the specific biomass objective function used in iNlan20. The insoluble carbohydrate component of the biomass was assumed to be solely chitin (27), and the amino acid composition of the protein component of the biomass was assumed to follow the amino acid distribution of the predicted genes (i.e., the predicted proteome). Similarly, the nucleotide composition was assumed to follow the composition of the genome (for the DNA nucleotides) and the transcriptome (for the RNA nucleotides). The lipid component was assumed to be composed of myristic, palmitic, and stearic acids, which were found to be the major fatty acid components of the lipid fraction of N. lanati (see Fig. S2). The growth-associated and non-growth-associated maintenance (GAM and NGAM, respectively) functions were estimated using experimental data (see Fig. S3).

View this table:
  • View inline
  • View popup
  • Download powerpoint
TABLE 2

Experimentally measured macromolecular constituents of N. lanati that were used to construct the biomass objective function for the genome-scale metabolic modela

FIG S2

Fatty acid composition of N. lanati reveals that the saturated C14:0, C16:0, and C18:0 fatty acids are the primary fatty acids produced by the gut fungus. The relative concentration of each fatty acid in lyophilized biomass samples of N. lanati after 4 days of growth in M2 at 39°C using cellobiose as a carbon source was determined by reconstituting the extracted lipid component into a 65% acetonitrile, 30% isopropanol, and 5% water solution. Mass spectrometry analysis, as detailed in Materials and Methods, was used to determine the concentrations. Only the relative concentrations were used in the biomass objective function. Download FIG S2, TIF file, 0.09 MB.
Copyright © 2021 Wilken et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

FIG S3

Growth- and non-growth-associated maintenance (GAM and NGAM, respectively) functions were estimated from experimental data for N. lanati. The GAM function was found to be 75.98 mmol ATP/g (dry weight)/h, while the NGAM function was found to be 2.27 mmol ATP/g (dry weight). While there is a large degree of uncertainty in the fit, these estimates are in line with those of other fungal GEMs (GAM, ∼15 to 86.6 mmol ATP/g [dry weight]/h; NGAM, ∼<1 to 18.7 mmol ATP/g [dry weight]) (P. Pan and Q. Hua, PLoS One, 7:e51535, 2012, https://doi.org/10.1371/journal.pone.0051535; M. Tomàs-Gamisans, P. Ferrer, and J. Albiol, PLoS One, 11:e0148031, 2016, https://doi.org/10.1371/journal.pone.0148031; J. Villadsen, J. Nielsen, and G. Liden (ed), Bioreaction engineering principles, 3rd edition, 2011, Springer, Berlin, Germany; W. Vongsangnak, A. Klanchui, I. Tawornsamretkit, et al., Gene, 583:121–129, 2016, https://doi.org/10.1016/j.gene.2016.02.028; C. Ye, W. Qiao, X. Yu, et al., BMC Genomics, 16:799, 2015, https://doi.org/10.1186/s12864-015-2042-y) and were assumed sufficiently informative for this model. Increasing concentrations of cellobiose were used to promote faster growth, inferred by pressure accumulation measurements. Metabolite measurements were taken 24 h apart, and the resultant fluxes were used to constrain the model. This constrained model was used to predict the maximum ATP yield given a specific growth rate. See Materials and Methods for more details. Download FIG S3, TIF file, 0.1 MB.
Copyright © 2021 Wilken et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

The model focuses on the primary metabolism but includes CAZymes as generalized cellulase and hemicellulase reactions. Of the 791 metabolic genes included in the model, 216 do not have gene assignments—reflecting the understudied nature of gut fungal metabolism. Despite this, the model is stoichiometrically consistent as well as mass and charged balanced (see Memote report or data set S1 in the iNlan20 GitHub repository available at https://github.com/stelmo/iNlan20). Table S2 explains the confidence rating assigned to the genes associated with each reaction in the model. In the energy-generating pathways, particular attention was paid to modeling the hydrogenosome (a mitochondrion-like organelle that functions completely anaerobically), which is discussed in greater detail in the following sections. More generally, the Embden-Meyerhof-Parnas variant of glycolysis is present in N. lanati as well as pathways for mixed-acid fermentation (succinate, acetate, lactate, formate, and ethanol), which are typically found in anaerobic gut bacteria (28). Interestingly, it was found that N. lanati possesses both the NAD+ and NADP+ variants of glyceraldehyde-3-phosphate dehydrogenase in glycolysis, with the latter used to conserve energy as NADPH instead of ATP. The pentose phosphate pathway of N. lanati is incomplete, with glucose-6-phosphate dehydrogenase and 6-phosphogluconate missing. These reactions regenerate NADPH and possibly explain the presence of the NADP+ variant of glyceraldehyde-3-phosphate dehydrogenase as a compensating mechanism (29, 30). Despite these missing genes, the model is still capable of producing nucleotide precursors. The xylose isomerase pathway is also present in N. lanati, as has been found in other sequenced anaerobic gut fungi (21).

TABLE S2

Confidence scoring system used for each reaction in the model suggests that there is genetic evidence for more than 70% of the reactions included in the model. The confidence score is annotated in the model under “notes”. Download Table S2, DOCX file, 0.01 MB.
Copyright © 2021 Wilken et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

The major components (amino acids, nucleotides, vitamins, fatty acids, and lipids) of the anabolic metabolism of N. lanati were found to be present, in agreement with its ability to grow in M2 medium without these components added. Specifically, the complete biosynthesis pathways for all the proteogenic amino acids and the modeled fatty acids were found. Most of the canonical vitamin and cofactor (vitamin B5, vitamin B6, riboflavin, and thiamine) biosynthesis pathways were found to be complete, with the exception of folate, where no synthesis mechanism of 4-aminobenzoate was found. The heme and biotin biosynthesis pathways appeared to be incomplete; however, the latter vitamin was not found to be essential (experimentally) in defined medium, suggesting that the pathway is either poorly annotated or that the fungus does not require it for growth. Finally, the model recapitulates the experimentally observed growth rate in defined medium using only the measured flux of glucose (1.5 mmol/g [dry weight]/h) as an input constraint (flux balance analysis predicted μ = 0.044 h−1 versus an experimentally measured μ = 0.045 ± 0.003 h−1).

iNlan20 includes an expanded model of the hydrogenosomal metabolism.Anaerobic gut fungi possess a variant of the hydrogenosome, with the core set of enzymes that catalyze the conversion of malate and pyruvate to acetate, H2 and formate already identified, as shown in Fig. 3 (2, 14, 16, 31). However, the metabolic pathway leading to H2 production is not resolved, with literature suggesting either pyruvate ferredoxin oxidoreductase (PFO) or pyruvate formate lyase (PFL) as possible routes (Fig. 3). Both enzymes were identified in the genome and transcriptome and are thus included in the model of the hydrogenosome.

FIG 3
  • Open in new tab
  • Download powerpoint
FIG 3

An expanded model of the hydrogenosome is included in the model based on genomic annotation, literature, and predicted localization data (14–16). Core hydrogenosome enzymes are colored in blue, while speculative enzymes are shown in black. PFL, pyruvate formate lyase; PFO, pyruvate ferredoxin oxidoreductase; Ac, acetate; SucCoA, succinyl coenzyme A; CoA, coenzyme A; AcCoA, acetyl coenzyme A; Frdx, ferredoxin.

By combining literature sources, gene annotation, transcriptomic expression, and subcellular localization data, we have included additional pathways in the model of the hydrogenosome for N. lanati (Table 3 and Fig. 3). Of note is the inclusion of a putative ATP synthase (32, 33), which was previously speculated to be present in other anaerobic gut fungal isolates (7, 16). Additionally, we also found evidence that complex 2 of the mitochondrial electron transport chain is present: homologs to all four subunits were found to be expressed and localized to the hydrogenosome (complex 2, subunits A, B, C, and D). We could not find any homologs of the membrane-bound subunits of complex 1 or the ATP synthase in the N. lanati genome, as was also reported previously for other anaerobic gut fungi (34). A possible explanation for this is that many of the membrane-bound subunits of the electron transport chain are encoded by mitochondrial DNA, which the fungal hydrogenosome appears to have lost. Despite this, it is perplexing that no homologs of the membrane-bound subunits of complex 1 were found, since these are used to shuttle electrons between the two complexes in the inner membrane of the mitochondrion, and it remains unclear how complex 2 could function without them. However, homologs of the soluble subunits of complex 1, nuoF and nuoE, are highly expressed relative to the other core enzymes of the hydrogenosome (Table 3). The presence of the soluble subunits, coupled with the absence of the membrane-associated subunits of complex 1, has also been observed in the hydrogenosomes of, e.g., Trichomonas vaginalis (35, 36) and Nyctotherus ovalis (37). This raises two possibilities. First, that N. lanati possesses a proton-pumping mechanism, but the missing membrane-bound genes are unannotated. While the membrane-bound subunits of complex 1 and the ATP synthase appear to be missing, we do find preliminary evidence of a pH gradient inside the hydrogenosome (see Fig. S4). This provides some support for a mechanism, possibly involving complex 1, 2, and the ATP synthase, that makes use of this gradient to generate energy, as found in other hydrogenosomes (32), if these genes are present but unannotated. Second, it is possible that its hydrogenase associates with the identified nuoF-like subunit of complex 1 to form a bifurcating hydrogenase, as has been speculated to occur in the hydrogenosome of T. vaginalis (32, 35). Indeed, we find high-homology sequences in the N. lanati genome to all three of the bifurcating hydrogenase subunits that have been enzymatically characterized in Thermotoga maritima (38) (see data set S2 in the iNlan20 GitHub repository available at https://github.com/stelmo/iNlan20). In this case, the function of the membrane-bound subunits of complex 2 and the soluble subunits of the ATP synthase are unknown. In either case, further experimental work is needed to clarify these questions.

View this table:
  • View inline
  • View popup
  • Download powerpoint
TABLE 3

Enzymes included in the model of the hydrogenosome metabolism

FIG S4

Zoospores were stained for the presence of intercellular organelles with electrical potential using the dye JC-1 as well as DAPI to illuminate the nuclear structure. Blue represents the DAPI stain, and orange/red represents the JC-1 stain, which indicates the presence of a pH gradient in the hydrogenosome (JC-1 is a mitochondrial selective stain, but here, the mitochondria is replaced by the hydrogenosome). Download FIG S4, TIF file, 0.7 MB.
Copyright © 2021 Wilken et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Taken together, our expanded model of the hydrogenosome includes the core enzymes previously reported in other gut fungal species as well as a speculative bifurcating hydrogenase, ATP synthase, and proton-pumping module composed of the complex 1 and complex 2 enzymes identified in the N. lanati genome (similar to what has been found in other H2-producing mitochondria [32]). Given the speculative nature of the proton-pumping mechanism, the ATP synthase, and the bifurcating hydrogenase, these reactions are constrained to carry zero flux in the base case model. These constraints were modified to investigate the consequences of this extended hydrogenosomal metabolism, as discussed later. Additionally, it was previously suggested that a hydrogen dehydrogenase [NAD(P)+ + H2 ↔ H+ + NAD(P)H] operates in the reverse direction in the hydrogenosome (7, 14, 15). Consequently, this hydrogen dehydrogenase simultaneously produces H2 and prevents the accumulation of NAD(P)H produced by the malic enzyme in the hydrogenosome. However, in this direction, the reaction is energetically very unfavorable, with a ΔG of ≈34 ± 5.9 kJ/mol assuming physiologically realistic conditions. Therefore, the flux bounds of this reaction in the hydrogenosome were set to reflect the assumption that the hydrogen dehydrogenase only carries flux in the forward, energetically favorable, direction.

iNlan20 accurately predicts substrate utilization and in vivo fluxes.The curated model was validated using a combination of growth curve, extracellular metabolite, and metabolic flux analysis (MFA) data. Substrate utilization tests were performed on 36 different carbon sources, focusing on metabolites that are present in the digestive tract environment of the anaerobic gut fungi (Table 4). The qualitative prediction accuracy of the model for the substrate utilization and vitamin essentiality validation tests is 89% (Matthews correlation coefficient, 0.79). Interestingly, despite the apparent presence of a full xylose isomerase pathway, N. lanati did not grow using xylose as its sole carbon source, as has been found in other gut fungi (21). In this case, the model’s predictions were incorrect. It was previously suggested that transport limitations may cause this issue (21). However, a xylose transporter was identified in the genome, suggesting that cellular regulation might explain this discrepancy better. Vitamin essentiality tests were also conducted (Table 4). It was found that both heme and 4-aminobenzoate were essential for growth, in agreement with the model’s predictions. In other gut fungi, heme has also been found to be essential (39), suggesting that its de novo biosynthesis pathway may be absent across the clade. It was found that only cysteine could be used as a sulfur source. However, it is not clear if this is a nutritional requirement, since every other reducing agent tested (Na2S, 2-mercaptoethanol, and dithiothreitol) appeared to be toxic to the fungus. Since cysteine was used to ensure anaerobicity of the medium, we could not test nitrogen source utilization.

View this table:
  • View inline
  • View popup
  • Download powerpoint
TABLE 4

Substrate utilization results suggest that the model accurately captures phenotypic behavior of N. lanatia

Metabolic flux analysis (MFA) was also used to experimentally verify the predicted intracellular fluxes of the GEM. A 1,2-13C-labeled glucose tracer was used in conjunction with a carbon atom transition model built from the N. lanati metabolic reconstruction. For the MFA model, metabolic degeneracy caused by the ability of the hydrogenosome to metabolize both malate and pyruvate resulted in large bounds on the fluxes involving these metabolites. To circumvent this, the MFA model was constrained to only import pyruvate into the hydrogenosome, based on previous observations (14). Extracellular metabolic product measurements (ethanol, formate, H2, acetate, succinate, and lactate) were also used to constrain the MFA model. This resulted in accurate internal metabolic flux measurements based on a statistically significant fit between measured and simulated proteinogenic amino acid labeling patterns (Fig. 4). These 13C measured fluxes were then compared to the fluxes predicted using the GEM with independently measured metabolite flux constraints (Table 5). Parsimonious flux-based analysis (pFBA) was then used to find unique flux predictions. Using these constraints, the coefficient of determination between the pFBA and MFA simulation was found to be 0.98 (linear regression fit P < 0.01) (see Fig. S5). This suggests that the constrained metabolic model accurately predicts the steady-state measured intracellular fluxes of N. lanati.

FIG 4
  • Open in new tab
  • Download powerpoint
FIG 4

The genome-scale metabolic model accurately predicts the in vivo carbon metabolism of N. lanati. Experimentally determined MFA fluxes and predicted pFBA fluxes (top and bottom, respectively) for glycolysis, the TCA cycle, and the hydrogenosome of N. lanati. Error estimates denote one standard deviation from the reported mean for the MFA measurements. Three serially passaged [1,2-13C]glucose tracer experiments, grown in M2 medium at 39°C and harvested during exponential phase, were used to measure the in vivo fluxes (see Materials and Methods for more details and the model).

View this table:
  • View inline
  • View popup
  • Download powerpoint
TABLE 5

Experimentally measured external fluxes of various metabolites produced by N. lanati growing on cellobiose in M2 medium during exponential phasea

FIG S5

Metabolic fluxes as measured by a 13C-labeled glucose tracer suggests that the model accurately captures internal cellular fluxes. Normalized parsimonious FBA (pFBA) fluxes and their corresponding measured MFA fluxes were used to compare the accuracy of the model’s predictions of the internal fluxes of N. lanati. The pFBA model was constrained by its exchange fluxes for glucose, H2, lactate, ethanol, formate, and acetate to reflect experimentally observed data. The MFA data were generated by growing N. lanati in M2 medium at 39°C and harvesting the cell pellet during exponential phase. The regression fit is significant (P < 0.01). Further measured flux analysis was conducted as detailed in Materials and Methods and shown in the included MFA model found in the GitHub repository of the model available at https://github.com/stelmo/iNlan20. Download FIG S5, TIF file, 0.06 MB.
Copyright © 2021 Wilken et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

The core hydrogenosome metabolism uses PFO to produce hydrogen, but PFL carries the most flux.There remains uncertainty regarding the presence of pyruvate ferredoxin oxidoreductase (PFO) and its relative importance in the hydrogenosomal metabolism of anaerobic fungi. Earlier enzymatic characterization of hydrogenosomal proteins in Neocallimastigomycota suggested that PFO is the primary route for H2 production through an associated ferredoxin hydrogenase, as found in the hydrogenosomes of other organisms (16, 31, 32, 40). However, more recent studies suggest that PFO is either absent or of only marginal importance in the gut fungal hydrogenosomal metabolism (14, 15). These later studies suggest that pyruvate formate lyase (PFL), which was likely acquired through horizontal gene transfer from bacteria (41), is significantly more active than PFO. It has been suggested that hydrogen evolution occurs through a hydrogen dehydrogenase working in an energetically infeasible reverse direction (7, 14). Both PFO and PFL were identified in all published gut fungal genomes as well as in N. lanati (Table 3). The model was used to reconcile the role and relative importance of these two enzymes to hydrogenosome function under steady-state growth conditions.

Due to the reaction stoichiometry of PFL, the molar ratio of formate to acetate and ethanol produced is expected to approach unity (1:1) if PFL is metabolically dominant. This is because ethanol and acetate are only produced from PFL (via acetyl coenzyme A) in the cytosol (15). Since PFO only produces acetate, and not formate, the ratio of formate to acetate and ethanol will not be unity if PFO carries significant metabolic flux. Figure S6A shows that the experimentally measured molar ratios of formate to acetate and ethanol are not significantly different (using the unequal variance test, P < 0.05) from unity for N. lanati. This is in agreement with earlier metabolite measurements for Piromyces sp. strain E2, suggesting that PFL is dominant (14). Figure S6B shows that the unconstrained model predicts a wide range of possible ratios, reflecting the metabolic degeneracy of the carbon metabolism of N. lanati. Since there is no energetic cost associated with using PFO versus PFL (both produce one ATP molecule per pyruvate) (Fig. 3), the model predicts that both could be used to maximize ATP production in the hydrogenosome. However, external metabolite flux measurements show only modest H2 production (Table 5), suggesting that cellular regulation may play a role in diverting flux to PFL instead of PFO. This can also be seen in the relative expression difference between PFL and PFO (an order of magnitude difference between them) in Table 3. When the model is constrained by the measured metabolite fluxes shown in Table 5, the range of possible ratios is reduced to those observed experimentally (Fig. S6B). Since PFO is the only (known) energetically feasible way to produce H2, this result is not surprising. Using this constraint, the model suggests that PFL carries the most flux in the hydrogenosome, but that PFO is used to produce H2.

FIG S6

The molar ratio of produced formate to acetate and ethanol (denoted by the flux ratio in this figure) is not significantly different from unity, confirming that PFL plays a dominant role in the metabolism of the hydrogenosome. (A) Experimentally measured extracellular metabolites indicate that the ratio of formate to acetate and ethanol is close to unity, as has been shown to be the case for other gut fungi (14). (B) The unconstrained genome-scale metabolic model prediction of the formate to acetate and ethanol ratio covers a broad range. This suggests that cellular regulation plays a role in directing flux away from PFO in vivo. When the H2 flux is constrained, the model accurately captures the experimental data. For panel A, N. lanati was grown in M2 medium at 39°C using the noted cellobiose initial concentration, and the end point metabolite concentrations were measured using the HPLC protocol detailed in Materials and Methods. The molar ratios were then compared as specified. For panel B, 2,000 feasible flux solutions were sampled from the model that was unconstrained and constrained in the H2 flux variable. Download FIG S6, TIF file, 0.08 MB.
Copyright © 2021 Wilken et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Electron bifurcation and proton pumping may form part of the hydrogenosomal metabolism.Electron bifurcation is an energy conservation mechanism that can be used to drive thermodynamically unfavorable reactions by coupling endergonic and exergonic reactions through an enzyme complex (42). In the simplest case, this phenomenon is used by anaerobes to increase the yield of ATP through their carbon metabolism by using H2 as an electron sink for the recycling of NADH to NAD+ (42). In the case of the anaerobic gut fungi, the hydrogenosome can be used to generate 2 extra moles of ATP for every mole of glucose that enters glycolysis. However, not all the glycolytic flux can be diverted to the hydrogenosome, because NAD+ needs to be regenerated from the NADH that is produced by glycolysis to maintain cellular redox balance. As mentioned before, NAD+ is unlikely to be produced by the hydrogen dehydrogenase since the redox potential of NADH/NAD+ is too electropositive to reduce H+ directly (32). On the other hand, the ferredoxin-based hydrogenase included in the model only recycles the oxidized ferredoxin, produced by PFO, to reduced ferredoxin and does not impact the NADH/NAD+ pools in the cell. Sequence homology suggests that the N. lanati hydrogenosome could potentially house a bifurcating hydrogenase, which would couple the reduction of H+ to the oxidation of NADH through the ferredoxin produced by PFO. This hydrogenase enzyme complex would allow more flux to be channeled into the hydrogenosome for energy production, since the hydrogenase would now generate NAD+ as well as H2 (Fig. 3).

These findings suggest that there is a significant energetic advantage associated with possessing a bifurcating hydrogenase. The model captures this benefit by predicting a 16% increase in growth rate associated with the use of the bifurcating hydrogenase, as opposed to the ferredoxin hydrogenase (μ = 0.051 h−1 versus μ = 0.044 h−1, respectively). The model also dictates that the production of NAD+ shifts from the cytosol to the hydrogenosome when the bifurcating hydrogenase is used, as shown in Fig. 5. However, this requires metabolic flux to be diverted from PFL to PFO in the hydrogenosome. Consequently, significantly more H2 is predicted to be produced, which is not observed experimentally (Table S3). This discrepancy could be due to metabolic regulation that is unaccounted for in the GEM. Further experimental work needs to be conducted to investigate the potential presence of the bifurcating hydrogenase in the anaerobic gut fungal hydrogenosome.

FIG 5
  • Open in new tab
  • Download powerpoint
FIG 5

The effects of including additional reactions in the hydrogenosome on NAD+ and ATP production show that the putative bifurcating hydrogenase has a large positive effect on NAD+ and ATP generation. Conversely, the putative ATP synthase has a negligible effect on both. The increase in growth rate caused by the putative bifurcating hydrogenase is due to NAD+ being regenerated in the hydrogenosome; this allows more flux to be channeled into the organelle, which in turn produces more ATP. Flux sampling was used to determine the fluxes associated with NAD+ and ATP production in each metabolic configuration. The base case model only includes the ferredoxin hydrogenase. The ferredoxin hydrogenase was replaced by a bifurcating hydrogenase to analyze its effect on the model. Finally, a complex 1, 2, and ATP synthase module was added to the base case model to investigate the consequences of this expanded metabolism. The model was constrained to produce biomass at 90% of the maximum yield; subsequently, 2,000 samples were drawn from each case. The average production of each metabolite in the hydrogenosome and cytosol are shown.

TABLE S3

H2 flux predicted by the unconstrained model using the ferredoxin and bifurcating hydrogenases as well as the experimentally measured values. Predicted flux bounds are based on 2,000 samples of the relevant model; the biomass objective function was constrained to carry at least 90% of the flux of the optimal solution. Download Table S3, DOCX file, 0.01 MB.
Copyright © 2021 Wilken et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Given that the hydrogenosome in anaerobic fungi is relatively understudied yet related to the mitochondrion (32, 43), we assumed that their metabolite trafficking machinery is similar. Specifically, we assumed that the hydrogenosome has an inner membrane that is impermeable to H+, like those of mitochondria (44). This implies that H+ can only enter and leave the hydrogenosome through the action of transporters. The H+ balance in the hydrogenosome has a direct effect on the ability of the putative ATP synthase to produce ATP. However, the impact of the ATP synthase on ATP production was found to be small, with it only supporting small fluxes (∼3% of the glucose flux into the model) (Fig. 5). This suggests that the putative hydrogenosomal proton gradient (Fig. S4) may not be important for the generation of ATP, as is also suggested by the low expression of the ATP synthase complex subunits (Table 3), or that the proton gradient mechanism is not yet fully understood. Without further experimental evidence, it remains an open question whether the identified ATP synthase components are localized to the hydrogenosome or to some other subcellular organelle (34). Furthermore, the mechanism by which the putative complex 2 operates without the membrane-bound subunits of complex 1 remains to be determined.

Metabolic degeneracy is related to the redox balance.As is typical of unconstrained GEMs, the modeled gut fungal metabolism displays significant degeneracy, as shown by the high degree of flux variability (Table S4). The degeneracy is primarily due to the ability of N. lanati to regulate how NAD+ is regenerated through its mixed acid fermentation pathways, i.e., through a combination of lactate dehydrogenase, acetaldehyde dehydrogenase, and alcohol dehydrogenase. Interestingly, the relative mean error between the predicted flux distributions and experimental measurements of the fermentation products is much more sensitive to constraints placed on acetate production than any other single measured external metabolite flux, as shown in Fig. 6. Likewise, constraints placed on lactic acid flux also narrow the deviation of the predicted flux distributions. This effect is due to the different yields of NAD+ that can be achieved per mole of pyruvate depending on which mixed-acid fermentation pathway, or combination thereof, is constrained. Without the bifurcating hydrogenase, H2 production does not significantly impact the overall redox balance of the cell. This possibly explains why its constraint has the smallest effect on the flux variability predicted by the model and may allow the cell to fine tune its metabolism to suit the environmental needs, e.g., sugar availability, by up- or downregulating the flux channeled to the hydrogenosome (2).

FIG 6
  • Open in new tab
  • Download powerpoint
FIG 6

The absolute relative error between the model predictions and the experimentally measured values suggest that constraining the flux of acetate production has the biggest impact on the model’s accuracy. The flux of acetate (Ac), ethanol (EtOH), formate (For), H2, and lactate (Lac) was constrained, individually, to their observed ranges (variables on the x axis). The resultant predicted fluxes of these metabolites (generated by sampling 2,000 possible solutions where the biomass objective function was within 90% of its optimal value and subject to the respective additional constraints as shown in the figure) were then compared to the experimental observations as shown in the legend.

TABLE S4

Flux variability analysis indicates that there is significant variability in the production rates of metabolic products that all lead to near optimal growth rates (the optimal growth rate, μ = 0.044 h−1, was lower bounded to 90% of its optimum). Download Table S4, DOCX file, 0.01 MB.
Copyright © 2021 Wilken et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

DISCUSSION

While anaerobic gut fungi specialize in lignocellulose decomposition, they are not as well understood as model microbes such as Escherichia coli and Saccharomyces cerevisiae (45). The most recent GEMs of E. coli and S. cerevisiae map 1,515 and 1,150 genes to reactions, respectively (46, 47). In comparison, iNlan20 maps 567 metabolic genes to reactions, mainly within the primary metabolism, with the balance being CAZymes. The large number of unannotated predicted genes (∼13,300) (Table S1) likely form part of the fungal secondary metabolism (48). Moreover, of the remaining ∼14,400 genes annotated with at least a single domain, only 2,761 were associated with EC numbers (Table S1), highlighting how undercharacterized this early-branching clade of fungi is. In iNlan20, annotation gaps are centered around transporters and lipid metabolism. With improvements in gene annotation and more experimental characterization, this coverage is likely to increase. Various simplifications were also made during the model construction that may affect the predictions. First, it was assumed that the carbohydrate content of N. lanati is wholly composed of chitin; this is a simplification, as gut fungi make use of other carbohydrate moieties as well (e.g., glycosylation sugars [12, 27]). Second, the biomass function only accounts for fatty acid synthesis and no other lipids, resulting in the bulk of the blocked reactions and dead-end metabolites being associated with lipid metabolism. While cell wall and membrane biosynthesis are clearly essential building blocks in cellular homeostasis, these metabolites were not included in the objective function of the model because the gut fungal lipid composition is not well studied. Thus, expanding knowledge of lipid and cell wall anabolism will likely improve the model as well as possibly contributing to a deeper understanding of how the gut fungi secrete very large proteins (e.g., the cellulosome, ∼1 MDa [4]).

Compared to the MFA estimates, the model-based flux predictions were the most accurate around the central glycolytic pathways (Fig. 4). The largest relative discrepancies were observed in the tricarboxylic acid (TCA) cycle, specifically in the reductive branch leading to succinate production: 12 ± 3 versus 1 (normalized flux units) for MFA and pFBA, respectively. In eukaryotes, the TCA cycle is typically involved in producing a proton gradient within the mitochondrion, which in the context of the gut fungi, would suggest that it may be localized to the hydrogenosome. Model-based analysis indicated that the proton-pumping module and associated ATP synthase were of negligible importance (Fig. 5). Yet, the experimental data were more inconclusive: while the expression of these genes was comparatively low (Table 3), the JC-1 staining suggested that the hydrogenosome did possess an electrochemical gradient (Fig. S4). The absence of any identified membrane-bound subunits of complex 1 and the ATP synthase makes it hard to reconcile these data. These ambiguous observations suggest that experimental characterization of the gut fungal hydrogenosome (e.g., through enzyme purification and proteomics) is critically important for improving the model and answering lingering questions surrounding the metabolic role of the hydrogenosome.

Under the growth conditions analyzed here, neither the model predictions nor the experimental data suggested that PFO carries high flux. Consequently, low H2 production was expected and observed. However, literature sources indicate that stable gut fungal/H2-consuming methanogen cocultures are readily formed synthetically and observed in nature (25, 49, 50). Given the tight association between H2-consuming methanogens and gut fungi, it is surprising that the H2 flux is low and not more tightly coupled to the central metabolism. It is possible that there may be more complex cellular regulation at play within the hydrogenosome regarding PFO and PFL than modeled here. Interestingly, there is both a large size difference between PFO and PFL (∼539 kDa versus 90 kDa, respectively) and a catalytic rate difference (∼48 ± 38 s−1 versus 8 ± 5 s−1, respectively; data derived from BRENDA [51]). It is possible that under high sugar availability conditions, as used here, N. lanati produces the smaller slower enzyme (PFL) without compromising growth. In contrast, under more challenging conditions, e.g., when the fungus is using lignocellulose as a carbon source, it might produce more PFO. Alternatively, recent high-resolution growth rate data suggest that gut fungi have highly variable growth rates that fluctuate with environmental conditions (22). Thus, it is also possible that the gut fungi modulate the expression of PFO and, consequently H2 production during their life cycle. High-temporal-resolution proteomics and transcriptomics that examine various life stages may shed more light on this question.

Despite these discussed caveats, the iNlan20 shows good agreement with experimental measurements. Due to their nonmodel nature, anaerobic gut fungi are vastly understudied, which presents unique challenges when constructing a metabolic model of these cryptic organisms. Model-based analysis represents a systematic framework that can be used to identify high-impact knowledge gaps and focus experimental attention, and in this case, the model indicates that focus should be directed toward a detailed characterization of the gut fungal hydrogenosome. Within the realm of nonmodel microbes, metabolic modeling is an appealing technique that can be used to speed up the biotechnological translation of anaerobic gut fungi and other nonmodel microbes.

Conclusion.Here, we have introduced a high-quality genome and transcriptome of a novel anaerobic gut fungus, N. lanati. While the genome is large, it is relatively unfragmented compared to the genomes of the other sequenced anaerobic gut fungi. Additionally, the genome encodes a large number and diversity of CAZymes, most of which are expressed in the transcriptome. This genome was used to construct the first genome-scale metabolic model of an anaerobic gut fungus. The model, iNlan20, accurately recapitulates the observed growth rate, in vivo fluxes, and substrate consumption and requirement profiles. The model refines and expands on our understanding of gut fungal hydrogenosomal metabolism. We confirm previous findings that suggested that PFL carries more flux than PFO in the hydrogenosome, but an energetically favorable route to hydrogen production still requires the action of PFO. The possible presence of a bifurcating hydrogenase and/or a proton-pumping mechanism suggests that anaerobic fungi may have evolved more complex energy conservation mechanisms that allow them to compete with faster-growing rumen bacteria. Experimental work, likely involving the isolation, purification, and enzymatic characterization (through assays and proteomic analysis) of the hydrogenosome, is necessary to further refine our understanding of its metabolism. This model is well poised to serve as a platform to build a better understanding of these nonmodel organisms. Moreover, the model will serve as a valuable tool to systematically guide future engineering efforts of gut fungi for converting lignocellulose into value-added products.

MATERIALS AND METHODS

Metabolic reconstruction, visualization, and simulation.All publicly available annotated genomes within the clade Neocallimastigomycota were downloaded from the Joint Genome Institute’s (JGI) MycoCosm database (48). This includes the high-quality PacBio-sequenced genomes of Anaeromyces robustus, Piromyces finnis, and Neocallimastix californiae (4) as well as the novel isolate Neocallimastix lanati introduced here. The genomes of Pecoramyces ruminatium, also known as Orpinomyces sp. strain C1A (7, 52), and Piromyces sp. strain E2 (4) were also included for completeness. The gene annotation data supplied by the JGI was combined with annotations derived from bidirectionally searching by blast (using BLASTp [53]) the predicted genes from the gut fungal genomes against the curated Swiss-Prot database from UniProt (54). Briefly, bidirectional blast searching annotates a predicted gut fungal gene if (i) the top hit using the fungal genome as the query and the reference collection as the database is the same as when (ii) the gut fungal genome is used as the database and the reference collection is used as the query. Furthermore, only matches with E values smaller than 1e−20 were considered for assigning Enzyme Commission (EC) annotations to genes. This information was collated into a master metabolic table (see data set S3 in the iNlan20 GitHub repository available at https://github.com/stelmo/iNlan20) and subsequently used to construct the model and assign genes to reactions. Enzyme complexes were assigned by using the “Subunit structure” field in the UniProt database. Protein localization was predicted using DeepLoc (55). Reaction directions were primarily inferred from MetaCyc (56), and specific Gibbs free energy change of reactions reported were calculated using eQuilibrator (57). Transcriptomic and expression experiments for N. lanati were conducted as part of this study (described later). These omics data sets were used to assign a confidence score to each gene in the model of N. lanati. Gaps in the model of N. lanati were filled by inspecting the EC assignments found for each other anaerobic fungus, as well as the GEMs of E. coli and S. cerevisiae, using the approach described above and looking for homologous genes in the genome of N. lanati (46, 58). The universal reactions and metabolites from the BiGG Models platform (59) was used to construct the in silico model where possible; if a reaction did not exist in that database it was manually added. The KEGG and MetaCyc databases were used as references to reconstruct the draft metabolic model based on the EC assignments of the metabolic annotation data (56, 60). The curated model for N. lanati was constructed by carefully following established genome-scale metabolic model construction protocols to refine the draft model (26). Specifically, each reaction was inspected to ensure consistency, mass, and charge balance where possible. Model quality was benchmarked by the Memote application (see data set S1 in the iNlan20 GitHub repository at the above-mentioned URL) (61). The curated N. lanati model as well as the entire reconstruction pipeline and all the data used in this work can be found in the model repository at https://github.com/stelmo/iNlan20. An experimentally measured flux of 1.5 mmol/g (dry weight)/h of glucose was used in all simulations. Flux balance analysis was used to simulate the genome-scale metabolic model of N. lanati using the COBRA Toolbox and COBRApy (62, 63). Flux samples (N = 2,000) were generated by sampling from the model and constraining the objective function to be within 90% of the optimum found by FBA. This threshold was set to reflect the assumption that the gut fungi need to maintain a high growth rate to compete with faster growing bacteria in their native microbiome (1). Escher was used to visualize the metabolism (64). Example code that can be used to run the model and computational experiments is supplied as an IPython notebook in the model repository available at the URL mentioned above.

Culturing conditions used for experiments.Standard anaerobic gut fungal culturing techniques were used (11) for all experiments. Briefly, N. lanati was grown at 39°C in sealed Hungate tubes (10-ml liquid volume) or 70-ml serum bottles (40-ml liquid volume) in both undefined complex medium C (MC) (65) and completely defined medium 2 (M2) (66), with 100% CO2 headspace unless otherwise specified. Pressure accumulation in the headspace (67) was used as a proxy for growth, and the fungus was serially passaged after 2 to 3 days of growth. The carbon source was cellobiose (5 g/liter) unless otherwise noted. The cultures were not shaken.

Genome and transcriptome isolation, sequencing, and analysis of N. lanati.N. lanati was isolated from the feces of a sheep located at the Santa Barbara Zoo according to an established protocol (3). Fungal cell pellets for genomic DNA (gDNA) isolation were grown by inoculating 20 ml from a serum bottle of fungi in exponential phase (2 to 3 days of growth given a 10% inoculation volume into the serum bottle) into a 1-liter bottle of medium C, using cellobiose as a carbon source. The serum bottle used to grow the inoculum was treated with chloramphenicol to reduce the risk of contamination. After 4 days of growth, the fungal cell mat was spun down and frozen at −80°C. Four of these frozen samples were subsequently shipped to the Arizona Genome Institute (University of Arizona, Tucson, AZ), where high-quality gDNA was isolated using a modified cetyltrimethylammonium bromide (CTAB) protocol (68). Briefly, these fungal cell mats were ground to a fine powder in a frozen mortar with liquid N2 followed by very gentle extraction in CTAB buffer, which included proteinase K, polyvinylpyrrolidone, molecular weight 40,000 (PVP-40), and 2-mercaptoethanol (Sigma, St. Louis, MO), for 1 h at 50°C. After centrifugation, the supernatant was gently extracted twice with 24:1 chloroform–iso-amyl alcohol. The upper phase was removed, adjusted to one-tenth volume with 3 M potassium acetate, and gently mixed, and the gDNA was precipitated with iso-propanol. Subsequently, the gDNA was collected by centrifugation, washed with 70% ethanol, air dried for 20 min, and dissolved thoroughly in 1× Tris-EDTA (TE) buffer at room temperature. The purified gDNA was shipped to the JGI, where it was sequenced and annotated. Briefly, 10 μg of genomic DNA was sheared to approximately 15 to 20 kb using Megaruptor3 (Diagenode). The sheared DNA was treated with DNA Prep to remove single-stranded ends and then with DNA damage repair mix, followed by end repair, A tailing, and ligation of PacBio overhang adapters using SMRTbell template prep kit 1.0 (Pacific Biosciences). The final library was size selected with BluePippin (Sage Science) at a 10-kb cutoff size and purified with AMPure PB beads. PacBio Sequencing primer v3 was then annealed to the SMRTbell template library, and sequencing polymerase was bound to them using a Sequel binding kit 3.0. The prepared SMRTbell template libraries were then sequenced on a Pacific Biosystem’s Sequel sequencer using 1M v3 SMRT cells, and version 3.0 sequencing chemistry with 10-h movie run times. Subsequently, the main assembly consisted of 62.05× of PacBio read coverage (7,821 bp average read size) and was assembled using MECAT version 1.8; the resulting sequence was polished using ARROW (version 2.2.3). The assembled genome was annotated using the JGI annotation pipeline. The genome is available at https://mycocosm.jgi.doe.gov/Neolan1.

RNA for transcriptome and expression analysis was isolated as previously described (3, 21), in the Biological Nanostructures Lab (University of California Santa Barbara, CA). For the transcriptome, the RNA was harvested from fungal cell pellets grown in serum bottles on a variety of substrates (cellobiose, filter paper, reed canary grass, and corn stover, solids loading 1% [wt/vol], in both medium C and medium 2) to capture as much transcript diversity as possible. For expression analysis, triplicate serum bottles of fungus grown on medium M2, using cellobiose as the sole carbon source, were used. The RNA was isolated and purified using an RNeasy kit (Qiagen, Germantown, MD). The concentration and quality of the RNA were measured on a Qubit (Qubit, New York, NY) and TapeStation 2200 (Agilent, Santa Clara, CA). The RNA used for the transcriptome was pooled in equal parts before sequencing. RNA libraries were made using NEBNext Ultra II directional RNA with mRNA purification beads (NEB, Ipswich, MA); these were subsequently sequenced on a NextSeq 500 (Illumina Inc., San Diego, CA) using high-output 300-cycle settings and 150-base pair paired-end reads (the resultant coverage is 470 and 364 for the transcriptome and expression analysis, respectively). The reads were assembled using Trinity (69). TransDecoder (https://github.com/TransDecoder/TransDecoder/wiki) was used to find the highest likelihood coding regions in the transcriptome. The transcript abundance was estimated using Kallisto (70). The raw assembled transcriptome and filtered output of TransDecoder may be found at https://github.com/stelmo/iNlan20. The mean expression data of transcripts mapped to genes may be found in the GitHub repository of the model available at the URL mentioned above.

A separate transcriptome was sequenced to aid the genome annotation; the same conditions as mentioned previously were used. Stranded cDNA libraries were generated using the Illumina TruSeq stranded mRNA library prep kit. mRNA was purified from 200 ng of total RNA using magnetic beads containing poly(T) oligonucleotides. mRNA was fragmented using divalent cations and high temperature. The fragmented RNA was reversed transcribed using random hexamers and SSII (Invitrogen) followed by second-strand synthesis. The fragmented cDNA was treated with end repair, A tailing, adapter ligation, and 10 cycles of PCR. The prepared libraries were quantified using KAPA Biosystem’s next-generation sequencing library quantitative PCR (qPCR) kit and run on a Roche LightCycler 480 real-time PCR instrument. The libraries were then multiplexed into pools, and sequencing was performed on the Illumina NovaSeq 6000 sequencer using NovaSeq XP v1 reagent kits and S4 flow cell and according to a 2-by-150 indexed run protocol. This transcriptome may be found at https://mycocosm.jgi.doe.gov/Neolan1.

High-performance liquid chromatography, gas chromatography, and liquid chromatography-mass spectrometry measurements.Liquid samples for high-performance liquid chromatography (HPLC) analysis were stored in microcentrifuge tubes at −20°C for batch analysis. Sulfuric acid (0.5 M) was added to the samples (1 in 100 volumes), vortexed, and allowed to mix at room temperature for 5 min. Thereafter, the samples were centrifuged for 5 min at 21,000 × g and filtered using a 0.22-μm syringe filter into HPLC vials. The samples were run on an Agilent 1260 Infinity (Agilent, Santa Clara, CA) using a Bio-Rad HPX-87H column (Bio-Rad, Hercules, CA). Samples were run at two column conditions to effectively separate all the fermentation products (71, 72). Succinate, lactate, cellobiose, glucose, and ethanol were measured at 50°C with a flow rate of 0.5 ml/min and run length of 30 min. Fumarate, formate, and acetate were measured at 25°C with a flow rate of 0.4 ml/min and a run length of 40 min. The mobile phase in both cases was 5 mM sulfuric acid, and the injection volume was 20 μl. Cellobiose, glucose, and ethanol were measured with a refractive index detector, and the other compounds were measured on a variable wavelength detector (λ = 210 nm). Standard curves for each compound were made at 3 concentrations bracketing the range expected in the samples. Gas samples were analyzed on a Thermo Fisher Scientific TRACE gas chromatograph according to a previously established protocol (25). Standard curves of H2 were made daily from supplier (Douglas Fluid & Integration Technology, Prosperity, SC) mixed gas at 1, 2, and 5% (mole basis). Fatty acid composition and extracellular amino acid concentrations were measured on a liquid chromatograph-mass spectrometer (LC-MS) using an Agilent Polaris 3 C18-ether 150-by-3.0 mm (part number [no.] A2021150X030) column. Total run time was 13 min and injection volume was 5 μl with a column temperature of 30°C and flow rate of 0.3 ml/min. The MS was run in negative ionization mode, with the gas temperature at 230°C and a flow rate of 12 liters/min. Standards were made to bracket the expected concentration ranges of the fatty acids.

Biomass composition measurements.The biomass composition of N. lanati was experimentally determined in triplicates for each major component (DNA, RNA, lipid, carbohydrate, and protein). Cultures grown in serum bottles on medium C, with cellobiose as the carbon source, were harvested during exponential phase for each measurement. DNA and RNA were immediately isolated from the wet cell pellet using a previously established CTAB protocol (73). A reference set of triplicate fungal mats was harvested at the same time and lyophilized to estimate the dry mass fraction of DNA and RNA. Lipids and total carbohydrates, which we assumed to be exclusively chitin, were isolated according to established protocols (74). Lipid composition was further refined by mass spectrometry using a fatty acid C18-ether column as described above. Protein extraction from fungal cell pellets was according to a previously developed method (75), and the concentration was measured on a Qubit (Q33327; Thermo Fisher Scientific). The amino acid composition of the protein fraction was assumed to follow the amino acid distribution of the predicted proteome.

Additionally, both the growth and non-growth-associated maintenance (GAM and NGAM, respectively) functions were estimated from experimental data. Briefly, five triplicate sets of Hungate tubes with various concentrations of cellobiose (1, 2, 3, 4, and 5 g/liter initially) in medium M2 were inoculated with 1 ml each (total liquid volume, 10 ml) from a single serum bottle of N. lanati growing at exponential phase (3 days postinoculation) in medium M2 with cellobiose as the carbon source. Pressure accumulation was measured twice daily to calculate the fungal growth rate (67). Liquid and gas samples from each triplicate set were harvested during two time points in exponential phase, 24 h apart. The gas samples were analyzed on a GC to determine the H2 fraction of the gas, and the liquid samples were analyzed by HPLC for organic acid concentration (see “High-performance liquid chromatography, gas chromatography, and liquid chromatography-mass spectrometry measurements” for details). After the last measurement, the fungal cell pellets were harvested by centrifugation, lyophilized, and weighed. The estimated growth rate for each sample was then used to extrapolate the dry cell mass at the respective time points (67). The fluxes of the fermentation products could then be estimated by the molar accumulation of each compound divided by the time between measurements and the difference in cell dry masses between these points. The difference in cell masses was taken because each mature cell lyses and dies; thus, its remaining biomass no longer contributes to metabolism. Finally, these estimated fluxes were used to constrain the model and maximize the ATP yield. The GAM and NGAM were then estimated by finding the line of best fit through the plot of maximum ATP yield predicted by the model and the growth rate associated with the fluxes previously measured (26).

13C metabolic flux analysis for N. lanati.Three serially passaged Hungate tubes using [1,2-13C]glucose as the sole carbon source in medium 2 at an initial concentration of 5 g/liter were used for the labeling experiment. Each Hungate tube was passaged during exponential growth phase, after which, the cell pellet and remaining media were frozen at −20°C for later processing. The medium was analyzed for glucose and fermentation products using the HPLC protocol described above. The pellets were lyophilized, after which GC-MS measurements were used to quantify the isotopic labeling of protein-bound amino acids, glycogen-bound glucose, and RNA-bound ribose as described previously (76). A carbon transition model for flux analysis was constructed using the genome-scale metabolic model of N. lanati as a basis for the flux reactions and biomass equation. Other carbon transition models were used to check the accuracy of the MFA model (77, 78). An ATP demand reaction was added to the carbon transition model and constrained to the maximum ATP flux capable of being generated by the system. This constraint ensured that the maximum feasible flux is routed through the hydrogenosome. INCA was used to perform the flux analysis and sensitivity calculations (79). The carbon transition model, constraints, and the GC-MS data can be found in the GitHub repository of the model available at https://github.com/stelmo/iNlan20.

Model validation experiments.Carbon utilization and vitamin essentiality experiments were conducted to test the predictive accuracy of the model. Carbon utilization was tested by growing N. lanati in medium 2 with each carbon source listed in Table 5 at 5 g/liter initial concentration instead of cellobiose. A carbon substrate was deemed able to support growth if the fungus could be passaged on it for 4 generations and still produce more than 8 lb/in2 gauge of accumulated pressure (no-carbon blanks produce <1 lb/in2 gauge of accumulated pressure). Similarly, the vitamin requirements of N. lanati were tested by individually removing each vitamin in medium 2 (listed in Table 5) and growing the fungus without it for 4 consecutive generations using cellobiose as the carbon source. Fluxes for comparing model predictions to experimental observations were measured similarly to how the fluxes for finding the GAM and NGAM functions were estimated; however, only 5 g/liter cellobiose loading was used. The total equivalent flux of glucose into the cell was calculated by measuring glucose accumulation and cellobiose depletion in the medium. It was assumed that N. lanati imports glucose, and not cellobiose, due to release of beta-glucosidases that decomposed the cellobiose in the medium.

Hydrogenosome staining protocol.The JC-1 dye was purchased from Invitrogen (part no. T3168; Carlsbad, CA, USA), and a standard protocol was used to visualize the presence of electrochemical gradients. Briefly, JC-1 was dissolved in dimethyl sulfoxide (DMSO; 1 mg/ml) and frozen until use. Dye aliquots were thawed and added to cultures of anaerobic gut fungal zoospores using final dye concentrations of 1 μg/ml. Zoospores were incubated with JC-1 for 30 min anaerobically in standard M2 medium at 39°C. After incubation, cultures were filtered onto 3-μm polycarbonate membranes (part no. TSTP02500; Millipore Sigma, Burlington, MA, USA) with a nitrocellulose backing filter (part no. HAWP04700; Millipore Sigma). Cells were counterstained with 4′,6-diamidino-2-phenylindole (DAPI; 2 μg/ml) and mounted on glass slides using an antifade mounting solution composed of 4:1 Citifluor (part no. AF1; Electron Microscopy Sciences, Hatfield PA, USA) to Vectashield (part no. H-1000; Vector Laboratories, Burlingame, CA, USA). Prepared slides were placed on ice and imaged immediately using a Zeiss Axiovert M200 fluorescence microscope (Carl Zeiss AG, Oberkochen, Germany).

Data availability.The full genome, raw sequencing data, assembly, predicted genes, and annotations are available at https://mycocosm.jgi.doe.gov/Neolan1/Neolan1.info.html.

ACKNOWLEDGMENTS

This work was supported by funding from the National Science Foundation (NSF) (MCB-1553721). This work was part of the DOE Joint BioEnergy Institute (http://www.jbei.org) supported by the Office of Biological and Environmental Research of the DOE Office of Science through contract DE-AC02–05CH11231 between Lawrence Berkeley National Laboratory and the DOE. Research was sponsored by the Army Research Office and was accomplished under grant number W911NF-19-1-0010. The work conducted by the U.S. Department of Energy Joint Genome Institute, a DOE Office of Science user facility, is supported by the Office of Science of the U.S. Department of Energy under contract no. DE-AC02-05CH11231. S. E. Wilken received funding support from the Dow Discovery Fellowship.

The views and conclusions contained in this document are those of the authors and should not be interpreted as representing the official policies, either expressed or implied, of the Army Research Office or the U.S. Government. The U.S. Government is authorized to reproduce and distribute reprints for Government purposes notwithstanding any copyright notation herein.

We further acknowledge the use of the Biological Nanostructures Laboratory within the California NanoSystems Institute, supported by the University of California, Santa Barbara and the University of California, Office of the President. We thank Jennifer Smith for help in sequencing the transcriptome of the gut fungus. Use was made of computational facilities purchased with funds from the National Science Foundation (CNS-1725797) and administered by the Center for Scientific Computing (CSC). The CSC is supported by the California NanoSystems Institute and the Materials Research Science and Engineering Center (MRSEC; NSF DMR 1720256) at UC Santa Barbara. We also thank David A. Kudrna and Shanmugam Rajasekar at the Arizona Genome Institute for their help in isolating and purifying high-quality gut fungal gDNA. We thank Maciek R. Antoniewicz, at the University of Michigan, Ann Arbor, MI, for help in generating the measurement data for the 13C tracer experiments. We also acknowledge the use of the UC Santa Barbara NRI-MCDB microscopy facility and thank Weiwei Li in the Keller Lab at UC Santa Barbara for the LC-MS fatty acid and amino acid analysis.

FOOTNOTES

    • Received January 11, 2021.
    • Accepted January 13, 2021.
  • The review history of this article can be read here.

  • Copyright © 2021 Wilken et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

REFERENCES

  1. 1.↵
    1. Gruninger RJ,
    2. Puniya AK,
    3. Callaghan TM,
    4. Edwards JE,
    5. Youssef N,
    6. Dagar SS,
    7. Fliegerova K,
    8. Griffith GW,
    9. Forster R,
    10. Tsang A,
    11. Mcallister T,
    12. Elshahed MS
    . 2014. Anaerobic fungi (phylum Neocallimastigomycota): advances in understanding their taxonomy, life cycle, ecology, role and biotechnological potential. FEMS Microbiol Ecol 90:1–17. doi:10.1111/1574-6941.12383.
    OpenUrlCrossRefPubMed
  2. 2.↵
    1. Hackstein JHP,
    2. Baker SE,
    3. van Hellemond JJ,
    4. Tielens AGM
    . 2019. Hydrogenosomes of anaerobic fungi: an alternative way to adapt to anaerobic environments, p 159–175. Springer, Cham, Switzerland.
  3. 3.↵
    1. Solomon KV,
    2. Haitjema CH,
    3. Henske JK,
    4. Gilmore SP,
    5. Borges-Rivera D,
    6. Lipzen A,
    7. Brewer HM,
    8. Purvine SO,
    9. Wright AT,
    10. Theodorou MK,
    11. Grigoriev IV,
    12. Regev A,
    13. Thompson DA,
    14. O'Malley MA
    . 2016. Early-branching gut fungi possess a large, comprehensive array of biomass-degrading enzymes. Science 351:1192–1196. doi:10.1126/science.aad1431.
    OpenUrlAbstract/FREE Full Text
  4. 4.↵
    1. Haitjema CH,
    2. Gilmore SP,
    3. Henske JK,
    4. Solomon KV,
    5. de Groot R,
    6. Kuo A,
    7. Mondo SJ,
    8. Salamov AA,
    9. LaButti K,
    10. Zhao Z,
    11. Chiniquy J,
    12. Barry K,
    13. Brewer HM,
    14. Purvine SO,
    15. Wright AT,
    16. Hainaut M,
    17. Boxma B,
    18. van Alen T,
    19. Hackstein JHP,
    20. Henrissat B,
    21. Baker SE,
    22. Grigoriev IV,
    23. O'Malley MA
    . 2017. A parts list for fungal cellulosomes revealed by comparative genomics. Nat Microbiol 2:17087. doi:10.1038/nmicrobiol.2017.87.
    OpenUrlCrossRef
  5. 5.↵
    1. Resch MG,
    2. Donohoe BS,
    3. Baker JO,
    4. Decker SR,
    5. Bayer EA,
    6. Beckham GT,
    7. Himmel ME
    . 2013. Fungal cellulases and complexed cellulosomal enzymes exhibit synergistic mechanisms in cellulose deconstruction. Energy Environ Sci 6:1858–1867. doi:10.1039/c3ee00019b.
    OpenUrlCrossRef
  6. 6.↵
    1. Seppälä S,
    2. Wilken SE,
    3. Knop D,
    4. Solomon KV,
    5. O'Malley MA
    . 2017. The importance of sourcing enzymes from non-conventional fungi for metabolic engineering and biomass breakdown. Metab Eng 44:45–59. doi:10.1016/j.ymben.2017.09.008.
    OpenUrlCrossRef
  7. 7.↵
    1. Youssef NH,
    2. Couger MB,
    3. Struchtemeyer CG,
    4. Liggenstoffer AS,
    5. Prade RA,
    6. Najar FZ,
    7. Atiyeh HK,
    8. Wilkins MR,
    9. Elshahed MS
    . 2013. The genome of the anaerobic fungus Orpinomyces sp. strain C1A reveals the unique evolutionary history of a remarkable plant biomass degrader. Appl Environ Microbiol 79:4620–4634. doi:10.1128/AEM.00821-13.
    OpenUrlAbstract/FREE Full Text
  8. 8.↵
    1. Sarkar N,
    2. Ghosh SK,
    3. Bannerjee S,
    4. Aikat K
    . 2012. Bioethanol production from agricultural wastes: an overview. Renew Energy 37:19–27. doi:10.1016/j.renene.2011.06.045.
    OpenUrlCrossRefWeb of Science
  9. 9.↵
    1. Himmel ME,
    2. Ding S-Y,
    3. Johnson DK,
    4. Adney WS,
    5. Nimlos MR,
    6. Brady JW,
    7. Foust TD
    . 2007. Biomass recalcitrance: engineering plants and enzymes for biofuels production. Science 315:804–807. doi:10.1126/science.1137016.
    OpenUrlAbstract/FREE Full Text
  10. 10.↵
    1. Liao JC,
    2. Mi L,
    3. Pontrelli S,
    4. Luo S
    . 2016. Fuelling the future: microbial engineering for the production of sustainable biofuels. Nat Rev Microbiol 14:288–304. doi:10.1038/nrmicro.2016.32.
    OpenUrlCrossRefPubMed
  11. 11.↵
    1. Haitjema CH,
    2. Solomon KV,
    3. Henske JK,
    4. Theodorou MK,
    5. O'Malley MA
    . 2014. Anaerobic gut fungi: advances in isolation, culture, and cellulolytic enzyme discovery for biofuel production. Biotechnol Bioeng 111:1471–1482. doi:10.1002/bit.25264.
    OpenUrlCrossRefPubMed
  12. 12.↵
    1. Wilken SE,
    2. Seppälä S,
    3. Lankiewicz TS,
    4. Saxena M,
    5. Henske JK,
    6. Salamov AA,
    7. Grigoriev IV,
    8. O’Malley MA
    . 2020. Genomic and proteomic biases inform metabolic engineering strategies for anaerobic fungi. Metab Eng Commun 10:e00107. doi:10.1016/j.mec.2019.e00107.
    OpenUrlCrossRef
  13. 13.↵
    1. Henske JK,
    2. Gilmore SP,
    3. Haitjema CH,
    4. Solomon KV,
    5. O'Malley MA
    . 2018. Biomass-degrading enzymes are catabolite repressed in anaerobic gut fungi. AIChE J 64:4263–4270. doi:10.1002/aic.16395.
    OpenUrlCrossRef
  14. 14.↵
    1. Boxma B,
    2. Voncken F,
    3. Jannink S,
    4. Van Alen T,
    5. Akhmanova A,
    6. Van Weelden SWH,
    7. Van Hellemond JJ,
    8. Ricard G,
    9. Huynen M,
    10. Tielens AGM,
    11. Hackstein JHP
    . 2004. The anaerobic chytridiomycete fungus Piromyces sp. E2 produces ethanol via pyruvate:formate lyase and an alcohol dehydrogenase E. Mol Microbiol 51:1389–1399. doi:10.1046/j.1365-2958.2003.03912.x.
    OpenUrlCrossRefPubMedWeb of Science
  15. 15.↵
    1. Akhmanova A,
    2. Voncken FGJ,
    3. Hosea KM,
    4. Harhangi H,
    5. Keltjens JT,
    6. Op den Camp HJM,
    7. Vogels GD,
    8. Hackstein JHP
    . 1999. A hydrogenosome with pyruvate formate-lyase: anaerobic chytrid fungi use an alternative route for pyruvate catabolism. Mol Microbiol 32:1103–1114. doi:10.1046/j.1365-2958.1999.01434.x.
    OpenUrlCrossRefPubMedWeb of Science
  16. 16.↵
    1. Marvin-Sikkema FD,
    2. Driessen AJM,
    3. Gottschal JC,
    4. Prins RA
    . 1994. Metabolic energy generation in hydrogenosomes of the anaerobic fungus Neocallimastix: evidence for a functional relationship with mitochondria. Mycol Res 98:205–212. doi:10.1016/S0953-7562(09)80187-1.
    OpenUrlCrossRef
  17. 17.↵
    1. Simeonidis E,
    2. Price ND
    . 2015. Genome-scale modeling for metabolic engineering. J Ind Microbiol Biotechnol 42:327–338. doi:10.1007/s10295-014-1576-3.
    OpenUrlCrossRef
  18. 18.↵
    1. Blazeck J,
    2. Alper H
    . 2010. Systems metabolic engineering: genome-scale models and beyond. Biotechnol J 5:647–659. doi:10.1002/biot.200900247.
    OpenUrlCrossRefPubMedWeb of Science
  19. 19.↵
    1. Aung HW,
    2. Henry SA,
    3. Walker LP
    . 2013. Revising the representation of fatty acid, glycerolipid, and glycerophospholipid metabolism in the consensus model of yeast metabolism. Ind Biotechnol (New Rochelle N Y) 9:215–228. doi:10.1089/ind.2013.0013.
    OpenUrlCrossRef
  20. 20.↵
    1. Wang TY,
    2. Chen HL,
    3. Lu MYJ,
    4. Chen YC,
    5. Sung HM,
    6. Mao CT,
    7. Cho HY,
    8. Ke HM,
    9. Hwa TY,
    10. Ruan SK,
    11. Hung KY,
    12. Chen CK,
    13. Li JY,
    14. Wu YC,
    15. Chen YH,
    16. Chou SP,
    17. Tsai YW,
    18. Chu TC,
    19. Shih CCA,
    20. Li WH,
    21. Shih MC
    . 2011. Functional characterization of cellulases identified from the cow rumen fungus Neocallimastix patriciarum W5 by transcriptomic and secretomic analyses. Biotechnol Biofuels 4:24. doi:10.1186/1754-6834-4-24.
    OpenUrlCrossRefPubMed
  21. 21.↵
    1. Henske JK,
    2. Wilken SE,
    3. Solomon KV,
    4. Smallwood CR,
    5. Shutthanandan V,
    6. Evans JE,
    7. Theodorou MK,
    8. O'Malley MA
    . 2018. Metabolic characterization of anaerobic fungi provides a path forward for bioprocessing of crude lignocellulose. Biotechnol Bioeng 115:874–884. doi:10.1002/bit.26515.
    OpenUrlCrossRef
  22. 22.↵
    1. Wilken SE,
    2. Leggieri PA,
    3. Kerdman‐Andrade C,
    4. Reilly M,
    5. Theodorou MK,
    6. O'Malley MA
    . 2020. An Arduino based automatic pressure evaluation system to quantify growth of non‐model anaerobes in culture. AIChE J 66:e16540. doi:10.1002/aic.16540.
    OpenUrlCrossRef
  23. 23.↵
    1. Nielsen BB,
    2. Zhu WY,
    3. Dhanoa MS,
    4. Trinci APJ,
    5. Theodorou MK
    . 2002. Comparison of the growth kinetics of anaerobic gut fungi on wheat straw in batch culture. Anaerobe 8:216–222. doi:10.1006/anae.2002.0432.
    OpenUrlCrossRef
  24. 24.↵
    1. Ranganathan A,
    2. Smith OP,
    3. Youssef NH,
    4. Struchtemeyer CG,
    5. Atiyeh HK,
    6. Elshahed MS
    . 2017. Utilizing anaerobic fungi for two-stage sugar extraction and biofuel production from lignocellulosic biomass. Front Microbiol 8:635. doi:10.3389/fmicb.2017.00635.
    OpenUrlCrossRef
  25. 25.↵
    1. Gilmore SP,
    2. Lankiewicz TS,
    3. Wilken SE,
    4. Brown JL,
    5. Sexton JA,
    6. Henske JK,
    7. Theodorou MK,
    8. Valentine DL,
    9. O’Malley MA
    . 2019. Top-down enrichment guides in formation of synthetic microbial consortia for biomass degradation. ACS Synth Biol 8:2174–2185. doi:10.1021/acssynbio.9b00271.
    OpenUrlCrossRef
  26. 26.↵
    1. Thiele I,
    2. Palsson B
    . 2010. A protocol for generating a high-quality genome-scale metabolic reconstruction. Nat Protoc 5:93–121. doi:10.1038/nprot.2009.203.
    OpenUrlCrossRefPubMedWeb of Science
  27. 27.↵
    1. Mountfort D,
    2. Orpin CG
    (ed). 1994. Anaerobic fungi: biology, ecology, and function, first ed. CRC Press, Boca Raton, FL.
  28. 28.↵
    1. Flint HJ,
    2. Bayer EA,
    3. Rincon MT,
    4. Lamed R,
    5. White BA
    . 2008. Polysaccharide utilization by gut bacteria: potential for new insights from genomic analysis. Nat Rev Microbiol 6:121–131. doi:10.1038/nrmicro1817.
    OpenUrlCrossRefPubMedWeb of Science
  29. 29.↵
    1. Martínez I,
    2. Zhu J,
    3. Lin H,
    4. Bennett GN,
    5. San KY
    . 2008. Replacing Escherichia coli NAD-dependent glyceraldehyde 3-phosphate dehydrogenase (GAPDH) with a NADP-dependent enzyme from Clostridium acetobutylicum facilitates NADPH dependent pathways. Metab Eng 10:352–359. doi:10.1016/j.ymben.2008.09.001.
    OpenUrlCrossRefPubMed
  30. 30.↵
    1. Boyd DA,
    2. Cvitkovitch DG,
    3. Hamilton IANR
    . 1995. Sequence, expression, and function of the gene for the phosphate dehydrogenase of Streptococcus mutans. J Bacteriol 177:2622–2627. doi:10.1128/jb.177.10.2622-2627.1995.
    OpenUrlAbstract/FREE Full Text
  31. 31.↵
    1. Yarlett N,
    2. Orpin CG,
    3. Munn EA,
    4. Yarlett NC,
    5. Greenwood CA
    . 1986. Hydrogenosomes in the rumen fungus Neocallimastix patriciarum. Biochem J 236:729–739. doi:10.1042/bj2360729.
    OpenUrlAbstract/FREE Full Text
  32. 32.↵
    1. Muller M,
    2. Mentel M,
    3. van Hellemond JJ,
    4. Henze K,
    5. Woehle C,
    6. Gould SB,
    7. Yu R-Y,
    8. van der Giezen M,
    9. Tielens AGM,
    10. Martin WF
    . 2012. Biochemistry and evolution of anaerobic energy metabolism in eukaryotes. Microbiol Mol Biol Rev 76:444–495. doi:10.1128/MMBR.05024-11.
    OpenUrlAbstract/FREE Full Text
  33. 33.↵
    1. Junge W,
    2. Nelson N
    . 2015. ATP synthase. Annu Rev Biochem 84:631–657. doi:10.1146/annurev-biochem-060614-034124.
    OpenUrlCrossRefPubMed
  34. 34.↵
    1. Seppälä S,
    2. Solomon KV,
    3. Gilmore SP,
    4. Henske JK,
    5. O'Malley MA
    . 2016. Mapping the membrane proteome of anaerobic gut fungi identifies a wealth of carbohydrate binding proteins and transporters. Microb Cell Fact 15:212. doi:10.1186/s12934-016-0611-7.
    OpenUrlCrossRef
  35. 35.↵
    1. Hrdy I,
    2. Hirt RP,
    3. Dolezal P,
    4. Bardonová L,
    5. Foster PG,
    6. Tachezy J,
    7. Embley TM
    . 2004. Trichomonas hydrogenosomes contain the NADH dehydrogenase module of mitochondrial complex I. Nature 432:618–622. doi:10.1038/nature03149.
    OpenUrlCrossRefPubMedWeb of Science
  36. 36.↵
    1. Schneider RE,
    2. Brown MT,
    3. Shiflett AM,
    4. Dyall SD,
    5. Hayes RD,
    6. Xie Y,
    7. Loo JA,
    8. Johnson PJ
    . 2011. The Trichomonas vaginalis hydrogenosome proteome is highly reduced relative to mitochondria, yet complex compared with mitosomes. Int J Parasitol 41:1421–1434. doi:10.1016/j.ijpara.2011.10.001.
    OpenUrlCrossRefPubMed
  37. 37.↵
    1. Boxma B,
    2. Ricard G,
    3. Van Hoek AHAM,
    4. Severing E,
    5. Moon-Van Der Staay SY,
    6. Van Der Staay GWM,
    7. Van Alen TA,
    8. De Graaf RM,
    9. Cremers G,
    10. Kwantes M,
    11. McEwan NR,
    12. Newbold CJ,
    13. Jouany JP,
    14. Michalowski T,
    15. Pristas P,
    16. Huynen MA,
    17. Hackstein JHP
    . 2007. The [FeFe] hydrogenase of Nyctotherus ovalis has a chimeric origin. BMC Evol Biol 7:230. doi:10.1186/1471-2148-7-230.
    OpenUrlCrossRefPubMed
  38. 38.↵
    1. Schut GJ,
    2. Adams MWW
    . 2009. The iron-hydrogenase of Thermotoga maritima utilizes ferredoxin and NADH synergistically: a new perspective on anaerobic hydrogen production. J Bacteriol 191:4451–4457. doi:10.1128/JB.01582-08.
    OpenUrlAbstract/FREE Full Text
  39. 39.↵
    1. Orpin CG,
    2. Greenwood Y
    . 1986. The role of haems and related compounds in the nutrition and zoosporogenesis of the rumen chytridiomycete Neocallimastix frontalis H8. J Gen Microbiol 132:2179–2185. doi:10.1099/00221287-132-8-2179.
    OpenUrlCrossRefWeb of Science
  40. 40.↵
    1. Marvin-Sikkema FD,
    2. Pedro Gomes TM,
    3. Grivet JP,
    4. Gottschal JC,
    5. Prins RA
    . 1993. Characterization of hydrogenosomes and their role in glucose metabolism of Neocallimastix sp. L2. Arch Microbiol 160:388–396. doi:10.1007/BF00252226.
    OpenUrlCrossRefPubMedWeb of Science
  41. 41.↵
    1. Stairs CW,
    2. Roger AJ,
    3. Hampl V
    . 2011. Eukaryotic pyruvate formate lyase and its activating enzyme were acquired laterally from a firmicute. Mol Biol Evol 28:2087–2099. doi:10.1093/molbev/msr032.
    OpenUrlCrossRefPubMedWeb of Science
  42. 42.↵
    1. Müller V,
    2. Chowdhury NP,
    3. Basen M
    . 2018. Electron bifurcation: a long-hidden energy-coupling mechanism. Annu Rev Microbiol 72:331–353. doi:10.1146/annurev-micro-090816-093440.
    OpenUrlCrossRefPubMed
  43. 43.↵
    1. Boxma B,
    2. de Graaf RM,
    3. van der Staay GWM,
    4. van Alen TA,
    5. Ricard G,
    6. Gabaldón T,
    7. van Hoek AHAM,
    8. Moon-van der Staay SY,
    9. Koopman WJH,
    10. van Hellemond JJ,
    11. Tielens AGM,
    12. Friedrich T,
    13. Veenhuis M,
    14. Huynen MA,
    15. Hackstein JHP
    . 2005. An anaerobic mitochondrion that produces hydrogen. Nature 434:74–79. doi:10.1038/nature03343.
    OpenUrlCrossRefPubMedWeb of Science
  44. 44.↵
    1. Zorova LD,
    2. Popkov VA,
    3. Plotnikov EY,
    4. Silachev DN,
    5. Pevzner IB,
    6. Jankauskas SS,
    7. Babenko VA,
    8. Zorov SD,
    9. Balakireva AV,
    10. Juhaszova M,
    11. Sollott SJ,
    12. Zorov DB
    . 2018. Mitochondrial membrane potential. Anal Biochem 552:50–59. doi:10.1016/j.ab.2017.07.009.
    OpenUrlCrossRefPubMed
  45. 45.↵
    1. Chubukov V,
    2. Mukhopadhyay A,
    3. Petzold CJ,
    4. Keasling JD,
    5. Martín HG
    . 2018. Synthetic and systems biology for microbial production of commodity chemicals. Npj Syst Biol Appl 2:1609. doi:10.1038/npjsba.2016.9.
    OpenUrlCrossRef
  46. 46.↵
    1. Monk JM,
    2. Lloyd CJ,
    3. Brunk E,
    4. Mih N,
    5. Sastry A,
    6. King Z,
    7. Takeuchi R,
    8. Nomura W,
    9. Zhang Z,
    10. Mori H,
    11. Feist AM,
    12. Palsson BO
    . 2017. iML1515, a knowledgebase that computes Escherichia coli traits. Nat Biotechnol 35:904–908. doi:10.1038/nbt.3956.
    OpenUrlCrossRef
  47. 47.↵
    1. Lu H,
    2. Li F,
    3. Sánchez BJ,
    4. Zhu Z,
    5. Li G,
    6. Domenzain I,
    7. Marcišauskas S,
    8. Anton PM,
    9. Lappa D,
    10. Lieven C,
    11. Beber ME,
    12. Sonnenschein N,
    13. Kerkhoven EJ,
    14. Nielsen J
    . 2019. A consensus S. cerevisiae metabolic model Yeast8 and its ecosystem for comprehensively probing cellular metabolism. Nat Commun 10:3586. doi:10.1038/s41467-019-11581-3.
    OpenUrlCrossRef
  48. 48.↵
    1. Grigoriev IV,
    2. Nikitin R,
    3. Haridas S,
    4. Kuo A,
    5. Ohm R,
    6. Otillar R,
    7. Riley R,
    8. Salamov A,
    9. Zhao X,
    10. Korzeniewski F,
    11. Smirnova T,
    12. Nordberg H,
    13. Dubchak I,
    14. Shabalov I
    . 2014. MycoCosm portal: gearing up for 1000 fungal genomes. Nucleic Acids Res 42:D699–D704. doi:10.1093/nar/gkt1183.
    OpenUrlCrossRefPubMedWeb of Science
  49. 49.↵
    1. Li Y,
    2. Li Y,
    3. Jin W,
    4. Sharpton TJ,
    5. Mackie RI,
    6. Cann I,
    7. Cheng Y,
    8. Zhu W
    . 2019. Combined genomic, transcriptomic, proteomic, and physiological characterization of the growth of Pecoramyces sp. F1 in monoculture and co-culture with a syntrophic methanogen. Front Microbiol 10:435. doi:10.3389/fmicb.2019.00435.
    OpenUrlCrossRef
  50. 50.↵
    1. Martinez-Fernandez G,
    2. Denman SE,
    3. Yang C,
    4. Cheung J,
    5. Mitsumori M,
    6. McSweeney CS
    . 2016. Methane inhibition alters the microbial community, hydrogen flow, and fermentation response in the rumen of cattle. Front Microbiol 7:1122. doi:10.3389/fmicb.2016.01122.
    OpenUrlCrossRef
  51. 51.↵
    1. Jeske L,
    2. Placzek S,
    3. Schomburg I,
    4. Chang A,
    5. Schomburg D
    . 2019. BRENDA in 2019: a European ELIXIR core data resource. Nucleic Acids Res 47:D542–D549. doi:10.1093/nar/gky1048.
    OpenUrlCrossRef
  52. 52.↵
    1. Hanafy RA,
    2. Elshahed MS,
    3. Liggenstoffer AS,
    4. Griffith GW,
    5. Youssef NH
    . 2017. Pecoramyces ruminantium, gen. nov., sp. nov., an anaerobic gut fungus from the feces of cattle and sheep. Mycologia 109:231–243. doi:10.1080/00275514.2017.1317190.
    OpenUrlCrossRef
  53. 53.↵
    1. Camacho C,
    2. Coulouris G,
    3. Avagyan V,
    4. Ma N,
    5. Papadopoulos J,
    6. Bealer K,
    7. Madden TL
    . 2009. BLAST+: architecture and applications. BMC Bioinformatics 10:421. doi:10.1186/1471-2105-10-421.
    OpenUrlCrossRefPubMed
  54. 54.↵
    UniProt Consortium. 2019. UniProt: a worldwide hub of protein knowledge. Nucleic Acids Res 47:D506–D515. doi:10.1093/nar/gky1049.
    OpenUrlCrossRefPubMed
  55. 55.↵
    1. Almagro Armenteros JJ,
    2. Sønderby CK,
    3. Sønderby SK,
    4. Nielsen H,
    5. Winther O
    . 2017. DeepLoc: prediction of protein subcellular localization using deep learning. Bioinformatics 33:3387–3395. doi:10.1093/bioinformatics/btx431.
    OpenUrlCrossRefPubMed
  56. 56.↵
    1. Caspi R,
    2. Billington R,
    3. Fulcher CA,
    4. Keseler IM,
    5. Kothari A,
    6. Krummenacker M,
    7. Latendresse M,
    8. Midford PE,
    9. Ong Q,
    10. Ong WK,
    11. Paley S,
    12. Subhraveti P,
    13. Karp PD
    . 2018. The MetaCyc database of metabolic pathways and enzymes. Nucleic Acids Res 46:D633–D639. doi:10.1093/nar/gkx935.
    OpenUrlCrossRefPubMed
  57. 57.↵
    1. Flamholz A,
    2. Noor E,
    3. Bar-Even A,
    4. Milo R
    . 2012. eQuilibrator - the biochemical thermodynamics calculator. Nucleic Acids Res 40:D770–D775. doi:10.1093/nar/gkr874.
    OpenUrlCrossRefPubMedWeb of Science
  58. 58.↵
    1. Mo ML,
    2. Palsson B,
    3. Herrgård MJ
    . 2009. Connecting extracellular metabolomic measurements to intracellular flux states in yeast. BMC Syst Biol 3:37. doi:10.1186/1752-0509-3-37.
    OpenUrlCrossRefPubMed
  59. 59.↵
    1. King ZA,
    2. Lu J,
    3. Dräger A,
    4. Miller P,
    5. Federowicz S,
    6. Lerman JA,
    7. Ebrahim A,
    8. Palsson BO,
    9. Lewis NE
    . 2016. BiGG Models: a platform for integrating, standardizing and sharing genome-scale models. Nucleic Acids Res 44:D515–D522. doi:10.1093/nar/gkv1049.
    OpenUrlCrossRefPubMed
  60. 60.↵
    1. Kanehisa M,
    2. Sato Y,
    3. Kawashima M,
    4. Furumichi M,
    5. Tanabe M
    . 2016. KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res 44:D457–D462. doi:10.1093/nar/gkv1070.
    OpenUrlCrossRefPubMed
  61. 61.↵
    1. Lieven C,
    2. Beber ME,
    3. Olivier BG,
    4. Bergmann FT,
    5. Ataman M,
    6. Babaei P,
    7. Bartell JA,
    8. Blank LM,
    9. Chauhan S,
    10. Correia K,
    11. Diener C,
    12. Dräger A,
    13. Ebert BE,
    14. Edirisinghe JN,
    15. Faria JP,
    16. Feist AM,
    17. Fengos G,
    18. Fleming RMT,
    19. García-Jiménez B,
    20. Hatzimanikatis V,
    21. van Helvoirt W,
    22. Henry CS,
    23. Hermjakob H,
    24. Herrgård MJ,
    25. Kaafarani A,
    26. Kim HU,
    27. King Z,
    28. Klamt S,
    29. Klipp E,
    30. Koehorst JJ,
    31. König M,
    32. Lakshmanan M,
    33. Lee DY,
    34. Lee SY,
    35. Lee S,
    36. Lewis NE,
    37. Liu F,
    38. Ma H,
    39. Machado D,
    40. Mahadevan R,
    41. Maia P,
    42. Mardinoglu A,
    43. Medlock GL,
    44. Monk JM,
    45. Nielsen J,
    46. Nielsen LK,
    47. Nogales J,
    48. Nookaew I,
    49. Palsson BO,
    50. Papin JA
    , et al. 2020. MEMOTE for standardized genome-scale metabolic model testing. Nat Biotechnol 38:272–276. doi:10.1038/s41587-020-0446-y.
    OpenUrlCrossRef
  62. 62.↵
    1. Heirendt L,
    2. Arreckx S,
    3. Pfau T,
    4. Mendoza SN,
    5. Richelle A,
    6. Heinken A,
    7. Haraldsdóttir HS,
    8. Wachowiak J,
    9. Keating SM,
    10. Vlasov V,
    11. Magnusdóttir S,
    12. Ng CY,
    13. Preciat G,
    14. Žagare A,
    15. Chan SHJ,
    16. Aurich MK,
    17. Clancy CM,
    18. Modamio J,
    19. Sauls JT,
    20. Noronha A,
    21. Bordbar A,
    22. Cousins B,
    23. El Assal DC,
    24. Valcarcel LV,
    25. Apaolaza I,
    26. Ghaderi S,
    27. Ahookhosh M,
    28. Ben Guebila M,
    29. Kostromins A,
    30. Sompairac N,
    31. Le HM,
    32. Ma D,
    33. Sun Y,
    34. Wang L,
    35. Yurkovich JT,
    36. Oliveira MAP,
    37. Vuong PT,
    38. El Assal LP,
    39. Kuperstein I,
    40. Zinovyev A,
    41. Hinton HS,
    42. Bryant WA,
    43. Aragón Artacho FJ,
    44. Planes FJ,
    45. Stalidzans E,
    46. Maass A,
    47. Vempala S,
    48. Hucka M,
    49. Saunders MA,
    50. Maranas CD
    , et al. 2019. Creation and analysis of biochemical constraint-based models using the COBRA Toolbox v.3.0. Nat Protoc 14:639–702. doi:10.1038/s41596-018-0098-2.
    OpenUrlCrossRef
  63. 63.↵
    1. Ebrahim A,
    2. Lerman JA,
    3. Palsson BO,
    4. Hyduke DR
    . 2013. COBRApy: COnstraints-Based Reconstruction and Analysis for Python. BMC Syst Biol 7:74. doi:10.1186/1752-0509-7-74.
    OpenUrlCrossRefPubMed
  64. 64.↵
    1. Rowe E,
    2. Palsson BO,
    3. King ZA
    . 2018. Escher-FBA: a web application for interactive flux balance analysis. BMC Syst Biol 12:84. doi:10.1186/s12918-018-0607-5.
    OpenUrlCrossRef
  65. 65.↵
    1. Davies DR,
    2. Theodorou MK,
    3. Lawrence MIG,
    4. Trinci APJ
    . 1993. Distribution of anaerobic fungi in the digestive tract of cattle and their survival in faeces. J Gen Microbiol 139:1395–1400. doi:10.1099/00221287-139-6-1395.
    OpenUrlCrossRefPubMedWeb of Science
  66. 66.↵
    1. Teunissen MJ,
    2. Op den Camp HJM,
    3. Orpin CG,
    4. Huis In 't Veld JHJ,
    5. Vogels GD
    . 1991. Comparison of growth characteristics of anaerobic fungi isolated from ruminant and non-ruminant herbivores during cultivation in a defined medium. J Gen Microbiol 137:1401–1408. doi:10.1099/00221287-137-6-1401.
    OpenUrlCrossRefPubMedWeb of Science
  67. 67.↵
    1. Theodorou MK,
    2. Davies DR,
    3. Nielsen BB,
    4. Lawrence MIG,
    5. Trinci APJ
    . 1995. Determination of growth of anaerobic fungi on soluble and cellulosic substrates using a pressure transducer. Microbiology 141:671–678. doi:10.1099/13500872-141-3-671.
    OpenUrlCrossRefWeb of Science
  68. 68.↵
    1. Doyle JJ,
    2. Doyle JL
    . 1987. A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochem Bull 19:11–15.
    OpenUrlCrossRef
  69. 69.↵
    1. Grabherr MG,
    2. Haas BJ,
    3. Yassour M,
    4. Levin JZ,
    5. Thompson DA,
    6. Amit I,
    7. Adiconis X,
    8. Fan L,
    9. Raychowdhury R,
    10. Zeng Q,
    11. Chen Z,
    12. Mauceli E,
    13. Hacohen N,
    14. Gnirke A,
    15. Rhind N,
    16. Di Palma F,
    17. Birren BW,
    18. Nusbaum C,
    19. Lindblad-Toh K,
    20. Friedman N,
    21. Regev A
    . 2011. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol 29:644–652. doi:10.1038/nbt.1883.
    OpenUrlCrossRefPubMed
  70. 70.↵
    1. Bray NL,
    2. Pimentel H,
    3. Melsted P,
    4. Pachter L
    . 2016. Near-optimal probabilistic RNA-seq quantification. Nat Biotechnol 34:525–527. doi:10.1038/nbt.3519.
    OpenUrlCrossRefPubMed
  71. 71.↵
    1. Chinnici F,
    2. Spinabelli U,
    3. Riponi C,
    4. Amati A
    . 2005. Optimization of the determination of organic acids and sugars in fruit juices by ion-exclusion liquid chromatography. J Food Compos Anal 18:121–130. doi:10.1016/j.jfca.2004.01.005.
    OpenUrlCrossRefWeb of Science
  72. 72.↵
    1. Qiu J,
    2. Jin X
    . 2002. Development and optimization of organic acid analysis in tobacco with ion chromatography and suppressed conductivity detection. J Chromatogr A 950:81–88. doi:10.1016/s0021-9673(02)00034-1.
    OpenUrlCrossRefPubMedWeb of Science
  73. 73.↵
    1. Lankiewicz TS,
    2. Cottrell MT,
    3. Kirchman DL
    . 2016. Growth rates and rRNA content of four marine bacteria in pure cultures and in the Delaware estuary. ISME J 10:823–832. doi:10.1038/ismej.2015.156.
    OpenUrlCrossRefPubMed
  74. 74.↵
    1. Beck AE,
    2. Hunt KA,
    3. Carlson RP
    . 2018. Measuring cellular biomass composition for computational biology applications. Processes 6:38. doi:10.3390/pr6050038.
    OpenUrlCrossRef
  75. 75.↵
    1. Soh L,
    2. Montazeri M,
    3. Haznedaroglu BZ,
    4. Kelly C,
    5. Peccia J,
    6. Eckelman MJ,
    7. Zimmerman JB
    . 2014. Evaluating microalgal integrated biorefinery schemes: empirical controlled growth studies and life cycle assessment. Bioresour Technol 151:19–27. doi:10.1016/j.biortech.2013.10.012.
    OpenUrlCrossRef
  76. 76.↵
    1. Long CP,
    2. Antoniewicz MR
    . 2019. High-resolution 13C metabolic flux analysis. Nat Protoc 14:2856–2877. doi:10.1038/s41596-019-0204-0.
    OpenUrlCrossRef
  77. 77.↵
    1. Crown SB,
    2. Long CP,
    3. Antoniewicz MR
    . 2016. Optimal tracers for parallel labeling experiments and 13C metabolic flux analysis: a new precision and synergy scoring system. Metab Eng 38:10–18. doi:10.1016/j.ymben.2016.06.001.
    OpenUrlCrossRefPubMed
  78. 78.↵
    1. Liu N,
    2. Qiao K,
    3. Stephanopoulos G
    . 2016. 13C metabolic flux analysis of acetate conversion to lipids by Yarrowia lipolytica. Metab Eng 38:86–97. doi:10.1016/j.ymben.2016.06.006.
    OpenUrlCrossRef
  79. 79.↵
    1. Young JD
    . 2014. INCA: a computational platform for isotopically non-stationary metabolic flux analysis. Bioinformatics 30:1333–1335. doi:10.1093/bioinformatics/btu015.
    OpenUrlCrossRefPubMedWeb of Science
PreviousNext
Back to top
Download PDF
Citation Tools
Experimentally Validated Reconstruction and Analysis of a Genome-Scale Metabolic Model of an Anaerobic Neocallimastigomycota Fungus
St. Elmo Wilken, Jonathan M. Monk, Patrick A. Leggieri, Christopher E. Lawson, Thomas S. Lankiewicz, Susanna Seppälä, Chris G. Daum, Jerry Jenkins, Anna M. Lipzen, Stephen J. Mondo, Kerrie W. Barry, Igor V. Grigoriev, John K. Henske, Michael K. Theodorou, Bernhard O. Palsson, Linda R. Petzold, Michelle A. O’Malley
mSystems Feb 2021, 6 (1) e00002-21; DOI: 10.1128/mSystems.00002-21

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print
Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this mSystems article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Experimentally Validated Reconstruction and Analysis of a Genome-Scale Metabolic Model of an Anaerobic Neocallimastigomycota Fungus
(Your Name) has forwarded a page to you from mSystems
(Your Name) thought you would be interested in this article in mSystems.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Experimentally Validated Reconstruction and Analysis of a Genome-Scale Metabolic Model of an Anaerobic Neocallimastigomycota Fungus
St. Elmo Wilken, Jonathan M. Monk, Patrick A. Leggieri, Christopher E. Lawson, Thomas S. Lankiewicz, Susanna Seppälä, Chris G. Daum, Jerry Jenkins, Anna M. Lipzen, Stephen J. Mondo, Kerrie W. Barry, Igor V. Grigoriev, John K. Henske, Michael K. Theodorou, Bernhard O. Palsson, Linda R. Petzold, Michelle A. O’Malley
mSystems Feb 2021, 6 (1) e00002-21; DOI: 10.1128/mSystems.00002-21
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS
    • DISCUSSION
    • MATERIALS AND METHODS
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • Review History
  • PDF

KEYWORDS

genome-scale metabolic model
13C metabolic flux analysis
nonmodel fungus
Neocallimastigomycota
flux balance analysis
Neocallimastix lanati
anaerobes
anaerobic fungi

Related Articles

Cited By...

About

  • About mSystems
  • Author Videos
  • Board of Editors
  • Policies
  • Overleaf Pilot
  • For Reviewers
  • For the Media
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Author Warranty
  • Types of Articles
  • Getting Started
  • Ethics
  • Contact Us

Follow #mSystemsJ

@ASMicrobiology

       

 

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

Copyright © 2021 American Society for Microbiology | Privacy Policy | Website feedback

Online ISSN: 2379-5077