Gene identification for risk of relapse in stage I lung adenocarcinoma patients: a combined methodology of gene expression profiling and computational gene network analysis

RNF152 Suppresses Fatty Acid Oxidation and Metastasis of Lung Adenocarcinoma by Inhibiting IRAK1-Mediated AKR1B10 Expression

Lung adenocarcinoma (LUAD) is a common subtype of primary lung cancer. Fatty acid oxidation plays a key role in LUAD development by providing energy for tumor cells. This study aimed to identify the role of ring finger protein 152 (RNF152) in LUAD. RNF152 was down-regulated in LUAD, and low RNF152 expression correlated with a poor prognosis in LUAD patients. RNF152 overexpression inhibited the proliferation and malignant phenotype of LUAD cells, whereas RNF152 knockdown exerted an opposite effect. Tumor cells overexpressing RNF152 showed less fatty acid oxidation compared with control cells, whereas RNF152 knockdown induced fatty acid uptake and oxidation. Further analysis revealed the binding reaction between RNF152 and interleukin-1 receptor-associated kinase 1 (IRAK1). RNF152 reduced the stability of IRAK1 in LUAD cells by promoting its ubiquitination. RNF152-overexpressed tumor cells exhibited a significantly lower level of Aldo-Keto reductase family 1 member 10 (AKR1B10), whereas up-regulation of IRAK1 restored the expression of AKR1B10 in RNF152-overexpressed cells. Furthermore, up-regulation of IRAK1 eliminated the antitumor effect of RNF152 in LUAD cells. Mouse xenograft models confirmed the inhibitory effect of RNF152 on the tumorigenesis and metastasis of LUAD. Taken together, RNF152 played a tumor suppressive role in LUAD by promoting IRAK1 ubiquitination and IRAK1-mediated down-regulation of AKR1B10, thereby reversing the malignant phenotype of LUAD.

Searching for New Molecular Targets for Oral Squamous Cell Carcinoma with a View to Clinical Implementation of Precision Medicine

Head and neck cancer, including oral squamous cell carcinoma (OSCC), is the eighth most common malignancy globally and is characterized by local invasiveness and high nodal metastatic potential. The OSCC incidence is also increasing, and the number of deaths is also rising steadily in Japan. The development of molecular markers to eradicate OSCC is an urgent issue for humankind. The increase in OSCC despite the declining smoking rate may be due to several viral infections through various sexual activities and the involvement of previously unfocused carcinogens, and genetic alterations in individual patients are considered to be more complicated. Given this situation, it is difficult to combat OSCC with conventional radiotherapy and chemotherapy using cell-killing anticancer drugs alone, and the development of precision medicine, which aims to provide tailor-made medicine based on the genetic background of each patient, is gaining attention. In this review article, the current status of the comprehensive search for driver genes and biomarkers in OSCC will be briefly described, and some of the candidates for novel markers of OSCC that were found will be outlined.

Neuropsychiatric Adverse Events of Montelukast: An Analysis of Real-World Datasets and drug-gene Interaction Network

Montelukast is a selective leukotriene receptor antagonist that is widely used to treat bronchial asthma and nasal allergy. To clarify the association between montelukast and neuropsychiatric adverse events (AEs), we evaluated case reports recorded between January 2004 and December 2018 in the Food and Drug Administration Adverse Event Reporting System (FAERS). Furthermore, we elucidated the potential toxicological mechanisms of montelukast-associated neuropsychiatric AEs through functional enrichment analysis of human genes interacting with montelukast. The reporting odds ratios of suicidal ideation and depression in the system organ class of psychiatric disorders were 21.5 (95% confidence interval (CI): 20.3–22.9) and 8.2 (95% CI: 7.8–8.7), respectively. We explored 1,144 human genes that directly or indirectly interact with montelukast. The molecular complex detection (MCODE) plug-in of Cytoscape detected 14 clusters. Functional analysis indicated that several genes were significantly enriched in the biological processes of “neuroactive ligand–receptor interaction.” “Mood disorders” and “major depressive disorder” were significant disease terms related to montelukast. Our retrospective analysis based on the FAERS demonstrated a significant association between montelukast and neuropsychiatric AEs. Functional enrichment analysis of montelukast-associated genes related to neuropsychiatric symptoms warrant further research on the underlying pharmacological mechanisms.

TFF-1 Functions to Suppress Multiple Phenotypes Associated with Lung Cancer Progression

Introduction

Trefoil Factor (TFF) is a member of a protein family comprised of three isoforms, of which TFF-1 exhibits antithetical functions; promotion or suppression of cell proliferation, survival and invasion, depending on the cancer type. However, the pathobiological function of TFF-1 in lung carcinoma has been still unclear.

Methods

We examined the expression and secretion of TFF-1 using cultured human lung carcinoma cells by immunoblotting, immunofluorescence, enzyme-linked immunosorbent assay and quantitative real-time PCR analyses. The effects of TFF-1 on various phenotypes were analyzed in two cell lines, including those transfected with cDNA encoding TFF-1. Cell proliferation and death were examined by hemocytometer cell counting and by colorimetric viability/cytotoxicity assay. Cell cycle profile, migration and invasion were also examined by flow cytometry, wound healing assay and Matrigel Transwell assay, respectively. The effect of TFF-1 overexpression was confirmed by additional transfection of TFF-1-specific siRNA.

Results

Endogenous TFF-1 protein expression and secretion into the media were observed exclusively in adenocarcinoma-derived cell lines. Forced overexpression of TFF-1 drove cell cycle transition, while the proliferation decreased by 19% to 25% due to increased cell death. This cell death was predominantly caused by apoptosis, as assessed by the activation of caspase 3/7. Cell migration was also suppressed by 71% to 82% in TFF-1-transfected cells. The suppressive effect of TFF-1 on proliferation and migration was restored by transfection of TFF-1 siRNA. Moreover, invasion was also suppressed to 77% to 83% in TFF-1-transfected cells.

Conclusion

These findings reveal that TFF-1 functions as a suppressor of cancer proliferation by induction of apoptosis, cell migration and invasion and thus may provide a synergistic target for potential treatment strategies for human lung carcinoma.

Gentamicin-induced hearing loss: A retrospective study using the Food and Drug Administration Adverse Event Reporting System and a toxicological study using drug-gene network analysis

The objectives of the study were to evaluate the relationship between gentamicin (GEN) and hearing loss using the Food and Drug Administration Adverse Event Reporting system (FAERS) database and elucidate the potential toxicological mechanism of GEN-induced hearing loss through a drug–gene network analysis.

Using the preferred terms and standardized queries from the Medical Dictionary for Regulatory Activities, we calculated the reporting odds ratios (RORs). We extracted GEN-associated genes (seed genes) and analyzed drug−gene interactions using the ClueGO plug-in in the Cytoscape software and the DIseAse MOdule Detection (DIAMOnD) algorithm.

The lower limit of the 95% confidence interval (CI) of the ROR for aminoglycosides (AG) antibacterials was over 1, and the ROR was 5.5 (5.1–6.0). We retrieved 17 seed genes related to GEN from the PharmGKB and Drug Gene Interaction databases. In total, 1018 human genes interacting with GEN were investigated using ClueGO. Through Molecular Complex Detection (MCODE) analysis, we identified 17 local gene clusters. The nodes and edges of the highest-ranked local gene cluster named “Cluster 1” were 30 and 433, respectively. According to the ClueGO analysis using the Kyoto Encyclopedia of Genes and Genomes (KEGG), Cluster 1 genes were highly enriched in “oxidative phosphorylation.” According to the ClueGO analysis using ClinVar, Cluster 1 genes were highly enriched in “mitochondrial diseases,” “mitochondrial complex I deficiency,” “hereditary hearing loss and deafness,” and “Leigh syndrome.” We identified 60 GEN-associated genes using the DIAMOnD algorithm. Several GEN-associated genes in the DIAMOnD algorithm were highly enriched in “PI3K-Akt signaling pathway,” “Ras signaling pathway,” “focal adhesion,” “MAPK signaling pathway,” “regulation of actin cytoskeleton,” “oxidative phosphorylation,” and “ECM-receptor interaction.”

Our analysis demonstrated an association between several AGs and hearing loss using the FAERS database. Drug−gene network analysis demonstrated that GEN may be associated with oxidative phosphorylation-associated genes and integrin genes, which may be associated with hearing loss.

The Role of AKR1B10 in Physiology and Pathophysiology

AKR1B10 is a human nicotinamide adenine dinucleotide phosphate (NADPH)-dependent reductase belonging to the aldo-keto reductase (AKR) 1B subfamily. It catalyzes the reduction of aldehydes, some ketones and quinones, and interacts with acetyl-CoA carboxylase and heat shock protein 90α. The enzyme is highly expressed in epithelial cells of the stomach and intestine, but down-regulated in gastrointestinal cancers and inflammatory bowel diseases. In contrast, AKR1B10 expression is low in other tissues, where the enzyme is upregulated in cancers, as well as in non-alcoholic fatty liver disease and several skin diseases. In addition, the enzyme's expression is elevated in cancer cells resistant to clinical anti-cancer drugs. Thus, growing evidence supports AKR1B10 as a potential target for diagnosing and treating these diseases. Herein, we reviewed the literature on the roles of AKR1B10 in a healthy gastrointestinal tract, the development and progression of cancers and acquired chemoresistance, in addition to its gene regulation, functions, and inhibitors.

Identification of molecular signatures associated with early relapse after complete resection of lung adenocarcinomas

The only potentially curative treatment for lung adenocarcinoma patients remains complete resection of early-stage tumors. However, many patients develop recurrence and die of their disease despite curative surgery. Underlying mechanisms leading to establishment of systemic disease after complete resection are mostly unknown. We therefore aimed at identifying molecular signatures of resected lung adenocarcinomas associated with the risk of an early relapse. The study comprised 89 patients with totally resected stage IA–IIIA lung adenocarcinomas. Patients suffering from an early relapse within two years after surgery were compared to patients without a relapse in two years. Patients were clinically and molecular pathologically characterized. Tumor tissues were immunohistochemically analyzed for the expression of Ki67, CD45, CD4, CD8, PD1, PD-L1, PD-L2 and CD34, by Nanostring nCounter PanCancer Immune Profiling Panel as well as a comprehensive methylome profiling using the Infinium MethylationEPIC BeadChip. We detected differential DNA methylation patterns as well as significantly differentially expressed genes associated with an early relapse after complete resection. Especially, CD1A was identified as a potential biomarker, whose reduced expression is associated with an early relapse. These findings might help to develop biomarkers improving risk assessment and patient selection for adjuvant therapy as well as establish novel targeted therapeutic strategies.

INSL4 as prognostic marker for proliferation and invasiveness in Non-Small-Cell Lung Cancer

Non-small-cell-lung cancer accounts for 80-85% of all forms of lung cancer as leading cause of cancer-related death in human. Despite remarkable advances in the diagnosis and therapy of lung cancer, no significant improvements have thus far been achieved in terms of patients' prognosis. Here, we investigated the role of INSL4 - a member of the relaxin-family - in NSCLC. We overexpressed INSL4 in NSCLC cells to analyse in vitro the growth rate and the tumourigenic features. We investigated the signalling pathways engaged in INSL4 overexpressing cells and the tumour growth ability by studying the tumour development in a patient derived tumour xenograft mouse model. We found an INSL4 cell growth promoting effect in vitro in H1299 cells and in vivo in NOD/SCID mice. Surprisingly, in NSCLC-A549 cells, INSL4 overexpression has not similar effect, despite huge basal INSL4-mRNA expression respect to H1299. The INSL4-mRNA analysis of eight different NSCLC-derived cell lines, revealed highly difference in the INSL4-mRNA amount. Transfection of NSCLC lines with INSL4-Myc showed huge level of INSL4-mRNA with a very low amount of protein expressed. Notably, similar discrepancy has been observed in NSCLC patients. However, in a cohort of NSCLC patients analysing a database, we found a significant inverse correlation between INSL4 expression and Overall Survival. By combining the in vitro and in vivo results, suggest that in patients whose NSCLC adenocarcinoma spontaneously expressed high levels of INSL4 post-transcriptional modifications affecting INSL4 do not allow to assess precision therapy in selected patients without consider protein INSL4 amount.

Identification of oral squamous cell carcinoma markers MUC2 and SPRR1B downstream of TANGO

PURPOSE

Transport and Golgi organization protein 1 (TANGO) promotes angiogenesis and lymphangiogenesis in oral squamous cell carcinoma (OSCC). To elucidate the underlying mechanisms, this study aims to identify and characterize elements downstream of TANGO that mediate its involvement in OSCC.

METHODS

In this study, microarray analysis compared gene expression between control and TANGO-repressed HSC3 cells. Protein expression in 213 OSCC tissue samples was analyzed immunohistochemically.

RESULTS

TANGO repression decreased or increased expression of Mucin 20 (MUC20) and small proline-rich protein 1B (SPRR1B), respectively. MUC20 increased the growth and invasiveness of OSCC cells via altered matrix metalloproteinase (MMP)-2 and E-cadherin expression and c-met phosphorylation. MUC20 induced angiogenesis and lymphangiogenesis by activating vascular endothelial growth factors A and C. In well-differentiated OSCC, SPRR1B expression was high (P = 0.0091) and correlated with keratinization markers and promoted proliferation by inducing mitogen-activated protein kinase p38 phosphorylation. MUC20 expression correlated significantly with clinical stage (P = 0.0024), lymph node metastasis (P = 0.0036), and number of blood and lymph vessels (P < 0.0001). MUC20-expressing cases had a significantly worse prognosis than non-expressing cases (P < 0.0001).

CONCLUSION

MUC20 and SPRR1B located downstream of TANGO may be useful molecular markers for OSCC.

High PD-L1/IDO-2 and PD-L2/IDO-1 Co-Expression Levels Are Associated with Worse Overall Survival in Resected Non-Small Cell Lung Cancer Patients

Programmed death ligand 1 (PD-L1) expression is a predictive biomarker of the success of PD-1/PD-L1 inhibitor therapy for patients with advanced non-small cell lung cancer (NSCLC) but its role as a prognostic marker for early-stage resectable NSCLC remains unclear. We studied gene expression levels of immune-related genes PD-1, PD-L1, PD-L2, IDO-1, IDO-2 and INFγ in tumor tissue of surgically resected NSCLC and correlated the finding with clinicopathological features and patient outcomes. A total of 191 consecutive early-stage NSCLC patients who underwent curative pulmonary resection were studied. The mRNA expression levels of immune-related genes were evaluated by quantitative reverse transcription polymerase chain reaction (qRT-PCR) using RT² Profiler PCR Arrays (Qiagen). PD-1, PD-L2 and IDO-2 gene expression levels were significantly higher in patients with squamous histology (p = 0.001, p = 0.021 and p < 0.001; respectively). PD-1, PD-L1 and IDO-2 gene expression levels were significantly higher in patients with higher stage (p = 0.005, p = 0.048 and p = 0.002, respectively). The univariate analysis for recurrence-free survival (RFS) and overall survival (OS) showed that patients with higher levels of three-genes (PD-L1/PD-L2/INFγ) (hazard ratio (HR)) 1.90 (95% confidence interval (CI), 1.13–3.21), p = 0.015) were associated with a worse RFS, while patients with higher levels of both genes (PD-L1/IDO-2) or (PD-L2/IDO-1) were associated with a worse OS (HR 1.63 95% CI, 1.06–2.51, p = 0.024; HR 1.54 95% CI, 1.02–2.33, p = 0.04; respectively). The multivariate interaction model adjusted for histology and stage confirmed that higher levels of three genes (PD-L1/PD-L2/INFγ) were significantly associated with worse RFS (HR 1.98, p = 0.031) and higher levels of both genes (PD-L1/IDO-2) and (PD-L2/IDO-1) with worse OS (HR 1.98, p = 0.042, HR 1.92, p = 0.022). PD-L1/IDO-2 and PD-L2/IDO-1 co-expression high levels are independent negative prognostic factors for survival in early NSCLC. These features may have important implications for future immune-checkpoint therapeutic approaches.

PRSS1 Upregulation Predicts Platinum Resistance in Ovarian Cancer Patients

Ovarian cancer is the most frequent cause of death among gynecologic malignancies. A total of 80% of patients who have completed platinum-based chemotherapy suffer from relapse and develop resistance within 2 years. In the present study, we obtained patients' complete platinum (cisplatin and carboplatin) medication information from The Cancer Genome Atlas database and then divided them into two categories: resistance and sensitivity. Difference analysis was performed to screen differentially expressed genes (DEgenes) related to platinum response. Subsequently, we annotated DEgenes into the protein–protein interaction network as seed nodes and analyzed them by random walk. Finally, second-ranking protease serine 1 gene (PRSS1) was selected as a candidate gene for verification analysis. PRSS1's expression pattern was continuously studied in Oncomine and cBio Cancer Genomic Portal databases, revealing the key roles of PRSS1 in ovarian cancer formation. Hereafter, we conducted in-depth explorations on PRSS1's platinum response to ovarian cancer through tissue and cytological experiments. Quantitative real-time polymerase chain reaction and Western blot assay results indicated that PRSS1 expression levels in platinum-resistant samples (tissue/cell) were significantly higher than in samples sensitive to platinum. By cell transfection assay, we observed that knockdown of PRSS1 reduced the resistance of ovarian cancer cells to cisplatin. Meanwhile, overexpression of PRSS1 increased the resistance to cisplatin. In conclusion, we identified a novel risk gene PRSS1 related to ovarian cancer platinum response and confirmed its key roles using multiple levels of low-throughput experiments, revealing a new treatment strategy based on a novel target factor for overcoming cisplatin resistance in ovarian cancer.

Transcriptomic analysis of micropapillary high grade T1 urothelial bladder cancer

No consensus currently exist on the optimal treatment of patients with high-risk nonmuscle invasive (HGT1) micropapillary variant of bladder cancer (MPBC). Transcripsome analysis may allow stratification of MPBC-HGT1 enabling prediction of recurrence and guide therapeutic management for individual patients. Whole transcriptome RNA-Sequencing of tumors from 23 patients with MPBC-HGT1 and 64 conventional urothelial carcinomas (cUC) (reference set) was performed. Differentially expressed genes between MPBC-HGT1 and cUC-HGT1 were explored. Cox proportional hazard models and Kapplan–Meier methods were used to assess the relation between time to progression (TTP) and individual gene expression adjusting for clinical covariates. Over 3000 genes were differentially expressed in MPBC-HGT1 as compared with cUC-HGT1 and a 26-gene signature is characteristic of MPBC within HGT1. A set of three genes; CD36, FAPB3 and RAETE1; were significantly associated with TTP. High expression of FABP3 and CD36 were associated with shorter TTP (p = 0.045 and p = 0.08) as was low expression of RAET1E (p = 0.01). Our study suggest that a 26-gene signature can define MPBC-HGT1 within conventional urothelial carcinomas. A prognostic risk index of three genes (FABP3, CD36 and RAET1E) was found to be associated with shorter TTP and may help classify a group of patients with MPBC-HGT1 with high-risk of early progression. These observations might have implications in terms of radical cystectomy recommendation in MPBC patients.

Upregulation of lncRNA ZFAS1 promotes lung adenocarcinoma progression by sponging miR-1271-5p and upregulating FRS2

Background

Nowadays, the important roles of long non‐coding RNAs (LncRNAs) in lung adenocarcinoma (LAD) is being increasingly recognized. The purpose of this study was to explore the regulatory mechanism of lncRNA ZFAS1 in LAD.

Methods

The expression and function of lncRNA ZFAS1 were assessed by RT‐qPCR, CCK‐8, transwell and dual luciferase reporter assays.

Results

Upregulation of lncRNA ZFAS1 was found in LAD tissues and cells. Knockdown of lncRNA ZFAS1 restrained cell proliferation, migration and invasion in LAD cells. In addition, we determined that lncRNA ZFAS1 could directly bind to miR‐1271‐5p. MiR‐1271‐5p functioned as a tumor suppressor in LAD, and lncRNA ZFAS1 promoted LAD development by downregulating miR‐1271‐5p. Furthermore, FRS2 was a direct target of miR‐1271‐5p. FRS2 promoted progression of LAD by mediating lncRNA ZFAS1/miR‐1271‐5p axis.

Conclusions

LncRNA ZFAS1 promotes cell proliferation, migration and invasion in LAD by downregulating miR‐1271‐5p or upregulating FRS2.

CLCA2 expression is associated with survival among African American women with triple negative breast cancer

Purpose

Black/African American (AA) women are twice as likely to be diagnosed with triple negative breast cancer (TNBC) compared to whites, an aggressive breast cancer subtype associated with poor prognosis. There are no routinely used targeted clinical therapies for TNBC; thus there is a clear need to identify prognostic markers and potential therapeutic targets.

Methods

We evaluated expression of 27,016 genes in 155 treatment-naïve TN tumors from AA women in Detroit. Associations with survival were evaluated using Cox proportional hazards models adjusting for stage and age at diagnosis, and p-values were corrected using a false discovery rate. Our validation sample consisted of 494 TN tumors using four publically available data sets. Meta-analyses were performed using summary statistics from the four validation results.

Results

In the Detroit AA cohort, CLCA2 [Hazard ratio (HR) = 1.56, 95% confidence interval (CI) 1.31–1.86, nominal p = 5.1x10-7, FDR p = 0.014], SPIC [HR = 1.47, 95%CI 1.26–1.73, nominal p = 1.8x10-6, FDR p = 0.022], and MIR4311 [HR = 1.57, 95% CI 1.31–1.92, nominal p = 2.5x10-5, FDR p = 0.022] expression were associated with overall survival. Further adjustment for treatment and breast cancer specific survival analysis did not substantially alter effect estimates. CLCA2 was also associated with increased risk of death in the validation cohorts [HR = 1.14, 95% CI 1.05–1.24, p = 0.038, p-heterogeneity = 0.88].

Conclusions

We identified CLCA2 as a potential prognostic marker for TNBC in AA women.

Silencing of LINC00461 enhances radiosensitivity of lung adenocarcinoma cells by down-regulating HOXA10 via microRNA-195

Lung adenocarcinoma is recognized as one of the most recurrent tumours in adults. Long non‐coding RNAs (lncRNAs) are non–protein‐coding transcripts and have been demonstrated to regulate biological functions during tumorigenesis. Our study aims to investigate the underlying molecular mechanisms of LINC00461/microRNA‐195 (miR‐195)/HOXA10 responsible for its involvement in lung adenocarcinoma. We firstly selected differentially expressed lncRNAs and genes by the Cancer Genome Atlas (TCGA) database and Gene Expression Omnibus (GEO). The functional role of LINC00461 in lung adenocarcinoma was then determined using ectopic expression, knockdown and reporter assay experiments. Besides, we detected the expression profiles of LINC00461, miR‐195, HOXA10 and apoptosis‐ and invasion‐related genes. Cell proliferation, migration and invasion were evaluated. In vivo tumour formation ability was analysed. Overexpressed LINC00461 and HOXA10 but down‐regulated miR‐195 were observed in primary and metastatic lung adenocarcinoma. LINC00461 negatively regulated miR‐195, while miR‐195 negatively regulated HOXA10. Forced LINC00461 expression decreased expression of miR‐195 and Bax, increased expression of HOXA10, MMP‐2, MMP‐9 and Bcl‐2, promoted cell proliferation, migration and invasion as well as tumour formation, and enhanced radiosensitivity of lung adenocarcinoma cells. However, these effects were reversed by lentivirus‐mediated miR‐195–forced expression, thereby suggesting that miR‐195 could antagonize the harmful effect of LINC00461 on lung adenocarcinoma cells. Collectively, the present study provides evidence supporting the inhibitory effect of LINC00461 silencing on lung adenocarcinoma, which suppresses lung adenocarcinoma cell migration, invasion and radiosensitivity via HOXA10 by binding to miR‐195, which provides a promising basis for the targeted intervention treatment for human lung adenocarcinoma.

Tumor-Promoting Activity of Long Noncoding RNA LINC00466 in Lung Adenocarcinoma via miR-144-Regulated HOXA10 Axis

Previous investigations have implicated long noncoding RNAs in lung adenocarcinoma, which is an aggressive disease with poor prognosis and high mortality. Through the alteration of lung adenocarcinoma-related long noncoding RNA and miRNA based on microarray analysis, our aim was to understand the role of LINC00466 and miR-144 in lung adenocarcinoma progression. The relationship among LINC00466, miR-144, and HOXA10 was also verified. Moreover, to examine whether the LINC00466/miR-144/HOXA10 axis contributed to the cellular processes in lung adenocarcinoma, A549 and XWLC-05 cells were transduced with siRNA LINC00466, siRNA HOXA10, or miR-144 mimic plasmids. Highly expressed LINC00466 and HOXA10 and lowly expressed miR-144 were eventually revealed in lung adenocarcinoma tissues. HOXA10 was down-regulated in response to the overexpression of miR-144, whereas inhibition of LINC00466 decreased its binding to miR-144, thereby up-regulating miR-144, which, in turn, halted the lung adenocarcinoma progression. LINC00466 silencing or miR-144 up-regulation exerted an inhibitory role in the tumorigenicity, invasion, migration, and proliferation, and it also promoted apoptosis of lung adenocarcinoma cells. Furthermore, tumor formation was inhibited by knockdown of LINC00466 or overexpression of miR-144. Taken together, LINC00466 could restrain the miR-144 expression to up-regulate HOXA10 and, therefore, promote lung adenocarcinoma.

Role of fatty acid binding proteins (FABPs) in cancer development and progression

Fatty acid binding proteins (FABPs) are small, water soluble proteins that bind long chain fatty acids and other biologically active ligands to facilitate intracellular localization. Twelve FABP family members have been identified to date, with 10 isoforms expressed in humans. Functionally, FABPs are important in fatty acid metabolism and transport, with distinct family members having the capacity to influence gene transcription. Expression of FABPs is usually cell/tissue specific to one predominant FABP family member. Dysregulation of FABP expression can occur through genetic mutation and/or environmental-lifestyle influences. In addition to intracellular function, exogenous, circulating FABP expression can occur and is associated with specific disease states such as insulin resistance. A role for FABPs is increasingly being reported in tumor biology with elevated exogenous FABP expression being associated with tumor progression and invasiveness. However, a less clear role has been appreciated for dysregulated FABP expression during cell transformation and early expansion.

Identification of potential diagnostic and therapeutic target genes for lung squamous cell carcinoma

The purpose of this study was to identify potential molecular markers of lung squamous cell carcinoma (LUSC). Three datasets containing LUSC mRNA sequencing data were downloaded from the Gene Expression Omnibus, The Cancer Genome Atlas and the Gene Expression Profiling Interactive Analysis databases. These datasets were used to identify significantly differentially expressed genes (DEGs) in LUSC. A protein-protein interaction network of the DEGs was constructed followed by Gene Ontology, Kyoto Encyclopedia of Genes and Genomes and overall survival analyses of the DEGs. A total of 37 DEGs between LUSC and normal tissues were identified, including 26 downregulated genes and 11 upregulated genes. Biological Process enrichment analysis revealed that the DEGs were mainly enriched in ‘cell adhesion’, ‘cell-matrix adhesion’, ‘anatomical structure morphogenesis’, ‘ECM-receptor interaction’ and ‘focal adhesion’. Overall survival analysis demonstrated that transcription factor 21, α-2-macroglobulin, acyl-CoA synthetase long chain family member 5, integrin subunit β8, meiotic nuclear divisions 1 and secretoglobin family 1A member 1 were significantly associated with the occurrence and development of lung cancer, and these genes were selected as hub genes. The results obtained in the present study may aid the elucidation of the molecular mechanisms involved in the development of LUSC and may provide potential targets for LUSC treatment.

Comprehensive clinical implications of homeobox A10 in 3,199 cases of non-small cell lung cancer tissue samples combining qRT-PCR, RNA sequencing and microarray data

In the current study, we proposed to explore the potential role and related signaling pathways of Homobox A10 (HOXA10) in non-small cell lung cancer (NSCLC). HOXA10 levels in lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC) were detected by qRT-PCR and the expression of HOXA10 was significantly up-regulated in the NSCLC tissue of all 55 pairs (P = 0.037). Overexpression of HOXA10 was closely correlated with the clinical stage of LUSC (P = 0.011). HOXA10 expression in RNA sequencing data based on 1, 077 cases exhibited concordant significant up-regulation in NSCLC, LUAD and LUSC (P < 0.001). In NSCLC, HOXA10 expression was closely correlated to patient T stage (P = 0.006). In LUAD, HOXA10 expression was compactly correlated to patient N stage (P = 0.02). Some of the microarrays from Gene Expression Omnibus (GEO) and ArrayExpress showed consistent over-expression of HOXA10 levels in NSCLCs. More importantly, the combined SMD value was 0.052 (95% CI: 0.29-0.75, P < 0.001) generated by meta-analysis from 47 datasets based on 4, 616 cases of NSCLC. The area under the curve (AUC) of SROC supported the over-expression of HOXA10 in NSCLC as being 0.88 (95% CI: 0.81-0.93), with sensitivity and specificity of 0.88 (95% CI: 0.81-0.93) and 0.56 (95% CI: 0.44-0.66), respectively. In addition, 111 co-expressed genes were collected from cBioPortal and enriched in "cell cycle", "cell adhesion molecules", "p53 signaling", and "adherens junction". Interestingly, an up-regulation trend of HOXA10 protein expression was also observed in NSCLC through tissue chips and immunohistochemistry. In conclusion, the overexpression of HOXA10 may play a pivotal role in the tumorigenesis of NSCLC, and this effect is observed more obviously in LUSC than in LUAD.

Targeted next-generation sequencing reveals recurrence-associated genomic alterations in early-stage non-small cell lung cancer

Purpose

The identification of genomic alterations related to recurrence in early-stage non-small cell lung cancer (NSCLC) patients may help better stratify high-risk individuals and guide treatment strategies. This study aimed to identify the molecular biomarkers of recurrence in early-stage NSCLC.

Results

Of the 42 tumors evaluable for genomic alterations, TP53 and EGFR were the most frequent alterations with population frequency 52.4% and 50.0%, respectively. Fusion genes were detected in four patients, which had lower mutational burden and relatively better genomic stability. EGFR mutation and fusion gene were mutually exclusive in this study. CDKN2A, FAS, SUFU and SMARCA4 genomic alterations were only observed in the relapsed patients. Increased copy number alteration index was observed in early relapsed patients. Among these genomic alterations, early-stage NSCLCs harboring CDKN2A, FAS, SUFU and SMARCA4 genomic alterations were found to be significantly associated with recurrence. Some of these new findings were validated using The Cancer Genome Atlas (TCGA) dataset.

Conclusions

The genomic alterations of CDKN2A, FAS, SUFU and SMARCA4 in early-stage NSCLC are found to be associated with recurrence, but confirmation in a larger independent cohort is required to define the clinical impact.

Materials and Methods

Paired primary tumor and normal lung tissue samples were collected for targeted next-generation sequencing analysis. A panel targets exons for 440 genes was used to assess the mutational and copy number status of selected genes in three clinically relevant groups of stage I/II NSCLC patients: 1) Early relapse; 2) Late relapse; and 3) No relapse.

Prognostic significance of AKR1B10 in patients with resected lung adenocarcinoma

Background

Aldo‐keto reductases (AKRs) modify carbonyl groups on aldehyde or ketones to form primary or secondary alcohols, which are then conjugated with sulfates or glucuronide for excretion. The AKR1B10 gene encodes a member of the AKR superfamily. Overexpression of AKR1B10 plays an important role in the tumorigenesis of lung cancer cells; however, the prognostic value of AKR1B10 expression in patients with lung adenocarcinoma has not been well demonstrated.

Methods

A total of 96 patients with resected lung adenocarcinoma were included in the study. AKR1B10 expression was determined by immunohistochemistry in tumor specimens. The prognostic value of AKR1B10 overexpression and its relationship with clinicopathological variables were investigated.

Results

AKR1B10 overexpression was identified in 22 (22.9%) of the 96 patients and tended to be significantly associated with N1 or N2 status (P = 0.055). AKR1B10 overexpression was not a significant prognostic factor for overall survival (P = 0.301) but was a significant prognostic factor for poor recurrence‐free survival (P = 0.015). T status (T3 or T4 vs. T1 or T2; P = 0.020), N1 or N2 (vs. N0; P = 0.019), predominant pattern group (lepidic/acinar/papillary vs. micropapillary/solid; P = 0.023), and AKR1B10 overexpression (P = 0.013) were significant prognostic factors for poor recurrence‐free survival in multivariate analysis.

Conclusions

AKR1B10 overexpression was a significant prognostic factor for poor recurrence‐free survival in patients with resected lung adenocarcinoma. This information is useful to stratify patients at high‐risk of recurrence after lung adenocarcinoma resection.

Prognostic implication of aquaporin 1 overexpression in resected lung adenocarcinoma

OBJECTIVES

Aquaporins (AQPs) are a group of transmembrane water-selective channel proteins thought to play a role in the regulation of water permeability for plasma membranes. Indeed, high AQP levels have been suggested to promote the progression, invasion and metastasis of tumours. Specifically, AQP1 and AQP5 overexpression in lung adenocarcinoma (AC) have been suggested to be involved in molecular mechanisms in lung cancer. The aim of this retrospective cohort single-centre study was to assess both the levels of expression and therein the prognostic significance, regarding outcome of AQP1 and AQP5 in resected AC patients.

METHODS

Patients with histological diagnoses of lung AC submitted to pulmonary resection were included in this cohort study. Tissue microarrays containing cores from 185 ACs were prepared. AQP1 and AQP5 expressions were assessed by immunohistochemistry. Results were scored as either low (Score 0-2) or high (Score 3-9). Clinical data, pathological tumour-node-metastasis staging and follow-up were recorded. Multivariate Cox survival analysis and Fisher's t-test were performed.

RESULTS

AQP1 overexpression was detected in 85 (46%) patients, while AQP5 overexpression was observed in 45 (24%) patients. AQP1 did not result being significantly correlated with clinical and pathological parameters, while AQP5 resulted more expressed in AC with mucinous and papillary predominant patterns. Patients with AQP1 overexpression had shorter disease-free survival (P = 0.001) compared with patients without AQP1 overexpression. Multivariate analysis confirmed that AQP1 overexpression was significantly associated with shorter disease-free survival (P = 0.001).

CONCLUSIONS

Our results evidenced that AQP1 overexpression resulted in a shorter disease-free survival in lung AC patients. Being so, AQP1 overexpression might be an important prognostic marker in lung AC.

Detection of oncogenic mutations in resected bronchial margins by next-generation sequencing indicates early relapse in stage IA lung adenocarcinoma patients

Stage I non-small cell lung cancer (NSCLC) patients experience a relatively high rate of recurrence, ranging from about 30-35%. We hypothesized that this elevated risk of recurrence is due to the presence of tumor cells at bronchial margins which was undetected by conventional light microscopy. Patients with clinical stage IA (T1N0M0) NSCLC were enrolled in this study, which included 8 early-relapse(ER) and 6 no-relapse(NR) patients. Primary tumor, bronchial margin, and normal lung tissues were collected and sent to a central site for targeted next-generation sequencing analysis. All of the patients were lung adenocarcinoma. Gene mutations were identified in all tumor tissue samples (100%). Oncogenic mutations were identified in 87.5%(7/8) bronchial margins of early relapse patients, whereas only 16.7%(1/6) no-relapse (NR) patient of marginal tissue had identified gene mutation. Additionally, concordance between primary tumor and bronchial margins was relatively high, with 4 of 8 (50%) ER patients having at least one identical mutation. Moreover, according to the gene mutation status in marginal tissue, 87.5% (7/8) of patients with at least one gene mutation in the bronchial margins had local recurrence or metastasis, whereas only 16.7% (1/6) of patients without any mutation detected had signs of relapse, the recurrence rate was significantly higher than that of the negative mutation margin group ((p (log-rank) = 0.023). The existence of oncogenic mutations in bronchial margins may represent occult residual tumor and elevated risk of recurrence in early stage NSCLC patients. Thus, assessing molecular status in bronchial margins may help identify patients who might benefit from extensive surgery or adjuvant treatment.

Differences in reproductive toxicology between alopecia drugs: an analysis on adverse events among female and male cases

Alopecia is a dermatological condition with limited therapeutic options. Only two drugs, finasteride and minoxidil, are approved by FDA for alopecia treatment. However, little is known about the differences in adverse effects between these two drugs. We examined the clinical reports submitted to the FDA Adverse Event Reporting System (FAERS) from 2004 to 2014. For both female and males, finasteride was found to be more associated with reproductive toxicity as compared to minoxidil. Among male alopecia cases, finasteride was significantly more concurrent with several forms of sexual dysfunction. Among female alopecia cases, finasteride was significantly more concurrent with harm to fetus and disorder of uterus. In addition, drug-gene network analysis indicated that finasteride could profoundly disturb pathways related to sex hormone signaling and oocyte maturation. These findings could provide clues for subsequent toxicological research. Taken together, this analysis suggested that finasteride could be more liable to various reproductive adverse effects. Some of these adverse effects have yet to be warned in FDA-approved drug label. This information can help improve the treatment regimen of alopecia and post-marketing regulation of drug products.

Prediction of potential drivers connecting different dysfunctional levels in lung adenocarcinoma via a protein-protein interaction network

Lung cancer is a serious disease that threatens an affected individual's life. Its pathogenesis has not yet to be fully described, thereby impeding the development of effective treatments and preventive measures. "Cancer driver" theory considers that tumor initiation can be associated with a number of specific mutations in genes called cancer driver genes. Four omics levels, namely, (1) methylation, (2) microRNA, (3) mutation, and (4) mRNA levels, are utilized to cluster cancer driver genes. In this study, the known dysfunctional genes of these four levels were used to identify novel driver genes of lung adenocarcinoma, a subtype of lung cancer. These genes could contribute to the initiation and progression of lung adenocarcinoma in at least two levels. First, random walk with restart algorithm was performed on a protein-protein interaction (PPI) network constructed with PPI information in STRING by using known dysfunctional genes as seed nodes for each level, thereby yielding four groups of possible genes. Second, these genes were further evaluated in a test strategy to exclude false positives and select the most important ones. Finally, after conducting an intersection operation in any two groups of genes, we obtained several inferred driver genes that contributed to the initiation of lung adenocarcinoma in at least two omics levels. Several genes from these groups could be confirmed according to recently published studies. The inferred genes reported in this study were also different from those described in a previous study, suggesting that they can be used as essential supplementary data for investigations on the initiation of lung adenocarcinoma. This article is part of a Special Issue entitled: Accelerating Precision Medicine through Genetic and Genomic Big Data Analysis edited by Yudong Cai & Tao Huang.

Loss of RUNX1 is associated with aggressive lung adenocarcinomas

The mammalian runt-related factor 1 (RUNX1) is a master transcription factor that regulates lineage specification of hematopoietic stem cells. RUNX1 translocations result in the development of myeloid leukemias. Recently, RUNX1 has been implicated as a tumor suppressor in other cancers. We postulated RUNX1 expression may be associated with lung adenocarcinoma etiology and/or progression. We evaluated the association of RUNX1 mRNA expression with overall survival data from The Cancer Genome Atlas (TCGA), a publically available database. Compared to high expression levels, Low RUNX1 levels from lung adenocarcinomas were associated with a worse overall survival (Hazard Ratio = 2.014 (1.042-3.730 95% confidence interval), log-rank p = 0.035) compared to those that expressed high RUNX1 levels. Further immunohistochemical examination of 85 surgical specimens resected at the University of Vermont Medical Center identified that low RUNX1 protein expression was associated with larger tumors (p = 0.038). Gene expression network analysis was performed on the same subset of TCGA cases that demonstrated differential survival by RUNX1 expression. This analysis, which reveals regulatory relationships, showed that reduced RUNX1 levels were closely linked to upregulation of the transcription factor E2F1. To interrogate this relationship, RUNX1 was depleted in a lung cancer cell line that expresses high levels of RUNX1. Loss of RUNX1 resulted in enhanced proliferation, migration, and invasion. RUNX1 depletion also resulted in increased mRNA expression of E2F1 and multiple E2F1 target genes. Our data implicate loss of RUNX1 as driver of lung adenocarcinoma aggression, potentially through deregulation of the E2F1 pathway.

Synthesis of Potent and Selective Inhibitors of Aldo-Keto Reductase 1B10 and Their Efficacy against Proliferation, Metastasis, and Cisplatin Resistance of Lung Cancer Cells

Aldo-keto reductase 1B10 (AKR1B10) is overexpressed in several extraintestinal cancers, particularly in non-small-cell lung cancer, where AKR1B10 is a potential diagnostic marker and therapeutic target. Selective AKR1B10 inhibitors are required because compounds should not inhibit the highly related aldose reductase that is involved in monosaccharide and prostaglandin metabolism. Currently, 7-hydroxy-2-(4-methoxyphenylimino)-2H-chromene-3-carboxylic acid benzylamide (HMPC) is known to be the most potent competitive inhibitor of AKR1B10, but it is nonselective. In this study, derivatives of HMPC were synthesized by removing the 4-methoxyphenylimino moiety and replacing the benzylamide with phenylpropylamide. Among them, 4c and 4e showed higher AKR1B10 inhibitory potency (IC₅₀ 4.2 and 3.5 nM, respectively) and selectivity than HMPC. The treatments with the two compounds significantly suppressed not only migration, proliferation, and metastasis of lung cancer A549 cells but also metastatic and invasive potentials of cisplatin-resistant A549 cells.

Reverse phase protein array (RPPA) combined with computational analysis to unravel relevant prognostic factors in non- small cell lung cancer (NSCLC): a pilot study

In this work high throughput technology and computational analysis were used to study two stage IV lung adenocarcinoma patients treated with standard chemotherapy with markedly different survival (128 months vs 6 months, respectively) and whose tumor samples exhibit a dissimilar protein activation pattern of the signal transduction. Tumor samples of the two patients were subjected to Reverse Phase Protein Microarray (RPPA) analysis to explore the expression/activation levels of 51 signaling proteins. We selected the most divergent proteins based on the ratio of their RPPA values in the two patients with short (s-OS) and long (l-OS) overall survival (OS) and tested them against a EGFR-IGF1R mathematical model. The model with RPPA data showed that the activation levels of 19 proteins were different in the two patients. The four proteins that most distinguished the two patients were BADS155/136 and c-KITY703/719 having a higher activation level in the patient with short survival and p70S6S371/T389 and b-RAFS445 that had a lower activation level in the s-OS patient. The final model describes the interactions between the MAPK and PI3K-mTOR pathways, including 21 nodes. According to our model mTOR and ERK activation levels were predicted to be lower in the s-OS patient than the l-OS patient, while the AMPK activation level was higher in the s-OS patient. Moreover, KRAS activation was predicted to be higher in the l-OS KRAS-mutated patient. In accordance with their different biological properties, the Moment Independent Robustness Indicator in s-OS and l-OS predicted the interaction of MAPK and mTOR and the crosstalk AKT with p90RSK as candidates to be prognostic factors and drug targets.

Digoxin-induced anemia among patients with atrial fibrillation and heart failure: clinical data analysis and drug-gene interaction network

Digoxin is widely used to treat various heart conditions. In order to clarify the association between digoxin and anemia adverse reaction, we inspected case reports submitted to the FDA Adverse Event Reporting System (FAERS) between January 2004 and December 2015. These reports involved 75618 atrial fibrillation patients and 15699 heart failure patients. Compared to other therapies, digoxin treatment was significantly more likely to be concurrent with anemia adverse reaction among both atrial fibrillation patients (pooled OR = 1.38, 95% CI 1.14-1.68, P-value = 0.001) and heart failure patients (pooled OR =1.50, 95% CI 1.33-1.59-, P =4.27×10^-5). We further explored previously published evidences and found 821 human genes directly or indirectly interacting with digoxin. Functional analysis indicated that these genes were significantly enriched in the biological processes of iron transport, which are closely related to iron deficiency anemia. Taken together, our retrospective analysis demonstrated the significant association between digoxin treatment and anemia adverse reaction, which should be seriously considered in clinical practice. Functional enrichment analysis on digoxin-related genes warranted subsequent research on the underlying toxicological mechanisms.