DSG2 expression is correlated with poor prognosis and promotes early-stage cervical cancer

Background The pathogenesis and developmental mechanism of early-stage (FIGO 2009 IA2-IIA2) cervical cancer (CC) remain unclear. Seeking novel molecular biomarkers based on The Cancer Genome Atlas (TCGA) will facilitate the understanding of CC pathogenesis and help evaluate early-stage CC prognosis. Methods To identify prognosis-related genes in early-stage CC, we analyzed TCGA mRNA-seq data and clinical data by univariate Cox and Kaplan–Meier plotter analyses. Differential expression analysis identified upregulated genes in early-stage CC. Combined with the genes correlated with unfavorable prognosis, we selected desmoglein-2 (DSG2) for further investigation. To detect DSG2 expression in early-stage CC, we used immunohistochemistry (IHC), quantitative real-time PCR (qRT-PCR) and western blotting. The relationship between the expression of DSG2 and clinical features was analyzed by the Chi square test. Cox analysis was applied to assess the relationship between CC overall survival (OS) and risk factors. The correlations between DSG2 expression and CC cell line proliferation and migration were investigated with Cell Counting Kit-8 (CCK-8) and migration assays. Results There were 416 prognosis-related genes in early-stage CC. DSG2, matrix metallopeptidase 1 (MMP1), carbonic anhydrase IX (CA9), homeobox A1 (HOXA1), and serine protease inhibitor B3 (SERPINB3) were upregulated in early-stage CC compared with adjacent noncancerous tissue (ANT) and correlated with unfavorable prognosis. Among them, DSG2 was most significantly correlated with patient survival. Coexpression analysis indicated that DSG2 was probably involved in cell division, positive regulation of transferase activity, positive regulation of cell migration, EGFR upregulation pathway and regulation of lymphangiogenesis. IHC, qRT-PCR and western blotting showed that DSG2 expression was higher in CC than in normal tissue. Significant correlations were identified between DSG2 expression and several aggressive clinical features, including pelvic lymph node metastasis (PLNM). Multivariate Cox analysis showed that DSG2 and PLNM were independent prognostic factors for OS. DSG2 knockdown inhibited CC cell proliferation and migration. Conclusions DSG2 is a biomarker that promotes tumor proliferation and metastasis and is correlated with poor prognosis in early-stage CC.

fourth among all female cancers [1]. In contrast to latestage CC patients, most early-stage (IA2-IIA2) CC patients have a significantly increased survival time after surgery and chemoradiotherapy. However, approximately 10-30% of early-stage patients were found to have pelvic lymph node metastasis (PLNM), and some of the patients eventually experienced adverse outcomes [2]. In earlystage CC, patients with moderately high-risk factors, including large tumor size (> 2 cm), poor differentiation, special pathologic types, deep stromal invasion, lymphovascular space invasion (LVSI), PLNM and parametrial infiltration, usually have relatively shorter survival times [2,3].
Currently, the pathogenesis and mechanism of CC metastasis remain unclear and probably involve the aberrant expression of numerous oncogenes and tumor suppressors. Rapid advances in molecular biotechnology revealed that some molecular biomarkers are related to the progression of CC [4]. Seeking novel molecular biomarkers of protein-coding genes would facilitate the understanding of CC pathogenesis and help us evaluate the prognosis of early-stage CC.
The Cancer Genome Atlas (TCGA) database has been developed in recent years. It is composed of a large amount of cancer mRNA-seq data as well as detailed clinical data, which makes bioinformatic data mining convenient and reliable [5]. We incorporated gene profiling, molecular signatures, and functional and pathway information with gene set enrichment analysis. Using bioinformatics analyses, we found a series of early-stage CC prognosis-related genes. Among all these genes, we found that desmoglein-2 (DSG2) was upregulated in early CC compared with normal samples and also predicted unfavorable prognosis in early CC.
DSG2 is a cell adhesion protein of the cadherin superfamily that is crucial for cardiomyocyte cohesion and function [6]. Its purpose is to regulate cell-cell contact with adjacent cells. The altered expression and function of desmosomal cadherins is associated with human tumorigenesis [7]. Brennan et al. [8] and Kurzen et al. [9] showed that DSG2 was more highly expressed in skin squamous cell carcinoma and basal cell carcinoma and that the positive rate was higher in high-risk patients. Kamekura et al. [10] showed that the downregulation of DSG2 inhibited the proliferation of colon cancer cells. Saaber et al. [11] showed that DSG2 was a novel biomarker of squamous cell lung carcinoma. Cai et al. [12] showed that DSG2 was more highly expressed in nonsmall cell lung cancer (NSCLC) and that the knockdown of DSG2 inhibited the progression of NSCLC. However, some studies have shown that DSG2 is expressed at lower levels in cancer and functions as a tumor suppressor. Yashiro et al. [13] showed that the high expression of DSG2 was correlated with a longer survival time among diffuse infiltrative carcinomas of the stomach. Ramani et al. [14] showed that the knockdown of DSG2 decreased the cell junction of pancreatic carcinoma cells and increased the rate of metastasis. Barber et al. [15] showed that the low expression of DSG2 was an independent prognostic factor for prostate cancer. Davies et al. [16] showed that lower-expressed DSG2 was correlated with poor differentiation, larger tumor size and lymph node metastasis in breast cancer.
However, the role of DSG2 in CC has never been explored. In the present study, we identified DSG2 as a novel CC prognosis-related gene using data mining. With clinical and cell line validation, we demonstrated that it probably increased the risk of PLNM and resulted in an unfavorable prognosis.

Datasets
The gene expression and clinicopathological data of 310 CC patients and 3 adjacent noncancerous tissues (ANTs) were downloaded from TCGA (https ://porta l.gdc.cance r.gov/) [17,18]. According to the TCGA publication guidelines (https ://www.cance r.gov/about -nci/organ izati on/ccg/resea rch/struc tural -genom ics/tcga), these mRNA sequencing data have no restrictions on publication, and no additional approval by an ethics committee was required to publish the use of the data.
With the Ensembl platform (http://www.ensem bl.org/), we separated the mRNAs from all the TCGA genes. Genes that had missing values in over 50% of the samples were removed. Finally, there were 12,084 genes included in the study. Samples without data on the survival state and survival time were also removed. Finally, 291 CC tissues, including 167 early-stage (FIGO 2009 IA2-IIA2) CC tissues and 3 ANTs, were included in the study. For the early-stage samples, any missing data on whether LVSI and corpus involvement occurred were all recorded as nonoccurrence (median of the available data).

Kaplan-Meier (KM), univariate Cox, Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) and protein-protein interaction (PPI) analyses
The prognostic value of each gene was calculated in the KM analyses and univariate Cox analyses for the earlystage cohort. A total of 416 genes with both P KM < 0.05 and P Cox < 0.05 were early-stage prognosis-related genes and were kept for further analyses. GO biological process, cellular component, and molecular function categories and KEGG pathway analyses and PPI network construction were conducted by the Metascape website (http://metas cape.org/gp/index .html), using false discovery rate (FDR) q-value < 0.05 as the standard for statistical significance.

Differential expression analyses (DEA)
To identify genes that are more highly expressed in early-stage CC than in ANTs, we performed a DEA of prognosis-related genes between 167 early-stage CC patients and 3 ANTs with the R package "DEseq 2". The differentially expressed mRNAs with log 2 |FC| > 1.5 and P-adjusted < 0.05 were considered to be significant. Hierarchical clustering analysis was applied to categorize the data into two groups with similar expression patterns between early-stage CC and ANTs.

Coexpression analyses
Coexpression analyses was conducted by the cBioPortal website (https ://www.cbiop ortal .org/). Using Spearman's correlation analyses, the genes with FDR q-value < 0.05 were regarded as coexpressed with DSG2. Then, GO biological process analysis and oncogenic signature analysis were conducted among the positively correlated genes (Spearman's correlation > 0) and negatively correlated genes (Spearman's correlation < 0) by the Metascape website.

Tissue sample collection
A total of 150 CC tissues, 6 ANTs and 30 normal cervical tissues (NCTs) collected from January 2006 to October 2012 were obtained from the archives of the Pathology Department and Gynecology Department of the First Affiliated Hospital of Sun Yat-sen University. All enrolled CC patients were matched from stage IA2 to IIA2 and underwent radical hysterectomy and lymphadenectomy. Only patients with no preoperative radiotherapy or chemotherapy and with available clinical follow-up data were enrolled. Thirty NCTs were collected from patients who underwent hysterectomy without malignant conditions. Written informed consent was obtained from each patient. All specimens were handled according to legal and ethical standards.

Cell lines and cell culture
In this study, SiHa, HeLa, C33A, CaSki, MS751 and ME180 cells were purchased from the American Type Culture Collection (ATCC, Rockville, MD, USA) and cultured according to their guidelines in a humidified atmosphere with 5% CO 2 at 37 °C. The SiHa, HeLa and ME180 cell lines were cultured in DMEM (Thermo Fisher, America). The CaSki cell line was cultured in RPMI 1640 medium (Thermo Fisher, America). The C33A and MS751 cell lines were cultured in Eagle's minimum essential medium (Thermo Fisher, America). The media were supplemented with 10% fetal bovine serum (Life Technology, America) and 1% antibiotics (100 U/ml penicillin and 100 µg/ml streptomycin) (Life Technology, America).

Immunohistochemistry (IHC)
For IHC, 4-µm paraffin-embedded sections were baked at 60 °C for 1 h, deparaffinized with xylene, rehydrated with a series of graded alcohols, and microwaved in EDTA antigen retrieval buffer. Then, the sections were blocked with 10% goat serum before incubation with a primary antibody at 4 °C overnight, followed by HRP-conjugated secondary antibody incubation for 30 min at room temperature. DAB was added to detect antibody binding. Once brown color appeared, the sections were immersed in distilled water to stop the reaction. The sections were counterstained with hematoxylin, dehydrated in graded alcohols and mounted. The primary antibodies were rabbit anti-human DSG2 monoclonal antibody (ab150372, Abcam, Britain) and mouse anti-human D2-40 monoclonal antibody (MAB-0567, MXB, China). The DSG2 staining results were scored based on the following criteria: (i) percentage of positive tumor cells in the tumor tissue: 0 (0%), 1 (1-10%), 2 (11-50%), 3 (51-70%) and 4 (71-100%); and (ii) staining intensity: 0 (none), 1 (weak), 2 (moderate), and 3 (strong). The staining index was calculated as the staining intensity score × the proportion of positive tumor cells (range from 0 to 12). The staining score of 6 was defined as the cutoff. Thus, patients with different positive staining levels of DSG2 expression were divided into low-and high-staining groups.

RNA extraction and quantitative real-time PCR (qRT-PCR)
Total RNA was extracted using Trizol reagent (TAKARA, Japan) according to the manufacturer's instructions, and the concentration of the RNA extracts of each sample was measured quantitatively by a NanoDrop ND-2000 spectrophotometer. RNA was reverse transcribed into cDNA by using PrimeScript RT Master Mix (TAKARA, Japan). cDNA was amplified and quantified using a 7500 Fast Real-Time PCR system (Applied Biosystems, USA) and SYBR Premix Ex Taq (TAKARA, Japan). The RT-PCR conditions for genes were set at 95 °C for 2 min, followed by 39 cycles at 95 °C for 20 s, 58 °C for 30 s and 72 °C for 30 s. The DSG2 sequences were 5′-CTC AGG TGT GCA GCC TAC TC-3′ (forward) and 5′-GTG GTG TTC CTA GCC GTC AT-3′ (reverse), while the GAPDH sequences were 5′-TGC ACC ACC AAC TGC TTA GC-3′ (forward) and 5′-GGC ATG GAC TGT GGT CAT GAG-3′ (reverse). qRT-PCR was repeated at least three times. mRNA expression was defined based on Ct, and relative expression levels were calculated using the comparative Ct (2 −ΔΔCt ) method after normalization with reference to the expression of the house-keeping gene GAPDH.

Western blot assay
Total protein was extracted with cold RIPA lysis buffer and fractionated by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and then transferred onto a 0.45-μm PVDF membrane (Millipore, America). The membranes were blocked with 5% skimmed milk and incubated with the primary antibody at 4 °C overnight, followed by secondary antibody incubation for 1 h at room temperature. Bound antibodies were detected with Immobilon Western Chemiluminescent HRP Substrate (Millipore, America). Rabbit antihuman DSG2 monoclonal antibody (ab150372, Abcam, Britain) and rabbit anti-human GAPDH antibody (XS20180808002, Bioworld, China) were used in this study.

Cell Counting Kit-8 (CCK-8) assay
For the CCK-8 assay, 5 × 10 3 SiHa and HeLa cells were seeded into each well of 96-well plates. The time calculation started when the cells adhered to the wall, and the wells were transfected with siRNA. Cell viability was measured at specific times by CCK-8 (CCK-8, DOJINDO, Japan). The absorbance value at 450 nm was read by a microplate reader (Tecan Sunrise, Tecan Group Ltd.).

Migration assay
The stable cell lines SiHa siRNA-NC, SiHa siRNA393, SiHa siRNA613, HeLa siRNA-NC, HeLa siRNA393 and HeLa siRNA613 were counted and then 10 × 10 4 stably infected SiHa cells and 20 × 10 4 stably infected HeLa cells in 250 µl of serum-free medium were separately plated into the upper chamber of 8-µm transwell inserts (BD Biosciences, Franklin Lakes, NJ), while 500 µl of medium containing 10% bovine serum albumin was added to the lower chamber. After 24 h of incubation at 37 °C, SiHa siRNA cells in the upper chamber were removed carefully. After 48 h of incubation at 37 °C, HeLa siRNA-NC and HeLa siRNA cells in the upper chamber were removed. Migrated cells on the lower membrane surface were fixed in 4% paraformaldehyde (Solarbio, Beijing, China) for 10 min and then stained with 0.1% crystal violet (KeyGEN biotech, Nanjing, China) for 10 min. The number of cells was counted in 5 randomly selected visual fields (100×) per well under an inverted microscope DMI4000B (Leica, Wetzlar, Germany).

Statistical analyses
Statistical analyses were performed using SPSS 22.0 statistical software (Chicago, IL, USA) and R version 3.6.0. The differences between two groups were analyzed by Student's t test. The differences among more than two groups were analyzed by ANOVA. The Chi square test and Fisher's exact test were used to analyze the relationship between DSG2 expression and the clinicopathological characteristics. Survival data were evaluated using univariate and multivariate Cox regression analyses. Survival curves were plotted by the KM method and compared using the log-rank test. In all cases, P < 0.05 was considered statistically significant.

Early-stage CC prognosis-related genes were identified by bioinformatic analyses
According to the KM plotter analyses and univariate Cox analyses, the TCGA data included 416 early-stage prognosis-related genes, including 217 protective (Cox coefficient < 0) and 199 hazardous (Cox coefficient > 0) genes. In the GO analyses using all survival-related genes, 4 biological process terms were significantly enriched (P-adjusted < 0.05), including the regulation of the mitotic cell cycle, the nucleobase-containing small molecule metabolic process, protein N-linked glycosylation, and organelle localization. GO cellular component analyses identified endoplasmic reticulum lumen as the significantly enriched signature (P-adjusted < 0.05). No signatures were significantly enriched in GO molecular function analysis. KEGG pathway analyses indicated that the significant pathways were purine metabolism, protein processing in the endoplasmic reticulum, and nucleotide excision repair (Fig. 1b, Additional file 1: Figure S1a). Each chromosome had different numbers of up-and downregulated prognosis-related genes (Additional file 1: Figure S1b). Additionally, we constructed a PPI network to interpret the potential biological roles of the prognosis-related mRNAs in early-stage CC (Additional file 1: Figure S1c). Five overlapping genes between hazardous (Cox coefficient > 0, P KM < 0.05 and P Cox < 0.05) genes and upregulated genes (Log 2 FC > 1.5, FDR < 0.05) in CC were identified. Except for DSG2, the other four genes have been explored in CC. Therefore, DSG2 was used for further validation in clinical samples and cells. The workflow for screening DSG2 and the 5 overlapping genes were shown in Fig. 1a and Table 1.

The potential functions of DSG2 in CC and other cancers were analyzed by bioinformatics
For further validation, we investigated the difference in DSG2 expression between normal tissue and CC based on Oncomine datasets (Fig. 2a). All datasets revealed that DSG2 was upregulated in the cancer group. Survival analyses of both the overall cohort and early-stage cohort showed that the expression of DSG2 predicted an unfavorable prognosis in CC (overall cohort: HR = 1.966, P = 0.006; early-stage cohort: HR = 2.122, P = 0.030) (Fig. 2c). Furthermore, the overall cohort survival analyses showed that the expression of DSG2 predicted an unfavorable prognosis in bladder urothelial carcinoma (BLCA), brain lower-grade glioma (LGG), lung adenocarcinoma (LUAD), pancreatic adenocarcinoma (PAAD) and uterine corpus endometrial carcinoma (UCEC), while predicting a favorable prognosis in colon adenocarcinoma (COAD), kidney renal clear cell carcinoma (KIRC) and kidney renal papillary cell carcinoma (KIRP) (all P < 0.05) (Fig. 2b). Genes that were coexpressed in conjunction with DSG2 were identified with cBioPortal analyses (P-adjusted < 0.05). There were 2610 positively correlated genes (Spearman's correlation > 0) and 2737 negatively correlated genes (Spearman's correlation < 0). The enriched GO biological process and oncogenic pathway items are shown in Fig. 2e and Additional file 1: Figure  S1d. According to the positively correlated gene enrichment, cell division, positive regulation of transferase activity, positive regulation of cell migration and the EGFR upregulation pathway were significantly enriched (P-adjusted < 0.05), revealing that DSG2 is involved in the process and metastasis of CC. Furthermore, two genes were significantly coexpressed with DSG2, CCBE1 and VASH1, which are genes that regulate lymphangiogenesis according to GO biological process analyses (Fig. 2d). DSG2 was positively correlated with CCBE1, which positively regulated lymphangiogenesis, while it was negatively correlated with VASH1, which negatively regulated lymphangiogenesis.
These results confirmed that DSG2 was important in the development of various cancers and was possibly an oncogenic gene in CC.

DSG2 expression is upregulated in CC tissues
IHC was performed on 150 early-stage CC samples and 30 NCTs, and the results revealed that DSG2 was more highly expressed in CC samples (Fig. 3b). Furthermore, DSG2 expression was upregulated in six early-stage CC samples compared to that in matched ANTs derived from the same patients (Fig. 3a).
Additionally, to validate the results above, the DSG2 mRNA expression level was determined in 20 NCTs and 20 early-stage CC tissues using qRT-PCR. Moreover, 3 NCT and 3 early-stage CC tissues were randomly selected from the abovementioned tissues for western blot analyses. A comparison of the results showed that DSG2 mRNA and protein levels were higher in earlystage CC tissues than in NCTs (Fig. 3c).

High expression of DSG2 is associated with poor clinical features and prognosis in early-stage CC
The correlation between DSG2 expression and clinicopathological features was analyzed according to the IHC score. High DSG2 expression was significantly correlated with several poor clinicopathological features, including Table 1 The up-regulated genes in early-stage CC tissue compared with ANT, correlated with unfavorable prognosis Gene Description  Table 2). No significant correlation was identified between DSG2 expression and age, FIGO stage, pathologic type, differentiation grade, stromal invasion, LVSI, vaginal involvement and parametrial infiltration ( Table 2). To verify the relationship between DSG2 expression and the prognosis of early-stage CC, univariate and multivariate Cox analyses were performed. Univariate analysis showed that DSG2 expression (P < 0.001), tumor size (P = 0.029), LVSI (P = 0.008) and PLNM (P < 0.001) were prognostic factors for overall survival (OS) ( Table 3). Multivariate analysis showed that DSG2 expression (P = 0.018) and PLNM (P = 0.006) were independent prognostic factors for OS (Table 4, Fig. 3e).

High DSG2 expression was correlated with the occurrence of PLNM
The IHC results showed that DSG2 expression was significantly correlated with PLNM ( Table 2, Fig. 3b). For further validation, qRT-PCR was performed to examine the mRNA levels of DSG2 in 20 PLNM and 20 non-PLNM tissues, while western blotting was performed to examine the protein levels of DSG2 in 3 PLNM and 3 non-PLNM tissues. Both the mRNA and protein levels of DSG2 in the PLNM group were higher than those in the non-PLNM group (P < 0.05) (Fig. 3d). This validation result was consistent with the IHC result. Moreover, to explore the mechanism of how DSG2 promoted PLNM, we detected the lymphatic microvessel density (LMVD) in the same IHC samples. We found that the high DSG2 expression group had higher LMVD than the low DSG2 group, indicating that DSG2 probably promoted PLNM by promoting lymphangiogenesis (Table 5, Fig. 3f ).

Knockdown of DSG2 expression decreased CC cell proliferation and migration
To determine the function of DSG2 in CC cell proliferation and migration, further investigation was performed using the CCK-8 assay and migration assay.
First, qRT-PCR and western blotting revealed that DSG2 expression in SiHa and HeLa cells was higher than that in other cells and NCT (Fig. 4a). Therefore, SiHa and HeLa were chosen for further experiments. DSG2 expression was downregulated in SiHa and HeLa cell lines by transfection of siRNA393 and siRNA613. The efficiencies of interference were confirmed by qRT-PCR and western blotting (Additional file 1: Figure S2). The CCK-8 assay showed that knockdown of DSG2 expression decreased the cell proliferative capacity (Fig. 4b). Migration assays showed that knockdown of DSG2 expression decreased cell migration (Fig. 4c).

Discussion
Using survival analyses and TCGA data, our study provided a series of prognosis-related genes of early-stage CC. With GO, KEGG and PPI analyses, we can determine the main function distribution and interaction of genes. We identified 5 genes that were upregulated in early CC compared with normal samples that were correlated with unfavorable prognosis, including DSG2, matrix metallopeptidase 1 (MMP1), carbonic anhydrase IX (CA9), homeobox A1 (HOXA1), and serine protease inhibitor B3 (SERPINB3). MMP1, CA9, HOXA1 and SERPINB3 had been explored. However, we could not find any studies exploring the relationship between DSG2 and CC. To the best of our knowledge, this is the first one.
Consistent with the TCGA data mining result, DSG2 was more highly expressed in CC tissue than in normal tissue in 4 Oncomine databases. By detecting DSG2 expression in tissues, we revealed that DSG2 was upregulated in CC tissue compared with ANT or NCT. Additionally, DSG2 was significantly correlated with tumor size, PLNM, recurrence and vital status in 5 years, but not FIGO stage, pathologic type, differentiation grade, stromal invasion, LVSI, vaginal involvement and parametrial infiltration. High DSG2 expression predicted an unfavorable prognosis in early-stage CC. As PLNM was the most important risk factor for CC development, we investigated the relationship between DSG2 and PLNM by IHC, qRT-PCR and western blot analyses. All experiments showed that DSG2 was more highly expressed in the PLNM group than in the non-PLNM group. Moreover, high DSG2 expression was associated with high LMVD. Furthermore, our in vitro studies demonstrated that knockdown of DSG2 inhibited the CC cell proliferative capacity and migration ability. In conclusion, DSG2 was a novel tumor promoter in CC, and probably promoted cancer development by promoting the occurrence of PLNM. However, some studies showed that the downregulation of DSG2 promoted the proliferation and metastasis of cancer cells because desmosome downregulation decreases adhesion junctions to drive tumor development and early invasion. The reasons why our results were contrary to some studies were probably as follows. First, as we found above, DSG2 played a different role in  DSG2 plays a different role in different kinds of cancer. With TCGA data mining, we found that high DSG2 expression was correlated with the unfavorable prognosis of BLCA, brain LGG, LUAD, PAAD and UCEC, while high DSG2 expression was correlated with the favorable prognosis of COAD, KIRC and KIRP. These findings were consistent with those of previous reports. DSG2 is probably a novel biomarker of cancers but has different functions in different cancers.
Desmosomal cadherins are a component in cell-cell junctions, which are involved in the process of intercellular communication, signal transduction and cell proliferation [20]. In addition to regulating cell adhesion, DSG2 influenced cell proliferation and invasion by regulating the signaling pathway. It could be upstream or downstream of a pathway. The coexpression analyses results showed that cell division, positive regulation of transferase activity, positive regulation of cell migration and the EGFR upregulation pathway were significantly enriched among the positively correlated genes, revealing that DSG2 is involved in the process and metastasis of CC. These results were consistent with those of previous reports. Cai et al. [12] showed that knockdown of DSG2 suppressed non-small cell lung cancer cell proliferation by targeting p27 and CDK2. Kamekura et al. [10] reported that DSG2 and DSC2 played opposite roles in colon cancer cell proliferation. The loss of DSG2 suppressed cell proliferation through the altered phosphorylation of EGFR, Src and Erk protein. Overmiller et al. [21] suggested that in skin squamous cell carcinoma, DSG2 stimulated cell growth and migration by positively regulating EGFR levels and signaling through a c-Src and Cav1-dependent mechanism using lipid rafts as signal modulatory platforms. Brennan-Crispi et al. [22] showed that in skin basal cell carcinoma and squamous cell carcinoma, DSG2 enhanced canonical hedgehog signaling downstream of Ptc1 to promote cancer development through the activation of phosphorylated Stat3 and regulation of Gli1 expression. Katharina et al. [23] identified a novel promigratory pathway of pancreatic cancer cells in which the loss of DSG2 reduces the levels of plakoglobin via deregulated MAPK signaling. All of the above results showed that DSG2 was involved in various signaling pathways, such as the EGFR and MAPK signaling pathways as well as cell cycle pathways, indicating its important function in signaling pathway regulation.
Our study was the first to investigate the relationship between DSG2 expression and lymphangiogenesis. Coexpression analyses showed that DSG2 was positively correlated with CCBE1, which positively regulated lymphangiogenesis, while it was negatively correlated with VASH1, which negatively regulated lymphangiogenesis. An experiment was conducted to detect LMVD in tissue, which has not been reported in previous studies, and high LMVD was found to be associated with high DSG2 expression, indicating that DSG2 probably increased the lymphangiogenesis of cancer.
In conclusion, our current study was the first to show that DSG2 was overexpressed in CC tumorigenesis and that DSG2 knockdown repressed CC cell proliferation and migration. However, further mechanisms and signaling pathways underlying the role of DSG2 in CC remain to be defined.

Conclusions
Based on the above data, we drew a conclusion that DSG2 was a biomarker that promotes CC cells proliferation and metastasis and is correlated with poor prognosis in early-stage CC. These findings facilitated us to discover novel targets for the therapy of patients with CC.

Availability of data and materials
The datasets used and analyzed during the current study are available from the corresponding author on reasonable request. c The effect of siRNA on the migration abilities of CC cells detected by migration assay. Original magnification: ×100. *P < 0.05; **P < 0.01; ***P < 0.001; ****P < 0.0001