A human mucosal melanoma organoid platform for modeling tumor heterogeneity and exploring immunotherapy combination options.

Mucosal melanoma (MM), an aggressive rare subtype of melanoma, is distinct from cutaneous melanoma and has poor prognoses. We addressed the lack of cell models for MM by establishing 30 organoids of human oral MM (OMM), which retained major histopathological and functional features of parental tumors. Organoid groups derived from chronologically or intratumorally distinct lesions within the same individual displayed heterogeneous genetics, expression profiles, and drug responses, indicating rapid tumor evolution and poor clinical response. Furthermore, transcriptome analysis revealed receptor tyrosine kinases (RTKs) signaling, particularly NGFR, a nerve growth factor receptor, was significantly up-regulated in OMMs and organoids from patients resistant to anti-programmed cell death protein 1 (anti-PD-1) therapy. Combining anti-PD-1 with anlotinib (a phase 2 multitarget RTK inhibitor for OMM) or NGFR knockdown enhanced the effective activity of immune cells in organoid-immune cell coculture systems. Together, our study suggested that OMM organoids serve as faithful models for exploring tumor evolution and immunotherapy combination strategies.

Predictive value of prognostic nutritional index (PNI) in recurrent or unresectable hepatocellular carcinoma received anti-PD1 therapy.

BACKGROUND: Clinical trials have shown that anti-PD1 therapy, either as a monotherapy or in combination, is effective and well-tolerated in patients with recurrent or unresectable hepatocellular carcinoma (HCC). In this study, we aimed to investigate the prognostic value of immune-nutritional biomarkers in measuring the effects of anti-PD1 therapy in these patients. METHODS: We enrolled and followed up with 85 patients diagnosed with advanced HCC who underwent anti-PD1 therapy at the First Medical Centre of Chinese People's Liberation Army (PLA) General Hospital between January 2016 and January 2021. The retrospective analysis aimed to determine whether immune-nutritional biomarkers could serve as promising prognostic indices in these patients. RESULTS: In this retrospective study, patients in the PNI-high group showed a better progression-free survival (PFS) compared to those in the PNI-low group (9.5 months vs. 4.2 months, P = 0.039). Similarly, the median overall survival (OS) was longer in the PNI-high group (23.9 months, 95%CI 17.45-30.35) than in the PNI-low group (11.7 months, 95%CI 9.27-14.13) (P = 0.002). These results were consistent with sub-analyses of the anti-PD1 therapy. Furthermore, both univariate and multivariate analyses indicated that a higher pre-treatment PNI ( > = 44.91) was a significant predictive factor for favorable outcomes in this patient cohort (HR = 0.411, P = 0.023). CONCLUSION: Our study suggests that pre-treatment PNI is a critical predictive factor in patients with recurrent or unresectable HCC undergoing anti-PD1 therapy.

Single-cell and spatial dissection of precancerous lesions underlying the initiation process of oral squamous cell carcinoma.

Precancerous lesions of the oral mucosa, especially those accompanied by moderate to severe dysplasia, contribute to the initiation of oral squamous cell carcinoma (OSCC). However, the cellular compositions and spatial organization of the precancerous stage and how these factors promote human OSCC initiation remain unclear. Here, we built a single-cell transcriptome atlas and a spatial transcriptome map after obtaining data from pairwise human oral mucosal biopsies of 9 individuals consisting of very early-stage OSCC, adjacent precancerous lesions with moderate to severe dysplasia, as well as a matched normal region. An altered epithelial gene-expression profile was identified which favored OSCC initiation. This observation was coupled with distinct fibroblast, monocytic, and regulatory T-cell subclusters involved in reshaping the microenvironment. In particular, a unique immune-inhibitory monocyte subtype and spatial-switching regulation of VEGF signaling were observed surrounding precancerous lesions, concertedly strengthening activities in promoting cancer initiation. Collectively, our work elucidated the cellular landscapes and roles of precancerous lesions underlying OSCC initiation, which is essential for understanding the entire OSCC initiation process and helps inform therapeutic strategies for cancer intervention.

A biomimetic nanoplatform for customized photothermal therapy of HNSCC evaluated on patient-derived xenograft models.

Cancer cell membrane (CCM) derived nanotechnology functionalizes nanoparticles (NPs) to recognize homologous cells, exhibiting translational potential in accurate tumor therapy. However, these nanoplatforms are majorly generated from fixed cell lines and are typically evaluated in cell line-derived subcutaneous-xenografts (CDX), ignoring the tumor heterogeneity and differentiation from inter- and intra- individuals and microenvironments between heterotopic- and orthotopic-tumors, limiting the therapeutic efficiency of such nanoplatforms. Herein, various biomimetic nanoplatforms (CCM-modified gold@Carbon, i.e., Au@C-CCM) were fabricated by coating CCMs of head and neck squamous cell carcinoma (HNSCC) cell lines and patient-derived cells on the surface of Au@C NP. The generated Au@C-CCMs were evaluated on corresponding CDX, tongue orthotopic xenograft (TOX), immune-competent primary and distant tumor models, and patient-derived xenograft (PDX) models. The Au@C-CCM generates a photothermal conversion efficiency up to 44.2% for primary HNSCC therapy and induced immunotherapy to inhibit metastasis via photothermal therapy-induced immunogenic cell death. The homologous CCM endowed the nanoplatforms with optimal targeting properties for the highest therapeutic efficiency, far above those with mismatched CCMs, resulting in distinct tumor ablation and tumor growth inhibition in all four models. This work reinforces the feasibility of biomimetic NPs combining modular designed CMs and functional cores for customized treatment of HNSCC, can be further extended to other malignant tumors therapy.

Pharmacogenomic landscape of head and neck squamous cell carcinoma informs precision oncology therapy.

Head and neck squamous cell carcinoma (HNSCC) is a common and frequently lethal cancer with few therapeutic options. In particular, there are few effective targeted therapies. Development of highly effective therapeutic strategies tailored to patients with HNSCC remains a pressing challenge. To address this, we present a pharmacogenomic study to facilitate precision treatments for patients with HNSCC. We established a large collection of 56 HNSCC patient-derived cells (PDCs), which recapitulated the molecular features of the original tumors. Pharmacological assessment of HNSCCs was conducted using a three-tiered high-throughput drug screening using 2248 compounds across these PDC models and an additional 18 immortalized cell lines. We integrated genomic, transcriptomic, and pharmacological analysis to predict biomarkers, gene-drug associations, and validated biomarkers. These results supported drug repurposing for multiple HNSCC subtypes, including the JAK2 inhibitor fedratinib, for low KRT18-expressing HNSCC cases, and the topoisomerase inhibitor mitoxantrone, for IL6R-activated HNSCC cases. Our results demonstrated concordance between susceptibility predictions from the PDCs and the matched patients' responses to standard clinical medication. Moreover, we identified and experimentally confirmed that high expression of ITGB1 elicited therapeutic resistance to docetaxel and high SOD1 expression conferred resistance to afatinib. We further validated ITGB1 as a predictive biomarker for the efficacy of docetaxel therapy in a phase 2 clinical trial. In summary, our study shows that this HNSCC cell resource, as well as the resulting pharmacogenomic profiles, is effective for biomarker discovery and for guiding precision oncology therapies in HNSCCs.

Identification of potential key molecules and signaling pathways for psoriasis based on weighted gene co-expression network analysis.

BACKGROUND: Psoriasis is a chronic inflammatory skin disease, the pathogenesis of which is more complicated and often requires long-term treatment. In particular, moderate to severe psoriasis usually requires systemic treatment. Psoriasis is also associated with many diseases, such as cardiometabolic diseases, malignant tumors, infections, and mood disorders. Psoriasis can appear at any age, and lead to a substantial burden for individuals and society. At present, psoriasis is still a treatable, but incurable, disease. Previous studies have found that microRNAs (miRNAs) play an important regulatory role in the progression of various diseases. Currently, miRNAs studies in psoriasis and dermatology are relatively new. Therefore, the identification of key miRNAs in psoriasis is helpful to elucidate the molecular mechanism of psoriasis. AIM: To identify key molecular markers and signaling pathways to provide potential basis for the treatment and management of psoriasis. METHODS: The miRNA and mRNA data were obtained from the Gene Expression Omnibus database. Then, differentially expressed mRNAs (DEmRNAs) and differentially expressed miRNAs (DEmiRNAs) were screened out by limma R package. Subsequently, DEmRNAs were analyzed for Gene Ontology and Kyoto Encyclopedia of Genes and Genomics functional enrichment. The "WGCNA" R package was used to analyze the co-expression network of all miRNAs. In addition, we constructed miRNA-mRNA regulatory networks based on identified hub miRNAs. Finally, in vitro validation was performed. All experimental procedures were approved by the ethics committee of Chinese PLA General Hospital (S2021-012-01). RESULTS: A total of 639 DEmRNAs and 84 DEmiRNAs were identified. DEmRNAs screening criteria were adjusted P (adj. P) value < 0.01 and |logFoldChange| (|logFC|) > 1. DEmiRNAs screening criteria were adj. P value < 0.01 and |logFC| > 1.5. KEGG functional analysis demonstrated that DEmRNAs were significantly enriched in immune-related biological functions, for example, toll-like receptor signaling pathway, cytokine-cytokine receptor interaction, and chemokine signaling pathway. In weighted gene co-expression network analysis, turquoise module was the hub module. Moreover, 10 hub miRNAs were identified. Among these 10 hub miRNAs, only 8 hub miRNAs predicted the corresponding target mRNAs. 97 negatively regulated miRNA-mRNA pairs were involved in the miRNA-mRNA regulatory network, for example, hsa-miR-21-5p-claudin 8 (CLDN8), hsa-miR-30a-3p-interleukin-1B (IL-1B), and hsa-miR-181a-5p/hsa-miR-30c-2-3p-C-X-C motif chemokine ligand 9 (CXCL9). Real-time polymerase chain reaction results showed that IL-1B and CXCL9 were up-regulated and CLDN8 was down-regulated in psoriasis with statistically significant differences. CONCLUSION: The identification of potential key molecular markers and signaling pathways provides potential research directions for further understanding the molecular mechanisms of psoriasis. This may also provide new research ideas for the prevention and treatment of psoriasis in the future.

Clinical utility of PDX cohorts to reveal biomarkers of intrinsic resistance and clonal architecture changes underlying acquired resistance to cetuximab in HNSCC.

Cetuximab is a widely used drug for treating head and neck squamous cell carcinomas (HNSCCs); however, it provides restricted clinical benefits, and its response duration is limited by drug resistance. Here, we conducted randomized "Phase II-like clinical trials" of 49 HNSCC PDX models and reveal multiple informative biomarkers for intrinsic resistance to cetuximab (e.g., amplification of ANKH, up-regulation of PARP3). After validating these intrinsic resistance biomarkers in another HNSCC PDX cohort (61 PDX models), we generated acquired cetuximab resistance PDX models and analyzed them to uncover resistance mechanisms. Whole exome sequencing and transcriptome sequencing revealed diverse patterns of clonal selection in acquired resistant PDXs, including the emergence of subclones with strongly activated RAS/MAPK. Extending these insights, we show that a combination of a RAC1/RAC3 dual-target inhibitor and cetuximab could overcome acquired cetuximab resistance in vitro and in vivo. Beyond revealing intrinsic resistance biomarkers, our PDX-based study shows how clonal architecture changes underlying acquired resistance can be targeted to expand the therapeutic utility of this important drug to more HNSCC patients.

Retinoic Acid-Induced 2 (RAI2) Is a Novel Antagonist of Wnt/ß-Catenin Signaling Pathway and Potential Biomarker of Chemosensitivity in Colorectal Cancer.

Objective: Aberrant activation of Wnt/ß-catenin signaling contributes to the maintenance of cancer stem cells and chemoresistance in colorectal cancer (CRC). Retinoic acid-induced 2 (RAI2) was proved to be a tumor suppressor in CRC in our previous report. In this study, the role of RAI2 in Wnt/ß-catenin signaling was further investigated. Methods: As a transcriptional co-regulator, C-terminal Binding Protein 2 (CtBP2) was reported to be involved in Wnt signaling in multiple and complex ways. The correlation of RAI2 and CtBP2 in CRC was analyzed by TCGA dataset, and the interaction between RAI2 and CtBP2 was explored by co-immunoprecipitation (Co-IP) in CRC cells. The effect of RAI2 on the activity of Wnt signaling and the location of ß-catenin was detected by Dual-Luciferase reporter assay and Immunofluorescence respectively. Western blotting analysis was performed to detect the expression of target genes involved in Wnt signaling. Sphere formation assay was employed to detect the effect of RAI2 on stem cell like properties. Cell viability assay was used to detect the chemosensitivity of cells before and after transfection of RAI2. Results: The interaction between RAI2 and CtBP2 was confirmed by Co-IP in CRC cells. Besides, the negative correlation of RAI2 and CtBP2 in CRC was found by analyzing the TCGA dataset. Re-expression of RAI2 in human colon cancer cells (HCT116 and LoVo) suppressed the fluorescent activity of Wnt signaling, increased the phosphorylation and inhibited nuclear translocation of ß-catenin, with down-regulation of target genes like c-Myc, CyclinD1, ASCL2, and LGR5. In contrast, the mutated RAI2, which can't interact with CtBP2, has no above effects. We observed low expression of RAI2 in 33.89% (101/298) of CRC patients, which was significantly associated with reduced phosphorylation of ß-catenin (r=0.8866, P<0.0001), poor 5-year relapse-free survival (RFS) (P = 0.0029) and overall survival (OS) (P = 0.0102). Restoration of RAI2 in HCT116 and LoVo cells inhibited stem cell-like properties of CRC cells and increased chemosensitivity of these cells to oxaliplatin and fluorouracil. Conclusion: Low expression of RAI2 can serve as an independent poor prognostic marker. RAI2 inhibits Wnt signaling by interacting with or down-regulating CtBP2, resulting in repression of stem cell-like properties and increased chemosensitivity of CRC cells.

A Retrospective Study of Biliary Drainage Strategies for Patients with Malignant Hilar Biliary Strictures.

PURPOSE: The main aim of this study was to compare the efficacy and safety of different biliary drainage strategies, including percutaneous transhepatic biliary drainage (PTBD) versus endoscopic biliary stenting (EBS) and unilateral versus bilateral stenting, in patients with unresectable malignant hilar biliary strictures (MHBSs). PATIENTS AND METHODS: This was a retrospective review of patients with inoperable MHBSs who underwent biliary drainage by either EBS or PTBD. Efficacy and safety were compared between the two pathways and between unilateral and bilateral stenting in the EBS group. The survival duration was analyzed with K-M curves and Log rank tests. RESULTS: From January 2015 to December 2019, a total of 206 (126: EBS and 80: PTBD) patients with MHBSs were enrolled in our study and underwent 270 procedures (173: EBS and 97: PTBD). Bilateral stenting was superior to unilateral stenting in terms of clinical success (69.6% vs 50.6%, p=0.039), especially for patients with Bismuth type IV (70.0% vs 30.3%, p=0.002). A higher decrease in bilirubin was seen with PTBD in patients with Bismuth types III-IV (66.9 vs 36.7, p=0.006). A survival advantage was seen in successful drainage (227 days vs 82 days, p<0.001), lower tumor-node-metastasis (TNM) scores (I-II) (195 days vs 139 days, p=0.012), and cholangiocarcinoma (184 days vs 84 days, p=0.001). CONCLUSION: For patients with advanced MHBSs, bilateral stenting may achieve a better drainage effect than unilateral stenting, and PTBD may have a better performance in relieving cholestasis than EBS. Successful drainage and cholangiocarcinoma may provide greater long-term survival benefits.

Erratum: Effects of HPV16E6 transfection on the biological behavior of Eca109 and Eca9706 cells.

[This corrects the article DOI: 10.3892/ol.2017.7469.].

Screening and identification of key genes between liver hepatocellular carcinoma (LIHC) and cholangiocarcinoma (CHOL) by bioinformatic analysis.

BACKGROUND: Liver hepatocellular carcinoma (LIHC) and cholangiocarcinoma (CHOL) are common primary liver cancers worldwide. Liver stem cells have biopotential to differentiate into either hepatocytes and cholangiocytes, the phenotypic overlap between LIHC and CHOL has been acceptable as a continuous liver cancer spectrum. However, few studies directly investigated the underlying molecular mechanisms between LIHC and CHOL. METHOD: To identify the candidate genes between LIHC and CHOL, three data series including GSE31370, GSE15765 and GSE40367 were downloaded from Gene Expression Omnibus (GEO) database. The differentially expressed genes (DEGs) were identified, and function enrichment analyses were performed. The protein-protein interaction network (PPI) was constructed and the module analysis was performed using STRING and Cytoscape. RESULTS: A total of 171 DEGs were identified, consisting of 49 downregulated genes and 122 upregulated genes. Compared with CHOL, the enriched functions of the DEGs mainly included steroid metabolic process, acute inflammatory response, coagulation. Meanwhile, the pathway of KEGG enrichment analyses showed that the upregulated gene(s) were mainly enriched complement and coagulation cascades, cholesterol metabolism and PPAR signaling pathway, while the downregulated gene(s) were mainly enriched in ECM-receptor interaction, focal adhesion, bile secretion. Similarly, the most significant module was identified and biological process analysis revealed that these genes were mainly enriched in regulation of blood coagulation, acute inflammatory response, complement and coagulation cascades. Finally, two (ITIH2 and APOA2) of 10 hub genes had been screened out to help differential diagnosis. CONCLUSION: 171 DEGs and two (ITIH2 and APOA2) of 10 hub genes identified in the present study help us understand the different molecular mechanisms between LIHC and CHOL, and provide candidate targets for differential diagnosis.

Immune-Stromal Score Signature: Novel Prognostic Tool of the Tumor Microenvironment in Lung Adenocarcinoma.

Background: Immune and stromal cells in the tumor microenvironment (TME) significantly contribute to the prognosis of lung adenocarcinoma; however, the TME-related immune prognostic signature is unknown. The aim of this study was to develop a novel immune prognostic model of the TME in lung adenocarcinoma. Methods: First, the immune and stromal scores among lung adenocarcinoma patients were determined using the ESTIMATE algorithm in accordance with The Cancer Genome Atlas (TCGA) database. Differentially expressed immune-related genes (IRGs) between high and low immune/stromal score groups were analyzed, and a univariate Cox regression analysis was performed to identify IRGs significantly correlated with overall survival (OS) among patients with lung adenocarcinoma. Furthermore, a least absolute shrinkage and selection operator (LASSO) regression analysis was performed to generate TME-related immune prognostic signatures. Gene set enrichment analysis was performed to analyze the mechanisms underlying these immune prognostic signatures. Finally, the functions of hub IRGs were further analyzed to delineate the potential prognostic mechanisms in comprehensive TCGA datasets. Results: In total, 702 intersecting differentially expressed IRGs (589 upregulated and 113 downregulated) were screened. Univariate Cox regression analysis revealed that 58 significant differentially expressed IRGs were correlated with patient prognosis in the training cohort, of which three IRGs (CLEC17A, INHA, and XIRP1) were identified through LASSO regression analysis. A robust prognostic model was generated on the basis of this three-IRG signature. Furthermore, functional enrichment analysis of the high-risk-score group was performed primarily on the basis of metabolic pathways, whereas analysis of the low-risk-score group was performed primarily on the basis of immunoregulation and immune cell activation. Finally, hub IRGs CLEC17A, INHA, and XIRP1 were considered novel prognostic biomarkers for lung adenocarcinoma. These hub genes had different mutation frequencies and forms in lung adenocarcinoma and participated in different signaling pathways. More importantly, these hub genes were significantly correlated with the infiltration of CD4+ T cells, CD8+ T cells, macrophages, B cells, and neutrophils. Conclusions: The robust novel TME-related immune prognostic signature effectively predicted the prognosis of patients with lung adenocarcinoma. Further studies are required to further elucidate the regulatory mechanisms of these hub IRGs in the TME and to develop new treatment strategies.

Identification of candidate genes and prognostic value analysis in patients with PDL1-positive and PDL1-negative lung adenocarcinoma.

BACKGROUND: Increasing bodies of evidence reveal that targeting a programmed cell death protein 1 (PD-1) monoclonal antibody is a promising immunotherapy for lung adenocarcinoma. Although PD receptor ligand 1 (PDL1) expression is widely recognized as the most powerful predictive biomarker for anti-PD-1 therapy, its regulatory mechanisms in lung adenocarcinoma remain unclear. Therefore, we conducted this study to explore differentially expressed genes (DEGs) and elucidate the regulatory mechanism of PDL1 in lung adenocarcinoma. METHODS: The GSE99995 data set was obtained from the Gene Expression Omnibus (GEO) database. Patients with and without PDL1 expression were divided into PDL1-positive and PDL1-negative groups, respectively. DEGs were screened using R. The Gene Ontology (GO) database and Kyoto Encyclopedia of Genes and Genomes (KEGG) were analyzed using the Database for Annotation, Visualization and Integrated Discovery. Protein-protein interaction (PPI) networks of DEGs was visualized using Cytoscape, and the MNC algorithm was applied to screen hub genes. A survival analysis involving Gene Expression Profiling Interactive Analysis was used to verify the GEO results. Mutation characteristics of the hub genes were further analyzed in a combined study of five datasets in The Cancer Genome Atlas (TCGA) database. RESULTS: In total, 869 DEGs were identified, 387 in the PDL1-positive group and 482 in the PDL1-negative group. GO and KEGG analysis results of the PDL1-positive group mainly exhibited enrichment of biological processes and pathways related to cell adhesion and the peroxisome proliferators-activated receptors (PPAR) signaling pathway, whereas biological process and pathways associated with cell division and repair were mainly enriched in the PDL1-negative group. The top 10 hub genes were screened during the PPI network analysis. Notably, survival analysis revealed BRCA1, mainly involved in cell cycle and DNA damage responses, to be a novel prognostic indicator in lung adenocarcinoma. Moreover, the prognosis of patients with different forms of lung adenocarcinoma was associated with differences in mutations and pathways in potential hub genes. CONCLUSIONS: PDL1-positive lung adenocarcinoma and PDL1-negative lung adenocarcinoma might be different subtypes of lung adenocarcinoma. The hub genes might play an important role in PDL1 regulatory pathways. Further studies on hub genes are warranted to reveal new mechanisms underlying the regulation of PDL1 expression. These results are crucial for understanding and applying precision immunotherapy for lung adenocarcinoma.

Clinical Analysis and Prognostic Significance of Hepatitis B Virus Infections with Diffuse Large B-Cell Lymphoma.

BACKGROUND: China is a high endemic area for the hepatitis B virus (HBV). The studies established the epidemiology between HBV and diffuse large B-cell lymphoma (DLBCL), but further research is necessary to clarify the potential link between HBV and DLBCL. PATIENTS AND METHODS: A total of 319 patients diagnosed with DLBCL were recruited as cases at First Medical Centre of Chinese PLA General Hospital from September 2010 to December 2018. During the same time, two age- and sex-matched controls were selected for each case, and the control groups comprised of 319 patients with non-hematological malignancy and 319 subjects with non-malignant conditions. Relative risk of developing DLBCL among individuals tested positive for hepatitis B surface antigen was calculated. After that, we retrospectively analyzed clinical data from DLBCL patients with different HBV infection statuses. RESULTS: The HBV infection rate of patients with DLBCL (11.60%) was significantly higher than the other two control groups (5.02% and 4.08%), indicating the risk of DLBCL may increased in HBV infections. Meanwhile, this study demonstrated an independent association between HBV infection and poorer clinical outcomes. CONCLUSION: Our study demonstrated that HBV infection may play an important role in the pathogenesis of DLBCL and show poor outcomes in HBV-endemic China.

Pseudo progression in heterogeneity of right-sided metastatic colon carcinoma during nivolumab therapy: A case report.

RATIONALE: Pseudo progression is a noted phenomenon of immune checkpoint inhibitors therapy, which has been defined as a response after an initial enlargement of the tumor followed by tumor reduction. In July 2017, the Food and Drug Administration granted accelerated approval of nivolumab for the treatment of metastatic colorectal cancer patients whose tumor harbors deficient mismatch repair. PATIENT CONCERNS AND DIAGNOSIS: We present a patient who received nivolumab for heterogeneity of right-sided metastatic colon carcinoma. INTERVENTION: The patient was treated with nivolumab combined with chemotherapy. OUTCOME: The computed tomography showed mass lesion in the left lobe of liver remained stable while metastasis tumors under envelop of liver were exacerbated after 6 cycles of nivolumab combined with chemotherapy, and later regressed. LESSONS: The status of mismatch repair in primary tumor and metastatic liver carcinoma is contradictory but using nivolumab demonstrated encouraging efficacy. This is the first case of pseudo progression undergoing immunotherapy for heterogeneity of right-sided metastatic colon carcinoma.

Effects of HPV16E6 transfection on the biological behavior of Eca109 and Eca9706 cells.

The aim of the present study was to investigate the effect of HPV16E6 gene integration on the biological behavior of Eca109 and Eca9706 cells. This was evaluated through positive liposome transfection of HPV16E6 into Eca109 cells and Eca9706 cells. The transfection efficiency was evaluated by calculating the ratio of fluorescent cells to total cells. After stable screening, reverse transcription-polymerase chain reaction (RT-PCR) was performed to detect the target gene and HPV16E6mRNA inside the cells. The distribution of HPV16E6 in esophageal carcinoma cells was observed by immunofluorescence and western blot analysis. A CCK-8 assay was performed to detect cell proliferation. The migration rate was measured by a wound healing assay, and a Transwell Matrigel invasion assay was used to detect invasive ability. The results of RT-PCR, immunofluorescence and western blot analyses indicated that HPV16E6 gene was expressed in the target group. The proliferation rates and clone group numbers were significantly higher in HPV16E6-transfected cell groups compared with nonsense-transfected (negative control) cell groups. The wound healing and Transwell invasion assays indicated that the migration rate and invasive ability were also significantly higher in the HPV16E6-transfected cell groups compared with negative control groups. In conclusion, Eca109 and Eca9706 cell lines with integration of HPV16E6 were successfully established in the present study. It was demonstrated that HPV16E6 expression enhanced the proliferation and migration of esophageal cancer cells. HPV16E6 may serve a key function in the occurrence and development of esophageal cancer.