1
|
Larsen FT, Hansen D, Terkelsen MK, Bendixen SM, Avolio F, Wernberg CW, Lauridsen MM, Grønkjaer LL, Jacobsen BG, Klinggaard EG, Mandrup S, Di Caterino T, Siersbæk MS, Indira Chandran V, Graversen JH, Krag A, Grøntved L, Ravnskjaer K. Stellate cell expression of SPARC-related modular calcium-binding protein 2 is associated with human non-alcoholic fatty liver disease severity. JHEP Rep 2023; 5:100615. [PMID: 36687468 PMCID: PMC9850195 DOI: 10.1016/j.jhepr.2022.100615] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 07/05/2022] [Revised: 09/30/2022] [Accepted: 10/15/2022] [Indexed: 11/07/2022] Open
Abstract
Background & Aims Histological assessment of liver biopsies is the gold standard for diagnosis of non-alcoholic steatohepatitis (NASH), the progressive form of non-alcoholic fatty liver disease (NAFLD), despite its well-established limitations. Therefore, non-invasive biomarkers that can offer an integrated view of the liver are needed to improve diagnosis and reduce sampling bias. Hepatic stellate cells (HSCs) are central in the development of hepatic fibrosis, a hallmark of NASH. Secreted HSC-specific proteins may, therefore, reflect disease state in the NASH liver and serve as non-invasive diagnostic biomarkers. Methods We performed RNA-sequencing on liver biopsies from a histologically characterised cohort of obese patients (n = 30, BMI >35 kg/m2) to identify and evaluate HSC-specific genes encoding secreted proteins. Bioinformatics was used to identify potential biomarkers and their expression at single-cell resolution. We validated our findings using single-molecule fluorescence in situ hybridisation (smFISH) and ELISA to detect mRNA in liver tissue and protein levels in plasma, respectively. Results Hepatic expression of SPARC-related modular calcium-binding protein 2 (SMOC2) was increased in NASH compared to no-NAFLD (p.adj <0.001). Single-cell RNA-sequencing data indicated that SMOC2 was primarily expressed by HSCs, which was validated using smFISH. Finally, plasma SMOC2 was elevated in NASH compared to no-NAFLD (p <0.001), with a predictive accuracy of AUROC 0.88. Conclusions Increased SMOC2 in plasma appears to reflect HSC activation, a key cellular event associated with NASH progression, and may serve as a non-invasive biomarker of NASH. Impact and implications Non-alcoholic fatty liver disease (NAFLD) and its progressive form, non-alcoholic steatohepatitis (NASH), are the most common forms of chronic liver diseases. Currently, liver biopsies are the gold standard for diagnosing NAFLD. Blood-based biomarkers to complement liver biopsies for diagnosis of NAFLD are required. We found that activated hepatic stellate cells, a cell type central to NAFLD pathogenesis, upregulate expression of the secreted protein SPARC-related modular calcium-binding protein 2 (SMOC2). SMOC2 was elevated in blood samples from patients with NASH and may hold promise as a blood-based biomarker for the diagnosis of NAFLD.
Collapse
Key Words
- AUROC, area under the receiver operating characteristic curve
- ECM, extracellular matrix
- HSC, hepatic stellate cells
- LSM, liver stiffness measurement
- MCP, matricellular protein
- NAFL, non-alcoholic fatty liver
- NAFLD
- NAFLD, non-alcoholic fatty liver disease
- NAS, NAFLD activity score
- NASH
- PCA, principal component analysis
- SAF, steatosis, activity, and fibrosis
- SE, sensitivity
- SMOC2
- SMOC2, SPARC-related modular calcium-binding protein 2
- SP, specificity
- SPARC, secreted protein acidic and cysteine-rich
- VSMCs, vascular smooth muscle cells
- WGCNA, weighted gene co-expression network analysis
- aHSC, activated HSC
- hepatic stellate cells
- non-invasive biomarker
- qHSC, quiescent HSC
- smFISH, single-molecule fluorescence in situ hybridisation
- transcriptomics
Collapse
Affiliation(s)
- Frederik T. Larsen
- Department of Biochemistry and Molecular Biology, University of Southern
Denmark, Odense, Denmark
- Center for Functional Genomics and Tissue Plasticity (ATLAS), University of
Southern Denmark, Odense, Denmark
| | - Daniel Hansen
- Department of Biochemistry and Molecular Biology, University of Southern
Denmark, Odense, Denmark
- Center for Functional Genomics and Tissue Plasticity (ATLAS), University of
Southern Denmark, Odense, Denmark
| | - Mike K. Terkelsen
- Department of Biochemistry and Molecular Biology, University of Southern
Denmark, Odense, Denmark
- Center for Functional Genomics and Tissue Plasticity (ATLAS), University of
Southern Denmark, Odense, Denmark
| | - Sofie M. Bendixen
- Department of Biochemistry and Molecular Biology, University of Southern
Denmark, Odense, Denmark
- Center for Functional Genomics and Tissue Plasticity (ATLAS), University of
Southern Denmark, Odense, Denmark
| | - Fabio Avolio
- Department of Biochemistry and Molecular Biology, University of Southern
Denmark, Odense, Denmark
- Center for Functional Genomics and Tissue Plasticity (ATLAS), University of
Southern Denmark, Odense, Denmark
| | - Charlotte W. Wernberg
- Center for Functional Genomics and Tissue Plasticity (ATLAS), University of
Southern Denmark, Odense, Denmark
- Department of Gastroenterology and Hepatology, University Hospital of
Southern Denmark, Esbjerg, Denmark
- Center for Liver Research (FLASH), Department of Gastroenterology and
Hepatology, Odense University Hospital, Odense, Denmark
| | - Mette M. Lauridsen
- Center for Functional Genomics and Tissue Plasticity (ATLAS), University of
Southern Denmark, Odense, Denmark
- Department of Gastroenterology and Hepatology, University Hospital of
Southern Denmark, Esbjerg, Denmark
| | - Lea L. Grønkjaer
- Department of Gastroenterology and Hepatology, University Hospital of
Southern Denmark, Esbjerg, Denmark
| | - Birgitte G. Jacobsen
- Department of Gastroenterology and Hepatology, University Hospital of
Southern Denmark, Esbjerg, Denmark
| | - Ellen G. Klinggaard
- Department of Biochemistry and Molecular Biology, University of Southern
Denmark, Odense, Denmark
- Center for Functional Genomics and Tissue Plasticity (ATLAS), University of
Southern Denmark, Odense, Denmark
| | - Susanne Mandrup
- Department of Biochemistry and Molecular Biology, University of Southern
Denmark, Odense, Denmark
- Center for Functional Genomics and Tissue Plasticity (ATLAS), University of
Southern Denmark, Odense, Denmark
| | - Tina Di Caterino
- Department of Pathology, Odense University Hospital, Odense,
Denmark
| | - Majken S. Siersbæk
- Department of Biochemistry and Molecular Biology, University of Southern
Denmark, Odense, Denmark
- Center for Functional Genomics and Tissue Plasticity (ATLAS), University of
Southern Denmark, Odense, Denmark
| | - Vineesh Indira Chandran
- Center for Functional Genomics and Tissue Plasticity (ATLAS), University of
Southern Denmark, Odense, Denmark
- Department of Molecular Medicine, University of Southern Denmark, Odense,
Denmark
| | - Jonas H. Graversen
- Center for Functional Genomics and Tissue Plasticity (ATLAS), University of
Southern Denmark, Odense, Denmark
- Department of Molecular Medicine, University of Southern Denmark, Odense,
Denmark
| | - Aleksander Krag
- Center for Functional Genomics and Tissue Plasticity (ATLAS), University of
Southern Denmark, Odense, Denmark
- Center for Liver Research (FLASH), Department of Gastroenterology and
Hepatology, Odense University Hospital, Odense, Denmark
| | - Lars Grøntved
- Department of Biochemistry and Molecular Biology, University of Southern
Denmark, Odense, Denmark
- Center for Functional Genomics and Tissue Plasticity (ATLAS), University of
Southern Denmark, Odense, Denmark
| | - Kim Ravnskjaer
- Department of Biochemistry and Molecular Biology, University of Southern
Denmark, Odense, Denmark
- Center for Functional Genomics and Tissue Plasticity (ATLAS), University of
Southern Denmark, Odense, Denmark
- Corresponding author. Address: Department of Biochemistry and Molecular
Biology, Campusvej 55, 5230 Odense M, Denmark. Tel.: +45 65508906/+45
93979317.
| |
Collapse
|
2
|
Gao L, Li GS, Li JD, He J, Zhang Y, Zhou HF, Kong JL, Chen G. Identification of the susceptibility genes for COVID-19 in lung adenocarcinoma with global data and biological computation methods. Comput Struct Biotechnol J 2021; 19:6229-6239. [PMID: 34840672 PMCID: PMC8605816 DOI: 10.1016/j.csbj.2021.11.026] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/06/2021] [Revised: 11/07/2021] [Accepted: 11/16/2021] [Indexed: 12/19/2022] Open
Abstract
Introduction The risk of infection with COVID-19 is high in lung adenocarcinoma (LUAD) patients, and there is a dearth of studies on the molecular mechanism underlying the high susceptibility of LUAD patients to COVID-19 from the perspective of the global differential expression landscape. Objectives To fill the research void on the molecular mechanism underlying the high susceptibility of LUAD patients to COVID-19 from the perspective of the global differential expression landscape. Methods Herein, we identified genes, specifically the differentially expressed genes (DEGs), correlated with the susceptibility of LUAD patients to COVID-19. These were obtained by calculating standard mean deviation (SMD) values for 49 SARS-CoV-2-infected LUAD samples and 24 non-affected LUAD samples, as well as 3931 LUAD samples and 3027 non-cancer lung samples from 40 pooled RNA-seq and microarray datasets. Hub susceptibility genes significantly related to COVID-19 were further selected by weighted gene co-expression network analysis. Then, the hub genes were further analyzed via an examination of their clinical significance in multiple datasets, a correlation analysis of the immune cell infiltration level, and their interactions with the interactome sets of the A549 cell line. Results A total of 257 susceptibility genes were identified, and these genes were associated with RNA splicing, mitochondrial functions, and proteasomes. Ten genes, MEA1, MRPL24, PPIH, EBNA1BP2, MRTO4, RABEPK, TRMT112, PFDN2, PFDN6, and NDUFS3, were confirmed to be the hub susceptibility genes for COVID-19 in LUAD patients, and the hub susceptibility genes were significantly correlated with the infiltration of multiple immune cells. Conclusion In conclusion, the susceptibility genes for COVID-19 in LUAD patients discovered in this study may increase our understanding of the high risk of COVID-19 in LUAD patients.
Collapse
Key Words
- CI, confidence interval
- COVID-19
- COVID-19, coronavirus disease 2019
- DEG
- DEG, differentially expressed genes
- FC, fold change
- FPKM, fragments per kilobase per million
- GTEx, Genotype-tissue Expression
- HPA, human protein atlas
- IHC, immunohistochemistry
- Immune infiltration
- LUAD
- LUAD, lung adenocarcinoma
- PPI, protein-to-protein interaction
- SARS-CoV-2, severe acute respiratory syndrome coronavirus 2
- SMD, standard mean difference
- SROC, summarized receiver’s operating characteristics
- Susceptibility
- TF, transcription factor
- TPM, transcripts per million reads
- WGCNA
- WGCNA, weighted gene co-expression network analysis
Collapse
Affiliation(s)
- Li Gao
- Department of Pathology, The First Affiliated Hospital of Guangxi Medical University, No. 6 Shuangyong Rd, Nanning, Guangxi Zhuang Autonomous Region 530021, PR China
| | - Guo-Sheng Li
- Department of Pathology, The First Affiliated Hospital of Guangxi Medical University, No. 6 Shuangyong Rd, Nanning, Guangxi Zhuang Autonomous Region 530021, PR China
| | - Jian-Di Li
- Department of Pathology, The First Affiliated Hospital of Guangxi Medical University, No. 6 Shuangyong Rd, Nanning, Guangxi Zhuang Autonomous Region 530021, PR China
| | - Juan He
- Department of Pathology, The First Affiliated Hospital of Guangxi Medical University, No. 6 Shuangyong Rd, Nanning, Guangxi Zhuang Autonomous Region 530021, PR China
| | - Yu Zhang
- Department of Pathology, Shandong Provincial Hospital Affiliated to Shandong First Medical University, No. 324. Jingwu Rd, Jinan, Shandong 250021, PR China
| | - Hua-Fu Zhou
- Department of Cardio-Thoracic Surgery, The First Affiliated Hospital of Guangxi Medical University, No. 6 Shuangyong Rd, Nanning, Guangxi Zhuang Autonomous Region 530021, PR China
| | - Jin-Liang Kong
- Ward of Pulmonary and Critical Care Medicine, Department of Respiratory Medicine, The First Affiliated Hospital of Guangxi Medical University, No. 6 Shuangyong Rd, Nanning, Guangxi Zhuang Autonomous Region 530021, PR China
| | - Gang Chen
- Department of Pathology, The First Affiliated Hospital of Guangxi Medical University, No. 6 Shuangyong Rd, Nanning, Guangxi Zhuang Autonomous Region 530021, PR China
| |
Collapse
|
3
|
Zhang H, Luo YB, Wu W, Zhang L, Wang Z, Dai Z, Feng S, Cao H, Cheng Q, Liu Z. The molecular feature of macrophages in tumor immune microenvironment of glioma patients. Comput Struct Biotechnol J 2021; 19:4603-18. [PMID: 34471502 DOI: 10.1016/j.csbj.2021.08.019] [Citation(s) in RCA: 66] [Impact Index Per Article: 22.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2021] [Revised: 08/11/2021] [Accepted: 08/12/2021] [Indexed: 12/12/2022] Open
Abstract
Background Gliomas are one of the most common types of primary tumors in central nervous system. Previous studies have found that macrophages actively participate in tumor growth. Methods Weighted gene co-expression network analysis was used to identify meaningful macrophage-related gene genes for clustering. Pamr, SVM, and neural network were applied for validating clustering results. Somatic mutation and methylation were used for defining the features of identified clusters. Differentially expressed genes (DEGs) between the stratified groups after performing elastic regression and principal component analyses were used for the construction of MScores. The expression of macrophage-specific genes were evaluated in tumor microenvironment based on single cell sequencing analysis. A total of 2365 samples from 15 glioma datasets and 5842 pan-cancer samples were used for external validation of MScore. Results Macrophages were identified to be negatively associated with the survival of glioma patients. Twenty-six macrophage-specific DEGs obtained by elastic regression and PCA were highly expressed in macrophages at single-cell level. The prognostic value of MScores in glioma was validated by the active proinflammatory and metabolic profile of infiltrating microenvironment and response to immunotherapies of samples with this signature. MScores managed to stratify patient survival probabilities in 15 external glioma datasets and pan-cancer datasets, which predicted worse survival outcome. Sequencing data and immunohistochemistry of Xiangya glioma cohort confirmed the prognostic value of MScores. A prognostic model based on MScores demonstrated high accuracy rate. Conclusion Our findings strongly support a modulatory role of macrophages, especially M2 macrophages in glioma progression and warrants further experimental studies.
Collapse
Key Words
- ACC, Adrenocortical carcinoma
- BBB, brain blood barrier
- BLCA, Bladder Urothelial Carcinoma
- BRCA, Breast invasive carcinoma
- CDF, cumulative distribution function
- CESC, Cervical squamous cell carcinoma and endocervical adenocarcinoma
- CGGA, Chinese Glioma Genome Atlas
- CHOL, Cholangiocarcinoma
- CNA, copy number alternations
- CNV, copy number variation
- COAD, Colon adenocarcinoma
- CSF-1, colony-stimulating factor-1
- DLBC, Lymphoid Neoplasm Diffuse Large B-cell Lymphoma
- DMP, differentially methylated position
- ESCA, Esophageal carcinoma
- GBM, glioblastoma
- GEO, Gene Expression Omnibus
- GO, gene ontology
- GSEA, gene set enrichment analysis
- GSVA, gene set variation analysis
- Glioma microenvironment
- HNSC, Head and Neck squamous cell carcinoma
- IGR, intergenic region
- IHC, immunohistochemistry
- IL, interleukin
- Immunotherapy
- KEGG, Kyoto Encyclopaedia of Genes and Genomes
- KICH, Kidney Chromophobe
- KIRC, Kidney renal clear cell carcinoma
- KIRP, Kidney renal papillary cell carcinoma
- LGG, low grade glioma
- LIHC, Liver hepatocellular carcinoma
- LUAD, Lung adenocarcinoma
- LUSC, Lung squamous cell carcinoma
- MMP-2, matrix metalloproteinase-2
- MT1, MMP membrane type 1 matrix metalloprotease
- Machine learning
- Macrophage
- OV, Ovarian serous cystadenocarcinoma
- PAAD, Pancreatic adenocarcinoma
- PAM, partition around medoids
- PCA, principal component analysis
- PCPG, Pheochromocytoma and Paraganglioma
- PRAD, Prostate adenocarcinoma
- Prognostic model
- READ, Rectum adenocarcinoma
- SARC, Sarcoma
- SKCM, Skin Cutaneous Melanoma
- SNP, single-nucleotide polymorphism
- SNV, single-nucleotide variant
- STAD, Stomach adenocarcinoma
- SVM, Support Vector Machines
- TAM, tumor associated macrophage
- TCGA, The Cancer Genome Atlas
- TGF-β, tumor growth factor-β
- THCA, Thyroid carcinoma
- THYM, Thymoma
- TIMP-2, tissue inhibitor of metalloproteinase-2
- TLR2, toll-like receptor 2
- TME, tumor microenvironment
- TNFα, tumor necrosis factor α
- TSS, transcription start site
- UCEC, Uterine Corpus Endometrial Carcinoma
- UCS, Uterine Carcinosarcoma
- WGCNA, weighted gene co-expression network analysis
- pamr, prediction analysis for microarrays
Collapse
|
4
|
Liu Q, Nie R, Li M, Li L, Zhou H, Lu H, Wang X. Identification of subtypes correlated with tumor immunity and immunotherapy in cutaneous melanoma. Comput Struct Biotechnol J 2021; 19:4472-4485. [PMID: 34471493 PMCID: PMC8379294 DOI: 10.1016/j.csbj.2021.08.005] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2021] [Revised: 08/04/2021] [Accepted: 08/04/2021] [Indexed: 01/15/2023] Open
Abstract
Because immune checkpoint inhibitors (ICIs) are effective for a subset of melanoma patients, identification of melanoma subtypes responsive to ICIs is crucial. We performed clustering analyses to identify immune subtypes of melanoma based on the enrichment levels of 28 immune cells using transcriptome datasets for six melanoma cohorts, including four cohorts not treated with ICIs and two cohorts treated with ICIs. We identified three immune subtypes (Im-H, Im-M, and Im-L), reproducible in these cohorts. Im-H displayed strong immune signatures, low stemness and proliferation potential, genomic stability, high immunotherapy response rate, and favorable prognosis. Im-L showed weak immune signatures, high stemness and proliferation potential, genomic instability, low immunotherapy response rate, and unfavorable prognosis. The pathways highly enriched in Im-H included immune, MAPK, apoptosis, calcium, VEGF, cell adhesion molecules, focal adhesion, gap junction, and PPAR. The pathways highly enriched in Im-L included Hippo, cell cycle, and ErbB. Copy number alterations correlated inversely with immune signatures in melanoma, while tumor mutation burden showed no significant correlation. The molecular features correlated with favorable immunotherapy response included immune-promoting signatures and pathways of PPAR, MAPK, VEGF, calcium, and glycolysis/gluconeogenesis. Our data recapture the immunological heterogeneity in melanoma and provide clinical implications for the immunotherapy of melanoma.
Collapse
Key Words
- Clustering analysis
- DMFS, distant-metastasis free survival
- DSS, disease-specific survival
- EMT, epithelial-mesenchymal transition
- FDR, false discovery rate
- GO, gene ontology
- GSEA, gene-set enrichment analysis
- HLA, human leukocyte antigen
- HRD, homologous recombination deficiency
- ICIs, immune checkpoint inhibitors
- Immune subtypes
- Immunotherapy
- MDSC, myeloid-derived suppressor cell
- Melanoma
- NK, natural killer
- OS, overall survival
- SCNAs, somatic copy number alterations
- TCGA, The Cancer Genome Atlas
- TIME, tumor immune microenvironment
- TMB, tumor mutation burden
- TME, tumor microenvironment
- Tumor immune microenvironment
- WGCNA, weighted gene co-expression network analysis
- ssGSEA, single-sample gene-set enrichment analysis
Collapse
Affiliation(s)
- Qian Liu
- Biomedical Informatics Research Lab, School of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing 211198, China
- Cancer Genomics Research Center, School of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing 211198, China
- Big Data Research Institute, China Pharmaceutical University, Nanjing 211198, China
| | - Rongfang Nie
- Biomedical Informatics Research Lab, School of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing 211198, China
- Cancer Genomics Research Center, School of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing 211198, China
- Big Data Research Institute, China Pharmaceutical University, Nanjing 211198, China
| | - Mengyuan Li
- Biomedical Informatics Research Lab, School of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing 211198, China
- Cancer Genomics Research Center, School of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing 211198, China
- Big Data Research Institute, China Pharmaceutical University, Nanjing 211198, China
| | - Lin Li
- Biomedical Informatics Research Lab, School of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing 211198, China
- Cancer Genomics Research Center, School of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing 211198, China
- Big Data Research Institute, China Pharmaceutical University, Nanjing 211198, China
| | - Haiying Zhou
- Department of Orthopedics, The First Affiliated Hospital, College of Medicine, Zhejiang University, Hangzhou 310003, China
| | - Hui Lu
- Department of Orthopedics, The First Affiliated Hospital, College of Medicine, Zhejiang University, Hangzhou 310003, China
| | - Xiaosheng Wang
- Biomedical Informatics Research Lab, School of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing 211198, China
- Cancer Genomics Research Center, School of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing 211198, China
- Big Data Research Institute, China Pharmaceutical University, Nanjing 211198, China
| |
Collapse
|
5
|
Zhang Z, Li L, Li M, Wang X. The SARS-CoV-2 host cell receptor ACE2 correlates positively with immunotherapy response and is a potential protective factor for cancer progression. Comput Struct Biotechnol J. 2020;18:2438-2444. [PMID: 32905022 PMCID: PMC7462778 DOI: 10.1016/j.csbj.2020.08.024] [Citation(s) in RCA: 30] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2020] [Revised: 08/23/2020] [Accepted: 08/26/2020] [Indexed: 02/06/2023] Open
Abstract
The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has infected more than 29 million people and has caused more than 900,000 deaths worldwide as of September 14, 2020. The SARS-CoV-2 human cell receptor ACE2 has recently received extensive attention for its role in SARS-CoV-2 infection. Many studies have also explored the association between ACE2 and cancer. However, a systemic investigation into associations between ACE2 and oncogenic pathways, tumor progression, and clinical outcomes in pan-cancer remains lacking. Using cancer genomics datasets from the Cancer Genome Atlas (TCGA) program, we performed computational analyses of associations between ACE2 expression and antitumor immunity, immunotherapy response, oncogenic pathways, tumor progression phenotypes, and clinical outcomes in 13 cancer cohorts. We found that ACE2 upregulation was associated with increased antitumor immune signatures and PD-L1 expression, and favorable anti-PD-1/PD-L1/CTLA-4 immunotherapy response. ACE2 expression levels inversely correlated with the activity of cell cycle, mismatch repair, TGF-β, Wnt, VEGF, and Notch signaling pathways. Moreover, ACE2 expression levels had significant inverse correlations with tumor proliferation, stemness, and epithelial-mesenchymal transition. ACE2 upregulation was associated with favorable survival in pan-cancer and in multiple individual cancer types. These results suggest that ACE2 is a potential protective factor for cancer progression. Our data may provide potential clinical implications for treating cancer patients infected with SARS-CoV-2.
Collapse
Key Words
- ACE2 expression
- ACE2, angiotensin-converting enzyme 2
- CESC, cervical squamous-cell carcinoma
- COAD, colon adenocarcinoma
- DFI, disease-free interval
- DSS, disease-specific survival
- EMT, epithelial-mesenchymal transition
- ESCA, esophageal carcinoma
- FDR, false discovery rate
- GO, gene ontology
- GSEA, gene set enrichment analysis
- HNSC, head and neck squamous cell carcinoma
- KIRC, kidney renal clear cell carcinoma
- KIRP, kidney renal papillary cell carcinoma
- LUAD, lung adenocarcinoma
- LUSC, lung squamous cell carcinoma
- OS, overall survival
- OV, ovarian carcinoma
- PAAD, pancreatic adenocarcinoma
- PFI, progression-free interval
- Pan-cancer
- SARS-CoV-2, severe acute respiratory syndrome coronavirus 2
- SKCM, skin cutaneous melanoma
- Survival prognosis
- TCGA, The Cancer Genome Atlas
- TF, transcription factor
- THYM, thymoma
- Tumor immunity and immunotherapy
- Tumor progression
- UCEC, uterine corpus endometrial carcinoma
- WGCNA, weighted gene co-expression network analysis
Collapse
|
6
|
Yuan L, Chen L, Qian K, Qian G, Wu CL, Wang X, Xiao Y. Co-expression network analysis identified six hub genes in association with progression and prognosis in human clear cell renal cell carcinoma (ccRCC). Genom Data 2017; 14:132-140. [PMID: 29159069 PMCID: PMC5683669 DOI: 10.1016/j.gdata.2017.10.006] [Citation(s) in RCA: 61] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/10/2017] [Revised: 10/12/2017] [Accepted: 10/25/2017] [Indexed: 12/21/2022]
Abstract
Human clear cell renal cell carcinoma (ccRCC) is one of the most common types of malignant adult kidney tumors. We constructed a weighted gene co-expression network to identify gene modules associated with clinical features of ccRCC (n = 97). Six hub genes (CCNB2, CDC20, CEP55, KIF20A, TOP2A and UBE2C) were identified in both co-expression and protein-protein interaction (PPI) networks, which were highly correlated with pathologic stage. The significance of expression of the hub genes in ccRCC was ranked top 4 among all cancers and correlated with poor prognosis. Functional analysis revealed that the hub genes were significantly enriched in cell cycle regulation and cell division. Gene set enrichment analysis suggested that the samples with highly expressed hub gene were correlated with cell cycle and p53 signaling pathway. Taken together, six hub genes were identified to be associated with progression and prognosis of ccRCC, and they might lead to poor prognosis by regulating p53 signaling pathway.
Collapse
Key Words
- Clear cell renal cell carcinoma (ccRCC)
- Co-expression network analysis
- DAVID, Database for Annotation, Visualization and Integrated Discovery
- DEG, differentially expressed gene
- DEGs, differentially expressed genes
- GS, gene significance
- GSEA, enrichment analysis and gene set enrichment
- HPA, human protein atlas
- Hub genes
- MEs, module eigengenes
- MS, module significance
- PPI, protein-protein interaction
- Prognosis
- Progression
- SAM, significance analysis of microarrays
- STRING, search tool for the retrieval of interacting genes
- TCGA, the cancer genome atlas
- TOM, topological overlap matrix
- WGCNA, weighted gene co-expression network analysis
- ccRCC, clear cell renal cell carcinoma
Collapse
Affiliation(s)
- Lushun Yuan
- Department of Urology, Zhongnan Hospital of Wuhan University, Wuhan, China
| | - Liang Chen
- Department of Urology, Zhongnan Hospital of Wuhan University, Wuhan, China
| | - Kaiyu Qian
- Department of Urology, Zhongnan Hospital of Wuhan University, Wuhan, China
- Department of Urology, The Fifth Hospital of Wuhan, Wuhan, China
| | - Guofeng Qian
- Department of Endocrinology, The First Affiliated Hospital of Zhejiang University, Hangzhou, China
| | - Chin-Lee Wu
- Department of Urology, Massachusetts General Hospital, Harvard Medical School, Boston, MA, USA
| | - Xinghuan Wang
- Department of Urology, Zhongnan Hospital of Wuhan University, Wuhan, China
- Corresponding author.
| | - Yu Xiao
- Department of Urology, Zhongnan Hospital of Wuhan University, Wuhan, China
- Laboratory of Precision Medicine, Zhongnan Hospital of Wuhan University, Wuhan, China
- Department of Biological Repositories, Zhongnan Hospital of Wuhan University, Wuhan, China
- Correspondence to: Y. Xiao, Department of Biological Repositories, Zhongnan Hospital of Wuhan University, Wuhan, China.Department of Biological RepositoriesZhongnan Hospital of Wuhan UniversityWuhanChina
| |
Collapse
|
7
|
Liu R, Guo CX, Zhou HH. Network-based approach to identify prognostic biomarkers for estrogen receptor-positive breast cancer treatment with tamoxifen. Cancer Biol Ther 2015; 16:317-24. [PMID: 25756514 DOI: 10.1080/15384047.2014.1002360] [Citation(s) in RCA: 45] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Abstract
This study aims to identify effective gene networks and prognostic biomarkers associated with estrogen receptor positive (ER+) breast cancer using human mRNA studies. Weighted gene coexpression network analysis was performed with a complex ER+ breast cancer transcriptome to investigate the function of networks and key genes in the prognosis of breast cancer. We found a significant correlation of an expression module with distant metastasis-free survival (HR = 2.25; 95% CI .21.03-4.88 in discovery set; HR = 1.78; 95% CI = 1.07-2.93 in validation set). This module contained genes enriched in the biological process of the M phase. From this module, we further identified and validated 5 hub genes (CDK1, DLGAP5, MELK, NUSAP1, and RRM2), the expression levels of which were strongly associated with poor survival. Highly expressed MELK indicated poor survival in luminal A and luminal B breast cancer molecular subtypes. This gene was also found to be associated with tamoxifen resistance. Results indicated that a network-based approach may facilitate the discovery of biomarkers for the prognosis of ER+ breast cancer and may also be used as a basis for establishing personalized therapies. Nevertheless, before the application of this approach in clinical settings, in vivo and in vitro experiments and multi-center randomized controlled clinical trials are still needed.
Collapse
Key Words
- CI, confidence interval
- ER+, estrogen receptor positive
- GS, gene significance
- HER2, human epidermal growth factor 2
- ME, module eigengene
- MS, module significance
- PCC, Pearson's correlation coefficient
- PR, progesterone receptor
- TOM, topologic overlap measure
- WGCNA, weighted gene co-expression network analysis
- biomarker
- breast cancer
- gene expression profiling
- k.in, intramodular connectivity
- k.total, Network connectivity
- systems biology
- tamoxifen resistance
Collapse
Affiliation(s)
- Rong Liu
- a Department of Clinical Pharmacology; Xiangya Hospital; Central South University ; Changsha , China
| | | | | |
Collapse
|