Thomas Muley

Identification of recurrent FGFR3 fusion genes in lung cancer through kinome‐centred RNA sequencing

Oncogenic fusion genes that involve kinases have proven to be effective targets for therapy in a wide range of cancers. Unfortunately, the diagnostic approaches required to identify these events are struggling to keep pace with the diverse array of genetic alterations that occur in cancer. Diagnostic screening in solid tumours is particularly challenging, as many fusion genes occur with a low frequency. To overcome these limitations, we developed a capture enrichment strategy to enable high‐throughput transcript sequencing of the human kinome. This approach provides a global overview of kinase fusion events, irrespective of the identity of the fusion partner. To demonstrate the utility of this system, we profiled 100 non‐small cell lung cancers and identified numerous genetic alterations impacting fibroblast growth factor receptor 3 (FGFR3) in lung squamous cell carcinoma and a novel ALK fusion partner in lung adenocarcinoma.

Non-small cell lung cancer is characterized by dramatic changes in phospholipid profiles

Non‐small cell lung cancer (NSCLC) is the leading cause of cancer death globally. To develop better diagnostics and more effective treatments, research in the past decades has focused on identification of molecular changes in the genome, transcriptome, proteome, and more recently also the metabolome. Phospholipids, which nevertheless play a central role in cell functioning, remain poorly explored. Here, using a mass spectrometry (MS)‐based phospholipidomics approach, we profiled 179 phospholipid species in malignant and matched non‐malignant lung tissue of 162 NSCLC patients (73 in a discovery cohort and 89 in a validation cohort). We identified 91 phospholipid species that were differentially expressed in cancer versus non‐malignant tissues. Most prominent changes included a decrease in sphingomyelins (SMs) and an increase in specific phosphatidylinositols (PIs). Also a decrease in multiple phosphatidylserines (PSs) was observed, along with an increase in several phosphatidylethanolamine (PE) and phosphatidylcholine (PC) species, particularly those with 40 or 42 carbon atoms in both fatty acyl chains together. 2D‐imaging MS of the most differentially expressed phospholipids confirmed their differential abundance in cancer cells. We identified lipid markers that can discriminate tumor versus normal tissue and different NSCLC subtypes with an AUC (area under the ROC curve) of 0.999 and 0.885, respectively. In conclusion, using both shotgun and 2D‐imaging lipidomics analysis, we uncovered a hitherto unrecognized alteration in phospholipid profiles in NSCLC. These changes may have important biological implications and may have significant potential for biomarker development.

Low Input Whole-Exome Sequencing to Determine the Representation of the Tumor Exome in Circulating DNA of Non-Small Cell Lung Cancer Patients.

Circulating cell-free DNA (cfDNA) released from cancerous tissues has been found to harbor tumor-associated alterations and to represent the molecular composition of the tumor. Recent advances in technologies, especially in next-generation sequencing, enable the analysis of low amounts of cfDNA from body fluids. We analyzed the exomes of tumor tissue and matched serum samples to investigate the molecular representation of the tumor exome in cfDNA. To this end, we implemented a workflow for sequencing of cfDNA from low serum volumes (200 μl) and performed whole-exome sequencing (WES) of serum and matched tumor tissue samples from six non-small cell lung cancer (NSCLC) patients and two control sera. Exomes, including untranslated regions (UTRs) of cfDNA were sequenced with an average coverage of 68.5x. Enrichment efficiency, target coverage, and sequencing depth of cfDNA reads were comparable to those from matched tissues. Discovered variants were compared between serum and tissue as well as to the COSMIC database of known mutations. Although not all tissue variants could be confirmed in the matched serum, up to 57% of the tumor variants were reflected in matched cfDNA with mutations in PIK3CA, ALK, and PTEN as well as variants at COSMIC annotated sites in all six patients analyzed. Moreover, cfDNA revealed a mutation in MTOR, which was not detected in the matched tissue, potentially from an untested region of the heterogeneous primary tumor or from a distant metastatic clone. WES of cfDNA may provide additional complementary molecular information about clinically relevant mutations and the clonal heterogeneity of the tumors.

Mutation analysis of circulating plasma DNA to determine response to EGFR tyrosine kinase inhibitor therapy of lung adenocarcinoma patients.

Long-lasting success in lung cancer therapy using tyrosine kinase inhibitors (TKIs) is rare since the tumors develop resistance due to the occurrence of molecularly altered subclones. The aim of this study was to monitor tumors over time based on the quantity of mutant plasma DNA and to identify early indications for therapy response and tumor progression. Serial plasma samples from lung adenocarcinoma patients treated with TKIs were used to quantify EGFR and KRAS mutations in circulating DNA by digital PCR. Mutant DNA levels were compared with the courses of responses to treatment with TKIs, conventional chemotherapy, radiotherapy, or combinations thereof. Variations in plasma DNA mutation levels over time were found in 15 patients. We categorize three major courses: First, signs of therapy response are associated with a fast clearing of plasma DNA mutations within a few days. Second, periods of stable disease are accompanied by either absence of mutations or fluctuation at low levels. Finally, dramatic increase of mutational load is followed by rapid tumor progression and poor patient survival. In summary, the serial assessment of EGFR mutations in the plasma of NSCLC patients allows conclusions about controlled disease and tumor progression earlier than currently available methods.

AURKA, DLGAP5, TPX2, KIF11 and CKAP5: Five specific mitosis-associated genes correlate with poor prognosis for non-small cell lung cancer patients

The growth of a tumor depends to a certain extent on an increase in mitotic events. Key steps during mitosis are the regulated assembly of the spindle apparatus and the separation of the sister chromatids. The microtubule-associated protein Aurora kinase A phosphorylates DLGAP5 in order to correctly segregate the chromatids. Its activity and recruitment to the spindle apparatus is regulated by TPX2. KIF11 and CKAP5 control the correct arrangement of the microtubules and prevent their degradation. In the present study, we investigated the role of these five molecules in non-small cell lung cancer (NSCLC). We analyzed the expression of the five genes in a large cohort of NSCLC patients (n=362) by quantitative real-time PCR. Each of the genes was highly overexpressed in the tumor tissues compared to corresponding normal lung tissue. The correlation of the expression of the individual genes depended on the histology. An increased expression of AURKA, DLGAP5, TPX2, KIF11 and CKAP5 was associated with poor overall survival (P=0.001–0.065). AURKA was a significant prognostic marker using multivariate analyses (P=0.006). Immunofluorescence studies demonstrated that the five mitosis-associated proteins co-localized with the spindle apparatus during cell division. Taken together, our data demonstrate that the expression of the mitosis-associated genes AURKA, DLGAP5, TPX2, KIF11 and CKAP5 is associated with the prognosis of NSCLC patients.

Gene signature driving invasive mucinous adenocarcinoma of the lung

Though invasive mucinous adenocarcinoma of the lung (IMA) is pathologically distinctive, the molecular mechanism driving IMA is not well understood, which hampers efforts to identify therapeutic targets. Here, by analyzing gene expression profiles of human and mouse IMA, we identified a Mucinous Lung Tumor Signature of 143 genes, which was unexpectedly enriched in mucin‐producing gastrointestinal, pancreatic, and breast cancers. The signature genes included transcription factors FOXA3, SPDEF, HNF4A, mucins MUC5AC, MUC5B, MUC3, and an inhibitory immune checkpoint VTCN1/B7‐H4 (but not PD‐L1/B7‐H1). Importantly, induction of FOXA3 or SPDEF along with mutant KRAS in lung epithelium was sufficient to develop benign or malignant mucinous lung tumors, respectively, in transgenic mice. FOXA3 and SPDEF induced MUC5AC and MUC5B, while HNF4A induced MUC3 in human mucinous lung cancer cells harboring a KRAS mutation. ChIP‐seq combined with CRISPR/Cas9 determined that upstream enhancer regions of the mucin genes MUC5AC and MUC5B, which were bound by SPDEF, were required for the expression of the mucin genes. Here, we report the molecular signature and gene regulatory network driving mucinous lung tumors.

Glycodelin: A New Biomarker with Immunomodulatory Functions in Non–Small Cell Lung Cancer

Purpose: In recent years, immune therapeutic strategies against non–small cell lung cancer (NSCLC) based on tissue-derived biomarkers, for example PD1/PD-L1 (CD274), have evolved as novel and promising treatment options. However, the crosstalk between tumor and immune cells is poorly understood. Glycodelin (gene name PAEP), initially described in the context of pregnancy and trophoblastic implantation, is a secreted immunosuppressive glycoprotein with an as-of-yet largely unknown function in lung cancer. Experimental Design: In this study, we characterized the expression and role of glycodelin in NSCLC through mRNA and protein expression analyses, functional knockdown experiments, and correlations with clinicopathologic parameters. Results: Glycodelin mRNA expression was significantly elevated in tumors (n = 336) compared with matched normal tissue (P < 0.0001). Overall survival (OS) was significantly reduced in NSCLC with high glycodelin mRNA levels in women but not in men. Glycodelin was detected in the sera of patients, and the levels correlated with recurrence and metastatic disease. Knockdown of glycodelin with siRNAs in NSCLC cell lines resulted in significant upregulation of immune system modulatory factors such as PDL1, CXCL5, CXCL16, MICA/B, and CD83 as well as proliferation stimulators EDN1 and HBEGF. Furthermore, decreased migration of tumor cells was observed. Conclusions: Altogether, the comprehensive characterization of glycodelin in NSCLC provides strong support for its use as a biomarker with immune modulatory function. Clin Cancer Res; 21(15); 3529–40. ©2015 AACR.

Concomitant expression of far upstream element (FUSE) binding protein (FBP) interacting repressor (FIR) and its splice variants induce migration and invasion of non‐small cell lung cancer (NSCLC) cells

Transcription factors integrate a variety of oncogenic input information, facilitate tumour growth and cell dissemination, and therefore represent promising therapeutic target structures. Because over‐expression of DNA‐interacting far upstream element binding protein (FBP) supports non‐small cell lung cancer (NSCLC) migration, we asked whether its repressor, FBP‐interacting repressor (FIR) is functionally inactivated and how FIR might affect NSCLC cell biology. Different FIR splice variants were highly expressed in the majority of NSCLCs, with the highest levels in tumours carrying genomic gains of chromosome 8q24.3, which contained the FIR gene locus. Nuclear FIR expression was significantly enriched at the invasion front of primary NSCLCs, but this did not correlate with tumour cell proliferation. FIR accumulation was associated with worse patient survival and tumour recurrence; in addition, FIR over‐expression significantly correlated with lymph node metastasis in squamous cell carcinomas (SCCs). In vitro, we applied newly developed methods and modelling approaches for the quantitative and time‐resolved description of the pro‐migratory and pro‐invasive capacities of SCC cells. siRNA‐mediated silencing of all FIR variants significantly reduced the speed and directional movement of tumour cells in all phases of migration. Furthermore, sprouting efficiency and single cell invasiveness were diminished following FIR inhibition. Interestingly, the silencing of FIR isoforms lacking exon 2 (FIRΔexon2) alone was sufficient to reduce lateral migration and invasion. In summary, by using scale‐spanning data derived from primary human tissues, quantitative cellular analyses and mathematical modelling, we have demonstrated that concomitant over‐expression of FIR and its splice variants drives NSCLC migration and dissemination. Copyright

Genetic changes of non-small cell lung cancer under neoadjuvant therapy

Background Large scale sequencing efforts defined common molecular alterations in non-small cell lung cancer (NSCLC) and revealed potentially druggable mutations. Yet, systematic data on the changes of the respective molecular profiles under standard therapy in NSCLC are limited. Results 14 out of 68 observed coding mutations (21%) and 6 out of 33 (18%) copy number variations (CNV) were lost or gained during therapy. Mutational and CNV changes clustered in 6/37 (16%) and 3/37 (8%) patients. Changes in clinically relevant mutations were rare but present in single cases for genes such as BRAF and PIK3CA. The type of radiochemotherapy but not the duration of therapy impacted on the frequency of mutational changes. Methods We established a lung cancer specific next-generation sequencing panel covering ~7500 hotspots of 41 genes frequently mutated in NSCLC and performed ultradeep multigene sequencing of 37 corresponding pre- and post-therapeutic formalin fixed paraffin-embedded specimens to discover mutational changes and copy number variations under neo-adjuvant radio- (RTX) and/or chemotherapy (CTX). Conclusion We unraveled changes in common driver gene candidates in NSCLC under neo-adjuvant therapy. Our data shed first light on the genetic changes of NSCLC under conventional therapy and might be taken into account when the relevance of sequential biopsy approaches is discussed.

Identification of lung cancer histology-specific variants applying Bayesian framework variant prioritization approaches within the TRICL and ILCCO consortia

Large-scale genome-wide association studies (GWAS) have likely uncovered all common variants at the GWAS significance level. Additional variants within the suggestive range (0.0001> P > 5×10(-8)) are, however, still of interest for identifying causal associations. This analysis aimed to apply novel variant prioritization approaches to identify additional lung cancer variants that may not reach the GWAS level. Effects were combined across studies with a total of 33456 controls and 6756 adenocarcinoma (AC; 13 studies), 5061 squamous cell carcinoma (SCC; 12 studies) and 2216 small cell lung cancer cases (9 studies). Based on prior information such as variant physical properties and functional significance, we applied stratified false discovery rates, hierarchical modeling and Bayesian false discovery probabilities for variant prioritization. We conducted a fine mapping analysis as validation of our methods by examining top-ranking novel variants in six independent populations with a total of 3128 cases and 2966 controls. Three novel loci in the suggestive range were identified based on our Bayesian framework analyses: KCNIP4 at 4p15.2 (rs6448050, P = 4.6×10(-7)) and MTMR2 at 11q21 (rs10501831, P = 3.1×10(-6)) with SCC, as well as GAREM at 18q12.1 (rs11662168, P = 3.4×10(-7)) with AC. Use of our prioritization methods validated two of the top three loci associated with SCC (P = 1.05×10(-4) for KCNIP4, represented by rs9799795) and AC (P = 2.16×10(-4) for GAREM, represented by rs3786309) in the independent fine mapping populations. This study highlights the utility of using prior functional data for sequence variants in prioritization analyses to search for robust signals in the suggestive range.