ORIGINAL RESEARCH article

Front. Genet., 23 June 2021

Sec. Systems Biology Archive

Volume 12 - 2021 | https://doi.org/10.3389/fgene.2021.595436

Association of CLDN6 and CLDN10 With Immune Microenvironment in Ovarian Cancer: A Study of the Claudin Family

  • 1. Cancer Biology Research Center (Key Laboratory of the Ministry of Education), Department of Obstetrics and Gynecology, Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China

  • 2. Department of Obstetrics and Gynecology, Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China

  • 3. Institute of Pathology, Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China

Abstract

Background:

The claudin family is a group of transmembrane proteins related to tight junctions. While their involvement in cancer has been studied extensively, their relationship with the tumor immune microenvironment remains poorly understood. In this research, we focused on genes related to the prognosis of ovarian cancer and explored their relationship with the tumor immune microenvironment.

Methods:

The cBioPortal for Cancer Genomics database was used to obtain the genetic variation pattern of the claudin family in ovarian cancer. The ONCOMINE and Gene Expression Profiling Interactive Analysis (GEPIA) databases were used to explore the mRNA expression of claudins in cancers. The prognostic potential of these genes was examined via the Kaplan-Meier plotter. The enrichment of immunological signatures was determined by gene set enrichment analysis (GSEA). The correlations between claudins and the tumor immune microenvironment in ovarian cancer were investigated via the Tumor Immune Estimation Resource (TIMER).

Results:

Claudin genes were altered in 363 (62%) of queried patients/samples. Abnormal expression levels of claudins were observed in various cancers. Among them, CLDN3, CLDN4, CLDN6, CLDN10, CLDN15, and CLDN16 were significantly correlated with overall survival in patients with ovarian cancer. GSEA revealed that CLDN6 and CLDN10 were significantly enriched in immunological signatures of B cell, CD4 T cell, and CD8 T cell. Furthermore, CLDN6 and CLDN10 were negatively correlated and positively correlated, respectively, with immune cell infiltration in ovarian cancer. The expression levels of CLDN6 and CLDN10 were also negatively correlated and positively correlated, respectively, with various gene markers of immune cells in ovarian cancer. Thus, CLDN6 and CLDN10 may participate in immune cell infiltration in ovarian cancer, and these mechanisms may be the reason for poor prognosis.

Conclusion:

Our study showed that CLDN6 and CLDN10 were prognostic biomarkers correlated with the immune microenvironment in ovarian cancer. These results reveal new roles for CLDN6 and CLDN10 as potential therapeutic targets in the treatment of ovarian cancer.

Introduction

Ovarian cancer is the most lethal gynecological cancer among women (Siegel et al., 2020). Although surgical techniques and combined chemotherapy applications have progressed since the 1970s, the 5 year survival rate of advanced ovarian cancer is only 40–45% (Henderson et al., 2018). Therefore, improved treatment of ovarian cancer remains an urgent issue. Immunotherapy is an emerging treatment for several solid tumors, which shows improved outcomes in patients. With the application of various immune-based interventions in ovarian cancer, immunotherapy has been proven useful in advanced disease (Bogani et al., 2020).

The claudin (CLDN) family consists of more than 20 transmembrane proteins, which are major components of tight junctions. They serve as a physical barrier to prevent molecules from passing freely through the paracellular space between epithelial or endothelial cell sheets and also play critical roles in maintaining cell polarity and signal transductions (Weinstein et al., 1976; Wodarz, 2000; Tsukita et al., 2001; Kirschner et al., 2013). Previous research has recognized various claudin gene expression patterns and identified several genes dysregulated in cancers (Hewitt et al., 2006). These genes play roles in the tumorigenesis of solid tumors (Swisshelm et al., 2005; Hagen, 2019) and represent promising targets for cancer detection, prognosis, and therapy (Morin, 2005). However, the relationship between claudins and the tumor immune microenvironment has not yet been elucidated. This study comprehensively analyzed claudin expression in ovarian cancer and further explored the relationship between claudins and the immune microenvironment.

Materials and Methods

cBioPortal

The cBioPortal database1 (Cerami et al., 2012; Gao et al., 2013) is an open platform for cancer genomics analysis. In total, 585 samples of ovarian serous cystadenocarcinoma (The Cancer Genome Atlas (TCGA), Pan-Cancer Atlas) were used for genetic variation analyses through the cBioPortal.

ONCOMINE Database Analysis

Claudin expression levels in various cancers were analyzed via the ONCOMINE database2 (Rhodes et al., 2007), which includes more than 35 types of cancer and normal samples.

Gene Expression Profiling Interactive Analysis (GEPIA)

GEPIA v23 (Tang et al., 2017) is used to analyze the RNA sequencing expression data of 9736 tumors and 8587 normal samples from the TCGA and GTEx projects using a standard processing pipeline. The expression profile of the claudins in ovarian cancer was explored via GEPIA v2. The p-value cutoff was 0.05 and | log2FC| cutoff was 1.5.

Kaplan-Meier Plotter Database Analysis

The Kaplan-Meier plotter4 (Gyorffy et al., 2012) assesses the effects of 54,000 genes on survival in 21 cancer types. The largest datasets include breast (n = 6234), ovarian (n = 2190), lung (n = 3452), and gastric (n = 1440) cancer. The system includes gene chip and RNA-seq data-sources from the Gene Expression Omnibus (GEO), European Genome-Phenome Archive (EGA), and TCGA databases. The prognostic significance of claudins in ovarian cancer was analyzed via the online database.

Tumor Immune Estimation Resource (TIMER)

TIMER5 (Li et al., 2017) allows comprehensive analysis of tumor-infiltrating immune cells. The correlation between claudin expression and immune cell infiltration was analyzed using this database. TIMER v2, an updated and enhanced version of TIMER, was used to analyze immune infiltration across diverse cancer types.

Statistical Analyses

The expression levels of claudins are presented as mean ± standard deviation (SD). Kaplan-Meier survival curves were established based on the log-rank test. The hazard ratio (HR) was determined using the Cox model. Spearman correlation was used for correlation analysis. A p-value of < 0.05 was considered to be significant.

Results

Gene Variation of Claudins in Ovarian Cancer

Twenty-four reviewed proteins of the claudin family were obtained from the UniProt Knowledgebase (UniProtKB)6 (Table 1) [an additional file shows this in more detail (see Table 1)]. Firstly, we investigated the genetic variation of the claudin family in ovarian cancer using the cBioPortal for Cancer Genomics. Twenty-four genes were queried in 585 samples of ovarian serous cystadenocarcinoma (TCGA, Pan-Cancer Atlas). Figure 1A shows the alteration frequency of genetic variation in serous ovarian cancer. As shown in Figure 1B, the queried genes were altered in 363 (62%) queried patients/samples. The top three gene variations were CLDN11 (24%), CLDN16 (22%), and CLDN1 (16%). Differences in overall survival (OS) between the altered and unaltered groups were compared using the Kruskal-Wallis test. We found that OS was reduced in the altered group compared to the unaltered group (p = 7.981e-3) (Figure 1C). Previous studies have shown that the claudin family is dysregulated in a variety of tumors and is involved in diagnosis, tumorigenesis, and prognosis (Zhang et al., 2013; Barros-Filho et al., 2015; Zhou et al., 2018). Thus, the claudin family is worthy of further research in ovarian cancer.

TABLE 1

EntryStatusGene namesProtein namesOrganism
O95832ReviewedCLDN1Claudin-1 (Senescence-associated epithelial membrane protein)Homo sapiens
P78369ReviewedCLDN10Claudin-10 (Oligodendrocyte-specific protein-like) (OSP-like)Homo sapiens
O75508ReviewedCLDN11Claudin-11 (Oligodendrocyte-specific protein)Homo sapiens
P56749ReviewedCLDN12Claudin-12Homo sapiens
O95500ReviewedCLDN14Claudin-14Homo sapiens
P56746ReviewedCLDN15Claudin-15Homo sapiens
Q9Y5I7ReviewedCLDN16Claudin-16 (Paracellin-1) (PCLN-1)Homo sapiens
P56750ReviewedCLDN17Claudin-17Homo sapiens
P56856ReviewedCLDN18Claudin-18Homo sapiens
Q8N6F1ReviewedCLDN19Claudin-19Homo sapiens
P57739ReviewedCLDN2Claudin-2 (SP82)Homo sapiens
P56880ReviewedCLDN20Claudin-20Homo sapiens
Q8N7P3ReviewedCLDN22Claudin-22Homo sapiens
Q96B33ReviewedCLDN23Claudin-23Homo sapiens
A6NM45ReviewedCLDN24/CLDN21Putative claudin-24 (Claudin-21)Homo sapiens
C9JDP6ReviewedCLDN25Putative claudin-25Homo sapiens
O15551ReviewedCLDN3Claudin-3 (CPE-receptor 2)Homo sapiens
H7C241ReviewedCLDN34Claudin-34Homo sapiens
O14493ReviewedCLDN4Claudin-4 (CPE-receptor)Homo sapiens
O00501ReviewedCLDN5Claudin-5 (Transmembrane protein deleted in VCFS) (TMDVCF)Homo sapiens
P56747ReviewedCLDN6Claudin-6 (Skullin)Homo sapiens
O95471ReviewedCLDN7Claudin-7Homo sapiens
P56748ReviewedCLDN8Claudin-8Homo sapiens
O95484ReviewedCLDN9Claudin-9Homo sapiens

Twenty-four reviewed proteins of claudin family from the UniProtKB.

FIGURE 1

Expression of Claudins Is Dysregulated in Various Cancers

To explore the mRNA expression of the claudin family, we investigated the expression profiles of claudins in various cancers via the ONCOMINE database. The thresholds were: p-value of 0.05, fold change of 1.5, and gene rank of all. Significant analyses are shown in Supplementary Figure 1 (those with < 3 significant analyses were not considered). Results showed that most claudins were dysregulated in various cancers. To verify the expression of claudins in ovarian cancer, GEPIA2 was used to analyze mRNA expression in TCGA and GTEx samples. The | Log2FC| cutoff was set to 1.5 and the p-value cutoff was set to 0.05. As shown in Figure 2, eight genes were overexpressed in ovarian cancer samples compared with normal tissue samples and included CLDN1, CLDN3, CLDN4, CLDN6, CLDN7, CLDN9, CLDN10, and CLDN16. Furthermore, three genes showed low expression in the ovarian cancer samples compared with normal tissue samples and included CLDN5, CLDN11, and CLDN15.

FIGURE 2

Correlation of Claudin Expression With Ovarian Cancer Prognosis

To identify genes with clinical significance, we studied the relationship between differentially expressed genes (DEGs) and ovarian cancer patient prognosis using the Kaplan-Meier plotter. As shown in Figure 3, overexpressed genes CLDN3, CLDN4, CLDN6, and CLDN16 were significantly correlated with poor OS (Figure 3A) and progression free survival (PFS) (Figure 3B). In addition, high expression of CLDN10 and CLDN15 were predictive of a good prognosis in ovarian cancer patients (Figures 3C,D). Surprisingly, CLDN10 was overexpressed in cancer, but patients with high expression of CLDN10 showed good OS (HR = 0.73, logrank P = 1.6e-06), PFS (HR = 0.83, logrank P = 0.0067), and post progression survival (PPS, HR = 0.73, logrank P = 0.00029). These results are somewhat counterintuitive, and the underlying mechanism requires further exploration.

FIGURE 3

TCGA projects have identified four molecular subtypes of high-grade serous ovarian carcinoma (HGSOC) (Cancer Genome Atlas Research Network, 2011): (i) the differentiated subtype; (ii) the immunoreactive subtype; (iii) the mesenchymal subtype; and (iv) the proliferative subtype. Among them, T-cell chemokine ligands, CXCL11 and CXCL10, and the receptor, CXCR3, characterized the immunoreactive subtype. Then, Thorsson et al. (2018) developed a global immune classification of solid tumors based on the transcriptomic profiles of 33 cancer types. They identified six distinct immune subtypes: C1 (Wound healing); C2 (IFN-γ dominant); C3 (Inflammatory); C4 (Lymphocyte depleted); C5 (Immunologically quiet); C6 (TGF-β dominant). These six categories represent features of the tumor microenvironment (Charoentong et al., 2017). In this research, we explored the relationships between the expression of differentially expressed genes related to prognosis and molecular subtypes or immune subtypes of ovarian cancer via the TISIDB (Ru et al., 2019). The Kruskal-Wallis test was used. As Supplementary Figure 2 shows, claudins including CLDN3, CLDN6, CLDN10, and CLDN15 are differentially expressed in different immune subtypes. And, claudins including CLDN3, CLDN4, CLDN6, CLDN10, and CLDN16 are differentially expressed in different molecular subtypes (Supplementary Figure 3). Among them, CLDN6 is relatively low expression, and CLDN10 is relatively high expression in the immunoreactive subtype.

GSEA of Immunological Signature Gene Sets

To characterize the potential function of claudins, GSEA was performed using gene expression data from TCGA ovarian cancer patients. Immunological signature gene sets were used. As shown in Figure 4, CLDN6 and CLDN10 were related to the effector differentiation of B cell, CD4 T cell, and CD8 T cell.

FIGURE 4

Correlation Analyses Between Claudins and Tumor Immune Microenvironment

To understand the role of claudins in immunity, we downloaded 379 RNA-seq FPKM (Fragments per kilobase per million) data of ovarian cancer from TCGA. Subsequently, the FPKM was converted to TPM (transcripts per million) (Li et al., 2010). The ESTIMATE algorithm (Yoshihara et al., 2013) was used to predict tumor purity based on TCGA ovarian cancer samples. Then, the relationship between claudin expression and the immune microenvironment was explored. As shown in Figure 5A, a significant negative correlation between CLDN6 expression and the immune score was observed (Spearman correlation = −0.23, p < 0.001). A significant positive correlation between CLDN10 expression and immune score (spearman correlation = 0.21, p < 0.001) was observed (Figure 5B). However, the expression levels of CLDN6 and CLDN10 were not correlated with the stromal score.

FIGURE 5

We next examined the relationship between immune cell infiltration and claudin expression. RNA-seq TPM data (n = 379) from TCGA ovarian cancer were used to assess 22 immune cells subtype concentrations with the CIBERSORT algorithm (Newman et al., 2019). TCGA samples were grouped by the median values of CLDN6 and CLDN10, respectively. Activated dendritic cells differed significantly between the CLDN6_high and CLDN6_low groups. Several cell types were significantly different between the CLDN10_high and CLDN10_low group, including naïve B cells, memory B cells, naïve CD4 T cells, CD4 memory-activated T cells, monocytes, M1 macrophages, and activated dendritic cells (Figure 5C).

The microarray expression values of ovarian cancer were used to calculate the abundances of six immune infiltrates (B cells, CD4+ T cells, CD8+ T cells, Neutrophils, Macrophages, and Dendritic cells) via the TIMER algorithm (Yoshihara et al., 2013). The gene expression levels correlated with tumor purity are displayed in the left-most panel (Figures 6A,B). Results showed that CLDN6 expression was negatively correlated with infiltration of B cell (partial correlation = −0.284, p = 2.21e-10), CD8+ T cells (partial correlation = −0.254, p = 1.64e-08), neutrophils (partial correlation = −0.152, p = 8.29e-04), and dendritic cells (partial correlation = −0.182, p = 6.31e-05) (Figure 6A). In contrast, there was a small but significant positive correlation between CLDN10 expression and infiltration of neutrophils (partial correlation = 0.185, p = 4.66e-05), and dendritic cells (partial correlation = 0.153, p = 7.74e-04) (Figure 6B).

FIGURE 6

To more accurately describe the relationship between gene expression and immune cell infiltration, we used the TIMER, CIBERSORT, quanTIseq, xCell, MCP-counter, and EPIC algorithms to assess the immune infiltration in tumor tissue (Sturm et al., 2019). TIMER2 provides a platform to analyze immune infiltrates across diverse cancer types based on available TCGA RNA-seq data (Li et al., 2016; Li T. et al., 2020). The correlations between claudin expression (CLDN6 and CLDN10) and immune cell infiltration in ovarian cancer are shown in Table 2. As seen in Figure 6C, CLDN6 was negatively correlated with immune cell infiltration, including that of B cells, CD8+ T cells, effector memory CD4+ T cells, M1 macrophages, and myeloid dendritic cells. In contrast, CLDN10 was positively correlated with immune cell infiltration, including that of B cells, CD8+ T cells, effector memory CD4+ T cells, M1 macrophages, and myeloid dendritic cells (Figure 6D). Relevant evidence suggests that cancer-associated fibroblasts (CAFs) play an important role in the progression of ovarian cancer (Mhawech-Fauceglia et al., 2014; Leung et al., 2018). Interestingly, here, CAFs also showed a positive correlation with CLDN6 expression, but a negative correlation with CLDN10 expression. In ovarian cancer, increased infiltration of tumor-infiltrating lymphocytes (TILs), and more specifically CD8+ T cells, have been proven to be associated with improved clinical outcomes (Sato et al., 2005; Hamanishi et al., 2007; Ovarian Tumor Tissue Analysis et al., 2017). These results suggest that CLDN6 and CLDN10 may participate in immune cell infiltration in ovarian cancer, and these mechanisms may be the reason for poor prognosis.

TABLE 2

CLDN6
CLDN10
CancerInfiltratesrhopadj.prhopadj.p
OV (n = 303)B cell memory_CIBERSORT−0.0180.7770.9214−0.1938***
OV (n = 303)B cell memory_CIBERSORT-ABS−0.01850.77130.9214−0.1795***
OV (n = 303)B cell memory_XCELL−0.03860.54460.78550.0910.15210.3381
OV (n = 303)B cell naive_CIBERSORT0.00530.93430.98950.255******
OV (n = 303)B cell naive_CIBERSORT-ABS−0.00580.92720.98950.2577******
OV (n = 303)B cell naive_XCELL0.09150.150.4803−0.142*0.0952
OV (n = 303)B cell plasma_CIBERSORT0.11640.06660.3075−0.03370.59630.7755
OV (n = 303)B cell plasma_CIBERSORT-ABS0.07410.24430.57680.00360.95520.9837
OV (n = 303)B cell plasma_XCELL0.040.53020.7759−0.120.05870.1821
OV (n = 303)B cell_EPIC0.0450.48010.7541−0.149*0.0782
OV (n = 303)B cell_MCPCOUNTER0.2482*****−0.08360.18880.3814
OV (n = 303)B cell_QUANTISEQ0.11530.06940.3139−0.11770.06360.1866
OV (n = 303)B cell_TIMER−0.3021******0.2164*****
OV (n = 303)B cell_XCELL−0.1283*0.26160.07560.23450.4401
OV (n = 303)Cancer associated fibroblast_EPIC0.1377*0.1353−0.09070.15370.4081
OV (n = 303)Cancer associated fibroblast_MCPCOUNTER0.1594*0.0746−0.09550.1330.3766
OV (n = 303)Cancer associated fibroblast_TIDE0.197***−0.178***
OV (n = 303)Cancer associated fibroblast_XCELL0.1913***−0.12010.05850.2122
OV (n = 303)Class-switched memory B cell_XCELL−0.10730.0910.37470.10940.0850.2267
OV (n = 303)Common lymphoid progenitor_XCELL−0.06280.32350.65960.07950.21120.4607
OV (n = 303)Common myeloid progenitor_XCELL−0.1444*0.1390.03330.60090.8165
OV (n = 303)Endothelial cell_EPIC0.0920.14780.4554−0.11350.07380.2627
OV (n = 303)Endothelial cell_MCPCOUNTER0.15*0.1218−0.11090.08070.2771
OV (n = 303)Endothelial cell_XCELL0.09230.14660.4554−0.08930.160.403
OV (n = 303)Eosinophil_CIBERSORT0.1312*0.1921−0.0060.92550.9687
OV (n = 303)Eosinophil_CIBERSORT-ABS0.1299*0.1983−0.00540.93230.9707
OV (n = 303)Eosinophil_XCELL0.04720.45880.7698−0.09080.15310.3919
OV (n = 303)Granulocyte-monocyte progenitor_XCELL0.04230.50610.78730.00610.92360.9687
OV (n = 303)Hematopoietic stem cell_XCELL0.07040.26850.6192−0.1648**0.0568
OV (n = 303)Macrophage M0_CIBERSORT0.120.05860.2045−0.1693***
OV (n = 303)Macrophage M0_CIBERSORT-ABS0.08540.17910.431−0.12190.05460.168
OV (n = 303)Macrophage M1_CIBERSORT−0.1565*0.08120.1868***
OV (n = 303)Macrophage M1_CIBERSORT-ABS−0.12010.05850.20450.1764***
OV (n = 303)Macrophage M1_QUANTISEQ−0.11150.07920.25410.1631***
OV (n = 303)Macrophage M1_XCELL−0.2436*****0.2096*****
OV (n = 303)Macrophage M2_CIBERSORT−0.1332*0.14810.09460.13660.3176
OV (n = 303)Macrophage M2_CIBERSORT-ABS−0.12010.05850.20450.1292*0.1388
OV (n = 303)Macrophage M2_QUANTISEQ−0.06320.32070.60290.12330.05210.1619
OV (n = 303)Macrophage M2_TIDE0.3074******−0.2819******
OV (n = 303)Macrophage M2_XCELL−0.2827******0.09920.11830.2886
OV (n = 303)Macrophage/Monocyte_MCPCOUNTER−0.1563*0.08120.06750.28840.5842
OV (n = 303)Macrophage/Monocyte_MCPCOUNTER−0.1563*0.11150.06750.28840.5244
OV (n = 303)Macrophage_EPIC−0.1983***0.1515*0.0698
OV (n = 303)Macrophage_TIMER0.03710.56020.7984−0.1785***
OV (n = 303)Macrophage_XCELL−0.2767******0.1879***
OV (n = 303)Mast cell activated_CIBERSORT0.01350.83250.9299−0.02710.66990.8355
OV (n = 303)Mast cell activated_CIBERSORT-ABS0.01180.85270.9352−0.02840.65550.8323
OV (n = 303)Mast cell resting_CIBERSORT−0.06450.31060.650.07650.22890.4775
OV (n = 303)Mast cell resting_CIBERSORT-ABS−0.07750.2230.56260.09790.12330.3433
OV (n = 303)Mast cell_XCELL−0.1516*0.1157−0.06980.27230.5282
OV (n = 303)MDSC_TIDE0.3588******−0.1393*0.1339
OV (n = 303)Monocyte_CIBERSORT0.04490.4810.77760.07390.24540.5578
OV (n = 303)Monocyte_CIBERSORT-ABS−0.00030.99660.99660.1240.05070.2355
OV (n = 303)Monocyte_MCPCOUNTER−0.1563*0.11150.06750.28840.5842
OV (n = 303)Monocyte_QUANTISEQ−0.3974******0.1651**0.0626
OV (n = 303)Monocyte_XCELL−0.11090.08070.33180.08240.1950.5043
OV (n = 303)Myeloid dendritic cell activated_CIBERSORT−0.1643**0.05590.1554*0.069
OV (n = 303)Myeloid dendritic cell activated_CIBERSORT-ABS−0.1626*0.05730.1618*0.0564
OV (n = 303)Myeloid dendritic cell activated_XCELL−0.2327*****0.1691***
OV (n = 303)Myeloid dendritic cell resting_CIBERSORT−0.03710.56050.7955−0.05460.39080.635
OV (n = 303)Myeloid dendritic cell resting_CIBERSORT-ABS−0.03670.56420.7962−0.04750.45510.6843
OV (n = 303)Myeloid dendritic cell_MCPCOUNTER−0.10320.10440.29890.02760.66520.8057
OV (n = 303)Myeloid dendritic cell_QUANTISEQ0.363******−0.1552*0.0693
OV (n = 303)Myeloid dendritic cell_TIMER−0.3143******0.2908******
OV (n = 303)Myeloid dendritic cell_XCELL−0.11960.05950.21380.1565*0.0675
OV (n = 303)Neutrophil_CIBERSORT−0.10290.10530.41270.11140.07930.2453
OV (n = 303)Neutrophil_CIBERSORT-ABS−0.09510.13450.46050.10720.09130.2681
OV (n = 303)Neutrophil_MCPCOUNTER−0.00170.97860.9929−0.03670.56390.7514
OV (n = 303)Neutrophil_QUANTISEQ0.1785**0.0595−0.02070.74470.863
OV (n = 303)Neutrophil_TIMER−0.07240.25520.610.06140.33480.5858
OV (n = 303)Neutrophil_XCELL−0.08690.17140.51220.08420.18510.418
OV (n = 303)NK cell activated_CIBERSORT−0.02630.67960.86630.02960.64230.8424
OV (n = 303)NK cell activated_CIBERSORT-ABS−0.04040.52560.77860.120.05870.2122
OV (n = 303)NK cell resting_CIBERSORT−0.10090.11240.3225−0.02460.69890.8788
OV (n = 303)NK cell resting_CIBERSORT-ABS−0.11090.08080.266−0.02260.72240.8908
OV (n = 303)NK cell_EPIC−0.1815***0.11490.07030.2474
OV (n = 303)NK cell_MCPCOUNTER−0.1553*0.08480.1402*0.12
OV (n = 303)NK cell_QUANTISEQ−0.05560.38210.67810.04110.5190.7789
OV (n = 303)NK cell_XCELL−0.08240.19510.44910.07990.20870.4765
OV (n = 303)Plasmacytoid dendritic cell_XCELL−0.208****0.2213*****
OV (n = 303)T cell CD4+ (non-regulatory)_QUANTISEQ−0.05360.39980.7259−0.06380.31560.5912
OV (n = 303)T cell CD4+ (non-regulatory)_XCELL0.00770.90320.9663−0.07230.25550.5347
OV (n = 303)T cell CD4+ central memory_XCELL0.04560.47360.78110.03440.58920.8122
OV (n = 303)T cell CD4+ effector memory_XCELL−0.1513*0.11090.1302*0.1625
OV (n = 303)T cell CD4+ memory activated_CIBERSORT−0.00470.94110.97980.05380.39820.6743
OV (n = 303)T cell CD4+ memory activated_CIBERSORT-ABS−0.00410.94850.97980.05260.4090.6835
OV (n = 303)T cell CD4+ memory resting_CIBERSORT0.10470.09940.3290.0150.81410.9242
OV (n = 303)T cell CD4+ memory resting_CIBERSORT-ABS0.00140.98270.9920.09430.13780.3757
OV (n = 303)T cell CD4+ memory_XCELL0.02530.69160.8970.06930.27620.5595
OV (n = 303)T cell CD4+ naive_CIBERSORT0.1349*0.1741−0.1428*0.1147
OV (n = 303)T cell CD4+ naive_CIBERSORT-ABS0.1349*0.1741−0.1428*0.1147
OV (n = 303)T cell CD4+ naive_XCELL−0.1611*0.08280.11010.0830.2652
OV (n = 303)T cell CD4+ Th1_XCELL−0.1385*0.16080.04990.43280.7009
OV (n = 303)T cell CD4+ Th2_XCELL0.06250.32630.65060.07660.22870.522
OV (n = 303)T cell CD4+ _EPIC0.04280.50140.8099−0.01480.81680.9242
OV (n = 303)T cell CD4+ _TIMER0.11490.07030.2735−0.00580.92730.9753
OV (n = 303)T cell CD8+ central memory_XCELL−0.1749***0.1568*0.0801
OV (n = 303)T cell CD8+ effector memory_XCELL0.08580.1770.46880.07960.21070.4441
OV (n = 303)T cell CD8+ naive_XCELL0.1924***−0.11910.06060.2145
OV (n = 303)T cell CD8+ _CIBERSORT−0.05340.40120.68290.03010.63660.8318
OV (n = 303)T cell CD8+ _CIBERSORT-ABS−0.04530.47650.70860.07020.26950.5033
OV (n = 303)T cell CD8+ _EPIC0.04340.49510.7166−0.05420.39440.6552
OV (n = 303)T cell CD8+ _MCPCOUNTER−0.03220.6130.79090.09250.14550.3528
OV (n = 303)T cell CD8+ _QUANTISEQ−0.2023***0.1851***
OV (n = 303)T cell CD8+ _TIMER−0.1707***0.1363*0.139
OV (n = 303)T cell CD8+ _XCELL−0.05440.39230.6765−0.00780.90280.9629
OV (n = 303)T cell follicular helper_CIBERSORT−0.0360.57160.82550.00320.96050.9889
OV (n = 303)T cell follicular helper_CIBERSORT-ABS−0.06180.33160.70460.0580.36180.6466
OV (n = 303)T cell gamma delta_CIBERSORT−0.02810.65910.8771−0.07380.24580.5578
OV (n = 303)T cell gamma delta_CIBERSORT-ABS−0.02760.66420.8771−0.07350.24810.5578
OV (n = 303)T cell gamma delta_XCELL−0.05450.39180.74310.030.63740.8533
OV (n = 303)T cell NK_XCELL−0.1745**0.064−0.0010.98690.9937
OV (n = 303)T cell regulatory (Tregs)_CIBERSORT−0.00560.92990.9886−0.04170.51230.7769
OV (n = 303)T cell regulatory (Tregs)_CIBERSORT-ABS−0.02780.66220.8546−0.0060.92480.971
OV (n = 303)T cell regulatory (Tregs)_QUANTISEQ−0.0010.98730.99980.1678***
OV (n = 303)T cell regulatory (Tregs)_XCELL0.06830.2830.5750.04020.52760.783

Correlation analysis between claudins and immune infiltration in ovarian cancer via TIMER2.0.

*P < 0.05;**P < 0.01;***P < 0.001.

Relationship Between Claudin Expression and Gene Markers of Immune Cells

To further illustrate the correlation between claudins and the immune microenvironment, we analyzed the relationship between CLDN6 and CLDN10 expression and gene markers of various immune cells in ovarian cancer (TIMER2 database), including B cells, T cells (general), CD8+ T cells, macrophages, dendritic cells, neutrophils, monocytes, natural killer (NK) cells, and regulatory T cells (Tregs) (Table 3). Purity-adjusted correlation heatmaps are shown in Supplementary Figure 4. After correlation adjustment by purity, CLDN6 expression was negatively correlated with most gene markers of dendritic cells, M1 macrophages, monocytes, NK cells, and tumor-associated macrophages (TAMs) in ovarian cancer. In contrast, CLDN10 expression was positively correlated with gene markers of dendritic cells, T cells (general), and TAMs in ovarian cancer.

TABLE 3

CLDN6
CLDN10
CancerImmune cellsGene markersrhopadj.prhopadj.p
OV (n = 303)B cellCD190.1220.0520.189−0.0750.2680.477
OV (n = 303)B cellCD79A0.0220.6920.853−0.0630.3070.521
OV (n = 303)CD8+ T cellCD8A−0.1030.1030.3050.09770.1210.298
OV (n = 303)CD8+ T cellCD8B−0.0320.6130.7980.09250.1450.336
OV (n = 303)DCCD1C−0.158*0.0980.08640.1720.467
OV (n = 303)DCHLA-DPA1−0.253*****0.2298*****
OV (n = 303)DCHLA-DPB1−0.3******0.2535******
OV (n = 303)DCHLA-DQB1−0.224*****0.2259*****
OV (n = 303)DCHLA-DRA−0.325******0.2428*****
OV (n = 303)DCITGAX−0.182***0.08590.1780.469
OV (n = 303)DCNRP10.1252*0.235−0.0040.9960.997
OV (n = 303)M1 MacrophageIRF5−0.186***0.08960.1570.342
OV (n = 303)M1 MacrophageNOS20.1436*0.106−0.0330.5450.757
OV (n = 303)M1 MacrophagePTGS20.09610.1350.3470.00930.8860.942
OV (n = 303)M2 MacrophageCD163−0.1060.0960.2890.06460.310.529
OV (n = 303)M2 MacrophageMS4A4A−0.1120.0720.2380.11470.0770.206
OV (n = 303)M2 MacrophageVSIG4−0.152*0.0860.07680.2240.433
OV (n = 303)MonocyteCD86−0.222*****0.1457*0.084
OV (n = 303)MonocyteCSF1R−0.196***0.07170.2560.473
OV (n = 303)NK cellKIR2DL1−0.0010.9240.9860.09910.1170.388
OV (n = 303)NK cellKIR2DL3−0.226*****0.1527*0.096
OV (n = 303)NK cellKIR2DL4−0.258*****0.1563*0.08
OV (n = 303)NK cellKIR2DS4−0.0970.1210.3910.08470.1850.475
OV (n = 303)NK cellKIR3DL10.0190.7640.9360.10370.1050.348
OV (n = 303)NK cellKIR3DL2−0.0630.3180.6310.1495*0.107
OV (n = 303)NK cellKIR3DL3−0.0440.4660.7510.05710.3680.682
OV (n = 303)NeutrophilCCR7−0.0680.3240.6330.09430.1380.421
OV (n = 303)NeutrophilCEACAM8−0.0650.3440.658−0.0340.6190.839
OV (n = 303)NeutrophilITGAM−0.185***0.05750.3670.682
OV (n = 303)T cell (general)CD2−0.157*0.0650.1651***
OV (n = 303)T cell (general)CD3D−0.142*0.1060.1524*0.077
OV (n = 303)T cell (general)CD3E−0.126*0.1750.1581*0.051
OV (n = 303)TAMCCL2−0.171***0.1709***
OV (n = 303)TAMCD68−0.203***0.1050.0930.258
OV (n = 303)TAMIL100.0460.4320.707−0.0070.9480.979
OV (n = 303)TfhIL21−0.128*0.164−0.0160.8440.934
OV (n = 303)TfhBCL6−0.195***0.1285*0.158
OV (n = 303)Th1IFNG−0.0880.1860.4380.1323*0.146
OV (n = 303)Th1STAT1−0.0770.2290.4890.08940.1580.384
OV (n = 303)Th1STAT4−0.0090.8730.9590.07680.2250.468
OV (n = 303)Th1TBX21−0.159*0.0780.1587*0.063
OV (n = 303)Th1TNF−0.0380.5680.7780.020.7590.882
OV (n = 303)Th17IL17A−0.0730.2650.5280.00430.9430.981
OV (n = 303)Th17STAT3−0.0420.4880.7310.01170.8570.938
OV (n = 303)Th2GATA3−0.0610.3050.562−0.0840.2060.436
OV (n = 303)Th2IL13−0.0170.7610.8930.06470.3030.575
OV (n = 303)Th2STAT5A−0.1150.0620.215−0.0510.4230.681
OV (n = 303)Th2STAT6−0.0460.4580.7110.08690.1760.396
OV (n = 303)TregCCR8−0.0040.8930.9680.02110.7410.882
OV (n = 303)TregFOXP3−0.0590.4150.6750.06350.3170.588
OV (n = 303)TregSTAT5B0.152*0.081−0.1677***
OV (n = 303)TregTGFB1−0.1270.0520.1810.01530.8150.924

Correlation analysis between claudins and markers of immune cells in ovarian cancer via TIMER2.0.

DC, Dendritic cell; NK cel, Natural killer cell; TAM, Tumor-associated macrophage; Tfh, Follicular helper T cell; Treg, Regulatory T cell; *P < 0.05; **P < 0.01; ***P < 0.001.

Studies have shown that the tumor-infiltrating immune cells mentioned above are related to the tumor immunotherapy response (Rodriguez et al., 2018). Immune cell-based immunotherapy (Baci et al., 2020), including NK Cells (Nersesian et al., 2019) and dendritic cells (Stiff et al., 2013), play important roles in the treatment of ovarian cancer. Taken together, these analyses and our research indicate that CLDN6 and CLDN10 may play important roles in immunotherapy in the future.

Discussion

CLDN6 and CLDN10 are important components of the claudin family related to tight junctions. Claudins were considered promising targets for diagnosis and therapy since they were involved in uncontrolled cancer growth and metastasis (Martin and Jiang, 2001; Morin, 2005; Bose and Mukhopadhyay, 2010). Moreover, studies have shown that they not only play a vital role in tumorigenesis (Swisshelm et al., 2005; Arabzadeh et al., 2007; Hagen, 2019), but also drug resistance (Gao et al., 2017).

CLDN6 had been demonstrated abnormal expression and can be a prognostic marker in cancers including ovarian cancer (Wang et al., 2013), endometrial cancer (Kojima et al., 2020), gastric cancer (Kohmoto et al., 2020), breast carcinoma (Liu et al., 2016; Jia et al., 2019), and lung cancer (Micke et al., 2014). Bioinformatic analysis has revealed that CLDN6 is regulated by a diverse set of transcription factors and promotes cancer cell behavior via the ASK1-p38/JNK MAPK secretory signaling pathway (Lin et al., 2017). A study revealed that CLDN6 may be a novel targeted therapy for ovarian cancer as a receptor for clostridium perfringens enterotoxin (Lal-Nag et al., 2012). In addition, 6PHU3, a T-cell-engaging bispecific single chain antibody with anti-CD3/anti-CLDN6 specificities, upregulated the cytotoxicity of T cells and made T cells acquire an effector phenotype (Stadler et al., 2016). Another recent study showed that CLDN6 as a chimeric antigen receptor target in solid tumors can be a strategy to overcome inefficient CAR-T cell stimulation in vivo (Reinhard et al., 2020). These studies suggested that CLDN6 has important research value in the treatment of cancer.

CLDN10, a glandular epithelial marker in epithelial ovarian cancer (Seo et al., 2010), was reported to be a key immune-related gene in clear cell renal cell carcinoma (Yang et al., 2021) and papillary thyroid carcinoma (Xiang et al., 2020). Furthermore, CLDN10 expression has proved to be a prognostic marker for ovarian cancer (Li Z. et al., 2020).

The present study combined and analyzed the prognostic potential of CLDN6 and CLDN10 with the tumor immune microenvironment. Consistent with previous reports, both CLDN6 and CLDN10 showed high expression in ovarian cancer. Prognostic analysis showed that the overexpression of CLDN6 was related to a poor prognosis for patients with ovarian cancer. However, CLDN10 overexpression predicted a better prognosis compared to the low CLDN10 expression group. We also found that CLDN6 overexpression was negatively related to immune cell infiltration, whereas CLDN10 overexpression was positively correlated with immune cell infiltration. Moreover, we found that CLDN6 and CLDN10 were related to gene markers of dendritic cells, NK cells, and TAMs. These results may explain why the overexpression of CLDN6 and low expression of CLDN10 predict poor OS in ovarian cancer. This study revealed that the prognostic potential of CLDN6 and CLDN10 is related to the tumor immune microenvironment in ovarian cancer.

Relevant evidence has emerged that immune-related gene expression and TILs are related to the prognosis, recurrence (Ojalvo et al., 2018), and chemotherapeutic response (Choi et al., 2020) of ovarian cancer. Furthermore, the presence of TILs may improve clinical outcomes in ovarian cancer patients (Odunsi, 2017). Immune cell-based immunotherapy (Baci et al., 2020), including NK Cells (Nersesian et al., 2019) and dendritic cells (Stiff et al., 2013), play an important role in the treatment of ovarian cancer. Previous studies and our analyses suggest that CLDN6 may be involved in immune evasion and that they could be an ideal candidate for immunotherapy in ovarian cancer. Future studies on the combined application of claudin-based molecular targeted therapy and immunotherapy are necessary.

Conclusion

CLDN6 and CLDN10 were identified as potential prognostic biomarkers and were correlated with immune cell infiltration in ovarian cancer. Our results revealed new roles for CLDN6 and CLDN10 in ovarian cancer and their potential as therapeutic targets in cancer treatment.

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

Author contributions

PeW was responsible for the study conception and design. PG, TP, CC, and SL were involved in data acquisition, data analysis, and interpretation. PG drafted the manuscript and took charge of supervising the manuscript. All authors read and approved the manuscript.

Funding

This work was supported by the Natural Science Foundation of China (81772775 to JW and 82072895 to PeW).

Acknowledgments

This manuscript was released as a pre-print at Research Square (https://www.researchsquare.com/article/rs-40048/v1) (DOI: 10.21203/rs.3.rs-40048/v1) (Gao et al., 2020).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.595436/full#supplementary-material

Supplementary Figure 1

The significant unique analyses of claudins expression in ONCOMINE. (Red: overexpression; Blue: low expression).

Supplementary Figure 2

Associations between the expression of claudins and immune subtypes of ovarian cancer.

Supplementary Figure 3

Associations between the expression of claudins and molecular subtypes of ovarian cancer.

Supplementary Figure 4

The correlations between claudins (CLDN6 and CLDN10) expression and gene markers of immune cells across gynecologic oncology (Red: positive correlation; Blue: negative correlation).

Abbreviations

  • GEPIA

    Gene Expression Profiling Interactive Analysis

  • TIMER

    Tumor Immune Estimation Resource

  • GSEA

    gene set enrichment analyses

  • TCGA

    The Cancer Genome Atlas

  • GEO

    Gene Expression Omnibus

  • EGA

    European Genome-Phenome Archive

  • FPKM

    Fragments per kilobase per million

  • TPM

    transcripts per million

  • TILs

    tumor-infiltrating lymphocytes

  • NK

    natural killer cells

  • Tregs

    regulatory T cells

  • CAFs

    cancer associated fibroblasts

  • TAMs

    Tumor-associated macrophages

  • CPE

    clostridium perfringens enterotoxin

  • OS

    overall survival

  • PFS

    progression free survival

  • PPS

    post progression survival.

References

  • 1

    ArabzadehA.TroyT. C.TurksenK. (2007). Changes in the distribution pattern of Claudin tight junction proteins during the progression of mouse skin tumorigenesis.BMC Cancer7:196. 10.1186/1471-2407-7-196

  • 2

    BaciD.BosiA.GallazziM.RizziM.NoonanD. M.PoggiA.et al (2020). The Ovarian Cancer Tumor Immune Microenvironment (TIME) as Target for Therapy: a Focus on Innate Immunity Cells as Therapeutic Effectors.Int. J. Mol. Sci.21:3125. 10.3390/ijms21093125

  • 3

    Barros-FilhoM. C.MarchiF. A.PintoC. A.RogattoS. R.KowalskiL. P. (2015). High Diagnostic Accuracy Based on CLDN10, HMGA2, and LAMB3 Transcripts in Papillary Thyroid Carcinoma.J. Clin. Endocrinol. Metab.100E890–E899.

  • 4

    BoganiG.LopezS.MantieroM.DucceschiM.BosioS.RuisiS.et al (2020). Immunotherapy for platinum-resistant ovarian cancer.Gynecol. Oncol.158484–488.

  • 5

    BoseC. K.MukhopadhyayA. (2010). Claudin and ovarian cancer.J. Turk. Ger. Gynecol. Assoc.1148–54.

  • 6

    Cancer Genome Atlas Research Network. (2011). Integrated genomic analyses of ovarian carcinoma.Nature474609–615. 10.1038/nature10166

  • 7

    CeramiE.GaoJ.DogrusozU.GrossB. E.SumerS. O.AksoyB. A.et al (2012). The cBio Cancer Genomics Portal: an Open Platform for Exploring Multidimensional Cancer Genomics Data.Cancer Discov.2401–404. 10.1158/2159-8290.cd-12-0095

  • 8

    CharoentongP.FinotelloF.AngelovaM.MayerC.EfremovaM.RiederD.et al (2017). Pan-cancer Immunogenomic Analyses Reveal Genotype-Immunophenotype Relationships and Predictors of Response to Checkpoint Blockade.Cell Rep.18248–262. 10.1016/j.celrep.2016.12.019

  • 9

    ChoiK. U.KimA.KimJ. Y.KimK. H.HwangC.LeeS. J.et al (2020). Differences in immune-related gene expressions and tumor-infiltrating lymphocytes according to chemotherapeutic response in ovarian high-grade serous carcinoma.J. Ovarian Res.13:65.

  • 10

    GaoJ.AksoyB. A.DogrusozU.DresdnerG.GrossB.SumerS. O.et al (2013). Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal.Sci. Signal.6:l1.

  • 11

    GaoP.PengT.CaoC.LinS.WuP.HuangX.et al (2020). CLDN6 and CLDN10 are Associated with Immune Infiltration of Ovarian Cancer: a Study of Claudin Family.Front. Genet.10.21203/rs.3.rs-40048/v1[Preprint].

  • 12

    GaoY.LiuX.LiT.WeiL.YangA.LuY.et al (2017). Cross-validation of genes potentially associated with overall survival and drug resistance in ovarian cancer.Oncol. Rep.373084–3092. 10.3892/or.2017.5534

  • 13

    GyorffyB.LanczkyA.SzallasiZ. (2012). Implementing an online tool for genome-wide validation of survival-associated biomarkers in ovarian-cancer using microarray data from 1287 patients.Endocr. Relat. Cancer19197–208. 10.1530/erc-11-0329

  • 14

    HagenS. J. (2019). Unraveling a New Role for Claudins in Gastric Tumorigenesis.Cell. Mol. Gastroenterol. Hepatol.8151–152. 10.1016/j.jcmgh.2019.04.004

  • 15

    HamanishiJ.MandaiM.IwasakiM.OkazakiT.TanakaY.YamaguchiK.et al (2007). Programmed cell death 1 ligand 1 and tumor-infiltrating CD8+ T lymphocytes are prognostic factors of human ovarian cancer.Proc. Natl. Acad. Sci. U. S. A.1043360–3365. 10.1073/pnas.0611533104

  • 16

    HendersonJ. T.WebberE. M.SawayaG. F. (2018). Screening for Ovarian Cancer: updated Evidence Report and Systematic Review for the US Preventive Services Task Force.JAMA319595–606. 10.1001/jama.2017.21421

  • 17

    HewittK. J.AgarwalR.MorinP. J. (2006). The claudin gene family: expression in normal and neoplastic tissues.BMC Cancer6:186. 10.1186/1471-2407-6-186

  • 18

    JiaH.ChaiX.LiS.WuD.FanZ. (2019). Identification of claudin-2, -6, -11 and -14 as prognostic markers in human breast carcinoma.Int. J. Clin. Exp. Pathol.122195–2204.

  • 19

    KirschnerN.RosenthalR.FuruseM.MollI.FrommM.BrandnerJ. M.et al (2013). Contribution of tight junction proteins to ion, macromolecule, and water barrier in keratinocytes.J. Invest. Dermatol.1331161–1169. 10.1038/jid.2012.507

  • 20

    KohmotoT.MasudaK.ShodaK.TakahashiR.UjiroS.TangeS.et al (2020). Claudin-6 is a single prognostic marker and functions as a tumor-promoting gene in a subgroup of intestinal type gastric cancer.Gastric Cancer23403–417. 10.1007/s10120-019-01014-x

  • 21

    KojimaM.SugimotoK.TanakaM.EndoY.KatoH.HondaT.et al (2020). Prognostic Significance of Aberrant Claudin-6 Expression in Endometrial Cancer.Cancers12:2748. 10.3390/cancers12102748

  • 22

    Lal-NagM.BattisM.SantinA. D.MorinP. J. (2012). Claudin-6: a novel receptor for CPE-mediated cytotoxicity in ovarian cancer.Oncogenesis1:e33. 10.1038/oncsis.2012.32

  • 23

    LeungC. S.YeungT. L.YipK. P.WongK. K.HoS. Y.MangalaL. S.et al (2018). Cancer-associated fibroblasts regulate endothelial adhesion protein LPP to promote ovarian cancer chemoresistance.J. Clin. Invest.128589–606. 10.1172/jci95200

  • 24

    LiB.RuottiV.StewartR. M.ThomsonJ. A.DeweyC. N. (2010). RNA-Seq gene expression estimation with read mapping uncertainty.Bioinformatics26493–500. 10.1093/bioinformatics/btp692

  • 25

    LiB.SeversonE.PignonJ. C.ZhaoH.LiT.NovakJ.et al (2016). Comprehensive analyses of tumor immunity: implications for cancer immunotherapy.Genome Biol.17:174.

  • 26

    LiT.FanJ.WangB.TraughN.ChenQ.LiuJ. S.et al (2017). TIMER: a Web Server for Comprehensive Analysis of Tumor-Infiltrating Immune Cells.Cancer Res.77e108–e110.

  • 27

    LiT.FuJ.ZengZ.CohenD.LiJ.ChenQ.et al (2020). TIMER2.0 for analysis of tumor-infiltrating immune cells.Nucleic Acids Res.48W509–W514.

  • 28

    LiZ.XuanW.HuangL.ChenN.HouZ.LuB.et al (2020). Claudin 10 acts as a novel biomarker for the prognosis of patients with ovarian cancer.Oncol. Lett.20373–381.

  • 29

    LinD.GuoY.LiY.RuanY.ZhangM.JinX.et al (2017). Bioinformatic analysis reveals potential properties of human Claudin-6 regulation and functions.Oncol. Rep.38875–885. 10.3892/or.2017.5756

  • 30

    LiuY.JinX.LiY.RuanY.LuY.YangM.et al (2016). DNA methylation of claudin-6 promotes breast cancer cell migration and invasion by recruiting MeCP2 and deacetylating H3Ac and H4Ac.J. Exp. Clin. Cancer Res.35:120.

  • 31

    MartinT. A.JiangW. G. (2001). Tight junctions and their role in cancer metastasis.Histol. Histopathol.161183–1195.

  • 32

    Mhawech-FaucegliaP.WangD.SamraoD.KimG.LawrensonK.MenesesT.et al (2014). Clinical Implications of Marker Expression of Carcinoma-Associated Fibroblasts (CAFs) in Patients with Epithelial Ovarian Carcinoma After Treatment with Neoadjuvant Chemotherapy.Cancer Microenviron733–39. 10.1007/s12307-013-0140-4

  • 33

    MickeP.MattssonJ. S.EdlundK.LohrM.JirströmK.BerglundA.et al (2014). Aberrantly activated claudin 6 and 18.2 as potential therapy targets in non-small-cell lung cancer.Int. J. Cancer1352206–2214. 10.1002/ijc.28857

  • 34

    MorinP. J. (2005). Claudin proteins in human cancer: promising new targets for diagnosis and therapy.Cancer Res.659603–9606. 10.1158/0008-5472.can-05-2782

  • 35

    NersesianS.GlazebrookH.ToulanyJ.GranthamS. R.BoudreauJ. E. (2019). Naturally Killing the Silent Killer: NK Cell-Based Immunotherapy for Ovarian Cancer.Front. Immunol.10:1782. 10.3389/fimmu.2019.01782

  • 36

    NewmanA. M.SteenC. B.LiuC. L.GentlesA. J.ChaudhuriA. A.SchererF.et al (2019). Determining cell type abundance and expression from bulk tissues with digital cytometry.Nat. Biotechnol.37773–782. 10.1038/s41587-019-0114-2

  • 37

    OdunsiK. (2017). Immunotherapy in ovarian cancer.Ann. Oncol.28viii1–viii7.

  • 38

    OjalvoL. S.ThompsonE. D.WangT. L.MeekerA. K.ShihI. M.FaderA. N.et al (2018). Tumor-associated macrophages and the tumor immune microenvironment of primary and recurrent epithelial ovarian cancer.Hum. Pathol.74135–147. 10.1016/j.humpath.2017.12.010

  • 39

    Ovarian Tumor Tissue AnalysisC.GoodeE. L.BlockM. S.KalliK. R.VierkantR. A.ChenW.et al (2017). Dose-Response Association of CD8+ Tumor-Infiltrating Lymphocytes and Survival Time in High-Grade Serous Ovarian Cancer.JAMA Oncol.3:e173290.

  • 40

    ReinhardK.RengstlB.OehmP.MichelK.BillmeierA.HaydukN.et al (2020). An RNA vaccine drives expansion and efficacy of claudin-CAR-T cells against solid tumors.Science367446–453. 10.1126/science.aay5967

  • 41

    RhodesD. R.Kalyana-SundaramS.MahavisnoV.VaramballyR.YuJ.BriggsB. B.et al (2007). Oncomine 3.0: genes, pathways, and networks in a collection of 18,000 cancer gene expression profiles.Neoplasia9166–180. 10.1593/neo.07112

  • 42

    RodriguezG. M.GalpinK. J. C.MccloskeyC. W.VanderhydenB. C. (2018). The Tumor Microenvironment of Epithelial Ovarian Cancer and Its Influence on Response to Immunotherapy.Cancers10:242. 10.3390/cancers10080242

  • 43

    RuB.WongC. N.TongY.ZhongJ. Y.ZhongS. S. W.WuW. C.et al (2019). TISIDB: an integrated repository portal for tumor-immune system interactions.Bioinformatics354200–4202. 10.1093/bioinformatics/btz210

  • 44

    SatoE.OlsonS. H.AhnJ.BundyB.NishikawaH.QianF.et al (2005). Intraepithelial CD8+ tumor-infiltrating lymphocytes and a high CD8+/regulatory T cell ratio are associated with favorable prognosis in ovarian cancer.Proc. Natl. Acad. Sci. U. S. A.10218538–18543. 10.1073/pnas.0509182102

  • 45

    SeoH. W.RengarajD.ChoiJ. W.AhnS. E.SongY. S.SongG.et al (2010). Claudin 10 is a glandular epithelial marker in the chicken model as human epithelial ovarian cancer.Int. J. Gynecol. Cancer201465–1473.

  • 46

    SiegelR. L.MillerK. D.JemalA. (2020). Cancer statistics, 2020.CA Cancer J. Clin.707–30.

  • 47

    StadlerC. R.Bähr-MahmudH.PlumL. M.SchmoldtK.KölschA. C.TüreciÖ, et al. (2016). Characterization of the first-in-class T-cell-engaging bispecific single-chain antibody for targeted immunotherapy of solid tumors expressing the oncofetal protein claudin 6.Oncoimmunology5:e1091555. 10.1080/2162402x.2015.1091555

  • 48

    StiffP. J.CzerlanisC.DrakesM. L. (2013). Dendritic cell immunotherapy in ovarian cancer.Expert Rev. Anticancer Ther.1343–53.

  • 49

    SturmG.FinotelloF.PetitprezF.ZhangJ. D.BaumbachJ.FridmanW. H.et al (2019). Comprehensive evaluation of transcriptome-based cell-type quantification methods for immuno-oncology.Bioinformatics35i436–i445.

  • 50

    SwisshelmK.MacekR.KubbiesM. (2005). Role of claudins in tumorigenesis.Adv. Drug Deliv. Rev.57919–928. 10.1016/j.addr.2005.01.006

  • 51

    TangZ.LiC.KangB.GaoG.LiC.ZhangZ.et al (2017). GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses.Nucleic Acids Res.45W98–W102.

  • 52

    ThorssonV.GibbsD. L.BrownS. D.WolfD.BortoneD. S.Ou YangT. H.et al (2018). The Immune Landscape of Cancer.Immunity48812–830.e14.

  • 53

    TsukitaS.FuruseM.ItohM. (2001). Multifunctional strands in tight junctions.Nat. Rev. Mol. Cell Biol.2285–293. 10.1038/35067088

  • 54

    WangL.JinX.LinD.LiuZ.ZhangX.LuY.et al (2013). Clinicopathologic significance of claudin-6, occludin, and matrix metalloproteinases -2 expression in ovarian carcinoma.Diagn. Pathol.8:190.

  • 55

    WeinsteinR. S.MerkF. B.AlroyJ. (1976). The structure and function of intercellular junctions in cancer.Adv. Cancer Res.2323–89. 10.1016/s0065-230x(08)60543-6

  • 56

    WodarzA. (2000). Tumor suppressors: linking cell polarity and growth control.Curr. Biol.10R624–R626.

  • 57

    XiangZ.ZhongC.ChangA.LingJ.ZhaoH.ZhouW.et al (2020). Immune-related key gene CLDN10 correlates with lymph node metastasis but predicts favorable prognosis in papillary thyroid carcinoma.Aging122825–2839. 10.18632/aging.102780

  • 58

    YangW.LiL.ZhangK.MaK.GongY.ZhouJ.et al (2021). CLDN10 associated with immune infiltration is a novel prognostic biomarker for clear cell renal cell carcinoma.Epigenomics1331–45. 10.2217/epi-2020-0256

  • 59

    YoshiharaK.ShahmoradgoliM.MartinezE.VegesnaR.KimH.Torres-GarciaW.et al (2013). Inferring tumour purity and stromal and immune cell admixture from expression data.Nat. Commun.4:2612.

  • 60

    ZhangZ.WangA.SunB.ZhanZ.ChenK.WangC.et al (2013). Expression of CLDN1 and CLDN10 in lung adenocarcinoma in situ and invasive lepidic predominant adenocarcinoma.J. Cardiothorac. Surg8:95.

  • 61

    ZhouY.XiangJ.BhandariA.GuanY.XiaE.ZhouX.et al (2018). CLDN10 is Associated with Papillary Thyroid Cancer Progression.J Cancer94712–4717. 10.7150/jca.28636

Summary

Keywords

ovarian cancer, CLDN6, CLDN10, prognosis, immune microenvironment

Citation

Gao P, Peng T, Cao C, Lin S, Wu P, Huang X, Wei J, Xi L, Yang Q and Wu P (2021) Association of CLDN6 and CLDN10 With Immune Microenvironment in Ovarian Cancer: A Study of the Claudin Family. Front. Genet. 12:595436. doi: 10.3389/fgene.2021.595436

Received

16 August 2020

Accepted

05 May 2021

Published

23 June 2021

Volume

12 - 2021

Edited by

Bairong Shen, Sichuan University, China

Reviewed by

Rajeev K. Singla, Sichuan University, China; Zhongqiu Xie, University of Virginia, United States

Updates

Copyright

*Correspondence: Peng Wu, Qin Yang,

This article was submitted to Systems Biology Archive, a section of the journal Frontiers in Genetics

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics