scholarly journals Cerebrospinal fluid cells immune landscape in multiple sclerosis

2021 ◽  
Vol 19 (1) ◽  
Author(s):  
Zijian Li ◽  
Yongchao Liu ◽  
Aili Jia ◽  
Yueran Cui ◽  
Juan Feng

Abstract Background Multiple Sclerosis (MS) is a potentially devastating autoimmune neurological disorder, which characteristically induces demyelination of white matter in the brain and spinal cord. Methods In this study, three characteristics of the central nervous system (CNS) immune microenvironment occurring during MS onset were explored; immune cell proportion alteration, differential gene expression profile, and related pathways. The raw data of two independent datasets were obtained from the ArrayExpress database; E-MTAB-69, which was used as a derivation cohort, and E-MTAB-2374 which was used as a validation cohort. Differentially expressed genes (DEGs) were identified by the false discovery rate (FDR) value of < 0.05 and |log2 (Fold Change)|> 1, for further analysis. Then, functional enrichment analyses were performed to explore the pathways associated with MS onset. The gene expression profiles were analyzed using CIBERSORT to identify the immune type alterations involved in MS disease. Results After verification, the proportion of five types of immune cells (plasma cells, monocytes, macrophage M2, neutrophils and eosinophils) in cerebrospinal fluid (CSF) were revealed to be significantly altered in MS cases compared to the control group. Thus, the complement and coagulation cascades and the systemic lupus erythematosus (SLE) pathways may play critical roles in MS. We identified NLRP3, LILRB2, C1QB, CD86, C1QA, CSF1R, IL1B and TLR2 as eight core genes correlated with MS. Conclusions Our study identified the change in the CNS immune microenvironment of MS cases by analysis of the in silico data using CIBERSORT. Our data may assist in providing directions for further research as to the molecular mechanisms of MS and provide future potential therapeutic targets in treatment.

2021 ◽  
Vol 9 (Suppl 3) ◽  
pp. A954-A955
Author(s):  
Jacob Kaufman ◽  
Doug Cress ◽  
Theresa Boyle ◽  
David Carbone ◽  
Neal Ready ◽  
...  

BackgroundLKB1 (STK11) is a commonly disrupted tumor suppressor in NSCLC. Its loss promotes an immune exclusion phenotype with evidence of low expression of interferon stimulated genes (ISG) and decreased microenvironment immune infiltration.1 2 Clinically, LKB1 loss induces primary immunotherapy resistance.3 LKB1 is a master regulator of a complex downstream kinase network and has pleiotropic effects on cell biology. Understanding the heterogeneous phenotypes associated with LKB1 loss and their influence on tumor-immune biology will help define and overcome mechanisms of immunotherapy resistance within this subset of lung cancer.MethodsWe applied multi-omic analyses across multiple lung adenocarcinoma datasets2 4–6 (>1000 tumors) to define transcriptional and genetic features enriched in LKB1-deficient lung cancer. Top scoring phenotypes exhibited heterogeneity across LKB1-loss tumors, and were further interrogated to determine association with increased or decreased markers of immune activity. Further, immune cell-types were estimated by Cibersort to identify effects of LKB1 loss on the immune microenvironment. Key conclusions were confirmed by blinded pathology review.ResultsWe show that LKB1 loss significantly affects differentiation patterns, with enrichment of ASCL1-expressing tumors with putative neuroendocrine differentiation. LKB1-deficient neuroendocrine tumors had lower expression of Interferon Stimulated Genes (ISG), MHC1 and MHC2 components, and immune infiltration compared to LKB1-WT and non-neuroendocrine LKB1-deficient tumors (figure 1).The abundances of 22 immune cell types assessed by Cibersort were compared between LKB1-deficient and LKB1-WT tumors. We observe skewing of immune microenvironmental composition by LKB1 loss, with lower abundance of dendritic cells, monocytes, and macrophages, and increased levels of neutrophils and plasma cells (table 1). These trends were most pronounced among tumors with neuroendocrine differentiation, and were concordant across three independent datasets. In a confirmatory subset of 20 tumors, plasma cell abundance was assessed by a blinded pathologist. Pathologist assessment was 100% concordant with Cibersort prediction, and association with LKB1 loss was confirmed (P=0.001).Abstract 909 Figure 1Immune-associated Gene Expression Profiles Affected by Neuroendocrine Differentiation within LKB1-Deficient Lung Adenocarcinomas. Gene expression profiles corresponding to five immune-associated phenotypes are shown with bars indicating average GEP scores for tumors grouped according to LKB1 and neuroendocrine status as indicated. P-values represent results from Student’s T-test between groups as indicated.Abstract 909 Table 1LKB1 Loss Affects Composition of Immune Microenvironment. Values indicate log10 P-values comparing LKB1-loss to LKB1-WT tumors. Positive (red) indicates increased abundance in LKB1 loss. Negative (blue) indicates decreased abundance.ConclusionsWe conclude that tumor differentiation patterns strongly influence the immune microenvironment and immune exclusion characteristics of LKB1-deficient tumors. Neuroendocrine differentiation is associated with the strongest immune exclusion characteristics and should be evaluated clinically for evidence of immunotherapy resistance. A novel observation of increased plasma cell abundance is observed across multiple datasets and confirmed by pathology. Causal mechanisms linking differentiation status to immune activity is not well understood, and the functional role of plasma cells in the immune biology of LKB1-deficient tumors is undefined. These questions warrant further study to inform precision immuno-oncology treatments for these patients.AcknowledgementsThis work was funded by SITC AZ Immunotherapy in Lung Cancer grant (SPS256666) and DOD Lung Cancer Research Program Concept Award (LC180633).ReferencesSkoulidis F, Byers LA, Diao L, et al. Co-occurring genomic alterations define major subsets of KRAS-mutant lung adenocarcinoma with distinct biology, immune profiles, and therapeutic vulnerabilities. Cancer Discov 2015;5:860–77.Schabath MB, Welsh EA, Fulp WJ, et al. Differential association of STK11 and TP53 with KRAS mutation-associated gene expression, proliferation and immune surveillance in lung adenocarcinoma. Oncogene 2016;35:3209–16.Skoulidis F, Goldberg ME, Greenawalt DM, et al. STK11/LKB1 mutations and PD-1 inhibitor resistance in KRAS-mutant lung adenocarcinoma. Cancer Discovery 2018;8:822-835.Cancer Genome Atlas Research Network. Comprehensive molecular profiling of lung adenocarcinoma. Nature 2014;511:543–50.Chitale D, Gong Y, Taylor BS, et al. An integrated genomic analysis of lung cancer reveals loss of DUSP4 in EGFR-mutant tumors. Oncogene 2009;28:2773–83.Shedden K, Taylor JM, Enkemann SA, et al. Gene expression-based survival prediction in lung adenocarcinoma: a multi-site, blinded validation study. Nat Med 2008;14:822–7.


2019 ◽  
Author(s):  
Yiling Cao ◽  
Weihao Tang ◽  
Wanxin Tang

Abstract Objects Lupus nephritis (LN) is a common complication of systemic lupus erythematosus that presents a high risk of end-stage renal disease. In the present study, we used CIBERSORT and gene set enrichment analysis (GSEA) of gene expression profiles to identify immune cell infiltration characteristics and related core genes in LN. Methods Datasets from the Gene Expression Omnibus, GSE32591 and GSE113342, were downloaded for further analysis. The GSE32591 dataset, which included 32 LN glomerular biopsy tissues and 14 glomerular tissues of living donors, was analyzed by CIBERSORT. Different immune cell types in LN were analyzed by the Limma software. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis based on GSEA were performed by clusterProfiler software. Lists of core genes were derived from Spearman correlation between the most significant GO term and differentially expressed immune cell gene from CIBERSORT. GSE113342 was employed to validate the association between selected core genes and clinical manifestation. Result Five types of immune cells revealed important associations with LN, and monocytes emerged to be the prominent differences. GO and KEGG analyses indicated that immune response pathways are significantly enriched in LN. The Spearman correlation indicated that 15 genes, including FCER1G, CLEC7A, MARCO, CLEC7A, PSMB9, and PSMB8, were closely related to clinical features. Conclusion This study is the first to identify immune cell infiltration with microarray data of glomeruli in LN by using CIBERSORT analysis and provides novel evidence and clues for further research of the molecular mechanisms of LN.


2020 ◽  
Author(s):  
Hui Xie ◽  
Xiao-hui Ding ◽  
Ce Yuan ◽  
Jin-jiang Li ◽  
Zhao-yang Li ◽  
...  

Abstract Background: To identify candidate key genes and pathways related to mast cells resting in meningioma and the underlying molecular mechanisms of meningioma.Methods: Gene expression profiles of GSE43290 and GSE16581 datasets were obtained from the Gene Expression Omnibus (GEO) database. GO and KEGG pathway enrichments of DEGs were analyzed using the ClusterProfiler package in R. The protein-protein interaction network (PPI), and TF-miRNA- mRNA co-expression networks were constructed. Further, the difference in immune infiltration was investigated using the CIBERSORT algorithm.Results: A total of 1499 DEGs were identified between tumor and normal controls. The analysis of the immune cell infiltration landscape showed that the probability of distribution of memory B cells, regulatory T cells (Tregs), and resting mast cells in tumor samples were significantly higher than those in the controls. Moreover, through WGCNA analysis, the module related to mast cells resting contained 158 DEGs, and KEGG pathway analysis revealed that the DEGs were dominant in the TNF signaling pathway, cytokine-cytokine receptor interaction, and IL-17 signaling pathway. Survival analysis of hub genes related to mast cells resting showed that the risk model was constructed based on 9 key genes. The TF-miRNA- mRNA co-regulation network, including MYC-miR-145-5p, TNFAIP3-miR-29c-3p, and TNFAIP3-hsa-miR-335-3p, were obtained. Further, 36 nodes and 197 interactions in the PPI network were identified. Conclusions: The results of this study revealed candidate key genes, miRNAs, and pathways related to mast cells resting involved in meningioma development, providing potential therapeutic targets for meningioma treatment.


Author(s):  
Qingchun Liang ◽  
Qin Zhou ◽  
Jinhe Li ◽  
Zhugui Chen ◽  
Zhihao Zhang ◽  
...  

Abstract Acute lung injury (ALI) is an inflammatory pulmonary disease that can easily develop into serious acute respiratory distress syndrome, which has high morbidity and mortality. However, the molecular mechanism of ALI remains unclear, and few molecular biomarkers for diagnosis and treatment have been identified. In this study, we aimed to identify novel molecular biomarkers using a bioinformatics approach. Gene expression data were obtained from the Gene Expression Omnibus database, co-expressed differentially expressed genes (CoDEGs) were identified using R software, and further functional enrichment analyses were conducted using the online tool Database for Annotation, Visualization, and Integrated Discovery. A protein–protein interaction network was established using the STRING database and Cytoscape software. Lipopolysaccharide (LPS)-induced ALI mouse model was constructed and verified. The hub genes were screened and validated in vivo. The transcription factors (TFs) and miRNAs associated with the hub genes were predicted using the NetworkAnalyst database. In total, 71 CoDEGs were screened and found to be mainly involved in the cytokine–cytokine receptor interactions, and the tumor necrosis factor and malaria signaling pathways. Animal experiments showed that the lung injury score, bronchoalveolar lavage fluid protein concentration, and wet-to-dry weight ratio were higher in the LPS group than those in the control group. Real-time polymerase chain reaction analysis indicated that most of the hub genes such as colony-stimulating factor 2 (Csf2) were overexpressed in the LPS group. A total of 20 TFs including nuclear respiratory factor 1 (NRF1) and two miRNAs were predicted to be regulators of the hub genes. In summary, Csf2 may serve as a novel diagnostic and therapeutic target for ALI. NRF1 and mmu-mir-122-5p may be key regulators in the development of ALI.


2021 ◽  
Vol 2021 ◽  
pp. 1-15
Author(s):  
Shunqiang Nong ◽  
Xiaohao Chen ◽  
Zechen Wang ◽  
Guidan Xu ◽  
Wujun Wei ◽  
...  

Background. Increasing evidence demonstrated that long noncoding RNA (lncRNA) could affect inflammatory tumor immune microenvironment by modulating gene expression and could be used as a biomarker for HBC-related hepatocellular carcinoma (HCC) but still needs further research. The aim of the present study was to determine an lncRNA signature for the diagnosis of HBV-related HCC. Methods. HBV-related HCC expression profiles (GSE55092, GSE19665, and GSE84402) were abstracted from the GEO (Gene Expression Omnibus) data resource, and R package limma and RobustRankAggreg were employed to identify common differentially expressed genes (DEGs). Using machine learning, optimal diagnostic lncRNA molecular markers for HBV-related HCC were identified. The expression of candidate lncRNAs was cross-validated in GSE121248, and an ROC (receiver operating characteristic) curve of lncRNA biomarkers was carried out. Additionally, a coexpression network and functional annotation was built, after which a PPI (protein-protein interaction) network along with module analysis were conducted with the Cytoscape open source software. Result. A total of 38 DElncRNAs and 543 DEmRNAs were identified with a fold change larger than 2.0 and a P value < 0.05. By machine learning, AL356056.2, AL445524.1, TRIM52-AS1, AC093642.1, EHMT2-AS1, AC003991.1, AC008040.1, LINC00844, and LINC01018 were screened out as optional diagnostic lncRNA biosignatures for HBV-related HCC. The AUC (areas under the curve) of the SVM (support vector machine) model and random forest model were 0.957 and 0.904, respectively, and the specificity and sensitivity were 95.7 and 100% and 94.3 and 86.5%, respectively. The results of functional enrichment analysis showed that the integrated coexpressed DEmRNAs shared common cascades in the p53 signaling pathway, retinol metabolism, PI3K-Akt signaling cascade, and chemical carcinogenesis. The integrated DEmRNA PPI network complex was found to be comprised of 87 nodes, and two vital modules with a high degree were selected with the MCODE app. Conclusion. The present study identified nine potential diagnostic biomarkers for HBV-related HCC, all of which could potentially modulated gene expression related to inflammatory conditions in the tumor immune microenvironment. The functional annotation of the target DEmRNAs yielded novel evidence in evaluating the precise functions of lncRNA in HBV-related HCC.


Cells ◽  
2021 ◽  
Vol 10 (12) ◽  
pp. 3484
Author(s):  
Jisun Lim ◽  
Yeon Bi Han ◽  
Soo Young Park ◽  
Soyeon Ahn ◽  
Hyojin Kim ◽  
...  

Many studies support a stepwise continuum of morphologic changes between atypical adenomatous hyperplasia (AAH) and lung adenocarcinoma (ADC). Here we characterized gene expression patterns and the association of differentially expressed genes and immune tumor microenvironment behaviors in AAH to ADC during ADC development. Tumor tissues from nine patients with ADC and synchronous multiple ground glass nodules/lesions (GGN/Ls) were analyzed using RNA sequencing. Using clustering, we identified genes differentially and sequentially expressed in AAH and ADC compared to normal tissues. Functional enrichment analysis using gene ontology terms was performed, and the fraction of immune cell types was estimated. We identified up-regulated genes (ACSL5 and SERINC2) with a stepwise change of expression from AAH to ADC and validated those expressions by quantitative PCR and immunohistochemistry. The immune cell profiles revealed increased B cell activities and decreased natural killer cell activities in AAH and ADC. A stepwise change of differential expression during ADC development revealed potential effects on immune function in synchronous precursors and in tumor lesions in patients with lung cancer.


2020 ◽  
Author(s):  
Junhong Li ◽  
Yang Zhai ◽  
Peng Wu ◽  
Yueqiang Hu ◽  
Wei Chen ◽  
...  

Abstract BACKGROUD: Microarray-based gene expression profiling is widely used in biomedical research. Weighted gene co-expression network analysis (WGCNA) links microarray data directly to clinical traits and identifies rules for predicting pathological stage and prognosis of disease.WGCNA is useful in understandingmany biological processes. Stroke is a common disease worldwide, however, molecular mechanisms of its pathogenesis are largely unknown. The aim of this study was to construct gene co-expression networks for identification of key modules and hub genes associated with stroke pathogenesis.METHODS: Gene microarray expression profiles of stroke samples were retrieved from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) were screened by the limma package in R software. WGCNA was used to construct free-scale gene co-expression networks to explore the associations between gene sets and clinical features, and to identify key modules and hub genes. Subsequently, functional enrichment analyses were performed. Further, receiver operating characteristic (ROC) curve analysis was carried out to validate expression of hub genes and literature validation was performed as well.RESULTS: A total of 11,747 most variant genes were used for co-expression network construction. Pink and yellow modules were significantly correlated to stroke pathogenesis. Functional enrichment analysis showed that the pink module was mainly involved in regulation of neuron regeneration, and repair of DNA damage.On the other hand, yellow module was mainly enriched in ion transport system dysfunction which was correlated with neuron death. A total of eight hub genes (PRR11, NEDD9, Notch2, RUNX1-IT1, ANP32A-IT1, ASTN2, SAMHD1 and STIM1) were identified and validated at transcriptional levels and through existing literature.CONCLUSION: The eight hub genes (PRR11, NEDD9, Notch2, RUNX1-IT1, ANP32A-IT1, ASTN2, SAMHD1 and STIM1) identified in the study are potentialbiomarkers and therapeutic targets for effective diagnosis and treatment of stroke.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Sheyda Khalilian ◽  
Zohreh Hojati ◽  
Fariba Dehghanian ◽  
Vahid Shaygannejad ◽  
Seyedeh Zahra Hosseini Imani ◽  
...  

AbstractAlterations in the regulatory mechanisms that control the process of myelination in the nervous system, may lead to the impaired myelination in the Multiple sclerosis. The Hippo pathway is an important mediator of myelination in the nervous system and might contribute to the pathophysiology of MS. This study examined via qPCR the RNA expression of YAP1, TAZ, and CRB3 as the key effectors of the Hippo pathway and also, VDR in the peripheral blood of 35 sporadic, 37 familial MS patients; and also 34 healthy first-degree relatives of the familial MS patients (HFR) and 40 healthy individuals without a family history of the disease (control). The results showed the increased expression of VDR in the sporadic group, as compared to other groups. There was also an increased expression of TAZ in the familial and HFR groups, as compared to the control group. The familial and sporadic patients displayed a significantly lower level of expression of YAP1 in comparison to the HFR group. The increased expression level in the sporadic patients and control group, as compared to the HFR group, was seen in CRB3. We also assessed different clinical parameters and MRI characteristics of the patients. Overall, these findings suggest that Hippo pathway effectors and also VDR gene may play a potential role in the pathophysiology of the sporadic and familial forms of MS. Confirmation of different gene expression patterns in sporadic and familial MS groups may have obvious implications for the personalization of therapies in the disease.


Author(s):  
Martina Reiter ◽  
Ales Tichopad ◽  
Irmgard Riedmaier ◽  
Michael W. Pfaffl ◽  
Heinrich H.D. Meyer

AbstractThe focus of this study was to evaluate data on the gene expression profiles induced by testosterone and a selective androgen receptor modulator (SARM, TAP Pharmaceutical Products Inc., Lake Forest, IL, USA) in androgen sensitive muscle tissue to obtain a better understanding on the molecular mechanisms of action and to identify biomarkers for SARM function in primate organs. A total of 24 male cyomolgus monkeys were divided into four groups: testosterone group, SARM1 group, SARM10 group, and control group, each consisting of six animals. The testosterone group was treated i.m. with 3.0 mg/kg Testostoviron


PeerJ ◽  
2019 ◽  
Vol 7 ◽  
pp. e8162
Author(s):  
Ying Le

Myelodysplastic syndrome (MDS) is a heterogeneous hematologic malignancy derived from hematopoietic stem cells and the molecular mechanism of MDS remains unclear. This study aimed to elucidate potential markers of diagnosis and prognosis of MDS. The gene expression profiles GSE19429 and GSE58831 were obtained and downloaded from the Gene Expression Omnibus (GEO) database. The differentially expressed genes (DEGs) in MDS were screened using GEO2R and overlapped DEGs were obtained with Venn Diagrams. Then, Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway functional enrichment analyses, protein–protein interaction network establishment and survival analyses were performed. Functional enrichment analysis indicated that these DEGs were significantly enriched in the interferon signaling pathway, immune response, hematopoietic cell lineage and the FOXO signaling pathway. Four hub genes and four significant modules including 25 module genes were obtained via Cytoscape MCODE. Survival analysis showed that the overall survival of MDS patients having BLNK, IRF4, IFITM1, IFIT1, ISG20, IFI44L alterations were worse than that without alterations. In conclusion, the identification of these genes and pathways helps understand the underlying molecular mechanisms of MDS and provides candidate targets for the diagnosis and prognosis of MDS.


Sign in / Sign up

Export Citation Format

Share Document