scholarly journals Integrative Analysis of Bulk and Single-Cell RNA Sequencing Data Reveals Cell Types Involved in Heart Failure

Author(s):  
Xin Shi ◽  
Li Zhang ◽  
Yi Li ◽  
Jieyuan Xue ◽  
Feng Liang ◽  
...  

Owing to the high mortality rates of heart failure (HF), a more detailed description of the HF becomes extremely urgent. Since the pathogenesis of HF remain elusive, a thorough identification of the genetic factors will provide novel insights into the molecular basis of this cardiac dysfunction. In our research, we performed publicly available transcriptome profiling datasets, including non-failure (NF), dilated cardiomyopathy (DCM) and ischemic cardiomyopathy (ICM) hearts tissues. Through principal component analysis (PCA), gene differential expression analysis, gene set enrichment analysis (GSEA), and gene Set Variation Analysis (GSVA), we figured out the candidate genes noticeably altered in HF, the specific biomarkers of endothelial cell (EC) and cardiac fibrosis, then validated the differences of the inflammation-related cell adhesion molecules (CAMs), extracellular matrix (ECM) genes, and immune responses. Taken together, our results suggested the EC and fibroblast could be activated in response to HF. DCM and ICM had both commonality and specificity in the pathogenesis of HF. Higher inflammation in ICM might related to autocrine CCL3/CCL4-CCR5 interaction induced chemokine signaling activation. Furthermore, the activities of neutrophil and macrophage were higher in ICM than DCM. These findings identified features of the landscape of previously underestimated cellular, transcriptomic heterogeneity between ICM and DCM.

F1000Research ◽  
2017 ◽  
Vol 6 ◽  
pp. 2010 ◽  
Author(s):  
Monther Alhamdoosh ◽  
Charity W. Law ◽  
Luyi Tian ◽  
Julie M. Sheridan ◽  
Milica Ng ◽  
...  

Gene set enrichment analysis is a popular approach for prioritising the biological processes perturbed in genomic datasets. The Bioconductor project hosts over 80 software packages capable of gene set analysis. Most of these packages search for enriched signatures amongst differentially regulated genes to reveal higher level biological themes that may be missed when focusing only on evidence from individual genes. With so many different methods on offer, choosing the best algorithm and visualization approach can be challenging. The EGSEA package solves this problem by combining results from up to 12 prominent gene set testing algorithms to obtain a consensus ranking of biologically relevant results.This workflow demonstrates how EGSEA can extend limma-based differential expression analyses for RNA-seq and microarray data using experiments that profile 3 distinct cell populations important for studying the origins of breast cancer. Following data normalization and set-up of an appropriate linear model for differential expression analysis, EGSEA builds gene signature specific indexes that link a wide range of mouse or human gene set collections obtained from MSigDB, GeneSetDB and KEGG to the gene expression data being investigated. EGSEA is then configured and the ensemble enrichment analysis run, returning an object that can be queried using several S4 methods for ranking gene sets and visualizing results via heatmaps, KEGG pathway views, GO graphs, scatter plots and bar plots. Finally, an HTML report that combines these displays can fast-track the sharing of results with collaborators, and thus expedite downstream biological validation. EGSEA is simple to use and can be easily integrated with existing gene expression analysis pipelines for both human and mouse data.


2021 ◽  
Author(s):  
HUA HUANG ◽  
Shanshan Xu ◽  
Youran Li ◽  
Yunfei Gu ◽  
Lijiang Ji

Abstract Background: Colorectal cancer (CRC), the commonly seen malignancy, ranks the 3rd place among the causes of cancer-associated mortality. As suggested by more and more studies, long coding RNAs (lncRNAs) have been considered as prognostic biomarkers for CRC. But the significance of hypoxic lncRNAs in predicting CRC prognosis remains unclear.Methods: The gene expressed profiles for CRC cases were obtained based on the Cancer Genome Atlas (TCGA) and applied to estimate the hypoxia score using a single-sample gene set enrichment analysis (ssGSEA) algorithm. Overall survival (OS) of high- and low-hypoxia score group was analyzed by the Kaplan–Meier (KM) plot. To identify differentially expressed lncRNAs (DELs) between two hypoxia score groups, this study carried out differential expression analysis, and then further integrated with the DELs between controls and CRC patients to generate the hypoxia-related lncRNAs for CRC. Besides, prognostic lncRNAs were screened by the univariate Cox regression, which were later utilized for constructing the prognosis nomogram for CRC by adopting the least absolute shrinkage and selection operator (LASSO) algorithm. In addition, both accuracy and specificity of the constructed prognostic signature were detected through the receiver operating characteristic (ROC) analysis. Moreover, our constructed prognosis signature also was validated in the internal testing test. This study operated gene set enrichment analysis (GSEA) for exploring potential biological functions associated with the prognostic signature. Finally, the ceRNA network of the prognostic lncRNAs was constructed.Results: Among 2299 hypoxia-related lncRNAs of CRC in total, LINC00327, LINC00163, LINC00174, SYNPR-AS1, and MIR31HG were identified as prognostic lncRNAs by the univariate Cox regression, and adopted for constructing the prognosis signature for CRC. ROC analysis showed the predictive power and accuracy of the prognostic signature. Additionally, the GSEA revealed that ECM-receptor interaction, PI3K-Akt pathway, phagosome, and Hippo pathway were mostly associated with the high-risk group. 352 miRNAs-mRNAs pairs and 177 lncRNAs-miRNAs were predicted.Conclusion: To conclude , we identified 5 hypoxia-related lncRNAs to establish an accurate prognostic signature for CRC, providing important prognostic markers and therapeutic target.


2021 ◽  
Vol 2021 ◽  
pp. 1-17
Author(s):  
Fei Liu ◽  
Xiaopeng Yu ◽  
Guijin He

Background. We analyzed the n6-methyladenosine (m6A) modification patterns of immune cells infiltrating the tumor microenvironment of breast cancer (BC) to provide a new perspective for the early diagnosis and treatment of BC. Methods. Based on 23 m6A regulatory factors, we identified m6A-related gene characteristics and m6A modification patterns in BC through unsupervised cluster analysis. To examine the differences in biological processes among various m6A modification modes, we performed genomic variation analysis. We then quantified the relative infiltration levels of different immune cell subpopulations in the tumor microenvironment of BC using the CIBERSORT algorithm and single-sample gene set enrichment analysis. Univariate Cox analysis was used to screen for m6A characteristic genes related to prognosis. Finally, we evaluated the m6A modification pattern of patients with a single BC by constructing the m6Ascore based on principal component analysis. Results. We identified three different m6A modification patterns in 2128 BC samples. A higher abundance of the immune infiltration of the m6Acluster C was indicated by the results of CIBERSORT and the single-sample gene set enrichment analysis. Based on the m6A characteristic genes obtained through screening, the m6Ascore was determined. The BC patients were segregated into m6Ascore groups of low and high categories, which revealed significant survival benefits among patients with low m6Ascores. Additionally, the high-m6Ascore group had a higher mutation frequency and was associated with low PD-L1 expression, and the m6Ascore and tumor mutation burden showed a positive correlation. In addition, treatment effects were better in patients in the high-m6Ascore group. Conclusions. In case of a single patient with BC, the immune cell infiltration characteristics of the tumor microenvironment and the m6A methylation modification pattern could be evaluated using the m6Ascore. Our results provide a foundation for improving personalized immunotherapy of BC.


2021 ◽  
Author(s):  
Rana Alghamdi ◽  
Maryam Al-Zahrani

Abstract Background: Claudin’s gene are associated with various aberrant physiological and cellular signaling. However, the association of claudins with survival prognosis, signaling pathways, and diagnostic efficacy in colon cancer remain lacking. Methods: We used various bioinformatics methods, including differential expression analysis, gene set enrichment analysis (GSEA), protein-protein interaction (PPI), survival analysis, single sample gene set enrichment analysis (ssGSEA), mutation analysis, and identifying receiver operating characteristic (ROC) curve of claudins in the TCGA colon adenocarcinoma (COAD). Results: We found that: CLDN2, CLDN1, CLDN14, CLDN16, CLDN18, CLDN9, CLDN12, and CLDN6 are elevated in COAD. In contrast, the CLDN8, CLDN23, CLDN5, CLDN11, CLDN7, and CLDN15 are downregulated in COAD. Various claudin’s genes are mutated and associated with diagnostic efficacy in the COAD. Conclusions: Claudin’s genes are associated with prognosis, immune regulation, signaling pathway regulations, and diagnosis. These findings may provide new molecular insight into the treatment of colon cancer.


2021 ◽  
Author(s):  
Hua Huang ◽  
Mingjia Gu ◽  
Shanshan Xu ◽  
Youran Li ◽  
Yunfei Gu ◽  
...  

Abstract Background:Colorectal cancer (CRC), the commonly seen malignancy, ranks 3rd place among the causes of cancer-associated mortality. As suggested by more and more studies, long non-coding RNAs (lncRNAs)have been considered as prognostic biomarkers for CRC. But the significance of hypoxic lncRNAs in predicting CRC prognosis remains unclear.Methods:The gene expressed profiles for CRC cases were obtained based on the Cancer Genome Atlas (TCGA) and applied to estimate the hypoxia score using a single-sample gene set enrichment analysis (ssGSEA) algorithm. Overall survival (OS) of the high- and low-hypoxia score group was analyzed by the Kaplan–Meier (KM) plot. To identify differentially expressed lncRNAs (DELs) between two hypoxia score groups, this study carried out differential expression analysis, and then further integrated with the DELs between controls and CRC patients to generate the hypoxia-related lncRNAs for CRC. Besides, prognostic lncRNAs were screened by the univariate Cox regression, which was later utilized for constructing the prognosis nomogram for CRC by adopting the least absolute shrinkage and selection operator (LASSO) algorithm. In addition, both accuracy and specificity of the constructed prognostic signature were detected through the receiver operating characteristic (ROC) analysis. Moreover, our constructed prognosis signature also was validated in the internal testing test. This study operated gene set enrichment analysis (GSEA) for exploring potential biological functions associated with the prognostic signature. Finally, the ceRNA network of the prognostic lncRNAs was constructed.Results:Among 2299 hypoxia-related lncRNAs of CRC in total, LINC00327, LINC00163, LINC00174, SYNPR-AS1, and MIR31HG were identified as prognostic lncRNAs by the univariate Cox regression and adopted for constructing the prognosis signature for CRC. ROC analysis showed the predictive power and accuracy of the prognostic signature. Additionally, the GSEA revealed that ECM-receptor interaction, PI3K-Akt pathway, phagosome, and Hippo pathway were mostly associated with the high-risk group. 352 miRNAs-mRNAs pairs and 177 lncRNAs-miRNAs were predicted.Conclusion:To conclude, we identified 5 hypoxia-related lncRNAs to establish an accurate prognostic signature for CRC, providing important prognostic markers and therapeutic targets.


2021 ◽  
Vol 12 ◽  
Author(s):  
Tao Feng ◽  
Dechao Wei ◽  
Qiankun Li ◽  
Xiaobing Yang ◽  
Yili Han ◽  
...  

Prostate cancer (PCa) is one of the most common malignancies for males, but very little is known about its pathogenesis. This study aimed to identify novel biomarkers associated with PCa prognosis and elucidate the underlying molecular mechanism. First, The Cancer Genome Atlas (TCGA) RNA-sequencing data were utilized to identify differentially expressed genes (DEGs) between tumor and normal samples. The DEGs were then applied to construct a co-expression and mined using structure network analysis. The magenta module that was highly related to the Gleason score (r = 0.46, p = 3e–26) and tumor stage (r = 0.38, p = 2e–17) was screened. Subsequently, all genes of the magenta module underwent function annotation. From the key module, CCNA2, CKAP2L, NCAPG, and NUSAP1 were chosen as the four candidate genes. Finally, internal (TCGA) and external data sets (GSE32571, GSE70770, and GSE141551) were combined to validate and predict the value of real hub genes. The results show that the above genes are up-regulated in PCa samples, and higher expression levels show significant association with higher Gleason scores and tumor T stage. Moreover, receiver operating characteristic curve and survival analysis validate the excellent value of hub genes in PCa progression and prognosis. In addition, the protein levels of these four genes also remain higher in tumor tissues when compared with normal tissues. Gene set enrichment analysis and gene set variation analysis for a single gene reveal the close relation with cell proliferation. Meanwhile, 11 small molecular drugs that have the potential to treat PCa were also screened. In conclusion, our research identified four potential prognostic genes and several candidate molecular drugs for treating PCa.


2021 ◽  
Vol 2021 ◽  
pp. 1-10
Author(s):  
Qianhong Yang ◽  
Xiaolu Bai ◽  
Xiang Li ◽  
Wei Hu

Purpose. Heart failure (HF) is a clinical syndrome caused by ventricular insufficiency. In order to further explore the biomarkers related to HF, we apply the high-throughput database. Materials and Methods. The GSE21610 was applied for the differentially expressed gene (DEG) analysis. The Database for Annotation, Visualization, and Integrated Discovery (DAVID) was performed to assess Gene ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses. The Gene Set Enrichment Analysis (GSEA) was used for gene expression profile GSE21610. The Protein-Protein Interaction (PPI) network and modules were also constructed for research. These hub gene function pathways were estimated in HF progression. Result. We have identified 434 DEGs in total, including 304 downregulated DEGs and 130 upregulated DEGs. GO and KEGG illustrated that DEGs in HF were significantly enriched in G protein-coupled receptor binding, peroxisome, and cAMP signaling pathway. GSEA results showed gene set GSE21610 was gathered in lipid digestion, defense response to fungus, and intestinal lipid absorption. Finally, through analyzing the PPI network, we screened hub genes CDH1, TFRC, CCL2, BUB1B, and CD19 by the Cytoscape software. Conclusion. This study uses a series of bioinformatics technologies to obtain hug genes and key pathways related to HF. These analysis results provide us with new ideas for finding biomarkers and treatment methods for HF.


2020 ◽  
Vol 15 ◽  
Author(s):  
Wei Han ◽  
Dongchen Lu ◽  
Chonggao Wang ◽  
Mengdi Cui ◽  
Kai Lu

Background: In the past decades, the incidence of thyroid cancer (TC) has been gradually increasing, owing to the widespread use of ultrasound scanning devices. However, the key mRNAs, miRNAs, and mRNA-miRNA network in papillary thyroid carcinoma (PTC) has not been fully understood. Material and Methods: In this study, multiple bioinformatics methods were employed, including differential expression analysis, gene set enrichment analysis, and miRNA-mRNA interaction network construction. Results: First, we investigated the key miRNAs that regulated significantly more differentially expressed genes based on GSEA method. Second, we searched for the key miRNAs based on the mRNA-miRNA interaction subnetwork involved in PTC. We identified hsa-mir-1275, hsa-mir-1291, hsa-mir-206 and hsa-mir-375 as the key miRNAs involved in PTC pathogenesis. Conclusion: The integrated analysis of the gene and miRNA expression data not only identified key mRNAs, miRNAs, and mRNA-miRNA network involved in papillary thyroid carcinoma, but also improved our understanding of the pathogenesis of PTC.


Genes ◽  
2021 ◽  
Vol 12 (1) ◽  
pp. 120
Author(s):  
Yiyun Sun ◽  
Dandan Xu ◽  
Chundong Zhang ◽  
Yitao Wang ◽  
Lian Zhang ◽  
...  

We previously demonstrated that proline-rich protein 11 (PRR11) and spindle and kinetochore associated 2 (SKA2) constituted a head-to-head gene pair driven by a prototypical bidirectional promoter. This gene pair synergistically promoted the development of non-small cell lung cancer. However, the signaling pathways leading to the ectopic expression of this gene pair remains obscure. In the present study, we first analyzed the lung squamous cell carcinoma (LSCC) relevant RNA sequencing data from The Cancer Genome Atlas (TCGA) database using the correlation analysis of gene expression and gene set enrichment analysis (GSEA), which revealed that the PRR11-SKA2 correlated gene list highly resembled the Hedgehog (Hh) pathway activation-related gene set. Subsequently, GLI1/2 inhibitor GANT-61 or GLI1/2-siRNA inhibited the Hh pathway of LSCC cells, concomitantly decreasing the expression levels of PRR11 and SKA2. Furthermore, the mRNA expression profile of LSCC cells treated with GANT-61 was detected using RNA sequencing, displaying 397 differentially expressed genes (203 upregulated genes and 194 downregulated genes). Out of them, one gene set, including BIRC5, NCAPG, CCNB2, and BUB1, was involved in cell division and interacted with both PRR11 and SKA2. These genes were verified as the downregulated genes via RT-PCR and their high expression significantly correlated with the shorter overall survival of LSCC patients. Taken together, our results indicate that GLI1/2 mediates the expression of the PRR11-SKA2-centric gene set that serves as an unfavorable prognostic indicator for LSCC patients, potentializing new combinatorial diagnostic and therapeutic strategies in LSCC.


2019 ◽  
Vol 8 (10) ◽  
pp. 1580 ◽  
Author(s):  
Kyoung Min Moon ◽  
Kyueng-Whan Min ◽  
Mi-Hye Kim ◽  
Dong-Hoon Kim ◽  
Byoung Kwan Son ◽  
...  

Ninety percent of patients with scrub typhus (SC) with vasculitis-like syndrome recover after mild symptoms; however, 10% can suffer serious complications, such as acute respiratory failure (ARF) and admission to the intensive care unit (ICU). Predictors for the progression of SC have not yet been established, and conventional scoring systems for ICU patients are insufficient to predict severity. We aimed to identify simple and robust indicators to predict aggressive behaviors of SC. We evaluated 91 patients with SC and 81 non-SC patients who were admitted to the ICU, and 32 cases from the public functional genomics data repository for gene expression analysis. We analyzed the relationships between several predictors and clinicopathological characteristics in patients with SC. We performed gene set enrichment analysis (GSEA) to identify SC-specific gene sets. The acid-base imbalance (ABI), measured 24 h before serious complications, was higher in patients with SC than in non-SC patients. A high ABI was associated with an increased incidence of ARF, leading to mechanical ventilation and worse survival. GSEA revealed that SC correlated to gene sets reflecting inflammation/apoptotic response and airway inflammation. ABI can be used to indicate ARF in patients with SC and assist with early detection.


Sign in / Sign up

Export Citation Format

Share Document