Validation of Immune Cell Modules in Multicellular Transcriptomic Data

Dysfunctional bronchoalveolar effector memory CD8+ T cells in tuberculosis-exposed people living with antiretroviral-naïve HIV infection

HIV causes susceptibility to respiratory pathogens, including tuberculosis (TB), but the underlying immunological mechanisms remain incompletely understood. We obtained whole blood and bronchoalveolar lavage (BAL) from TB-exposed people in the presence or absence of antiretroviral-naïve HIV co-infection. Bulk transcriptional profiling demonstrated compartment-specific enrichment of immunological processes. Systems-level deconvolution of whole blood from people living with HIV identified elevated type I and type II interferon cytokine activity and T cell proliferation. Transcriptional modules derived from both peripheral blood and sorted BAL immune cells demonstrated an increased frequency of effector memory CD8 T cells in whole BAL samples. Both compartments displayed reduced induction of CD8 T-cell-derived interleukin-17A (IL-17A) in people with HIV, associated with elevated T cell regulatory molecule expression. The data suggest that dysfunctional CD8 T cell responses in uncontrolled HIV may contribute to compromised respiratory immunity to pathogens, a process that could be modulated by host-directed therapies that target CD8 T cell effector functions.

Immune cell type signature discovery and random forest classification for analysis of single cell gene expression datasets

Background

Robust immune cell gene expression signatures are central to the analysis of single cell studies. Nearly all known sets of immune cell signatures have been derived by making use of only single gene expression datasets. Utilizing the power of multiple integrated datasets could lead to high-quality immune cell signatures which could be used as superior inputs to machine learning-based cell type classification approaches.

Results

We established a novel workflow for the discovery of immune cell type signatures based primarily on gene-versus-gene expression similarity. It leverages multiple datasets, here seven single cell expression datasets from six different cancer types and resulted in eleven immune cell type-specific gene expression signatures. We used these to train random forest classifiers for immune cell type assignment for single-cell RNA-seq datasets. We obtained similar or better prediction results compared to commonly used methods for cell type assignment in independent benchmarking datasets. Our gene signature set yields higher prediction scores than other published immune cell type gene sets in random forest-based cell type classification. We further demonstrate how our approach helps to avoid bias in downstream statistical analyses by re-analysis of a published IFN stimulation experiment.

Discussion and conclusion

We demonstrated the quality of our immune cell signatures and their strong performance in a random forest-based cell typing approach. We argue that classifying cells based on our comparably slim sets of genes accompanied by a random forest-based approach not only matches or outperforms widely used published approaches. It also facilitates unbiased downstream statistical analyses of differential gene expression between cell types for significantly more genes compared to previous cell classification algorithms.

Dysfunctional effector memory CD8 T cells in the bronchoalveolar compartment of people living with HIV

Mechanisms by which HIV causes susceptibility to respiratory pathogens remain incompletely understood. We obtained whole blood and bronchoalveolar lavage (BAL) from people with latent TB infection in the presence or absence of antiretroviral-naïve HIV co-infection. Transcriptomic and flow cytometric analyses demonstrated HIV-associated cell proliferation plus type I interferon activity in blood and effector memory CD8 T-cells in BAL. Both compartments displayed reduced induction of CD8 T-cell-derived IL-17A in people with HIV, associated with elevated T-cell regulatory molecule expression. The data suggest that dysfunctional CD8 T-cell responses in uncontrolled HIV contribute to susceptibility to secondary bacterial infections, including tuberculosis.

Parvimonas micra is associated with tumour immune profiles in molecular subtypes of colorectal cancer

The importance of the tumour microbiome in different aspects of colorectal cancer (CRC) has been increasingly recognised, but many questions remain. The aim of this study was to explore the effect of specific CRC associated microbes on the tumour immune response, which has a considerable prognostic value in CRC. We applied specific qPCR to detect Parvimonas micra and Fusobacterium nucleatum in tumour tissues from an immunologically well-characterised cohort of 69 CRC patients. This cohort included detailed analyses of immune profiles based on flow cytometry and transcriptomics in tumour tissue and blood, along with comprehensive analyses of molecular subtypes. P. micra and F. nucleatum were detected in 24% and 64% of tumour tissues, respectively. We found a significant association of P. micra with high-grade tumours and tumours of CMS1 subtype. F. nucleatum was significantly associated with right-sided tumours, microsatellite instability, and CMS1 tumours. The immunological analyses revealed significant associations of P. micra with activated CD69⁺ T lymphocytes and increased antigen-presenting HLA-DR⁺ B lymphocytes. P. micra was also positively associated with M1 and M2 macrophage traits. The impact of P. micra tumour colonisation on the immune response was further assessed using transcriptomics in validation of our findings. No associations were found between F. nucleatum and immune profiles in this study. Our findings support novel associations between P. micra and the immune response in CRC. A better understanding of these interactions might help to identify important predictive and prognostic tools as well as new targets for therapy.

Rapid synchronous type 1 IFN and virus-specific T cell responses characterize first wave non-severe SARS-CoV-2 infections

Effective control of SARS-CoV-2 infection on primary exposure may reveal correlates of protective immunity to future variants, but we lack insights into immune responses before or at the time virus is first detected. We use blood transcriptomics, multiparameter flow cytometry, and T cell receptor (TCR) sequencing spanning the time of incident non-severe infection in unvaccinated virus-naive individuals to identify rapid type 1 interferon (IFN) responses common to other acute respiratory viruses and cell proliferation responses that discriminate SARS-CoV-2 from other viruses. These peak by the time the virus is first detected and sometimes precede virus detection. Cell proliferation is most evident in CD8 T cells and associated with specific expansion of SARS-CoV-2-reactive TCRs, in contrast to virus-specific antibodies, which lag by 1-2 weeks. Our data support a protective role for early type 1 IFN and CD8 T cell responses, with implications for development of universal T cell vaccines.

Exaggerated IL-17A activity in human in vivo recall responses discriminates active tuberculosis from latent infection and cured disease

Host immune responses at the site of Mycobacterium tuberculosis (Mtb) infection can mediate pathogenesis of tuberculosis (TB) and onward transmission of infection. We hypothesized that pathological immune responses would be enriched at the site of host-pathogen interactions modelled by a standardized tuberculin skin test (TST) challenge in patients with active TB compared to those without disease, and interrogated immune responses by genome-wide transcriptional profiling. We show exaggerated interleukin (IL)-17A and Th17 responses among 48 individuals with active TB compared to 191 with latent TB infection, associated with increased neutrophil recruitment and matrix metalloproteinase-1 expression, both involved in TB pathogenesis. Curative antimicrobial treatment reversed these observed changes. Increased IL-1β and IL-6 responses to mycobacterial stimulation were evident in both circulating monocytes and in molecular changes at the site of TST in individuals with active TB, supporting a model in which monocyte-derived IL-1β and IL-6 promote Th17 differentiation within tissues. Modulation of these cytokine pathways may provide a rational strategy for host-directed therapy in active TB.

Persistent T Cell Repertoire Perturbation and T Cell Activation in HIV After Long Term Treatment

Objective

In people living with HIV (PLHIV), we sought to test the hypothesis that long term anti-retroviral therapy restores the normal T cell repertoire, and investigate the functional relationship of residual repertoire abnormalities to persistent immune system dysregulation.

Methods

We conducted a case-control study in PLHIV and HIV-negative volunteers, of circulating T cell receptor repertoires and whole blood transcriptomes by RNA sequencing, complemented by metadata from routinely collected health care records.

Results

T cell receptor sequencing revealed persistent abnormalities in the clonal T cell repertoire of PLHIV, characterized by reduced repertoire diversity and oligoclonal T cell expansion correlated with elevated CD8 T cell counts. We found no evidence that these expansions were driven by cytomegalovirus or another common antigen. Increased frequency of long CDR3 sequences and reduced frequency of public sequences among the expanded clones implicated abnormal thymic selection as a contributing factor. These abnormalities in the repertoire correlated with systems level evidence of persistent T cell activation in genome-wide blood transcriptomes.

Conclusions

The diversity of T cell receptor repertoires in PLHIV on long term anti-retroviral therapy remains significantly depleted, and skewed by idiosyncratic clones, partly attributable to altered thymic output and associated with T cell mediated chronic immune activation. Further investigation of thymic function and the antigenic drivers of T cell clonal selection in PLHIV are critical to efforts to fully re-establish normal immune function.

Transcriptional response modules characterize IL-1β and IL-6 activity in COVID-19

Dysregulated IL-1β and IL-6 responses have been implicated in the pathogenesis of severe Coronavirus Disease 2019 (COVID-19). Innovative approaches for evaluating the biological activity of these cytokines in vivo are urgently needed to complement clinical trials of therapeutic targeting of IL-1β and IL-6 in COVID-19. We show that the expression of IL-1β or IL-6 inducible transcriptional signatures (modules) reflects the bioactivity of these cytokines in immunopathology modelled by juvenile idiopathic arthritis (JIA) and rheumatoid arthritis. In COVID-19, elevated expression of IL-1β and IL-6 response modules, but not the cytokine transcripts themselves, is a feature of infection in the nasopharynx and blood but is not associated with severity of COVID-19 disease, length of stay, or mortality. We propose that IL-1β and IL-6 transcriptional response modules provide a dynamic readout of functional cytokine activity in vivo, aiding quantification of the biological effects of immunomodulatory therapies in COVID-19.

Transcriptional profiling of pediatric cholestatic livers identifies three distinct macrophage populations

Background & aims

Limited understanding of the role for specific macrophage subsets in the pathogenesis of cholestatic liver injury is a barrier to advancing medical therapy. Macrophages have previously been implicated in both the mal-adaptive and protective responses in obstructive cholestasis. Recently two macrophage subsets were identified in non-diseased human liver; however, no studies to date fully define the heterogeneous macrophage subsets during the pathogenesis of cholestasis. Here, we aim to further characterize the transcriptional profile of macrophages in pediatric cholestatic liver disease.

Methods

We isolated live hepatic immune cells from patients with biliary atresia (BA), Alagille syndrome (ALGS), and non-cholestatic pediatric liver by fluorescence activated cell sorting. Through single-cell RNA sequencing analysis and immunofluorescence, we characterized cholestatic macrophages. We next compared the transcriptional profile of pediatric cholestatic and non-cholestatic macrophage populations to previously published data on normal adult hepatic macrophages.

Results

We identified 3 distinct macrophage populations across cholestatic liver samples and annotated them as lipid-associated macrophages, monocyte-like macrophages, and adaptive macrophages based on their transcriptional profile. Immunofluorescence of liver tissue using markers for each subset confirmed their presence across BA (n = 6) and ALGS (n = 6) patients. Cholestatic macrophages demonstrated reduced expression of immune regulatory genes as compared to normal hepatic macrophages and were distinct from macrophage populations defined in either healthy adult or pediatric non-cholestatic liver.

Conclusions

We are the first to perform single-cell RNA sequencing on human pediatric cholestatic liver and identified three macrophage subsets with distinct transcriptional signatures from healthy liver macrophages. Further analyses will identify similarities and differences in these macrophage sub-populations across etiologies of cholestatic liver disease. Taken together, these findings may allow for future development of targeted therapeutic strategies to reprogram macrophages to an immune regulatory phenotype and reduce cholestatic liver injury.

Vitamin D3 replacement enhances antigen-specific immunity in older adults

Introduction

Ageing is associated with increased number of infections, decreased vaccine efficacy and increased systemic inflammation termed inflammageing. These changes are reflected by reduced recall responses to varicella zoster virus (VZV) challenge in the skin of older adults. Vitamin D deficiency is more common in the old and has been associated with frailty and increased inflammation. In addition, vitamin D increases immunoregulatory mechanisms and therefore has the potential to inhibit inflammageing.

Objectives

We investigated the use of vitamin D3 replacement to enhance cutaneous antigen-specific immunity in older adults (≥65 years).

Methods

Vitamin D insufficient older adults (n = 18) were administered 6400IU of vitamin D3/day orally for 14 weeks. Antigen-specific immunity to VZV was assessed by clinical score assessment of the injection site and transcriptional analysis of skin biopsies collected from challenged injection sites pre- and post-vitamin D3 replacement.

Results

We showed that older adults had reduced VZV-specific cutaneous immune response and increased non-specific inflammation as compared to young. Increased non-specific inflammation observed in the skin of older adults negatively correlated with vitamin D sufficiency. We showed that vitamin D3 supplementation significantly increased the response to cutaneous VZV antigen challenge in older adults. This enhancement was associated with a reduction in inflammatory monocyte infiltration with a concomitant enhancement of T cell recruitment to the site of antigen challenge in the skin.

Conclusion

Vitamin D3 replacement can boost antigen-specific immunity in older adults with sub-optimal vitamin D status.

Transcriptional response modules characterise IL-1β and IL-6 activity in COVID-19

Dysregulated IL-1β and IL-6 responses have been implicated in the pathogenesis of severe Coronavirus Disease 2019 (COVID-19). Innovative approaches for evaluating the biological activity of these cytokines in vivo are urgently needed to complement clinical trials of therapeutic targeting of IL-1β and IL-6 in COVID-19. We show that the expression of IL-1β or IL-6 inducible transcriptional signatures (modules) reflects the bioactivity of these cytokines in immunopathology modelled by juvenile idiopathic arthritis (JIA) and rheumatoid arthritis. In COVID-19, elevated expression of IL-1β and IL-6 response modules, but not the cytokine transcripts themselves, is a feature of infection in the nasopharynx and blood, but is not associated with severity of COVID-19 disease, length of stay or mortality. We propose that IL-1β and IL-6 transcriptional response modules provide a dynamic readout of functional cytokine activity in vivo, aiding quantification of the biological effects of immunomodulatory therapies in COVID-19.

A Detailed Flow Cytometric Analysis of Immune Activity Profiles in Molecular Subtypes of Colorectal Cancer

Simple Summary

Colorectal cancer is one of the deadliest cancers worldwide, with around 40% of patients dying from distant metastasis. Tumour immune cell infiltration has powerful positive prognostic value in this disease, suggesting immunotherapy as a potential treatment modality. The aim of this explorative study was to assess in detail the local and systemic immune response in different molecular subgroups of colorectal cancer. An improved molecular understanding of the disease may lead to important advances in personalised medicine, identifying prognostic and predictive tools, in addition to new therapeutic targets.

Abstract

The local anti-tumour immune response has important prognostic value in colorectal cancer (CRC). In the era of immunotherapy, a better understanding of the immune response in molecular subgroups of CRC may lead to significant advances in personalised medicine. On this note, microsatellite instable (MSI) tumours have been characterised by increased immune infiltration, suggesting MSI as a marker for immune inhibitor checkpoint therapy. Here, we used flow cytometry to perform a comprehensive analysis of immune activity profiles in tumour tissues, adjacent non-malignant tissues and blood, from a cohort of 69 CRC patients. We found several signs of immune suppression in tumours compared to adjacent non-malignant tissues, including T cells more often expressing the immune checkpoint molecules programmed cell death protein (PD-1) and cytotoxic T lymphocyte-associated protein 4 (CTLA-4). We further analysed immune cell infiltration in molecular subgroups of CRC. MSI tumours were indeed found to be associated with increased immune infiltration, including increased fractions of PD-1⁺ T cells. No correlation was, however, found between MSI and the fraction of CTLA-4⁺ T cells. Interestingly, within the group of patients with microsatellite stable (MSS) tumours, some also presented with increased immune infiltration, including comparably high portions of PD-1⁺ T cells, but also CTLA-4⁺ T cells. Furthermore, no correlation was found between PD-1⁺ and CTLA-4⁺ T cells, suggesting that different tumours may, to some extent, be regulated by different immune checkpoints. We further evaluated the distribution of immune activity profiles in the consensus molecular subtypes of CRC. In conclusion, our findings suggest that different immune checkpoint inhibitors may be beneficial for selected CRC patients irrespective of MSI status. Improved predictive tools are required to identify these patients.

Assessing the Impact of Sample Heterogeneity on Transcriptome Analysis of Human Diseases Using MDP Webtool

Transcriptome analyses have increased our understanding of the molecular mechanisms underlying human diseases. Most approaches aim to identify significant genes by comparing their expression values between healthy subjects and a group of patients with a certain disease. Given that studies normally contain few samples, the heterogeneity among individuals caused by environmental factors or undetected illnesses can impact gene expression analyses. We present a systematic analysis of sample heterogeneity in a variety of gene expression studies relating to inflammatory and infectious diseases and show that novel immunological insights may arise once heterogeneity is addressed. The perturbation score of samples is quantified using nonperturbed subjects (i.e., healthy subjects) as a reference group. Such a score allows us to detect outlying samples and subgroups of diseased patients and even assess the molecular perturbation of single cells infected with viruses. We also show how removal of outlying samples can improve the "signal" of the disease and impact detection of differentially expressed genes. The method is made available via the mdp Bioconductor R package and as a user-friendly webtool, webMDP, available at http://mdp.sysbio.tools.

Spatial heterogeneity of the T cell receptor repertoire reflects the mutational landscape in lung cancer

Somatic mutations together with immunoediting drive extensive heterogeneity within non-small-cell lung cancer (NSCLC). Herein we examine heterogeneity of the T cell antigen receptor (TCR) repertoire. The number of TCR sequences selectively expanded in tumors varies within and between tumors and correlates with the number of nonsynonymous mutations. Expanded TCRs can be subdivided into TCRs found in all tumor regions (ubiquitous) and those present in a subset of regions (regional). The number of ubiquitous and regional TCRs correlates with the number of ubiquitous and regional nonsynonymous mutations, respectively. Expanded TCRs form part of clusters of TCRs of similar sequence, suggestive of a spatially constrained antigen-driven process. CD8+ tumor-infiltrating lymphocytes harboring ubiquitous TCRs display a dysfunctional tissue-resident phenotype. Ubiquitous TCRs are preferentially detected in the blood at the time of tumor resection as compared to routine follow-up. These findings highlight a noninvasive method to identify and track relevant tumor-reactive TCRs for use in adoptive T cell immunotherapy.

An Efficient and Flexible Method for Deconvoluting Bulk RNA-Seq Data with Single-Cell RNA-Seq Data

Estimating cell type compositions for complex diseases is an important step to investigate the cellular heterogeneity for understanding disease etiology and potentially facilitate early disease diagnosis and prevention. Here, we developed a computationally statistical method, referring to Multi-Omics Matrix Factorization (MOMF), to estimate the cell-type compositions of bulk RNA sequencing (RNA-seq) data by leveraging cell type-specific gene expression levels from single-cell RNA sequencing (scRNA-seq) data. MOMF not only directly models the count nature of gene expression data, but also effectively accounts for the uncertainty of cell type-specific mean gene expression levels. We demonstrate the benefits of MOMF through three real data applications, i.e., Glioblastomas (GBM), colorectal cancer (CRC) and type II diabetes (T2D) studies. MOMF is able to accurately estimate disease-related cell type proportions, i.e., oligodendrocyte progenitor cells and macrophage cells, which are strongly associated with the survival of GBM and CRC, respectively.

Spatial Network Mapping of Pulmonary Multidrug-Resistant Tuberculosis Cavities Using RNA Sequencing

Rationale: There is poor understanding about protective immunity and the pathogenesis of cavitation in patients with tuberculosis.Objectives: To map pathophysiological pathways at anatomically distinct positions within the human tuberculosis cavity.Methods: Biopsies were obtained from eight predetermined locations within lung cavities of patients with multidrug-resistant tuberculosis undergoing therapeutic surgical resection (n = 14) and healthy lung tissue from control subjects without tuberculosis (n = 10). RNA sequencing, immunohistochemistry, and bacterial load determination were performed at each cavity position. Differentially expressed genes were normalized to control subjects without tuberculosis, and ontologically mapped to identify a spatially compartmentalized pathophysiological map of the cavity. In silico perturbation using a novel distance-dependent dynamical sink model was used to investigate interactions between immune networks and bacterial burden, and to integrate these identified pathways.Measurements and Main Results: The median (range) lung cavity volume on positron emission tomography/computed tomography scans was 50 cm3 (15-389 cm3). RNA sequence reads (31% splice variants) mapped to 19,049 annotated human genes. Multiple proinflammatory pathways were upregulated in the cavity wall, whereas a downregulation "sink" in the central caseum-fluid interface characterized 53% of pathways including neuroendocrine signaling, calcium signaling, triggering receptor expressed on myeloid cells-1, reactive oxygen and nitrogen species production, retinoic acid-mediated apoptosis, and RIG-I-like receptor signaling. The mathematical model demonstrated that neuroendocrine, protein kinase C-θ, and triggering receptor expressed on myeloid cells-1 pathways, and macrophage and neutrophil numbers, had the highest correlation with bacterial burden (r > 0.6), whereas T-helper effector systems did not.Conclusions: These data provide novel insights into host immunity to Mycobacterium tuberculosis-related cavitation. The pathways defined may serve as useful targets for the design of host-directed therapies, and transmission prevention interventions.

Microinvasion by Streptococcus pneumoniae induces epithelial innate immunity during colonisation at the human mucosal surface

Control of Streptococcus pneumoniae colonisation at human mucosal surfaces is critical to reducing the burden of pneumonia and invasive pneumococcal disease, interrupting transmission, and achieving herd protection. Here, we use an experimental human pneumococcal carriage model (EHPC) to show that S. pneumoniae colonisation is associated with epithelial surface adherence, micro-colony formation and invasion, without overt disease. Interactions between different strains and the epithelium shaped the host transcriptomic response in vitro. Using epithelial modules from a human epithelial cell model that recapitulates our in vivo findings, comprising of innate signalling and regulatory pathways, inflammatory mediators, cellular metabolism and stress response genes, we find that inflammation in the EHPC model is most prominent around the time of bacterial clearance. Our results indicate that, rather than being confined to the epithelial surface and the overlying mucus layer, the pneumococcus undergoes micro-invasion of the epithelium that enhances inflammatory and innate immune responses associated with clearance.

Streptococcus pneumoniae is a common coloniser of the human nasopharynx, but it also causes severe diseases. Here, Weight et al. use an experimental human pneumococcal carriage model to show that bacterial colonisation is associated with invasion of the epithelium and enhancement of immune responses.

Immune Cell Gene Signatures for Profiling the Microenvironment of Solid Tumors

The immune composition of the tumor microenvironment regulates processes including angiogenesis, metastasis, and the response to drugs or immunotherapy. To facilitate the characterization of the immune component of tumors from transcriptomics data, a number of immune cell transcriptome signatures have been reported that are made up of lists of marker genes indicative of the presence a given immune cell population. The majority of these gene signatures have been defined through analysis of isolated blood cells. However, blood cells do not reflect the differentiation or activation state of similar cells within tissues, including tumors, and consequently markers derived from blood cells do not necessarily transfer well to tissues. To address this issue, we generated a set of immune gene signatures derived directly from tissue transcriptomics data using a network-based deconvolution approach. We define markers for seven immune cell types, collectively named ImSig, and demonstrate how these markers can be used for the quantitative estimation of the immune cell content of tumor and nontumor tissue samples. The utility of ImSig is demonstrated through the stratification of melanoma patients into subgroups of prognostic significance and the identification of immune cells with the use of single-cell RNA-sequencing data derived from tumors. Use of ImSig is facilitated by an R package (imsig). Cancer Immunol Res; 6(11); 1388-400. ©2018 AACR.

Tumor Necrosis Factor (TNF) Bioactivity at the Site of an Acute Cell-Mediated Immune Response Is Preserved in Rheumatoid Arthritis Patients Responding to Anti-TNF Therapy

The impact of anti-tumor necrosis factor (TNF) therapies on inducible TNF-dependent activity in humans has never been evaluated in vivo. We aimed to test the hypothesis that patients responding to anti-TNF treatments exhibit attenuated TNF-dependent immune responses at the site of an immune challenge. We developed and validated four context-specific TNF-inducible transcriptional signatures to quantify TNF bioactivity in transcriptomic data. In anti-TNF treated rheumatoid arthritis (RA) patients, we measured the expression of these biosignatures in blood, and in skin biopsies from the site of tuberculin skin tests (TSTs) as a human experimental model of multivariate cell-mediated immune responses. In blood, anti-TNF therapies attenuated TNF bioactivity following ex vivo stimulation. However, at the site of the TST, TNF-inducible gene expression and genome-wide transcriptional changes associated with cell-mediated immune responses were comparable to that of RA patients receiving methotrexate only. These data demonstrate that anti-TNF agents in RA patients do not inhibit inducible TNF activity at the site of an acute inflammatory challenge in vivo, as modeled by the TST. We hypothesize instead that their therapeutic effects are limited to regulating TNF activity in chronic inflammation or by alternative non-canonical pathways.