Treatment of non-small cell lung cancer (NSCLC) with radiotherapy or chemoradiotherapy is often accompanied by the development of esophagitis and pneumonitis. Identifying patients who might be at increased risk for normal tissue toxicity would help in determination of the optimal radiation dose to avoid these events. We profiled 59 single nucleotide polymorphisms (SNPs) from 37 inflammation-related genes in 173 NSCLC patients with stage IIIA/IIIB (dry) disease who were treated with definitive radiation or chemoradiation. For esophagitis risk, nine SNPs were associated with a 1.5- to 4-fold increase in risk, including three PTGS2 (COX2) variants: rs20417 (HR:1.93, 95% CI:1.10-3.39), rs5275 (HR:1.58, 95% CI:1.09-2.27), and rs689470 (HR:3.38, 95% CI:1.09-10.49). Significantly increased risk of pneumonitis was observed for patients with genetic variation in the proinflammatory genes IL1A, IL8, TNF, TNFRSF1B, and MIF. In contrast, NOS3:rs1799983 displayed a protective effect with a 45% reduction in pneumonitis risk (HR:0.55, 95% CI:0.31-0.96). Pneumonitis risk was also modulated by polymorphisms in anti-inflammatory genes, including genetic variation in IL13. rs20541 and rs180925 each resulted in increased risk (HR:2.95, 95% CI:1.14-7.63 and HR:3.23, 95% CI:1.03-10.18, respectively). The cumulative effect of these SNPs on risk was dose-dependent, as evidenced by a significantly increased risk of either toxicity with an increasing number of risk genotypes (P<0.001). These results suggest that genetic variations among inflammation pathway genes may modulate the development of radiation-induced toxicity and, ultimately, help in identifying patients who are at an increased likelihood for such events.

Non-small cell lung cancer (NSCLC) represents the majority (85%) of lung cancers and is comprised mainly of adenocarcinomas and squamous cell carcinomas (SCCs). The sequential pathogenesis of lung adenocarcinomas and SCCs occurs through dissimilar phases as the former tumors typically arise in the lung periphery whereas the latter normally arise near the central airway.

Activating mutations in one allele of an oncogene (heterozygous mutations) are widely believed to be sufficient for tumorigenesis. However, mutant allele specific imbalance (MASI) has been observed in tumors and cell lines harboring mutations of oncogenes.

Traditionally, non-small cell lung cancer is treated as a single disease entity in terms of systemic therapy. Emerging evidence suggests the major subtypes--adenocarcinoma (AC) and squamous cell carcinoma (SqCC)--respond differently to therapy. Identification of the molecular differences between these tumor types will have a significant impact in designing novel therapies that can improve the treatment outcome.

Generating a balanced network of inhibitory and excitatory neurons during development requires precise transcriptional control. In the dorsal spinal cord, Ptf1a, a basic helix-loop-helix (bHLH) transcription activator, maintains this delicate balance by inducing homeodomain (HD) transcription factors such as Pax2 to specify the inhibitory lineage while suppressing HD factors such as Tlx1/3 that specify the excitatory lineage. We uncover the mechanism by which Ptf1a represses excitatory cell fate in the inhibitory lineage. We identify Prdm13 as a direct target of Ptf1a and reveal that Prdm13 actively represses excitatory cell fate by binding to regulatory sequences near the Tlx1 and Tlx3 genes to silence their expression. Prdm13 acts through multiple mechanisms, including interactions with the bHLH factor Ascl1, to repress Ascl1 activation of Tlx3. Thus, Prdm13 is a key component of a highly coordinated transcriptional network that determines the balance of inhibitory versus excitatory neurons in the dorsal spinal cord.

KRAS is one of the most commonly mutated oncogenes in human cancer. Mutant KRAS aberrantly regulates metabolic networks. However, the contribution of cellular metabolism to mutant KRAS tumorigenesis is not completely understood. We report that mutant KRAS regulates intracellular fatty acid metabolism through Acyl-coenzyme A (CoA) synthetase long-chain family member 3 (ACSL3), which converts fatty acids into fatty Acyl-CoA esters, the substrates for lipid synthesis and β-oxidation. ACSL3 suppression is associated with depletion of cellular ATP and causes the death of lung cancer cells. Furthermore, mutant KRAS promotes the cellular uptake, retention, accumulation, and β-oxidation of fatty acids in lung cancer cells in an ACSL3-dependent manner. Finally, ACSL3 is essential for mutant KRAS lung cancer tumorigenesis in vivo and is highly expressed in human lung cancer. Our data demonstrate that mutant KRAS reprograms lipid homeostasis, establishing a metabolic requirement that could be exploited for therapeutic gain.

It is well established that lung tumors induce the formation of lymphatic vessels. However, the molecular mechanisms controlling tumor lymphangiogenesis in lung cancer have not been fully delineated. In the present study, we identify a panel of non-small cell lung cancer (NSCLC) cell lines that induce lymphangiogenesis and use genome-wide mRNA expression to characterize the molecular mechanisms regulating tumor lymphangiogenesis. We show that Calu-1, H1993, HCC461, HCC827, and H2122 NSCLC cell lines form tumors that induce lymphangiogenesis whereas Calu-3, H1155, H1975, and H2073 NSCLC cell lines form tumors that do not induce lymphangiogenesis. By analyzing genome-wide mRNA expression data, we identify a 17-gene expression signature that distinguishes lymphangiogenic from non-lymphangiogenic NSCLC cell lines. Importantly, VEGF-C is the only lymphatic growth factor in this expression signature and is approximately 50-fold higher in the lymphangiogenic group than in the non-lymphangiogenic group. We show that forced expression of VEGF-C by H1975 cells induces lymphangiogenesis and that knockdown of VEGF-C in H1993 cells inhibits lymphangiogenesis. Additionally, we demonstrate that the triple angiokinase inhibitor, nintedanib (small molecule that blocks all FGFRs, PDGFRs, and VEGFRs), suppresses tumor lymphangiogenesis in H1993 tumors. Together, these data suggest that VEGF-C is the dominant driver of tumor lymphangiogenesis in NSCLC and reveal a specific therapy that could potentially block tumor lymphangiogenesis in NSCLC patients.

Lung cancer is a highly heterogeneous disease in terms of both underlying genetic lesions and response to therapeutic treatments. We performed deep whole-genome sequencing and transcriptome sequencing on 19 lung cancer cell lines and three lung tumor/normal pairs. Overall, our data show that cell line models exhibit similar mutation spectra to human tumor samples. Smoker and never-smoker cancer samples exhibit distinguishable patterns of mutations. A number of epigenetic regulators, including KDM6A, ASH1L, SMARCA4, and ATAD2, are frequently altered by mutations or copy number changes. A systematic survey of splice-site mutations identified 106 splice site mutations associated with cancer specific aberrant splicing, including mutations in several known cancer-related genes. RAC1b, an isoform of the RAC1 GTPase that includes one additional exon, was found to be preferentially up-regulated in lung cancer. We further show that its expression is significantly associated with sensitivity to a MAP2K (MEK) inhibitor PD-0325901. Taken together, these data present a comprehensive genomic landscape of a large number of lung cancer samples and further demonstrate that cancer-specific alternative splicing is a widespread phenomenon that has potential utility as therapeutic biomarkers. The detailed characterizations of the lung cancer cell lines also provide genomic context to the vast amount of experimental data gathered for these lines over the decades, and represent highly valuable resources for cancer biology.

Neural networks are balanced by inhibitory and excitatory neuronal activity. The formation of these networks is initially generated through neuronal subtype specification controlled by transcription factors. The basic helix-loop-helix (bHLH) transcription factor Ptf1a is essential for the generation of GABAergic inhibitory neurons in the dorsal spinal cord, cerebellum, and retina. The transcription factor Rbpj is a transducer of the Notch signaling pathway that functions to maintain neural progenitor cells. Here we demonstrate Ptf1a and Rbpj interact in a complex that is required in vivo for specification of the GABAergic neurons, a function that cannot be substituted by the classical form of the bHLH heterodimer with E-protein or Notch signaling through Rbpj. We show that a mutant form of Ptf1a without the ability to bind Rbpj, while retaining its ability to interact with E-protein, is incapable of inducing GABAergic (Pax2)- and suppressing glutamatergic (Tlx3)-expressing cells in the chick and mouse neural tube. Moreover, we use an Rbpj conditional mutation to demonstrate that Rbpj function is essential for GABAergic specification, and that this function is independent of the Notch signaling pathway. Together, these findings demonstrate the requirement for a Ptf1a-Rbpj complex in controlling the balanced formation of inhibitory and excitatory neurons in the developing spinal cord, and point to a novel Notch-independent function for Rbpj in nervous system development.

RUVBL1 and RUVBL2 (collectively RUVBL1/2) are essential AAA+ ATPases that function as co-chaperones and have been implicated in cancer. Here we investigated the molecular and phenotypic role of RUVBL1/2 ATPase activity in non-small cell lung cancer (NSCLC). We find that RUVBL1/2 are overexpressed in NSCLC patient tumors, with high expression associated with poor survival. Utilizing a specific inhibitor of RUVBL1/2 ATPase activity, we show that RUVBL1/2 ATPase activity is necessary for the maturation or dissociation of the PAQosome, a large RUVBL1/2-dependent multiprotein complex. We also show that RUVBL1/2 have roles in DNA replication, as inhibition of its ATPase activity can cause S-phase arrest, which culminates in cancer cell death via replication catastrophe. While in vivo pharmacological inhibition of RUVBL1/2 results in modest antitumor activity, it synergizes with radiation in NSCLC, but not normal cells, an attractive property for future preclinical development.

The mechanics of the cellular microenvironment continuously modulates cell functions such as growth, survival, apoptosis, differentiation and morphogenesis via cytoskeletal remodelling and actomyosin contractility1-3. Although all of these processes consume energy4,5, it is unknown whether and how cells adapt their metabolic activity to variable mechanical cues. Here we report that the transfer of human bronchial epithelial cells from stiff to soft substrates causes a downregulation of glycolysis via proteasomal degradation of the rate-limiting metabolic enzyme phosphofructokinase (PFK). PFK degradation is triggered by the disassembly of stress fibres, which releases the PFK-targeting E3 ubiquitin ligase tripartite motif (TRIM)-containing protein 21 (TRIM21). Transformed non-small-cell lung cancer cells, which maintain high glycolytic rates regardless of changing environmental mechanics, retain PFK expression by downregulating TRIM21, and by sequestering residual TRIM21 on a stress-fibre subset that is insensitive to substrate stiffness. Our data reveal a mechanism by which glycolysis responds to architectural features of the actomyosin cytoskeleton, thus coupling cell metabolism to the mechanical properties of the surrounding tissue. These processes enable normal cells to tune energy production in variable microenvironments, whereas the resistance of the cytoskeleton in response to mechanical cues enables the persistence of high glycolytic rates in cancer cells despite constant alterations of the tumour tissue.

Small cell lung cancer (SCLC) is classified as a high-grade neuroendocrine (NE) tumor, but a subset of SCLC has been termed "variant" due to the loss of NE characteristics. In this study, we computed NE scores for patient-derived SCLC cell lines and xenografts, as well as human tumors. We aligned NE properties with transcription factor-defined molecular subtypes. Then we investigated the different immune phenotypes associated with high and low NE scores. We found repression of immune response genes as a shared feature between classic SCLC and pulmonary neuroendocrine cells of the healthy lung. With loss of NE fate, variant SCLC tumors regain cell-autonomous immune gene expression and exhibit higher tumor-immune interactions. Pan-cancer analysis revealed this NE lineage-specific immune phenotype in other cancers. Additionally, we observed MHC I re-expression in SCLC upon development of chemoresistance. These findings may help guide the design of treatment regimens in SCLC.

Detection of somatic mosaicism in non-proliferative cells is a new challenge in genome research, however, the accuracy of current detection strategies remains uncertain due to the lack of a ground truth. Herein, we sought to present a set of ultra-deep sequenced WES data based on reference standards generated by cell line mixtures, providing a total of 386,613 mosaic single-nucleotide variants (SNVs) and insertion-deletion mutations (INDELs) with variant allele frequencies (VAFs) ranging from 0.5% to 56%, as well as 35,113,417 non-variant and 19,936 germline variant sites as a negative control. The whole reference standard set mimics the cumulative aspect of mosaic variant acquisition such as in the early developmental stage owing to the progressive mixing of cell lines with established genotypes, ultimately unveiling 741 possible inter-sample relationships with respect to variant sharing and asymmetry in VAFs. We expect that our reference data will be essential for optimizing the current use of mosaic variant detection strategies and for developing algorithms to enable future improvements.

DNA methylation is a key regulator of gene expression and a clinical therapeutic predictor. We examined global DNA methylation beyond the generally used promoter areas in human small cell lung cancer (SCLC) and find that gene body methylation is a robust positive predictor of gene expression. Combining promoter and gene body methylation better predicts gene expression than promoter methylation alone including genes involved in the neuroendocrine classification of SCLC and the expression of therapeutically relevant genes including MGMT, SLFN11, and DLL3. Importantly, for super-enhancer (SE) covered genes such as NEUROD1 or MYC, using H3K27ac and NEUROD1, ASCL1, and POU2F3 ChIP-seq data, we show that genic methylation is inversely proportional to expression, thus providing a new approach to identify potential SE regulated genes involved in SCLC pathogenesis. To advance SCLC transitional research, these data are integrated into our web portal (https://discover.nci.nih.gov/SclcCellMinerCDB/) for open and easy access to basic and clinical investigators.

Autophagy is elevated in metastatic tumors and is often associated with active epithelial-to-mesenchymal transition (EMT). However, the extent to which EMT is dependent on autophagy is largely unknown. This study aimed to identify the mechanisms by which autophagy facilitates EMT.

Atezolizumab (anti-PD-L1), combined with carboplatin and etoposide (CE), is now a standard of care for extensive-stage small-cell lung cancer (ES-SCLC). A clearer understanding of therapeutically relevant SCLC subsets could identify rational combination strategies and improve outcomes. We conduct transcriptomic analyses and non-negative matrix factorization on 271 pre-treatment patient tumor samples from IMpower133 and identify four subsets with general concordance to previously reported SCLC subtypes (SCLC-A, -N, -P, and -I). Deeper investigation into the immune heterogeneity uncovers two subsets with differing neuroendocrine (NE) versus non-neuroendocrine (non-NE) phenotypes, demonstrating immune cell infiltration hallmarks. The NE tumors with low tumor-associated macrophage (TAM) but high T-effector signals demonstrate longer overall survival with PD-L1 blockade and CE versus CE alone than non-NE tumors with high TAM and high T-effector signal. Our study offers a clinically relevant approach to discriminate SCLC patients likely benefitting most from immunotherapies and highlights the complex mechanisms underlying immunotherapy responses.

Phosphatidylserine (PS) is an anionic phospholipid maintained on the inner-leaflet of the cell membrane and is externalized in malignant cells. We previously launched a careful unbiased selection targeting biomolecules (e.g. protein, lipid or carbohydrate) distinct to cancer cells by exploiting HCC4017 lung cancer and HBEC30KT normal epithelial cells derived from the same patient, identifying HCC4017 specific peptide-peptoid hybrid PPS1. In this current study, we identified PS as the target of PPS1. We validated direct PPS1 binding to PS using ELISA-like assays, lipid dot blot and liposome based binding assays. In addition, PPS1 recognized other negatively charged and cancer specific lipids such as phosphatidic acid, phosphatidylinositol and phosphatidylglycerol. PPS1 did not bind to neutral lipids such as phosphatidylethanolamine found in cancer and phosphatidylcholine and sphingomyelin found in normal cells. Further we found that the dimeric version of PPS1 (PPS1D1) displayed strong cytotoxicity towards lung cancer cell lines that externalize PS, but not normal cells. PPS1D1 showed potent single agent anti-tumor activity and enhanced the efficacy of docetaxel in mice bearing H460 lung cancer xenografts. Since PS and anionic phospholipid externalization is common across many cancer types, PPS1 may be an alternative to overcome limitations of protein targeted agents.

Lung cancer is the leading cause of cancer-related death. Despite a number of studies that have provided prognostic biomarkers for lung cancer, a paucity of reliable markers and therapeutic targets exist to diagnose and treat this aggressive disease. In this study we investigated the potential of nuclear receptors (NRs), many of which are well-established drug targets, as therapeutic markers in lung cancer. Using quantitative real-time PCR, we analyzed the expression of the 48 members of the NR superfamily in a human panel of 55 normal and lung cancer cell lines. Unsupervised cluster analysis of the NR expression profile segregated normal from tumor cell lines and grouped lung cancers according to type (i.e. small vs. non-small cell lung cancers). Moreover, we found that the NR signature was 79% accurate in diagnosing lung cancer incidence in smokers (n = 129). Finally, the evaluation of a subset of NRs (androgen receptor, estrogen receptor, vitamin D receptor, and peroxisome proliferator-activated receptor-γ) demonstrated the therapeutic potential of using NR expression to predict ligand-dependent growth responses in individual lung cancer cells. Preclinical evaluation of one of these receptors (peroxisome proliferator activated receptor-γ) in mouse xenografts confirmed that ligand-dependent inhibitory growth responses in lung cancer can be predicted based on a tumor's receptor expression status. Taken together, this study establishes NRs as theragnostic markers for predicting lung cancer incidence and further strengthens their potential as therapeutic targets for individualized treatment.

NSCLC (non-small cell lung cancer) often exhibits resistance to paclitaxel treatment. Identifying the elements regulating paclitaxel response will advance efforts to overcome such resistance in NSCLC therapy. Using in vitro approaches, we demonstrated that over-expression of the microRNA miR-337-3p sensitizes NCI-H1155 cells to paclitaxel, and that miR-337-3p mimic has a general effect on paclitaxel response in NSCLC cell lines, which may provide a novel adjuvant strategy to paclitaxel in the treatment of lung cancer. By combining in vitro and in silico approaches, we identified STAT3 and RAP1A as direct targets that mediate the effect of miR-337-3p on paclitaxel sensitivity. Further investigation showed that miR-337-3p mimic also sensitizes cells to docetaxel, another member of the taxane family, and that STAT3 levels are significantly correlated with taxane resistance in lung cancer cell lines, suggesting that endogenous STAT3 expression is a determinant of intrinsic taxane resistance in lung cancer. The identification of a miR-337-3p as a modulator of cellular response to taxanes, and STAT3 and RAP1A as regulatory targets which mediate that response, defines a novel regulatory pathway modulating paclitaxel sensitivity in lung cancer cells, which may provide novel adjuvant strategies along with paclitaxel in the treatment of lung cancer and may also provide biomarkers for predicting paclitaxel response in NSCLC.

TUSC2-defective gene expression is detected in the majority of lung cancers and is associated with worse overall survival. We analyzed the effects of TUSC2 re-expression on tumor cell sensitivity to the AKT inhibitor, MK2206, and explored their mutual signaling connections, in vitro and in vivo. TUSC2 transient expression in three LKB1-defective non-small cell lung cancer (NSCLC) cell lines combined with MK2206 treatment resulted in increased repression of cell viability and colony formation, and increased apoptotic activity. In contrast, TUSC2 did not affect the response to MK2206 treatment for two LKB1-wild type NSCLC cell lines. In vivo, TUSC2 systemic delivery, by nanoparticle gene transfer, combined with MK2206 treatment markedly inhibited growth of tumors in a human LKB1-defective H322 lung cancer xenograft mouse model. Biochemical analysis showed that TUSC2 transient expression in LKB1-defective NSCLC cells significantly stimulated AMP-activated protein kinase (AMPK) phosphorylation and enzymatic activity. More importantly, AMPK gene knockdown abrogated TUSC2-MK2206 cooperation, as evidenced by reduced sensitivity to the combined treatment. Together, TUSC2 re-expression and MK2206 treatment was more effective in inhibiting the phosphorylation and kinase activities of AKT and mTOR proteins than either single agent alone. In conclusion, these findings support the hypothesis that TUSC2 expression status is a biological variable that potentiates MK2206 sensitivity in LKB1-defective NSCLC cells, and identifies the AMPK/AKT/mTOR signaling axis as an important regulator of this activity.