Shuzhao Li

Molecular signatures of antibody responses derived from a systems biology study of five human vaccines.

Many vaccines induce protective immunity via antibodies. Systems biology approaches have been used to determine signatures that can be used to predict vaccine-induced immunity in humans, but whether there is a universal signature that can be used to predict antibody responses to any vaccine is unknown. Here we did systems analyses of immune responses to the polysaccharide and conjugate vaccines against meningococcus in healthy adults, in the broader context of published studies of vaccines against yellow fever virus and influenza virus. To achieve this, we did a large-scale network integration of publicly available human blood transcriptomes and systems-scale databases in specific biological contexts and deduced a set of transcription modules in blood. Those modules revealed distinct transcriptional signatures of antibody responses to different classes of vaccines, which provided key insights into primary viral, protein recall and anti-polysaccharide responses. Our results elucidate the early transcriptional programs that orchestrate vaccine immunity in humans and demonstrate the power of integrative network modeling.

Autophagy is essential for effector CD8 + T cell survival and memory formation

The importance of autophagy in the generation of memory CD8+ T cells in vivo is not well defined. We report here that autophagy was dynamically regulated in virus-specific CD8+ T cells during acute infection of mice with lymphocytic choriomeningitis virus. In contrast to the current paradigm, autophagy decreased in activated proliferating effector CD8+ T cells and was then upregulated when the cells stopped dividing just before the contraction phase. Consistent with those findings, deletion of the gene encoding either of the autophagy-related molecules Atg5 or Atg7 had little to no effect on the proliferation and function of effector cells, but these autophagy-deficient effector cells had survival defects that resulted in compromised formation of memory T cells. Our studies define when autophagy is needed during effector and memory differentiation and warrant reexamination of the relationship between T cell activation and autophagy.

Predicting Network Activity from High Throughput Metabolomics

The functional interpretation of high throughput metabolomics by mass spectrometry is hindered by the identification of metabolites, a tedious and challenging task. We present a set of computational algorithms which, by leveraging the collective power of metabolic pathways and networks, predict functional activity directly from spectral feature tables without a priori identification of metabolites. The algorithms were experimentally validated on the activation of innate immune cells.

The amino acid sensor GCN2 controls gut inflammation by inhibiting inflammasome activation

The integrated stress response (ISR) is a homeostatic mechanism by which eukaryotic cells sense and respond to stress-inducing signals, such as amino acid starvation. General controlled non-repressed (GCN2) kinase is a key orchestrator of the ISR, and modulates protein synthesis in response to amino acid starvation. Here we demonstrate in mice that GCN2 controls intestinal inflammation by suppressing inflammasome activation. Enhanced activation of ISR was observed in intestinal antigen presenting cells (APCs) and epithelial cells during amino acid starvation, or intestinal inflammation. Genetic deletion of Gcn2 (also known as Eif2ka4) in CD11c+ APCs or intestinal epithelial cells resulted in enhanced intestinal inflammation and T helper 17 cell (TH17) responses, owing to enhanced inflammasome activation and interleukin (IL)-1β production. This was caused by reduced autophagy in Gcn2−/− intestinal APCs and epithelial cells, leading to increased reactive oxygen species (ROS), a potent activator of inflammasomes. Thus, conditional ablation of Atg5 or Atg7 in intestinal APCs resulted in enhanced ROS and TH17 responses. Furthermore, in vivo blockade of ROS and IL-1β resulted in inhibition of TH17 responses and reduced inflammation in Gcn2−/− mice. Importantly, acute amino acid starvation suppressed intestinal inflammation via a mechanism dependent on GCN2. These results reveal a mechanism that couples amino acid sensing with control of intestinal inflammation via GCN2.

Effects of age, sex, and genotype on high-sensitivity metabolomic profiles in the fruit fly, Drosophila melanogaster.

Researchers have used whole‐genome sequencing and gene expression profiling to identify genes associated with age, in the hope of understanding the underlying mechanisms of senescence. But there is a substantial gap from variation in gene sequences and expression levels to variation in age or life expectancy. In an attempt to bridge this gap, here we describe the effects of age, sex, genotype, and their interactions on high‐sensitivity metabolomic profiles in the fruit fly, Drosophila melanogaster. Among the 6800 features analyzed, we found that over one‐quarter of all metabolites were significantly associated with age, sex, genotype, or their interactions, and multivariate analysis shows that individual metabolomic profiles are highly predictive of these traits. Using a metabolomic equivalent of gene set enrichment analysis, we identified numerous metabolic pathways that were enriched among metabolites associated with age, sex, and genotype, including pathways involving sugar and glycerophospholipid metabolism, neurotransmitters, amino acids, and the carnitine shuttle. Our results suggest that high‐sensitivity metabolomic studies have excellent potential not only to reveal mechanisms that lead to senescence, but also to help us understand differences in patterns of aging among genotypes and between males and females.

Hybrid feature detection and information accumulation using high-resolution LC- MS metabolomics data

Feature detection is a critical step in the preprocessing of liquid chromatography-mass spectrometry (LC-MS) metabolomics data. Currently, the predominant approach is to detect features using noise filters and peak shape models based on the data at hand alone. Databases of known metabolites and historical data contain information that could help boost the sensitivity of feature detection, especially for low-concentration metabolites. However, utilizing such information in targeted feature detection may cause large number of false positives because of the high levels of noise in LC-MS data. With high-resolution mass spectrometry such as liquid chromatograph-Fourier transform mass spectrometry (LC-FTMS), high-confidence matching of peaks to known features is feasible. Here we describe a computational approach that serves two purposes. First it boosts feature detection sensitivity by using a hybrid procedure of both untargeted and targeted peak detection. New algorithms are designed to reduce the chance of false-positives by nonparametric local peak detection and filtering. Second, it can accumulate information on the concentration variation of metabolites over large number of samples, which can help find rare features and/or features with uncommon concentration in future studies. Information can be accumulated on features that are consistently found in real data even before their identities are found. We demonstrate the value of the approach in a proof-of-concept study. The method is implemented as part of the R package apLCMS at http://www.sph.emory.edu/apLCMS/ .

Metabolic Phenotypes of Response to Vaccination in Humans

Herpes zoster (shingles) causes significant morbidity in immune compromised hosts and older adults. Whereas a vaccine is available for prevention of shingles, its efficacy declines with age. To help to understand the mechanisms driving vaccinal responses, we constructed a multiscale, multifactorial response network (MMRN) of immunity in healthy young and older adults immunized with the live attenuated shingles vaccine Zostavax. Vaccination induces robust antigen-specific antibody, plasmablasts, and CD4+ Txa0cells yet limited CD8+ Txa0cell and antiviral responses. The MMRN reveals striking associations between orthogonal datasets, such as transcriptomic and metabolomics signatures, cell populations, and cytokine levels, and identifies immune and metabolic correlates of vaccine immunity. Networks associated with inositol phosphate, glycerophospholipids, and sterol metabolism are tightly coupled with immunity. Critically, the sterol regulatory binding protein 1 and its targets are key integrators of antibody and T follicular cell responses. Our approach is broadly applicable to study human immunity and can help to identify predictors of efficacy as well as mechanisms controlling immunity to vaccination.

Correlation of the lung microbiota with metabolic profiles in bronchoalveolar lavage fluid in HIV infection

BackgroundWhile 16S ribosomal RNA (rRNA) sequencing has been used to characterize the lung’s bacterial microbiota in human immunodeficiency virus (HIV)-infected individuals, taxonomic studies provide limited information on bacterial function and impact on the host. Metabolic profiles can provide functional information on host-microbe interactions in the lungs. We investigated the relationship between the respiratory microbiota and metabolic profiles in the bronchoalveolar lavage fluid of HIV-infected and HIV-uninfected outpatients.ResultsTargeted sequencing of the 16S rRNA gene was used to analyze the bacterial community structure and liquid chromatography-high-resolution mass spectrometry was used to detect features in bronchoalveolar lavage fluid. Global integration of all metabolic features with microbial species was done using sparse partial least squares regression. Thirty-nine HIV-infected subjects and 20 HIV-uninfected controls without acute respiratory symptoms were enrolled. Twelve mass-to-charge ratio (m/z) features from C18 analysis were significantly different between HIV-infected individuals and controls (false discovery rate (FDR)u2009=u20090.2); another 79 features were identified by network analysis. Further metabolite analysis demonstrated that four features were significantly overrepresented in the bronchoalveolar lavage (BAL) fluid of HIV-infected individuals compared to HIV-uninfected, including cystine, two complex carbohydrates, and 3,5-dibromo-l-tyrosine. There were 231u2009m/z features significantly associated with peripheral blood CD4 cell counts identified using sparse partial least squares regression (sPLS) at a variable importance on projection (VIP) threshold of 2. Twenty-five percent of these 91u2009m/z features were associated with various microbial species. Bacteria from families Caulobacteraceae, Staphylococcaceae, Nocardioidaceae, and genus Streptococcus were associated with the greatest number of features. Glycerophospholipid and lineolate pathways correlated with these bacteria.ConclusionsIn bronchoalveolar lavage fluid, specific metabolic profiles correlated with bacterial organisms known to play a role in the pathogenesis of pneumonia in HIV-infected individuals. These findings suggest that microbial communities and their interactions with the host may have functional metabolic impact in the lung.

Detailed Mitochondrial Phenotyping by High Resolution Metabolomics

Mitochondrial phenotype is complex and difficult to define at the level of individual cell types. Newer metabolic profiling methods provide information on dozens of metabolic pathways from a relatively small sample. This pilot study used “top-down” metabolic profiling to determine the spectrum of metabolites present in liver mitochondria. High resolution mass spectral analyses and multivariate statistical tests provided global metabolic information about mitochondria and showed that liver mitochondria possess a significant phenotype based on gender and genotype. The data also show that mitochondria contain a large number of unidentified chemicals.

MetaboAnalyst 4.0: towards more transparent and integrative metabolomics analysis

Abstract We present a new update to MetaboAnalyst (version 4.0) for comprehensive metabolomic data analysis, interpretation, and integration with other omics data. Since the last major update in 2015, MetaboAnalyst has continued to evolve based on user feedback and technological advancements in the field. For this years update, four new key features have been added to MetaboAnalyst 4.0, including: (1) real-time R command tracking and display coupled with the release of a companion MetaboAnalystR package; (2) a MS Peaks to Pathways module for prediction of pathway activity from untargeted mass spectral data using the mummichog algorithm; (3) a Biomarker Meta-analysis module for robust biomarker identification through the combination of multiple metabolomic datasets and (4) a Network Explorer module for integrative analysis of metabolomics, metagenomics, and/or transcriptomics data. The user interface of MetaboAnalyst 4.0 has been reengineered to provide a more modern look and feel, as well as to give more space and flexibility to introduce new functions. The underlying knowledgebases (compound libraries, metabolite sets, and metabolic pathways) have also been updated based on the latest data from the Human Metabolome Database (HMDB). A Docker image of MetaboAnalyst is also available to facilitate download and local installation of MetaboAnalyst. MetaboAnalyst 4.0 is freely available at http://metaboanalyst.ca.