Genetic architecture of alcohol consumption identified by a genotype-stratified GWAS and impact on esophageal cancer risk in Japanese people.

An East Asian-specific variant on aldehyde dehydrogenase 2 (ALDH2 rs671, G>A) is the major genetic determinant of alcohol consumption. We performed an rs671 genotype-stratified genome-wide association study meta-analysis of alcohol consumption in 175,672 Japanese individuals to explore gene-gene interactions with rs671 behind drinking behavior. The analysis identified three genome-wide significant loci (GCKR, KLB, and ADH1B) in wild-type homozygotes and six (GCKR, ADH1B, ALDH1B1, ALDH1A1, ALDH2, and GOT2) in heterozygotes, with five showing genome-wide significant interaction with rs671. Genetic correlation analyses revealed ancestry-specific genetic architecture in heterozygotes. Of the discovered loci, four (GCKR, ADH1B, ALDH1A1, and ALDH2) were suggested to interact with rs671 in the risk of esophageal cancer, a representative alcohol-related disease. Our results identify the genotype-specific genetic architecture of alcohol consumption and reveal its potential impact on alcohol-related disease risk.

Long non-coding RNA TILR constitutively represses TP53 and apoptosis in lung cancer.

Non-coding RNAs have an integral regulatory role in numerous functions related to lung cancer development. Here, we report identification of a novel lncRNA, termed TP53-inhibiting lncRNA (TILR), which was found to function as a constitutive negative regulator of p53 expression, including activation of downstream genes such as p21 and MDM2, and induction of apoptosis. A proteomic search for TILR-associated proteins revealed an association with PCBP2, while the mid-portion of TILR was found to be required for both PCBP2 and p53 mRNA binding. In addition, depletion of PCBP2 resulted in phenocopied effects of TILR silencing. TILR was also shown to suppress p53 expression in a post-transcriptional manner, as well as via a positive feedback loop involving p53 and Fanconi anemia pathway genes. Taken together, the present findings clearly demonstrate that TILR constitutively inhibits p53 expression in cooperation with PCBP2, thus maintaining p53 transcriptional activity at a level sufficiently low for avoidance of spurious apoptosis induction.

Meclozine Attenuates the MARK Pathway in Mammalian Chondrocytes and Ameliorates FGF2-Induced Bone Hyperossification in Larval Zebrafish.

Meclozine has been developed as an inhibitor of fibroblast growth factor receptor 3 (FGFR3) to treat achondroplasia (ACH). Extracellular signal regulated kinase (ERK) phosphorylation was attenuated by meclozine in FGF2-treated chondrocyte cell line, but the site of its action has not been elucidated. Although orally administered meclozine promoted longitudinal bone growth in a mouse model of ACH, its effect on craniofacial bone development during the early stage remains unknown. Herein, RNA-sequencing analysis was performed using murine chondrocytes from FGF2-treated cultured tibiae, which was significantly elongated by meclozine treatment. Gene set enrichment analysis demonstrated that FGF2 significantly increased the enrichment score of mitogen-activated protein kinase (MAPK) family signaling cascades in chondrocytes; however, meclozine reduced this enrichment. Next, we administered meclozine to FGF2-treated larval zebrafish from 8 h post-fertilization (hpf). We observed that FGF2 significantly increased the number of ossified vertebrae in larval zebrafish at 7 days post-fertilization (dpf), while meclozine delayed vertebral ossification in FGF2-induced zebrafish. Meclozine also reversed the FGF2-induced upregulation of ossified craniofacial bone area, including ceratohyal, hyomandibular, and quadrate. The current study provided additional evidence regarding the inhibitory effect of meclozine on the FGF2-induced upregulation of MAPK signaling in chondrocytes and FGF2-induced development of craniofacial and vertebral bones.

Genome-wide association meta-analysis identifies GP2 gene risk variants for pancreatic cancer.

Pancreatic cancer is the fourth leading cause of cancer-related deaths in Japan. To identify risk loci, we perform a meta-analysis of three genome-wide association studies comprising 2,039 pancreatic cancer patients and 32,592 controls in the Japanese population. Here, we identify 3 (13q12.2, 13q22.1, and 16p12.3) genome-wide significant loci (P < 5.0 × 10-8), of which 16p12.3 has not been reported in the Western population. The lead single nucleotide polymorphism (SNP) at 16p12.3 is rs78193826 (odds ratio = 1.46, 95% confidence interval = 1.29-1.66, P = 4.28 × 10-9), an Asian-specific, nonsynonymous glycoprotein 2 (GP2) gene variant. Associations between selected GP2 gene variants and pancreatic cancer are replicated in 10,822 additional cases and controls of East Asian origin. Functional analyses using cell lines provide supporting evidence of the effect of rs78193826 on KRAS activity. These findings suggest that GP2 gene variants are probably associated with pancreatic cancer susceptibility in populations of East Asian ancestry.

PTEN Loss Mediates Clinical Cross-Resistance to CDK4/6 and PI3Kα Inhibitors in Breast Cancer.

The combination of CDK4/6 inhibitors with antiestrogen therapies significantly improves clinical outcomes in ER-positive advanced breast cancer. To identify mechanisms of acquired resistance, we analyzed serial biopsies and rapid autopsies from patients treated with the combination of the CDK4/6 inhibitor ribociclib with letrozole. This study revealed that some resistant tumors acquired RB loss, whereas other tumors lost PTEN expression at the time of progression. In breast cancer cells, ablation of PTEN, through increased AKT activation, was sufficient to promote resistance to CDK4/6 inhibition in vitro and in vivo. Mechanistically, PTEN loss resulted in exclusion of p27 from the nucleus, leading to increased activation of both CDK4 and CDK2. Because PTEN loss also causes resistance to PI3Kα inhibitors, currently approved in the post-CDK4/6 setting, these findings provide critical insight into how this single genetic event may cause clinical cross-resistance to multiple targeted therapies in the same patient, with implications for optimal treatment-sequencing strategies. SIGNIFICANCE: Our analysis of serial biopsies uncovered RB and PTEN loss as mechanisms of acquired resistance to CDK4/6 inhibitors, utilized as first-line treatment for ER-positive advanced breast cancer. Importantly, these findings have near-term clinical relevance because PTEN loss also limits the efficacy of PI3Kα inhibitors currently approved in the post-CDK4/6 setting.This article is highlighted in the In This Issue feature, p. 1.

Increased Synthesis of MCL-1 Protein Underlies Initial Survival of EGFR-Mutant Lung Cancer to EGFR Inhibitors and Provides a Novel Drug Target.

Purpose: EGFR inhibitors (EGFRi) are effective against EGFR-mutant lung cancers. The efficacy of these drugs, however, is mitigated by the outgrowth of resistant cells, most often driven by a secondary acquired mutation in EGFR, T790M We recently demonstrated that T790M can arise de novo during treatment; it follows that one potential therapeutic strategy to thwart resistance would be identifying and eliminating these cells [referred to as drug-tolerant cells (DTC)] prior to acquiring secondary mutations like T790M Experimental Design: We have developed DTCs to EGFRi in EGFR-mutant lung cancer cell lines. Subsequent analyses of DTCs included RNA-seq, high-content microscopy, and protein translational assays. Based on these results, we tested the ability of MCL-1 BH3 mimetics to combine with EGFR inhibitors to eliminate DTCs and shrink EGFR-mutant lung cancer tumors in vivo Results: We demonstrate surviving EGFR-mutant lung cancer cells upregulate the antiapoptotic protein MCL-1 in response to short-term EGFRi treatment. Mechanistically, DTCs undergo a protein biosynthesis enrichment resulting in increased mTORC1-mediated mRNA translation of MCL-1, revealing a novel mechanism in which lung cancer cells adapt to short-term pressures of apoptosis-inducing kinase inhibitors. Moreover, MCL-1 is a key molecule governing the emergence of early EGFR-mutant DTCs to EGFRi, and we demonstrate it can be effectively cotargeted with clinically emerging MCL-1 inhibitors both in vitro and in vivo Conclusions: Altogether, these data reveal that this novel therapeutic combination may delay the acquisition of secondary mutations, therefore prolonging therapy efficacy. Clin Cancer Res; 24(22); 5658-72. ©2018 AACR.

Analysis of the androgen receptor-regulated lncRNA landscape identifies a role for ARLNC1 in prostate cancer progression.

The androgen receptor (AR) plays a critical role in the development of the normal prostate as well as prostate cancer. Using an integrative transcriptomic analysis of prostate cancer cell lines and tissues, we identified ARLNC1 (AR-regulated long noncoding RNA 1) as an important long noncoding RNA that is strongly associated with AR signaling in prostate cancer progression. Not only was ARLNC1 induced by the AR protein, but ARLNC1 stabilized the AR transcript via RNA-RNA interaction. ARLNC1 knockdown suppressed AR expression, global AR signaling and prostate cancer growth in vitro and in vivo. Taken together, these data support a role for ARLNC1 in maintaining a positive feedback loop that potentiates AR signaling during prostate cancer progression and identify ARLNC1 as a novel therapeutic target.

Oncogenic Role of THOR, a Conserved Cancer/Testis Long Non-coding RNA.

Large-scale transcriptome sequencing efforts have vastly expanded the catalog of long non-coding RNAs (lncRNAs) with varying evolutionary conservation, lineage expression, and cancer specificity. Here, we functionally characterize a novel ultraconserved lncRNA, THOR (ENSG00000226856), which exhibits expression exclusively in testis and a broad range of human cancers. THOR knockdown and overexpression in multiple cell lines and animal models alters cell or tumor growth supporting an oncogenic role. We discovered a conserved interaction of THOR with IGF2BP1 and show that THOR contributes to the mRNA stabilization activities of IGF2BP1. Notably, transgenic THOR knockout produced fertilization defects in zebrafish and also conferred a resistance to melanoma onset. Likewise, ectopic expression of human THOR in zebrafish accelerated the onset of melanoma. THOR represents a novel class of functionally important cancer/testis lncRNAs whose structure and function have undergone positive evolutionary selection.

Inactivating mutations and hypermethylation of the NKX2-1/TTF-1 gene in non-terminal respiratory unit-type lung adenocarcinomas.

The major driver mutations of lung cancer, EGFR mutations and EML4-ALK fusion, are mainly detected in terminal respiratory unit (TRU)-type lung adenocarcinomas, which typically show lepidic and/or papillary patterns, but are rarely associated with a solid or invasive mucinous morphology. In order to elucidate the key genetic events in non-TRU-type lung cancer, we carried out whole-exome sequencing on 43 non-TRU-type lung adenocarcinomas based on morphology (17 acinar, nine solid, and two enteric adenocarcinomas, and 15 adenocarcinomas with a mucinous morphology). Our analysis identified mutations in TP53 (16/43, 37.2%), KRAS (13/43, 30.2%), and NKX2-1/TTF-1 (7/43; 16.3%) as the top three significantly mutated genes, while the EGFR mutation was rare (1/43, 2.3%) in this cohort. Eight NKX2-1/TTF-1 mutations (five frameshift, two nonsense, and one missense) were identified, with one case harboring two distinct NKX2-1/TTF-1 mutations (one missense and one frameshift). Functional assays with the NK2 homeobox 1 (NKX2-1)/thyroid transcription factor 1 (TTF-1) mutants revealed that none of them retain the activity as a transcriptional factor. Histologically, invasive mucinous adenocarcinomas accounted for most of the NKX2-1/TTF-1 mutations (five cases), as well as one enteric and one acinar adenocarcinoma. Immunohistochemistry showed that the cohort was largely divided into TTF-1-postive/hepatocyte nuclear factor 4-α (HNF4-α)-negative and TTF-1-negative/HNF4-α-positive groups. NKX2-1/TTF-1 mutations were exclusively found in the latter, in which the gastrointestinal markers, mucin 5AC and cytokeratin 20, were frequently expressed. Bisulfite sequencing revealed that the NKX2-1/TTF-1 gene body was highly methylated in NKX2-1/TTF-1-negative cases, including those without the NKX2-1/TTF-1 mutations. The genetic or epigenetic inactivation of NKX2-1/TTF-1 may play an essential role in the development and aberrant differentiation of non-TRU-type lung adenocarcinomas.

The lncRNA landscape of breast cancer reveals a role for DSCAM-AS1 in breast cancer progression.

Molecular classification of cancers into subtypes has resulted in an advance in our understanding of tumour biology and treatment response across multiple tumour types. However, to date, cancer profiling has largely focused on protein-coding genes, which comprise <1% of the genome. Here we leverage a compendium of 58,648 long noncoding RNAs (lncRNAs) to subtype 947 breast cancer samples. We show that lncRNA-based profiling categorizes breast tumours by their known molecular subtypes in breast cancer. We identify a cohort of breast cancer-associated and oestrogen-regulated lncRNAs, and investigate the role of the top prioritized oestrogen receptor (ER)-regulated lncRNA, DSCAM-AS1. We demonstrate that DSCAM-AS1 mediates tumour progression and tamoxifen resistance and identify hnRNPL as an interacting protein involved in the mechanism of DSCAM-AS1 action. By highlighting the role of DSCAM-AS1 in breast cancer biology and treatment resistance, this study provides insight into the potential clinical implications of lncRNAs in breast cancer.

Blood-borne miRNA profile-based diagnostic classifier for lung adenocarcinoma.

Accumulated evidence indicates that various types of miRNA are aberrantly expressed in lung cancer and secreted into the bloodstream. For this study, we constructed a serum diagnostic classifier based on detailed bioinformatics analysis of miRNA profiles from a training cohort of 143 lung adenocarcinoma patients and 49 healthy subjects, resulting in a 20 miRNA-based classifier. Validation performed with an independent cohort of samples from lung adenocarcinoma patients (n = 110), healthy subjects (n = 52), and benign pulmonary disease patients (n = 47) showed a sensitivity of 89.1% and specificity of 94.9%, with an area under the curve value of 0.958. Notably, 90.8% of Stage I lung adenocarcinoma cases were correctly diagnosed. Interestingly, this classifier also detected squamous and large cell lung carcinoma cases at relatively high rates (70.4% and 70.0%, respectively), which appears to be consistent with organ site-dependent miRNA expression in cancer tissues. In contrast, we observed significantly lower rates (0-35%) using samples from 96 cases of cancer in other major organs, with breast cancer the lowest. These findings warrant a future study to realize its clinical application as a part of diagnostic procedures for lung cancers, for which early detection and surgical removal is presently the only hope for eventual cure.

KRAS Engages AGO2 to Enhance Cellular Transformation.

Oncogenic mutations in RAS provide a compelling yet intractable therapeutic target. Using co-immunoprecipitation mass spectrometry, we uncovered an interaction between RAS and Argonaute 2 (AGO2). Endogenously, RAS and AGO2 co-sediment and co-localize in the endoplasmic reticulum. The AGO2 N-terminal domain directly binds the Switch II region of KRAS, agnostic of nucleotide (GDP/GTP) binding. Functionally, AGO2 knockdown attenuates cell proliferation in mutant KRAS-dependent cells and AGO2 overexpression enhances KRAS(G12V)-mediated transformation. Using AGO2-/- cells, we demonstrate that the RAS-AGO2 interaction is required for maximal mutant KRAS expression and cellular transformation. Mechanistically, oncogenic KRAS attenuates AGO2-mediated gene silencing. Overall, the functional interaction with AGO2 extends KRAS function beyond its canonical role in signaling.

The landscape of long noncoding RNAs in the human transcriptome.

Long noncoding RNAs (lncRNAs) are emerging as important regulators of tissue physiology and disease processes including cancer. To delineate genome-wide lncRNA expression, we curated 7,256 RNA sequencing (RNA-seq) libraries from tumors, normal tissues and cell lines comprising over 43 Tb of sequence from 25 independent studies. We applied ab initio assembly methodology to this data set, yielding a consensus human transcriptome of 91,013 expressed genes. Over 68% (58,648) of genes were classified as lncRNAs, of which 79% were previously unannotated. About 1% (597) of the lncRNAs harbored ultraconserved elements, and 7% (3,900) overlapped disease-associated SNPs. To prioritize lineage-specific, disease-associated lncRNA expression, we employed non-parametric differential expression testing and nominated 7,942 lineage- or cancer-associated lncRNA genes. The lncRNA landscape characterized here may shed light on normal biology and cancer pathogenesis and may be valuable for future biomarker development.

NKX2-1/TTF-1: an enigmatic oncogene that functions as a double-edged sword for cancer cell survival and progression.

Emerging evidence indicates that NKX2-1, a homeobox-containing transcription factor also known as TTF-1, plays a role as a "lineage-survival" oncogene in lung adenocarcinomas. In T cell acute lymphoblastic leukemia, gene rearrangements lead to aberrant expression of NKX2-1/TTF-1. Despite accumulating evidence supporting its oncogenic role, it has become apparent that NKX2-1/TTF-1 expression also has biological and clinical functions in the opposite direction that act against tumor progression. Herein, we review recent findings showing these enigmatic double-edged characteristics, with special attention given to the roles of NKX2-1/TTF-1 in lung development and carcinogenesis.

MYBPH inhibits NM IIA assembly via direct interaction with NMHC IIA and reduces cell motility.

Actomyosin filament assembly is a critical step in tumor cell migration. We previously found that myosin binding protein H (MYBPH) is directly transactivated by the TTF-1 lineage-survival oncogene in lung adenocarcinomas and inhibits phosphorylation of the myosin regulatory light chain (RLC) of non-muscle myosin IIA (NM IIA) via direct interaction with Rho kinase 1 (ROCK1). Here, we report that MYBPH also directly interacts with an additional molecule, non-muscle myosin heavy chain IIA (NMHC IIA), which was found to occur between MYBPH and the rod portion of NMHC IIA. MYBPH inhibited NMHC IIA assembly and reduced cell motility. Conversely, siMYBPH-induced increased motility was partially, yet significantly, suppressed by blebbistatin, a non-muscle myosin II inhibitor, while more profound effects were attained by combined treatment with siROCK1 and blebbistatin. Electron microscopy observations showed well-ordered paracrystals of NMHC IIA reflecting an assembled state, which were significantly less frequently observed in the presence of MYBPH. Furthermore, an in vitro sedimentation assay showed that a greater amount of NMHC IIA was in an unassembled state in the presence of MYBPH. Interestingly, treatment with a ROCK inhibitor that impairs transition of NM IIA from an assembly-incompetent to assembly-competent state reduced the interaction between MYBPH and NMHC IIA, suggesting that MYBPH has higher affinity to assembly-competent NM IIA. These results suggest that MYBPH inhibits RLC and NMHC IIA, independent components of NM IIA, and negatively regulates actomyosin organization at 2 distinct steps, resulting in firm inhibition of NM IIA assembly.

NKX2-1/TITF1/TTF-1-Induced ROR1 is required to sustain EGFR survival signaling in lung adenocarcinoma.

We and others previously identified NKX2-1, also known as TITF1 and TTF-1, as a lineage-survival oncogene in lung adenocarcinomas. Here we show that NKX2-1 induces the expression of the receptor tyrosine kinase-like orphan receptor 1 (ROR1), which in turn sustains a favorable balance between prosurvival PI3K-AKT and pro-apoptotic p38 signaling, in part through ROR1 kinase-dependent c-Src activation, as well as kinase activity-independent sustainment of the EGFR-ERBB3 association, ERBB3 phosphorylation, and consequential PI3K activation. Notably, ROR1 knockdown effectively inhibited lung adenocarcinoma cell lines, irrespective of their EGFR status, including those with resistance to the EGFR tyrosine kinase inhibitor gefitinib. Our findings thus identify ROR1 as an "Achilles' heel" in lung adenocarcinoma, warranting future development of therapeutic strategies for this devastating cancer.

MYBPH, a transcriptional target of TTF-1, inhibits ROCK1, and reduces cell motility and metastasis.

Cell migration driven by actomyosin filament assembly is a critical step in tumour invasion and metastasis. Herein, we report identification of myosin binding protein H (MYBPH) as a transcriptional target of TTF-1 (also known as NKX2-1 and TITF1), a master regulator of lung development that also plays a role as a lineage-survival oncogene in lung adenocarcinoma development. MYBPH inhibited assembly competence-conferring phosphorylation of the myosin regulatory light chain (RLC) as well as activating phosphorylation of LIM domain kinase (LIMK), unexpectedly through its direct physical interaction with Rho kinase 1 (ROCK1) rather than with RLC. Consequently, MYBPH inhibited ROCK1 and negatively regulated actomyosin organization, which in turn reduced single cell motility and increased collective cell migration, resulting in decreased cancer invasion and metastasis. Finally, we also show that MYBPH is epigenetically inactivated by promoter DNA methylation in a fraction of TTF-1-positive lung adenocarcinomas, which appears to be in accordance with its deleterious functions in lung adenocarcinoma invasion and metastasis, as well as with the paradoxical association of TTF-1 expression with favourable prognosis in lung adenocarcinoma patients.

A novel network profiling analysis reveals system changes in epithelial-mesenchymal transition.

Patient-specific analysis of molecular networks is a promising strategy for making individual risk predictions and treatment decisions in cancer therapy. Although systems biology allows the gene network of a cell to be reconstructed from clinical gene expression data, traditional methods, such as bayesian networks, only provide an averaged network for all samples. Therefore, these methods cannot reveal patient-specific differences in molecular networks during cancer progression. In this study, we developed a novel statistical method called NetworkProfiler, which infers patient-specific gene regulatory networks for a specific clinical characteristic, such as cancer progression, from gene expression data of cancer patients. We applied NetworkProfiler to microarray gene expression data from 762 cancer cell lines and extracted the system changes that were related to the epithelial-mesenchymal transition (EMT). Out of 1732 possible regulators of E-cadherin, a cell adhesion molecule that modulates the EMT, NetworkProfiler, identified 25 candidate regulators, of which about half have been experimentally verified in the literature. In addition, we used NetworkProfiler to predict EMT-dependent master regulators that enhanced cell adhesion, migration, invasion, and metastasis. In order to further evaluate the performance of NetworkProfiler, we selected Krueppel-like factor 5 (KLF5) from a list of the remaining candidate regulators of E-cadherin and conducted in vitro validation experiments. As a result, we found that knockdown of KLF5 by siRNA significantly decreased E-cadherin expression and induced morphological changes characteristic of EMT. In addition, in vitro experiments of a novel candidate EMT-related microRNA, miR-100, confirmed the involvement of miR-100 in several EMT-related aspects, which was consistent with the predictions obtained by NetworkProfiler.