Genetic mapping of microbial and host traits reveals production of immunomodulatory lipids by Akkermansia muciniphila in the murine gut.

The molecular bases of how host genetic variation impacts the gut microbiome remain largely unknown. Here we used a genetically diverse mouse population and applied systems genetics strategies to identify interactions between host and microbe phenotypes including microbial functions, using faecal metagenomics, small intestinal transcripts and caecal lipids that influence microbe-host dynamics. Quantitative trait locus (QTL) mapping identified murine genomic regions associated with variations in bacterial taxa; bacterial functions including motility, sporulation and lipopolysaccharide production and levels of bacterial- and host-derived lipids. We found overlapping QTL for the abundance of Akkermansia muciniphila and caecal levels of ornithine lipids. Follow-up in vitro and in vivo studies revealed that A. muciniphila is a major source of these lipids in the gut, provided evidence that ornithine lipids have immunomodulatory effects and identified intestinal transcripts co-regulated with these traits including Atf3, which encodes for a transcription factor that plays vital roles in modulating metabolism and immunity. Collectively, these results suggest that ornithine lipids are potentially important for A. muciniphila-host interactions and support the role of host genetics as a determinant of responses to gut microbes.

Defining intermediates and redundancies in coenzyme Q precursor biosynthesis.

Coenzyme Q (CoQ), a redox-active lipid essential for oxidative phosphorylation, is synthesized by virtually all cells, but how eukaryotes make the universal CoQ head group precursor 4-hydroxybenzoate (4-HB) from tyrosine is unknown. The first and last steps of this pathway have been defined in Saccharomyces cerevisiae, but the intermediates and enzymes involved in converting 4-hydroxyphenylpyruvate (4-HPP) to 4-hydroxybenzaldehyde (4-HBz) have not been described. Here, we interrogate this pathway with genetic screens, targeted LC-MS, and chemical genetics. We identify three redundant aminotransferases (Bna3, Bat2, and Aat2) that support CoQ biosynthesis in the absence of the established pathway tyrosine aminotransferases, Aro8 and Aro9. We use isotope labeling to identify bona fide tyrosine catabolites, including 4-hydroxyphenylacetate (4-HPA) and 4-hydroxyphenyllactate (4-HPL). Additionally, we find multiple compounds that rescue this pathway when exogenously supplemented, most notably 4-hydroxyphenylacetaldehyde (4-HPAA) and 4-hydroxymandelate (4-HMA). Finally, we show that the Ehrlich pathway decarboxylase Aro10 is dispensable for 4-HB production. These results define new features of 4-HB synthesis in yeast, demonstrate the redundant nature of this pathway, and provide a foundation for further study.

A large-scale genome-lipid association map guides lipid identification.

Despite the crucial roles of lipids in metabolism, we are still at the early stages of comprehensively annotating lipid species and their genetic basis. Mass spectrometry-based discovery lipidomics offers the potential to globally survey lipids and their relative abundances in various biological samples. To discover the genetics of lipid features obtained through high-resolution liquid chromatography-tandem mass spectrometry, we analysed liver and plasma from 384 diversity outbred mice, and quantified 3,283 molecular features. These features were mapped to 5,622 lipid quantitative trait loci and compiled into a public web resource termed LipidGenie. The data are cross-referenced to the human genome and offer a bridge between genetic associations in humans and mice. Harnessing this resource, we used genome-lipid association data as an additional aid to identify a number of lipids, for example gangliosides through their association with B4galnt1, and found evidence for a group of sex-specific phosphatidylcholines through their shared locus. Finally, LipidGenie's ability to query either mass or gene-centric terms suggests acyl-chain-specific functions for proteins of the ABHD family.

Assessing the Viability of Recovery of Hydroxycinnamic Acids from Lignocellulosic Biorefinery Alkaline Pretreatment Waste Streams.

Invited for this month's cover is the research team from the D.O.E. Great Lake Bioenergy Research Center (GLBRC) at the University of Wisconsin-Madison. The cover image shows how a diverse team with expertise in many different fields works together in an integrated fashion to address complex problems. Only when the whole system, from field to the liquid fuels and co-products, is assessed, can we identify the key parameters needed to design an economically viable biorefinery-based economy. Cover art by Chelsea Mamott. The Full Paper itself is available at 10.1002/cssc.201903345.

Assessing the Viability of Recovery of Hydroxycinnamic Acids from Lignocellulosic Biorefinery Alkaline Pretreatment Waste Streams.

The hydroxycinnamic acids p-coumaric acid (pCA) and ferulic acid (FA) add diversity to the portfolio of products produced by using grass-fed lignocellulosic biorefineries. The level of lignin-bound pCA in Zea mays was modified by the alteration of p-coumaroyl-CoA monolignol transferase expression. The biomass was processed in a lab-scale alkaline-pretreatment biorefinery process and the data were used for a baseline technoeconomic analysis to determine where to direct future research efforts to couple plant design to biomass utilization processes. It is concluded that future plant engineering efforts should focus on strategies that ramp up accumulation of one type of hydroxycinnamate (pCA or FA) predominantly and suppress that of the other. Technoeconomic analysis indicates that target extraction titers of one hydroxycinnamic acid need to be >50 g kg-1 biomass, at least five times higher than observed titers for the impure pCA/FA product mixture from wild-type maize. The technical challenge for process engineers is to develop a viable process that requires more than 80 % reduction of the isolation costs.

Multiomic Fermentation Using Chemically Defined Synthetic Hydrolyzates Revealed Multiple Effects of Lignocellulose-Derived Inhibitors on Cell Physiology and Xylose Utilization in Zymomonas mobilis.

Utilization of both C5 and C6 sugars to produce biofuels and bioproducts is a key goal for the development of integrated lignocellulosic biorefineries. Previously we found that although engineered Zymomonas mobilis 2032 was able to ferment glucose to ethanol when fermenting highly concentrated hydrolyzates such as 9% glucan-loading AFEX-pretreated corn stover hydrolyzate (9% ACSH), xylose conversion after glucose depletion was greatly impaired. We hypothesized that impaired xylose conversion was caused by lignocellulose-derived inhibitors (LDIs) in hydrolyzates. To investigate the effects of LDIs on the cellular physiology of Z. mobilis during fermentation of hydrolyzates, including impacts on xylose utilization, we generated synthetic hydrolyzates (SynHs) that contained nutrients and LDIs at concentrations found in 9% ACSH. Comparative fermentations of Z. mobilis 2032 using SynH with or without LDIs were performed, and samples were collected for end product, transcriptomic, metabolomic, and proteomic analyses. Several LDI-specific effects were observed at various timepoints during fermentation including upregulation of sulfur assimilation and cysteine biosynthesis, upregulation of RND family efflux pump systems (ZMO0282-0285) and ZMO1429-1432, downregulation of a Type I secretion system (ZMO0252-0255), depletion of reduced glutathione, and intracellular accumulation of mannose-1P and mannose-6P. Furthermore, when grown in SynH containing LDIs, Z. mobilis 2032 only metabolized ∼50% of xylose, compared to ∼80% in SynH without LDIs, recapitulating the poor xylose utilization observed in 9% ACSH. Our metabolomic data suggest that the overall flux of xylose metabolism is reduced in the presence of LDIs. However, the expression of most genes involved in glucose and xylose assimilation was not affected by LDIs, nor did we observe blocks in glucose and xylose metabolic pathways. Accumulations of intracellular xylitol and xylonic acid was observed in both SynH with and without LDIs, which decreased overall xylose-to-ethanol conversion efficiency. Our results suggest that xylose metabolism in Z. mobilis 2032 may not be able to support the cellular demands of LDI mitigation and detoxification during fermentation of highly concentrated lignocellulosic hydrolyzates with elevated levels of LDIs. Together, our findings identify several cellular responses to LDIs and possible causes of impaired xylose conversion that will enable future strain engineering of Z. mobilis.

Genetic determinants of gut microbiota composition and bile acid profiles in mice.

The microbial communities that inhabit the distal gut of humans and other mammals exhibit large inter-individual variation. While host genetics is a known factor that influences gut microbiota composition, the mechanisms underlying this variation remain largely unknown. Bile acids (BAs) are hormones that are produced by the host and chemically modified by gut bacteria. BAs serve as environmental cues and nutrients to microbes, but they can also have antibacterial effects. We hypothesized that host genetic variation in BA metabolism and homeostasis influence gut microbiota composition. To address this, we used the Diversity Outbred (DO) stock, a population of genetically distinct mice derived from eight founder strains. We characterized the fecal microbiota composition and plasma and cecal BA profiles from 400 DO mice maintained on a high-fat high-sucrose diet for ~22 weeks. Using quantitative trait locus (QTL) analysis, we identified several genomic regions associated with variations in both bacterial and BA profiles. Notably, we found overlapping QTL for Turicibacter sp. and plasma cholic acid, which mapped to a locus containing the gene for the ileal bile acid transporter, Slc10a2. Mediation analysis and subsequent follow-up validation experiments suggest that differences in Slc10a2 gene expression associated with the different strains influences levels of both traits and revealed novel interactions between Turicibacter and BAs. This work illustrates how systems genetics can be utilized to generate testable hypotheses and provide insight into host-microbe interactions.

Accelerating Lipidomic Method Development through in Silico Simulation.

Judicious selection of mass spectrometry (MS) acquisition parameters is essential for effectively profiling the broad diversity and dynamic range of biomolecules. Typically, acquisition parameters are individually optimized to maximally characterize analytes from each new sample matrix. This time-consuming process often ignores the synergistic relationship between MS method parameters, producing suboptimal results. Here we detail the creation of an algorithm which accurately simulates LC-MS/MS lipidomic data acquisition performance for a benchtop quadrupole-Orbitrap MS system. By coupling this simulation tool with a genetic algorithm for constrained parameter optimization, we demonstrate the efficient identification of LC-MS/MS method parameter sets individually suited for specific sample matrices. Finally, we utilize the in silico simulation to examine how continued developments in MS acquisition speed and sensitivity will further increase the power of MS lipidomics as a vital tool for impactful biochemical analysis.

Mapping Lipid Fragmentation for Tailored Mass Spectral Libraries.

Libraries of simulated lipid fragmentation spectra enable the identification of hundreds of unique lipids from complex lipid extracts, even when the corresponding lipid reference standards do not exist. Often, these in silico libraries are generated through expert annotation of spectra to extract and model fragmentation rules common to a given lipid class. Although useful for a given sample source or instrumental platform, the time-consuming nature of this approach renders it impractical for the growing array of dissociation techniques and instrument platforms. Here, we introduce Library Forge, a unique algorithm capable of deriving lipid fragment mass-to-charge (m/z) and intensity patterns directly from high-resolution experimental spectra with minimal user input. Library Forge exploits the modular construction of lipids to generate m/z transformed spectra in silico which reveal the underlying fragmentation pathways common to a given lipid class. By learning these fragmentation patterns directly from observed spectra, the algorithm increases lipid spectral matching confidence while reducing spectral library development time from days to minutes. We embed the algorithm within the preexisting lipid analysis architecture of LipiDex to integrate automated and robust library generation within a comprehensive LC-MS/MS lipidomics workflow. Graphical Abstract.

S-adenosylmethionine biosynthesis is a targetable metabolic vulnerability of cancer stem cells.

PURPOSE: Many transformed cells and embryonic stem cells are dependent on the biosynthesis of the universal methyl-donor S-adenosylmethionine (SAM) from methionine by the enzyme MAT2A to maintain their epigenome. We hypothesized that cancer stem cells (CSCs) rely on SAM biosynthesis and that the combination of methionine depletion and MAT2A inhibition would eradicate CSCs. METHODS: Human triple (ER/PR/HER2)-negative breast carcinoma (TNBC) cell lines were cultured as CSC-enriched mammospheres in control or methionine-free media. MAT2A was inhibited with siRNAs or cycloleucine. The effects of methionine restriction and/or MAT2A inhibition on the formation of mammospheres, the expression of CSC markers (CD44hi/C24low), MAT2A and CSC transcriptional regulators, apoptosis induction and histone modifications were determined. A murine model of metastatic TNBC was utilized to evaluate the effects of dietary methionine restriction, MAT2A inhibition and the combination. RESULTS: Methionine restriction inhibited mammosphere formation and reduced the CD44hi/C24low CSC population; these effects were partly rescued by SAM. Methionine depletion induced MAT2A expression (mRNA and protein) and sensitized CSCs to inhibition of MAT2A (siRNAs or cycloleucine). Cycloleucine enhanced the effects of methionine depletion on H3K4me3 demethylation and suppression of Sox9 expression. Dietary methionine restriction induced MAT2A expression in mammary tumors, and the combination of methionine restriction and cycloleucine was more effective than either alone at suppressing primary and lung metastatic tumor burden in a murine TNBC model. CONCLUSIONS: Our findings point to SAM biosynthesis as a unique metabolic vulnerability of CSCs that can be targeted by combining methionine depletion with MAT2A inhibition to eradicate drug-resistant CSCs.

Coenzyme Q biosynthetic proteins assemble in a substrate-dependent manner into domains at ER-mitochondria contacts.

Coenzyme Q (CoQ) lipids are ancient electron carriers that, in eukaryotes, function in the mitochondrial respiratory chain. In mitochondria, CoQ lipids are built by an inner membrane-associated, multicomponent, biosynthetic pathway via successive steps of isoprenyl tail polymerization, 4-hydroxybenzoate head-to-tail attachment, and head modification, resulting in the production of CoQ. In yeast, we discovered that head-modifying CoQ pathway components selectively colocalize to multiple resolvable domains in vivo, representing supramolecular assemblies. In cells engineered with conditional ON or OFF CoQ pathways, domains were strictly correlated with CoQ production and substrate flux, respectively, indicating that CoQ lipid intermediates are required for domain formation. Mitochondrial CoQ domains were also observed in human cells, underscoring their conserved functional importance. CoQ domains within cells were highly enriched adjacent to ER-mitochondria contact sites. Together, our data suggest that CoQ domains function to facilitate substrate accessibility for processive and efficient CoQ production and distribution in cells.

LipiDex: An Integrated Software Package for High-Confidence Lipid Identification.

State-of-the-art proteomics software routinely quantifies thousands of peptides per experiment with minimal need for manual validation or processing of data. For the emerging field of discovery lipidomics via liquid chromatography-tandem mass spectrometry (LC-MS/MS), comparably mature informatics tools do not exist. Here, we introduce LipiDex, a freely available software suite that unifies and automates all stages of lipid identification, reducing hands-on processing time from hours to minutes for even the most expansive datasets. LipiDex utilizes flexible in silico fragmentation templates and lipid-optimized MS/MS spectral matching routines to confidently identify and track hundreds of lipid species and unknown compounds from diverse sample matrices. Unique spectral and chromatographic peak purity algorithms accurately quantify co-isolation and co-elution of isobaric lipids, generating identifications that match the structural resolution afforded by the LC-MS/MS experiment. During final data filtering, ionization artifacts are removed to significantly reduce dataset redundancy. LipiDex interfaces with several LC-MS/MS software packages, enabling robust lipid identification to be readily incorporated into pre-existing data workflows.

Multi-omic Mitoprotease Profiling Defines a Role for Oct1p in Coenzyme Q Production.

Mitoproteases are becoming recognized as key regulators of diverse mitochondrial functions, although their direct substrates are often difficult to discern. Through multi-omic profiling of diverse Saccharomyces cerevisiae mitoprotease deletion strains, we predicted numerous associations between mitoproteases and distinct mitochondrial processes. These include a strong association between the mitochondrial matrix octapeptidase Oct1p and coenzyme Q (CoQ) biosynthesis-a pathway essential for mitochondrial respiration. Through Edman sequencing and in vitro and in vivo biochemistry, we demonstrated that Oct1p directly processes the N terminus of the CoQ-related methyltransferase, Coq5p, which markedly improves its stability. A single mutation to the Oct1p recognition motif in Coq5p disrupted its processing in vivo, leading to CoQ deficiency and respiratory incompetence. This work defines the Oct1p processing of Coq5p as an essential post-translational event for proper CoQ production. Additionally, our data visualization tool enables efficient exploration of mitoprotease profiles that can serve as the basis for future mechanistic investigations.

Ptc7p Dephosphorylates Select Mitochondrial Proteins to Enhance Metabolic Function.

Proper maintenance of mitochondrial activity is essential for metabolic homeostasis. Widespread phosphorylation of mitochondrial proteins may be an important element of this process; yet, little is known about which enzymes control mitochondrial phosphorylation or which phosphosites have functional impact. We investigate these issues by disrupting Ptc7p, a conserved but largely uncharacterized mitochondrial matrix PP2C-type phosphatase. Loss of Ptc7p causes respiratory growth defects concomitant with elevated phosphorylation of select matrix proteins. Among these, Δptc7 yeast exhibit an increase in phosphorylation of Cit1p, the canonical citrate synthase of the tricarboxylic acid (TCA) cycle, that diminishes its activity. We find that phosphorylation of S462 can eliminate Cit1p enzymatic activity likely by disrupting its proper dimerization, and that Ptc7p-driven dephosphorylation rescues Cit1p activity. Collectively, our work connects Ptc7p to an essential TCA cycle function and to additional phosphorylation events that may affect mitochondrial activity inadvertently or in a regulatory manner.

Ongoing resolution of duplicate gene functions shapes the diversification of a metabolic network.

The evolutionary mechanisms leading to duplicate gene retention are well understood, but the long-term impacts of paralog differentiation on the regulation of metabolism remain underappreciated. Here we experimentally dissect the functions of two pairs of ancient paralogs of the GALactose sugar utilization network in two yeast species. We show that the Saccharomyces uvarum network is more active, even as over-induction is prevented by a second co-repressor that the model yeast Saccharomyces cerevisiae lacks. Surprisingly, removal of this repression system leads to a strong growth arrest, likely due to overly rapid galactose catabolism and metabolic overload. Alternative sugars, such as fructose, circumvent metabolic control systems and exacerbate this phenotype. We further show that S. cerevisiae experiences homologous metabolic constraints that are subtler due to how the paralogs have diversified. These results show how the functional differentiation of paralogs continues to shape regulatory network architectures and metabolic strategies long after initial preservation.

NANOG is multiply phosphorylated and directly modified by ERK2 and CDK1 in vitro.

NANOG is a divergent homeobox protein and a core component of the transcriptional circuitry that sustains pluripotency and self-renewal. Although NANOG has been extensively studied on the transcriptional level, little is known regarding its posttranslational regulation, likely due to its low abundance and challenging physical properties. Here, we identify eleven phosphorylation sites on endogenous human NANOG, nine of which mapped to single amino acids. To screen for the signaling molecules that impart these modifications, we developed the multiplexed assay for kinase specificity (MAKS). MAKS simultaneously tests activity for up to ten kinases while directly identifying the substrate and exact site of phosphorylation. Using MAKS, we discovered site-specific phosphorylation by ERK2 and CDK1/CyclinA2, providing a putative link between key signaling pathways and NANOG.

Multipurpose dissociation cell for enhanced ETD of intact protein species.

We describe and characterize an improved implementation of ETD on a modified hybrid linear ion trap-Orbitrap instrument. Instead of performing ETD in the mass-analyzing quadrupole linear ion trap (A-QLT), the instrument collision cell was modified to enable ETD. We partitioned the collision cell into a multi-section rf ion storage and transfer device to enable injection and simultaneous separate storage of precursor and reagent ions. Application of a secondary (axial) confinement voltage to the cell end lens electrodes enables charge-sign independent trapping for ion-ion reactions. The approximately 2-fold higher quadrupole field frequency of this cell relative to that of the A-QLT enables higher reagent ion densities and correspondingly faster ETD reactions, and, with the collision cell's longer axial dimensions, larger populations of precursor ions may be reacted. The higher ion capacity of the collision cell permits the accumulation and reaction of multiple full loads of precursor ions from the A-QLT followed by FT Orbitrap m/z analysis of the ETD product ions. This extends the intra-scan dynamic range by increasing the maximum number of product ions in a single MS/MS event. For analyses of large peptide/small protein precursor cations, this reduces or eliminates the need for spectral averaging to achieve acceptable ETD product ion signal-to-noise levels. Using larger ion populations, we demonstrate improvements in protein sequence coverage and aggregate protein identifications in LC-MS/MS analysis of intact protein species as compared to the standard ETD implementation.

Characterization and quantification of intact 26S proteasome proteins by real-time measurement of intrinsic fluorescence prior to top-down mass spectrometry.

Quantification of gas-phase intact protein ions by mass spectrometry (MS) is impeded by highly-variable ionization, ion transmission, and ion detection efficiencies. Therefore, quantification of proteins using MS-associated techniques is almost exclusively done after proteolysis where peptides serve as proxies for estimating protein abundance. Advances in instrumentation, protein separations, and informatics have made large-scale sequencing of intact proteins using top-down proteomics accessible to the proteomics community; yet quantification of proteins using a top-down workflow has largely been unaddressed. Here we describe a label-free approach to determine the abundance of intact proteins separated by nanoflow liquid chromatography prior to MS analysis by using solution-phase measurements of ultraviolet light-induced intrinsic fluorescence (UV-IF). UV-IF is measured directly at the electrospray interface just prior to the capillary exit where proteins containing at least one tryptophan residue are readily detected. UV-IF quantification was demonstrated using commercially available protein standards and provided more accurate and precise protein quantification than MS ion current. We evaluated the parallel use of UV-IF and top-down tandem MS for quantification and identification of protein subunits and associated proteins from an affinity-purified 26S proteasome sample from Arabidopsis thaliana. We identified 26 unique proteins and quantified 13 tryptophan-containing species. Our analyses discovered previously unidentified N-terminal processing of the ß6 (PBF1) and ß7 (PBG1) subunit - such processing of PBG1 may generate a heretofore unknown additional protease active site upon cleavage. In addition, our approach permitted the unambiguous identification and quantification both isoforms of the proteasome-associated protein DSS1.

Comprehensive mass spectrometric mapping of the hydroxylated amino acid residues of the α1(V) collagen chain.

BACKGROUND: α1(V) is an extensively modified collagen chain important in disease. RESULTS: Comprehensive mapping of α1(V) post-translational modifications reveals unexpectedly large numbers of X-position hydroxyprolines in Gly-X-Y amino acid triplets. CONCLUSION: The unexpected abundance of X-position hydroxyprolines suggests a mechanism for differential modification of collagen properties. SIGNIFICANCE: Positions, numbers, and occupancy of modified sites can provide insights into α1(V) biological properties. Aberrant expression of the type V collagen α1(V) chain can underlie the connective tissue disorder classic Ehlers-Danlos syndrome, and autoimmune responses against the α1(V) chain are linked to lung transplant rejection and atherosclerosis. The α1(V) collagenous COL1 domain is thought to contain greater numbers of post-translational modifications (PTMs) than do similar domains of other fibrillar collagen chains, PTMs consisting of hydroxylated prolines and lysines, the latter of which can be glycosylated. These types of PTMs can contribute to epitopes that underlie immune responses against collagens, and the high level of PTMs may contribute to the unique biological properties of the α1(V) chain. Here we use high resolution mass spectrometry to map such PTMs in bovine placental α1(V) and human recombinant pro-α1(V) procollagen chains. Findings include the locations of those PTMs that vary and those PTMs that are invariant between these α1(V) chains from widely divergent sources. Notably, an unexpectedly large number of hydroxyproline residues were mapped to the X-positions of Gly-X-Y triplets, contrary to expectations based on previous amino acid analyses of hydrolyzed α1(V) chains from various tissues. We attribute this difference to the ability of tandem mass spectrometry coupled to nanoflow chromatographic separations to detect lower-level PTM combinations with superior sensitivity and specificity. The data are consistent with the presence of a relatively large number of 3-hydroxyproline sites with less than 100% occupancy, suggesting a previously unknown mechanism for the differential modification of α1(V) chain and type V collagen properties.