Integrative genomic and transcriptomic profiling of pulmonary sarcomatoid carcinoma identifies molecular subtypes associated with distinct immune features and clinical outcomes.

Background: Pulmonary sarcomatoid carcinoma (PSC) is a rare and aggressive subtype of non-small cell lung cancer (NSCLC), characterized by the presence of epithelial and sarcoma-like components. The molecular and immune landscape of PSC has not been well defined. Methods: Multiomics profiling of 21 pairs of PSCs with matched normal lung tissues was performed through targeted high-depth DNA panel, whole-exome, and RNA sequencing. We describe molecular and immune features that define subgroups of PSC with disparate genomic and immunogenic features as well as distinct clinical outcomes. Results: In total, 27 canonical cancer gene mutations were identified, with TP53 the most frequently mutated gene, followed by KRAS. Interestingly, most TP53 and KRAS mutations were earlier genomic events mapped to the trunks of the tumors, suggesting branching evolution in most PSC tumors. We identified two distinct molecular subtypes of PSC, driven primarily by immune infiltration and signaling. The Immune High (IM-H) subtype was associated with superior survival, highlighting the impact of immune infiltration on the biological and clinical features of localized PSCs. Conclusions: We provided detailed insight into the mutational landscape of PSC and identified two molecular subtypes associated with prognosis. IM-H tumors were associated with favorable recurrence-free survival and overall survival, highlighting the importance of tumor immune infiltration in the biological and clinical features of PSCs.

Impact of Cancer Therapy on Clonal Hematopoiesis Mutations and Subsequent Clinical Outcomes.

Exposure to cancer therapies is associated with an increased risk of clonal hematopoiesis (CH). The objective of our study was to investigate the genesis and evolution of CH following cancer therapy. In this prospective study, we undertook error-corrected duplex DNA sequencing in blood samples collected prior to and at two timepoints following chemoradiation in patients with esophageal or lung cancer recruited from 2013-2018. We applied a customized workflow to identify the earliest changes in CH mutation count and clone size and determine their association with clinical outcomes. Our study included 29 patients (87 samples). Their median age was 67 years, 76% (n = 22) were male; the median follow-up period was 3.9 years. The most mutated genes were DNMT3A, TET2, TP53, and ASXL1. We observed a two-fold increase in the number of mutations from before to after treatment in TP53, which differed from all other genes examined (P < .001). Among mutations detected before and after treatment, we observed an increased clone size in 38% and a decreased clone size in 5% of TP53 mutations (odds ratio = 3.7; 95% CI = 1.75-7.84; P < .001). Changes in mutation count and clone size were not observed in other genes. Individuals with an increase in the number of TP53 mutations following chemoradiation experienced shorter overall survival (hazard ratio = 7.07; 95% CI = 1.50-33.46; P = .014). In summary, we found an increase in the number and size of TP53 CH clones following chemoradiation that were associated with clinical outcomes.

Cytokine/Chemokine Release Patterns and Transcriptomic Profiles of LPS/IFNγ-Activated Human Macrophages Differentiated with Heat-Killed Mycobacterium obuense, M-CSF, or GM-CSF.

Macrophages (Mφs) are instrumental regulators of the immune response whereby they acquire diverse functional phenotypes following their exposure to microenvironmental cues that govern their differentiation from monocytes and their activation. The complexity and diversity of the mycobacterial cell wall have empowered mycobacteria with potent immunomodulatory capacities. A heat-killed (HK) whole-cell preparation of Mycobacterium obuense (M. obuense) has shown promise as an adjunctive immunotherapeutic agent for the treatment of cancer. Moreover, HK M. obuense has been shown to trigger the differentiation of human monocytes into a monocyte-derived macrophage (MDM) type named Mob-MDM. However, the transcriptomic profile and functional properties of Mob-MDMs remain undefined during an activation state. Here, we characterized cytokine/chemokine release patterns and transcriptomic profiles of lipopolysaccharide (LPS)/interferon γ (IFNγ)-activated human MDMs that were differentiated with HK M. obuense (Mob-MDM(LPS/IFNγ)), macrophage colony-stimulating factor M-MDM(LPS/IFNγ)), or granulocyte/macrophage colony-stimulating factor (GM-MDM(LPS/IFNγ)). Mob-MDM(LPS/IFNγ) demonstrated a unique cytokine/chemokine release pattern (interleukin (IL)-10low, IL-12/23p40low, IL-23p19/p40low, chemokine (C-x-C) motif ligand (CXCL)9low) that was distinct from those of M-MDM(LPS/IFNγ) and GM-MDM(LPS/IFNγ). Furthermore, M-MDM(LPS/IFNγ) maintained IL-10 production at significantly higher levels compared to GM-MDM(LPS/IFNγ) and Mob-MDM(LPS/IFNγ) despite being activated with M1-Mφ-activating stimuli. Comparative RNA sequencing analysis pointed to a distinct transcriptome profile for Mob-MDM(LPS/IFNγ) relative to both M-MDM(LPS/IFNγ) and GM-MDM(LPS/IFNγ) that comprised 417 transcripts. Functional gene-set enrichment analysis revealed significant overrepresentation of signaling pathways and biological processes that were uniquely related to Mob-MDM(LPS/IFNγ). Our findings lay a foundation for the potential integration of HK M. obuense in specific cell-based immunotherapeutic modalities such as adoptive transfer of Mφs (Mob-MDM(LPS/IFNγ)) for cancer treatment.

Augmented Lipocalin-2 Is Associated with Chronic Obstructive Pulmonary Disease and Counteracts Lung Adenocarcinoma Development.

Rationale: Early pathogenesis of lung adenocarcinoma (LUAD) remains largely unknown. We found that, relative to wild-type littermates, the innate immunomodulator Lcn2 (lipocalin-2) was increased in normal airways from mice with knockout of the airway lineage gene Gprc5a (Gprc5a-/-) and that are prone to developing inflammation and LUAD. Yet, the role of LCN2 in lung inflammation and LUAD is poorly understood.Objectives: Delineate the role of Lcn2 induction in LUAD pathogenesis.Methods: Normal airway brushings, uninvolved lung tissues, and tumors from Gprc5a-/- mice before and after tobacco carcinogen exposure were analyzed by RNA sequencing. LCN2 mRNA was analyzed in public and in-house data sets of LUAD, lung squamous cancer (LUSC), chronic obstructive pulmonary disease (COPD), and LUAD/LUSC with COPD. LCN2 protein was immunohistochemically analyzed in a tissue microarray of 510 tumors. Temporal lung tumor development, gene expression programs, and host immune responses were compared between Gprc5a-/- and Gprc5a-/-/Lcn2-/- littermates.Measurements and Main Results:Lcn2 was progressively elevated during LUAD development and positively correlated with proinflammatory cytokines and inflammation gene sets. LCN2 was distinctively elevated in human LUADs, but not in LUSCs, relative to normal lungs and was associated with COPD among smokers and patients with LUAD. Relative to Gprc5a-/- mice, Gprc5a-/-/Lcn2-/- littermates exhibited significantly increased lung tumor development concomitant with reduced T-cell abundance (CD4+) and richness, attenuated antitumor immune gene programs, and increased immune cell expression of protumor inflammatory cytokines.Conclusions: Augmented LCN2 expression is a molecular feature of COPD-associated LUAD and counteracts LUAD development in vivo by maintaining antitumor immunity.

Driver Mutations in Normal Airway Epithelium Elucidate Spatiotemporal Resolution of Lung Cancer.

Rationale: Uninvolved normal-appearing airway epithelium has been shown to exhibit specific mutations characteristic of nearby non-small cell lung cancers (NSCLCs). Yet, its somatic mutational landscape in patients with early-stage NSCLC is unknown.Objectives: To comprehensively survey the somatic mutational architecture of the normal airway epithelium in patients with early-stage NSCLC.Methods: Multiregion normal airways, comprising tumor-adjacent small airways, tumor-distant large airways, nasal epithelium and uninvolved normal lung (collectively airway field), matched NSCLCs, and blood cells (n = 498) from 48 patients were interrogated for somatic single-nucleotide variants by deep-targeted DNA sequencing and for chromosomal allelic imbalance events by genome-wide genotype array profiling. Spatiotemporal relationships between the airway field and NSCLCs were assessed by phylogenetic analysis.Measurements and Main Results: Genomic airway field carcinogenesis was observed in 25 cases (52%). The airway field epithelium exhibited a total of 269 somatic mutations in most patients (n = 36) including key drivers that were shared with the NSCLCs. Allele frequencies of these acquired variants were overall higher in NSCLCs. Integrative analysis of single-nucleotide variants and allelic imbalance events revealed driver genes with shared "two-hit" alterations in the airway field (e.g., TP53, KRAS, KEAP1, STK11, and CDKN2A) and those with single hits progressing to two in the NSCLCs (e.g., PIK3CA and NOTCH1).Conclusions: Tumor-adjacent and tumor-distant normal-appearing airway epithelia exhibit somatic driver alterations that undergo selection-driven clonal expansion in NSCLC. These events offer spatiotemporal insights into the development of NSCLC and, thus, potential targets for early treatment.

Genomic landscape of allelic imbalance in premalignant atypical adenomatous hyperplasias of the lung.

BACKGROUND: Genomic investigation of atypical adenomatous hyperplasia (AAH), the only known precursor lesion to lung adenocarcinomas (LUAD), presents challenges due to the low mutant cell fractions. This necessitates sensitive methods for detection of chromosomal aberrations to better study the role of critical alterations in early lung cancer pathogenesis and the progression from AAH to LUAD. METHODS: We applied a sensitive haplotype-based statistical technique to detect chromosomal alterations leading to allelic imbalance (AI) from genotype array profiling of 48 matched normal lung parenchyma, AAH and tumor tissues from 16 stage-I LUAD patients. To gain insights into shared developmental trajectories among tissues, we performed phylogenetic analyses and integrated our results with point mutation data, highlighting significantly-mutated driver genes in LUAD pathogenesis. FINDINGS: AI was detected in nine AAHs (56%). Six cases exhibited recurrent loss of 17p. AI and the enrichment of 17p events were predominantly identified in patients with smoking history. Among the nine AAH tissues with detected AI, seven exhibited evidence for shared chromosomal aberrations with matched LUAD specimens, including losses harboring tumor suppressors on 17p, 8p, 9p, 9q, 19p, and gains encompassing oncogenes on 8q, 12p and 1q. INTERPRETATION: Chromosomal aberrations, particularly 17p loss, appear to play critical roles early in AAH pathogenesis. Genomic instability in AAH, as well as truncal chromosomal aberrations shared with LUAD, provide evidence for mutation accumulation and are suggestive of a cancerized field contributing to the clonal selection and expansion of these premalignant lesions. FUND: Supported in part by Cancer Prevention and Research Institute of Texas (CPRIT) grant RP150079 (PS and HK), NIH grant R01HG005859 (PS) and The University of Texas MD Anderson Cancer Center Core Support Grant.

Genome-Wide Gene Expression Changes in the Normal-Appearing Airway during the Evolution of Smoking-Associated Lung Adenocarcinoma.

Smoking perpetuates in cytologically normal airways a molecular "field of injury" that is pertinent to lung cancer and early detection. The evolution of airway field changes prior to lung oncogenesis is poorly understood largely due to the long latency of lung cancer in smokers. Here, we studied airway expression changes prior to lung cancer onset in mice with knockout of the Gprc5a gene (Gprc5a-/-) and tobacco carcinogen (NNK) exposure and that develop the most common type of lung cancer, lung adenocarcinoma, within 6 months following exposure. Airway epithelial brushings were collected from Gprc5a-/- mice before exposure and at multiple times post-NNK until time of lung adenocarcinoma development and then analyzed by RNA sequencing. Temporal airway profiles were identified by linear models and analyzed by comparative genomics in normal airways of human smokers with and without lung cancer. We identified significantly altered profiles (n = 926) in the NNK-exposed mouse normal airways relative to baseline epithelia, a subset of which were concordantly modulated with smoking status in the human airway. Among airway profiles that were significantly modulated following NNK, we found that expression changes (n = 22) occurring as early as 2 months following exposure were significantly associated with lung cancer status when examined in airways of human smokers. Furthermore, a subset of a recently reported human bronchial gene classifier (Percepta; n = 56) was enriched in the temporal mouse airway profiles. We underscore evolutionarily conserved profiles in the normal-appearing airway that develop prior to lung oncogenesis and that comprise viable markers for early lung cancer detection in suspect smokers. Cancer Prev Res; 11(4); 237-48. ©2018 AACR.

Strategies for identification of somatic variants using the Ion Torrent deep targeted sequencing platform.

BACKGROUND: 'Next-generation' (NGS) sequencing has wide application in medical genetics, including the detection of somatic variation in cancer. The Ion Torrent-based (IONT) platform is among NGS technologies employed in clinical, research and diagnostic settings. However, identifying mutations from IONT deep sequencing with high confidence has remained a challenge. We compared various computational variant-calling methods to derive a variant identification pipeline that may improve the molecular diagnostic and research utility of IONT. RESULTS: Using IONT, we surveyed variants from the 409-gene Comprehensive Cancer Panel in whole-section tumors, intra-tumoral biopsies and matched normal samples obtained from frozen tissues and blood from four early-stage non-small cell lung cancer (NSCLC) patients. We used MuTect, Varscan2, IONT's proprietary Ion Reporter, and a simple subtraction we called "Poor Man's Caller." Together these produced calls at 637 loci across all samples. Visual validation of 434 called variants was performed, and performance of the methods assessed individually and in combination. Of the subset of inspected putative variant calls (n=223) in genomic regions that were not intronic or intergenic, 68 variants (30%) were deemed valid after visual inspection. Among the individual methods, the Ion Reporter method offered perhaps the most reasonable tradeoffs. Ion Reporter captured 83% of all discovered variants; 50% of its variants were visually validated. Aggregating results from multiple packages offered varied improvements in performance. CONCLUSIONS: Overall, Ion Reporter offered the most attractive performance among the individual callers. This study suggests combined strategies to maximize sensitivity and positive predictive value in variant calling using IONT deep sequencing.

TBX2 subfamily suppression in lung cancer pathogenesis: a high-potential marker for early detection.

The TBX2 subfamily (TBXs 2, 3, 4 and 5) transactivates or represses genes involved in lung organogenesis. Yet TBX2 subfamily expression in pathogenesis of non-small cell lung cancer (NSCLC), the most common lung malignancy, remains elusive. We sought to probe the expression profile of the TBX2 subfamily in early phases of NSCLC. Expression of TBX2 subfamily was analyzed in datasets of pan-normal specimens as well as NSCLCs and normal lung tissues. TBX2 subfamily expression in matched normal lungs, premalignant hyperplasias and NSCLCs was profiled by transcriptome sequencing. TBX2 subfamily expression was evaluated in the cancerization field consisting of matched NSCLCs and adjacent cytologically-normal airways relative to distant normal lungs and in a dataset of normal bronchial samples from smokers with indeterminate nodules suspicious for malignancy. Statistical analysis was performed using R. TBX2 subfamily expression was markedly elevated in normal lungs relative to other organ-specific normal tissues. Expression of the TBXs was significantly suppressed in NSCLCs relative to normal lungs (P < 10-9). TBX2 subfamily was significantly progressively decreased across premalignant lesions and NSCLCs relative to normal lungs (P < 10-4). The subfamily was significantly suppressed in NSCLCs and adjacent normal-appearing airways relative to distant normal lung tissues (P < 10-15). Further, suppressed TBX2 subfamily expression in normal bronchi was associated with lung cancer status (P < 10-5) in smokers. Our findings suggest that the TBX2 subfamily is notably suppressed in human NSCLC pathogenesis and may serve as a high-potential biomarker for early lung cancer detection in high-risk smokers.

Defining Genome-Wide Expression and Phenotypic Contextual Cues in Macrophages Generated by Granulocyte/Macrophage Colony-Stimulating Factor, Macrophage Colony-Stimulating Factor, and Heat-Killed Mycobacteria.

Heat-killed (HK) Mycobacterium obuense (NCTC13365) is currently being evaluated in the clinic as an immunotherapeutic agent for cancer treatment. Yet, the molecular underpinnings underlying immunomodulatory properties of HK M. obuense are still largely undefined. To fill this void, we sought to perform immunophenotyping, chemokine/cytokine release analysis and genome-wide characterization of monocyte-derived macrophages (MDM) in which monocytes were originally isolated from healthy donors and differentiated by HK M. obuense (Mob-MDM) relative to macrophage colony-stimulating factor (M-MDM) and granulocyte/macrophage colony-stimulating factor (GM-MDM). Immunophenotyping and cytokine release analysis revealed downregulated surface expression of CD36, decreased spontaneous release of CCL2 and increased spontaneous secretion of CCL5, CXCL8/IL-8, IL-6, and TNF-α in Mob-MDM relative to M-MDM and GM-MDM. Analysis of cytostatic activity showed that Mob-MDM exhibited similar growth inhibitory effects on immortalized and malignant epithelial cells compared with GM-MDM but at an elevated rate relative to M-MDM. To understand global cues in Mob-MDM, we performed comparative RNA-sequencing (RNA-Seq) analysis of Mob-MDM relative to GM-MDM and M-MDM (n = 4 donors). Clustering analysis underscored expression profiles (n = 256) that were significantly modulated in Mob-MDM versus both M-MDM and GM-MDM including, among others, chemokines/cytokines and their receptors, enzymes and transcriptions factors. Topological functional analysis of these profiles identified pathways and gene sets linked to Mob-MDM phenotype including nitric oxide production, acute phase response signaling and microbe recognition pathways as well as signaling cues mediated by the proinflammatory cytokine, interferon-gamma, and the intracellular pattern recognition receptor, nucleotide-binding oligomerization domain-containing protein 2. Taken together, our study highlights molecular immune phenotypes and global signaling cues in Mob-MDM that may underlie immunomodulatory properties of HK M. obuense. Such properties could be of valuable use in immunotherapy approaches such as adoptive cell therapy against cancer.

Genomic Landscape of Atypical Adenomatous Hyperplasia Reveals Divergent Modes to Lung Adenocarcinoma.

There is a dearth of knowledge about the pathogenesis of premalignant lung lesions, especially for atypical adenomatous hyperplasia (AAH), the only known precursor for the major lung cancer subtype adenocarcinoma (LUAD). In this study, we performed deep DNA and RNA sequencing analyses of a set of AAH, LUAD, and normal tissues. Somatic BRAF variants were found in AAHs from 5 of 22 (23%) patients, 4 of 5 of whom had matched LUAD with driver EGFR mutations. KRAS mutations were present in AAHs from 4 of 22 (18%) of patients. KRAS mutations in AAH were only found in ever-smokers and were exclusive to BRAF-mutant cases. Integrative analysis revealed profiles expressed in KRAS-mutant cases (UBE2C, REL) and BRAF-mutant cases (MAX) of AAH, or common to both sets of cases (suppressed AXL). Gene sets associated with suppressed antitumor (Th1; IL12A, GZMB) and elevated protumor (CCR2, CTLA-4) immune signaling were enriched in AAH development and progression. Our results reveal potentially divergent BRAF or KRAS pathways in AAH as well as immune dysregulation in the pathogenesis of this premalignant lung lesion. Cancer Res; 77(22); 6119-30. ©2017 AACR.

Development of Kras mutant lung adenocarcinoma in mice with knockout of the airway lineage-specific gene Gprc5a.

Despite the urgency for prevention and treatment of lung adenocarcinoma (LUAD), we still do not know drivers in pathogenesis of the disease. Earlier work revealed that mice with knockout of the G-protein coupled receptor Gprc5a develop late onset lung tumors including LUADs. Here, we sought to further probe the impact of Gprc5a expression on LUAD pathogenesis. We first surveyed GPRC5A expression in human tissues and found that GPRC5A was markedly elevated in human normal lung relative to other normal tissues and was consistently downregulated in LUADs. In sharp contrast to wild-type littermates, Gprc5a-/- mice treated chronically with the nicotine-specific carcinogen NNK developed LUADs by 6 months following NNK exposure. Immunofluorescence analysis revealed that the LUADs exhibited abundant expression of surfactant protein C and lacked the clara cell marker Ccsp, suggesting that these LUADs originated from alveolar type II cells. Next, we sought to survey genome-wide alterations in the pathogenesis of Gprc5a-/- LUADs. Using whole exome sequencing, we found that carcinogen-induced LUADs exhibited markedly higher somatic mutation burdens relative to spontaneous tumors. All LUADs were found to harbor somatic mutations in the Kras oncogene (p. G12D or p. Q61R). In contrast to spontaneous lesions, carcinogen-induced Gprc5a-/- LUADs exhibited mutations (variants and copy number gains) in additional drivers (Atm, Kmt2d, Nf1, Trp53, Met, Ezh2). Our study underscores genomic alterations that represent early events in the development of Kras mutant LUAD following Gprc5a loss and tobacco carcinogen exposure and that may constitute targets for prevention and early treatment of this disease.