iLipidome: enhancing statistical power and interpretability using hidden biosynthetic interdependencies in the lipidome.

Numerous biological processes and diseases are influenced by lipid composition. Advances in lipidomics are elucidating their roles, but analyzing and interpreting lipidomics data at the systems level remain challenging. To address this, we present iLipidome, a method for analyzing lipidomics data in the context of the lipid biosynthetic network, thus accounting for the interdependence of measured lipids. iLipidome enhances statistical power, enables reliable clustering and lipid enrichment analysis, and links lipidomic changes to their genetic origins. We applied iLipidome to investigate mechanisms driving changes in cellular lipidomes following supplementation of docosahexaenoic acid (DHA) and successfully identified the genetic causes of alterations. We further demonstrated how iLipidome can disclose enzyme-substrate specificity and pinpoint prospective glioblastoma therapeutic targets. Finally, iLipidome enabled us to explore underlying mechanisms of cardiovascular disease and could guide the discovery of early lipid biomarkers. Thus, iLipidome can assist researchers studying the essence of lipidomic data and advance the field of lipid biology.

LeGenD: determining N-glycoprofiles using an explainable AI-leveraged model with lectin profiling.

Glycosylation affects many vital functions of organisms. Therefore, its surveillance is critical from basic science to biotechnology, including biopharmaceutical development and clinical diagnostics. However, conventional glycan structure analysis faces challenges with throughput and cost. Lectins offer an alternative approach for analyzing glycans, but they only provide glycan epitopes and not full glycan structure information. To overcome these limitations, we developed LeGenD, a lectin and AI-based approach to predict N-glycan structures and determine their relative abundance in purified proteins based on lectin-binding patterns. We trained the LeGenD model using 309 glycoprofiles from 10 recombinant proteins, produced in 30 glycoengineered CHO cell lines. Our approach accurately reconstructed experimentally-measured N-glycoprofiles of bovine Fetuin B and IgG from human sera. Explanatory AI analysis with SHapley Additive exPlanations (SHAP) helped identify the critical lectins for glycoprofile predictions. Our LeGenD approach thus presents an alternative approach for N-glycan analysis.

Mechanotransduction-induced interplay between phospholamban and yes-activated protein induces smooth muscle cell hypertrophy.

The gastrointestinal system is a hollow organ affected by fibrostenotic diseases that cause volumetric compromise of the lumen via smooth muscle hypertrophy and fibrosis. Many of the driving mechanisms remain unclear. Yes-associated protein-1 (YAP) is a critical mechanosensory transcriptional regulator that mediates cell hypertrophy in response to elevated extracellular rigidity. In the type 2 inflammatory disorder, eosinophilic esophagitis (EoE), phospholamban (PLN) can induce smooth muscle cell hypertrophy. We used EoE as a disease model for understanding a mechanistic pathway in which PLN and YAP interact in response to rigid extracellular substrate to induce smooth muscle cell hypertrophy. PLN-induced YAP nuclear sequestration in a feed-forward loop caused increased cell size in response to a rigid substrate. This mechanism of rigidity sensing may have previously unappreciated clinical implications for PLN-expressing hollow systems such as the esophagus and heart.

LipidSIM: Inferring mechanistic lipid biosynthesis perturbations from lipidomics with a flexible, low-parameter, Markov modeling framework.

Lipid metabolism is a complex and dynamic system involving numerous enzymes at the junction of multiple metabolic pathways. Disruption of these pathways leads to systematic dyslipidemia, a hallmark of many pathological developments, such as nonalcoholic steatohepatitis and diabetes. Recent advances in computational tools can provide insights into the dysregulation of lipid biosynthesis, but limitations remain due to the complexity of lipidomic data, limited knowledge of interactions among involved enzymes, and technical challenges in standardizing across different lipid types. Here, we present a low-parameter, biologically interpretable framework named Lipid Synthesis Investigative Markov model (LipidSIM), which models and predicts the source of perturbations in lipid biosynthesis from lipidomic data. LipidSIM achieves this by accounting for the interdependency between the lipid species via the lipid biosynthesis network and generates testable hypotheses regarding changes in lipid biosynthetic reactions. This feature allows the integration of lipidomics with other omics types, such as transcriptomics, to elucidate the direct driving mechanisms of altered lipidomes due to treatments or disease progression. To demonstrate the value of LipidSIM, we first applied it to hepatic lipidomics following Keap1 knockdown and found that changes in mRNA expression of the lipid pathways were consistent with the LipidSIM-predicted fluxes. Second, we used it to study lipidomic changes following intraperitoneal injection of CCl4 to induce fast NAFLD/NASH development and the progression of fibrosis and hepatic cancer. Finally, to show the power of LipidSIM for classifying samples with dyslipidemia, we used a Dgat2-knockdown study dataset. Thus, we show that as it demands no a priori knowledge of enzyme kinetics, LipidSIM is a valuable and intuitive framework for extracting biological insights from complex lipidomic data.

GlycoMME, a Markov modeling platform for studying N-glycosylation biosynthesis from glycomics data.

Variations in N-glycosylation, which is crucial to glycoprotein functions, impact many diseases and the safety and efficacy of biotherapeutic drugs. Here, we present a protocol for using GlycoMME (Glycosylation Markov Model Evaluator) to study N-glycosylation biosynthesis from glycomics data. We describe steps for annotating glycomics data and quantifying perturbations to N-glycan biosynthesis with interpretable models. We then detail procedures to predict the impact of mutations in disease or potential glycoengineering strategies in drug development. For complete details on the use and execution of this protocol, please refer to Liang et al. (2020).1.

Preparing glycomics data for robust statistical analysis with GlyCompareCT.

GlyCompareCT is a portable command-line tool to facilitate downstream glycomic data analyses, by addressing data inherent sparsity and non-independence. Inputting glycan abundances, users can run GlyCompareCT with one line of code to obtain the abundances of a minimal substructure set, named glycomotif, thereby quantifying hidden biosynthetic relationships between measured glycans. Optional parameters tuning and annotation are supported for personal preference. For complete details on the use and execution of this protocol, please refer to Bao et al. (2021).1.

CHOGlycoNET: Comprehensive glycosylation reaction network for CHO cells.

Chinese hamster ovary (CHO) cells are extensively used for the production of glycoprotein therapeutics proteins, for which N-linked glycans are a critical quality attribute due to their influence on activity and immunogenicity. Manipulation of protein glycosylation is commonly achieved through cell or process engineering, which are often guided by mathematical models. However, each study considers a unique glycosylation reaction network that is tailored around the cell line and product at hand. Herein, we use 200 glycan datasets for both recombinantly produced and native proteins from different CHO cell lines to reconstruct a comprehensive reaction network, CHOGlycoNET, based on the individual minimal reaction networks describing each dataset. CHOGlycoNET is used to investigate the distribution of mannosidase and glycosyltransferase enzymes in the Golgi apparatus and identify key network reactions using machine learning and dimensionality reduction techniques. CHOGlycoNET can be used for accelerating glycomodel development and predicting the effect of glycoengineering strategies. Finally, CHOGlycoNET is wrapped in a SBML file to be used as a standalone model or in combination with CHO cell genome scale models.

A predictive ensemble classifier for the gene expression diagnosis of ASD at ages 1 to 4 years.

Autism Spectrum Disorder (ASD) diagnosis remains behavior-based and the median age of diagnosis is ~52 months, nearly 5 years after its first-trimester origin. Accurate and clinically-translatable early-age diagnostics do not exist due to ASD genetic and clinical heterogeneity. Here we collected clinical, diagnostic, and leukocyte RNA data from 240 ASD and typically developing (TD) toddlers (175 toddlers for training and 65 for test). To identify gene expression ASD diagnostic classifiers, we developed 42,840 models composed of 3570 gene expression feature selection sets and 12 classification methods. We found that 742 models had AUC-ROC ≥ 0.8 on both Training and Test sets. Weighted Bayesian model averaging of these 742 models yielded an ensemble classifier model with accurate performance in Training and Test gene expression datasets with ASD diagnostic classification AUC-ROC scores of 85-89% and AUC-PR scores of 84-92%. ASD toddlers with ensemble scores above and below the overall ASD ensemble mean of 0.723 (on a scale of 0 to 1) had similar diagnostic and psychometric scores, but those below this ASD ensemble mean had more prenatal risk events than TD toddlers. Ensemble model feature genes were involved in cell cycle, inflammation/immune response, transcriptional gene regulation, cytokine response, and PI3K-AKT, RAS and Wnt signaling pathways. We additionally collected targeted DNA sequencing smMIPs data on a subset of ASD risk genes from 217 of the 240 ASD and TD toddlers. This DNA sequencing found about the same percentage of SFARI Level 1 and 2 ASD risk gene mutations in TD (12 of 105) as in ASD (13 of 112) toddlers, and classification based only on the presence of mutation in these risk genes performed at a chance level of 49%. By contrast, the leukocyte ensemble gene expression classifier correctly diagnostically classified 88% of TD and ASD toddlers with ASD risk gene mutations. Our ensemble ASD gene expression classifier is diagnostically predictive and replicable across different toddler ages, races, and ethnicities; out-performs a risk gene mutation classifier; and has potential for clinical translation.

Non-protective immune imprint underlies failure of Staphylococcus aureus IsdB vaccine.

Humans frequently encounter Staphylococcus aureus (SA) throughout life. Animal studies have yielded SA candidate vaccines, yet all human SA vaccine trials have failed. One notable vaccine "failure" targeted IsdB, critical for host iron acquisition. We explored a fundamental difference between humans and laboratory animals-natural SA exposure. Recapitulating the failed phase III IsdB vaccine trial, mice previously infected with SA do not mount protective antibody responses to vaccination, unlike naive animals. Non-protective antibodies exhibit increased α2,3 sialylation that blunts opsonophagocytosis and preferentially targets a non-protective IsdB domain. IsdB vaccination of SA-infected mice recalls non-neutralizing humoral responses, further reducing vaccine efficacy through direct antibody competition. IsdB vaccine interference was overcome by immunization against the IsdB heme-binding domain. Purified human IsdB-specific antibodies also blunt IsdB passive immunization, and additional SA vaccines are susceptible to SA pre-exposure. Thus, failed anti-SA immunization trials could be explained by non-protective imprint from prior host-SA interaction.

Artificial intelligence in the analysis of glycosylation data.

Glycans are complex, yet ubiquitous across biological systems. They are involved in diverse essential organismal functions. Aberrant glycosylation may lead to disease development, such as cancer, autoimmune diseases, and inflammatory diseases. Glycans, both normal and aberrant, are synthesized using extensive glycosylation machinery, and understanding this machinery can provide invaluable insights for diagnosis, prognosis, and treatment of various diseases. Increasing amounts of glycomics data are being generated thanks to advances in glycoanalytics technologies, but to maximize the value of such data, innovations are needed for analyzing and interpreting large-scale glycomics data. Artificial intelligence (AI) provides a powerful analysis toolbox in many scientific fields, and here we review state-of-the-art AI approaches on glycosylation analysis. We further discuss how models can be analyzed to gain mechanistic insights into glycosylation machinery and how the machinery shapes glycans under different scenarios. Finally, we propose how to leverage the gained knowledge for developing predictive AI-based models of glycosylation. Thus, guiding future research of AI-based glycosylation model development will provide valuable insights into glycosylation and glycan machinery.

Elucidating Human Milk Oligosaccharide biosynthetic genes through network-based multi-omics integration.

Human Milk Oligosaccharides (HMOs) are abundant carbohydrates fundamental to infant health and development. Although these oligosaccharides were discovered more than half a century ago, their biosynthesis in the mammary gland remains largely uncharacterized. Here, we use a systems biology framework that integrates glycan and RNA expression data to construct an HMO biosynthetic network and predict glycosyltransferases involved. To accomplish this, we construct models describing the most likely pathways for the synthesis of the oligosaccharides accounting for >95% of the HMO content in human milk. Through our models, we propose candidate genes for elongation, branching, fucosylation, and sialylation of HMOs. Our model aggregation approach recovers 2 of 2 previously known gene-enzyme relations and 2 of 3 empirically confirmed gene-enzyme relations. The top genes we propose for the remaining 5 linkage reactions are consistent with previously published literature. These results provide the molecular basis of HMO biosynthesis necessary to guide progress in HMO research and application with the goal of understanding and improving infant health and development.

LIGHT controls distinct homeostatic and inflammatory gene expression profiles in esophageal fibroblasts via differential HVEM and LTßR-mediated mechanisms.

Fibroblasts mediate tissue remodeling in eosinophilic esophagitis (EoE), a chronic allergen-driven inflammatory pathology. Diverse fibroblast subtypes with homeostasis-regulating or inflammatory profiles have been recognized in various tissues, but which mediators induce these alternate differentiation states remain largely unknown. We recently identified that TNFSF14/LIGHT promotes an inflammatory esophageal fibroblast in vitro. Herein we used esophageal biopsies and primary fibroblasts to investigate the role of the LIGHT receptors, herpes virus entry mediator (HVEM) and lymphotoxin-beta receptor (LTßR), and their downstream activated pathways, in EoE. In addition to promoting inflammatory gene expression, LIGHT down-regulated homeostatic factors including WNTs, BMPs and type 3 semaphorins. In vivo, WNT2B+ fibroblasts were decreased while ICAM-1+ and IL-34+ fibroblasts were expanded in EoE, suggesting that a LIGHT-driven gene signature was imprinted in EoE versus normal esophageal fibroblasts. HVEM and LTßR overexpression and deficiency experiments demonstrated that HVEM regulates a limited subset of LIGHT targets, whereas LTßR controls all transcriptional effects. Pharmacologic blockade of the non-canonical NIK/p100/p52-mediated NF-κB pathway potently silenced LIGHT's transcriptional effects, with a lesser role found for p65 canonical NF-κB. Collectively, our results show that LIGHT promotes differentiation of esophageal fibroblasts toward an inflammatory phenotype and represses homeostatic gene expression via a LTßR-NIK-p52 NF-κB dominant pathway.

Model-based assessment of mammalian cell metabolic functionalities using omics data.

Omics experiments are ubiquitous in biological studies, leading to a deluge of data. However, it is still challenging to connect changes in these data to changes in cell functions because of complex interdependencies between genes, proteins, and metabolites. Here, we present a framework allowing researchers to infer how metabolic functions change on the basis of omics data. To enable this, we curated and standardized lists of metabolic tasks that mammalian cells can accomplish. Genome-scale metabolic networks were used to define gene sets associated with each metabolic task. We further developed a framework to overlay omics data on these sets and predict pathway usage for each metabolic task. We demonstrated how this approach can be used to quantify metabolic functions of diverse biological samples from the single cell to whole tissues and organs by using multiple transcriptomic datasets. To facilitate its adoption, we integrated the approach into GenePattern (www.genepattern.org-CellFie).

Correcting for sparsity and interdependence in glycomics by accounting for glycan biosynthesis.

Glycans are fundamental cellular building blocks, involved in many organismal functions. Advances in glycomics are elucidating the essential roles of glycans. Still, it remains challenging to properly analyze large glycomics datasets, since the abundance of each glycan is dependent on many other glycans that share many intermediate biosynthetic steps. Furthermore, the overlap of measured glycans can be low across samples. We address these challenges with GlyCompare, a glycomic data analysis approach that accounts for shared biosynthetic steps for all measured glycans to correct for sparsity and non-independence in glycomics, which enables direct comparison of different glycoprofiles and increases statistical power. Using GlyCompare, we study diverse N-glycan profiles from glycoengineered erythropoietin. We obtain biologically meaningful clustering of mutant cell glycoprofiles and identify knockout-specific effects of fucosyltransferase mutants on tetra-antennary structures. We further analyze human milk oligosaccharide profiles and find mother's fucosyltransferase-dependent secretor-status indirectly impact the sialylation. Finally, we apply our method on mucin-type O-glycans, gangliosides, and site-specific compositional glycosylation data to reveal tissues and disease-specific glycan presentations. Our substructure-oriented approach will enable researchers to take full advantage of the growing power and size of glycomics data.

Dysregulation of the secretory pathway connects Alzheimer's disease genetics to aggregate formation.

Amyloid disorders such as Alzheimer's disease (AD) involve the aggregation of secreted proteins. However, it is largely unclear how secretory-pathway proteins contribute to amyloid formation. We developed a systems biology framework integrating expression data with protein-protein interaction networks to estimate a tissue's fitness for producing specific secreted proteins and analyzed the fitness of the secretory pathway of various brain regions and cell types for synthesizing the AD-associated amyloid precursor protein (APP). While key amyloidogenic pathway components were not differentially expressed in AD brains, we found Aß deposition correlates with systemic down- and upregulation of the secretory-pathway components proximal to APP and amyloidogenic secretases, respectively, in AD. Our analyses suggest that perturbations from three AD risk loci cascade through the APP secretory-support network and into the endocytosis pathway, connecting amyloidogenesis to dysregulation of secretory-pathway components supporting APP and suggesting novel therapeutic targets for AD. A record of this paper's transparent peer review process is included in the supplemental information.

Systems glycobiology for discovering drug targets, biomarkers, and rational designs for glyco-immunotherapy.

Cancer immunotherapy has revolutionized treatment and led to an unprecedented wave of immuno-oncology research during the past two decades. In 2018, two pioneer immunotherapy innovators, Tasuku Honjo and James P. Allison, were awarded the Nobel Prize for their landmark cancer immunotherapy work regarding "cancer therapy by inhibition of negative immune regulation" -CTLA4 and PD-1 immune checkpoints. However, the challenge in the coming decade is to develop cancer immunotherapies that can more consistently treat various patients and cancer types. Overcoming this challenge requires a systemic understanding of the underlying interactions between immune cells, tumor cells, and immunotherapeutics. The role of aberrant glycosylation in this process, and how it influences tumor immunity and immunotherapy is beginning to emerge. Herein, we review current knowledge of miRNA-mediated regulatory mechanisms of glycosylation machinery, and how these carbohydrate moieties impact immune cell and tumor cell interactions. We discuss these insights in the context of clinical findings and provide an outlook on modulating the regulation of glycosylation to offer new therapeutic opportunities. Finally, in the coming age of systems glycobiology, we highlight how emerging technologies in systems glycobiology are enabling deeper insights into cancer immuno-oncology, helping identify novel drug targets and key biomarkers of cancer, and facilitating the rational design of glyco-immunotherapies. These hold great promise clinically in the immuno-oncology field.

Genome-wide screens uncover KDM2B as a modifier of protein binding to heparan sulfate.

Heparan sulfate (HS) proteoglycans bind extracellular proteins that participate in cell signaling, attachment and endocytosis. These interactions depend on the arrangement of sulfated sugars in the HS chains generated by well-characterized biosynthetic enzymes; however, the regulation of these enzymes is largely unknown. We conducted genome-wide CRISPR-Cas9 screens with a small-molecule ligand that binds to HS. Screening of A375 melanoma cells uncovered additional genes and pathways impacting HS formation. The top hit was the epigenetic factor KDM2B, a histone demethylase. KDM2B inactivation suppressed multiple HS sulfotransferases and upregulated the sulfatase SULF1. These changes differentially affected the interaction of HS-binding proteins. KDM2B-deficient cells displayed decreased growth rates, which was rescued by SULF1 inactivation. In addition, KDM2B deficiency altered the expression of many extracellular matrix genes. Thus, KDM2B controls proliferation of A375 cells through the regulation of HS structure and serves as a master regulator of the extracellular matrix.

A unique esophageal extracellular matrix proteome alters normal fibroblast function in severe eosinophilic esophagitis.

BACKGROUND: Eosinophilic esophagitis (EoE) is a chronic TH2 disorder complicated by tissue fibrosis and loss of esophageal luminal patency. The fibrostenotic esophagus does not respond well to therapy, but profibrotic therapeutic targets are largely unclear. OBJECTIVE: Our aim was to utilize proteomics and primary cells as a novel approach to determine relevant profibrotic factors. METHODS: We utilized primary esophageal EoE and normal fibroblasts, their derivative extracellular matrixes (ECMs), an approach of fibroblast culture on autologous versus nonautologous ECM, and proteomics to elucidate EoE ECM proteins that dysregulate cellular function. RESULTS: We cultured esophageal fibroblasts from normal esophagi and esophagi from patients with severe EoE on autologous versus nonautologous ECM. The EoE ECM proteome shifted normal esophageal fibroblast protein expression. Proteomic analysis demonstrated that thrombospondin-1 is detected only in the EoE ECM, is central in the EoE ECM protein-protein interactome, is found at significantly elevated levels in biopsy specimens from patients with active EoE, and induces fibroblast collagen I production. CONCLUSION: Fibroblasts from patients with EoE secrete a unique ECM proteome that reflects their in vivo state and induces collagen I and α-smooth muscle actin protein expression from normal fibroblasts. Thrombospondin-1 is a previously unappreciated profibrotic molecule in EoE.

Type 2 Immunity and Age Modify Gene Expression of Coronavirus-induced Disease 2019 Receptors in Eosinophilic Gastrointestinal Disorders.

ABSTRACT: Infection with severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) can lead to coronavirus-induced disease 2019 (COVID-19). The gastrointestinal (GI) tract is now an appreciated portal of infection. SARS-CoV-2 enters host cells via angiotensin-converting enzyme-2 (ACE2) and the serine protease TMPRSS2. Eosinophilic gastrointestinal disorders (EGIDs) are inflammatory conditions caused by chronic type 2 (T2) inflammation. the effects of the T2 atopic inflammatory milieu on SARS-COV-2 viral entry gene expression in the GI tract is poorly understood. We analyzed tissue ACE2 and TMPRSS2 gene expression in pediatric eosinophilic esophagitis (EoE), eosinophilic gastritis (EG), and in normal adult esophagi using publicly available RNA-sequencing datasets. Similar to findings evaluating the airway, there was no difference in tissue ACE2/TMPRSS2 expression in EoE or EG when compared with control non-EoE/EG esophagus/stomach. ACE2 gene expression was significantly lower in esophagi from children with or without EoE and from adults with EoE as compared with normal adult esophagi. Type 2 immunity and pediatric age could be protective for infection by SARS-CoV-2 in the gastrointestinal tract because of decreased expression of ACE2.