- 1Department of Pharmacology, Hebei University of Chinese Medicine, Shijiazhuang, China
- 2Department of Hematology, The Fourth Hospital of Hebei Medical University, Shijiazhuang, China
- 3College of Natural Science, Michigan State University, East Lansing, MI, United States
- 4Shijiazhuang People’s Hospital, Shijiazhuang, China
Background: Claudins (CLDNs) are a family of closely related transmembrane proteins that have been linked to oncogenic transformation and metastasis across a range of cancers, suggesting that they may be valuable diagnostic and/or prognostic biomarkers that can be used to evaluate patient outcomes. However, CLDN expression patterns associated with colorectal cancer (CRC) remain to be defined.
Methods: The mRNA levels of 21 different CLDN family genes were assessed across 20 tumor types using the Oncomine database. Correlations between these genes and patient clinical outcomes, immune cell infiltration, clinicopathological staging, lymph node metastasis, and mutational status were analyzed using the GEPIA, UALCAN, Human Protein Atlas, Tumor Immune Estimation Resource, STRING, Genenetwork, cBioportal, and DAVID databases in an effort to clarify the potential functional roles of different CLDN protein in CRC. Molecular docking analyses were used to probe potential interactions between CLDN4 and TGFβ1. Levels of CLDN4 and CLDN11 mRNA expression in clinical CRC patient samples and in the HT29 and HCT116 cell lines were assessed via qPCR. CLDN4 expression levels in these 2 cell lines were additionally assessed following TGFβ1 inhibitor treatment.
Results: These analyses revealed that COAD and READ tissues exhibited the upregulation of CLDN1, CLDN2, CLDN3, CLDN4, CLDN7, and CLDN12 as well as the downregulation of CLDN5 and CLDN11 relative to control tissues. Higher CLDN11 and CLDN14 expression as well as lower CLDN23 mRNA levels were associated with poorer overall survival (OS) outcomes. Moreover, CLDN2 and CLDN3 or CLDN11 mRNA levels were significantly associated with lymph node metastatic progression in COAD or READ lower in COAD and READ tissues. A positive correlation between the expression of CLDN11 and predicted macrophage, dendritic cell, and CD4+ T cell infiltration was identified in CRC, with CLDN12 expression further being positively correlated with CD4+ T cell infiltration whereas a negative correlation was observed between such infiltration and the expression of CLDN3 and CLDN15. A positive correlation between CLDN1, CLDN16, and neutrophil infiltration was additionally detected, whereas neutrophil levels were negatively correlated with the expression of CLDN3 and CLDN15. Molecular docking suggested that CLDN4 was able to directly bind via hydrogen bond with TGFβ1. Relative to paracancerous tissues, clinical CRC tumor tissue samples exhibited CLDN4 and CLDN11 upregulation and downregulation, respectively. LY364947 was able to suppress the expression of CLDN4 in both the HT29 and HCT116 cell lines.
Conclusion: Together, these results suggest that the expression of different CLDN family genes is closely associated with CRC tumor clinicopathological staging and immune cell infiltration. Moreover, CLDN4 expression is closely associated with TGFβ1 in CRC, suggesting that it and other CLDN family members may represent viable targets for antitumor therapeutic intervention.
Introduction
Colorectal cancer (CRC) is the third most common and fourth deadliest cancer (Mármol et al., 2017). As CRC is often diagnosed at an early stage, it generally has a good 5-year overall survival (OS) rate, but patients diagnosed at a later stage with metastatic disease often fair poorly (Yamagishi et al., 2016). Diagnostic and prognostic biomarkers of CRC thus offer clear clinical value for identifying and monitoring affected patients.
Claudins (CLDNs) are members of a 27 + gene family of transmembrane proteins with four transmembrane helical domains, two extracellular loops, and short N- and C-terminal domains (Suzuki et al., 2014). These CLDNs play important roles in tumorigenesis and can influence aggressive growth and motility owing to their role as regulators of intercellular adhesion. Indeed, there is growing evidence that CLDN dysregulation is common across many cancer types such as gastric, lung, breast, ovarian, and colorectal cancer (Yu et al., 2009; Tabariès and Siegel, 2017).
Herein, we surveyed patterns of CLDN mRNA and protein expression in different CRC patient subsets in order to elucidate the links between these different genes and outcomes associated with different CRC stages. In addition, relationships between CLDN expression, mutational status, and immune cell infiltration in CRC tumors were assessed to better clarify the mechanistic role of these CLDNs and to guide the selection of future targets for therapeutic intervention when treating patients with this cancer type.
Materials and Methods
Clinical Samples
CRC patient tumor and paracancerous tissues were collected from the First Department of General Surgery, Shijiazhuang People’s Hospital (Hebei Province, China), and Hebei Medical University Fourth Hospital (China). Pathologists confirmed the diagnosis and staging of all patients, and the Ethical Committee of Shijiazhuang People’s Hospital approved this study. All patients provided written informed consent to participate.
Cell Culture
Human HT29 and HCT116 cell lines were grown in McCoy’s 5A media (Gibco, CA, USA) containing 10% FBS and penicillin/streptomycin at 37°C in a 5% CO₂ incubator (Yang et al., 2021).
Oncomine Database Analysis of CLDN mRNA Profiles
The Oncomine database (http://www.oncomine.com) (Rhodes et al., 2004) was utilized to evaluate the expression of 21 different CLDN family members in 20 cancer types, using the following criteria for differential expression: p = 0.01, Fold-change > 1.5, gene rank ≤10%.
GEPIA Analyses
The GEPIA database (http://gepia.cancer-pku.cn/) (Tang et al., 2017) was used to analyze colon adenocarcinoma (COAD) and rectal adenocarcinoma (READ) tissue samples in order to explore the relationships between different CLDN family members and key clinical outcomes including staging, overall survival (OS), and disease-free survival (DFS). Default GEPIA parameters were utilized for these analyses.
Human Protein Atlas Analyses
IHC staining data pertaining to the protein level expression of different CLDN family members were assessed with the Human Protein Atlas database (https://www.proteinatlas.org/) in CRC patient tumors and normal tissue samples (Pontén et al., 2008).
Immune Cell Infiltration Analysis
The association between different CLDN family members and the infiltration of immune cells into CRC tumors was assessed with the TIMER database (https://cistrome.shinyapps.io/timer/), with scatter plots for different CLDN genes being generated to demonstrate purity-corrected partial Spearman’s r values and corresponding significance metrics (Li et al., 2021).
UALCAN Database Analyses
Staging and nodal metastasis status for COAD (n = 324) and READ (n = 172) patient clinicopathologic parameters were evaluated using the UALCAN database (Lyu et al., 2020). According to multivariate Kaplan-Meier survival analysis, with p-values calculated using the log-rank test (log-rank test). The Ualcan database uses TCGA RNA-seq and clinical data for 31 cancer types.
Analysis of CLDN Mutational Status in CRC
The mutational status of different CLDN family genes in CRC was assessed using the cBioPortal database (http://cbioportal.org) (Cerami et al., 2012). Assessed mutations included deep deletions, missense mutations, copy number amplifications, and mRNA upregulation.
CLDN Matrix Interaction Analyses
An interaction matrix for CLDN family genes was generated using the “matrix” command in the GeneNetwork database (Mulligan et al., 2017), enabling a correlation analysis between different CLDN family members.
Protein-Protein Interaction Network
The STRING database (http://string-db.org; version 11.0) was used to construct a CLDN family gene PPI network incorporating 24 co-expressed genes with a score of >0.4 (Franceschini et al., 2013). The network was visualized using Cytoscape (v 3.8.2).
GO and KEGG Enrichment Analysis
GO analyses were used to assess the enrichment of particular genes in specific functional categories including biological processes (BPs), cellular components (CCs), and molecular functions (MFs), enabling efficient analyses of transcriptomic datasets (Thomas et al., 2019). The Kyoto Encyclopedia of Genes and Genomes (KEGG) database compiles functional data pertaining to a diverse array of regulatory pathways (Kanehisa et al., 2016). The online DAVID database (https://david.ncifcrf.gov/; version 6.8) was used for GO annotation and KEGG pathway analyses in this study (Huang et al., 2007). R was used to visualize the enrichment data.
Molecular Docking Analysis
Full-length wild-type protein sequences were obtained from UniProt (https://www.uniprot.org/) (UniProt Consortium, 2021), with Uniprot ID. O14493 and Uniprot ID. P01137 being used for CLDN4 and TGFβ1, respectively. Three-dimensional structures for these proteins were obtained from the RCSB PDB database (https://www.rcsb.org/) (Berman et al., 2000), using PDB IDs of 7KP4 and 5VQP, respectively, for docking analyses. The CLDN4 and TGFβ1 proteins were docked with the Zdock server 3.0.2 (https://zdock.umassmed.edu/) (Pierce et al., 2014). Prior to docking, PyMOL was utilized to remove water molecules, heteroatoms, and repeated subunits, with the best docked complex being assessed for hydrogen bonding.
qRT-PCR
The ISOGEN reagent (Nippon Gene Co. Ltd., Kokyo, Japan) was utilized to isolate RNA from tissue and cell line samples, after which qPCR was conducted by initially reverse transcribing 1 µg of total RNA per sample using a Revert Aid first strand cDNA synthesis kit (ThermoFisher Scientific). Primers used for qPCR are compiled in Supplementary Table S1. A Real-Time PCR system (BIOER Co. Ltd., Kokyo, Japan) was used for all qPCR analyses, and the relative expression of CLDN11, CLDN4, and TGFβ1 was assessed via the ΔCT method as in prior reports (Yang et al., 2021). For appropriate experiments, the HT29 and HCT116 cells were treated for 48 h with the TGFβ1 inhibitor LY364947 (5.00 or 10.0uM).
Western Blot
HT29 and HCT116 cells were treated for 48 h with LY364947 at concentrations of 5, and 10 μM. After being washed three times with PBS, the cells were lysed for 30 min on ice before being centrifuged at 10000 g for 5 min at 4°C for 5 min. The protein concentrations in cell lysates were determined using an ND-1000 Spectrophotometer (National Instruments) (NanoDrop, ThermoFisher Scientific). SDS-PAGE was used to separate out equal amounts of total protein, which was then transferred to PVDF membranes (Millipore, Billerica, MA, USA). After blocking with 5% milk for 2 h, we did overnight immunoblotting with primary antibodies at 4°C. The primary antibodies were raised against TGFβ1 (1:1,000; Abcam, Cambridge, MA, USA), CLDN4 (1:500; abways technology) and β-actin (AC026, at 1 : 10000 dilution). A secondary antibody of goat anti-rabbit IG (KPL074-1,506, at 1: 5,000 dilution) was then used to incubate the membranes for 1 h at 37°C. Fusion FX5 Spectra was used to determine the intensity of protein bands (Fusion, France). The membranes were then incubated for 1 h at 37°C with secondary antibodies of goat anti-rabbit IG (KPL074-1,506, at 1:5,000 dilution). The intensity of protein bands was estimated via executing Fusion FX5 Spectra (Fusion, France).
Results
The Association Between CLDN Family Gene Expression and CRC Pathological Type
We began by using the Oncomine database to evaluate CLDN family member expression in CRC. In total, 19 CLDNs were found to be differentially expressed between normal tissues and COAD and READ patient samples. At the mRNA level, CLDN1, CLDN2, CLDN12, and CLDN14 were upregulated whereas CLDN3, CLDN5, CLDN7, CLDN8, CLDN11, CLDN15, and CLDN23 were downregulated (FC > 1.5) in patients with CRC (Figure 1). Specifically, fold-change values for CLDN1, CLDN2, and CLDN12 in COAD tissues were 6.755, 3.006, and 2.096, respectively (Supplementary Table S1), while in READ tissues these respective fold-change values were 4.301, 1.849, and 1.583 in the dataset analysis of Kaiser et al. (2007). CLDN19 expression was altered by 1.117 and 1.158-fold in COAD and READ tissues from this same dataset, while CLDN20 expression was changed by 1.075- and 1.091-fold, respectively. Gaedcke et al. (2010) found fold-change values for the expression of CLDN1, CLDN10, CLDN11, CLDN15, and CLDN22 of 19.563, 1.583, 1.718, 1.417, and 1.05, respectively, in READ tissues. Skrzypczak et al. (2010) further reported CLDN1 and CLDN2 expression levels in CRC tissues that were altered by 6.351- and 7.754-fold, respectively. In the TCGA dataset, CLDN14 expression in READ and COAD was altered by 4.876-fold and 4.368-fold, respectively, while it was altered by 3.471-fold in the dataset generated by Gaedke et al. CLDN18 expression was altered by 1.298-fold and 1.184-fold in COAD data from TCGA and Kaiser et al., and by 1.216-fold and 1.25-fold in READ samples from these two respective data sources.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/783016/fgene-13-783016-HTML/image_m/fgene-13-783016-g001.jpg)
FIGURE 1. CLDN family member expression levels across human cancers. Cut-off values for these mRNA analyses were as follows: p-value = 0.01, Fold-change = 1.5, gene rank = 10%. Red and blue correspond to overexpression and underexpression, respectively, with the strength of the color be proportional to the expression level for that gene.
The GEPIA database was next utilized to compare the expression of different CLDN family members in CRC and in normal colon tissue samples, revealing CLDN1, CLDN2, CLDN3, CLDN4, CLDN7, and CLDN12 to have been upregulated and CLDN5 and CLDN11 to have been downregulated in COAD and READ tissues relative to control samples (Figure 2). Significant differences in CLDN3, CLDN4, CLDN5, CLDN6, CLDN9, CLDN11, and CLDN12 expression were evident as a function of tumor stage (Figure 3). CLDN3, CLDN4 and CLDN6 were related with clinicopathological stage in CRC. The more CLDN5, CLDN9, and CLDN11 are expressed in CRC, the worse the stage. CLDN12 expression in CRC associated with pathological stage. We additionally examined the link between CLDN family gene mRNA levels and CRC patient lymph node metastasis, revealing a significant relationship between the expression of CLDN1, CLDN2, CLDN3, CLDN7, CLDN8, CLDN9, CLDN11, CLDN14, CLDN16 and CLDN23 nodal metastasis in COAD and READ (Figures 4A, B, C, G, H, I, K, M, N, R) (Figures 5A, B, C, E, F, G, I, K, M, O).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/783016/fgene-13-783016-HTML/image_m/fgene-13-783016-g002.jpg)
FIGURE 2. CLDN family gene expression in COAD and READ tissues. Tumors and normal tissue controls are respectively shown in red and grey (A-I). *p < 0.05.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/783016/fgene-13-783016-HTML/image_m/fgene-13-783016-g004.jpg)
FIGURE 4. The UALCAN database shows a link between CLDN gene expression and nodal metastases in COAD patients. Box plots represent the CLDN mRNA expression levels in normal tissues or in COAD patients with N0–N2 disease (A–R). *p < 0.05; **p < 0.01; ***p < 0.001.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/783016/fgene-13-783016-HTML/image_m/fgene-13-783016-g005.jpg)
FIGURE 5. The UALCAN database shows a link between CLDN gene expression and nodal metastases in READ patients. Box plots represent the CLDN mRNA expression levels in normal tissues or in READ patients with N0–N2 disease (A–O). *p < 0.05; **p < 0.01; ***p < 0.001.
The Association Between CLDN Expression and CRC Patient Survival Outcomes
Kaplan-Meier curves and log-rank tests were next performed, revealing a relationship between the upregulation of CLDN11, CLDN14, and CLDN23 at the mRNA level and COAD and READ patient OS (Figure 6). Specifically, high levels of CLDN11 and CLDN14 expression were correlated with worse patient OS, whereas CLDN23 overexpression was linked to better OS outcomes.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/783016/fgene-13-783016-HTML/image_m/fgene-13-783016-g006.jpg)
FIGURE 6. Kaplan–Meier survival analyses assing the link between CLDN expression and CRC patient survival in GEPIA database (A-C) (p < 0.05).
Analysis of CLDN Protein Levels in CRC
The Human Protein Atlas yielded findings that were somewhat consistent with the mRNA level data above, with higher CLDN1 and CLDN12 protein levels and lower CLDN11 levels in CRC tumors relative to samples of normal colon tissue (Figure 7).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/783016/fgene-13-783016-HTML/image_m/fgene-13-783016-g007.jpg)
FIGURE 7. IHC staining for CLDN protein expression levels in tumor samples from CRC patients and normal colon tissue samples (A-G).
The Relationship Between CLDN Expression and CRC Tumor Immune Cell Infiltration
Using the TIMER database, we detected a negative correlation between the mRNA level expression of CLDN5, CLDN8, CLDN11, and CLDN18 and CRC tumor purity (Figures 8E, H, J, O). Moreover, CLDN12 expression was correlated with B cell infiltration and positively correlated with CD8+ T cells (Figures 8K). Additionally, CLDN5, CLDN9, CLDN11, CLDN15, and CLDN20 were significantly positively correlated with CD4+ T cell infiltration into CRC tumors, while CLDN7 expression was negatively correlated with such infiltration (Figures 8E, H, J, M, Q). Additionally, a clear positive correlation was observed between macrophage infiltration and the expression of CLDN1, CLDN9, CLDN11, and CLDN16 (Figures 8A,H,J,N). Neutrophil infiltration was negatively correlated with CLDN15 expression and positively correlated with CLDN1 and CLDN16 expression (Figures 8A, M, N), while CLDN11 expression was positively correlated with dendritic cell infiltration (Figures 8J).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/783016/fgene-13-783016-HTML/image_m/fgene-13-783016-g008.jpg)
FIGURE 8. The relationship between CLDN mRNA level expression and CRC tumor immune cell infiltration (A-S).
CLDN Mutation and the Relationships in CLDN Family in COAD Patients
Next, we assessed CLDN family gene mutations in 640 COAD samples using the online cBioPortal tool, revealing mutations in 283 samples (44%) (Figures 9A). Mutations in CLDN10 (6%) were attributable to missense mutations and gene amplification.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/783016/fgene-13-783016-HTML/image_m/fgene-13-783016-g009.jpg)
FIGURE 9. (A) Cancer type summary and Oncoprint for CLDN family genes in BioPortal. (B) Heatmap demonstrating correlations among CLDN family genes with respect to mRNA expression. Colors correspond to Pearson correlation coefficients, which span from blue to red (negative and positively correlated, respectively). (C) A PPI network of CLDNs was generated using the STRING database to clarify interactions among these proteins.
The Genenetwork database was also used to explore relationships among CLDN gene expression levels in COAD samples, revealing the following correlative relationships: CLDN1 with CLDN2; CLDN3 with CLDN4, CLDN7, CLDN8, CLDN14, and CLDN23; CLDN4 with CLDN7, CLDN8, and CLDN14; CLDN5 with CLDN10, CLDN11, CLDN18, and CLDN24; CLDN7 with CLDN8 and CLDN14; CLDN8 with CLDN14; CLDN9 with CLDN3, CLDN4, CLDN7, CLDN8, CLDN14, CLDN15, and CLDN23; CLDN10 with CLDN18 and CLDN24; CLDN11 with CLDN10, CLDN18, and CLDN24; CLDN15 with CLDN3, CLDN4, CLDN7, CLDN8, CLDN14, and CLDN23; CLDN18 with CLDN24; and CLDN23 with CLDN4, CLDN7, CLDN8, CLDN14, and CLDN23. A co-expression network for these CLDNs is shown in Figures 9B.
A PPI network for these CLDN family genes was generated composed of 24 nodes and 254 edges (Figures 9C).
GO Function and KEGG Pathway Enrichment Analyses of CLDN Genes in CRC
The DAVID database was utilized to predict the functional roles of different CLDNs in key CRC-related biological processes, molecular functions, and cellular components (Figures 10A). KEGG pathway analyses further revealed a close relationship between these CLDN family members and tight junctions, cell adhesion molecules, leukocyte transendothelial migration, and Hepatitis C (Supplementary Table S2).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/783016/fgene-13-783016-HTML/image_m/fgene-13-783016-g010.jpg)
FIGURE 10. (A) DAVID was used for the functional enrichment analysis of CLDN family genes, with results shown using a bubble chart.
Molecular Docking Analyses
A prior study of inflammatory bowel disease demonstrated the ability of TGFβ1 signaling inhibition to suppress CLDN4 expression (Marincola Smith et al., 2021). A molecular docking analysis was therefore conducted to explore putative interactions between CLDN4 and TGFβ1. For this approach, ZDOCK scores were used to estimate protein-protein binding affinity, with higher scores being indicating of greater binding affinity (Pierce et al., 2011). This approach revealed that nine hydrogen bonds were predicted to form between CLDN4 and TGFβ1, with interactions between Glu109 of CLDN4 and Gln177 of TGFβ1, Cys107 of CLDN4 and Gly173 and Arg181 of TGFβ1, Val95 of CLDN4 rand Ser244 of TGFβ1, Ser98 of CLDN4 and His247 and Arg249 of TGFβ1, Trp18 of CLDN4 and Leu242 of TGFβ1, Ala20 of CLDN4 and Glu261 of TGFβ1, and Leu173 of CLDN4 and Thr260 of TGFβ1 (Figures 11A).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/783016/fgene-13-783016-HTML/image_m/fgene-13-783016-g011.jpg)
FIGURE 11. (A) A docking analysis of CLDN4 and TGFβ1. (B) (C) The protein of CLDN4 levels was increased in human CRC tissues (n = 5, *p < 0.05; **p < 0.01). (D) CLDN4 and CLDN11 mRNA levels were respectively increased and decreased in human CRC tissues (n = 10, *p < 0.05; **p < 0.01). (E) (F) (G) Western blot and qRT-PCR showed TGFβ1 and CLDN4 were remarkably downregulated at both protein and mRNA levels after treatment with LY364947 at 5.0 and 10.0 uM/ml, respectively, on HT29 cells. (H) (I) (J) Western blot and qRT-PCR showed TGFβ1 and CLDN4 were remarkably downregulated at both protein and mRNA levels after treatment with LY364947 at 5.0 and 10.0 uM/ml, respectively, on HCT116 cells. TGFβ1 reduced CLDN4 expression in colorectal cancer cell lines.
Assessment of CLDN4 and CLDN11 Expression and the Responsiveness of CLDN4 to TGFβ1 Inhibitor Treatment
At the mRNA level, we found that CLDN11 and CLDN4 expression were respectively reduced and increased in CRC tumor samples relative to paracancerous controls as measured via qPCR (n = 10, p < 0.05). At the protein level, CLDN4 expression was increased by WB (n = 5, p < 0.05) (Figures 11B,C,D) (Supplementary Table S3). Treatment with LY364947 (TGFβ1 inhibitor) at a dose of five or 10 uM was sufficient to suppress CLDN4 mRNA and protein expression at 48 h post-treatment in both HT29 and HCT116 cells (Figure 11E, F, G, H, I, J).
Discussion
Claudin (CLDN) genes encode a family of proteins that play a critical role in the generation and functioning of tight junctions between epithelial and endothelial cells (Hewitt et al., 2006). Many malignancies rely on this protein family for their development, progression, and metastasis. Although some studies have demonstrated the importance of claudins in tumorigenesis and the prognosis of human cancers (Arabzadeh et al., 2007; Morel et al., 2012), no extensive bioinformatic analyses of this gene family have been performed in colorectal cancer to date. In this study, each gene in the claudin family was evaluated to examine its expression and prognostic relevance in colorectal cancer using bioinformatics tools, offering an opportunity to better understand how claudin gene dysregulation affects colorectal cancer, as well as how treatment can be optimized to improve patient prognosis.
Our findings revealed that the mRNA expression levels of CLDN1, CLDN2, CLDN12, and CLDN14 were upregulated in colorectal cancer tissues relative to normal tissues in the Oncomine database, whereas CLDN3, CLDN5, CLDN7, CLDN8, CLDN11, CLDN15, and CLDN23 were downregulated, as previously reported. However, when comparing colorectal cancer tissues from the GEPIA database to normal colonic mucosa, the mRNA level expression of CLDN3, CLDN7, and CLDN23 was shown to be significantly higher in these tumor tissues. Given the inconsistencies in claudin gene expression in these two databases, we focused on those claudins exhibiting comparable gene expression patterns in both databases. In addition, three genes were shown to be related with advanced colorectal cancer: CLDN11 (p = 0.007), and CLDN14 (p = 0.03) and CLDN23 (p = 0.021). Moreover, we found that CLDN11 and CLDN18 protein levels were decreased in CRC tumors relative to normal tissues. At the same time, high levels of CLDN11 expression were associated with an increase in the proportion of CD4+ T cells, macrophages, neutrophils, and dendritic cells within tumors. We further determined that CLDN4 and TGFβ1 were capable of forming hydrogen bonds, and that CLDN4 was overexpressed in CRC cell lines and inhibited by TGFβ inhibitor treatment.
Claudin expression is related to tumor proliferation, migration, and invasion (Mori et al., 2011; Pope et al., 2014; Wang and Zöller, 2019; Li et al., 2020a). Many claudins are dysregulated in a range of cancers, as in the case of the CLDN1 and CLDN2 genes, which are overexpressed in colorectal cancer (Dhawan et al., 2011; Mori et al., 2011; Pope et al., 2014; Wang et al., 2019). CLDN3 is similarly overexpressed at the mRNA level in CRC tumors that belong to the consensus molecular subgroup (CMS)-CMS2 and CMS3, which correspond to a poor prognosis (Tang et al., 2011; Perez et al., 2020). Colorectal cancers have significant levels of CLDN4 expression (Georges et al., 2012). It has been shown that combination treatment with 5-fluorouracil (FU) and an anti-CLDN4 extracellular domain antibody (4D3) enhances antitumor efficacy against CRC (Fujiwara-Tani et al., 2018). CLDN5 expression has been shown to be reduced in squamous cell carcinomas of the lung (Akizuki et al., 2017), cervical cancer (Zhu et al., 2015), hepatocellular carcinoma (Sakaguchi et al., 2008), oral squamous cell carcinoma (Phattarataratip and Sappayatosok, 2016), breast cancer (Escudero-Esparza et al., 2012), and pancreatic cancer (Jakab et al., 2011). Similarly, compared to normal mucosal tissues, CLDN5 expression was downregulated in CRC tumors (Cherradi et al., 2019), although the mechanistic basis for this change remains poorly understood. CLDN6 is a cellular adhesion protein that is abundantly expressed in gastric tumors tissues and cell lines, and is associated with a poor prognosis. CLDN6 enhances gastric cancer cell proliferation and invasiveness via the YAP1 and YAP1-Snail1 axis (Li et al., 2020a). High expression of CLDN7 has been shown to accelerate pancreatic and colon cancer progression (Mu et al., 2019; Wang and Zöller, 2019). Cancer initiation cells (CICs) and nasopharyngeal carcinoma metastases in gastrointestinal tumors can be detected using CLDN7 as a marker protein (Liu et al., 2017; Kyuno et al., 2019). The MAPK/ERK signaling pathway is activated in colorectal cancers when CLDN8 is overexpressed (Cheng et al., 2019). In contrast, low levels of CLDN8 expression slow the progression of renal clear cell carcinoma via the EMT/AKT signaling pathway (Zhu et al., 2020). Colonic epithelial barrier function in mice is influenced by TGF family signaling, which regulates CLDN2, CLDN4, CLDN7, and CLDN8 in colonic epithelial cells. We additionally confirmed that TGFβ1 and CLDN4 engage in molecular structural interactions and that TGFβ1 inhibits CLDN4 transcription and protein expression in human colorectal cancer cell lines. Increased levels of TGFβ1 may affect the expression of CLDN4, potentially influencing CRC development. CLDN9 can be employed as a prognostic biomarker for gastric and esophageal cancer (Kang et al., 2020; Yu et al., 2020). CLDN10 overexpression promotes carcinogenesis in papillary thyroid cancer, which has a poor prognosis (Zhou et al., 2018). However, infiltration of B cells, CD8 + T cells, and macrophages as a consequence of increased CLDN10 expression improves the prognosis of patients with papillary thyroid carcinoma (Xiang et al., 2020). TGFβ or WNT/β-catenin-induced epithelial-mesenchymal transition (EMT) and tumor resistance may be linked to CLDN10, which offers promise as a prognostic biomarker in ovarian cancer (Gao et al., 2017; Li et al., 2020b). Hypermethylated CLDN11 can be utilized as a biomarker for the early detection of CRC melanoma, and gastric cancer (Agarwal et al., 2009; Karagiannis et al., 2014; Guo et al., 2019). While CLDN11 is perhaps the best-studied claudin in oncogenic contexts to date, it has yet to be identified as either a promoter or repressor of tumor growth (Bhat et al., 2020). CLDN12 is involved in the formation of Ca2+ channels in intestinal epithelial cells and is involved in Ca2+ uptake in intestinal epithelial cells (Fujita et al., 2008). According to prior studies, CLDN12 is the most uniformly expressed of these genes across murine organs (Hwang et al., 2014). By activating PI3K/AKT/mTOR, CLDN14 promotes CRC development (Qiao et al., 2021). CLDN15 has been shown to be a positive indicator associated with malignant pleural mesothelioma in clinical contexts (Watanabe et al., 2021), while CLDN16 is a predictor of oral squamous cell carcinoma, breast cancer, and thyroid cancer (Chen et al., 2010; Gomez-Rueda et al., 2016; Ribeiro et al., 2021). In gastric cancer, CLDN18 can be employed as a tumor marker (Matsusaka et al., 2016). Moreover, CRC patients with low levels of CLDN23 expression may benefit from its use as a prognostic, diagnostic, or therapeutic target (Oh et al., 2017; Liu et al., 2020).
Previous studies exploring the effects of CLDN1, CLDN2, CLDN4, CLDN11, CLDN12, CLDN14, and CLDN23 expression on CRC tumor growth have yielded findings consistent with our results in the present study. However, the influence of additional genes expressed in CRC on tumor growth remains unknown. In addition, as many of these studies are subject to various limitations, the same genes may not be consistently expressed in both databases analyzed in the present study, and the results may thus be skewed. Secondly, as we used an online database to analyze gene expression patterns, our findings remain to be verified in a large-scale colorectal cancer clinical study.
In conclusion, we herein found that many claudin family genes and proteins were dysregulated in CRC by querying multiple publicly available databases. Claudin family proteins have also been implicated in various phases of cancer development, lymph node metastasis, and immune cell infiltration. In addition, TGFβ1 was found to suppress the expression of CLDN4 in CRC cell lines. In general, claudin family genes contribute to the progression and prognosis of colorectal cancer, and they may thus be prospective targets for the development of novel pharmaceutical preparations capable of treating CRC.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics Statement
The studies involving human participants were reviewed and approved by The Ethic Committee of Hebei University of Chinese Medicine. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
WZ, PC and ZZ designed the study. WZ, PC and ZZ supervised the study and wrote the manuscript. LY and WZ performed the experiments. ML, KH, YW, ZQ and TS collected and analyzed the data. All authors contributed to the article and approved the submitted version.
Funding
The work was supported by the Key Basic Applied Project of Hebei Provincial Department of Science & Technology (grant NO. 15967730D) and Hebei Provincial Administration of Traditional Chinese Medicine (2022077).
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.
Publisher’s Note
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.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.783016/full#supplementary-material
References
Agarwal, R., Mori, Y., Cheng, Y., Jin, Z., Olaru, A. V., Hamilton, J. P., et al. (2009). Silencing of Claudin-11 Is Associated with Increased Invasiveness of Gastric Cancer Cells. PLoS One 4 (11), e8002. doi:10.1371/journal.pone.0008002
Akizuki, R., Shimobaba, S., Matsunaga, T., Endo, S., and Ikari, A. (2017). Claudin-5, -7, and -18 Suppress Proliferation Mediated by Inhibition of Phosphorylation of Akt in Human Lung Squamous Cell Carcinoma. Biochim. Biophys. Acta (Bba) - Mol. Cel Res. 1864, 293–302. doi:10.1016/j.bbamcr.2016.11.018
Arabzadeh, A., Troy, T.-C., and Turksen, K. (2007). Changes in the Distribution Pattern of Claudin Tight junction Proteins during the Progression of Mouse Skin Tumorigenesis. BMC Cancer 7, 196. doi:10.1186/1471-2407-7-196
Berman, H. M., Westbrook, J., Feng, Z., and Bourne, P. E. (2000). The Protein Data Bank. Nucleic Acids Res. 28 (1), 235–242. doi:10.1093/nar/28.1.235
Bhat, A. A., Syed, N., Therachiyil, L., Nisar, S., Hashem, S., Macha, M. A., et al. (2020). Claudin-1, A Double-Edged Sword in Cancer. Ijms 21 (2), 569. doi:10.3390/ijms21020569
Cerami, E., Gao, J., Dogrusoz, U., Gross, B. E., Sumer, S. O., Aksoy, B. A., et al. (2012). The cBio Cancer Genomics Portal: An Open Platform for Exploring Multidimensional Cancer Genomics Data: Figure 1. Cancer Discov. 2 (5), 401–404. doi:10.1158/2159-8290.CD-12-0095
Chen, S. J., Chien, S. Y., Lin, C., and Chen, D. R. (2010). Significant Elevation of CLDN16 and HAPLN3 Gene Expression in Human Breast Cancer. Oncol. Rep. 24 (3), 759–766. doi:10.3892/or_00000918
Cheng, B., Rong, A., Zhou, Q., and Li, W. (2019). CLDN8 Promotes Colorectal Cancer Cell Proliferation, Migration, and Invasion by Activating MAPK/ERK Signaling. Cmar 11, 3741–3751. doi:10.2147/CMAR.S189558
Cherradi, S., Martineau, P., Gongora, C., and Del Rio, M. (2019). Claudin Gene Expression Profiles and Clinical Value in Colorectal Tumors Classified According to Their Molecular Subtype. Cmar 11, 1337–1348. doi:10.2147/CMAR.S188192
Dhawan, P., Ahmad, R., Chaturvedi, R., Smith, J. J., Midha, R., Mittal, M. K., et al. (2011). Claudin-2 Expression Increases Tumorigenicity of colon Cancer Cells: Role of Epidermal Growth Factor Receptor Activation. Oncogene 30 (29), 3234–3247. doi:10.1038/onc.2011.43
Escudero-Esparza, A., Jiang, W. G., and Martin, T. A. (2012). Claudin-5 Is Involved in Breast Cancer Cell Motility through the N-WASP and ROCK Signalling Pathways. J. Exp. Clin. Cancer Res. 31, 43. doi:10.1186/1756-9966-31-43
Franceschini, A., Szklarczyk, D., Frankild, S., Kuhn, M., Simonovic, M., Roth, A., et al. (2013). STRING v9.1: Protein-Protein Interaction Networks, with Increased Coverage and Integration. Nucleic Acids Res. 41, D808–D815. doi:10.1093/nar/gks1094
Fujita, H., Sugimoto, K., Inatomi, S., Maeda, T., Osanai, M., Uchiyama, Y., et al. (2008). Tight junction Proteins Claudin-2 and -12 Are Critical for Vitamin D-dependent Ca2+ Absorption between Enterocytes. MBoC 19 (5), 1912–1921. doi:10.1091/mbc.e07-09-0973
Fujiwara-Tani, R., Sasaki, T., Luo, Y., Goto, K., Kawahara, I., Nishiguchi, Y., et al. (2018). Anti-claudin-4 Extracellular Domain Antibody Enhances the Antitumoral Effects of Chemotherapeutic and Antibody Drugs in Colorectal Cancer. Oncotarget 9 (100), 37367–37378. doi:10.18632/oncotarget.26427
Gaedcke, J., Grade, M., Jung, K., Camps, J., Jo, P., Emons, G., et al. (2010). Mutated KRAS Results in Overexpression of DUSP4, a MAP-Kinase Phosphatase, and SMYD3, a Histone Methyltransferase, in Rectal Carcinomas. Genes Chromosom. Cancer 49 (11), 1024–1034. doi:10.1002/gcc.20811
Gao, Y., Liu, X., Li, T., Wei, L., Yang, A., Lu, Y., et al. (2017). Cross-validation of Genes Potentially Associated with Overall Survival and Drug Resistance in Ovarian Cancer. Oncol. Rep. 37 (5), 3084–3092. doi:10.3892/or.2017.5534
Georges, R., Bergmann, F., Hamdi, H., Zepp, M., Eyol, E., Hielscher, T., et al. (2012). Sequential Biphasic Changes in Claudin1 and Claudin4 Expression Are Correlated to Colorectal Cancer Progression and Liver Metastasis. J. Cel Mol Med 16 (2), 260–272. doi:10.1111/j.1582-4934.2011.01289.x
Gomez-Rueda, H., Palacios-Corona, R., Gutiérrez-hermosillo, H., and Trevino, V. (2016). A Robust Biomarker of Differential Correlations Improves the Diagnosis of Cytologically Indeterminate Thyroid Cancers. Int. J. Mol. Med. 37 (5), 1355–1362. doi:10.3892/ijmm.2016.2534
Guo, Y., Long, J., and Lei, S. (2019). Promoter Methylation as Biomarkers for Diagnosis of Melanoma: A Systematic Review and Meta‐analysis. J. Cel Physiol 234 (5), 7356–7367. doi:10.1002/jcp.27495
Hewitt, K. J., Agarwal, R., and Morin, P. J. (2006). The Claudin Gene Family: Expression in normal and Neoplastic Tissues. BMC Cancer 6, 186. doi:10.1186/1471-2407-6-186
Huang, D., Sherman, B. T., Tan, Q., Collins, J. R., Alvord, W. G., Roayaei, J., et al. (2007). The DAVID Gene Functional Classification Tool: a Novel Biological Module-Centric Algorithm to Functionally Analyze Large Gene Lists. Genome Biol. 8 (9), R183. doi:10.1186/gb-2007-8-9-r183
Hwang, I., Yang, H., Kang, H.-S., Ahn, C.-H., Lee, G.-S., Hong, E.-J., et al. (2014). Spatial Expression of Claudin Family Members in Various Organs of Mice. Mol. Med. Rep. 9 (5), 1806–1812. doi:10.3892/mmr.2014.2031
Jakab, C., Rusvai, M., Gálfi, P., Halász, J., and Kulka, J. (2011). Expression of Claudin-5 in Canine Pancreatic Acinar Cell Carcinoma - an Immunohistochemical Study. Acta Vet. Hung 59 (1), 87–98. doi:10.1556/AVet.59.2011.1.8
Kaiser, S., Park, Y.-K., Franklin, J. L., Halberg, R. B., Yu, M., Jessen, W. J., et al. (2007). Transcriptional Recapitulation and Subversion of Embryonic colon Development by Mouse colon Tumor Models and Human colon Cancer. Genome Biol. 8 (7), R131. doi:10.1186/gb-2007-8-7-r131
Kanehisa, M., Sato, Y., Kawashima, M., Furumichi, M., and Tanabe, M. (2016). KEGG as a Reference Resource for Gene and Protein Annotation. Nucleic Acids Res. 44 (D1), D457–D462. doi:10.1093/nar/gkv1070
Kang, H., Wang, N., Wang, X., Zhang, Y., Lin, S., Mao, G., et al. (2020). A Glycolysis-Related Gene Signature Predicts Prognosis of Patients with Esophageal Adenocarcinoma. Aging 12 (24), 25828–25844. doi:10.18632/aging.104206
Karagiannis, G. S., Schaeffer, D. F., Cho, C.-K. J., Musrap, N., Saraon, P., Batruch, I., et al. (2014). Collective Migration of Cancer-Associated Fibroblasts Is Enhanced by Overexpression of Tight junction-associated Proteins Claudin-11 and Occludin. Mol. Oncol. 8 (2), 178–195. doi:10.1016/j.molonc.2013.10.008
Kyuno, D., Zhao, K., Schnölzer, M., Provaznik, J., Hackert, T., and Zöller, M. (2019). Claudin7‐dependent Exosome‐promoted Reprogramming of Nonmetastasizing Tumor Cells. Int. J. Cancer 145 (8), 2182–2200. doi:10.1002/ijc.32312
Li, D., Zhao, W., Zhang, X., Lv, H., Li, C., and Sun, L. (2021). NEFM DNA Methylation Correlates with Immune Infiltration and Survival in Breast Cancer. Clin. Epigenet 13 (1), 112. doi:10.1186/s13148-021-01096-4
Li, X., Li, H., Liu, C., Leng, X., Liu, T., Zhang, X., et al. (2020). CLDN6-mediates SB431542 Action through MMPs to Regulate the Invasion, Migration, and EMT of Breast Cancer Cells. Int. J. Clin. Exp. Pathol. 13 (7), 1590–1600.
Li, Z., Xuan, W., Huang, L., Chen, N., Hou, Z., Lu, B., et al. (2020). Claudin 10 Acts as a Novel Biomarker for the Prognosis of Patients with Ovarian Cancer. Oncol. Lett. 20 (1), 373–381. doi:10.3892/ol.2020.11557
Liu, H., Jiang, F., Jia, X., Lan, J., Guo, H., Li, E., et al. (2017). Cycling Hypoxia Affects Cell Invasion and Proliferation through Direct Regulation of Claudin1/Claudin7 Expression, and Indirect Regulation of P18 through Claudin7. Oncotarget 8 (6), 10298–10311. doi:10.18632/oncotarget.14397
Liu, X., Bing, Z., Wu, J., Zhang, J., Zhou, W., Ni, M., et al. (2020). Integrative Gene Expression Profiling Analysis to Investigate Potential Prognostic Biomarkers for Colorectal Cancer. Med. Sci. Monit. 26, e918906. doi:10.12659/MSM.918906
Lyu, L., Xiang, W., Zheng, F., Huang, T., Feng, Y., Yuan, J., et al. (2020). Significant Prognostic Value of the Autophagy-Related Gene P4HB in Bladder Urothelial Carcinoma. Front. Oncol. 10, 1613. doi:10.3389/fonc.2020.01613
Marincola Smith, P., Choksi, Y. A., Markham, N. O., Hanna, D. N., Zi, J., Weaver, C. J., et al. (2021). Colon Epithelial Cell TGFβ Signaling Modulates the Expression of Tight junction Proteins and Barrier Function in Mice. Am. J. Physiology-Gastrointestinal Liver Physiol. 320 (6), G936–G957. doi:10.1152/ajpgi.00053.2021
Mármol, I., Sánchez-de-Diego, C., Pradilla Dieste, A., Cerrada, E., and Rodriguez Yoldi, M. (2017). Colorectal Carcinoma: A General Overview and Future Perspectives in Colorectal Cancer. Ijms 18 (1), 197. doi:10.3390/ijms18010197
Matsusaka, K., Ushiku, T., Urabe, M., Fukuyo, M., Abe, H., Ishikawa, S., et al. (2016). Coupling CDH17 and CLDN18 Markers for Comprehensive Membrane-Targeted Detection of Human Gastric Cancer. Oncotarget 7 (39), 64168–64181. doi:10.18632/oncotarget.11638
Morel, A.-P., Hinkal, G. W., Thomas, C., Fauvet, F., Courtois-Cox, S., Wierinckx, A., et al. (2012). EMT Inducers Catalyze Malignant Transformation of Mammary Epithelial Cells and Drive Tumorigenesis towards Claudin-Low Tumors in Transgenic Mice. Plos Genet. 8 (5), e1002723. doi:10.1371/journal.pgen.1002723
Mori, M., Miyoshi, N., Ishii, H., and Mori, M. (2011). Expression of CLDN1 in Colorectal Cancer: a Novel Marker for Prognosis. Int. J. Oncol. 39 (4), 791–796. doi:10.3892/ijo.2011.1102
Mu, W., Wang, Z., and Zöller, M. (2019). Ping-Pong-Tumor and Host in Pancreatic Cancer Progression. Front. Oncol. 9, 1359. doi:10.3389/fonc.2019.01359
Mulligan, M. K., Mozhui, K., Prins, P., and Williams, R. W. (2017). GeneNetwork: A Toolbox for Systems Genetics. Methods Mol. Biol. 1488, 75–120. doi:10.1007/978-1-4939-6427-7_4
Oh, B. Y., Cho, J., Hong, H. K., Bae, J. S., Park, W.-Y., Joung, J.-G., et al. (2017). Exome and Transcriptome Sequencing Identifies Loss of PDLIM2 in Metastatic Colorectal Cancers. Cmar 9, 581–589. doi:10.2147/cmar.s149002
Perez, A., Andrade-Da-Costa, J. s., De Souza, W., De Souza Ferreira, M., Boroni, M., De Oliveira, I., et al. (2020). N-glycosylation and R-eceptor T-yrosine K-inase S-ignaling A-ffect C-laudin-3 L-evels in C-olorectal C-ancer C-ells. Oncol. Rep. 44 (4), 1649–1661. doi:10.3892/or.2020.7727
Phattarataratip, E., and Sappayatosok, K. (2016). Expression of Claudin-5, Claudin-7 and Occludin in Oral Squamous Cell Carcinoma and Their Clinico-Pathological Significance. J. Clin. Exp. Dent 8 (3), e299–306. doi:10.4317/jced.52801
Pierce, B. G., Hourai, Y., and Weng, Z. (2011). Accelerating Protein Docking in ZDOCK Using an Advanced 3D Convolution Library. PLoS One 6 (9), e24657. doi:10.1371/journal.pone.0024657
Pierce, B. G., Wiehe, K., Hwang, H., Kim, B.-H., Vreven, T., and Weng, Z. (2014). ZDOCK Server: Interactive Docking Prediction of Protein-Protein Complexes and Symmetric Multimers. Bioinformatics 30 (12), 1771–1773. doi:10.1093/bioinformatics/btu097
Pontén, F., Jirström, K., and Uhlen, M. (2008). The Human Protein Atlas-A Tool for Pathology. J. Pathol. 216 (4), 387–393. doi:10.1002/path.2440
Pope, J. L., Ahmad, R., Bhat, A. A., Washington, M. K., Singh, A. B., and Dhawan, P. (2014). Claudin-1 Overexpression in Intestinal Epithelial Cells Enhances Susceptibility to Adenamatous Polyposis Coli-Mediated colon Tumorigenesis. Mol. Cancer 13, 167. doi:10.1186/1476-4598-13-167
Qiao, T.-Y., Yuan, Z.-M., Ma, T.-Y., Hu, H.-Q., Zhu, Y.-H., Zhang, W.-Y., et al. (2021). Claudin14 Promotes Colorectal Cancer Progression via the PI3K/AKT/mTOR Pathway. neo 68, 947–954. doi:10.4149/neo_2021_210210N203
Rhodes, D. R., Yu, J., Shanker, K., Deshpande, N., Varambally, R., Ghosh, D., et al. (2004). ONCOMINE: a Cancer Microarray Database and Integrated Data-Mining Platform. Neoplasia 6 (1), 1–6. doi:10.1016/s1476-5586(04)80047-2
Ribeiro, I. P., Esteves, L., Santos, A., Barroso, L., Marques, F., Caramelo, F., et al. (2021). A Seven-Gene Signature to Predict the Prognosis of Oral Squamous Cell Carcinoma. Oncogene 40 (22), 3859–3869. doi:10.1038/s41388-021-01806-5
Sakaguchi, T., Suzuki, S., Higashi, H., Inaba, K., Nakamura, S., Baba, S., et al. (2008). Expression of Tight junction Protein Claudin-5 in Tumor Vessels and Sinusoidal Endothelium in Patients with Hepatocellular Carcinoma. J. Surg. Res. 147 (1), 123–131. doi:10.1016/j.jss.2007.07.013
Skrzypczak, M., Goryca, K., Rubel, T., Paziewska, A., Mikula, M., Jarosz, D., et al. (2010). Modeling Oncogenic Signaling in colon Tumors by Multidirectional Analyses of Microarray Data Directed for Maximization of Analytical Reliability. PLoS One 5 (10), e13091. doi:10.1371/journal.pone.0013091
Suzuki, H., Nishizawa, T., Tani, K., Yamazaki, Y., Tamura, A., Ishitani, R., et al. (2014). Crystal Structure of a Claudin Provides Insight into the Architecture of Tight Junctions. Science 344 (6181), 304–307. doi:10.1126/science.1248571
Tabariès, S., and Siegel, P. M. (2017). The Role of Claudins in Cancer Metastasis. Oncogene 36 (9), 1176–1190. doi:10.1038/onc.2016.289
Tang, W., Dou, T., Zhong, M., and Wu, Z. (2011). Dysregulation of Claudin Family Genes in Colorectal Cancer in a Chinese Population. Biofactors 37 (1), 65–73. doi:10.1002/biof.138
Tang, Z., Li, C., Kang, B., Gao, G., Li, C., and Zhang, Z. (2017). GEPIA: a Web Server for Cancer and normal Gene Expression Profiling and Interactive Analyses. Nucleic Acids Res. 45 (W1), W98–W102. doi:10.1093/nar/gkx247
Thomas, P. D., Hill, D. P., Mi, H., Osumi-Sutherland, D., Van Auken, K., Carbon, S., et al. (2019). Gene Ontology Causal Activity Modeling (GO-CAM) Moves beyond GO Annotations to Structured Descriptions of Biological Functions and Systems. Nat. Genet. 51 (10), 1429–1433. doi:10.1038/s41588-019-0500-1
UniProt Consortium Nucleic Acids Res. UniProt: the Universal Protein Knowledgebase in 2021. UniProt Consortium.2021 ;49(D1):D480–D489. doi:10.1093/nar/gkaa1100
Wang, Y.-B., Shi, Q., Li, G., Zheng, J.-H., Lin, J., and Qiu, W. (2019). MicroRNA-488 Inhibits Progression of Colorectal Cancer via Inhibition of the Mitogen-Activated Protein Kinase Pathway by Targeting Claudin-2. Am. J. Physiology-Cell Physiol. 316 (1), C33–C47. doi:10.1152/ajpcell.00047.2018
Wang, Z., and Zöller, M. (2019). Exosomes, Metastases, and the Miracle of Cancer Stem Cell Markers. Cancer Metastasis Rev. 38 (1-2), 259–295. doi:10.1007/s10555-019-09793-6
Watanabe, M., Higashi, T., Ozeki, K., Higashi, A. Y., Sugimoto, K., Mine, H., et al. (2021). CLDN15 Is a Novel Diagnostic Marker for Malignant Pleural Mesothelioma. Sci. Rep. 11 (1), 12554. doi:10.1038/s41598-021-91464-0
Xiang, Z., Zhong, C., Chang, A., Ling, J., Zhao, H., Zhou, W., et al. (2020). Immune-related Key Gene CLDN10 Correlates with Lymph Node Metastasis but Predicts Favorable Prognosis in Papillary Thyroid Carcinoma. Aging 12 (3), 2825–2839. doi:10.18632/aging.102780
Yamagishi, H., Kuroda, H., Imai, Y., Hiraishi, H., Yamagishi, H., Kuroda, H., et al. (2016). Molecular Pathogenesis of Sporadic Colorectal Cancers. Chin. J. Cancer 35, 4. doi:10.1186/s40880-015-0066-y
Yang, T., Wang, H., Li, M., Yang, L., Han, Y., Liu, C., et al. (2021). CD151 Promotes Colorectal Cancer Progression by a Crosstalk Involving CEACAM6, LGR5 and Wnt Signaling via TGFβ1. Int. J. Biol. Sci. 17 (3), 848–860. doi:10.7150/ijbs.53657
Yu, A. S. L., Cheng, M. H., Angelow, S., Günzel, D., Kanzawa, S. A., Schneeberger, E. E., et al. (2009). Molecular Basis for Cation Selectivity in Claudin-2-Based Paracellular Pores: Identification of an Electrostatic Interaction Site. J. Gen. Physiol. 133 (1), 111–127. doi:10.1085/jgp.200810154
Yu, S., Hu, C., Cai, L., Du, X., Lin, F., Yu, Q., et al. (2020). Seven-Gene Signature Based on Glycolysis Is Closely Related to the Prognosis and Tumor Immune Infiltration of Patients with Gastric Cancer. Front. Oncol. 10, 1778. doi:10.3389/fonc.2020.01778
Zhou, Y., Xiang, J., Bhandari, A., Guan, Y., Xia, E., Zhou, X., et al. (2018). CLDN10 Is Associated with Papillary Thyroid Cancer Progression. J. Cancer 9 (24), 4712–4717. doi:10.7150/jca.28636
Zhu, J., Wang, R., Cao, H., Zhang, H., Xu, S., Wang, A., et al. (2015). Expression of Claudin-5, -7, -8 and -9 in Cervical Carcinoma Tissues and Adjacent Non-neoplastic Tissues. Int. J. Clin. Exp. Pathol. 8 (8), 9479–9486.
Keywords: colorectal cancer, CLDNs, public databases, TGFβ1, molecular docking, prognostic value
Citation: Yang L, Zhang W, Li M, Dam J, Huang K, Wang Y, Qiu Z, Sun T, Chen P, Zhang Z and Zhang W (2022) Evaluation of the Prognostic Relevance of Differential Claudin Gene Expression Highlights Claudin-4 as Being Suppressed by TGFβ1 Inhibitor in Colorectal Cancer. Front. Genet. 13:783016. doi: 10.3389/fgene.2022.783016
Received: 25 September 2021; Accepted: 25 January 2022;
Published: 24 February 2022.
Edited by:
Xiao Zhu, Guangdong Medical University, ChinaReviewed by:
Jenny Paredes, Memorial Sloan Kettering Cancer Center, United StatesMohamed Abd Elhamid Tantawy, National Research Centre, Egypt
Copyright © 2022 Yang, Zhang, Li, Dam, Huang, Wang, Qiu, Sun, Chen, Zhang and Zhang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Wei Zhang, zhangwei@hebcm.edu.cn; Pingping Chen, chenpp@hebcm.edu.cn; Zhenduo Zhang, 289230116@qq.com