The expression profile and clinic significance of the SIX family in non-small cell lung cancer

All eyes on Eya: A unique transcriptional co-activator and phosphatase in cancer

The Eya family of proteins (consisting of Eyas1-4 in mammals) play vital roles in embryogenesis by regulating processes such as proliferation, migration/invasion, cellular survival and pluripotency/plasticity of epithelial and mesenchymal states. Eya proteins carry out such diverse functions through a unique combination of transcriptional co-factor, Tyr phosphatase, and PP2A/B55α-mediated Ser/Thr phosphatase activities. Since their initial discovery, re-expression of Eyas has been observed in numerous tumor types, where they are known to promote tumor progression through a combination of their transcriptional and enzymatic activities. Eya proteins thus reinstate developmental processes during malignancy and represent a compelling class of therapeutic targets for inhibiting tumor progression.

NIK-mediated reactivation of SIX2 enhanced the CSC-like traits of hepatocellular carcinoma cells through suppressing ubiquitin-proteasome system

The critical roles of NF-κB Inducing Kinase (NIK) in tumor progression have been elucidated in various tumors; however, its effects on hepatocellular carcinoma (HCC) progression are still confusing. Here, we found that NIK level was upregulated in HCC tissues compared to that of normal tissues, and positively correlated with the levels of cancer stem cell (CSC) markers. Then we established HCC cells with NIK-stable knockdown and found that NIK knockdown suppressed the CSC-like traits of HCC cells through in vivo and in vitro experiments. Mechanistically, we revealed that SIX2 protein level, but not its mRNA level, was significantly reduced in HCC cells with NIK knockdown, which was rescued by MG132 treatment. Furthermore, NIK knockdown promoted the ubiquitination level of SIX2 and decreased its protein stability. Moreover, Six2 overexpression partially reversed the inhibition of NIK knockdown on the CSC-like traits of HCC cells. This study identified a novel NIK/SIX2 axis conferring HCC stemness.

S100A8/A9 as a risk factor for breast cancer negatively regulated by DACH1

BACKGROUND

S100A8 and S100A9 are members of Ca2+-binding EF-hand superfamily, mainly expressed by macrophages and neutrophils. Limited by the poor stability of homodimers, they commonly exist as heterodimers. Beyond acting as antibacterial cytokines, S100A8/A9 is also associated with metabolic and autoimmune diseases such as obesity, diabetes, and rheumatoid arthritis. While the involvement of S100A8/A9 in breast cancer development has been documented, its prognostic significance and the precise regulatory mechanisms remain unclear.

METHODS

S100A8/A9 protein in breast cancer samples was evaluated by immunohistochemistry staining with tumor tissue microarrays. The serum S100A8 concentration in patients was measured by enzyme-linked immunosorbent assay (ELISA). The S100A8 secreted by breast cancer cells was detected by ELISA as well. Pooled analyses were conducted to explore the relationships between S100A8/A9 mRNA level and clinicopathological features of breast cancer patients. Besides, the effects of S100A8/A9 and DACH1 on patient outcomes were analyzed by tissue assays. Finally, xenograft tumor assays were adopted to validate the effects of DACH1 on tumor growth and S100A8/A9 expression.

RESULTS

The level of S100A8/A9 was higher in breast cancer, relative to normal tissue. Increased S100A8/A9 was related to poor differentiation grade, loss of hormone receptors, and Her2 positive. Moreover, elevated S100A8/A9 predicted a worse prognosis for breast cancer patients. Meanwhile, serum S100A8 concentration was upregulated in Grade 3, basal-like, and Her2-overexpressed subtypes. Additionally, the results of public databases showed S100A8/A9 mRNA level was negatively correlated to DACH1. Stable overexpressing DACH1 in breast cancer cells significantly decreased the generation of S100A8. The survival analysis demonstrated that patients with high S100A8/A9 and low DACH1 achieved the shortest overall survival. The xenograft models indicated that DACH1 expression significantly retarded tumor growth and downregulated S100A8/A9 protein abundance.

CONCLUSION

S100A8/A9 is remarkedly increased in basal-like and Her2-overexpressed subtypes, predicting poor prognosis of breast cancer patients. Tumor suppressor DACH1 inhibits S100A8/A9 expression. The combination of S100A8/A9 and DACH1 predicted the overall survival of breast cancer patients more preciously.

Forkhead box proteins as the critical regulators of cisplatin response in tumor cells

Cisplatin (CDDP) is one of the most common chemotherapy drugs used in a wide range of cancer patients; however, there is a high rate of CDDP resistance among cancer patients. Considering the side effects of cisplatin in normal tissues, it is necessary to predict the CDDP response in cancer patients. Therefore, identifying the molecular mechanisms involved in CDDP resistance can help to introduce the prognostic markers. Several molecular mechanisms such as apoptosis inhibition, drug efflux, drug detoxification, and increased DNA repair are involved in CDDP resistance. Regarding the key role of transcription factors in regulation of many cellular processes related to drug resistance, in the present review, we discussed the role of Forkhead box (FOX) protein family in CDDP response. It has been reported that FOX proteins mainly promote CDDP resistance through the regulation of DNA repair, autophagy, epithelial-mesenchymal transition (EMT), and signaling pathways. Therefore, FOX proteins can be introduced as the prognostic markers to predict CDDP response in cancer patients. In addition, considering that oncogenic role of FOX proteins, the CDDP treatment along with FOX inhibition can be used as a therapeutic strategy in cancer patients.

Comprehensive DNA methylation profiling of COVID-19 and hepatocellular carcinoma to identify common pathogenesis and potential therapeutic targets

BACKGROUND & AIMS

The effects of SARS-CoV-2 infection can be more complex and severe in patients with hepatocellular carcinoma (HCC) as compared to other cancers. This is due to several factors, including pre-existing conditions such as viral hepatitis and cirrhosis, which are commonly associated with HCC.

METHODS

We conducted an analysis of epigenomics in SARS-CoV-2 infection and HCC patients, and identified common pathogenic mechanisms using weighted gene co-expression network analysis (WGCNA) and other analyses. Hub genes were identified and analyzed using LASSO regression. Additionally, drug candidates and their binding modes to key macromolecular targets of COVID-19 were identified using molecular docking.

RESULTS

The epigenomic analysis of the relationship between SARS-CoV-2 infection and HCC patients revealed that the co-pathogenesis was closely linked to immune response, particularly T cell differentiation, regulation of T cell activation and monocyte differentiation. Further analysis indicated that CD4⁺ T cells and monocytes play essential roles in the immunoreaction triggered by both conditions. The expression levels of hub genes MYLK2, FAM83D, STC2, CCDC112, EPHX4 and MMP1 were strongly correlated with SARS-CoV-2 infection and the prognosis of HCC patients. In our study, mefloquine and thioridazine were identified as potential therapeutic agents for COVID-19 in combined with HCC.

CONCLUSIONS

In this research, we conducted an epigenomics analysis to identify common pathogenetic processes between SARS-CoV-2 infection and HCC patients, providing new insights into the pathogenesis and treatment of HCC patients infected with SARS-CoV-2.

Plasma Exosome Gene Signature Differentiates Colon Cancer from Healthy Controls

BACKGROUND

Liquid biopsies have become an integral part of cancer management as minimally invasive options to detect molecular and genetic changes. However, current options show poor sensitivity in peritoneal carcinomatosis (PC). Novel exosome-based liquid biopsies may provide critical information on these challenging tumors. In this initial feasibility analysis, we identified an exosome gene signature of 445 genes (ExoSig445) from colon cancer patients, including those with PC, that is distinct from healthy controls.

METHODS

Plasma exosomes from 42 patients with metastatic and non-metastatic colon cancer and 10 healthy controls were isolated and verified. RNAseq analysis of exosomal RNA was performed and differentially expressed genes (DEGs) were identified by the DESeq2 algorithm. The ability of RNA transcripts to discriminate control and cancer cases was assessed by principal component analysis (PCA) and Bayesian compound covariate predictor classification. An exosomal gene signature was compared with tumor expression profiles of The Cancer Genome Atlas.

RESULTS

Unsupervised PCA using exosomal genes with greatest expression variance showed stark separation between controls and patient samples. Using separate training and test sets, gene classifiers were constructed capable of discriminating control and patient samples with 100% accuracy. Using a stringent statistical threshold, 445 DEGs fully delineated control from cancer samples. Furthermore, 58 of these exosomal DEGs were found to be overexpressed in colon tumors.

CONCLUSIONS

Plasma exosomal RNAs can robustly discriminate colon cancer patients, including patients with PC, from healthy controls. ExoSig445 can potentially be developed as a highly sensitive liquid biopsy test in colon cancer.

Retinal determination gene networks: from biological functions to therapeutic strategies

The retinal determinant gene network (RDGN), originally discovered as a critical determinator in Drosophila eye specification, has become an important regulatory network in tumorigenesis and progression, as well as organogenesis. This network is not only associated with malignant biological behaviors of tumors, such as proliferation, and invasion, but also regulates the development of multiple mammalian organs. Three members of this conservative network have been extensively investigated, including DACH, SIX, and EYA. Dysregulated RDGN signaling is associated with the initiation and progression of tumors. In recent years, it has been found that the members of this network can be used as prognostic markers for cancer patients. Moreover, they are considered to be potential therapeutic targets for cancer. Here, we summarize the research progress of RDGN members from biological functions to signaling transduction, especially emphasizing their effects on tumors. Additionally, we discuss the roles of RDGN members in the development of organs and tissue as well as their correlations with the pathogenesis of chronic kidney disease and coronary heart disease. By summarizing the roles of RDGN members in human diseases, we hope to promote future investigations into RDGN and provide potential therapeutic strategies for patients.

Silencing of SiX-4 enhances the chemosensitivity of melanoma cells to Cisplatin

Melanoma is the riskiest type of skin cancer. Its prevalence has been rapidly increased over the last three decades. SIX1, SIX2, SIX3, SIX4, SIX5, and SIX6 are members of the sine oculis homeobox (SIX) homolog family. It is imperative to identify new melanoma biomarkers to improve the predictive value for melanoma prognosis, which could enhance our understanding of carcinogenesis and tumor progression. In this study, we investigated whether silencing of SIX4 in a melanoma cell line (A375 cells) in combination with Cisplatin can affect the apoptosis and suppression of cell cycle progression, migration of the melanoma cells. MTT test and colony formation assay was applied to determine the IC50 of Cisplatin and the combined effect of SIX4 siRNA and Cisplatin on the viability and clonogenesis of the A-375 cells. qRT-PCR was performed to determine the c-myc, BCL-2, BAX, MMP-9, CXCR4, and Rock genes expression. Furthermore, flow cytometry was applied to evaluate apoptosis, autophagy, and the cell cycle status in different groups. Finally, wound healing assay was employed to evaluate the effect of this combination therapy on migratory capacity. SIX4 suppression increased the chemosensitivity of A-375 cells to Cisplatin and decreased its efficient dose. Furthermore, SIX4 suppression alongside Cisplatin reduced cell migration rate, arrested the cell cycle at the G1 phase, induced apoptosis by modulating the expression of apoptotic target genes, induced autophagy, and also significantly inhibits clonogenesis of A-375 cells. SIX4 plays a significant role in the chemosensitivity and pathogenesis of melanoma. Therefore, SIX4 suppression, in combination with Cisplatin, may be a promising therapeutic approach in treating melanoma.

Sine oculis homeobox homolog 1 plays a critical role in pulmonary fibrosis

Idiopathic pulmonary fibrosis (IPF) is a fatal disease with limited treatment options. The role of the developmental transcription factor Sine oculis homeobox homolog 1 (SIX1) in the pathophysiology of lung fibrosis is not known. IPF lung tissue samples and IPF-derived alveolar type II cells (AT2) showed a significant increase in SIX1 mRNA and protein levels, and the SIX1 transcriptional coactivators EYA1 and EYA2 were elevated. Six1 was also upregulated in bleomycin-treated (BLM-treated) mice and in a model of spontaneous lung fibrosis driven by deletion of Telomeric Repeat Binding Factor 1 (Trf1) in AT2 cells. Conditional deletion of Six1 in AT2 cells prevented or halted BLM-induced lung fibrosis, as measured by a significant reduction in histological burden of fibrosis, reduced fibrotic mediator expression, and improved lung function. These effects were associated with increased macrophage migration inhibitory factor (MIF) in lung epithelial cells in vivo following SIX1 overexpression in BLM-induced fibrosis. A MIF promoter-driven luciferase assay demonstrated direct binding of Six1 to the 5'-TCAGG-3' consensus sequence of the MIF promoter, identifying a likely mechanism of SIX1-driven MIF expression in the pathogenesis of lung fibrosis and providing a potentially novel pathway for targeting in IPF therapy.

FOXA1 Leads to Aberrant Expression of SIX4 Affecting Cervical Cancer Cell Growth and Chemoresistance

Cervical cancer (CC) is among the most prevalent cancers among female populations with high recurrence rates all over the world. Cisplatin (DDP) is the first-line treatment for multiple cancers, including CC. The main problem associated with its clinical application is drug resistance. This study is aimed at investigating the function and downstream regulation mechanism of forkhead-box A1 (FOXA1) in CC, which was verified as an oncogene in several cancers. Using GEO database and bioinformatics analysis, we identified FOXA1 as a possible oncogene in CC. Silencing of FOXA1 inhibited CC cell growth, invasion, and chemoresistance. Afterwards, the downstream gene of FOXA1 was predicted using a bioinformatics website and validated using ChIP and dual-luciferase assays. SIX4, a possible target of FOXA1, promoted CC cell malignant aggressiveness and chemoresistance. In addition, overexpression of SIX4 promoted phosphorylation of PI3K and AKT proteins and activated the PI3K/AKT signaling pathway. Further overexpression of SIX4 reversed the repressive effects of FOXA1 knockdown on CC cell growth, invasion, and chemoresistance in DDP-resistant cells. FOXA1-induced SIX4 facilitates CC progression and chemoresistance, highlighting a strong potential for FOXA1 to serve as a promising therapeutic target in CC.

NIK-SIX1 signalling axis regulates high glucose-induced endothelial cell dysfunction and inflammation

Endothelial dysfunction and inflammation are the main manifestations of diabetes-associated atherosclerosis. This paper studied the roles of NF-κB-inducing kinase (NIK) and sine oculis homeobox homolog 1 (SIX1) in regulating high glucose-induced endothelial dysfunction and inflammation. The expression of NIK and SIX1 in human umbilical vein endothelial cells (HUVECs) was silenced by transfection with the specific shRNAs. HUVECs exposed to high glucose were considered as a cell model of endothelial dysfunction. Expression of NIK and SIX1 following transfection was measured by qRT-PCR and western blotting analysis. The proliferation, migration, and inflammation of HUVECs were evaluated by EdU staining, scratch test, ELISA, and western blotting. High glucose (30 mM) significantly decreased the proliferation and migration of HUVECs. High glucose-induced the expression of adhesion molecules VCAM-1 and ICAM-1. Moreover, high glucose increased the release of IL-1β, IL-6, TNF-α, and MCP-1. Transfection of cells with NIK shRNA significantly reversed the toxic effects of high glucose on HUVECs. Of contrast, SIX1 shRNA accelerated the effects of high glucose on HUVECs. NIK shRNA inhibited the accumulation of RelA, RelB, and p52. Meanwhile, NIK shRNA led to SIX1 downregulation which further induced the activation of the NF-κB pathway. NIK-SIX1 signalling axis was suggested to be critical in the regulation of high glucose-induced endothelial dysfunction and inflammation. SIX1 may function as an immunological gatekeeper to control the excessive inflammation mediated by NIK in diabetes-associated atherosclerosis.

NFAT as a Biomarker and Therapeutic Target in Non-Small Cell Lung Cancer-Related Brain Metastasis

Background

Brain metastases (BMs) are associated with poor prognosis and significant mortality, and approximately 25% of patients with non-small cell lung cancer (NSCLC) develop BMs. The present study was aimed to understand the relationships between BM and NSCLC and reveal potential biomarkers and therapeutic targets in NSCLC-related BM.

Methods

The differentially expressed genes (DEGs) expressed during NSCLC and BM development were predicted by bioinformatics analysis, and the expression of the upstream transcription factor nuclear factor of activated T cells (NFAT) was confirmed as a differential gene expressed in both NSCLC and BM. In addition, the expression of proteins encoded by these DEGs was verified by immunohistochemical experiments examining normal lung tissue, lung cancer tissue, and brain metastasis tissue from 30 patients with NSCLC related BM.

Results

The co-DEGs interleukin (IL)-11, cadherin 5 (CDH5) and C-C motif chemokine 2 (CCL2) link NSCLC and BM in the Gene Expression Omnibus (GEO) database, and NFAT may target the expression of these co-DEGs. In the GEO database, NFATc1 and NFATc3 were significantly downregulated in NSCLC tissues (P <.05), whereas NFATc1, NFATc2, NFATc3, and NFATc4 were significantly downregulated in BMs (P <.05). Consistent findings were observed in the immunohistochemical analysis.

Conclusion

NFATc1 and NFATc3 may play important roles in the occurrence of NSCLC and BM by regulating IL-11, CDH5, and CCL2.

A promising prognostic signature for lung adenocarcinoma (LUAD) patients basing on 6 hypoxia-related genes

BACKGROUND

Hypoxia signaling plays a critical role in the development of lung adenocarcinoma (LUAD). We herein aimed to explore the prognostic value of hypoxia-related genes and construct the hypoxia-related prognostic signature for LUAD patients.

METHODS

A total of 26 hypoxia-related genes were collected. Five hundred thirteen and 246 LUAD samples were obtained from the Cancer Genome Atlas and Gene Expression Omnibus databases, respectively. Univariate Cox regression and LASSO Cox regression analyses were conducted to screen the hypoxia-related genes associated with the prognosis of LUAD patients, which would be used for constructing prognosis predictive model for LUAD patients. Multivariate Cox regression analysis was done to determine the independent prognostic factors. The Nomogram model was constructed to predict the prognosis of LUAD patients.

RESULTS

Based on 26 hypoxia-related genes, LUAD samples could be divided into 4 clusters with different prognoses. Among which, 6 genes were included to construct the Risk Score and the LUAD patients with higher Risk Score had worse prognosis. Besides, the Nomogram based on all the independent risk factors could relatively reliably predict the survival probability. And 9 types of immune cells' infiltration was significantly differential between high and low risk LUAD patients.

CONCLUSION

The Risk Score model based on the 6 crucial hypoxia-related genes could relatively reliably predict the prognosis of LUAD patients.

Homeoprotein SIX1 compromises antitumor immunity through TGF-β-mediated regulation of collagens

The tumor microenvironment (TME), including infiltrated immune cells, is known to play an important role in tumor growth; however, the mechanisms underlying tumor immunogenicity have not been fully elucidated. Here, we discovered an unexpected role for the transcription factor SIX1 in regulating the tumor immune microenvironment. Based on analyses of patient datasets, we found that SIX1 was upregulated in human tumor tissues and that its expression levels were negatively correlated with immune cell infiltration in the TME and the overall survival rates of cancer patients. Deletion of Six1 in cancer cells significantly reduced tumor growth in an immune-dependent manner with enhanced antitumor immunity in the TME. Mechanistically, SIX1 was required for the expression of multiple collagen genes via the TGFBR2-dependent Smad2/3 activation pathway, and collagen deposition in the TME hampered immune cell infiltration and activation. Thus, our study uncovers a crucial role for SIX1 in modulating tumor immunogenicity and provides proof-of-concept evidence for targeting SIX1 in cancer immunotherapy.

The Significance of SIX1 as a Prognostic Biomarker for Survival Outcome in Various Cancer Patients: A Systematic Review and Meta-Analysis

Sine Oculis Homeobox Homolog 1 (SIX1) is reported to promote cancer initiation and progression in many preclinical models and is demonstrated in human cancer tissues. However, the correlation between SIX1 and cancer patients’ prognosis has not yet been systematically evaluated. Therefore, we performed a systematic review and meta-analysis in various human cancer types and extracted some data from TCGA datasets for further verification and perfection. We constructed 27 studies and estimated the association between SIX1 expression in various cancer patients’ overall survival and verified with TCGA datasets. Twenty-seven studies with 4899 patients are include in the analysis of overall, and disease-free survival, most of them were retrospective. The pooled hazard ratios (HRs) for overall and disease-free survival in high SIX1 expression patients were 1.54 (95% CI: 1.32-1.80, P<0.00001) and 1.83 (95% CI: 1.31-2.55, P=0.0004) respectively. On subgroup analysis classified in cancer type, high SIX1 expression was associated with poor overall survival in patients with hepatocellular carcinoma (HR 1.50; 95% CI: 1.17-1.93, P =0.001), breast cancer (HR 1.31; 95% CI: 1.10-1.55, P =0.002) and esophageal squamous cell carcinoma (HR 1.89; 95% CI: 1.42-2.52, P<0.0001). Next, we utilized TCGA online datasets, and the consistent results were verified in various cancer types. SIX1 expression indicated its potential to serve as a cancer biomarker and deliver prognostic information in various cancer patients. More works still need to improve the understandings of SIX1 expression and prognosis in different cancer types.

A Novel Seven Gene Signature-Based Prognostic Model to Predict Distant Metastasis of Lymph Node-Negative Triple-Negative Breast Cancer

Background

The prognosis of lymph node-negative triple-negative breast cancer (TNBC) is still worse than that of other subtypes despite adjuvant chemotherapy. Reliable prognostic biomarkers are required to identify lymph node-negative TNBC patients at a high risk of distant metastasis and optimize individual treatment.

Methods

We analyzed the RNA sequencing data of primary tumor tissue and the clinicopathological data of 202 lymph node-negative TNBC patients. The cohort was randomly divided into training and validation sets. Least absolute shrinkage and selection operator Cox regression and multivariate Cox regression were used to construct the prognostic model.

Results

A clinical prognostic model, seven-gene signature, and combined model were constructed using the training set and validated using the validation set. The seven-gene signature was established based on the genomic variables associated with distant metastasis after shrinkage correction. The difference in the risk of distant metastasis between the low- and high-risk groups was statistically significant using the seven-gene signature (training set: P < 0.001; validation set: P = 0.039). The combined model showed significance in the training set (P < 0.001) and trended toward significance in the validation set (P = 0.071). The seven-gene signature showed improved prognostic accuracy relative to the clinical signature in the training data (AUC value of 4-year ROC, 0.879 vs. 0.699, P = 0.046). Moreover, the composite clinical and gene signature also showed improved prognostic accuracy relative to the clinical signature (AUC value of 4-year ROC: 0.888 vs. 0.699, P = 0.029; AUC value of 5-year ROC: 0.882 vs. 0.693, P = 0.038). A nomogram model was constructed with the seven-gene signature, patient age, and tumor size.

Conclusions

The proposed signature may improve the risk stratification of lymph node-negative TNBC patients. High-risk lymph node-negative TNBC patients may benefit from treatment escalation.

The SIX Family of Transcription Factors: Common Themes Integrating Developmental and Cancer Biology

The sine oculis (SIX) family of transcription factors are key regulators of developmental processes during embryogenesis. Members of this family control gene expression to promote self-renewal of progenitor cell populations and govern mechanisms of cell differentiation. When the function of SIX genes becomes disrupted, distinct congenital defects develops both in animal models and humans. In addition to the embryonic setting, members of the SIX family have been found to be critical regulators of tumorigenesis, promoting cell proliferation, epithelial-to-mesenchymal transition, and metastasis. Research in both the fields of developmental biology and cancer research have provided an extensive understanding of SIX family transcription factor functions. Here we review recent progress in elucidating the role of SIX family genes in congenital disease as well as in the promotion of cancer. Common themes arise when comparing SIX transcription factor function during embryonic and cancer development. We highlight the complementary nature of these two fields and how knowledge in one area can open new aspects of experimentation in the other.

Whole-Genome Methylation Study of Congenital Lung Malformations in Children

Background and Objectives

The treatment of asymptomatic patients with congenital pulmonary malformations (CPMs) remains controversial, partially because the relationship between congenital lung malformations and malignancy is still undefined. Change in methylation pattern is a crucial event in human cancer, including lung cancer. We therefore studied all differentially methylated regions (DMRs) in a series of CPMs in an attempt to find methylation anomalies in genes already described in association with malignancy.

Methods

The DNA extracted from resected congenital lung malformations and control lung tissue was screened using Illumina MethylationEPIC arrays. Comparisons between the group of malformed samples or the malformed samples of same histology or each malformed sample and the controls and between a pleuropulmonary blastoma (PPB) and controls were performed. Moreover, each malformed sample was pairwise compared with its respective control. All differentially methylated regions (DMRs) with an adjusted p-value <0,05 were studied.

Results

Every comparison highlighted a number of DMRs closed to genes involved either in cell proliferation or in embryonic development or included in the Cancer Gene Census. Their abnormal methylation had been already described in lung tumors.

Conclusions

Methylation anomalies already described in lung tumors and also shared by the PPB were found in congenital lung malformations, regardless the histology. The presence of methylation abnormalities is suggestive of a correlation between congenital lung malformations and some step of malignant transformation.

miR-4306 Suppresses Proliferation of Esophageal Squamous Cell Carcinoma Cell by Targeting SIX3

Esophageal squamous cell carcinoma (ESCC) is one of the most aggressive cancers and the primary cause of cancer-related mortality in China. micoRNA plays a vital role during tumor initiation and malignant progression. miR-4306 has been reported to negatively regulate aggressive cell phenotypes in triple-negative breast cancer (TNBC). Nevertheless, the function of miR-4306 in ESCC was still not clear. In this study, we detected miR-4306 expression by quantitative real-time reverse transcription-PCR (qRT-PCR) and found that miR-4306 expression was downregulated in human ESCC tissue samples and cell lines. Moreover, miR-4306 overexpression could restrain ESCC cell proliferation, migratory and invasive ability and epithelial-mesenchymal transition (EMT), promote cell apoptosis after treatment with or without cisplatin. In contrast, inhibiting the expression of miR-4306 showed the opposing results. Furthermore, we explored the molecular mechanism of effects of miR-4306 and found that miR-4306 inhibited the expression of SIX3 by interaction with SIX3 3'UTR in ESCC cells, and SIX3 overexpression significantly reversed the effect of miR-4306-mediated ESCC cells proliferation. The current study provided evidence of miR-4306 as a tumor suppression gene in ESCC.

DNA Methylation Haplotype Block Markers Efficiently Discriminate Follicular Thyroid Carcinoma from Follicular Adenoma

CONTEXT

Follicular thyroid carcinoma (FTC) is the second most common type of thyroid carcinoma and must be pathologically distinguished from benign follicular adenoma (FA). Additionally, the clinical assessment of thyroid tumors with uncertain malignant potential (TT-UMP) demands effective indicators.

OBJECTIVE

We aimed to identify discriminating DNA methylation markers between FA and FTC.

METHODS

DNA methylation patterns were investigated in 33 FTC and 33 FA samples using reduced representation bisulfite sequencing and methylation haplotype block-based analysis. A prediction model was constructed and validated in an independent cohort of 13 FTC and 13 FA samples. Moreover, 36 TT-UMP samples were assessed using this model.

RESULTS

A total of 70 DNA methylation markers, approximately half of which were located within promoters, were identified to be significantly different between the FTC and FA samples. All the Gene Ontology terms enriched among the marker-associated genes were related to "DNA binding," implying that the inactivation of DNA binding played a role in FTC development. A random forest model with an area under the curve of 0.994 was constructed using those markers for discriminating FTC from FA in the validation cohort. When the TT-UMP samples were scored using this model, those with fewer driver mutations also exhibited lower scores.

CONCLUSION

An FTC-predicting model was constructed using DNA methylation markers, which distinguished between FA and FTC tissues with a high degree of accuracy. This model can also be used to help determine the potential of malignancy in TT-UMP.

Silencing SIX1 by miR-7160 inhibits non-small cell lung cancer cell growth

The homeoprotein SIX1 is upregulated in non-small cell lung cancer (NSCLC) and associated with NSCLC tumorigenesis and progression. We identified microRNA-7160 (miR-7160) as a SIX1-targeting miRNA. RNA immunoprecipitation results confirmed a direct binding between miR-7160 and SIX1 mRNA in NSCLC cells. In the primary and established NSCLC cells, forced overexpression of miR-7160 downregulated SIX1 and inhibited cancer cell growth, proliferation, migration and invasion. Furthermore, miR-7160 overexpression induced apoptosis activation in NSCLC cells. Conversely, miR-7160 inhibition elevated SIX1 expression and enhanced NSCLC cell progression in vitro. Restoring SIX1 expression, by an untranslated region-depleted SIX1 expression construct, reversed miR-7160-induced anti-NSCLC cell activity. CRISPR/Cas9-inudced knockout of SIX1 mimicked miR-7160-induced actions and produced anti-NSCLC cell activity. In vivo, intratumoral injection of miR-7160-expressing lentivirus downregulated SIX1 mRNA and inhibited NSCLC xenograft growth in severe combined immunodeficient mice. Significantly, miR-7160 expression is downregulated in human NSCLC tissues and is correlated with SIX1 mRNA upregulation. Collectively, miR-7160 silenced SIX1 and inhibited NSCLC cell growth in vitro and in vivo.

VEGF-C mediates tumor growth and metastasis through promoting EMT-epithelial breast cancer cell crosstalk

It is well established that a subset of cells within primary breast cancers can undergo an epithelial-to-mesenchymal transition (EMT), although the role of EMT in metastasis remains controversial. We previously demonstrated that breast cancer cells that had undergone an oncogenic EMT could increase metastasis of neighboring cancer cells via non-canonical paracrine-mediated activation of GLI activity that is dependent on SIX1 expression in the EMT cancer cells. However, the mechanism by which these SIX1-expressing EMT cells activate GLI signaling remained unclear. In this study, we demonstrate a novel mechanism for activation of GLI-mediated signaling in epithelial breast tumor cells via EMT cell-induced production and secretion of VEGF-C. We show that VEGF-C, secreted by breast cancer cells that have undergone an EMT, promotes paracrine-mediated increases in proliferation, migration, and invasion of epithelial breast cancer cells, via non-canonical activation of GLI-signaling. We further show that the aggressive phenotypes, including metastasis, imparted by EMT cells on adjacent epithelial cancer cells can be disrupted by either inhibiting VEGF-C in EMT cells or by knocking down NRP2, a receptor which interacts with VEGF-C, in neighboring epithelial cancer cells. Interrogation of TCGA and GEO public datasets supports the relevance of this pathway in human breast cancer, demonstrating that VEGF-C strongly correlates with activation of Hedgehog signaling and EMT in the human disease. Our study suggests that the VEGF-C/NRP2/GLI axis is a novel and conserved paracrine means by which EMT cells enhance metastasis, and provides potential targets for therapeutic intervention in this heterogeneous disease.

Upregulation of SIX4 indicates poor clinical outcome and promotes tumor growth and cell metastasis in esophageal squamous cell carcinoma

Background

The role of sine oculis homeobox 4 (SIX4) has been found in some malignant tumors. However, there have been few studies on the function of SIX4 in esophageal squamous cell carcinoma (ESCC). This study aimed to explore the regulatory mechanism of SIX4 in ESCC.

Methods

RT‐qPCR and Western blot analysis were used to measure mRNA and protein expression. The function of SIX4 was investigated using CCK‐8, colony formation, flow cytometry, wound healing and transwell assays. A mouse xenograft tumor assay was designed to perform in vivo experiments.

Results

SIX4 was upregulated in ESCC and indicated poor clinical outcomes in ESCC patients. Functionally, knockdown of SIX4 inhibited cell proliferation and induced apoptosis in ESCC. In addition, the silencing of SIX4 inhibited cell migration, invasion and EMT in ESCC. More importantly, upregulation of SIX4 could activate the PI3K/AKT pathway in ESCC cells and promote tumor growth in vivo.

Conclusions

Upregulation of SIX4 indicates poor clinical outcomes in ESCC patients and promotes tumor growth and cell metastasis in ESCC.

TRIM27 acts as an oncogene and regulates cell proliferation and metastasis in non-small cell lung cancer through SIX3-β-catenin signaling

The Wnt/β-catenin pathway plays vital roles in diverse biological processes, including cell differentiation, proliferation, migration, and insulin sensitivity. A recent study reported that the DNA-binding transcriptional factor SIX3 is essential during embryonic development in vertebrates and capable of downregulating target genes of the Wnt/β-catenin pathway in lung cancer, indicating negative regulation of Wnt/β-catenin activation. However, regulation of the SIX3-Wnt/β-catenin pathway axis remains unknown. We measured the expression of TRIM27 and SIX3 as well as investigated whether there was a correlation between them in lung cancer tissue samples. Herein, we found that the E3 ubiquitin ligase, TRIM27, ubiquitinates, and degrades SIX3. TRIM27 induces non-small cell lung cancer (NSCLC) cell proliferation and metastasis, and the expression of β-catenin, S100P, TGFB3, and MMP-9 were significantly inhibited by SIX3. Furthermore, XAV939 is a selective β-catenin-mediated transcription inhibitor that inhibited TRIM27- and SIX3-mediated NSCLC cell proliferation, migration, and invasion. Clinically, lung tissue samples of cancer patients showed increased TRIM27 expression and decreased SIX3 expression. Taken together, these data demonstrate that TRIM27 acts as an oncogene regulating cell proliferation and metastasis in NSCLC through SIX3-β-catenin signaling.

SIX6 is a TAL1-regulated transcription factor in T-ALL and associated with inferior outcome

T-cell acute lymphoblastic leukemia (T-ALL) is a hematological malignancy driven by abnormal activity of transcription factors. Here we report an aberrant expression of the developmental transcription factor SIX6 in the TAL1-subtype of T-ALL. Our results demonstrate that the binding of TAL1 and GATA3 transcription factors into an upstream enhancer element directly regulates SIX6 expression. High expression of SIX6 was associated with inferior event-free survival within three independent patient cohorts. At a functional level, CRISPR-Cas9-mediated knockout of the SIX6 gene in TAL1 positive Jurkat cells induced changes in genes associated with the mTOR-, K-RAS-, and TNFα-related molecular signatures but did not impair cell proliferation or viability. There was also no acceleration of T-ALL development within a Myc driven zebrafish tumor model in vivo. Taken together, our results show that SIX6 belongs to the TAL1 regulatory gene network in T-ALL but is alone insufficient to influence the development or maintenance of T-ALL.

Integrated Analysis of lncRNAs and mRNAs Identifies a Potential Driver lncRNA FENDRR in Lung Cancer in Xuanwei, China

This study was to screen out potential driver long non-coding RNAs (lncRNAs) in lung cancer in Xuanwei (LCXW) differently expressed mRNAs and lncRNAs were detected by gene expression microarrays in 23 paired lung adenocarcinoma and adjacent tissues. Combined bioinformatics analysis was performed to identify potential driver lncRNAs and their potential regulatory relationships. Transcriptome and clinical data in TCGA-LUAD were used as comparison and validation dataset. The comparison of LCXW and TCGA-LUAD revealed significant differences in expression of some genes, signaling pathways affected by differentially expressed genes, and the 5-year survival rate of patients. We identified 14 consistently deregulated mRNAs and 5 lncRNAs as candidate genes, which affected multiple cancer-related pathways and influenced patients' overall survival. By combined bioinformatics analysis, we further identified a potential driver lncRNA fetal-lethal non-coding developmental regulatory RNA (FENDRR) and proposed its possible regulation mechanism. The low expression of FENDRR was positively correlated with Krüppel-like factor4 (KLF4), KLF4 down-regulation may loss the activation function of cyclin-dependent kinase inhibitor 1A (CDKN1A) and cyclin-dependent kinase inhibitor 1C (CDKN1C) and the inhibition function of CyclinB1 (CCNB1), eventually cause excessive cell cycle activation and lead to lung cancer. This study revealed a potential FENDRR-KLF4-cell cycle regulation axis. These results lay an important foundation for further research on the pathogenesis of LCXW and identification of potential novel biomarkers or therapeutic targets.

RDGN-based predictive model for the prognosis of breast cancer

Background

Breast cancer is the most diagnosed malignancy in females in the United States. The members of retinal determination gene network (RDGN) including DACH, EYA, as well as SIX families participate in the proliferation, apoptosis, and metastasis of multiple tumors including breast cancer. A comprehensive predictive model of RDGN might be helpful to herald the prognosis of breast cancer patients.

Methods

In this study, the Gene Expression Ominibus (GEO) and Gene Set Expression Analysis (GSEA) algorithm were used to investigate the effect of RDGN members on downstream signaling pathways. Besides, based on The Cancer Genome Atlas (TCGA) database, we explored the expression patterns of RDGN members in tumors, normal tissues, and different breast cancer subtypes. Moreover, we estimated the relationship between RDGN members and the outcomes of breast cancer patients. Lastly, we constructed a RDGN-based predictive model by Cox proportional hazard regression and verified the model in two separate GEO datasets.

Results

The results of GSEA showed that the expression of DACH1 was negatively correlated with cell cycle and DNA replication pathways. On the contrary, the levels of EYA2 and SIX1 were significantly positively correlated with DNA replication, mTOR, and Wnt pathways. Further investigation in TCGA database indicated that DACH1 expression was lower in breast cancers especially basal-like subtype. In the meanwhile, SIX1 was remarkably upregulated in breast cancers while EYA2 level was increased in Basal-like and Her-2 enriched subtypes. Survival analyses demonstrated that DACH1 was a favorable factor while EYA2 and SIX1 were risk factors for breast cancer patients. Given the results of Cox proportional hazard regression analysis, two members of RDGN were involved in the present predictive model and patients with high model index had poorer outcomes.

Conclusion

This study showed that aberrant RDGN expression was an unfavorable factor for breast cancer. This RDGN-based comprehensively framework was meaningful for predicting the prognosis of breast cancer patients.

Autoantibody profiling identifies predictive biomarkers of response to anti-PD1 therapy in cancer patients

Background: Programmed cell death protein 1 (PD1) inhibitors have revolutionized cancer therapy, yet many patients fail to respond. Thus, the identification of accurate predictive biomarkers of therapy response will improve the clinical benefit of anti-PD1 therapy.

Method: We assessed the baseline serological autoantibody (AAb) profile against ~2300 proteins in 10 samples and ~4600 proteins in 35 samples with alveolar soft part sarcoma (ASPS), non-small-cell lung cancer (NSCLC) and lymphoma using Nucleic Acid Programmable Protein Arrays (NAPPA). 23 selected potential AAb biomarkers were verified using simple, affordable and rapid enzyme linked immune sorbent assay (ELISA) technology with baseline plasma samples from 12 ASPS, 16 NSCLC and 46 lymphoma patients. SIX2 and EIF4E2 AAbs were further validated in independent cohorts of 17 NSCLC and 43 lymphoma patients, respectively, using ELISA. The IgG subtypes in response to therapy were also investigated.

Results: Distinct AAb profiles between ASPS, NSCLC and lymphoma were observed. In ASPS, the production of P53 and PD1 AAbs were significantly increased in non-responders (p=0.037). In NSCLC, the SIX2 AAb was predictive of response with area under the curve (AUC) of 0.87, 0.85 and 0.90 at 3 months, 4.5 months, 6 months evaluation time points, respectively. In the validation cohort, the SIX2 AAb was consistently up-regulated in non-responders (p=0.024). For lymphoma, the EIF4E2 AAb correlated with a favorable response with AUCs of 0.68, 0.70, and 0.70 at 3 months, 4.5 months, and 6 months, respectively. In the validation cohort, the AUCs were 0.74, 0.75 and 0.66 at 3 months, 4.5 months, and 6 months, respectively. The PD1 and PD-L1 IgG2 AAbs were highly produced in ~20% of lymphoma responders. Furthermore, bioinformatics analysis revealed antigen functions of these AAb biomarkers.

Conclusion: This study provides the first evidence that AAb biomarkers selected using high-throughput protein microarrays can predict anti-PD1 therapeutic response and guide anti-PD1 therapy.

CXCL1 as an Unfavorable Prognosis Factor Negatively Regulated by DACH1 in Non-small Cell Lung Cancer

Background: Interaction between cancer cells with microenvironment is essential for cancer progression, therapeutic resistance and prognosis. Chemokine CXCL1 shows variable roles in the development of cancers. DACH1 has been considered as a tumor suppressor and represses the expressions of several chemokines. The relationship between CXCL1 and DACH1 in non-small cell lung cancer (SCLC) deserves further investigation.

Methods: Immunohistochemistry staining was performed on tumor tissue microarrays from lung cancer patients to detect CXCL1 protein. The CXCL1 concentration in the serum of adenocarcinoma patients was measured by ELISA. The CXCL1 protein secreted by cancer cell lines was detected by SearchLight proteome array and human cytokine antibody array. The meta-analysis of CXCL1 expression form public databases was performed and correlation between CXCL1 and DACH1 was analyzed. Moreover, the association between clinicopathological features and prognosis with CXCL1 and DACH1 was analyzed by tissue array and KM-plotter from public database.

Results: The protein abundance of CXCL1 in lung cancer tissues was significantly higher than that in adjacent normal tissues. CXCL1 was closely related to TNM stage, tumor size, and lymph node metastasis and predicted worse overall survival in adenocarcinoma. The level of CXCL1 in the peripheral blood of adenocarcinoma patients also significantly elevated and positively related with clinical stage. The meta-analysis demonstrated that CXCL1 mRNA level was increased in lung cancer tissues and high level of CXCL1 indicated tumor progression in lung adenocarcinoma. In addition, public database analyses showed that CXCL1 negatively correlated with DACH1. Stable overexpressing DACH1 in cultured lung cancer cells remarkably decreased CXCL1 protein. Moreover, ectopic expression of DACH1 significantly inhibited the expression of CXCL1, Ki67, and cyclin D1 in tumor tissues compared with A549 cells with empty vector. Survival analysis showed that high CXCL1 and low DACH1 indicated poor overall survival and progression-free survival.

Conclusion: CXCL1 is closely associated with tumor progression and poor survival. DACH1 significantly inhibits the expression of CXCL1 and indicates good prognosis. Therefore, combined detection of CXCL1 and DACH1 could more precisely predict prognosis of lung adenocarcinoma.

SIX4 promotes metastasis through STAT3 activation in breast cancer

Sine oculis homeobox homolog 4 (SIX4), a member of the SIX family, play important role in the development and construction of vertebrate tissues and organs. There is very little known about the function of SIX4 in cancer cells. Herein, we investigated whether SIX4 promote cancer metastasis in addition to its direct role in breast cancer cells. Our study showed that the expression of SIX4 was profoundly increased in breast cancer tissues, and the high expression of SIX4 correlated strongly with distant metastasis and poor prognosis. Functional experiments demonstrated that SIX4 obviously promoted the cell migration and invasion of breast cancer cells, and up-regulated the expression of EMT mesenchymal marker, down-regulated the epithelial molecules by Snai1 induction in vitro. Moreover, SIX4 knockdown significant suppressed breast cancer lung metastasis in vivo. Mechanistically, SIX4 directly interacted with STAT3 and promoted the phosphorylated STAT3 nuclear translocation, thereby inducing EMT program activation. Collectively, our findings demonstrated that SIX4 may be a promising target for intervention to prevent cancer metastasis.

Distinct Roles of VEGFA and ANGPT2 in Lung Adenocarcinoma and Squamous Cell Carcinoma

Background: Vascular endothelial growth factor A (VEGFA) and angiopoietin 2 (ANGPT2) are key mediators in angiogenesis. The expression and clinical significance of VEGFA and ANGPT2 have been investigated in lung cancer, but the results are controversial. The specific roles of VEGFA and ANGPT2 in adenocarcinoma (ADC) and squamous cell carcinoma (SQC) are still not fully understood. To characterize it, we conducted the current study.

Materials and Methods: The relationships between clinic-pathological characteristics and the protein expressions of VEGFA and ANGPT2 were analyzed on tissue microarrays by immunohistochemistry (IHC) staining. Then public databases were used to evaluate the association of VEGFA and ANGPT2 mRNA expressions with clinic-pathological parameters and prognosis. Cobalt chloride (CoCl₂) was adopted to mimic a hypoxic microenvironment and western blot was used to detect the expression of hypoxia inducible factor-1α (HIF-1α), VEGFA and ANGPT2 in lung cancer cell lines.

Results: IHC staining revealed that the expressions of VEGFA and ANGPT2 were enriched in lung cancer tissues compared with normal tissues. Additionally, both VEGFA and ANGPT2 protein levels were significantly associated with the tumor size and lymph node metastasis only in ADC, not SQC. More importantly, increased VEGFA and ANGPT2 protein levels were negatively correlated with overall survival (OS) of ADC individuals. Meta-analyses of 22 gene expression omnibus (GEO) databases of lung cancer implicated that patients with higher VEGFA and ANGPT2 mRNA expressions tended to have advanced stage in ADC rather than SQC. Kaplan-Meier plot analyses further verified that high levels of VEGFA and ANGPT2 mRNA were associated with poor survival only in ADC. Moreover, the combination of VEGFA and ANGPT2 could more precisely predict prognosis in ADC. In hypoxia-mimicking conditions, induced expression of HIF-1α unregulated VEGFA and ANGPT2 proteins abundance.

Conclusion: Our results showed hypoxia upregulated the protein levels of VEGFA and ANGPT2 in lung cancer cell lines, and the roles of VEGFA and ANGPT2 were distinct in ADC and SQC. Combined detections of VEGFA and ANGPT2 may be valuable prognostic biomarkers for ADC and double block of VEGFA and ANGPT2 may improve therapeutic outcome.

Competing endogenous RNA (ceRNA) hypothetic model based on comprehensive analysis of long non-coding RNA expression in lung adenocarcinoma

Background

Non-small cell lung cancer (NSCLC) is a major subtype of lung cancer with high malignancy and bad prognosis, consisted of lung adenocarcinomas (LUAD) and lung squamous cell carcinomas (LUSC) chiefly. Multiple studies have indicated that competing endogenous RNA (ceRNA) network centered long noncoding RNAs (lncRNAs) can regulate gene expression and the progression of various cancers. However, the research about lncRNAs-mediated ceRNA network in LUAD is still lacking.

Methods

In this study, we analyzed the RNA-seq database from The Cancer Genome Atlas (TCGA) and obtained dysregulated lncRNAs in NSCLC, then further identified survival associated lncRNAs through Kaplan–Meier analysis. Quantitative real time PCR (qRT-PCR) was performed to confirm their expression in LUAD tissues and cell lines. The ceRNA networks were constructed based on DIANA-TarBase and TargetScan databases and visualized with OmicShare tools. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed to investigate the potential function of ceRNA networks.

Results

In total, 1,437 and 1,699 lncRNAs were found to be up-regulated in LUAD and LUSC respectively with 895 lncRNAs overlapping (|log2FC| > 3, adjusted P value <0.01). Among which, 222 lncRNAs and 46 lncRNAs were associated with the overall survival (OS) of LUAD and LUSC, and 18 out of 222 up-regulated lncRNAs were found to have inverse correlation with LUAD patients’ OS (|log2FC| > 3, adjusted P value < 0.02). We selected 3 lncRNAs (CASC8, LINC01842 and VPS9D1-AS1) out of these 18 lncRNAs and confirmed their overexpression in lung cancer tissues and cells. CeRNA networks were further constructed centered CASC8, LINC01842 and VPS9D1-AS1 with 3 miRNAs and 100 mRNAs included respectively.

Conclusion

Through comprehensively analyses of TCGA, our study identified specific lncRNAs as candidate diagnostic and prognostic biomarkers for LUAD. The novel ceRNA network we created provided more insights into the regulatory mechanisms underlying LUAD.

A New Switch for TGFβ in Cancer

The TGFβ cytokine plays dichotomous roles during tumor progression. In normal and premalignant cancer cells, the TGFβ signaling pathway inhibits proliferation and promotes cell-cycle arrest and apoptosis. However, the activation of this pathway in late-stage cancer cells could facilitate the epithelial-to-mesenchymal transition, stemness, and mobile features to enhance tumorigenesis and metastasis. The opposite functions of TGFβ signaling during tumor progression make it a challenging target to develop anticancer interventions. Nevertheless, the recent discovery of cellular contextual determinants, especially the binding partners of the transcription modulators Smads, is critical to switch TGFβ responses from proapoptosis to prometastasis. In this review, we summarize the recently identified contextual determinants (such as PSPC1, KLF5, 14-3-3ζ, C/EBPβ, and others) and the mechanisms of how tumor cells manage the context-dependent autonomous TGFβ responses to potentiate tumor progression. With the altered expression of some contextual determinants and their effectors during tumor progression, the aberrant molecular prometastatic switch might serve as a new class of theranostic targets for developing anticancer strategies.

Tumor-suppressive activity of Hnf1β in Wilms' tumor

Hepatocyte nuclear factor 1β (HNF1β) is a transcription factor belonging to the HNF-1 family and has been implicated in a number of cancers, but its role in Wilms' tumor (nephroblastoma) has not been addressed. Here, we compared its expression between Wilms' tumor patient kidney tissue and adjacent tissue based on the Oncomine database ( www.oncomine.com ). Cell proliferation, apoptosis, migration, and HNF1β expression level were analyzed in Wilms' tumor-derived G401 cells. Using a variety of mouse tissues (lung, heart, kidney, etc.), we found that HNF1β is the highest expression in the kidneys. Oncomine analysis further demonstrated that HNF1β has a lower expression in Wilms' tumor tissue than in paracancerous tissues. Overexpression of HNF1β decreased cell proliferation and migration, but promoted cell apoptosis. Knockdown of HNF1β produced the opposite results. These results indicated that HNF1β may play important roles in kidney development and function, and its activation may negatively regulate Wilms' tumor progression.

SNPs in LncRNA genes are associated with non-small cell lung cancer in a Chinese population

Background

It has indicated that single nuclear polymorphisms (SNPs) in the regions encoding non‐coding transcripts are associated with lung cancer susceptibility. In a previous microarray study, we identified 13 differentially expressed long non‐coding RNAs (lncRNAs) in non‐small cell lung cancer (NSCLC) and associations of SNPs in these lncRNA genes with lung cancer were unknown. We conducted a case‐control study to address this issue.

Methods

Using the TaqMan method, we genotyped 17 SNPs located in the 13 lncRNA genes in 1294 cases with NSCLC and 1729 healthy controls. Unconditional logistic regression and Cox proportional hazards regression were used to analyze the associations of these SNPs with NSCLC risk and patient survival, respectively. These analyses were also repeated in subgroups of cases and controls stratified by gender, age group, smoking status, disease stage, and histological type.

Results

We identified three SNPs associated with NSCLC risk. For SNP rs498238, CC genotype was associated with lower risk compared to TT genotype (adjusted OR = 0.33, 95%CI: 0.11‐0.97, P = 0.043). For rs16901995, CT/TT genotypes were associated with lower risk compared to CC genotype in non‐smokers (adjusted OR = 0.78, 95%CI: 0.62‐0.98, P = 0.035). Variant genotypes in rs219741 were associated with NSCLC risk in young patients, and the adjusted OR was 1.47 (95%CI: 1.03‐2.10, P = 0.033) when compared to the wild genotype. No SNPs were found to be associated with patient overall survival in the study.

Conclusion

The study suggests that some genetic polymorphisms in the lncRNA genes may influence the risk of NSCLC among Chinese.

Transcriptional factor six2 promotes the competitive endogenous RNA network between CYP4Z1 and pseudogene CYP4Z2P responsible for maintaining the stemness of breast cancer cells

Background

The expression of CYP4Z1 and the pseudogene CYP4Z2P has been shown to be specifically increased in breast cancer by our group and others. Additionally, we previously revealed the roles of the competitive endogenous RNA (ceRNA) network mediated by these genes (ceRNET_CC) in breast cancer angiogenesis, apoptosis, and tamoxifen resistance. However, the roles of ceRNET_CC in regulating the stemness of breast cancer cells and the mechanisms through which ceRNET_CC is regulated remain unclear.

Methods

Transcriptional factor six2, CYP4Z1-3′UTR, and CYP4Z2P-3′UTR were stably overexpressed or knocked down in breast cancer cells via lentivirus infection. ChIP-sequencing and RNA-sequencing analysis were performed to reveal the mechanism through which ceRNET_CC is regulated and the transcriptome change mediated by ceRNET_CC. Clinical samples were used to validate the correlation between six2 and ceRNET_CC. Finally, the effects of the six2/ceRNET_CC axis on the stemness of breast cancer cells and chemotherapy sensitivity were evaluated by in vitro and in vivo experiments.

Results

We revealed that ceRNET_CC promoted the stemness of breast cancer cells. Mechanistically, six2 activated ceRNET_CC by directly binding to their promoters, thus activating the downstream PI3K/Akt and ERK1/2 pathways. Finally, we demonstrated that the six2/ceRNET_CC axis was involved in chemoresistance.

Conclusions

Our results uncover the mechanism through which ceRNET_CC is regulated, identify novel roles for the six2/ceRNET_CC axis in regulating the stemness of breast cancer cells, and propose the possibility of targeting the six2/ceRNET_CC axis to inhibit breast cancer stem cell (CSC) traits.

Electronic supplementary material

The online version of this article (10.1186/s13045-019-0697-6) contains supplementary material, which is available to authorized users.

Encapsulated miR-200c and Nkx2.1 in a nuclear/mitochondria transcriptional regulatory network of non-metastatic and metastatic lung cancer cells

BACKGROUND

MicroRNAs are noncoding RNA molecules of ~ 22 nucleotides with diagnostic and therapeutic action [Curr Drug Targets, 2015. 16(12): p. 1381-403], affecting the expression of mRNAs involved in invasion, migration, and development [Oncotarget, 2015. 6(9): p. 6472-98, Cancer Manag Res, 2014. 6: p. 205-16]. miR-200c is part of the miR-200c/141 cluster on chromosome 12p13. Its mechanism of action when encapsulated is critical in lung cancer when patients express changes in miRNAs. miR-200c be a potential biomarkers for various lung diseases. As a potential therapy, miR-200c can impacts lives as target lung cancer is a leading cause of death with about 234,000 cases annually, high heterogeneity, complex screening, and a 5-year survival rate of 16% [CA Cancer J Clin, 2016.66(1): p. 7-30]. Encapsulated miR-200c efficiently enhances bioavailability, pharmacokinetics of therapeutics and targeting to cells, improves efficacy and provides potential cure.

METHODS

The functions of miR-200c were determined in non-metastatic KW-634 and metastatic 821-T4 and 821-LN mouse lung cancer cell lines after various Nano vehicle treatments. Viability and cytotoxicity were determined by cell cycle and quantitative real-time PCR analyses were used to quantify levels of miR-200c and its target genes. In situ hybridization was used to visualize patterns of expression of miR-200c and others in the lung and many organs. Next-generation sequencing accession number GSE125000, invasion and migration assays using transwell chambers, and ActivSignal were used to elucidate the activation and inhibition profiles and perform direct expression measurements and modification of cellular components.

RESULTS

Due to their effectiveness as intracellular vesicles transporting miR-200c into, out, and between parts of the cells, miR-200c is encapsulated with cholesterol, an integral part of the biological membranes with very important physical properties of the vehicle. Nano miR-200c showed efficient cellular uptake in KW-634, 821-T4, and 821-LN cells with important changes in gene expression and new isoforms. In KW-634, when treated with encapsulated miR-200c and compare to the non-encapsulated control; miR-29b increased by 5261-fold, and in 821-T4/LN, miR-1247 increased by 150-fold. Conversely, miR-1247 and miR-675 decreased by 348 and 1029.5-fold, respectively. miR-189 decreased by 34-fold in treated 821-T4 cells. A reduction of growth was observed only after 48 h of treatment with Nano miR-200c. Moreover, labeling the vehicle with carboxy-fluorescein showed that the encapsulated particles enter the nucleus and mitochondria. Encapsulated miR-200c by entering the cells, the nucleus and mitochondria, trigger changes in cell cycle phases with 4 up to 12 fold percentage in G2 and S phase respectively compare to miR-200c. Endogenous expression of Nkx2.1, miR-200c, and their targets Myb, Nfib, Six4 and Six1 showed an inverse correlation, as observed in development.

CONCLUSIONS

Little is known about miR-200c involvement in regulatory processes. Nano miR-200c affects invasion and migration mechanisms. The expression of encapsulated miR-200c contributes to the inhibition/activation of Kras, EMT, Hippo, regulatory pathways and blockers of metastasis. Delivery of miR-200c increases the expression of miR-29b, an EMY regulator, and miR-1247, an inhibitor of cancer genes, both tumor suppressors involved in lung metastasis. Encapsulated miR-200c act on different proteins that regulates cell cycle pathways. These findings represent a part of a regulatory network providing new insights towards improvement of therapy.

Transcriptome profiling of caspase-2 deficient EμMyc and Th-MYCN mouse tumors identifies distinct putative roles for caspase-2 in neuronal differentiation and immune signaling

Caspase-2 is a highly conserved cysteine protease with roles in apoptosis and tumor suppression. Our recent findings have also demonstrated that the tumor suppression function of caspase-2 is context specific. In particular, while caspase-2 deficiency augments lymphoma development in the EμMyc mouse model, it leads to delayed neuroblastoma development in Th-MYCN mice. However, it is unclear how caspase-2 mediates these differential outcomes. Here we utilized RNA sequencing to define the transcriptomic changes caused by caspase-2 (Casp2^−/−) deficiency in tumors from EμMyc and Th-MYCN mice. We describe key changes in both lymphoma and neuroblastoma-associated genes and identified differential expression of the EGF-like domain-containing gene, Megf6, in the two tumor types that may contribute to tumor outcome following loss of Casp2. We identified a panel of genes with altered expression in Th-MYCN/Casp2^−/− tumors that are strongly associated with neuroblastoma outcome, with roles in melanogenesis, Wnt and Hippo pathway signaling, that also contribute to neuronal differentiation. In contrast, we found that key changes in gene expression in the EμMyc/Casp2^−/− tumors, are associated with increased immune signaling and T-cell infiltration previously associated with more aggressive lymphoma progression. In addition, Rap1 signaling pathway was uniquely enriched in Casp2 deficient EμMyc tumors. Our findings suggest that Casp2 deficiency augments immune signaling pathways that may be in turn, enhance lymphomagenesis. Overall, our study has identified new genes and pathways that contribute to the caspase-2 tumor suppressor function and highlight distinct roles for caspase-2 in different tissues.

Identification of nine microRNAs as potential biomarkers for lung adenocarcinoma

Lung cancer is a leading global cause of cancer‐related death, and lung adenocarcinoma (LUAD) accounts for ~ 50% of lung cancer. Here, we screened for novel and specific biomarkers of LUAD by searching for differentially expressed mRNAs (DEmRNAs) and microRNAs (DEmiRNAs) in LUAD patient expression data within The Cancer Genome Atlas (TCGA). The identified optimal diagnostic miRNA biomarkers were used to establish classification models (including support vector machine, decision tree, and random forest) to distinguish between LUAD and adjacent tissues. We then predicted the targets of identified optimal diagnostic miRNA biomarkers, functionally annotated these target genes, and performed receiver operating characteristic curve analysis of the respective DEmiRNA biomarkers, their target DEmRNAs, and combinations of DEmiRNA biomarkers. We validated the expression of selected DEmiRNA biomarkers by quantitative real‐time PCR (qRT‐PCR). In all, we identified a total of 13 DEmiRNAs, 2301 DEmRNAs and 232 DEmiRNA–target DEmRNA pairs between LUAD and adjacent tissues and selected nine DEmiRNAs (hsa‐mir‐486‐1, hsa‐mir‐486‐2, hsa‐mir‐153, hsa‐mir‐210, hsa‐mir‐9‐1, hsa‐mir‐9‐2, hsa‐mir‐9‐3, hsa‐mir‐577, and hsa‐mir‐4732) as optimal LUAD‐specific biomarkers with great diagnostic value. The predicted targets of these nine DEmiRNAs were significantly enriched in transcriptional misregulation in cancer and central carbon metabolism. Our qRT‐PCR results were generally consistent with our integrated analysis. In summary, our study identified nine DEmiRNAs that may serve as potential diagnostic biomarkers of LUAD. Functional annotation of their target DEmRNAs may provide information on their roles in LUAD.

Progression and prognostic value of ECT2 in non-small-cell lung cancer and its correlation with PCNA

Purpose

Epithelial cell transforming sequence 2 (ECT2) is a guanine nucleotide exchange factor, which is involved in cell division regulation and cell cycle modulation. Recent evidence indicates that ECT2 is overexpressed in many human cancers. However, the exact prognostic value of ECT2 in lung cancer has not been elucidated.

Patients and methods

In the current study, we performed correlation and prognosis analyses using public databases and conducted immunohistochemical staining in tissue microarrays, using samples from 204 lung cancer patients with survival data.

Results

We found that the expression of ECT2 was markedly increased in lung cancer tissues compared with normal tissues. Moreover, we demonstrated that the expression of ECT2 was related to tumor cell differentiation degree, TNM stage, lymph node metastasis, and prognosis in non-small-cell lung cancer (NSCLC). A correlation analysis indicated that ECT2 levels were significantly correlated with proliferating cell nuclear antigen (PCNA) levels in NSCLC. Furthermore, Kaplan–Meier analyses revealed that high ECT2 expression was associated with unfavorable overall survival (OS) and progression-free survival (PFS) in NSCLC patients.

Conclusion

Taken together, these results indicate that the overexpression of ECT2 contributes to tumor invasion and progression, suggesting that ECT2 is a potential prognostic marker for NSCLC patients.

Molecular Changes of Lung Malignancy in HIV Infection

Malignancy of the lung is a major source of morbidity and mortality in persons with human immunodeficiency virus infection; as the most prevalent non-acquired immunodeficiency syndrome-defining malignancy, it represents an important and growing problem confronting HIV-infected patients. To evaluate the molecular changes of lung malignancy in HIV infection, we analyzed differential gene expression profiles and screened for early detection biomarkers of HIV-associated lung cancer using Affymetrix arrays and IPA analysis. A total of 59 patients were diagnosed with HIV-associated lung cancer from Jan 2010 to May 2018. The primary outcome was a significant difference in survival outcome between stages III-IV (10.46 ± 1.87 months) and I-II (17.66 ± 2.88 months). We identified 758 differentially expressed genes in HIV-associated lung cancer. The expression levels of SIX1 and TFAP2A are specifically increased in HIV-associated lung cancer and are associated with poorly differentiated tumor tissue. We also found decreased ADH1B, INMT and SYNPO2 mRNA levels in HIV lung cancer. A comprehensive network and pathway analysis of the dysregulated genes revealed that these genes were associated with four network functions and six canonical pathways relevant to the development of HIV-associated lung cancer. The molecular changes in lung malignancy may help screen the growing population of HIV patients who have or will develop this malignancy.

miR-409 Inhibits Human Non-Small-Cell Lung Cancer Progression by Directly Targeting SPIN1

Lung cancers, the leading cause of cancer mortality worldwide, are characterized by a high metastatic potential. Growing evidence reveals that Spindlin 1 (SPIN1) is involved in tumor progression and carcinogenesis. However, the role of SPIN1 in non-small-cell lung cancer (NSCLC) and the molecular mechanisms underlying SPIN1 in human NSCLC remain undetermined. Here we examined the function of SPIN1 in human NSCLC and found that the expression of SPIN1 was closely correlated with the overall survival and poor prognosis of NSCLC patients. Aberrant regulation of microRNAs (miRNAs) has an important role in cancer progression. We revealed that miR-409 inhibits the expression of SPIN1 by binding directly to the 3′ UTR of SPIN1 using dual-luciferase reporter assays. Overexpression of miR-409 significantly suppressed cell migration, growth, and proliferation by inhibiting SPIN1 in vitro and in vivo. SPIN1 overexpression in miR-409-transfected NSCLC cells effectively rescued the suppression of cell migration, growth, and proliferation regulated by miR-409. miR-409 regulates the PI3K/AKT (protein kinase B) pathway in NSCLC. Moreover, clinical data showed that NSCLC patients with high levels of miR-409 experienced significantly better survival. miR-409 expression was also negatively associated with SPIN1 expression. Taken together, these findings highlight that the miR-409/SPIN1 axis is a useful pleiotropic regulatory network and could predict the metastatic potential in NSCLC patients early, indicating the possibility that miR-409 and SPIN1 might be attractive prognostic markers for treating NSCLC patients.

DACH1 antagonizes CXCL8 to repress tumorigenesis of lung adenocarcinoma and improve prognosis

BACKGROUND

C-X-C motif ligand 8 (CXCL8), known as a proinflammatory chemokine, exerts multiple effects on the proliferation, invasion, and migration of tumor cells via the autocrine or paracrine manner. Conversely, the human Dachshund homologue 1 (DACH1) is recognized as a tumor suppressor which retards the progression of various cancers. In prostate cancer, it has been demonstrated that DACH1 was negatively correlated with the expression of CXCL8 and able to antagonize the effects of CXCL8 on cellular migration. Herein, we explored the mechanisms by which DACH1 regulated the CXCL8 in non-small cell lung cancer (NSCLC).

METHODS

Public microarray and Kaplan-Meier plotter datasets were analyzed. Blood serum samples from lung adenocarcinoma (ADC) patients were collected for enzyme-linked immunosorbent assay (ELISA) analysis. Immunohistochemical staining was conducted on tissue microarray. Cell lines with stable expression of DACH1 were established, and relative gene expression was measured by Western blot, ELISA, real-time PCR, and human cytokine array. Correspondingly, cell lines transfected with shDACH1 were established, and relative gene expression was measured by real-time PCR and immunofluorescence array. Functional studies were performed by transwell and xenograft mice models. Luciferase reporter gene assay was applied to measure the regulation of DACH1 on CXCL8.

RESULTS

Our study indicated that CXCL8 both at the mRNA and protein level was associated with the high tumor burden of ADC. Correlational analyses in ADC cell lines and ADC tissues showed that DACH1 was inversely correlated with CXCL8. Meanwhile, patients with high DACH1 expression and low CXCL8 expression had prolonged time to death and recurrence. Moreover, we verified the inhibitory effects of DACH1 on CXCL8 both in vitro and in vivo. Mechanism studies proved that DACH1 transcriptionally repressed CXCL8 promoter activity through activator protein-1 (AP-1) and nuclear transcription factor-kappa B (NF-κB) sites.

CONCLUSIONS

Our study proved that CXCL8 acted as an unfavorable factor promoting to tumor progression and poor prognosis of ADC, while DACH1 antagonized CXCL8 to provide a favorable survival of ADC patients. Double detection of DACH1 and CXCL8 may provide a precise information for further evaluating the prognosis of ADC patients.

MAT1 correlates with molecular subtypes and predicts poor survival in breast cancer

Objective

Menage a trois 1 (MAT1) is a targeting subunit of cyclin-dependent kinase-activating kinase and general transcription factor IIH kinase, which modulates cell cycle, transcription and DNA repair. Its dysregulation is responsible for diseases including cancers. To further explore the role of MAT1 in breast cancer, we investigated the pathways in which MAT1 might be involved, the association between MAT1 and molecular subtypes, and the role of MAT1 in clinical outcomes of breast cancer patients.

Methods

We conducted immunohistochemistry staining on tissue microarray and immunofluorescence staining on sections of MAT1 stable breast cancer cells. Also, we performed Kyoto Encyclopedia of Genes and Genomes pathway analysis, correlation analysis and prognosis analysis on public databases.

Results

MAT1 was involved in multiple pathways including normal physiology signaling and disease-related signaling. Furthermore, MAT1 positively correlated with the protein status of estrogen receptor and progesterone receptor, and was enriched in luminal-type and human epidermal growth factor receptor 2-enriched breast cancer in comparison with basal-like subtype at both mRNA and protein levels. Correlation analysis revealed significant association between MAT1 mRNA amount and epithelial markers, mesenchymal markers, cancer stem cell markers, apoptosis markers, transcription markers and oncogenes. Consistently, the results of immunofluorescence stain indicated that MAT1 overexpression enhanced the protein abundance of epidermal growth factor receptor, vimentin, sex determining region Y-box 2 and sine oculis homeobox homolog 1. Importantly, Kaplan-Meier Plotter analysis reflected that MAT1 could serve as a prognostic biomarker predicting worse relapse-free survival and metastasis-free survival.

Conclusions

MAT1 is correlated with molecular subtypes and is associated with unfavorable prognosis for breast cancer patients.

The Homeotic Protein SIX3 Suppresses Carcinogenesis and Metastasis through Recruiting the LSD1/NuRD(MTA3) Complex

The homeodomain transcription factor SIX3 was recently reported to be a negative regulator of the Wnt pathway and has an emerging role in cancer. However, how SIX3 contributes to tumorigenesis and metastasis is poorly understood.

METHODS

We employed affinity purification and mass spectrometry (MS) to identify the proteins physically associated with SIX3. Genome-wide analysis of the SIX3/LSD1/NuRD(MTA3) complex using a chromatin immunoprecipitation-on-chip approach identified a cohort of target genes including WNT1 and FOXC2, which are critically involved in cell proliferation and epithelial-to-mesenchymal transition. Also, we used flow cytometry, growth curve analysis, EdU incorporation assay, colony formation assays, trans-well invasion assays, immunohistochemical staining and in vivo bioluminescence assay to investigate the function of SIX3 in tumorigenesis.

RESULTS

We demonstrate that the SIX3/LSD1/NuRD(MTA3) complex inhibits carcinogenesis in breast cancer cells and suppresses metastasis in breast cancer. SIX3 expression is downregulated in various human cancers and high SIX3 is correlated with improved prognosis.

CONCLUSION

Our study revealed an important mechanistic link between the loss of function of SIX3 and tumor progression, identified a molecular basis for the opposing actions of MTA1 and MTA3, and may provide new potential prognostic indicators and targets for cancer therapy.

The clinical significance of CXCL5 in non-small cell lung cancer

As a CXC-type chemokine, ENA78/CXCL5 is an important attractant for granulocytes by binding to its receptor CXCR2. Recent studies proved that CXCL5/CXCR2 axis plays an oncogenic role in many human cancers. However, the exact clinical significance of CXCL5 in lung cancer has not been well defined. Here, we found that the serum protein expression of CXCL5 was significantly increased in non-small cell lung cancer (NSCLC) compared with that in healthy volunteers. Immunohistochemistry staining revealed that CXCL5 protein was higher in various lung cancer tissues compared with normal tissues. Moreover, CXCL5 expression correlated with histological grade, tumor size, and TNM stage in NSCLC. Elevated CXCL5 protein abundance predicted poor overall survival in adenocarcinoma patients. Further meta-analysis demonstrated that CXCL5 mRNA expression was also positively associated with tumor stage, lymph node metastasis, and worse survival. Kaplan–Meier plot analyses indicated high CXCL5 was associated with short overall survival and progression-free survival. Together, these results indicated that CXCL5 may be a potential biomarker for NSCLC.

DACH1 suppresses breast cancer as a negative regulator of CD44

Dachshund homolog 1 (DACH1), a key cell fate determination factor, contributes to tumorigenesis, invasion, metastasis of human breast neoplasm. However, the exact molecular mechanisms for the anti-tumor roles of DACH1 in breast carcinoma are still lack of extensive understanding. Herein, we utilized immunohistochemistry (IHC) staining and public microarray data analysis showing that DACH1 was higher in normal breast, low-grade and luminal-type cancer in comparison with breast carcinoma, high-grade and basal-like tumors respectively. Additionally, both correlation analysis of public databases of human breast carcinoma and IHC analysis of mice xenograft tumors demonstrated that DACH1 inversely related to cancer stem cells (CSCs) markers, epithelial-mesenchymal transition (EMT) inducers and basal-enriched molecules, while cluster of differentiation 44 (CD44) behaved in an opposite manner. Furthermore, mice transplanted tumor model indicated that breast cancer cells Met-1 with up-regulation of DACH1 were endowed with remarkably reduced potential of tumorigenesis. Importantly, meta-analysis of 19 Gene Expression Omnibus (GEO) databases of breast cancer implicated that patients with higher DACH1 expression had prolonged time to death, recurrence and metastasis, while CD44 was a promising biomarker predicting worse overall survival (OS) and metastasis-free survival (MFS). Collectively, our study indicated that CD44 might be a novel target of DACH1 in breast carcinoma.