Nkx2-1: a novel tumor biomarker of lung cancer

A STAG2-PAXIP1/PAGR1 axis suppresses lung tumorigenesis

The cohesin complex is a critical regulator of gene expression. STAG2 is the most frequently mutated cohesin subunit across several cancer types and is a key tumor suppressor in lung cancer. Here, we coupled somatic CRISPR-Cas9 genome editing and tumor barcoding with an autochthonous oncogenic KRAS-driven lung cancer model and showed that STAG2 is uniquely tumor-suppressive among all core and auxiliary cohesin components. The heterodimeric complex components PAXIP1 and PAGR1 have highly correlated effects with STAG2 in human lung cancer cell lines, are tumor suppressors in vivo, and are epistatic to STAG2 in oncogenic KRAS-driven lung tumorigenesis in vivo. STAG2 inactivation elicits changes in gene expression, chromatin accessibility, and 3D genome conformation that impact the cancer cell state. Gene expression and chromatin accessibility similarities between STAG2- and PAXIP1-deficient neoplastic cells further relate STAG2-cohesin to PAXIP1/PAGR1. These findings reveal a STAG2-PAXIP1/PAGR1 tumor-suppressive axis and uncover novel PAXIP1-dependent and PAXIP1-independent STAG2-cohesin-mediated mechanisms of lung tumor suppression.

Neutrophils Recruited by NKX2-1 Suppression via Activation of CXCLs/CXCR2 Axis Promote Lung Adenocarcinoma Progression

NK2 Homeobox 1 (NKX2-1) is a well-characterized pathological marker that delineates lung adenocarcinoma (LUAD) progression. The advancement of LUAD is influenced by the immune tumor microenvironment through paracrine signaling. However, the involvement of NKX2-1 in modeling the tumor immune microenvironment is still unclear. Here, the downregulation of NKX2-1 is observed in high-grade LUAD. Meanwhile, single-cell RNA sequencing and Visium in situ capturing profiling revealed the recruitment and infiltration of neutrophils in orthotopic syngeneic tumors exhibiting strong cell-cell communication through the activation of CXCLs/CXCR2 signaling. The depletion of NKX2-1 triggered the expression and secretion of CXCL1, CXCL2, CXCL3, and CXCL5 in LUAD cells. Chemokine secretion is analyzed by chemokine array and validated by qRT-PCR. ATAC-seq revealed the restrictive regulation of NKX2-1 on the promoters of CXCL1, CXCL2, and CXCL5 genes. This phenomenon led to increased tumor growth, and conversely, tumor growth decreased when inhibited by the CXCR2 antagonist SB225002. This study unveils how NKX2-1 modulates the infiltration of tumor-promoting neutrophils by inhibiting CXCLs/CXCR2-dependent mechanisms. Hence, targeting CXCR2 in NKX2-1-low tumors is a potential antitumor therapy that may improve LUAD patient outcomes.

Significance of NKX2-1 as a biomarker for clinical prognosis, immune infiltration, and drug therapy in lung squamous cell carcinoma

Background

This study was performed to determine the biological processes in which NKX2-1 is involved and thus its role in the development of lung squamous cell carcinoma (LUSC) toward improving the prognosis and treatment of LUSC.

Methods

Raw RNA sequencing (RNA-seq) data of LUSC from The Cancer Genome Atlas (TCGA) were used in bioinformatics analysis to characterize NKX2-1 expression levels in tumor and normal tissues. Survival analysis of Kaplan-Meier curve, the time-dependent receiver operating characteristic (ROC) curve, and a nomogram were used to analyze the prognosis value of NKX2-1 for LUSC in terms of overall survival (OS) and progression-free survival (PFS). Then, differentially expressed genes (DEGs) were identified, and Kyoto Encyclopedia of Genes and Genomes (KEGG), Gene Ontology (GO), and Gene Set Enrichment Analysis (GSEA) were used to clarify the biological mechanisms potentially involved in the development of LUSC. Moreover, the correlation between the NKX2-1 expression level and tumor mutation burden (TMB), tumor microenvironment (TME), and immune cell infiltration revealed that NKX2-1 participates in the development of LUSC. Finally, we studied the effects of NKX2-1 on drug therapy. To validate the protein and gene expression levels of NKX2-1 in LUSC, we employed immunohistochemistry(IHC) datasets, The Gene Expression Omnibus (GEO) database, and qRT-PCR analysis.

Results

NKX2-1 expression levels were significantly lower in LUSC than in normal lung tissue. It significantly differed in gender, stage and N classification. The survival analysis revealed that high expression of NKX2-1 had shorter OS and PFS in LUSC. The multivariate Cox regression hazard model showed the NKX2-1 expression as an independent prognostic factor. Then, the nomogram predicted LUSC prognosis. There are 51 upregulated DEGs and 49 downregulated DEGs in the NKX2-1 high-level groups. GO, KEGG and GSEA analysis revealed that DEGs were enriched in cell cycle and DNA replication.The TME results show that NKX2-1 expression was positively associated with mast cells resting, neutrophils, monocytes, T cells CD4 memory resting, and M2 macrophages but negatively associated with M1 macrophages. The TMB correlated negatively with NKX2-1 expression. The pharmacotherapy had great sensitivity in the NKX2-1 low-level group, the immunotherapy is no significant difference in the NKX2-1 low-level and high-level groups. The analysis of GEO data demonstrated concurrence with TCGA results. IHC revealed NKX2-1 protein expression in tumor tissues of both LUAD and LUSC. Meanwhile qRT-PCR analysis indicated a significantly lower NKX2-1 expression level in LUSC compared to LUAD. These qRT-PCR findings were consistent with co-expression analysis of NKX2-1.

Conclusion

We conclude that NKX2-1 is a potential biomarker for prognosis and treatment LUSC. A new insights of NKX2-1 in LUSC is still needed further research.

Modeling RET-Rearranged Non-Small Cell Lung Cancer (NSCLC): Generation of Lung Progenitor Cells (LPCs) from Patient-Derived Induced Pluripotent Stem Cells (iPSCs)

REarranged during Transfection (RET) oncogenic rearrangements can occur in 1-2% of lung adenocarcinomas. While RET-driven NSCLC models have been developed using various approaches, no model based on patient-derived induced pluripotent stem cells (iPSCs) has yet been described. Patient-derived iPSCs hold great promise for disease modeling and drug screening. However, generating iPSCs with specific oncogenic drivers, like RET rearrangements, presents challenges due to reprogramming efficiency and genotypic variability within tumors. To address this issue, we aimed to generate lung progenitor cells (LPCs) from patient-derived iPSCs carrying the mutation RETC634Y, commonly associated with medullary thyroid carcinoma. Additionally, we established a RETC634Y knock-in iPSC model to validate the effect of this oncogenic mutation during LPC differentiation. We successfully generated LPCs from RETC634Y iPSCs using a 16-day protocol and detected an overexpression of cancer-associated markers as compared to control iPSCs. Transcriptomic analysis revealed a distinct signature of NSCLC tumor repression, suggesting a lung multilineage lung dedifferentiation, along with an upregulated signature associated with RETC634Y mutation, potentially linked to poor NSCLC prognosis. These findings were validated using the RETC634Y knock-in iPSC model, highlighting key cancerous targets such as PROM2 and C1QTNF6, known to be associated with poor prognostic outcomes. Furthermore, the LPCs derived from RETC634Y iPSCs exhibited a positive response to the RET inhibitor pralsetinib, evidenced by the downregulation of the cancer markers. This study provides a novel patient-derived off-the-shelf iPSC model of RET-driven NSCLC, paving the way for exploring the molecular mechanisms involved in RET-driven NSCLC to study disease progression and to uncover potential therapeutic targets.

ANKRD49 promotes the metastasis of NSCLC via activating JNK-ATF2/c-Jun-MMP-2/9 axis

BACKGROUND

Ankyrin repeat domain 49 (ANKRD49) has been found to be highly expressed in multiple cancer including lung adenocarcinoma (LUAD) and lung squamous carcinoma (LUSC). However, the function of ANKRD49 in the pathogenesis of NSCLC still remains elusive. Previously, ANKRD49 has been demonstrated to promote the invasion and metastasis of A549 cells, a LUAD cell line, via activating the p38-ATF-2-MMP2/MMP9 pathways. Considering the heterogeneity of tumor cells, the function and mechanism of ANKRD49 in NSCLC need more NSCLC-originated cells to clarify.

METHODS

Real-time qPCR was employed to test ANKRD49 expression levels in nine pairs of fresh NSCLC tissues and the corresponding adjacent normal tissues. The function of ANKRD49 was investigated using overexpression and RNA interference assays in lung adenocarcinoma cell line (NCI-H1299) and lung squamous carcinoma cell line (NCI-H1703) through gelatin zymography, cell counting kit-8, colony formation, wound healing, migration and invasion assays mmunoprecipitation was performed to in vitro. Immunoprecipitation was performed to test the interaction of c-Jun and ATF2. Chromatin immunoprecipitation was conducted to assess the transcriptional regulation of ATF2/c-Jun on MMP-2/9. Moreover, the tumorigenicity of ANKRD49 was evaluated in nude mice models and the involved signal molecular was also measured by immunohistochemical method.

RESULTS

We found that the levels of ANKRD49 in cancerous tissues were higher than those in adjacent normal tissues. in vitro assay showed that ANKRD49 promoted the migration and invasion of NCI-H1299 and NCI-H1703 cells via enhancing the levels of MMP-2 and MMP-9. Furthermore, ANKRD49 elevated phosphorylation of JNK and then activated c-Jun and ATF2 which interact in nucleus to promote the binding of ATF2:c-Jun with the promoter MMP-2 or MMP-9. In vivo assay showed that ANKRD49 promoted lung metastasis of injected-NSCLC cells and the high metastatic rate was positively correlated with the high expression of ANKRD49, MMP-2, MMP-9, p-JNK, p-c-Jun and p-ATF2.

CONCLUSION

The present study indicated that ANKRD49 accelerated the invasion and metastasis of NSCLC cells via JNK-mediated transcription activation of c-Jun and ATF2 which regulated the expression of MMP-2/MMP-9. The molecular mechanisms of ANKRD49's function is different from those found in A549 cells. The current study is a supplement and improvement to the previous research.

Visual Exploration of Relationships and Structure in Low-Dimensional Embeddings

In this work, we propose an interactive visual approach for the exploration and formation of structural relationships in embeddings of high-dimensional data. These structural relationships, such as item sequences, associations of items with groups, and hierarchies between groups of items, are defining properties of many real-world datasets. Nevertheless, most existing methods for the visual exploration of embeddings treat these structures as second-class citizens or do not take them into account at all. In our proposed analysis workflow, users explore enriched scatterplots of the embedding, in which relationships between items and/or groups are visually highlighted. The original high-dimensional data for single items, groups of items, or differences between connected items and groups are accessible through additional summary visualizations. We carefully tailored these summary and difference visualizations to the various data types and semantic contexts. During their exploratory analysis, users can externalize their insights by setting up additional groups and relationships between items and/or groups. We demonstrate the utility and potential impact of our approach by means of two use cases and multiple examples from various domains.

[Association of ventricular septal defect with rare variations of the HAND2 gene]

OBJECTIVES

To study the association of ventricular septal defect (VSD) with rare variations in the promoter region of HAND2 gene, as well as related molecular mechanisms.

METHODS

Blood samples were collected from 349 children with VSD and 345 healthy controls. The target fragments were amplified by polymerase chain reaction and sequenced to identify the rare variation sites in the promoter region of the HAND2 gene. Dual-luciferase reporter assay was used to perform a functional analysis of the variation sites. Electrophoretic mobility shift assay (EMSA) was used to investigate related molecular mechanisms. TRANSFAC and JASPAR databases were used to predict transcription factors.

RESULTS

Sequencing revealed that three variation sites (g.173530852A>G, g.173531173A>G, and g.173531213C>G) were only observed in the promoter region of the HAND2 gene in 10 children with VSD, among whom 4 children had only one variation site. The dual-luciferase reporter assay revealed that g.173531213C>G reduced the transcriptional activity of the HAND2 gene promoter. EMSA and transcription factor prediction revealed that g.173531213C>G created a binding site for transcription factor.

CONCLUSIONS

The rare variation, g.173531213C>G, in the promoter region of the HAND2 gene participates in the development and progression of VSD possibly by affecting the binding of transcription factors.

Zebrafish as an emerging tool for drug discovery and development for thyroid diseases

Zebrafish is a useful model for understanding human genetics and diseases and has evolved into a prominent scientific research model. The genetic structure of zebrafish is 70% identical to that of humans. Its small size, low cost, and transparent embryo make it a valuable tool in experimentation. Zebrafish and mammals possess the same molecular mechanism of thyroid organogenesis and development. Thus, thyroid hormone signaling, embryonic development, thyroid-related disorders, and novel genes involved in early thyroid development can all be studied using zebrafish as a model. Here in this review, we emphasize the evolving role of zebrafish as a possible tool for studying the thyroid gland in the context of physiology and pathology. The transcription factors nkx2.1a, pax2a, and hhex which contribute a pivotal role in the differentiation of thyroid primordium are discussed. Further, we have described the role of zebrafish as a model for thyroid cancer, evaluation of defects in thyroid hormone transport, thyroid hormone (TH) metabolism, and as a screening tool to study thyrotoxins. Hence, the present review highlights the role of zebrafish as a novel approach to understand thyroid development and organogenesis.

ANKRD49 promotes the invasion and metastasis of lung adenocarcinoma via a P38/ATF-2 signalling pathway

Lung adenocarcinoma (LUAD) is the most challenging neoplasm to treat in clinical practice. Ankyrin repeat domain 49 protein (ANKRD49) is highly expressed in several carcinomas; however, its pattern of expression and role in LUAD are not known. Tissue microarrays, immunohistochemistry, χ² test, Spearman correlation analysis, Kaplan–Meier, log‐rank test, and Cox's proportional hazard model were used to analyse the clinical cases. The effect of ANKRD49 on the LUAD was investigated using CCK‐8, clonal formation, would healing, transwell assays, and nude mice experiment. Expressions of ANKRD49 and its associated downstream protein molecules were verified by real‐time PCR, Western blot, immunohistochemistry, and/or immunofluorescence analyses. ANKRD49 expression was highly elevated in LUAD. The survival rate and Cox's modelling analysis indicated that there may be an independent prognostic indicator for LUAD patients. We also found that ANKRD49 promoted the invasion and migration in both in in vitro and in vivo assays, through upregulating matrix metalloproteinase (MMP)‐2 and MMP‐9 activities via the P38/ATF‐2 signalling pathway Our findings suggest that ANKRD49 is a latent biomarker for evaluating LUAD prognosis and promotes the metastasis of A549 cells via upregulation of MMP‐2 and MMP‐9 in a P38/ATF‐2 pathway‐dependent manner.

[Eccrine adenoma in a cat]

Tumors originating from eccrine glands are rare findings in dogs and cats. In most cases, the tumors are malignant, while adenomas are only reported anecdotally. In the present case, a one-year-old female, spayed cat was presented with a swelling of the footpad of the right forelimb. Initially, the mass possessed a diameter of 2 cm which progressed to 4 cm within the following two months. At the latter time point the tumor was ulcerated. After surgical removal, histological and immunohistochemical analyses were performed. Histologically, a well demarcated, nodular, multilobular mass was present. The cuboidal to columnar neoplastic cells were arranged in tubular and acinar structures. Tumor cells possessed large, round to oval nuclei with moderately distinct nucleoli. Mitotic figures averaged 0-1 per high power field. Additionally, large areas of chondroid metaplasia were evident. Immunohistochemically, neoplastic cells were positive for pan-cytokeratin AE1/AE3 whereas thyroid transcription factor 1 (TTF1) was not expressed. Based on the histological and immunohistochemical findings an adenoma of the eccrine glands was diagnosed.

Functional and Molecular Characters of Cancer Stem Cells Through Development to Establishment

Cancer stem cells (CSCs) are small subpopulation sharing similar properties like normal stem such as self-renewal and differentiation potential to direct tumor growth. Last few years, scientists considered CSCs as the cause of phenotypic heterogeneity in diverse cancer types. Also, CSCs contribute to cancer metastasis and recurrence. The cellular and molecular regulators influence on the CSCs' phenotype changing their behaviors in different stages of cancer progression. CSC markers play significance roles in cancer diagnosis and characterization. We delineate the cross-talks between CSCs and the tumor microenvironment that supports their intrinsic properties including survival, stemness, quiescence and their cellular and molecular adaptation. An insight into the markers of CSCs specific to organs is described.

NKX6-1 mediates cancer stem-like properties and regulates sonic hedgehog signaling in leiomyosarcoma

Background

Leiomyosarcoma (LMS), the most common soft tissue sarcoma, exhibits heterogeneous and complex genetic karyotypes with severe chromosomal instability and rearrangement and poor prognosis.

Methods

Clinical variables associated with NKX6-1 were obtained from The Cancer Genome Atlas (TCGA). NKX6-1 mRNA expression was examined in 49 human uterine tissues. The in vitro effects of NXK6-1 in LMS cells were determined by reverse transcriptase PCR, western blotting, colony formation, spheroid formation, and cell viability assays. In vivo tumor growth was evaluated in nude mice.

Results

Using The Cancer Genome Atlas (TCGA) and human uterine tissue datasets, we observed that NKX6-1 expression was associated with poor prognosis and malignant potential in LMS. NKX6-1 enhanced in vitro tumor cell aggressiveness via upregulation of cell proliferation and anchorage-independent growth and promoted in vivo tumor growth. Moreover, overexpression and knockdown of NKX6-1 were associated with upregulation and downregulation, respectively, of stem cell transcription factors, including KLF8, MYC, and CD49F, and affected sphere formation, chemoresistance, NOTCH signaling and Sonic hedgehog (SHH) pathways in human sarcoma cells. Importantly, treatment with an SHH inhibitor (RU-SKI 43) but not a NOTCH inhibitor (DAPT) reduced cell survival in NKX6-1-expressing cancer cells, indicating that an SHH inhibitor could be useful in treating LMS. Finally, using the TCGA dataset, we demonstrated that LMS patients with high expression of NKX6-1 and HHAT, an SHH pathway acyltransferase, had poorer survival outcomes compared to those without.

Conclusions

Our findings indicate that NKX6-1 and HHAT play critical roles in the pathogenesis of LMS and could be promising diagnostic and therapeutic targets for LMS patients.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12929-021-00726-6.

Nkx2.5 Functions as a Conditional Tumor Suppressor Gene in Colorectal Cancer Cells via Acting as a Transcriptional Coactivator in p53-Mediated p21 Expression

NK2 homeobox 5 (Nkx2.5), a homeobox-containing transcription factor, is associated with a spectrum of congenital heart diseases. Recently, Nkx2.5 was also found to be differentially expressed in several kinds of tumors. In colorectal cancer (CRC) tissue and cells, hypermethylation of Nkx2.5 was observed. However, the roles of Nkx2.5 in CRC cells have not been fully elucidated. In the present study, we assessed the relationship between Nkx2.5 and CRC by analyzing the expression pattern of Nkx2.5 in CRC samples and the adjacent normal colonic mucosa (NCM) samples, as well as in CRC cell lines. We found higher expression of Nkx2.5 in CRC compared with NCM samples. CRC cell lines with poorer differentiation also had higher expression of Nkx2.5. Although this expression pattern makes Nkx2.5 seem like an oncogene, in vitro and in vivo tumor suppressive effects of Nkx2.5 were detected in HCT116 cells by establishing Nkx2.5-overexpressed CRC cells. However, Nkx2.5 overexpression was incapacitated in SW480 cells. To further assess the mechanism, different expression levels and mutational status of p53 were observed in HCT116 and SW480 cells. The expression of p21^WAF1/CIP1, a downstream antitumor effector of p53, in CRC cells depends on both expression level and mutational status of p53. Overexpressed Nkx2.5 could elevate the expression of p21^WAF1/CIP1 only in CRC cells with wild-type p53 (HCT116), rather than in CRC cells with mutated p53 (SW480). Mechanistically, Nkx2.5 could interact with p53 and increase the transcription of p21^WAF1/CIP1 without affecting the expression of p53. In conclusion, our findings demonstrate that Nkx2.5 could act as a conditional tumor suppressor gene in CRC cells with respect to the mutational status of p53. The tumor suppressive effect of Nkx2.5 could be mediated by its role as a transcriptional coactivator in wild-type p53-mediated p21^WAF1/CIP1 expression.

DeePaN: deep patient graph convolutional network integrating clinico-genomic evidence to stratify lung cancers for immunotherapy

Immuno-oncology (IO) therapies have transformed the therapeutic landscape of non-small cell lung cancer (NSCLC). However, patient responses to IO are variable and influenced by a heterogeneous combination of health, immune, and tumor factors. There is a pressing need to discover the distinct NSCLC subgroups that influence response. We have developed a deep patient graph convolutional network, we call "DeePaN", to discover NSCLC complexity across data modalities impacting IO benefit. DeePaN employs high-dimensional data derived from both real-world evidence (RWE)-based electronic health records (EHRs) and genomics across 1937 IO-treated NSCLC patients. DeePaN demonstrated effectiveness to stratify patients into subgroups with significantly different (P-value of 2.2 × 10-11) overall median survival of 20.35 months and 9.42 months post-IO therapy. Significant differences in IO outcome were not seen from multiple non-graph-based unsupervised methods. Furthermore, we demonstrate that patient stratification from DeePaN has the potential to augment the emerging IO biomarker of tumor mutation burden (TMB). Characterization of the subgroups discovered by DeePaN indicates potential to inform IO therapeutic insight, including the enrichment of mutated KRAS and high blood monocyte count in the IO beneficial and IO non-beneficial subgroups, respectively. Our work has proven the concept that graph-based AI is feasible and can effectively integrate high-dimensional genomic and EHR data to meaningfully stratify cancer patients on distinct clinical outcomes, with potential to inform precision oncology.

Predicting Deep Learning Based Multi-Omics Parallel Integration Survival Subtypes in Lung Cancer Using Reverse Phase Protein Array Data

Mortality attributed to lung cancer accounts for a large fraction of cancer deaths worldwide. With increasing mortality figures, the accurate prediction of prognosis has become essential. In recent years, multi-omics analysis has emerged as a useful survival prediction tool. However, the methodology relevant to multi-omics analysis has not yet been fully established and further improvements are required for clinical applications. In this study, we developed a novel method to accurately predict the survival of patients with lung cancer using multi-omics data. With unsupervised learning techniques, survival-associated subtypes in non-small cell lung cancer were first detected using the multi-omics datasets from six categories in The Cancer Genome Atlas (TCGA). The new subtypes, referred to as integration survival subtypes, clearly divided patients into longer and shorter-surviving groups (log-rank test: p = 0.003) and we confirmed that this is independent of histopathological classification (Chi-square test of independence: p = 0.94). Next, an attempt was made to detect the integration survival subtypes using only one categorical dataset. Our machine learning model that was only trained on the reverse phase protein array (RPPA) could accurately predict the integration survival subtypes (AUC = 0.99). The predicted subtypes could also distinguish between high and low risk patients (log-rank test: p = 0.012). Overall, this study explores novel potentials of multi-omics analysis to accurately predict the prognosis of patients with lung cancer.

LncRNA DSCAM-AS1 interacts with YBX1 to promote cancer progression by forming a positive feedback loop that activates FOXA1 transcription network

Rationale: The forkhead box A1 (FOXA1) is a crucial transcription factor in initiation and development of breast, lung and prostate cancer. Previous studies about the FOXA1 transcriptional network were mainly focused on protein-coding genes. Its regulatory network of long non-coding RNAs (lncRNAs) and their role in FOXA1 oncogenic activity remains unknown.

Methods: The Cancer Genome Atlas (TCGA) data, RNA-seq and ChIP-seq data were used to analyze FOXA1 regulated lncRNAs. RT-qPCR was used to detect the expression of DSCAM-AS1, RT-qPCR and Western blotting were used to determine the expression of FOXA1, estrogen receptor α (ERα) and Y box binding protein 1 (YBX1). RNA pull-down and RIP-qPCR were employed to investigate the interaction between DSCAM-AS1 and YBX1. The effect of DSCAM-AS1 on malignant phenotypes was examined through in vitro and in vivo assays.

Results: In this study, we conducted a global analysis of FOXA1 regulated lncRNAs. For detailed analysis, we chose lncRNA DSCAM-AS1, which is specifically expressed in lung adenocarcinoma, breast and prostate cancer. The expression level of DSCAM-AS1 is regulated by two super-enhancers (SEs) driven by FOXA1. High expression levels of DSCAM-AS1 was associated with poor prognosis. Knockout experiments showed DSCAM-AS1 was essential for the growth of xenograft tumors. Moreover, we demonstrated DSCAM-AS1 can regulate the expression of the master transcriptional factor FOXA1. In breast cancer, DSCAM-AS1 was also found to regulate ERα. Mechanistically, DSCAM-AS1 interacts with YBX1 and influences the recruitment of YBX1 in the promoter regions of FOXA1 and ERα.

Conclusion: Our study demonstrated that lncRNA DSCAM-AS1 was transcriptionally activated by super-enhancers driven by FOXA1 and exhibited lineage-specific expression pattern. DSCAM-AS1 can promote cancer progression by interacting with YBX1 and regulating expression of FOXA1 and ERα.

NKX6.1 Represses Tumorigenesis, Metastasis, and Chemoresistance in Colorectal Cancer

Accumulating evidence suggests that NKX6.1 (NK homeobox 1) plays a role in various types of cancer. In our previous studies, we identified NKX6.1 hypermethylation as a promising marker and demonstrated that the NKX6.1 gene functions as a metastasis suppressor through the epigenetic regulation of the epithelial-to-mesenchymal transition (EMT) in cervical cancer. More recently, we have demonstrated that NKX6.1 methylation is related to the chemotherapy response in colorectal cancer (CRC). Nevertheless, the biological function of NKX6.1 in the tumorigenesis of CRC remains unclear. In this study, we showed that NKX6.1 suppresses tumorigenic and metastatic ability both in vitro and in vivo. NKX6.1 represses cell invasion partly through the modulation of EMT. The overexpression of NKX6.1 enhances chemosensitivity in CRC cells. To further explore how NKX6.1 exerts its tumor-suppressive function, we used RNA sequencing technology for comprehensive analysis. The results showed that differentially expressed genes (DEGs) were mainly related to cell migration, response to drug, transcription factor activity, and growth factor activity, suggesting that these DEGs are involved in the function of NKX6.1 suppressing cancer invasion and metastasis. Our results demonstrated that NKX6.1 functions as a tumor suppressor partly by repressing EMT and enhancing chemosensitivity in CRC, making it a potential therapeutic target.

Thyroid transcription factor 1 and p63 expression is associated with survival outcome in patients with non-small cell lung cancer treated with erlotinib

While erlotinib is primarily administered to patients with non-small cell lung cancer with sensitizing epidermal growth factor receptor (EGFR) mutations, it is also prescribed to patients with wild type (wt) EGFR in higher lines of treatment. However, there is no predictive marker for erlotinib efficacy in patients with EGFR wt. Certain immunohistochemical (IHC) parameters, including thyroid transcription factor 1 (TTF1) and p63, have been reported to indicate predictive power in patients with EGFR wt. The present study focused on retrospective data from the University Hospital in Pilsen using the TULUNG register. TTF1 and p63 expression data were extracted from the hospital information system and merged with registry data to calculate progression-free survival (PFS) and overall survival (OS) rates. A cohort of 345 patients with adenocarcinoma (ADC) or squamous cell carcinoma (SCC) exhibited similar erlotinib efficacies when TTF1 and p63 were ignored. However, significant differences were reported in PFS and OS rates of a subgroup of 126 patients where TTF1 and p63 parameters were known. In a univariate analysis, group A (ADC TTF1+/p63-) achieved PFS of 2.6 months, group B (SSC TTF1-/p63+) 1.9 months and group C (did not fit into groups A or B, i.e., ADC TTF1-/p63+ or SCC TTF1+/p63-) 1.4 months (P=0.006). Median OS was 14.2, 19.1 and 5.3 months for A, B and C, respectively (P=0.002). Furthermore, a multivariate analysis demonstrated IHC markers to be the only significant parameters for PFS and OS. Group C had a negative prognostic factor for PFS [hazard ratio (HR), 1.812; P=0.02] and OS (HR=2.367; P=0.01). In conclusion, patients with EGFR wt and lung carcinomas without TTF1 and p63 expression typical for ADC (TTF1+/p633-) or SCC (TTF1-/p63+) do not appear to be suitable candidates for erlotinib treatment.

TTF-1/Nkx2.1 functional connection with mutated EGFR relies on LRIG1 and β-catenin pathways in lung cancer cells

In non-small lung cancer, the expression of the transcription factor TTF-1/Nkx2.1 correlates with the presence of EGFR mutations, therefore TTF-1/Nkx2.1 expression is used to optimize an EGFR testing strategy and to guide clinical treatment. We investigate the molecular mechanisms underlying the functional connection between EGFR and TTF-1/Nkx2.1 gene expression in lung adenocarcinoma. Using the H1975 cell line as a non-small cell lung cancer model system and short hairpin RNA, we have selected clones with TTF-1/Nkx2.1 silenced expression. We have found that Leucine-rich immunoglobulin repeats-1 (LRIG1) gene is a direct target of TTF-1/Nkx2.1 and the transcription factor binding to the LRIG1 genomic sequence inhibits its gene expression. In TTF-1/Nkx2.1 depleted clones, we have found high levels of LRIG1 and decreased presence of EGFR protein. Furthermore, in TTF-1/Nkx2.1 depleted clones we detected a reduced β-catenin level and we provide experimental evidence indicating that TTF-1/Nkx2.1 gene expression is regulated by β-catenin. Published studies indicate that LRIG1 triggers EGFR degradation and that mutated EGFR induces β-catenin activity. Hence, with the present study we show that mutated EGFR, enhancing β-catenin, stimulates TTF-1/Nkx2.1 gene expression and, at the same time, TTF-1/Nkx2.1, down-regulating LRIG1, sustains EGFR pathway. Therefore, LRIG1 and β-catenin mediate the functional connection between TTF-1/Nkx2.1 and mutated EGFR.

Clinicopathological indicators of survival among patients with pulmonary carcinoid tumor

BACKGROUND

Pulmonary carcinoids (PC) are rare malignant neoplasms that cover approximately 1% of all lung cancers. PCs are classified by histological criteria as either typical (TC) or atypical (AC). Histological subtype is the most studied prognostic factor. The aim of this study was to evaluate if other tissue or clinical features are associated with patient outcomes.

MATERIAL AND METHODS

We retrospectively reviewed clinical records of 133 PC patients who underwent operation in the Helsinki University Hospital between 1990 and 2013. Tissue specimens were re-evaluated, processed into tissue microarray format and stained immunohistochemically with serotonin, calcitonin, adrenocorticotropic hormone (ACTH), thyroid transcription factor-1 (TTF-1) and Ki-67. Survival and risk analyses were performed.

RESULTS

Based on histology, 75% (n = 100) of the tumors were TCs and 25% (n = 33) ACs. TCs had higher 10-year disease-specific survival (DSS) rate than ACs (99% (95% CI, 93-100%) for TCs vs. 82% (95% CI, 61-92%) for ACs). Hormonally active tumors expressing serotonin, calcitonin or ACTH were noted in 53% of the specimens but hormonal expression was not associated with DSS. TTF-1 was positive in 78% of the specimens but was not associated with DSS. Ki-67 index varied between <1% and 15%. Ki-67 ≥ 2.5% was associated with shorter DSS (p = .004). The presence of metastatic disease (p = .001), tumor size ≥30 mm (p = .021) and atypical histology (p = .011) were also associated with disease-specific mortality.

CONCLUSIONS

We conclude that PCs are uncommon tumors. When resected, the long-term survival is in general favorable. In this consecutive, single-institution cohort of patients, presence of metastatic disease, tumor size, histological subtype and Ki-67 index were associated with shorter disease-specific survival. As TC and AC have different clinical behaviors, the correct tumor classification at the time of diagnosis is a necessity.

Nkx2-3-A Slippery Slope From Development Through Inflammation Toward Hematopoietic Malignancies

The development of peripheral lymphoid tissues from the mesoderm is the result of a complex convergence combining lymphohematopoietic differentiation with the local specification of nonhematopoietic mesenchymal components. Although the various transcriptional regulators with fate-determining effects in diversifying the mobile leukocyte subsets have been thoroughly studied and identified, the tissue-specific determinants promoting the regional differentiation of resident mesenchyme are less understood. Of these factors, various members of the NK-class Nkx paralogues have emerged as key regulators for the organogenesis of spleen and mucosal lymphoid tissues, and recent data have also indicated their involvement in various pathological events, including gut inflammation and hematopoietic malignancies. Here, we summarize available data on the roles of Nkx2-3 in lymphoid tissue development and discuss its possible value as a developmental marker and disease-associated pathogenic trait.

Network modeling of kinase inhibitor polypharmacology reveals pathways targeted in chemical screens

Small molecule screens are widely used to prioritize pharmaceutical development. However, determining the pathways targeted by these molecules is challenging, since the compounds are often promiscuous. We present a network strategy that takes into account the polypharmacology of small molecules in order to generate hypotheses for their broader mode of action. We report a screen for kinase inhibitors that increase the efficacy of gemcitabine, the first-line chemotherapy for pancreatic cancer. Eight kinase inhibitors emerge that are known to affect 201 kinases, of which only three kinases have been previously identified as modifiers of gemcitabine toxicity. In this work, we use the SAMNet algorithm to identify pathways linking these kinases and genetic modifiers of gemcitabine toxicity with transcriptional and epigenetic changes induced by gemcitabine that we measure using DNaseI-seq and RNA-seq. SAMNet uses a constrained optimization algorithm to connect genes from these complementary datasets through a small set of protein-protein and protein-DNA interactions. The resulting network recapitulates known pathways including DNA repair, cell proliferation and the epithelial-to-mesenchymal transition. We use the network to predict genes with important roles in the gemcitabine response, including six that have already been shown to modify gemcitabine efficacy in pancreatic cancer and ten novel candidates. Our work reveals the important role of polypharmacology in the activity of these chemosensitizing agents.

The NANCI-Nkx2.1 gene duplex buffers Nkx2.1 expression to maintain lung development and homeostasis

A subset of long noncoding RNAs (lncRNAs) is spatially correlated with transcription factors (TFs) across the genome, but how these lncRNA–TF gene duplexes regulate tissue development and homeostasis is unclear. Here, Herriges et al. identified a feedback loop within the NANCI–Nkx2.1 gene duplex that is essential for buffering Nkx2.1 expression, lung epithelial cell identity, and tissue homeostasis.

A subset of long noncoding RNAs (lncRNAs) is spatially correlated with transcription factors (TFs) across the genome, but how these lncRNA–TF gene duplexes regulate tissue development and homeostasis is unclear. We identified a feedback loop within the NANCI (Nkx2.1-associated noncoding intergenic RNA)–Nkx2.1 gene duplex that is essential for buffering Nkx2.1 expression, lung epithelial cell identity, and tissue homeostasis. Within this locus, Nkx2.1 directly inhibits NANCI, while NANCI acts in cis to promote Nkx2.1 transcription. Although loss of NANCI alone does not adversely affect lung development, concurrent heterozygous mutations in both NANCI and Nkx2.1 leads to persistent Nkx2.1 deficiency and reprogramming of lung epithelial cells to a posterior endoderm fate. This disruption in the NANCI–Nkx2.1 gene duplex results in a defective perinatal innate immune response, tissue damage, and progressive degeneration of the adult lung. These data point to a mechanism in which lncRNAs act as rheostats within lncRNA–TF gene duplex loci that buffer TF expression, thereby maintaining tissue-specific cellular identity during development and postnatal homeostasis.

A homeobox protein, NKX6.1, up-regulates interleukin-6 expression for cell growth in basal-like breast cancer cells

Among breast cancer subtypes, basal-like breast cancer is particularly aggressive, and research on the molecules involved in its pathology might contribute to therapy. In this study, we found that expression of NKX6.1, a homeobox transcription factor, is higher in basal-like breast cancer than in other subtypes. In loss-of-function experiments on basal-like breast cancer cell lines, NKX6.1-depleted cells exhibited reduced cell growth. Because cytokine interleukin-6 (IL-6) is expressed in basal-like breast cancer, and increases cell growth, we analyzed expression levels of IL6, an IL-6 gene, and observed reduced IL6 expression in NKX6.1-depleted cells. In a reporter assay, IL6 promoter activity was reduced by loss of NKX6.1 function. A pull-down assay showed that NKX6.1 binds to the proximal region in IL6 promoter. These results indicate that NKX6.1 directly up-regulates IL6 expression. To investigate further, we established cells with forced expression of IL-6. We observed that exogenous IL-6 expression restored the reduced cell growth of NKX6.1-depleted cells. Furthermore, orthotopic xenografts showed that NKX6.1-depleted cells lost the capacity for tumor formation. We therefore conclude that NKX6.1 is a factor for IL-6-regulated growth and tumor formation in basal-like breast cancer. Our findings facilitate profound understanding of basal-like breast cancer, and the development of suitable therapy.

Informative gene selection and the direct classification of tumors based on relative simplicity

BACKGROUND

Selecting a parsimonious set of informative genes to build highly generalized performance classifier is the most important task for the analysis of tumor microarray expression data. Many existing gene pair evaluation methods cannot highlight diverse patterns of gene pairs only used one strategy of vertical comparison and horizontal comparison, while individual-gene-ranking method ignores redundancy and synergy among genes.

RESULTS

Here we proposed a novel score measure named relative simplicity (RS). We evaluated gene pairs according to integrating vertical comparison with horizontal comparison, finally built RS-based direct classifier (RS-based DC) based on a set of informative genes capable of binary discrimination with a paired votes strategy. Nine multi-class gene expression datasets involving human cancers were used to validate the performance of new method. Compared with the nine reference models, RS-based DC received the highest average independent test accuracy (91.40%), the best generalization performance and the smallest informative average gene number (20.56). Compared with the four reference feature selection methods, RS also received the highest average test accuracy in three classifiers (Naïve Bayes, k-Nearest Neighbor and Support Vector Machine), and only RS can improve the performance of SVM.

CONCLUSIONS

Diverse patterns of gene pairs could be highlighted more fully while integrating vertical comparison with horizontal comparison strategy. DC core classifier can effectively control over-fitting. RS-based feature selection method combined with DC classifier can lead to more robust selection of informative genes and classification accuracy.

TGF-β regulates the proliferation of lung adenocarcinoma cells by inhibiting PIK3R3 expression

PIK3R3, an isoform of class IA phosphoinositide 3-kinase (PI3K), specifically interacts with cell proliferation regulators, such as retinoblastoma and proliferation cell nuclear antigen, to promote cell proliferation. However, the mechanisms behind the upstream signaling pathway of PIK3R3 remain unclear to date. This study showed that PIK3R3 expression was regulated by transforming growth factor-β (TGF-β) signaling and that PIK3R3 mediated the TGF-β-induced inhibition of lung adenocarcinoma cell proliferation. TGF-β down-regulated PIK3R3 expression in lung adenocarcinoma cells. However, this TGF-β-induced inhibition of cell proliferation can be attenuated by PIK3R3 overexpression. In addition, TGF-β can attenuate the transcriptional activity of NKX2.1, a transcription factor that binds to the promoter of PIK3R3. This result indicated that TGF-β regulated PIK3R3 expression by targeting NKX2.1. We confirmed the correlation between NKX2.1 and PIK3R3 in clinical samples. Therefore, the TGF-β/NKX2.1/PIK3R3 axis is crucial in the TGF-β-induced inhibition of cell proliferation, and the NKX2.1/PIK3R3 axis might become a target in TGF-β receptor-repressed lung adenocarcinoma.

Association between gene expression profiles and clinical outcome of pemetrexed-based treatment in patients with advanced non-squamous non-small cell lung cancer: exploratory results from a phase II study

INTRODUCTION

We report exploratory gene-expression profiling data from a single-arm Phase-II-study in patients with non-squamous (ns)NSCLC treated with pemetrexed and cisplatin. Previously disclosed results indicated a significant association of low thymidylate-synthase (TS)-expression with longer progression-free and overall survival (PFS/OS).

METHODS

Treatment-naïve nsNSCLC patients (IIIB/IV) received 4 cycles of pemetrexed/cisplatin; non-progressing patients continued on pemetrexed-maintenance. Diagnostic tissue-samples were used to assess TS-expression by immunohistochemistry (IHC) and mRNA-expression array-profiling (1,030 lung cancer-specific genes). Cox proportional-hazard models were applied to explore the association between each gene and PFS/OS. Genes significantly correlated with PFS/OS were further correlated with TS-protein expression (Spearman-rank). Unsupervised clustering was applied to all evaluable samples (n = 51) for all 1,030 genes and an overlapping 870-gene subset associated with adenocarcinoma (ADC, n = 47).

RESULTS

51/70 tissue-samples (72.9%) were evaluable; 9 of 1,030 genes were significantly associated with PFS/OS (unadjusted p < 0.01, genes: Chromosome 16 open reading frame 89, napsin A, surfactant protein B, aquaporin 4, TRAF2- and Nck-interacting kinase, Lysophosphatidylcholine acyltransferase 1, Interleukin 1 receptor type II, NK2 homeobox 1, ABO glycosyl-transferase); expression for all except IL1R2 correlated negatively with nuclear TS-expression (statistically significant for 5/8 genes, unadjusted p<0.01). Cluster-analysis based on 1,030 genes revealed no clear trend regarding PFS/OS; the ADC-based cluster analysis identified 3 groups (n = 21/11/15) with median (95%CI) PFS of 8.1(6.9,NE)/2.4(1.2,NE)/4.4(1.2,NE) months and OS of 20.3(17.5,NE)/4.3(1.4,NE)/8.3(3.9,NE) months, respectively.

CONCLUSIONS

These exploratory gene-expression profiling results describe genes potentially linked to low TS-expression. Nine genes were significantly associated with PFS/OS but could not be differentiated as prognostic or predictive as this was a single-arm study. Although these hypotheses-generating results are interesting, they provide no evidence to change the current histology-based treatment approach with pemetrexed.