Search

Supplementary Data from Integrative Profiling of T790M-Negative EGFR-Mutated NSCLC Reveals Pervasive Lineage Transition and Therapeutic Opportunities

Purpose:Despite the established role of EGFR tyrosine kinase inhibitors (TKIs) in EGFR-mutated NSCLC, drug resistance inevitably ensues, with a paucity of treatment options especially in EGFRT790M-negative resistance.Experimental Design:We performed whole-exome and transcriptome analysis of 59 patients with first- and second-generation EGFR TKI-resistant metastatic EGFR-mutated NSCLC to characterize and compare molecular alterations mediating resistance in T790M-positive (T790M+) and -negative (T790M−) disease.Results:Transcriptomic analysis revealed ubiquitous loss of adenocarcinoma lineage gene expression in T790M− tumors, orthogonally validated using multiplex IHC. There was enrichment of genomic features such as TP53 alterations, 3q chromosomal amplifications, whole-genome doubling and nonaging mutational signatures in T790M− tumors. Almost half of resistant tumors were further classified as immunehot, with clinical outcomes conditional on immune cell-infiltration state and T790M status. Finally, using a Bayesian statistical approach, we explored how T790M− and T790M+ disease might be predicted using comprehensive genomic and transcriptomic profiles of treatment-naïve patients.Conclusions:Our results illustrate the interplay between genetic alterations, cell lineage plasticity, and immune microenvironment in shaping divergent TKI resistance and outcome trajectories in EGFR-mutated NSCLC. Genomic and transcriptomic profiling may facilitate the design of bespoke therapeutic approaches tailored to a tumor's adaptive potential.

Study Objectives Genetics impacts sleep, yet, the molecular mechanisms underlying sleep regulation remain elusive. In this study, we built machine learning models to predict sleep genes based on their similarity to genes that are known to regulate sleep. Methods We trained a prediction model on thousands of published datasets, representing circadian, immune, sleep deprivation, and many other processes, using a manually curated list of 109 sleep genes. Results Our predictions fit with prior knowledge of sleep regulation and identified key genes and pathways to pursue in follow-up studies. As an example, we focused on the NF-κB pathway and showed that chronic activation of NF-κB in a genetic mouse model impacted the sleep-wake patterns. Conclusion Our study highlights the power of machine learning in integrating prior knowledge and genome-wide data to study genetic regulation of complex behaviors such as sleep.

SignificanceWe successfully identified theduperallele as a null mutation of Cryptochrome 1 in Syrian hamsters. Here, we have shown the use of fast homozygosity mapping as an effective approach to identify causal mutations in mammals, despite lacking chromosomal genome information. In the course of this work, we improved the draft Syrian hamster genome and generated datasets necessary to exploit Syrian hamsters as a modern genetic research model. The unique physiological features of Syrian hamsters make them a desirable model to investigate human diseases, including circadian disorders, cancer, heart function, metabolism, and infectious diseases (e.g., severe acute respiratory syndrome coronavirus 2).

Introduction Intermittent hypoxia (IH) is a condition that occurs in people with obstructive sleep apnea (OSA), in which breathing is repeatedly interrupted during sleep. This can cause a lack of oxygen (hypoxia) in the body. OSA is associated with chronic systemic inflammation, thought to increase the risk of comorbid conditions, including cardiovascular disease. However, the mechanisms behind these pathologic consequences are not well understood. We hypothesize that IH dysregulates the circadian system and heightens inflammatory responses through hypoxia-inducible factors (HIFs), contributing to the end organ damage among patients with OSA. Here, we tested this hypothesis using a mouse model of OSA. Methods C57BL/6J background mice were exposed to normoxic or IH conditions for either seven days (short-term) or six weeks (long-term) under 12-hour light and 12-hour dark cycles. After exposure, cardiopulmonary tissues were collected and processed to evaluate circadian transcriptomes, including among HIFs. We also measured protein levels of inflammatory mediators in the heart and lung to compare the effects of short- and long-term exposure to IH. In addition, we measured the real-time circadian physiology of blood pressure, heart rate, activity, and core body temperature using a DSI-telemetry system. Results We observed dysregulated circadian systems of the lung and heart following IH exposure. IH also had a negative impact on immune responses in the heart and lung tissue, resulting in a decrease in inflammatory mediators following short-term exposure. However, long-term exposure to IH had the opposite effect, causing an overall increase in inflammatory mediators that contribute to comorbid cardiopulmonary pathologies. In addition, we also observed altered circadian physiology, particularly a significant increase in core body temperature. Conclusion Our results suggest that a dysregulated circadian system is associated with altered inflammatory responses in the heart and lung following IH exposure. Our study provides novel insight into the pathophysiological mechanisms associated with IH and potential biomarkers for the identification of comorbid cardiopulmonary pathology. Support (if any) National Institutes of Health grants 5K08HL148551; Triological Society; and the CCHMC Procter Scholar Award.

At least one-third of adults in the United States experience intermittent hypoxia (IH) due to health or living conditions. The majority of these adults suffer with sleep breathing conditions and associated circadian rhythm disorders. The impact of IH on the circadian clock is not well characterized. In the current study, we used an IH mouse model to understand the impact of IH on the circadian gene expression of the canonical clock genes in the central (the brain) and peripheral (the liver) tissues. Gene expression was measured using a Quantitative Reverse Transcription Polymerase Chain Reaction (qRT-PCR). CircaCompare was used to evaluate the differential rhythmicity between normoxia and IH. Our observations suggested that the circadian clock in the liver was less sensitive to IH compared to the circadian clock in the brain.

Lung cancer is the world’s leading cause of cancer death and shows strong ancestry disparities. By sequencing and assembling a large genomic and transcriptomic dataset of lung adenocarcinoma (LUAD) in individuals of East Asian ancestry (EAS; n = 305), we found that East Asian LUADs had more stable genomes characterized by fewer mutations and fewer copy number alterations than LUADs from individuals of European ancestry. This difference is much stronger in smokers as compared to nonsmokers. Transcriptomic clustering identified a new EAS-specific LUAD subgroup with a less complex genomic profile and upregulated immune-related genes, allowing the possibility of immunotherapy-based approaches. Integrative analysis across clinical and molecular features showed the importance of molecular phenotypes in patient prognostic stratification. EAS LUADs had better prediction accuracy than those of European ancestry, potentially due to their less complex genomic architecture. This study elucidated a comprehensive genomic landscape of EAS LUADs and highlighted important ancestry differences between the two cohorts. Genomic and transcriptomic analysis of lung adenocarcinoma (LUAD) in Asia indicates that Asian LUADs have fewer mutations, lower driver prevalence and fewer copy number alterations than European LUADs.

Genetics impacts sleep, yet, the molecular mechanisms underlying sleep regulation remain elusive. We built machine learning (ML) models to predict genes based on their similarity to known sleep genes. Using a manually curated list of 109 labeled sleep genes, we trained a prediction model on thousands of published datasets, representing circadian, immune, sleep deprivation, and many other processes. Our predictions fit with prior knowledge of sleep regulation and also identify several key genes/pathways to pursue in follow-up studies. We tested one of our findings, the NF-κB pathway, and showed that its genetic alteration affects sleep duration in mice. Our study highlights the power of ML to integrate prior knowledge and genome-wide data to study genetic regulation of sleep and other complex behaviors.

Obstructive sleep apnea (OSA) results from episodes of airway collapse and intermittent hypoxia and is associated with a host of health complications including dementia, diabetes, heart failure, and stroke. Although the lung is the first organ to sense changes in inspired oxygen levels, little is known about the consequences of IH to the lung hypoxia-inducible factor (HIF)-responsive pathways. Furthermore, cellular mechanisms causing disease progression across multiple systems in OSA are unknown. We hypothesized that exposure to IH would lead to up- and down-regulation of diverse expression pathways and that individual cell populations would show distinctive responses to IH. We identify changes in circadian and immune pathways in lungs from mice exposed to IH. Among all cell types, endothelial cells show the most prominent transcriptional changes. Interestingly, up-regulated genes in endothelial, fibroblast, and myofibroblast cells were enriched for genes associated with pulmonary fibrosis and pulmonary hypertension. These genes include targets of several drugs currently used to treat chronic pulmonary diseases. Our results reveal potential candidates for cell-targeted therapy seeking to minimize pulmonary effects of OSA. A better understanding of the pathophysiologic mechanisms underlying diseases associated with OSA could improve our therapeutic approaches, directing therapies to the most relevant cells and molecular pathways.

Various forms of immunotherapy, such as checkpoint blockade immunotherapy, are proving to be effective at restoring T cell-mediated immune responses that can lead to marked and sustained clinical responses, but only in some patients and cancer types1–4. Patients and tumours may respond unpredictably to immunotherapy partly owing to heterogeneity of the immune composition and phenotypic profiles of tumour-infiltrating lymphocytes (TILs) within individual tumours and between patients5,6. Although there is evidence that tumour-mutation-derived neoantigen-specific T cells play a role in tumour control2,4,7–10, in most cases the antigen specificities of phenotypically diverse tumour-infiltrating T cells are largely unknown. Here we show that human lung and colorectal cancer CD8+ TILs can not only be specific for tumour antigens (for example, neoantigens), but also recognize a wide range of epitopes unrelated to cancer (such as those from Epstein–Barr virus, human cytomegalovirus or influenza virus). We found that these bystander CD8+ TILs have diverse phenotypes that overlap with tumour-specific cells, but lack CD39 expression. In colorectal and lung tumours, the absence of CD39 in CD8+ TILs defines populations that lack hallmarks of chronic antigen stimulation at the tumour site, supporting their classification as bystanders. Expression of CD39 varied markedly between patients, with some patients having predominantly CD39− CD8+ TILs. Furthermore, frequencies of CD39 expression among CD8+ TILs correlated with several important clinical parameters, such as the mutation status of lung tumour epidermal growth factor receptors. Our results demonstrate that not all tumour-infiltrating T cells are specific for tumour antigens, and suggest that measuring CD39 expression could be a straightforward way to quantify or isolate bystander T cells. Human lung and colorectal tumours contain a population of tumour-infiltrating lymphocytes that are specific for tumour-unrelated antigens and, unlike tumour-antigen-specific tumour-infiltrating lymphocytes, do not express CD39.

Purpose Mesenchymal epithelial transition factor ( MET) activation has been implicated as an oncogenic driver in epidermal growth factor receptor ( EGFR)–mutant non–small-cell lung cancer (NSCLC) and can mediate primary and secondary resistance to EGFR tyrosine kinase inhibitors (TKI). High copy number thresholds have been suggested to enrich for response to MET inhibitors. We examined the clinical relevance of MET copy number gain (CNG) in the setting of treatment-naive metastatic EGFR-mutant–positive NSCLC. Patients and Methods MET fluorescence in situ hybridization was performed in 200 consecutive patients identified as metastatic treatment-naïve EGFR-mutant–positive. We defined MET-high as CNG greater than or equal to 5, with an additional criterion of MET/centromeric portion of chromosome 7 ratiο greater than or equal to 2 for amplification. Time-to-treatment failure (TTF) to EGFR TKI in patients identified as MET-high and -low was estimated by Kaplan-Meier method and compared using log-rank test. Multiregion single-nucleotide polymorphism array analysis was performed on 13 early-stage resected EGFR-mutant–positive NSCLC across 59 sectors to investigate intratumoral heterogeneity of MET CNG. Results Fifty-two (26%) of 200 patients in the metastatic cohort were MET-high at diagnosis; 46 (23%) had polysomy and six (3%) had amplification. Median TTF was 12.2 months (95% CI, 5.7 to 22.6 months) versus 13.1 months (95% CI, 10.6 to 15.0 months) for MET-high and -low, respectively ( P = .566), with no significant difference in response rate regardless of copy number thresholds. Loss of MET was observed in three of six patients identified as MET-high who underwent postprogression biopsies, which is consistent with marked intratumoral heterogeneity in MET CNG observed in early-stage tumors. Suboptimal response (TTF, 1.0 to 6.4 months) to EGFR TKI was observed in patients with coexisting MET amplification (five [3.2%] of 154). Conclusion Although up to 26% of TKI-naïve EGFR-mutant–positive NSCLC harbor high MET CNG by fluorescence in situ hybridization, this did not significantly affect response to TKI, except in patients identified as MET-amplified. Our data underscore the limitations of adopting arbitrary copy number thresholds and the need for cross-assay validation to define therapeutically tractable MET pathway dysregulation in EGFR-mutant–positive NSCLC.

Lung cancer is the world's leading cause of cancer death and shows strong ancestry disparities. By sequencing and assembling a large genomic and transcriptomic dataset of lung adenocarcinoma (LUAD) in individuals of East Asian ancestry (EAS; n = 305), we found that East Asian LUADs had more stable genomes characterized by fewer mutations and fewer copy number alterations than LUADs from individuals of European ancestry. This difference is much stronger in smokers as compared to nonsmokers. Transcriptomic clustering identified a new EAS-specific LUAD subgroup with a less complex genomic profile and upregulated immune-related genes, allowing the possibility of immunotherapy-based approaches. Integrative analysis across clinical and molecular features showed the importance of molecular phenotypes in patient prognostic stratification. EAS LUADs had better prediction accuracy than those of European ancestry, potentially due to their less complex genomic architecture. This study elucidated a comprehensive genomic landscape of EAS LUADs and highlighted important ancestry differences between the two cohorts.

The discovery that half of the mammalian protein-coding genome is regulated by the circadian clock has clear implications for medicine. Recent studies demonstrated that the circadian clock influences therapeutic outcomes in human heart disease and cancer. However, biological time is rarely given clinical consideration. A key barrier is the absence of information on tissue-specific molecular rhythms in the human body. We have applied the cyclic ordering by periodic structure (CYCLOPS) algorithm, designed to reconstruct sample temporal order in the absence of time-of-day information, to the gene expression collection of 13 tissues from 632 human donors. We identified rhythms in gene expression across the body; nearly half of protein-coding genes were shown to be cycling in at least 1 of the 13 tissues analyzed. One thousand of these cycling genes encode proteins that either transport or metabolize drugs or are themselves drug targets. These results provide a useful resource for studying the role of circadian rhythms in medicine and support the idea that biological time might play a role in determining drug response.

Lung cancer has the highest cancer associated mortality rate in many countries across the world. In contrast to western populations, approximately half of lung adenocarcinoma cases in Singapore harbour activating epidermal growth factor receptor (EGFR) mutations, with preponderance for never smokers and female gender. Although EGFR tyrosine kinase inhibitors (TKIs) confer high response rates of up to 70%, drug resistance invariably ensues - most commonly through the “acquisition” of EGFR T790M mutation. While the extent and pattern of intratumoral heterogeneity (ITH) in non-small cell lung cancer (NSCLC) were recently described, these studies examined histologically and molecularly diverse cohort of patients, majority being current or ex-smokers. Here we report ITH in eight never-smoker EGFR mutant lung adenocarcinoma cases of Asian ethnicity. All eight patients had no prior treatment history and harboured an activating EGFR mutation (5 L858R , 2 exon 19 deletion, 1 exon 20 insertion). They underwent lobectomy for Stage IA, IB NSCLC. Tumors were harvested using a systematic sectoring protocol according to standard operation procedures, with tissue banked for exome sequencing, RNA-sequencing and SNP array. A total of 46 tumor sectors (at least 4 regions from each of the 8 tumors) were subject to whole exome sequencing, with matched normal samples. With an average sequencing depth of 100x, we identified 860 somatic exonic SNVs (601 being non-synonymous) and 49 indels across all samples. The median number of SNVs per patient was 112 and per sector was 49. Notably, activating EGFR mutations were identified across all tumor sectors of all but two patients (for whom it was identified in 3 of 5 and 5 of 7 sectors respectively). In addition, we did not identify the EGFR T790M mutation in any of the sequenced tumor sectors, suggesting that, this resistance mutation is not present at detectable frequencies even as a branch or subclonal event in a treatment naïve scenario. Of 20 genes that were significantly mutated across 46 individual tumor sectors, only two overlapped with published recurrently mutated genes in NSCLC. In conclusion, we show that activating EGFR mutations are ubiquitous truncal events in 6 of 8 Asian never-smoker lung adenocarcinoma – consistent with its role as a therapeutically tractable driver gene. The T790M mutation commonly associated with TKI resistance was not detected as a subclonal event in treatment naïve patients. Further, our study reveals the unique mutation spectra of Asian EGFR mutant lung adenocarcinoma, highlighting the value of multi-region sequencing in characterising the genomic architecture of defined molecular subsets of NSCLC from different ethnic backgrounds. Citation Format: Rahul Nahar, Weiwei Zhai, Angela Takano, Alexis Jiaying Khng, Xingliang Liu, Chong Hee Lim, Audrey S.M. Teo, Cheryl Xueli Chan, Apoorva Gogna, Kiat-Hon Lim, Tina Koh, Huay Mei Poh, Yin Yeng Lee, Liang He Chen, Tong Zhang, Vidhya Gomathi Krishnan, N Gopalakrishna Iyer, Pauline Ng, Wan Teck Lim, Bing Lim, Eng-Huat Tan, Daniel S.W. Tan, Axel M. Hillmer. Intratumor heterogeneity in never-smoker Asian EGFR mutant lung adenocarcinoma. [abstract]. In: Proceedings of the AACR Special Conference on Computational and Systems Biology of Cancer; Feb 8-11 2015; San Francisco, CA. Philadelphia (PA): AACR; Cancer Res 2015;75(22 Suppl 2):Abstract nr B2-54.

Lung cancer has the highest cancer associated mortality rate in many countries across the world. In contrast to western populations, approximately half of lung adenocarcinoma cases in Singapore harbour activating epidermal growth factor receptor (EGFR) mutations, with preponderance for never smokers and female gender. Although EGFR tyrosine kinase inhibitors (TKIs) confer high response rates of up to 70%, drug resistance invariably ensues - most commonly through the “acquisition” of EGFR T790M mutation. While the extent and pattern of intratumoral heterogeneity (ITH) in non-small cell lung cancer (NSCLC) were recently described, these studies examined histologically and molecularly diverse cohort of patients, majority being current or ex-smokers. Here we report ITH in eight never-smoker EGFR mutant lung adenocarcinoma cases of Asian ethnicity. All eight patients had no prior treatment history and harboured an activating EGFR mutation (5 L858R, 2 exon 19 deletion, 1 exon 20 insertion). They underwent lobectomy for Stage IA, IB NSCLC. Tumors were harvested using a systematic sectoring protocol according to standard operation procedures, with tissue banked for exome sequencing, RNA-sequencing and SNP array. A total of 46 tumor sectors (at least 4 regions from each of the 8 tumors) were subject to whole exome sequencing, with matched normal samples. With an average sequencing depth of 100x, we identified 860 somatic exonic SNVs (601 being non-synonymous) and 49 indels across all samples. The median number of SNVs per patient was 112 and per sector was 49. Notably, activating EGFR mutations were identified across all tumor sectors of all but two patients (for whom it was identified in 3 of 5 and 5 of 7 sectors respectively). In addition, we did not identify the EGFR T790M mutation in any of the sequenced tumor sectors, suggesting that, this resistance mutation is not present at detectable frequencies even as a branch or subclonal event in a treatment naïve scenario. Of 20 genes that were significantly mutated across 46 individual tumor sectors, only two overlapped with published recurrently mutated genes in NSCLC. In conclusion, we show that activating EGFR mutations are ubiquitous truncal events in 6 of 8 Asian never-smoker lung adenocarcinoma – consistent with its role as a therapeutically tractable driver gene. The T790M mutation commonly associated with TKI resistance was not detected as a subclonal event in treatment naïve patients. Further, our study reveals the unique mutation spectra of Asian EGFR mutant lung adenocarcinoma, highlighting the value of multi-region sequencing in characterising the genomic architecture of defined molecular subsets of NSCLC from different ethnic backgrounds. Citation Format: Rahul Nahar, Weiwei Zhai, Angela Takano, Alexis Jiaying Khng, Xingliang Liu, Chong Hee Lim, Audrey S.M. Teo, Cheryl Xueli Chan, Apoorva Gogna, Kiat-Hon Lim, Tina Koh, Huay Mei Poh, Yin Yeng Lee, Liang He Chen, Tong Zhang, Vidhya Gomathi Krishnan, N Gopalakrishna Iyer, Pauline Ng, Wan Teck Lim, Bing Lim, Eng-Huat Tan, Daniel S.W. Tan, Axel M. Hillmer. Intratumor heterogeneity in never-smoker Asian EGFR mutant lung adenocarcinoma. [abstract]. In: Proceedings of the AACR Special Conference on Translation of the Cancer Genome; Feb 7-9, 2015; San Francisco, CA. Philadelphia (PA): AACR; Cancer Res 2015;75(22 Suppl 1):Abstract nr A1-25.

Various forms of immunotherapy, such as checkpoint blockade immunotherapy, are proving to be effective at restoring T cell-mediated immune responses that can lead to marked and sustained clinical responses, but only in some patients and cancer types

Gastrodia elata (Tianma) is a traditional Chinese medicine often used for the treatment of headache, convulsions, hypertension, and cardiovascular diseases. The vasodilatory actions of Tianma led us to investigate its specific effects on memory and learning as well as on Alzheimer's disease (AD)-related signaling. We conducted a radial arm water maze analysis and the novel object recognition test to assess the cognitive functions of Tianma-treated mice. Our data show that Tianma enhances cognitive functions in mice. Further investigations revealed that Tianma enhances the α-secretase-mediated proteolytic processing of the amyloid precursor protein (App) that precludes the amyloid-β peptide production and supports the non-amyloidogenic processing of App which is favorable in AD treatment. We hypothesize that Tianma promotes cognitive functions and neuronal survival by inhibiting β-site App-cleaving enzyme 1 activity and promoting the neuroprotective α-secretase activity.

SUMMARY Taming cell-to-cell variability in gene expression is critical for precise pattern formation during embryonic development. To investigate the source and buffering mechanism of expression variability, we studied a biological clock, the vertebrate segmentation clock, controlling the precise spatiotemporal patterning of the vertebral column. By counting single transcripts of segmentation clock genes in zebrafish, we show that clock genes have low RNA amplitudes and expression variability is primarily driven by gene extrinsic sources, which is suppressed by Notch signaling. We further show that expression noise surprisingly increases from the posterior progenitor zone to the anterior segmentation and differentiation zone. Our computational model reproduces the spatial noise profile by incorporating spatially increasing time delays in gene expression. Our results, suggesting that expression variability is controlled by the balance of time delays and cell signaling in a vertebrate tissue, will shed light on the accuracy of natural clocks in multi-cellular systems and inspire engineering of robust synthetic oscillators., In Brief Keskin et al. show that segmentation clock transcription levels display low amplitude and high heterogeneity. Clock gene expression variability is primarily driven by gene extrinsic sources, which is suppressed by Notch signaling. Gene expression noise increases along the posteroanterior axis. Spatial gradients of time delays contribute to the noise gradient along the axis.