127
Views
2
CrossRef citations to date
0
Altmetric
Original Research

A Pan-Cancer Analysis of the Prognostic Value and Expression of Adenylate Cyclase 7 (ADCY7) in Human Tumors

ORCID Icon, , , , , , , , & show all
Pages 5415-5429 | Published online: 11 Sep 2021

Abstract

Background

The role of adenylate cyclase 7 (ADCY7) in cancer is still unclear. This study analyzed the interrelationship between the expression and immune function of ADCY7.

Methods

ADCY7 expression in multiple human cancers was analyzed using the databases of Genotype-Tissue Expression Project (GTEx), Cancer Cell Line Encyclopedia (CCLE), and The Cancer Genome Atlas (TCGA). Correlations among ADCY7 gene expression, mismatch repair (MMR) gene expression, and DNA methyltransferase (DNMT) expression were assessed using Spearman correlation analysis. Univariate survival analysis and Kaplan–Meier (KM) curve were used to examine the effect of ADCY7 expression on prognosis. The Tumor Immune Estimation Resource (TIMER) database was used to evaluate the relationship between ADCY7 gene expression and tumor immune invasion or immune checkpoint gene (ICG) expression.

Results

ADCY7 was abnormally expressed in multiple human cancers and was correlated with MMR genes and DNMT expression. Univariate survival analysis and KM curve showed that ADCY7 expression influences the overall survival (OS) of six types of cancer, disease-specific survival (DSS) of eight, and progression-free interval (PFI) of three. The high expression of ADCY7 in OS, DSS, and PFI was strongly associated with poor outcomes in patients with breast cancer and lung squamous cell carcinoma. ADCY7 expression was strongly associated with immune cell infiltration and ICG expression.

Conclusion

The results of this study indicated that ADCY7 may be a prognostic biomarker of tumorigenesis. The study may also provide a new perspective on the role of ADCY7 in human cancers.

Introduction

Adenylate cyclase 7 (ADCY7) encodes a membrane protein that is a member of the adenylate cyclase family. It promotes transduction of extracellular signals into intracellular reactions.Citation1 Studies have reported that ADCY7 can catalyze the generation of cyclic adenosine monophosphate (cAMP) from adenosine triphosphate (ATP).Citation2,Citation3 ADCY7 plays a critical role in nervous system diseases, inflammatory responses and immune responses.Citation4Citation6 Studies have shown that the adenylate cyclase gene family plays a key role in tumor development. For example, ADCY2 is associated with the prognosis of malignant melanoma and colorectal cancer,Citation7 while ADCY4 and ADCY7 are associated with poorer prognosis of acute myeloid leukemia.Citation8 Among the members of the adenylate cyclase family, the function of ADCY1 in tumors have been widely investigated. ADCY1 was highly expressed and demonstrated a poorer prognosis in most cancers.Citation9 However, no pan-cancer analysis was performed on the association of ADCY7 with multiple tumor types.

Cancer is one of the major public problems threatening global health. This threat may be attributed to the increasing incidence of malignant tumors in patients with cancer. This increase is influenced by factors such as unhealthy lifestyle, genetic makeup, and living environment of the patients. Even in countries with very advanced medical technology, the treatment success rate for patients with cancer is very low,Citation10 which is attributed to genetic alteration, tumor immunity, tumor recurrence, and prognosis. Thus, investigating their mechanisms was thought to improve cancer diagnosis and treatment. Given the commonality and uniqueness of tumor formation and development, analyzing the pan-cancer expression of any genes that may influence cancer occurrence and development is necessary, to evaluate the association of such genes with clinical prognosis and molecular mechanisms.

In this study, The Cancer Genome Atlas (TCGA) and other tumor research databases, such as the Genotype-Tissue Expression Project (GTEx), Cancer Cell Line Encyclopedia (CCLE), cBioPortal, and Tumor IMmune Estimation Resource (TIMER), were used for the first time to access a pan-cancer analysis of ADCY7. Various factors were considered, including genetic alterations, DNA methylation, gene expression, immune infiltration, survival status, and related signaling pathways, in exploring the specific mechanisms of ADCY7 in the tumorigenesis and progression of different cancers. The results of this study revealed that ADCY7 might be significant in tumorigenesis and immune infiltration and may serve as a prognostic biomarker in many cancers.

Methods

Analysis of ADCY7 Expression in Human Cancer

The GTEx database contains the RNA-sequencing (RNA-seq) data from normal tissues.Citation11 GTEx shares its data resource on the differences in gene regulation between tissues and between individuals, which primarily involve tissue-specific gene expression and identification of the associations of gene expression levels in different tissues. In this study, the GTEx database was used to extract the ADCY7 expression data of 31 normal tissues. In addition, the CCLE database contains information on gene expression, methylation, gene mutation, and copy number variation in cancer cells. This database was used to obtain the ADCY7 expression data of 21 tumor cell lines. TCGA database contains bioinformatics data, such as mRNA/microRNA expression profiles, copy number variation, and mutations of related cancer genes. In this study, GTEx and TCGA data were combined to perform the differential analysis of ADCY7 expression between cancer tissues and normal tissues. The expression data were normalized through log2 conversion.

Genetic Alteration Analysis of ADCY7 in Human Cancer

The cBioPortal is a cancer research database that can help in cancer tissue and cytology analysis regrading molecular data recognition and understanding of genetics, epigenetics, gene expression, and proteomics.Citation12 This database was used to analyze genetic alteration characteristics, including missense mutation with uncharted significance, copy number amplification, deep deletion, and mRNA upregulation of ADCY7.

Relationship Between ADCY7 Expression and the Expression of Mismatch Repair (MMR) Gene and DNA Methyltransferase (DNMT) Expressions

The MMR system is a safety guarantee system used in repairing DNA base pair mismatch in human cells.Citation13 It helps maintain the integrity and stability of the genetic material and prevents genetic mutation. TCGA database was used to extract the expressions of MMR genes, including MLH1, MSH2, MSH6, PMS2, and EPCAM, to analyze the relationship between MMR gene expression and ADCY7 expression using Spearman correlation analysis. DNA methylation, a relatively stable modification state, is an important epigenetic mechanism that can be inherited by the new offspring DNA during replication through DNMT.Citation14 DNA methylation is a common and normal modification of eukaryotic cells and a major epigenetic form of gene expression regulation in mammals.Citation15,Citation16 The expressions of four DNMTs, namely, DNMT1, DNMT2, DNMT3A, and DNMT3B, were downloaded from TCGA database to perform the relationship analysis between DNMT expression and ADCY7 expression by Spearman correlation analysis.

Relationship Between ADCY7 Expression and Tumor Mutational Burden (TMB) and Microsatellite Instability (MSI) in Cancers

TMB is associated with overall survival (OS) after immunotherapy for various cancers, suggesting that TMB can be used as a biomarker to predict the efficacy of immune checkpoint inhibitor therapy.Citation17,Citation18 MSI is defined as the occurrence of new microsatellite alleles at a microsatellite locus in tumors as a result of the insertion or deletion of duplicate units, compared with normal tissues.Citation19 MSI associated with DNA MMR defects is an important clinical tumor marker. In this study, the relationship between ADCY7 expression and TMB and MSI was investigated using Spearman correlation analysis.

Immune Correlation Analysis

The TIMER database is a free and publicly accessible database where immune cell infiltration in tumor tissues is detected using RNA-seq data.Citation20 Tumor-infiltrating immune cells are critical to the spread of cancer. In rumors macrophages and monocytes may serve as specific biomarkers for cancer, helping physicians develop novel cancer diagnostic and therapeutic strategies. In this study, the TIMER database was used to investigate the relationship between ADCY7 expression and immune infiltrates based on all TCGA cancers. The relationship between ADCY7 expression and six types of immune cells was evaluated using Spearman correlation analysis.

Survival Prognosis Analysis

TCGA database also contains clinical data on a variety of human cancers, including information on the patients’ prognosis. This information, including OS time, disease-specific survival (DSS) time, and progression-free interval (PFI) time, was downloaded from TCGA database. The relationship between prognosis and ADCY7 expression was investigated in 33 types of cancer. The forest plots and Kaplan–Meier (KM) curves were used to present the results.

ADCY7-Correlated Gene Functional Enrichment Analyses

The STRING database is another free and publicly accessible database that can be used to explore protein-protein interactions.Citation21 First, the top 50 genes closely correlated with ADCY7 were screened by the STRING database. The thresholds set were set as follows: the organism was Homo sapiens, the minimum interaction score was 0.9, the meaning of network edges was evidence, the maximum number of interactors to show was less than 50. Moreover, ADCY7 and its correlated genes were enriched and analyzed using the “clusterProfiler” package in R software for Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, respectively.Citation22,Citation23

Statistical Analysis

Kruskal–Wallis test was used to analyze ADCY7 expression levels in different normal tissues and different tumor cell lines. Differences in ADCY7 expression levels in tumor tissues and normal tissues were evaluated by t-test. Univariate survival analysis was used to analyze the correlation of ADCY7 expression and patients’ survival. KM methods were employed to compare survival by different levels of BRAP expression. Spearman correlation analysis was performed between BRAP expression and MMR gene mutation level, methyltransferases level, and immune checkpoint marker level. Correlations were considered significant and positive at P <0.05 and r >0.20. P <0.05 was considered significant for all statistical analyses.

Results

Pan-Cancer Expression Analysis of ADCY7

First, ADCY7 expression in 31 human tissues was analyzed using the GTEx database. shows that ADCY7 had high expressions in the bone marrow, lung, ovary, and spleen, but low expressions in the muscle, pancreas, and kidney tissues. Then, ADCY7 expressions in 21 human tumor cells were analyzed using the CCLE database. As shown in , ADCY7 was expressed in all 21 tumor cells. Moreover, the differential expression analysis of ADCY7 between normal tissues and 20 types of tumors was performed using TCGA database. The mRNA expression of ADCY7 was higher in bladder cancer (BLCA), cholangiocarcinoma (CHOL), esophageal carcinoma (ESCA), head and neck squamous cell carcinoma (HNSC), kidney renal clear cell carcinoma (KIRC), kidney renal papillary cell carcinoma (KIRP), liver hepatocellular carcinoma (LIHC), stomach adenocarcinoma (STAD), and uterine corpus endometrial carcinoma (UCEC) tissues than in normal tissues but lower in lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC) tissues than in normal tissues (). Since few normal tissues were included in TCGA database, TCGA and GTEx databases were merged for differential expression analysis. As shown in , ADCY7 was differentially expressed in 22 tumor tissues, including up-expression in CHOL, glioblastoma multiforme (GBM), HNSC, KIRC, KIRP, brain lower-grade glioma (LGG), pancreatic adenocarcinoma (PAAD), STAD, testicular germ cell tumors (TGCT), and UCEC tumor tissues, and down-expression in adenoid cystic carcinoma (ACC), breast cancer (BRCA), colon adenocarcinoma (COAD), LUAD, lung squamous cell carcinoma (LUSC), LIHC, thyroid carcinoma (THCA), ovarian serous cystadenocarcinoma (OV), uterine carcinosarcoma (UCS), prostate adenocarcinoma (PRAD), and skin cutaneous melanoma (SKCM) tumor tissues compared with normal tissues. ADCY7 is abnormally expressed in numerous tumors.

Figure 1 The expression of ADCY7 in cancers. (A) The expression of ADCY7 in 31 tissues (GTEx). (B) The expression of ADCY7 in twenty-one kinds of tumor cells (CCLE). (C) Differential expression of ADCY7 in cancers and normal tissues from TCGA database. (D) ADCY7 is differentially expressed in 22 tumor tissues from the GTEx database and TCGA database. *P-value < 0.05, **P-value < 0.01, ***P-value < 0.001, - No significant.

Abbreviations: ADCY7, Adenylate cyclase 7; GTEx, Genotype-Tissue Expression Project; CCLE, Cancer Cell Line Encyclopedia; TCGA, The Cancer Genome Atlas.
Figure 1 The expression of ADCY7 in cancers. (A) The expression of ADCY7 in 31 tissues (GTEx). (B) The expression of ADCY7 in twenty-one kinds of tumor cells (CCLE). (C) Differential expression of ADCY7 in cancers and normal tissues from TCGA database. (D) ADCY7 is differentially expressed in 22 tumor tissues from the GTEx database and TCGA database. *P-value < 0.05, **P-value < 0.01, ***P-value < 0.001, - No significant.

Genetic Alteration Analysis

A genomic mutation is the main driving factor in the occurrence and development of tumors. The cBioportal database was used to analyze the characteristics of ADCY7 genetic alterations (). ADCY7 mutation is one of the most important factors because of its expression alteration in the following tumor tissues: SKCM, UCEC, stomach adenocarcinoma, LUAD, colorectal adenocarcinoma, bladder urothelial carcinoma, LUSC, esophageal adenocarcinoma, glioblastoma multiforme, HNSCC, adrenocortical carcinoma, CSCC, thymoma, KIRC, KIRP, brain lower grade glioma, and thyroid carcinoma. The deep deletion of ADCY7 is one of the most important factors due to its expression alteration in the following tumor tissues: uterine carcinosarcoma, diffuse large B-cell lymphoma, sarcoma, prostate adenocarcinoma, breast-invasive carcinoma, mesothelioma, and liver hepatocellular carcinoma.

Figure 2 The kinds of genetic alterations and frequency of ADCY7 in different tumors (cBioPortal). The cBioPortal database was applied to study the ADCY7 mutation in cancers. The results are displayed as a histogram of the alteration frequencies of ADCY7 across cancer studies. Color images are available online.

Abbreviation: ADCY7, Adenylate cyclase 7.
Figure 2 The kinds of genetic alterations and frequency of ADCY7 in different tumors (cBioPortal). The cBioPortal database was applied to study the ADCY7 mutation in cancers. The results are displayed as a histogram of the alteration frequencies of ADCY7 across cancer studies. Color images are available online.

Relationship Between ADCY7 Expression, MMR Genes, and DNA Methyltransferase Expression

To study the role of ADCY7 in the genesis and development of tumors, the relationship between ADCY7 expression, MMR genes, and DNMTs expression were analyzed. The expression of ADCY7 was associated with the expression of five MMR genes in most tumors, except for ACC and GBM (). The connection analysis between ADCY7 expression and 4 DNMTs expression was performed. shows a strong relationship between ADCY7 expression and expressions of four DNMTs in cancers, except for ESCA, GBM, and ACC. The result suggested that ADCY7 may regulate the DNA damage process and repair of the methylation status of related genes so that they can participate in the development of human cancers.

Figure 3 The expression of ADCY7 is related to the expression of five MMR genes and DNMTs expression in cancers. (A) The relationship between the expression of ADCY7 and expression of MMR genes in cancers. The expression of ADCY7 was associated with the expression of five MMR genes in most tumors, except for ACC and GBM. (B) The relationship between the expression of ADCY7 and expression of DNMTs in cancers. Strong relationship between ADCY7 expression and expressions of four DNMTs in cancers, except for ESCA, GBM, and ACC. *P-value < 0.05, **P-value < 0.01, ***P-value < 0.001.

Abbreviations: ADCY7, Adenylate cyclase 7; MMR, Mismatch repair; DNMTs, DNA methyltransferases.
Figure 3 The expression of ADCY7 is related to the expression of five MMR genes and DNMTs expression in cancers. (A) The relationship between the expression of ADCY7 and expression of MMR genes in cancers. The expression of ADCY7 was associated with the expression of five MMR genes in most tumors, except for ACC and GBM. (B) The relationship between the expression of ADCY7 and expression of DNMTs in cancers. Strong relationship between ADCY7 expression and expressions of four DNMTs in cancers, except for ESCA, GBM, and ACC. *P-value < 0.05, **P-value < 0.01, ***P-value < 0.001.

Relationship Between ADCY7 Expression and TMB and MSI in Cancers

In the study, the relationship between ADCY7 expression and TMB and MSI were analyzed. shows that ADCY7 expression had positive relevance to TMB in BRCA, COAD, and UCEC, and negative relevance to TMB in LIHC and THCA. In addition, ADCY7 expression was positively related to MSI in CESC, COAD, and LUSC, and negatively related to DLBC, HNSC, and SKCM (). These results indicated that ADCY7 expression was related to immunity in different types of cancer.

Figure 4 Correlation between ADCY7 expression and TMB and MSI in human cancers. (A) Radar map of correlation between ADCY7 expression and TMB. ADCY7 expression had positive relevance to TMB in BRCA, COAD, and UCEC, and negative relevance to TMB in LIHC and THCA. (B) Radar map of correlation between ADCY7 expression and MSI. ADCY7 expression was positively related to MSI in CESC, COAD, and LUSC, and negatively related to DLBC, HNSC, and SKCM.

Abbreviations: ADCY7, Adenylate cyclase 7; TMB, Tumor mutational burden; MSI, Microsatellite instability.
Figure 4 Correlation between ADCY7 expression and TMB and MSI in human cancers. (A) Radar map of correlation between ADCY7 expression and TMB. ADCY7 expression had positive relevance to TMB in BRCA, COAD, and UCEC, and negative relevance to TMB in LIHC and THCA. (B) Radar map of correlation between ADCY7 expression and MSI. ADCY7 expression was positively related to MSI in CESC, COAD, and LUSC, and negatively related to DLBC, HNSC, and SKCM.

Immune Correlation Analysis

The immune cells’ infiltration of tumor tissue is closely in matter of clinical outcomes and can act as therapeutic targets to improve the survival time. Therefore, the study explored the dependency of the immune infiltration level and ADCY7 expression level in 33 cancers using the TIMER database. Results show that ADCY7 expression was related to the infiltration of six types of immune cells in ACC, BRCA, and CHOL ().

Figure 5 Correlation of ADCY7 expression with immune infiltration level in ACC, BRCA, and CHOL. ADCY7 expression is positively correlated with immune infiltration in ACC, BRCA, and CHOL.

Abbreviation: ADCY7, Adenylate cyclase 7.
Figure 5 Correlation of ADCY7 expression with immune infiltration level in ACC, BRCA, and CHOL. ADCY7 expression is positively correlated with immune infiltration in ACC, BRCA, and CHOL.

Survival Prognosis Analysis

The study also explored how ADCY7 expression affected the prognosis (OS, DSS, DFI, and PFI) of patients with cancer. According to the univariate survival analysis of 33 types of cancer, the results revealed that ADCY7 expression was correlated with OS in patients with BLCA, LUSC, MESO, OV, PCPG, and SKCM (). KM curves presented that the overexpression of ADCY7 indicated poorer prognosis in BLCA (P<0.0001, HR=1.07), LUSC (P=0.004, HR=1.04), MESO (P<0.0001, HR=1.07), OV (P=0.0073, HR=1.06), and PCPG (P=0.0035, HR=1.93), but better prognosis in SKCM (P=0.0018, HR=0.96) (). The study also found that ADCY7 expression was associated with patients’ DSS in eight types of cancer, seven of which were associated with worse prognosis, including BLCA (P<0.0001, HR=1.07), COAD (P=0.0071, HR=1.09), LUSC (P=0.00032, HR=1.06), MESO (P=0.0004, HR=1.07), PCPG (P=0.0012, HR=2.38), STAD (P=0.017, HR=1.03), and UVM (P=0.0011, HR=1.49), while one was related to better prognosis in SKCM (P=0.0078, HR=0.95) ( ). Regarding PFI, results revealed that ADCY7 expression was associated with patients’ PFI in three types of cancer, including BLCA, LUSC, and UVM (). KM curves showed that the overexpression of ADCY7 was associated with better prognosis in BLCA (P<0.0001, HR=1.07), LUSC (P<0.0001, HR=1.05), and UVM (P=0.024, HR=1.81) ( ).

Figure 6 The relationship between ADCY7 expression and OS. (A) Forest plots of ADCY7 in thirty-three kinds of cancers. (B) KM curves showed that ADCY7 expression was correlated with OS in 6 kinds of cancers (BLCA, LUSC, MESO, OV, PCPG, and SKCM).

Abbreviations: ADCY7, Adenylate cyclase 7; OS, Overall survival.
Figure 6 The relationship between ADCY7 expression and OS. (A) Forest plots of ADCY7 in thirty-three kinds of cancers. (B) KM curves showed that ADCY7 expression was correlated with OS in 6 kinds of cancers (BLCA, LUSC, MESO, OV, PCPG, and SKCM).

Figure 7 The relationship between ADCY7 expression and DSS. (A) Forest plots of hazard ratios of ADCY7 in thirty-three kinds of cancers. (B) KM curves showed that ADCY7 expression was correlated with DSS in 8 kinds of cancers (BLCA, COAD, PCPG, SKCM, UVM, STAD, MESO, and LUSC).

Abbreviations: ADCY7, Adenylate cyclase 7; DSS, Disease-specific survival.
Figure 7 The relationship between ADCY7 expression and DSS. (A) Forest plots of hazard ratios of ADCY7 in thirty-three kinds of cancers. (B) KM curves showed that ADCY7 expression was correlated with DSS in 8 kinds of cancers (BLCA, COAD, PCPG, SKCM, UVM, STAD, MESO, and LUSC).

Figure 8 The relationship between ADCY7 expression and patients’ PFI. (A) Forest plots of hazard ratios of ADCY7 in thirty-three kinds of cancers. (B) KM curves showed that ADCY7 expression was correlated with PFI in 3 kinds of cancers (BLCA, UVM, and LUSC).

Abbreviations: ADCY7, Adenylate cyclase 7; PFI, Progression-free interval.
Figure 8 The relationship between ADCY7 expression and patients’ PFI. (A) Forest plots of hazard ratios of ADCY7 in thirty-three kinds of cancers. (B) KM curves showed that ADCY7 expression was correlated with PFI in 3 kinds of cancers (BLCA, UVM, and LUSC).

ADCY7 and Its Correlated Gene Functional Enrichment Analyses

To investigate the molecular mechanism of ADCY7 in the occurrence and development of diseases, ADCY7 and its correlated genes were enriched and analyzed by the “clusterProfiler” package in R software (version 3.6.3) for GO and KEGG pathways, respectively. In GO terms, ADCY7 and its correlated gene were mainly enriched in the following: cAMP-mediated signaling, cyclic-nucleotide-mediated signaling, second-messenger-mediated signaling, cAMP-dependent protein kinase complex, heterotrimeric G-protein complex, GTPase complex, ciliary base, catecholamine binding, G-protein beta/gamma−subunit complex binding, peptide binding, and dopamine neurotransmitter receptor activity (). In KEGG terms, ADCY7 and its correlated gene were mainly enriched in the following: neuroactive ligand-receptor interaction, cAMP signaling pathway, dopaminergic synapse, morphine addiction, renin secretion, circadian entrainment, regulation of lipolysis in adipocytes, glutamatergic synapse, gap junction, and serotonergic synapse ().

Figure 9 GO analysis of ADCY7 and its correlated genes.

Abbreviations: ADCY7, Adenylate cyclase 7; GO, Gene Ontology.
Figure 9 GO analysis of ADCY7 and its correlated genes.

Figure 10 KEGG analysis of ADCY7 and its correlated genes.

Abbreviations: ADCY7, Adenylate cyclase 7; KEGG, Kyoto Encyclopedia of Genes and Genomes.
Figure 10 KEGG analysis of ADCY7 and its correlated genes.

Discussion

cAMP serves as a second messenger in the immune system. It regulates cell proliferation, differentiation, and apoptosis.Citation24,Citation25 Adenylate cyclase regulates intracellular cAMP concentration. Studies have shown that ADCY7 can catalyze the generation of cAMP from ATP.Citation26,Citation27 Therefore, ADCY7 may be involved in cell proliferation, differentiation, and apoptosis by regulating cAMP. It has been reported that targeting cAMP signaling pathway can be used in the treatment of chronic lymphocytic leukemia.Citation28 Cancer has many common biological characteristics, such as abnormal cell proliferation and differentiation, invasion, and metastasis. Thus, it can be hypothesized that ADCY7 may be involved in the genesis and progression of cancer. Although the role of ADCY7 in acute myeloid leukemia has been demonstrated, its role in pervasive cancer and its potential to be used as a biomarker are unclear.Citation29 In this study, TCGA and other tumor research databases (including GTEx, CCLE, cBioportal, and TIMER) were used in the pan-cancer analysis of ADCY7. According to the study’s results, it was found that ADCY7 was expressed differently in numerous cancers, and the number of ADCY7 alternations was high in most cancer types. According to the Cox and KM survival analyses, abnormal expression of ADCY7 may be a prognostic factor for some types of cancer. The study also predicted that ADCY7 dysregulation was related to cancer immunity, MMR, and DNA methylation. Therefore, ADCY7 plays a significant role in tumor immunity and might serve as a potential biomarker for prognosis.

Precision medicine has allowed us to identify specific tumor targets or characteristics for individualized treatment and cure.Citation30 The pan-cancer analysis of ADCY7 is significant in investigating differential expressions of ADCY7 and their relative function in cancers.Citation31 Abnormal gene expression is closely related to tumorigenesis and invasion. To the best of our knowledge, this is the first study to find abnormal expressions of ADCY7 in human cancers, including CHOL, GBM, HNSC, KIRC, KIRP, LGG, PAAD, STAD, TGCT, UCEC, ACC, BRCA, COAD, LIHC, LUAD, LUSC, OV, PRAD, SKCM, THCA, and UCS. This suggests that ADCY7 may be related to the occurrence and progression of these cancers.

MMR is closely related to MSI, and the loss of MMR expression leads to the numerous errors in gene replication. This will further affect the MSI and change the DNA repeat sequence, leading to the occurrence of tumors.Citation32 Defective DNA mismatch repair (dMMR) and MSI high (MSI-H) are important indicators of the effectiveness of immune checkpoint inhibitors.Citation33 For patients with dMMR/MSI-H solid tumor, immunotherapy can benefit regardless of the tissue type of tumor origin.Citation34,Citation35 Studies have shown that TMB levels can influence a patient’s response to immune checkpoint inhibitors.Citation36,Citation37 MSI-H was also associated with increased TMB levels.Citation38 Therefore, further investigation of the relationship between ADCY7 expression level and MMR gene expression, MSI, and TMB level in patients with tumors is necessary. In this study, the expression of ADCY7 was related to the expression of 5 MMR genes (MLH1, MSH2, MSH6, PMS2, and EPCAM) in most of the human cancers, with ACC, GBM as the exceptions. Moreover, ADCY7 expression was related to both MSI and TMB in COAD. More studies are needed to confirm whether ADCY7 can be used as a predictor of the efficacy of COAD immunotherapy.

Cancer immunotherapy has attracted the attention of many researchers. Studies have found that the expression level of ADCY7 is related to cancer immunity. Tumor immune microenvironment is mainly composed of tumor-infiltrating lymphocytes and other immune cells, which inhibit or penetrate tumor tension.Citation39,Citation40 Tumor-infiltrating lymphocytes are a new generation of anti-tumor adaptive immune effector cells that can be used as a prognostic marker of cancer.Citation40 Immune escape is involved in the development of cancer.Citation41 Macrophages are involved in anti-tumor immune response.Citation42 These findings suggest that tumor-infiltrating immune cells are involved in tumor progression. In this study, ADCY7 expression was significantly related to six immune-infiltrating cells in ACC, BRCA, and CHOL. Therefore, further investigation of the molecular mechanism of ADCY7 in tumor immunity is warranted.

This study also found that the expression of ADCY7 was related to poor prognosis in several types of tumors especially in BLCA and LUSC. Regardless of OS, PFI, or DSS analysis, the survival time of low ADCY7 expression in patients with BLCA and LUSC was longer than that in patients with BLCA and LUSC patients having high ADCY7 expression. This indicated that ADCY7 is a potential prognostic biomarker for BLCA and LUSC.

Conclusion

This study found that ADCY7 is a potential prognostic biomarker for many tumors, and the expression level of ADCY7 is related to cancer immunity. These results are based on online data, which means that further experiments are needed to investigate the mechanism of ADCY7 in cancer. This study has provided a new perspective on the role of ADCY7 in human cancers.

Ethical Statement

The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. The authors confirm that the data is available from the corresponding author upon reasonable request.

Author Contributions

All authors made a significant contribution to the work reported, whether that is in the conception, study design, execution, acquisition of data, analysis and interpretation, or in all these areas; took part in drafting, revising or critically reviewing the article; gave final approval of the version to be published; have agreed on the journal to which the article has been submitted; and agree to be accountable for all aspects of the work.

Acknowledgments

We would like to thank BulletEdits for English language editing. Yu Zeng and Nanhong Li are co-first authors for this study.

Disclosure

The authors have no conflicts of interest to declare.

Additional information

Funding

The work was supported by the Natural Science Foundation of Guangdong Province (2021A1515011373).

References

  • Hanoune J, Defer N. Regulation and role of adenylyl cyclase isoforms. Annu Rev Pharmacol Toxicol. 2001;41:145–174. doi:10.1146/annurev.pharmtox.41.1.145
  • Baron L, Gombault A, Fanny M, et al. The NLRP3 inflammasome is activated by nanoparticles through ATP, ADP and adenosine. Cell Death Dis. 2015;6:e1629–e.
  • Ludwig M-G, Seuwen K. Characterization of the human adenylyl cyclase gene family: cDNA, gene structure, and tissue distribution of the nine isoforms. J Receptors Signal Transd. 2002;22:79–110. doi:10.1081/RRS-120014589
  • Li C, Xie J, Lu Z, et al. ADCY7 supports development of acute myeloid leukemia. Biochem Biophys Res Commun. 2015;465:47–52.
  • Price T, Brust TF. Adenylyl cyclase 7 and neuropsychiatric disorders: a new target for depression? Pharmacol Res. 2019;143:106–112.
  • Jiang LI, Sternweis PC, Wang JE. Zymosan activates protein kinase A via adenylyl cyclase VII to modulate innate immune responses during inflammation. Mol Immunol. 2013;54:14–22. doi:10.1016/j.molimm.2012.10.027
  • Masugi Y, Tanese K, Emoto K, et al. Overexpression of adenylate cyclase-associated protein 2 is a novel prognostic marker in malignant melanoma. Pathol Int. 2015;65:627–634. doi:10.1111/pin.12351
  • Chen SL, Hu F, Wang DW, et al. Prognosis and regulation of an adenylyl cyclase network in acute myeloid leukemia. Aging. 2020;12:11864–11877. doi:10.18632/aging.103357
  • Sethna F, Feng W, Ding Q, et al. Enhanced expression of ADCY1 underlies aberrant neuronal signalling and behaviour in a syndromic autism model. Nat Commun. 2017;8:14359. doi:10.1038/ncomms14359
  • Bray F, Ferlay J, Soerjomataram I, et al. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68:394–424.
  • Consortium G. The GTEx Consortium atlas of genetic regulatory effects across human tissues. Science. 2020;369:1318–1330.
  • Cerami E, Gao J, Dogrusoz U, et al. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov. 2012;2:401–404. doi:10.1158/2159-8290.CD-12-0095
  • Csiszar A, Balasubramanian P, Tarantini S, et al. Chemically induced carcinogenesis in rodent models of aging: assessing organismal resilience to genotoxic stressors in geroscience research. Geroscience. 2019;41:209–227. doi:10.1007/s11357-019-00064-4
  • Kim S, Wyckoff J, Morris AT, et al. DNA methylation associated with healthy aging of elderly twins. Geroscience. 2018;40:469–484. doi:10.1007/s11357-018-0040-0
  • Masser DR, Hadad N, Porter H, et al. Analysis of DNA modifications in aging research. Geroscience. 2018;40:11–29. doi:10.1007/s11357-018-0005-3
  • Tiffen J, Gallagher SJ, Filipp F, et al. EZH2 cooperates with DNA methylation to downregulate key tumor suppressors and IFN gene signatures in melanoma. J Invest Dermatol. 2020;140(12):2442–2454.e5. doi:10.1016/j.jid.2020.02.042
  • Le DT, Durham JN, Smith KN, et al. Mismatch repair deficiency predicts response of solid tumors to PD-1 blockade. Science. 2017;357:409–413. doi:10.1126/science.aan6733
  • Asaoka Y, Ijichi H, Koike K. PD-1 blockade in tumors with mismatch-repair deficiency. N Engl J Med. 2015;373:1979.
  • Hause RJ, Pritchard CC, Shendure J, et al. Classification and characterization of microsatellite instability across 18 cancer types. Nat Med. 2016;22:1342–1350. doi:10.1038/nm.4191
  • Li T, Fan J, Wang B, et al. TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res. 2017;77:e108–e10. doi:10.1158/0008-5472.CAN-17-0307
  • Szklarczyk D, Gable AL, Lyon D, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47:D607–d13. doi:10.1093/nar/gky1131
  • Yu G, Wang LG, Han Y, et al. clusterProfiler: an R package for comparing biological themes among gene clusters. Omics. 2012;16:284–287. doi:10.1089/omi.2011.0118
  • Liu S, Xie X, Lei H, et al. Identification of key circRNAs/lncRNAs/miRNAs/mRNAs and pathways in preeclampsia using bioinformatics analysis. Med Sci Monit. 2019;25:1679–1693. doi:10.12659/MSM.912801
  • Kammer GM. The adenylate cyclase-cAMP-protein kinase A pathway and regulation of the immune response. Immunol Today. 1988;9:222–229. doi:10.1016/0167-5699(88)91220-0
  • Castro A, Jerez MJ, Gil C, et al. Cyclic nucleotide phosphodiesterases and their role in immunomodulatory responses: advances in the development of specific phosphodiesterase inhibitors. Med Res Rev. 2005;25:229–244. doi:10.1002/med.20020
  • Patel TB, Du Z, Pierre S, et al. Molecular biological approaches to unravel adenylyl cyclase signaling and function. Gene. 2001;269:13–25. doi:10.1016/S0378-1119(01)00448-6
  • Baron L, Gombault A, Fanny M, et al. The NLRP3 inflammasome is activated by nanoparticles through ATP, ADP and adenosine. Cell Death Dis. 2015;6:e1629.
  • Murray F, Insel PA. Targeting cAMP in chronic lymphocytic leukemia: a pathway-dependent approach for the treatment of leukemia and lymphoma. Expert Opin Ther Targets. 2013;17:937–949. doi:10.1517/14728222.2013.798304
  • Li C, Xie J, Lu Z, et al. ADCY7 supports development of acute myeloid leukemia. Biochem Biophys Res Commun. 2015;465:47–52. doi:10.1016/j.bbrc.2015.07.123
  • Lancet T. Moving toward precision medicine. Lancet. 2011;378:1678. doi:10.1016/S0140-6736(11)61725-X
  • Cao Z, Zhang S. An integrative and comparative study of pan-cancer transcriptomes reveals distinct cancer common and specific signatures. Sci Rep. 2016;6:33398. doi:10.1038/srep33398
  • Armaghany T, Wilson JD, Chu Q, et al. Genetic alterations in colorectal cancer. Gastrointest Cancer Res. 2012;5:19–27.
  • Baretti M, Le DT. DNA mismatch repair in cancer. Pharmacol Ther. 2018;189:45–62. doi:10.1016/j.pharmthera.2018.04.004
  • Li K, Luo H, Huang L, et al. Microsatellite instability: a review of what the oncologist should know. Cancer Cell Int. 2020;20:16. doi:10.1186/s12935-019-1091-8
  • Herbst RS, Baas P, Perez-Gracia JL, et al. Use of archival versus newly collected tumor samples for assessing PD-L1 expression and overall survival: an updated analysis of KEYNOTE-010 trial. Ann Oncol. 2019;30:281–289. doi:10.1093/annonc/mdy545
  • Wu HX, Wang ZX, Zhao Q, et al. Tumor mutational and indel burden: a systematic pan-cancer evaluation as prognostic biomarkers. Ann Transl Med. 2019;7:640. doi:10.21037/atm.2019.10.116
  • Havel JJ, Chowell D, Chan TA. The evolving landscape of biomarkers for checkpoint inhibitor immunotherapy. Nat Rev Cancer. 2019;19:133–150.
  • Chalmers ZR, Connelly CF, Fabrizio D, et al. Analysis of 100,000 human cancer genomes reveals the landscape of tumor mutational burden. Genome Med. 2017;9:34. doi:10.1186/s13073-017-0424-2
  • Junttila MR, de Sauvage FJ. Influence of tumour micro-environment heterogeneity on therapeutic response. Nature. 2013;501:346–354. doi:10.1038/nature12626
  • Woo JR, Liss MA, Muldong MT, et al. Tumor infiltrating B-cells are increased in prostate cancer tissue. J Transl Med. 2014;12:30. doi:10.1186/1479-5876-12-30
  • Ju Q, Li XM, Zhang H, et al. BRCA1-associated protein is a potential prognostic biomarker and is correlated with immune infiltration in liver hepatocellular carcinoma: a pan-cancer analysis. Front Mol Biosci. 2020;7:573619. doi:10.3389/fmolb.2020.573619
  • Shan X, Zhang C, Wang Z, et al. Prognostic value of a nine-gene signature in glioma patients based on tumor-associated macrophages expression profiling. Clin Immunol. 2020;216:108430. doi:10.1016/j.clim.2020.108430