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

Intersection of Immunity, Metabolism, and Muscle Regeneration in an Autoimmune-Prone MRL Mouse Model

Due to the limited capacity of mammals to regenerate complex tissues, researchers have worked to understand the mechanisms of tissue regeneration in organisms that maintain that capacity. One example is the MRL/MpJ mouse strain with unique regenerative capacity in ear pinnae that is absent from other strains, such as the common C57BL/6 strain. The MRL/MpJ mouse has also been associated with an autoimmune phenotype even in the absence of the mutant Fas gene described in its parent strain MRL/lpr. Due to these findings, the differences between the responses of MRL/MpJ versus C57BL/6 strain are evaluated in volumetric muscle injury and subsequent material implantation. One salient feature of the MRL/MpJ response to injury is robust adipogenesis within the muscle. This is associated with a decrease in M2-like polarization in response to biologically derived extracellular matrix scaffolds. In pro-fibrotic materials, such as polyethylene, there are fewer foreign body giant cells in the MRL/MpJ mice. As there are reports of both positive and negative influences of adipose tissue and adipogenesis on wound healing, this model can provide an important lens to investigate the interplay between stem cells, adipose tissue, and immune responses in trauma and material implantation.

Ectopic adipogenesis in response to injury and material implantation in an autoimmune mouse model

Due to the limited capacity of mammals to regenerate complex tissues, researchers have worked to understand the mechanisms of tissue regeneration in organisms that maintain that capacity. One example is the MRL/MpJ mouse strain with unique regenerative capacity in ear pinnae that is absent from other strains, such as the common C57BL/6 strain. The MRL/MpJ mouse has also been associated with an autoimmune phenotype even in the absence of the mutant Fas gene described in its parent strain MRL/lpr. Due to these findings, we evaluated the differences between the responses of MRL/MpJ versus C57BL/6 strain in traumatic muscle injury and subsequent material implantation. One salient feature of the MRL/MpJ response to injury was a robust adipogenesis within the muscle. This was associated with a decrease in M2-like polarization in response to biologically derived extracellular matrix scaffolds. In pro-fibrotic materials, such as polyethylene, there were fewer foreign body giant cells in the MRL/MpJ mice. As there are reports of both positive and negative influences of adipose tissue and adipogenesis on wound healing, this model could provide an important lens to investigate the interplay between stem cells, adipose tissue, and immune responses in trauma and materials implantation.

Characterization of virus-mediated autoimmunity and the consequences for pathological process in patients with systemic lupus erythematosus

INTRODUCTION/OBJECTIVES

This study aimed to identify differentially expressed genes (DEGs) of systemic lupus erythematosus (SLE) using gene expression-based computational methodologies to analyze disease-immune interactions, which affect the development and progression of SLE.

METHOD

Twenty-six patients with SLE and 46 healthy controls were selected from the Gene Expression Omnibus (GEO) database. The significantly enriched immune and virus-related gene lists were computed and visualized by using the DEGs from the gene set enrichment analysis (GSEA). Quantification of 38 immune cells was performed in determining the impact of immune cells on the virus mediated immunity in SLE by using ImmQuant algorithm.

RESULTS

Thirty-nine upregulated and 57 downregulated were identified in SLE patient compared to the healthy controls. Upregulated genes were significantly implicated in Gene Ontology gene sets as cytokine mediated signaling, secretion, and exocytosis in immune response pathways in 26 female SLE patients. In addition, these genes were enriched in hepatitis C, influenza A, measles, Epstein-Barr virus, and herpes simplex virus 1 infection in Kyoto Encyclopedia of Genes and Genomes pathways. Especially, FCGR1A, IRF7, OAS2, CAMP, MX1, OAS3, OAS1, DEFA3, ISG15, and RSAD2 were involved in virus mediated SLE mechanism, and the expression for OAS1, OAS2, and IRF7 was closely associated with the quantities of colony forming unit-monocyte and colony forming unit-granulocyte.

CONCLUSIONS

Identifying virus-mediated SLE genes and quantifies of immune cells were used to understand the pathological process and perform early diagnosis of female SLE, and will lead to clinical tools for treating SLE in patients. Key Points • Using gene expression-based computational methodologies, the 57 immune and viral genes were significantly upregulated in 26 SLE patients. • The identified three key  viral genes such as OAS1, OAS2, and IF7 were closely associated with colony-forming unit-monocytes and colony-forming unit-granulocytes, which affect the virus mediated immunity in SLE. • The viral genes and quantifies of immune cells are useful in understanding pathogenesis of SLE, and this will provide clinical strategies of potential treatment choices in SLE patients.

Analysis of shared ceRNA networks and related-hub genes in rats with primary and secondary photoreceptor degeneration

Introduction

Photoreceptor degenerative diseases are characterized by the progressive death of photoreceptor cells, resulting in irreversible visual impairment. However, the role of competing endogenous RNA (ceRNA) in photoreceptor degeneration is unclear. We aimed to explore the shared ceRNA regulation network and potential molecular mechanisms between primary and secondary photoreceptor degenerations.

Methods

We established animal models for both types of photoreceptor degenerations and conducted retina RNA sequencing to identify shared differentially expressed long non-coding RNAs (lncRNAs), microRNAs (miRNAs), and messenger RNAs (mRNAs). Using ceRNA regulatory principles, we constructed a shared ceRNA network and performed function enrichment and protein-protein interaction (PPI) analyses to identify hub genes and key pathways. Immune cell infiltration and drug-gene interaction analyses were conducted, and hub gene expression was validated by quantitative real-time polymerase chain reaction (qRT-PCR).

Results

We identified 37 shared differentially expressed lncRNAs, 34 miRNAs, and 247 mRNAs and constructed a ceRNA network consisting of 3 lncRNAs, 5 miRNAs, and 109 mRNAs. Furthermore, we examined 109 common differentially expressed genes (DEGs) through functional annotation, PPI analysis, and regulatory network analysis. We discovered that these diseases shared the complement and coagulation cascades pathway. Eight hub genes were identified and enriched in the immune system process. Immune infiltration analysis revealed increased T cells and decreased B cells in both photoreceptor degenerations. The expression of hub genes was closely associated with the quantities of immune cell types. Additionally, we identified 7 immune therapeutical drugs that target the hub genes.

Discussion

Our findings provide new insights and directions for understanding the common mechanisms underlying the development of photoreceptor degeneration. The hub genes and related ceRNA networks we identified may offer new perspectives for elucidating the mechanisms and hold promise for the development of innovative treatment strategies.

Measuring the transcriptome-wide effects of aging on murine adipocytes using RNAseq

Adipose tissue plays a central role in age-related diseases. While RNAseq protocols exist for many tissues, few data have been generated with this technology to explore gene expression in adipocytes, particularly during aging. Here, we present a protocol to analyze the transcriptional changes that occur in adipose tissue during normal and accelerated aging in mouse models. We describe steps for genotyping, diet control, euthanasia, and dissection. We then detail RNA purification and genome-wide data generation and analysis. For complete details on the use and execution of this protocol, please refer to De Cauwer et al. (2022) iScience. Sep 16;25(10):105149.

Dicer1 deficient mice exhibit premature aging and metabolic perturbations in adipocytes

Age-related diseases are major concern in developed countries. To avoid disabilities that accompany increased lifespan, pharmaceutical approaches are considered. Therefore, appropriate animal models are required for a better understanding of aging processes and potential in vivo assays to evaluate the impact of molecules that may delay the occurrence of age-related diseases. Few mouse models exhibiting pathological aging exist, but currently, none of them reproducibly mimics human diseases like osteoporosis, cognitive dysfunctions or sarcopenia that can be seen in some, but not all, elders. Here, we describe the premature aging phenotypes of Dicer-deficient mature animals, which exhibit an overall deterioration of many organs and tissues (skin, heart, and adipose tissue) ultimately leading to a significant reduction of their lifespan. Molecular characterization of transcriptional responses focused on the adipose tissue suggested that both canonical and non-canonical functions of DICER are involved in this process and highlight potential actionable pathways to revert it.

•

Dicer1-deficient mice have a reduced lifespan with early age-related symptoms

•

Mutant mice are resistant to high fat diet-induced disorders

•

Myokines FGF21 and GDF15 are likely key regulators of adipocytes metabolism

CX3CR1-Expressing Myeloid Cells Regulate Host-Helminth Interaction and Lung Inflammation

Many helminth life cycles, including hookworm, involve a mandatory lung phase, where myeloid and granulocyte subsets interact with the helminth and respond to infection-induced lung injury. To evaluate these innate subsets in Nippostrongylus brasiliensis infection, reporter mice for myeloid cells (CX3CR1^GFP ) and granulocytes (PGRP^dsRED ) are employed. Nippostrongylus infection induces lung infiltration of reporter cells, including CX3CR1⁺ myeloid cells and PGRP⁺ eosinophils. Strikingly, CX3CR1^GFP/GFP mice, which are deficient in CX3CR1, are protected from Nippostrongylus infection with reduced weight loss, lung leukocyte infiltration, and worm burden compared to CX3CR1^+/+ mice. This protective effect is specific for CX3CR1 as CCR2-deficient mice do not exhibit reduced worm burdens. Nippostrongylus co-culture with lung Ly6C⁺ monocytes or CD11c⁺ cells demonstrates that CX3CR1^GFP/GFP monocytes secrete more pro-inflammatory cytokines and actively bind the parasites causing reduced motility. RNA sequencing of Ly6C⁺ or CD11c⁺ cells shows Nippostrongylus-induced gene expression changes, particularly in monocytes, associated with inflammation, chemotaxis, and extracellular matrix remodeling pathways. Analysis reveals cytotoxic and adhesion molecules as potential effectors against the parasite, such as Gzma and Gzmb, which are elevated in CX3CR1^GFP/GFP monocytes. These studies validate a dual innate cell reporter for lung helminth infection and demonstrate that CX3CR1 impairs monocyte-helminth interaction.

Multidimensional analysis of the host response reveals prognostic and pathogen-driven immune subtypes among adults with sepsis in Uganda

Background

The global burden of sepsis is concentrated in sub-Saharan Africa, where severe infections disproportionately affect young, HIV-infected adults and high-burden pathogens are unique. In this context, poor understanding of sepsis immunopathology represents a crucial barrier to development of locally-effective treatment strategies. We sought to determine inter-individual immunologic heterogeneity among adults hospitalized with sepsis in a sub-Saharan African setting, and characterize associations between immune subtypes, infecting pathogens, and clinical outcomes.

Methods

Among a prospective observational cohort of 288 adults hospitalized with suspected sepsis in Uganda, we applied machine learning methods to 14 soluble host immune mediators, reflective of key domains of sepsis immunopathology (innate and adaptive immune activation, endothelial dysfunction, fibrinolysis), to identify immune subtypes in randomly-split discovery (N = 201) and internal validation (N = 87) sub-cohorts. In parallel, we applied similar methods to whole-blood RNA-sequencing data from a consecutive subset of patients (N = 128) to identify transcriptional subtypes, which we characterized using biological pathway and immune cell-type deconvolution analyses.

Results

Unsupervised clustering consistently identified two immune subtypes defined by differential activation of pro-inflammatory innate and adaptive immune pathways, with transcriptional evidence of concomitant CD56(-)/CD16( +) NK-cell expansion, T-cell exhaustion, and oxidative-stress and hypoxia-induced metabolic and cell-cycle reprogramming in the hyperinflammatory subtype. Immune subtypes defined by greater pro-inflammatory immune activation, T-cell exhaustion, and metabolic reprogramming were consistently associated with a high-prevalence of severe and often disseminated HIV-associated tuberculosis, as well as more extensive organ dysfunction, worse functional outcomes, and higher 30-day mortality.

Conclusions

Our results highlight unique host- and pathogen-driven features of sepsis immunopathology in sub-Saharan Africa, including the importance of severe HIV-associated tuberculosis, and reinforce the need to develop more biologically-informed treatment strategies in the region, particularly those incorporating immunomodulation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13054-022-03907-3.

Corticosteroid treatment in COVID-19 modulates host inflammatory responses and transcriptional signatures of immune dysregulation

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of coronavirus disease-2019 (COVID-19), a respiratory disease that varies in severity from mild to severe/fatal. Several risk factors for severe disease have been identified, notably age, male sex, and pre-existing conditions such as diabetes, obesity, and hypertension. Several advancements in clinical care have been achieved over the past year, including the use of corticosteroids (e.g., corticosteroids) and other immune-modulatory treatments that have now become standard of care for patients with acute severe COVID-19. While the understanding of the mechanisms that underlie increased disease severity with age has improved over the past few months, it remains incomplete. Furthermore, the molecular impact of corticosteroid treatment on host response to acute SARS-CoV-2 infection has not been investigated. In this study, a cross-sectional and longitudinal analysis of Ab, soluble immune mediators, and transcriptional responses in young (65 ≤ years) and aged (≥ 65 years) diabetic males with obesity hospitalized with acute severe COVID-19 was conducted. Additionally, the transcriptional profiles in samples obtained before and after corticosteroids became standard of care were compared. The analysis indicates that severe COVID-19 is characterized by robust Ab responses, heightened systemic inflammation, increased expression of genes related to inflammatory and pro-apoptotic processes, and reduced expression of those important for adaptive immunity regardless of age. In contrast, COVID-19 patients receiving steroids did not show high levels of systemic immune mediators and lacked transcriptional indicators of heightened inflammatory and apoptotic responses. Overall, these data suggest that inflammation and cell death are key drivers of severe COVID-19 pathogenesis in the absence of corticosteroid therapy.

Systematic evaluation of transcriptomics-based deconvolution methods and references using thousands of clinical samples

Estimating cell type composition of blood and tissue samples is a biological challenge relevant in both laboratory studies and clinical care. In recent years, a number of computational tools have been developed to estimate cell type abundance using gene expression data. Although these tools use a variety of approaches, they all leverage expression profiles from purified cell types to evaluate the cell type composition within samples. In this study, we compare 12 cell type quantification tools and evaluate their performance while using each of 10 separate reference profiles. Specifically, we have run each tool on over 4000 samples with known cell type proportions, spanning both immune and stromal cell types. A total of 12 of these represent in vitro synthetic mixtures and 300 represent in silico synthetic mixtures prepared using single-cell data. A final 3728 clinical samples have been collected from the Framingham cohort, for which cell populations have been quantified using electrical impedance cell counting. When tools are applied to the Framingham dataset, the tool Estimating the Proportions of Immune and Cancer cells (EPIC) produces the highest correlation, whereas Gene Expression Deconvolution Interactive Tool (GEDIT) produces the lowest error. The best tool for other datasets is varied, but CIBERSORT and GEDIT most consistently produce accurate results. We find that optimal reference depends on the tool used, and report suggested references to be used with each tool. Most tools return results within minutes, but on large datasets runtimes for CIBERSORT can exceed hours or even days. We conclude that deconvolution methods are capable of returning high-quality results, but that proper reference selection is critical.

Transcriptional Analysis of Infection With Early or Late Isolates From the 2013-2016 West Africa Ebola Virus Epidemic Does Not Suggest Attenuated Pathogenicity as a Result of Genetic Variation

The 2013-2016 West Africa Ebola virus (EBOV) epidemic caused by the EBOV-Makona isolate is the largest and longest recorded to date. It incurred over 28,000 infections and ∼11,000 deaths. Early in this epidemic, several mutations in viral glycoprotein (A82V), nucleoprotein (R111C), and polymerase L (D759G) emerged and stabilized. In vitro studies of these new EBOV-Makona isolates showed enhanced fitness and viral replication capacity. However, in vivo studies in mice and rhesus macaques did not provide any evidence of enhanced viral fitness or shedding. Infection with late isolates carrying or early isolates lacking (early) these mutations resulted in uniformly lethal disease in nonhuman primates (NHPs), albeit with slightly delayed kinetics with late isolates. The recent report of a possible reemergence of EBOV from a persistent infection in a survivor of the epidemic highlights the urgency for understanding the impact of genetic variation on EBOV pathogenesis. However, potential molecular differences in host responses remain unknown. To address this gap in knowledge, we conducted the first comparative analysis of the host responses to lethal infection with EBOV-Mayinga and EBOV-Makona isolates using bivariate, longitudinal, regression, and discrimination transcriptomic analyses. Our analysis shows a conserved core of differentially expressed genes (DEGs) involved in antiviral defense, immune cell activation, and inflammatory processes in response to EBOV-Makona and EBOV-Mayinga infections. Additionally, EBOV-Makona and EBOV-Mayinga infections could be discriminated based on the expression pattern of a small subset of genes. Transcriptional responses to EBOV-Makona isolates that emerged later during the epidemic, specifically those from Mali and Liberia, lacked signatures of profound lymphopenia and excessive inflammation seen following infection with EBOV-Mayinga and early EBOV-Makona isolate C07. Overall, these findings provide novel insight into the mechanisms underlying the lower case fatality rate (CFR) observed with EBOV-Makona compared to EBOV-Mayinga.

Growth faltering regardless of chronic diarrhea is associated with mucosal immune dysfunction and microbial dysbiosis in the gut lumen

Despite the impact of childhood diarrhea on morbidity and mortality, our understanding of its sequelae has been significantly hampered by the lack of studies that examine samples across the entire intestinal tract. Infant rhesus macaques are naturally susceptible to human enteric pathogens and recapitulate the hallmarks of diarrheal disease such as intestinal inflammation and growth faltering. Here, we examined intestinal biopsies, lamina propria leukocytes, luminal contents, and fecal samples from healthy infants and those experiencing growth faltering with distant acute or chronic active diarrhea. We show that growth faltering in the presence or absence of active diarrhea is associated with a heightened systemic and mucosal pro-inflammatory state centered in the colon. Moreover, polyclonal stimulation of colonic lamina propria leukocytes resulted in a dampened cytokine response, indicative of immune exhaustion. We also detected a functional and taxonomic shift in the luminal microbiome across multiple gut sites including the migration of Streptococcus and Prevotella species between the small and large intestine, suggesting a decompartmentalization of gut microbial communities. Our studies provide valuable insight into the outcomes of diarrheal diseases and growth faltering not attainable in humans and lays the groundwork to test interventions in a controlled and reproducible setting.

Transcriptomic landscape of skin lesions in cutaneous leishmaniasis reveals a strong CD8+ T cell immunosenescence signature linked to immunopathology

The severity of lesions that develop in patients infected by Leishmania braziliensis is mainly associated with a highly cytotoxic and inflammatory cutaneous environment. Recently, we demonstrated that senescent T and NK cells play a role in the establishment and maintenance of this tissue inflammation. Here, we extended those findings using transcriptomic analyses that demonstrate a strong co-induction of senescence and pro-inflammatory gene signatures in cutaneous leishmaniasis (CL) lesions. The senescence-associated signature was characterized by marked expression of key genes such as ATM, Sestrin 2, p16, p21 and p38. The cell type identification from deconvolution of bulk sequencing data showed that the senescence signature was linked with CD8⁺ effector memory and T_EMRA subsets and also senescent NK cells. A key observation was that the senescence markers in the skin lesions are age-independent of patients and were correlated with lesion size. Moreover, a striking expression of the senescence-associated secretory phenotype (SASP), pro-inflammatory cytokine and chemokines genes was found within lesions that were most strongly associated with the senescent CD8 T_EMRA subset. Collectively, our results confirm that there is a senescence transcriptomic signature in CL lesions and supports the hypothesis that lesional senescent cells have a major role in mediating immunopathology of the disease.

Gene Expression Differences Between Young Adults Based on Trauma History and Post-traumatic Stress Disorder

Background: The purpose of this study was to identify gene expression differences associated with post-traumatic stress disorder (PTSD) and trauma exposure (TE) in a three-group study design comprised of those with and without trauma exposure and PTSD. Methods: We conducted gene expression and gene network analyses in a sample (n = 45) composed of female subjects of European Ancestry (EA) with PTSD, TE without PTSD, and controls. Results: We identified 283 genes differentially expressed between PTSD-TE groups. In an independent sample of Veterans (n = 78) a small minority of these genes were also differentially expressed. We identified 7 gene network modules significantly associated with PTSD and TE (Bonferroni corrected p ≤ 0.05), which at a false discovery rate (FDR) of q ≤ 0.2, were significantly enriched for biological pathways involved in focal adhesion, neuroactive ligand receptor interaction, and immune related processes among others. Conclusions: This study uses gene network analyses to identify significant gene modules associated with PTSD, TE, and controls. On an individual gene level, we identified a large number of differentially expressed genes between PTSD-TE groups, a minority of which were also differentially expressed in the independent sample. We also demonstrate a lack of network module preservation between PTSD and TE, suggesting that the molecular signature of PTSD and trauma are likely independent of each other. Our results provide a basis for the identification of likely disease pathways and biomarkers involved in the etiology of PTSD.

High-fat diet induces a predisposition to follicular hyperkeratosis and neutrophilic folliculitis in mice

BACKGROUND

Neutrophilic folliculitis is an inflammatory condition of hair follicles. In some neutrophilic folliculitis, such as in patients with acne and hidradenitis suppurativa, follicular hyperkeratosis is also observed. Neutrophilic folliculitis is often induced and/or exacerbated by a high-fat diet (HFD). However, the molecular mechanisms by which an HFD affects neutrophilic folliculitis are not fully understood.

OBJECTIVE

Our aim was to elucidate how an HFD promotes the development of neutrophilic folliculitis.

METHODS

Mice were fed an HFD, and their skin was subjected to histologic, RNA sequencing, and imaging mass spectrometry analyses. To examine the effect of an HFD on neutrophil accumulation around the hair follicles, phorbol 12-myristate 13-acetate (PMA) was used as an irritant to the skin.

RESULTS

Histologic analysis revealed follicular hyperkeratosis in the skin of HFD-fed mice. RNA sequencing analysis showed that genes related to keratinization, especially in upper hair follicular keratinocytes, were significantly upregulated in HFD-fed mice. Application of PMA to the skin induced neutrophilic folliculitis in HFD-fed mice but not in mice fed a normal diet. Accumulation of neutrophils in the skin and around hair follicles was dependent on CXCR2 signaling, and CXCL1 (a CXCR2 ligand) was produced mainly by hair follicular keratinocytes. Imaging mass spectrometry analysis revealed an increase in fatty acids in the skin of HFD-fed mice. Application of these fatty acids to the skin induced follicular hyperkeratosis and caused PMA-induced neutrophilic folliculitis even in mice fed a normal diet.

CONCLUSION

An HFD can facilitate the development of neutrophilic folliculitis with the induction of hyperkeratosis of hair follicles and increased neutrophil infiltration around the hair follicles via CXCR2 signaling.

Progress and challenge for computational quantification of tissue immune cells

Tissue immune cells have long been recognized as important regulators for the maintenance of balance in the body system. Quantification of the abundance of different immune cells will provide enhanced understanding of the correlation between immune cells and normal or abnormal situations. Currently, computational methods to predict tissue immune cell compositions from bulk transcriptomes have been largely developed. Therefore, summarizing the advantages and disadvantages is appropriate. In addition, an examination of the challenges and possible solutions for these computational models will assist the development of this field. The common hypothesis of these models is that the expression of signature genes for immune cell types might represent the proportion of immune cells that contribute to the tissue transcriptome. In general, we grouped all reported tools into three groups, including reference-free, reference-based scoring and reference-based deconvolution methods. In this review, a summary of all the currently reported computational immune cell quantification tools and their applications, limitations, and perspectives are presented. Furthermore, some critical problems are found that have limited the performance and application of these models, including inadequate immune cell type, the collinearity problem, the impact of the tissue environment on the immune cell expression level, and the deficiency of standard datasets for model validation. To address these issues, tissue specific training datasets that include all known immune cells, a hierarchical computational framework, and benchmark datasets including both tissue expression profiles and the abundances of all the immune cells are proposed to further promote the development of this field.

Gene expression profiling reveals insights into infant immunological and febrile responses to group B meningococcal vaccine

Neisseria meningitidis is a major cause of meningitis and septicaemia. A MenB vaccine (4CMenB) was licensed by the European Medicines Agency in January 2013. Here we describe the blood transcriptome and proteome following infant immunisations with or without concomitant 4CMenB, to gain insight into the molecular mechanisms underlying post-vaccination reactogenicity and immunogenicity. Infants were randomised to receive control immunisations (PCV13 and DTaP-IPV-Hib) with or without 4CMenB at 2 and 4 months of age. Blood gene expression and plasma proteins were measured prior to, then 4 h, 24 h, 3 days or 7 days post-vaccination. 4CMenB vaccination was associated with increased expression of ENTPD7 and increased concentrations of 4 plasma proteins: CRP, G-CSF, IL-1RA and IL-6. Post-vaccination fever was associated with increased expression of SELL, involved in neutrophil recruitment. A murine model dissecting the vaccine components found the concomitant regimen to be associated with increased gene perturbation compared with 4CMenB vaccine alone with enhancement of pathways such as interleukin-3, -5 and GM-CSF signalling. Finally, we present transcriptomic profiles predictive of immunological and febrile responses following 4CMenB vaccine.

Metastasis-Specific Gene Expression in Autochthonous and Allograft Mouse Mammary Tumor Models: Stratification and Identification of Targetable Signatures

Breast cancer metastasis is a leading cause of cancer-related death of women in the United States. A hurdle in advancing metastasis-targeted intervention is the phenotypic heterogeneity between primary and secondary lesions. To identify metastasis-specific gene expression profiles we performed RNA-sequencing of breast cancer mouse models; analyzing metastases from models of various drivers and routes. We contrasted the models and identified common, targetable signatures. Allograft models exhibited more mesenchymal-like gene expression than genetically engineered mouse models (GEMM), and primary culturing of GEMM-derived metastatic tissue induced mesenchymal-like gene expression. In addition, metastasis-specific transcriptomes differed between tail vein and orthotopic injection of the same cell line. Gene expression common to models of spontaneous metastasis included sildenafil response and nicotine degradation pathways. Strikingly, in vivo sildenafil treatment significantly reduced metastasis by 54%, while nicotine significantly increased metastasis by 46%. These data suggest that (i) actionable metastasis-specific pathways can be readily identified, (ii) already available drugs may have great potential to alleviate metastatic incidence, and (iii) metastasis may be influenced greatly by lifestyle choices such as the choice to consume nicotine products. In summary, while mouse models of breast cancer metastasis vary in ways that must not be ignored, there are shared features that can be identified and potentially targeted therapeutically. IMPLICATIONS: The data we present here exposes critical variances between preclinical models of metastatic breast cancer and identifies targetable pathways integral to metastatic spread. VISUAL OVERVIEW: http://mcr.aacrjournals.org/content/molcanres/18/9/1278/F1.large.jpg.

Bone morphogenetic protein 8B promotes the progression of non-alcoholic steatohepatitis

Non-alcoholic steatohepatitis (NASH) is characterized by lipotoxicity, inflammation and fibrosis, ultimately leading to end-stage liver disease. The molecular mechanisms promoting NASH are poorly understood, and treatment options are limited. Here, we demonstrate that hepatic expression of bone morphogenetic protein 8B (BMP8B), a member of the transforming growth factor beta (TGFβ)-BMP superfamily, increases proportionally to disease stage in people and animal models with NASH. BMP8B signals via both SMAD2/3 and SMAD1/5/9 branches of the TGFβ-BMP pathway in hepatic stellate cells (HSCs), promoting their proinflammatory phenotype. In vivo, the absence of BMP8B prevents HSC activation, reduces inflammation and affects the wound-healing responses, thereby limiting NASH progression. Evidence is featured in primary human 3D microtissues modelling NASH, when challenged with recombinant BMP8. Our data show that BMP8B is a major contributor to NASH progression. Owing to the near absence of BMP8B in healthy livers, inhibition of BMP8B may represent a promising new therapeutic avenue for NASH treatment.

Early Transcriptional Changes within Liver, Adrenal Gland, and Lymphoid Tissues Significantly Contribute to Ebola Virus Pathogenesis in Cynomolgus Macaques

Ebola virus (EBOV) continues to pose a significant threat to human health, as evidenced by the 2013-2016 epidemic in West Africa and the ongoing outbreak in the Democratic Republic of the Congo. EBOV causes hemorrhagic fever, organ damage, and shock culminating in death, with case fatality rates as high as 90%. This high lethality combined with the paucity of licensed medical countermeasures makes EBOV a critical human pathogen. Although EBOV infection results in significant damage to the liver and the adrenal glands, little is known about the molecular signatures of injury in these organs. Moreover, while changes in peripheral blood cells are becoming increasingly understood, the host responses within organs and lymphoid tissues remain poorly characterized. To address this knowledge gap, we tracked longitudinal transcriptional changes in tissues collected from EBOV-Makona-infected cynomolgus macaques. Following infection, both liver and adrenal glands exhibited significant and early downregulation of genes involved in metabolism, coagulation, hormone synthesis, and angiogenesis; upregulated genes were associated with inflammation. Analysis of lymphoid tissues showed early upregulation of genes that play a role in innate immunity and inflammation and downregulation of genes associated with cell cycle and adaptive immunity. Moreover, transient activation of innate immune responses and downregulation of humoral immune responses in lymphoid tissues were confirmed with flow cytometry. Together, these data suggest that the liver, adrenal gland, and lymphatic organs are important sites of EBOV infection and that dysregulating the function of these vital organs contributes to the development of Ebola virus disease.IMPORTANCE Ebola virus (EBOV) remains a high-priority pathogen since it continues to cause outbreaks with high case fatality rates. Although it is well established that EBOV results in severe organ damage, our understanding of tissue injury in the liver, adrenal glands, and lymphoid tissues remains limited. We begin to address this knowledge gap by conducting longitudinal gene expression studies in these tissues, which were collected from EBOV-infected cynomolgus macaques. We report robust and early gene expression changes within these tissues, indicating they are primary sites of EBOV infection. Furthermore, genes involved in metabolism, coagulation, and adaptive immunity were downregulated, while inflammation-related genes were upregulated. These results indicate significant tissue damage consistent with the development of hemorrhagic fever and lymphopenia. Our study provides novel insight into EBOV-host interactions and elucidates how host responses within the liver, adrenal glands, and lymphoid tissues contribute to EBOV pathogenesis.

Immune correlates of postexposure vaccine protection against Marburg virus

Postexposure immunization can prevent disease and reduce transmission following pathogen exposure. The rapid immunostimulatory properties of recombinant vesicular stomatitis virus (rVSV)-based vaccines make them suitable postexposure treatments against the filoviruses Ebola virus and Marburg virus (MARV); however, the mechanisms that drive this protection are undefined. Previously, we reported 60–75% survival of rhesus macaques treated with rVSV vectors expressing MARV glycoprotein (GP) 20–30 minutes after a low dose exposure to the most pathogenic variant of MARV, Angola. Survival in this model was linked to production of GP-specific antibodies and lower viral load. To confirm these results and potentially identify novel correlates of postexposure protection, we performed a similar experiment, but analyzed plasma cytokine levels, frequencies of immune cell subsets, and the transcriptional response to infection in peripheral blood. In surviving macaques (80–89%), we observed induction of genes mapping to antiviral and interferon-related pathways early after treatment and a higher percentage of T helper 1 (Th1) and NK cells. In contrast, the response of non-surviving macaques was characterized by hypercytokinemia; a T helper 2 signature; recruitment of low HLA-DR expressing monocytes and regulatory T-cells; and transcription of immune checkpoint (e.g., PD-1, LAG3) genes. These results suggest dysregulated immunoregulation is associated with poor prognosis, whereas early innate signaling and Th1-skewed immunity are important for survival.

Comparative analysis of lung and blood transcriptomes in mice exposed to multi-walled carbon nanotubes

Pulmonary exposure to multi-walled carbon nanotubes (MWCNT) causes inflammation, fibroproliferation, immunotoxicity, and systemic responses in rodents. However, the search for representative biomarkers of exposure is an ongoing endeavor. Whole blood gene expression profiling is a promising new approach for the identification of novel disease biomarkers. We asked if the whole blood transcriptome reflects pathology-specific changes in lung gene expression caused by MWCNT. To answer this question, we performed mRNA sequencing analysis of the whole blood and lung in mice administered MWCNT or vehicle solution via pharyngeal aspiration and sacrificed 56 days later. The pattern of lung mRNA expression as determined using Ingenuity Pathway Analysis (IPA) was indicative of continued inflammation, immune cell trafficking, phagocytosis, and adaptive immune responses. Simultaneously, innate immunity-related transcripts (Plunc, Bpifb1, Reg3g) and cancer-related pathways were downregulated. IPA analysis of the differentially expressed genes in the whole blood suggested increased hematopoiesis, predicted activation of cancer/tumor development pathways, and atopy. There were several common upregulated genes between whole blood and lungs, important for adaptive immune responses: Cxcr1, Cd72, Sharpin, and Slc11a1. Trim24, important for T_H2 cell effector function, was downregulated in both datasets. Hla-dqa1 mRNA was upregulated in the lungs and downregulated in the blood, as was Lilrb4, which controls the reactivity of immune response. "Cancer" disease category had opposing activation status in the two datasets, while the only commonality was "Hypersensitivity". Transcriptome changes occurring in the lungs did not produce a completely replicable pattern in whole blood; however, specific systemic responses may be shared between transcriptomic profiles.

A mouse model of MSU-induced acute inflammation in vivo suggests imiquimod-dependent targeting of Il-1β as relevant therapy for gout patients

Rationale: The role of Monosodium Urate (MSU) crystals in gout pathophysiology is well described, as is the major impact of IL-1β in the inflammatory reaction that constitutes the hallmark of the disease. However, despite the discovery of the NLRP3 inflammasome and its role as a Pattern Recognition Receptor linking the detection of a danger signal (MSU) to IL-1β secretion in vitro, the precise mechanisms leading to joint inflammation in gout patients are still poorly understood. Methods: Acute urate crystal inflammation was obtained by subcutaneous injections of MSU crystals in mice. Symptoms were followed by scoring, cytokine quantification by ELISA and western blot, gene expression by RT-qPCR and RNAseq; Magnetic Resonance Imaging was also used to assess inflammation. Results: We provide an extensive clinical, biological and molecular characterization of an acute uratic inflammation mouse model which accurately mimics human gout. We report the efficacy of topical imiquimod treatment and its impact on Interferon-dependent down modulation of Il-1β gene expression in this experimental model. Conclusion: Our work reveals several key features of MSU-dependent inflammation and identifies novel therapeutic opportunities for gout patients.

Transcriptional Profiles and Stromal Changes Reveal Bone Marrow Adaptation to Early Breast Cancer in Association with Deregulated Circulating microRNAs

The presence of a growing tumor establishes a chronic state of inflammation that acts locally and systemically. Bone marrow responds to stress signals by expanding myeloid cells endowed with immunosuppressive functions, further fostering tumor growth and dissemination. How early in transformation the cross-talk with the bone marrow begins and becomes detectable in blood is unknown. Here, gene expression profiling of the bone marrow along disease progression in a spontaneous model of mammary carcinogenesis demonstrates that transcriptional modifications in the hematopoietic compartment occurred as early as preinvasive disease stages. The transcriptional profile showed downregulation of adaptive immunity and induction of programs related to innate immunity and response to danger signals triggered by activating transcription factor 3. Transcriptional reprogramming was paralleled by the expansion of myeloid populations at the expense of erythroid and B lymphoid fractions. Hematopoietic changes were associated with modifications of the bone marrow stromal architecture through relocalization and increased density in the interstitial area of Nestin⁺ mesenchymal cells expressing CXCL12 and myeloid cells expressing CXCL12 receptor CXCR4. These early events were concomitant with deregulation of circulating miRNAs, which were predicted regulators of transcripts downregulated in the bone marrow and involved in lymphoid differentiation and activation. These data provide a link between sensing of peripheral cancer initiation by the bone marrow and hematopoietic adaptation to distant noxia through transcriptional rewiring toward innate/inflammatory response programs. SIGNIFICANCE: The bone marrow senses distant tissue transformation at premalignant/preinvasive stages, suggesting that circulating messengers, intercepted in the blood, could serve as early diagnostic markers.

Variable gene expression and parasite load predict treatment outcome in cutaneous leishmaniasis

Patients infected with Leishmania braziliensis develop chronic lesions that often fail to respond to treatment with antiparasite drugs. To determine whether genes whose expression is highly variable in lesions between patients might influence disease outcome, we obtained biopsies of lesions from patients before treatment with pentavalent antimony and performed transcriptomic profiling on these clinical samples. We identified genes that were highly variably expressed between patients, and the variable expression of these genes correlated with treatment outcome. Among the most variable genes in all the patients were components of the cytolytic pathway, and the expression of these genes correlated with parasite load in the skin. We demonstrated that treatment failure was linked to the cytolytic pathway activated during infection. Using a host-pathogen marker profile of as few as three genes, we showed that eventual treatment outcome could be predicted before the start of treatment in two separate cohorts of patients with cutaneous leishmaniasis (n = 21 and n = 25). These findings raise the possibility of point-of-care diagnostic screening to identify patients at high risk of treatment failure and provide a rationale for a precision medicine approach to drug selection in cutaneous leishmaniasis. This work more broadly demonstrates the value of identifying genes of high variability in other diseases to better understand and predict diverse clinical outcomes.

Macrophage-associated wound healing contributes to African green monkey SIV pathogenesis control

Natural hosts of simian immunodeficiency virus (SIV) avoid AIDS despite lifelong infection. Here, we examined how this outcome is achieved by comparing a natural SIV host, African green monkey (AGM) to an AIDS susceptible species, rhesus macaque (RM). To asses gene expression profiles from acutely SIV infected AGMs and RMs, we developed a systems biology approach termed Conserved Gene Signature Analysis (CGSA), which compared RNA sequencing data from rectal AGM and RM tissues to various other species. We found that AGMs rapidly activate, and then maintain, evolutionarily conserved regenerative wound healing mechanisms in mucosal tissue. The wound healing protein fibronectin shows distinct tissue distribution and abundance kinetics in AGMs. Furthermore, AGM monocytes exhibit an embryonic development and repair/regeneration signature featuring TGF-β and concomitant reduced expression of inflammatory genes compared to RMs. This regenerative wound healing process likely preserves mucosal integrity and prevents inflammatory insults that underlie immune exhaustion in RMs.

Computational deconvolution of transcriptomics data from mixed cell populations

Summary

Gene expression analyses of bulk tissues often ignore cell type composition as an important confounding factor, resulting in a loss of signal from lowly abundant cell types. In this review, we highlight the importance and value of computational deconvolution methods to infer the abundance of different cell types and/or cell type-specific expression profiles in heterogeneous samples without performing physical cell sorting. We also explain the various deconvolution scenarios, the mathematical approaches used to solve them and the effect of data processing and different confounding factors on the accuracy of the deconvolution results.

Contact

katleen.depreter@ugent.be.

Supplementary information

Supplementary data are available at Bioinformatics online.

Predicting bacterial infection outcomes using single cell RNA-sequencing analysis of human immune cells

Complex interactions between different host immune cell types can determine the outcome of pathogen infections. Advances in single cell RNA-sequencing (scRNA-seq) allow probing of these immune interactions, such as cell-type compositions, which are then interpreted by deconvolution algorithms using bulk RNA-seq measurements. However, not all aspects of immune surveillance are represented by current algorithms. Here, using scRNA-seq of human peripheral blood cells infected with Salmonella, we develop a deconvolution algorithm for inferring cell-type specific infection responses from bulk measurements. We apply our dynamic deconvolution algorithm to a cohort of healthy individuals challenged ex vivo with Salmonella, and to three cohorts of tuberculosis patients during different stages of disease. We reveal cell-type specific immune responses associated not only with ex vivo infection phenotype but also with clinical disease stage. We propose that our approach provides a predictive power to identify risk for disease, and human infection outcomes.

Complex interactions between different host immune cell types can determine the outcome of pathogen infections. Here, Avraham and colleagues present a deconvolution algorithm that uses single-cell RNA and bulk RNA sequencing measurements of pathogen-infected cells to predict disease risk outcomes.

Antiviral Innate Responses Induced by VSV-EBOV Vaccination Contribute to Rapid Protection

Ebola virus (EBOV) is a single-stranded RNA virus that causes Ebola virus disease (EVD), characterized by excessive inflammation, lymphocyte apoptosis, hemorrhage, and coagulation defects leading to multiorgan failure and shock. Recombinant vesicular stomatitis virus expressing the EBOV glycoprotein (VSV-EBOV), which is highly efficacious against lethal challenge in nonhuman primates, is the only vaccine that successfully completed a phase III clinical trial. Additional studies showed VSV-EBOV provides complete and partial protection to macaques immunized 7 and 3 days before EBOV challenge, respectively. However, the mechanisms by which this live-attenuated vaccine elicits rapid protection are only partially understood. To address this, we carried out a longitudinal transcriptome analysis of host responses in whole-blood samples collected from cynomolgus macaques vaccinated with VSV-EBOV 28, 21, 14, 7, and 3 days before EBOV challenge. Our findings indicate the transcriptional response to the vaccine peaks 7 days following vaccination and contains signatures of both innate antiviral immunity as well as B-cell activation. EBOV challenge 1 week after vaccination resulted in large gene expression changes suggestive of a recall adaptive immune response 14 days postchallenge. Lastly, the timing and magnitude of innate immunity and interferon-stimulated gene expression correlated with viral burden and disease outcome in animals vaccinated 3 days before challenge.IMPORTANCE Ebola virus (EBOV) is the causative agent of Ebola virus disease (EVD), a deadly disease and major public health threat worldwide. A safe and highly efficacious vesicular stomatitis virus-based vaccine against EBOV is the only platform that has successfully completed phase III clinical trials and has been used in recent and ongoing outbreaks. Earlier studies showed that antibodies are the main mode of protection when this vaccine is administered 28 days before EBOV challenge. Recently, we showed this vaccine can provide protection when administered as early as 3 days before challenge and before antibodies are detected. This study seeks to identify the mechanisms of rapid protection, which in turn will pave the way for improved vaccines and therapeutics. Additionally, this study provides insight into host gene expression signatures that could provide early biomarkers to identify infected individuals who are at highest risk of poor outcomes.

Aicardi-Goutières syndrome gene Rnaseh2c is a metastasis susceptibility gene in breast cancer

Breast cancer is the second leading cause of cancer-related deaths in the United States, with the majority of these deaths due to metastatic lesions rather than the primary tumor. Thus, a better understanding of the etiology of metastatic disease is crucial for improving survival. Using a haplotype mapping strategy in mouse and shRNA-mediated gene knockdown, we identified Rnaseh2c, a scaffolding protein of the heterotrimeric RNase H2 endoribonuclease complex, as a novel metastasis susceptibility factor. We found that the role of Rnaseh2c in metastatic disease is independent of RNase H2 enzymatic activity, and immunophenotyping and RNA-sequencing analysis revealed engagement of the T cell-mediated adaptive immune response. Furthermore, the cGAS-Sting pathway was not activated in the metastatic cancer cells used in this study, suggesting that the mechanism of immune response in breast cancer is different from the mechanism proposed for Aicardi-Goutières Syndrome, a rare interferonopathy caused by RNase H2 mutation. These results suggest an important novel, non-enzymatic role for RNASEH2C during breast cancer progression and add Rnaseh2c to a panel of genes we have identified that together could determine patients with high risk for metastasis. These results also highlight a potential new target for combination with immunotherapies and may contribute to a better understanding of the etiology of Aicardi-Goutières Syndrome autoimmunity.

Chronic ethanol consumption alters lamina propria leukocyte response to stimulation in a region-dependent manner

Chronic heavy alcohol consumption, also referred to as chronic heavy drinking (CHD), results in intestinal injury characterized by increased permeability, dysbiosis, nutrient malabsorption, potentially higher susceptibility to infection, and increased risk of colorectal cancer. However, our understanding of the mechanisms by which CHD results in intestinal damage remains incomplete. Here, we investigated the impact of chronic drinking on transcriptional and functional responses of lamina propria leukocytes (LPLs) isolated from the 4 major gut sections. Although no significant differences were detected between LPLs isolated from the ethanol and control groups at resting state within each major gut section, our analysis uncovered key regional differences in composition and function of LPLs independent of alcohol consumption. However, in response to phorbol myristate acetate and ionomycin, duodenal LPLs from ethanol-drinking animals generated a dampened response, whereas jejunal and ileal LPLs from ethanol-drinking animals produced a heightened response. Transcriptional responses following stimulation were pronounced in ileal and duodenal LPLs from the ethanol-drinking group but less evident in jejunal and colonic LPLs compared with controls, suggesting a more significant impact of alcohol on these gut regions. The altered intestinal LPL function detected in our study reveals remarkable region specificity and novel insight into potential mechanisms of intestinal injury associated with CHD.-Barr, T., Lewis, S. A., Sureshchandra, S., Doratt, B., Grant, K. A., Messaoudi, I. Chronic ethanol consumption alters lamina propria leukocyte response to stimulation in a region-dependent manner.

Cell composition analysis of bulk genomics using single-cell data

Single-cell expression profiling (scRNA-seq) is a rich resource of cellular heterogeneity. While profiling every sample under study would be advantageous, it is time-consuming and costly. Here we introduce Cell Population Mapping (CPM), a deconvolution algorithm in which the composition of cell types and states is inferred from the bulk transcriptome using reference scRNA-seq profiles ('scBio' CRAN R-package). Analysis of individual variations in lungs of influenza virus-infected mice, using CPM, revealed 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 was confirmed in subsequent experiments and was 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.

Altered Immunity and Microbial Dysbiosis in Aged Individuals With Long-Term Controlled HIV Infection

The introduction of highly active antiretroviral therapy (HAART) resulted in a significant increase in life expectancy for HIV patients. Indeed, in 2015, 45% of the HIV+ individuals in the United States were ≥55 years of age. Despite improvements in diagnosis and treatment of HIV infection, geriatric HIV+ patients suffer from higher incidence of comorbidities compared to age-matched HIV- individuals. Both chronic inflammation and dysbiosis of the gut microbiome are believed to be major contributors to this phenomenon, however carefully controlled studies investigating the impact of long-term (>10 years) controlled HIV (LTC-HIV) infection are lacking. To address this question, we profiled circulating immune cells, immune mediators, and the gut microbiome from elderly (≥55 years old) LTC-HIV+ and HIV- gay men living in the Palm Springs area. LTC-HIV+ individuals had lower frequency of circulating monocytes and CD4+ T-cells, and increased frequency CD8+ T-cells. Moreover, levels of systemic INFγ and several growth factors were increased while levels of IL-2 and several chemokines were reduced. Upon stimulation, immune cells from LTC-HIV+ individuals produced higher levels of pro-inflammatory cytokines. Last but not least, the gut microbiome of LTC-HIV+ individuals was enriched in bacterial taxa typically found in the oral cavity suggestive of loss of compartmentalization, while levels of beneficial butyrate producing taxa were reduced. Additionally, prevalence of Prevotella negatively correlated with CD4+ T-cells numbers in LTC-HIV+ individuals. These results indicate that despite long-term adherence and undetectable viral loads, LTC-HIV infection results in significant shifts in immune cell frequencies and gut microbial communities.

Comparative Transcriptomics in Ebola Makona-Infected Ferrets, Nonhuman Primates, and Humans

The domestic ferret is a uniformly lethal model of infection for 3 species of Ebolavirus known to be pathogenic in humans. Reagents to systematically analyze the ferret host response to infection are lacking; however, the recent publication of a draft ferret genome has opened the potential for transcriptional analysis of ferret models of disease. In this work, we present comparative analysis of longitudinally sampled blood taken from ferrets and nonhuman primates infected with lethal doses of the Makona variant of Zaire ebolavirus. Strong induction of proinflammatory and prothrombotic signaling programs were present in both ferrets and nonhuman primates, and both transcriptomes were similar to previously published datasets of fatal cases of human Ebola virus infection.

Quantifying tumor-infiltrating immune cells from transcriptomics data

By exerting pro- and anti-tumorigenic actions, tumor-infiltrating immune cells can profoundly influence tumor progression, as well as the success of anti-cancer therapies. Therefore, the quantification of tumor-infiltrating immune cells holds the promise to unveil the multi-faceted role of the immune system in human cancers and its involvement in tumor escape mechanisms and response to therapy. Tumor-infiltrating immune cells can be quantified from RNA sequencing data of human tumors using bioinformatics approaches. In this review, we describe state-of-the-art computational methods for the quantification of immune cells from transcriptomics data and discuss the open challenges that must be addressed to accurately quantify immune infiltrates from RNA sequencing data of human bulk tumors.

The immunosuppressive face of sepsis early on intensive care unit-A large-scale microarray meta-analysis

Background

Sepsis is defined as a life-threatening condition, resulting from a dysregulated and harmful response of the hosts’ immune system to infection. Apart from this, the (over-)compensating mechanisms counterbalancing the inflammatory response have been proven to render the host susceptible to further infections and increase delayed mortality. Our study aimed to unravel the heterogeneity of immune response in early sepsis and to explain the biology behind it.

Methods

A systematic search of public repositories yielded 949 microarray samples from patients with sepsis of different infectious origin and early after clinical manifestation. These were merged into a meta-expression set, and after applying sequential conservative bioinformatics filtering, an in-deep analysis of transcriptional heterogeneity, as well as a comparison to samples of healthy controls was performed.

Results

We can identify two distinct clusters of patients (cluster 1: 655 subjects, cluster 2: 294 subjects) according to their global blood transcriptome. While both clusters exhibit only moderate differences in direct comparison, a comparison of both clusters individually to healthy controls yielded strong expression changes of genes involved in immune responses. Both comparisons found similar regulated genes, with a stronger dysregulation occurring in the larger patient cluster and implicating a loss of monocyte and T cell function, co-occurring with an activation of neutrophil granulocytes.

Conclusion

We propose a consistent—but in its extent varying—presence of immunosuppression, occurring as early in sepsis as its clinical manifestation and irrespective of the infectious origin. While certain cell types possess contradictory activation states, our finding underlines the urgent need for an early host-directed therapy of sepsis side-by-side with antibiotics.

Cell-type deconvolution from DNA methylation: a review of recent applications

Recent advances in cell-type deconvolution approaches are adding to our understanding of the biology underlying disease development and progression. DNA methylation (DNAm) can be used as a biomarker of cell types, and through deconvolution approaches, to infer underlying cell type proportions. Cell-type deconvolution algorithms have two main categories: reference-based and reference-free. Reference-based algorithms are supervised methods that determine the underlying composition of cell types within a sample by leveraging differentially methylated regions (DMRs) specific to cell type, identified from DNAm measures of purified cell populations. Reference-free algorithms are unsupervised methods for use when cell-type specific DMRs are not available, allowing scientists to estimate putative cellular proportions or control for potential confounding from cell type. Reference-based deconvolution is typically applied to blood samples and has potentiated our understanding of the relation between immune profiles and disease by allowing estimation of immune cell proportions from archival DNA. Bioinformatic analyses using DNAm to infer immune cell proportions, part of a new field known as Immunomethylomics, provides a new direction for consideration in epigenome wide association studies (EWAS).

Computational de novo discovery of distinguishing genes for biological processes and cell types in complex tissues

Bulk tissue samples examined by gene expression studies are usually heterogeneous. The data gained from these samples display the confounding patterns of mixtures consisting of multiple cell types or similar cell types in various functional states, which hinders the elucidation of the molecular mechanisms underlying complex biological phenomena. A realistic approach to compensate for the limitations of experimentally separating homogenous cell populations from mixed tissues is to computationally identify cell-type specific patterns from bulk, heterogeneous measurements. We designed the CellDistinguisher algorithm to analyze the gene expression data of mixed samples, identifying genes that best distinguish biological processes and cell types. Coupled with a deconvolution algorithm that takes cell type specific gene lists as input, we show that CellDistinguisher performs as well as partial deconvolution algorithms in predicting cell type composition without the need for prior knowledge of cell type signatures. This approach is also better in predicting cell type signatures than the one-step traditional complete deconvolution methods. To illustrate its wide applicability, the algorithm was tested on multiple publicly available data sets. In each case, CellDistinguisher identified genes reflecting biological processes typical for the tissues and development stages of interest and estimated the sample compositions accurately.

Sequence variation in PPP1R13L results in a novel form of cardio-cutaneous syndrome

Dilated cardiomyopathy (DCM) is a life-threatening disorder whose genetic basis is heterogeneous and mostly unknown. Five Arab Christian infants, aged 4-30 months from four families, were diagnosed with DCM associated with mild skin, teeth, and hair abnormalities. All passed away before age 3. A homozygous sequence variation creating a premature stop codon at PPP1R13L encoding the iASPP protein was identified in three infants and in the mother of the other two. Patients' fibroblasts and PPP1R13L-knocked down human fibroblasts presented higher expression levels of pro-inflammatory cytokine genes in response to lipopolysaccharide, as well as Ppp1r13l-knocked down murine cardiomyocytes and hearts of Ppp1r13l-deficient mice. The hypersensitivity to lipopolysaccharide was NF-κB-dependent, and its inducible binding activity to promoters of pro-inflammatory cytokine genes was elevated in patients' fibroblasts. RNA sequencing of Ppp1r13l-knocked down murine cardiomyocytes and of hearts derived from different stages of DCM development in Ppp1r13l-deficient mice revealed the crucial role of iASPP in dampening cardiac inflammatory response. Our results determined PPP1R13L as the gene underlying a novel autosomal-recessive cardio-cutaneous syndrome in humans and strongly suggest that the fatal DCM during infancy is a consequence of failure to regulate transcriptional pathways necessary for tuning cardiac threshold response to common inflammatory stressors.