scholarly journals Integrative transcriptomic analysis of SLE reveals IFN-driven cross-talk between immune cells

2020 ◽  
Author(s):  
Bharat Panwar ◽  
Benjamin J. Schmiedel ◽  
Shu Liang ◽  
Brandie White ◽  
Enrique Rodriguez ◽  
...  

ABSTRACTThe systemic lupus erythematosus (SLE) is an incurable autoimmune disease disproportionately affecting women and may lead to damage in multiple different organs. The marked heterogeneity in its clinical manifestations is a major obstacle in finding targeted treatments and involvement of multiple immune cell types further increases this complexity. Thus, identifying molecular subtypes that best correlate with disease heterogeneity and severity as well as deducing molecular cross-talk among major immune cell types that lead to disease progression are critical steps in the development of more informed therapies for SLE. Here we profile and analyze gene expression of six major circulating immune cell types from patients with well-characterized SLE (classical monocytes (n=64), T cells (n=24), neutrophils (n=24), B cells (n=20), conventional (n=20) and plasmacytoid (n=22) dendritic cells) and from healthy control subjects. Our results show that the interferon (IFN) response signature was the major molecular feature that classified SLE patients into two distinct groups: IFN-signature negative (IFNneg) and positive (IFNpos). We show that the gene expression signature of IFN response was consistent (i) across all immune cell types, (ii) all single cells profiled from three IFNpos donors using single-cell RNA-seq, and (iii) longitudinal samples of the same patient. For a better understanding of molecular differences of IFNpos versus IFNneg patients, we combined differential gene expression analysis with differential Weighted Gene Co-expression Network Analysis (WGCNA), which revealed a relatively small list of genes from classical monocytes including two known immune modulators, one the target of an approved therapeutic for SLE (TNFSF13B/BAFF: belimumab) and one itself a therapeutic for Rheumatoid Arthritis (IL1RN: anakinra). For a more integrative understanding of the cross-talk among different cell types and to identify potentially novel gene or pathway connections, we also developed a novel gene co-expression analysis method for joint analysis of multiple cell types named integrated WGNCA (iWGCNA). This method revealed an interesting cross-talk between T and B cells highlighted by a significant enrichment in the expression of known markers of T follicular helper cells (Tfh), which also correlate with disease severity in the context of IFNpos patients. Interestingly, higher expression of BAFF from all myeloid cells also shows a strong correlation with enrichment in the expression of genes in T cells that may mark circulating Tfh cells or related memory cell populations. These cell types have been shown to promote B cell class-switching and antibody production, which are well-characterized in SLE patients. In summary, we generated a large-scale gene expression dataset from sorted immune cell populations and present a novel computational approach to analyze such data in an integrative fashion in the context of an autoimmune disease. Our results reveal the power of a hypothesis-free and data-driven approach to discover drug targets and reveal novel cross-talk among multiple immune cell types specific to a subset of SLE patients. This approach is immediately useful for studying autoimmune diseases and is applicable in other contexts where gene expression profiling is possible from multiple cell types within the same tissue compartment.

2019 ◽  
Author(s):  
Anne-Marie Madore ◽  
Lucile Pain ◽  
Anne-Marie Boucher-Lafleur ◽  
Jolyane Meloche ◽  
Andréanne Morin ◽  
...  

AbstractBackgroundThe 17q12-21 locus is the most replicated association with asthma. However, no study had described the genetic mechanisms underlying this association considering all genes of the locus in immune cell samples isolated from asthmatic and non-asthmatic individuals.ObjectiveThis study takes benefit of samples from naïve CD4+ T cells and eosinophils isolated from the same 200 individuals to describe specific interactions between genetic variants, gene expression and DNA methylation levels for the 17q12-21 asthma locus.Methods and ResultsAfter isolation of naïve CD4+ T cells and eosinophils from blood samples, next generation sequencing was used to measure DNA methylation levels and gene expression counts. Genetic interactions were then evaluated considering genetic variants from imputed genotype data. In naïve CD4+ T cells but not eosinophils, 20 SNPs in the fourth and fifth haplotype blocks modulated both GSDMA expression and methylation levels, showing an opposite pattern of allele frequencies and expression counts in asthmatics compared to controls. Moreover, negative correlations have been measured between methylation levels of CpG sites located within the 1.5 kb region from the transcription start site of GSDMA and its expression counts.ConclusionAvailability of sequencing data from two key cell types isolated from asthmatic and non-asthmatic individuals allowed identifying a new gene in naïve CD4+ T cells that drives the association with the 17q12-21 locus, leading to a better understanding of the genetic mechanisms taking place in it.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Ikko Mito ◽  
Hideyuki Takahashi ◽  
Reika Kawabata-Iwakawa ◽  
Shota Ida ◽  
Hiroe Tada ◽  
...  

AbstractHead and neck squamous carcinoma (HNSCC) is highly infiltrated by immune cells, including tumor-infiltrating lymphocytes and myeloid lineage cells. In the tumor microenvironment, tumor cells orchestrate a highly immunosuppressive microenvironment by secreting immunosuppressive mediators, expressing immune checkpoint ligands, and downregulating human leukocyte antigen expression. In the present study, we aimed to comprehensively profile the immune microenvironment of HNSCC using gene expression data obtained from public database. We calculated enrichment scores of 33 immune cell types based on gene expression data of HNSCC tissues and adjacent non-cancer tissues. Based on these scores, we performed non-supervised clustering and identified three immune signatures—cold, lymphocyte, and myeloid/dendritic cell (DC)—based on the clustering results. We then compared the clinical and biological features of the three signatures. Among HNSCC and non-cancer tissues, human papillomavirus (HPV)-positive HNSCCs exhibited the highest scores in various immune cell types, including CD4+ T cells, CD8+ T cells, B cells, plasma cells, basophils, and their subpopulations. Among the three immune signatures, the proportions of HPV-positive tumors, oropharyngeal cancers, early T tumors, and N factor positive cases were significantly higher in the lymphocyte signature than in other signatures. Among the three signatures, the lymphocyte signature showed the longest overall survival (OS), especially in HPV-positive patients, whereas the myeloid/DC signature demonstrated the shortest OS in these patients. Gene set enrichment analysis revealed the upregulation of several pathways related to inflammatory and proinflammatory responses in the lymphocyte signature. The expression of PRF1, IFNG, GZMB, CXCL9, CXCL10, PDCD1, LAG3, CTLA4, HAVCR2, and TIGIT was the highest in the lymphocyte signature. Meanwhile, the expression of PD-1 ligand genes CD274 and PDCD1LG2 was highest in the myeloid/DC signature. Herein, our findings revealed the transcriptomic landscape of the immune microenvironment that closely reflects the clinical and biological significance of HNSCC, indicating that molecular profiling of the immune microenvironment can be employed to develop novel biomarkers and precision immunotherapies for HNSCC.


2017 ◽  
Vol 313 (5) ◽  
pp. E528-E539 ◽  
Author(s):  
Katya B. Rubinow ◽  
Jing H. Chao ◽  
Derek Hagman ◽  
Mario Kratz ◽  
Brian Van Yserloo ◽  
...  

Male hypogonadism results in changes in body composition characterized by increases in fat mass. Resident immune cells influence energy metabolism in adipose tissue and could promote increased adiposity through paracrine effects. We hypothesized that manipulation of circulating sex steroid levels in healthy men would alter adipose tissue immune cell populations. Subjects ( n = 44 men, 19–55 yr of age) received 4 wk of treatment with the gonadotropin-releasing hormone receptor antagonist acyline with daily administration of 1) placebo gel, 2) 1.25 g testosterone gel (1.62%), 3) 5 g testosterone gel, or 4) 5 g testosterone gel with an aromatase inhibitor. Subcutaneous adipose tissue biopsies were performed at baseline and end-of-treatment, and adipose tissue immune cells, gene expression, and intra-adipose estrogen levels were quantified. Change in serum total testosterone level correlated inversely with change in the number of CD3+ (β = −0.36, P = 0.04), CD4+ (β = −0.34, P = 0.04), and CD8+ (β = −0.33, P = 0.05) T cells within adipose tissue. Change in serum 17β-estradiol level correlated inversely with change in the number of adipose tissue macrophages (ATMs) (β = −0.34, P = 0.05). A negative association also was found between change in serum testosterone and change in CD11c+ ATMs (β = −0.41, P = 0.01). Overall, sex steroid deprivation was associated with increases in adipose tissue T cells and ATMs. No associations were found between changes in serum sex steroid levels and changes in adipose tissue gene expression. Circulating sex steroid levels may regulate adipose tissue immune cell populations. These exploratory findings highlight a possible novel mechanism that could contribute to increased metabolic risk in hypogonadal men.


2019 ◽  
Vol 9 (1) ◽  
Author(s):  
Konstantin Carlberg ◽  
Marina Korotkova ◽  
Ludvig Larsson ◽  
Anca I. Catrina ◽  
Patrik L. Ståhl ◽  
...  

AbstractLately it has become possible to analyze transcriptomic profiles in tissue sections with retained cellular context. We aimed to explore synovial biopsies from rheumatoid arthritis (RA) and spondyloarthritis (SpA) patients, using Spatial Transcriptomics (ST) as a proof of principle approach for unbiased mRNA studies at the site of inflammation in these chronic inflammatory diseases. Synovial tissue biopsies from affected joints were studied with ST. The transcriptome data was subjected to differential gene expression analysis (DEA), pathway analysis, immune cell type identification using Xcell analysis and validation with immunohistochemistry (IHC). The ST technology allows selective analyses on areas of interest, thus we analyzed morphologically distinct areas of mononuclear cell infiltrates. The top differentially expressed genes revealed an adaptive immune response profile and T-B cell interactions in RA, while in SpA, the profiles implicate functions associated with tissue repair. With spatially resolved gene expression data, overlaid on high-resolution histological images, we digitally portrayed pre-selected cell types in silico. The RA displayed an overrepresentation of central memory T cells, while in SpA effector memory T cells were most prominent. Consequently, ST allows for deeper understanding of cellular mechanisms and diversity in tissues from chronic inflammatory diseases.


2020 ◽  
Vol 22 (Supplement_3) ◽  
pp. iii312-iii312
Author(s):  
Timothy Ritzmann ◽  
Anbarasu Lourdusamy ◽  
Andrew Jackson ◽  
Lisa Storer ◽  
Andrew Donson ◽  
...  

Abstract Ependymoma is the third commonest childhood brain tumour. Relapse is frequent, often fatal and current therapeutic strategies are inadequate. Previous ependymoma research describes an immunosuppressive environment with T-cell exhaustion, indicating a lack of response to T-cell directed immunotherapy. Understanding the immune microenvironment is therefore critical. We present a computational analysis of ependymoma, gene expression derived, immune profiles. Using 465 ependymoma samples from gene expression datasets (GSE64415, GSE50385, GSE100240) and two RNA-seq databases from UK ependymomas, we applied bulk tumour deconvolution methods (CIBERSORT and xCell) to infer immune cell populations. Additionally, we measured checkpoint blockade related mRNAs and used immunohistochemistry to investigate cell populations in ependymoma sections. CIBERSORT indicated high proportions of M2-like macrophages and smaller proportions of activated natural killer (NK) cells, T follicular helper cells, CD4+ memory T-cells and B-cells. xCell overlapped with the M2-like macrophage and CD4+ memory T-cell signatures seen in CIBERSORT. On immunohistochemistry, T and B cells were scarce, with small numbers of CD8+, CD4+ and CD20+ cells in the parenchyma but greater numbers in surrounding regions. CD68 was more highly expressed in the parenchyma. Analysis of nine checkpoint ligands and receptors demonstrated only the TIM3/GAL9 combination was reliably detectable. GAL9 is implicated in tumour interactions with T-cells and macrophages elsewhere, possibly contributing to poorer outcomes. Our study supports the presence of myeloid cells being leading contributors to the ependymoma immune microenvironment. Further work will delineate the extent of myeloid contribution to immunosuppression across molecular subtypes. Modulation of tumour immunity may contribute to better clinical outcomes.


Cancers ◽  
2021 ◽  
Vol 13 (15) ◽  
pp. 3854
Author(s):  
Saira Khalique ◽  
Sarah Nash ◽  
David Mansfield ◽  
Julian Wampfler ◽  
Ayoma Attygale ◽  
...  

Ovarian clear cell carcinoma (OCCC) is a rare subtype of epithelial ovarian cancer characterised by a high frequency of loss-of-function ARID1A mutations and a poor response to chemotherapy. Despite their generally low mutational burden, an intratumoural T cell response has been reported in a subset of OCCC, with ARID1A purported to be a biomarker for the response to the immune checkpoint blockade independent of micro-satellite instability (MSI). However, assessment of the different immune cell types and spatial distribution specifically within OCCC patients has not been described to date. Here, we characterised the immune landscape of OCCC by profiling a cohort of 33 microsatellite stable OCCCs at the genomic, gene expression and histological level using targeted sequencing, gene expression profiling using the NanoString targeted immune panel, and multiplex immunofluorescence to assess the spatial distribution and abundance of immune cell populations at the protein level. Analysis of these tumours and subsequent independent validation identified an immune-related gene expression signature associated with risk of recurrence of OCCC. Whilst histological quantification of tumour-infiltrating lymphocytes (TIL, Salgado scoring) showed no association with the risk of recurrence or ARID1A mutational status, the characterisation of TILs via multiplexed immunofluorescence identified spatial differences in immunosuppressive cell populations in OCCC. Tumour-associated macrophages (TAM) and regulatory T cells were excluded from the vicinity of tumour cells in low-risk patients, suggesting that high-risk patients have a more immunosuppressive microenvironment. We also found that TAMs and cytotoxic T cells were also excluded from the vicinity of tumour cells in ARID1A-mutated OCCCs compared to ARID1A wild-type tumours, suggesting that the exclusion of these immune effectors could determine the host response of ARID1A-mutant OCCCs to therapy. Overall, our study has provided new insights into the immune landscape and prognostic associations in OCCC and suggest that tailored immunotherapeutic approaches may be warranted for different subgroups of OCCC patients.


Author(s):  
Jiaqi Li ◽  
Lifang Li ◽  
Yimeng Wang ◽  
Gan Huang ◽  
Xia Li ◽  
...  

To date, nearly 100 autoimmune diseases have been an area of focus, and these diseases bring health challenges to approximately 5% of the population worldwide. As a type of disease caused by tolerance breakdown, both environmental and genetic risk factors contribute to autoimmune disease development. However, in most cases, there are still gaps in our understanding of disease pathogenesis, diagnosis, and treatment. Therefore, more detailed knowledge of disease pathogenesis and potential therapies is indispensable. DNA methylation, which does not affect the DNA sequence, is one of the key epigenetic silencing mechanisms and has been indicated to play a key role in gene expression regulation and to participate in the development of certain autoimmune diseases. Potential epigenetic regulation via DNA methylation has garnered more attention as a disease biomarker in recent years. In this review, we clarify the basic function and distribution of DNA methylation, evaluate its effects on gene expression and discuss related key enzymes. In addition, we summarize recent aberrant DNA methylation modifications identified in the most important cell types related to several autoimmune diseases and then provide potential directions for better diagnosing and monitoring disease progression driven by epigenetic control, which may broaden our understanding and contribute to further epigenetic research in autoimmune diseases.


2021 ◽  
Vol 12 (1) ◽  
Author(s):  
Benjamin J. Schmiedel ◽  
Job Rocha ◽  
Cristian Gonzalez-Colin ◽  
Sourya Bhattacharyya ◽  
Ariel Madrigal ◽  
...  

AbstractCommon genetic polymorphisms associated with COVID-19 illness can be utilized for discovering molecular pathways and cell types driving disease pathogenesis. Given the importance of immune cells in the pathogenesis of COVID-19 illness, here we assessed the effects of COVID-19-risk variants on gene expression in a wide range of immune cell types. Transcriptome-wide association study and colocalization analysis revealed putative causal genes and the specific immune cell types where gene expression is most influenced by COVID-19-risk variants. Notable examples include OAS1 in non-classical monocytes, DTX1 in B cells, IL10RB in NK cells, CXCR6 in follicular helper T cells, CCR9 in regulatory T cells and ARL17A in TH2 cells. By analysis of transposase accessible chromatin and H3K27ac-based chromatin-interaction maps of immune cell types, we prioritized potentially functional COVID-19-risk variants. Our study highlights the potential of COVID-19 genetic risk variants to impact the function of diverse immune cell types and influence severe disease manifestations.


2016 ◽  
Author(s):  
Santiago J. Carmona ◽  
Sarah A. Teichmann ◽  
Lauren Ferreira ◽  
Iain C. Macaulay ◽  
Michael J.T. Stubbington ◽  
...  

AbstractThe immune system of vertebrate species consists of many different cell types that have distinct functional roles and are subject to different evolutionary pressures. Here, we first analysed gene conservation of all major immune cell types in human and mouse. Our results revealed higher gene turnover and faster evolution of trans-membrane proteins in NK cells compared to other immune cell populations, and especially T cells, but similar conservation of nuclear and cytoplasmic protein coding genes. To validate these findings in a distant vertebrate species, we used single-cell RNA-Sequencing of lck:GFP cells in zebrafish to obtain the first transcriptome of specific immune cell types in a non-mammalian species. Unsupervised clustering and single-cell TCR locus reconstruction identified three cell populations, T-cells, a novel type of NK-like cells and a smaller population of myeloid-like cells. Differential expression analysis uncovered new immune cell specific genes, including novel immunoglobulin-like receptors, and neofunctionalization of recently duplicated paralogs. Evolutionary analyses confirmed a higher gene turnover and lower conservation of trans-membrane proteins in NK cells compared to T cells in fish species, suggesting that this is a general property of immune cell types across all vertebrates.


Sign in / Sign up

Export Citation Format

Share Document