An IFNγ-dependent immune-endocrine circuit lowers blood glucose to potentiate the innate antiviral immune response.

Viral infection makes us feel sick as the immune system alters systemic metabolism to better fight the pathogen. The extent of these changes is relative to the severity of disease. Whether blood glucose is subject to infection-induced modulation is mostly unknown. Here we show that strong, nonlethal infection restricts systemic glucose availability, which promotes the antiviral type I interferon (IFN-I) response. Following viral infection, we find that IFNγ produced by γδ T cells stimulates pancreatic ß cells to increase glucose-induced insulin release. Subsequently, hyperinsulinemia lessens hepatic glucose output. Glucose restriction enhances IFN-I production by curtailing lactate-mediated inhibition of IRF3 and NF-κB signaling. Induced hyperglycemia constrained IFN-I production and increased mortality upon infection. Our findings identify glucose restriction as a physiological mechanism to bring the body into a heightened state of responsiveness to viral pathogens. This immune-endocrine circuit is disrupted in hyperglycemia, possibly explaining why patients with diabetes are more susceptible to viral infection.

The pathophysiology of sepsis and precision-medicine-based immunotherapy.

Sepsis remains a major cause of morbidity and mortality in both low- and high-income countries. Antibiotic therapy and supportive care have significantly improved survival following sepsis in the twentieth century, but further progress has been challenging. Immunotherapy trials for sepsis, mainly aimed at suppressing the immune response, from the 1990s and 2000s, have largely failed, in part owing to unresolved patient heterogeneity in the underlying immune disbalance. The past decade has brought the promise to break this blockade through technological developments based on omics-based technologies and systems medicine that can provide a much larger data space to describe in greater detail the immune endotypes in sepsis. Patient stratification opens new avenues towards precision medicine approaches that aim to apply immunotherapies to sepsis, on the basis of precise biomarkers and molecular mechanisms defining specific immune endotypes. This approach has the potential to lead to the establishment of immunotherapy as a successful pillar in the treatment of sepsis for future generations.

Systematic discovery of TLR signaling components delineates viral-sensing circuits.

Deciphering the signaling networks that underlie normal and disease processes remains a major challenge. Here, we report the discovery of signaling components involved in the Toll-like receptor (TLR) response of immune dendritic cells (DCs), including a previously unkown pathway shared across mammalian antiviral responses. By combining transcriptional profiling, genetic and small-molecule perturbations, and phosphoproteomics, we uncover 35 signaling regulators, including 16 known regulators, involved in TLR signaling. In particular, we find that Polo-like kinases (Plk) 2 and 4 are essential components of antiviral pathways in vitro and in vivo and activate a signaling branch involving a dozen proteins, among which is Tnfaip2, a gene associated with autoimmune diseases but whose role was unknown. Our study illustrates the power of combining systematic measurements and perturbations to elucidate complex signaling circuits and discover potential therapeutic targets.

Pleural macrophages translocate to the lung during infection to promote improved influenza outcomes.

Seasonal influenza results in 3 to 5 million cases of severe disease and 250,000 to 500,000 deaths annually. Macrophages have been implicated in both the resolution and progression of the disease, but the drivers of these outcomes are poorly understood. We probed mouse lung transcriptomic datasets using the Digital Cell Quantifier algorithm to predict immune cell subsets that correlated with mild or severe influenza A virus (IAV) infection outcomes. We identified a unique lung macrophage population that transcriptionally resembled small serosal cavity macrophages and whose presence correlated with mild disease. Until now, the study of serosal macrophage translocation in the context of viral infections has been neglected. Here, we show that pleural macrophages (PMs) migrate from the pleural cavity to the lung after infection with IAV. We found that the depletion of PMs increased morbidity and pulmonary inflammation. There were increased proinflammatory cytokines in the pleural cavity and an influx of neutrophils within the lung. Our results show that PMs are recruited to the lung during IAV infection and contribute to recovery from influenza. This study expands our knowledge of PM plasticity and identifies a source of lung macrophages independent of monocyte recruitment and local proliferation.

A physical and regulatory map of host-influenza interactions reveals pathways in H1N1 infection.

During the course of a viral infection, viral proteins interact with an array of host proteins and pathways. Here, we present a systematic strategy to elucidate the dynamic interactions between H1N1 influenza and its human host. A combination of yeast two-hybrid analysis and genome-wide expression profiling implicated hundreds of human factors in mediating viral-host interactions. These factors were then examined functionally through depletion analyses in primary lung cells. The resulting data point to potential roles for some unanticipated host and viral proteins in viral infection and the host response, including a network of RNA-binding proteins, components of WNT signaling, and viral polymerase subunits. This multilayered approach provides a comprehensive and unbiased physical and regulatory model of influenza-host interactions and demonstrates a general strategy for uncovering complex host-pathogen relationships.

Cell composition analysis of bulk genomics using single-cell data.

Single-cell RNA sequencing (scRNA-seq) is a rich resource of cellular heterogeneity, opening new avenues in the study of complex tissues. We introduce Cell Population Mapping (CPM), a deconvolution algorithm in which reference scRNA-seq profiles are leveraged to infer the composition of cell types and states from bulk transcriptome data ('scBio' CRAN R-package). Analysis of individual variations in lungs of influenza-virus-infected mice reveals that the relationship between cell abundance and clinical symptoms is a cell-state-specific property that varies gradually along the continuum of cell-activation states. The gradual change is confirmed in subsequent experiments and is further explained by a mathematical model in which clinical outcomes relate to cell-state dynamics along the activation process. Our results demonstrate the power of CPM in reconstructing the continuous spectrum of cell states within heterogeneous tissues.

Inferring cellular heterogeneity of associations from single cell genomics.

MOTIVATION: Cell-to-cell variation has uncovered associations between cellular phenotypes. However, it remains challenging to address the cellular diversity of such associations. RESULTS: Here, we do not rely on the conventional assumption that the same association holds throughout the entire cell population. Instead, we assume that associations may exist in a certain subset of the cells. We developed CEllular Niche Association (CENA) to reliably predict pairwise associations together with the cell subsets in which the associations are detected. CENA does not rely on predefined subsets but only requires that the cells of each predicted subset would share a certain characteristic state. CENA may therefore reveal dynamic modulation of dependencies along cellular trajectories of temporally evolving states. Using simulated data, we show the advantage of CENA over existing methods and its scalability to a large number of cells. Application of CENA to real biological data demonstrates dynamic changes in associations that would be otherwise masked. AVAILABILITY AND IMPLEMENTATION: CENA is available as an R package at Github: https://github.com/mayalevy/CENA and is accompanied by a complete set of documentations and instructions. CONTACT: iritgv@gmail.com. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Anemia measurements to distinguish between viral and bacterial infections in the emergency department.

The clinical diagnosis of acute infections in the emergency department is a challenging task due to the similarity in symptom presentation between virally and bacterially infected individuals, while the use of routine laboratory tests for pathogen identification is often time-consuming and may contain contaminants. We investigated the ability of various anemia-related parameters, including hemoglobin, red cell distribution width (RDW), and iron, to differentiate between viral and bacterial infection in a retrospective study of 3883 patients admitted to the emergency department with a confirmed viral (n = 1238) or bacterial (n = 2645) infection based on either laboratory tests or microbiological cultures. The ratio between hemoglobin to RDW was found to be significant in distinguishing between virally and bacterially infected patients and outperformed other anemia measurements. Moreover, the predictive value of the ratio was high even in patients presenting with low C-reactive protein values (< 21 mg/L). We followed the dynamics of hemoglobin, RDW, and the ratio between them up to 72 h post emergency department admission, and observed a consistent discrepancy between virally and bacterially infected patients over time. Additional analysis demonstrated higher levels of ferritin and lower levels of iron in bacterially infected compared with virally infected patients. The anemia measurements were associated with length of hospital stay, where all higher levels, except for RDW, corresponded to a shorter hospitalization period. We highlighted the importance of various anemia measurements as an additional host-biomarker to discern virally from bacterially infected patients.

ImmQuant: a user-friendly tool for inferring immune cell-type composition from gene-expression data.

: The composition of immune-cell subsets is key to the understanding of major diseases and pathologies. Computational deconvolution methods enable researchers to investigate immune cell quantities in complex tissues based on transcriptome data. Here we present ImmQuant, a software tool allowing immunologists to upload transcription profiles of multiple tissue samples, apply deconvolution methodology to predict differences in cell-type quantities between the samples, and then inspect the inferred cell-type alterations using convenient visualization tools. ImmQuant builds on the DCQ deconvolution algorithm and allows a user-friendly utilization of this method by non-bioinformatician researchers. Specifically, it enables investigation of hundreds of immune cell subsets in mouse tissues, as well as a few dozen cell types in human samples. AVAILABILITY AND IMPLEMENTATION: ImmQuant is available for download at http://csgi.tau.ac.il/ImmQuant/ CONTACT: iritgv@post.tau.ac.ilSupplementary information: Supplementary data are available at Bioinformatics online.

Exploiting Gene-Expression Deconvolution to Probe the Genetics of the Immune System.

Sequence variation can affect the physiological state of the immune system. Major experimental efforts targeted at understanding the genetic control of the abundance of immune cell subpopulations. However, these studies are typically focused on a limited number of immune cell types, mainly due to the use of relatively low throughput cell-sorting technologies. Here we present an algorithm that can reveal the genetic basis of inter-individual variation in the abundance of immune cell types using only gene expression and genotyping measurements as input. Our algorithm predicts the abundance of immune cell subpopulations based on the RNA levels of informative marker genes within a complex tissue, and then provides the genetic control on these predicted immune traits as output. A key feature of the approach is the integration of predictions from various sets of marker genes and refinement of these sets to avoid spurious signals. Our evaluation of both synthetic and real biological data shows the significant benefits of the new approach. Our method, VoCAL, is implemented in the freely available R package ComICS.

A statistical framework for revealing signaling pathways perturbed by DNA variants.

Much of the inter-individual variation in gene expression is triggered via perturbations of signaling networks by DNA variants. We present a novel probabilistic approach for identifying the particular pathways by which DNA variants perturb the signaling network. Our procedure, called PINE, relies on a systematic integration of established biological knowledge of signaling networks with data on transcriptional responses to various experimental conditions. Unlike previous approaches, PINE provides statistical aspects that are critical for prioritizing hypotheses for followup experiments. Using simulated data, we show that higher accuracy is attained with PINE than with existing methods. We used PINE to analyze transcriptional responses of immune dendritic cells to several pathogenic stimulations. PINE identified statistically significant genetic perturbations in the pathogen-sensing signaling network, suggesting previously uncharacterized regulatory mechanisms for functional DNA variants.

Mapping liver fat female-dependent quantitative trait loci in collaborative cross mice.

Non-alcoholic fatty liver disease (NAFLD) is the most common cause of chronic liver disease in the western world, with spectrum from simple steatosis to non-alcoholic steatohepatitis, which can progress to cirrhosis. NAFLD developments are known to be affected by host genetic background. Herein we emphasize the power of collaborative cross (CC) mouse for dissecting this complex trait and revealing quantitative trait loci (QTL) controlling hepatic fat accumulation in mice. 168 female and 338 male mice from 24 and 37 CC lines, respectively, of 18-20 weeks old, maintained on standard rodent diet, since weaning. Hepatic fat content was assessed, using dual DEXA scan in the liver. Using the available high-density genotype markers of the CC line, QTL mapping associated with percentage liver fat accumulation was performed. Our results revealed significant fatty liver accumulation QTL that were specifically, mapped in females. Two significant QTLs on chromosomes 17 and 18, with genomic intervals 3 and 2 Mb, respectively, were mapped. A third QTL, with a less significant P value, was mapped to chromosome 4, with genomic interval of 2 Mb. These QTLs were named Flal1-Flal3, referring to Fatty Liver Accumulation Locus 1-3, for the QTLs on chromosomes 17, 18, and 4, respectively. Unfortunately, no QTL was mapped with males. Searching the mouse genome database suggested several candidate genes involved in hepatic fat accumulation. Our results show that susceptibility to hepatic fat accumulations is a complex trait, controlled by multiple genetic factors in female mice, but not in male.

CoD: inferring immune-cell quantities related to disease states.

MOTIVATION: The immune system comprises a complex network of genes, cells and tissues, coordinated through signaling pathways and cell-cell communications. However, the orchestrated role of the multiple immunological components in disease is still poorly understood. Classifications based on gene-expression data have revealed immune-related signaling pathways in various diseases, but how such pathways describe the immune cellular physiology remains largely unknown. RESULTS: We identify alterations in cell quantities discriminating between disease states using ' Cell type of Disease' (CoD), a classification-based approach that relies on computational immune-cell decomposition in gene-expression datasets. CoD attains significantly higher accuracy than alternative state-of-the-art methods. Our approach is shown to recapitulate and extend previous knowledge acquired with experimental cell-quantification technologies. CONCLUSIONS: The results suggest that CoD can reveal disease-relevant cell types in an unbiased manner, potentially heralding improved diagnostics and treatment. AVAILABILITY AND IMPLEMENTATION: The software described in this article is available at http://www.csgi.tau.ac.il/CoD/.

Proteomics-level analysis of myelin formation and regeneration in a mouse model for Vanishing White Matter disease.

Vanishing white matter (VWM) is a recessive neurodegenerative disease caused by mutations in translation initiation factor eIF2B and leading to progressive brain myelin deterioration, secondary axonal damage, and death in early adolescence. Eif2b5(R132H/R132H) mice exhibit delayed developmental myelination, mild early neurodegeneration and a robust remyelination defect in response to cuprizone-induced demyelination. In the current study we used Eif2b5(R132H/R132H) mice for mass-spectrometry analyses, to follow the changes in brain protein abundance in normal- versus cuprizone-diet fed mice during the remyelination recovery phase. Analysis of proteome profiles suggested that dysregulation of mitochondrial functions, altered proteasomal activity and impaired balance between protein synthesis and degradation play a role in VWM pathology. Consistent with these findings, we detected elevated levels of reactive oxygen species in mutant-derived primary fibroblasts and reduced 20S proteasome activity in mutant brain homogenates. These observations highlight the importance of tight translational control to precise coordination of processes involved in myelin formation and regeneration and point at cellular functions that may contribute to VWM pathology. Eif2b5(R132H/R132H) mouse model for vanishing white matter (VWM) disease was used for mass spectrometry of brain proteins at two time points under normal conditions and along recovery from cuprizone-induced demyelination. Comparisons of proteome profiles revealed the importance of mitochondrial function and tight coordination between protein synthesis and degradation to myelination formation and regeneration, pointing at cellular functions that contribute to VWM pathology.

Digital cell quantification identifies global immune cell dynamics during influenza infection.

Hundreds of immune cell types work in coordination to maintain tissue homeostasis. Upon infection, dramatic changes occur with the localization, migration, and proliferation of the immune cells to first alert the body of the danger, confine it to limit spreading, and finally extinguish the threat and bring the tissue back to homeostasis. Since current technologies can follow the dynamics of only a limited number of cell types, we have yet to grasp the full complexity of global in vivo cell dynamics in normal developmental processes and disease. Here, we devise a computational method, digital cell quantification (DCQ), which combines genome-wide gene expression data with an immune cell compendium to infer in vivo changes in the quantities of 213 immune cell subpopulations. DCQ was applied to study global immune cell dynamics in mice lungs at ten time points during 7 days of flu infection. We find dramatic changes in quantities of 70 immune cell types, including various innate, adaptive, and progenitor immune cells. We focus on the previously unreported dynamics of four immune dendritic cell subtypes and suggest a specific role for CD103(+) CD11b(-) DCs in early stages of disease and CD8(+) pDC in late stages of flu infection.

Dissecting dynamic genetic variation that controls temporal gene response in yeast.

Inter-individual variation in regulatory circuits controlling gene expression is a powerful source of functional information. The study of associations among genetic variants and gene expression provides important insights about cell circuitry but cannot specify whether and when potential variants dynamically alter their genetic effect during the course of response. Here we develop a computational procedure that captures temporal changes in genetic effects, and apply it to analyze transcription during inhibition of the TOR signaling pathway in segregating yeast cells. We found a high-order coordination of gene modules: sets of genes co-associated with the same genetic variant and sharing a common temporal genetic effect pattern. The temporal genetic effects of some modules represented a single state-transitioning pattern; for example, at 10-30 minutes following stimulation, genetic effects in the phosphate utilization module attained a characteristic transition to a new steady state. In contrast, another module showed an impulse pattern of genetic effects; for example, in the poor nitrogen sources utilization module, a spike up of a genetic effect at 10-20 minutes following stimulation reflected inter-individual variation in the timing (rather than magnitude) of response. Our analysis suggests that the same mechanism typically leads to both inter-individual variation and the temporal genetic effect pattern in a module. Our methodology provides a quantitative genetic approach to studying the molecular mechanisms that shape dynamic changes in transcriptional responses.

The host transcriptional response to superinfection by influenza A virus and Streptococcus pneumoniae.

Secondary bacterial challenges during influenza virus infection "superinfection") cause excessive mortality and hospitalization. Here, we present a longitudinal study of bulk gene expression changes in murine lungs during superinfection, with an initial influenza A virus infection and a subsequent Streptococcus pneumoniae infection. In addition to the well-characterized impairment of the host response, we identified superinfection-specific alterations in the global transcriptional program that are linked to the host's ability to resist the pathogens. Particularly, whereas superinfected mice manifested an excessive rapid induction of the resistance-to-infection program, there was a substantial tissue-level rewiring of this program: upon superinfection, interferon-regulated genes were switched from positive to negative correlations with the host's resistance state, whereas genes of fatty acid metabolism switched from negative to positive correlations with resistance states. Thus, the transcriptional resistance state in superinfection is reprogrammed toward repressed interferon signaling and induced fatty acid metabolism. Our findings suggest new insights into a tissue-level remodeling of the host defense upon superinfection, providing promising targets for future therapeutic interventions. IMPORTANCE: Secondary bacterial infections are the most frequent complications during influenza A virus (IAV) pandemic outbreaks, contributing to excessive morbidity and mortality in the human population. Most IAV-related deaths are attributed to Streptococcus pneumoniae (SP) infections, which usually begin within the first week of IAV infection in the respiratory tracts. Here, we focused on longitudinal transcriptional responses during a superinfection model consisting of an SP infection that follows an initial IAV infection, comparing superinfection to an IAV-only infection, an SP-only infection, and control treatments. Our longitudinal data allowed a fine analysis of gene expression changes during superinfection. For instance, we found that superinfected mice exhibited rapid gene expression induction or reduction within the first 12 h after encountering the second pathogen. Cell proliferation and immune response activation processes were upregulated, while endothelial processes, vasculogenesis, and angiogenesis were downregulated, providing promising targets for future therapeutic interventions. We further analyzed the longitudinal transcriptional responses in the context of a previously defined spectrum of the host's resistance state, revealing superinfection-specific reprogramming of resistance states, such as reprogramming of fatty acid metabolism and interferon signaling. The reprogrammed functions are compelling new targets for switching the pathogenic superinfection state into a single-infection state.

Identifying genetic susceptibility to Aspergillus fumigatus infection using collaborative cross mice and RNA-Seq approach.

BACKGROUND: Aspergillus fumigatus (Af) is one of the most ubiquitous fungi and its infection potency is suggested to be strongly controlled by the host genetic background. The aim of this study was to search for candidate genes associated with host susceptibility to Aspergillus fumigatus (Af) using an RNAseq approach in CC lines and hepatic gene expression. METHODS: We studied 31 male mice from 25 CC lines at 8 weeks old; the mice were infected with Af. Liver tissues were extracted from these mice 5 days post-infection, and next-generation RNA-sequencing (RNAseq) was performed. The GENE-E analysis platform was used to generate a clustered heat map matrix. RESULTS: Significant variation in body weight changes between CC lines was observed. Hepatic gene expression revealed 12 top prioritized candidate genes differentially expressed in resistant versus susceptible mice based on body weight changes. Interestingly, three candidate genes are located within genomic intervals of the previously mapped quantitative trait loci (QTL), including Gm16270 and Stox1 on chromosome 10 and Gm11033 on chromosome 8. CONCLUSIONS: Our findings emphasize the CC mouse model's power in fine mapping the genetic components underlying susceptibility towards Af. As a next step, eQTL analysis will be performed for our RNA-Seq data. Suggested candidate genes from our study will be further assessed with a human cohort with aspergillosis.

Stratification of Atherosclerosis based on Plasma Metabolic States.

CONTEXT: Atherosclerosis is a dominant cause of cardiovascular disease (CVD), including myocardial infarction and stroke. OBJECTIVE: To investigate metabolic states that are associated with the development of atherosclerosis. METHODS: Cross-sectional cohort study at a university hospital in the Netherlands. A total of 302 adult subjects with a body mass index (BMI) ≥ 27 kg/m2 were included. We integrated plasma metabolomics with clinical metadata to quantify the "atherogenic state" of each individual, providing a continuous spectrum of atherogenic states that ranges between nonatherogenic states to highly atherogenic states. RESULTS: Analysis of groups of individuals with different clinical conditions-such as metabolically healthy individuals with obesity, and individuals with metabolic syndrome-confirmed the generalizability of this spectrum; revealed a wide variation of atherogenic states within each condition; and allowed identification of metabolites that are associated with the atherogenic state regardless of the particular condition, such as gamma-glutamyl-glutamic acid and homovanillic acid sulfate. The analysis further highlighted metabolic pathways such as catabolism of phenylalanine and tyrosine and biosynthesis of estrogens and phenylpropanoids. Using validation cohorts, we confirmed variation in atherogenic states in healthy subjects (before atherosclerosis plaques become visible), and showed that metabolites associated with the atherogenic state were also associated with future CVD. CONCLUSION: Our results provide a global view of atherosclerosis risk states using plasma metabolomics.

Understanding gene sequence variation in the context of transcription regulation in yeast.

DNA sequence polymorphism in a regulatory protein can have a widespread transcriptional effect. Here we present a computational approach for analyzing modules of genes with a common regulation that are affected by specific DNA polymorphisms. We identify such regulatory-linkage modules by integrating genotypic and expression data for individuals in a segregating population with complementary expression data of strains mutated in a variety of regulatory proteins. Our procedure searches simultaneously for groups of co-expressed genes, for their common underlying linkage interval, and for their shared regulatory proteins. We applied the method to a cross between laboratory and wild strains of S. cerevisiae, demonstrating its ability to correctly suggest modules and to outperform extant approaches. Our results suggest that middle sporulation genes are under the control of polymorphism in the sporulation-specific tertiary complex Sum1p/Rfm1p/Hst1p. In another example, our analysis reveals novel inter-relations between Swi3 and two mitochondrial inner membrane proteins underlying variation in a module of aerobic cellular respiration genes. Overall, our findings demonstrate that this approach provides a useful framework for the systematic mapping of quantitative trait loci and their role in gene expression variation.