Small cell lung cancer (SCLC) tumors comprise heterogeneous mixtures of cell states, categorized into neuroendocrine (NE) and non-neuroendocrine (non-NE) transcriptional subtypes. NE to non-NE state transitions, fueled by plasticity, likely underlie adaptability to treatment and dismal survival rates. Here, we apply an archetypal analysis to model plasticity by recasting SCLC phenotypic heterogeneity through multi-task evolutionary theory. Cell line and tumor transcriptomics data fit well in a five-dimensional convex polytope whose vertices optimize tasks reminiscent of pulmonary NE cells, the SCLC normal counterparts. These tasks, supported by knowledge and experimental data, include proliferation, slithering, metabolism, secretion, and injury repair, reflecting cancer hallmarks. SCLC subtypes, either at the population or single-cell level, can be positioned in archetypal space by bulk or single-cell transcriptomics, respectively, and characterized as task specialists or multi-task generalists by the distance from archetype vertex signatures. In the archetype space, modeling single-cell plasticity as a Markovian process along an underlying state manifold indicates that task trade-offs, in response to microenvironmental perturbations or treatment, may drive cell plasticity. Stifling phenotypic transitions and plasticity may provide new targets for much-needed translational advances in SCLC. A record of this paper’s Transparent Peer Review process is included in the supplemental information.
2021
J Thor Onc
Beyond Programmed Death-Ligand 1: B7-H6 Emerges as a Potential Immunotherapy Target in SCLC
Portia L. Thomas, Sarah M. Groves, Yun-Kai Zhang, and 16 more authors
Introduction The programmed death-ligand 1 (PD-L1) immune checkpoint inhibitors, atezolizumab and durvalumab, have received regulatory approval for the first-line treatment of patients with extensive-stage SCLC. Nevertheless, when used in combination with platinum-based chemotherapy, these PD-L1 inhibitors only improve overall survival by 2 to 3 months. This may be due to the observation that less than 20% of SCLC tumors express PD-L1 at greater than 1%. Evaluating the composition and abundance of checkpoint molecules in SCLC may identify molecules beyond PD-L1 that are amenable to therapeutic targeting. Methods We analyzed RNA-sequencing data from SCLC cell lines (n = 108) and primary tumor specimens (n = 81) for expression of 39 functionally validated inhibitory checkpoint ligands. Furthermore, we generated tissue microarrays containing SCLC cell lines and patient with SCLC specimens to confirm expression of these molecules by immunohistochemistry. We annotated patient outcomes data, including treatment response and overall survival. Results The checkpoint protein B7-H6 (NCR3LG1) exhibited increased protein expression relative to PD-L1 in cell lines and tumors (p < 0.05). Higher B7-H6 protein expression correlated with longer progression-free survival (p = 0.0368) and increased total immune infiltrates (CD45+) in patients. Furthermore, increased B7-H6 gene expression in SCLC tumors correlated with a decreased activated natural killer cell gene signature, suggesting a complex interplay between B7-H6 expression and immune signature in SCLC. Conclusions We investigated 39 inhibitory checkpoint molecules in SCLC and found that B7-H6 is highly expressed and associated with progression-free survival. In addition, 26 of 39 immune checkpoint proteins in SCLC tumors were more abundantly expressed than PD-L1, indicating an urgent need to investigate additional checkpoint targets for therapy in addition to PD-L1.
Cancer Cell
Patterns of transcription factor programs and immune pathway activation define four major subtypes of SCLC with distinct therapeutic vulnerabilities
Carl M. Gay, C. Allison Stewart, Elizabeth M. Park, and 27 more authors
Despite molecular and clinical heterogeneity, small cell lung cancer (SCLC) is treated as a single entity with predictably poor results. Using tumor expression data and non-negative matrix factorization, we identify four SCLC subtypes defined largely by differential expression of transcription factors ASCL1, NEUROD1, and POU2F3 or low expression of all three transcription factor signatures accompanied by an Inflamed gene signature (SCLC-A, N, P, and I, respectively). SCLC-I experiences the greatest benefit from the addition of immunotherapy to chemotherapy, while the other subtypes each have distinct vulnerabilities, including to inhibitors of PARP, Aurora kinases, or BCL-2. Cisplatin treatment of SCLC-A patient-derived xenografts induces intratumoral shifts toward SCLC-I, supporting subtype switching as a mechanism of acquired platinum resistance. We propose that matching baseline tumor subtype to therapy, as well as manipulating subtype switching on therapy, may enhance depth and duration of response for SCLC patients.
Genes & Dev
ASCL1 represses a SOX9+ neural crest stem-like state in small cell lung cancer
Rachelle R. Olsen, Abbie S. Ireland, David W. Kastner, and 12 more authors
ASCL1 is a neuroendocrine lineage-specific oncogenic driver of small cell lung cancer (SCLC), highly expressed in a significant fraction of tumors. However, ∼25% of human SCLC are ASCL1-low and associated with low neuroendocrine fate and high MYC expression. Using genetically engineered mouse models (GEMMs), we show that alterations in Rb1/Trp53/Myc in the mouse lung induce an ASCL1+ state of SCLC in multiple cells of origin. Genetic depletion of ASCL1 in MYC-driven SCLC dramatically inhibits tumor initiation and progression to the NEUROD1+ subtype of SCLC. Surprisingly, ASCL1 loss promotes a SOX9+ mesenchymal/neural crest stem-like state and the emergence of osteosarcoma and chondroid tumors, whose propensity is impacted by cell of origin. ASCL1 is critical for expression of key lineage-related transcription factors NKX2-1, FOXA2, and INSM1 and represses genes involved in the Hippo/Wnt/Notch developmental pathways in vivo. Importantly, ASCL1 represses a SOX9/RUNX1/RUNX2 program in vivo and SOX9 expression in human SCLC cells, suggesting a conserved function for ASCL1. Together, in a MYC-driven SCLC model, ASCL1 promotes neuroendocrine fate and represses the emergence of a SOX9+ nonendodermal stem-like fate that resembles neural crest.
2019
PLOS Comp Bio
Systems-level network modeling of Small Cell Lung Cancer subtypes identifies master regulators and destabilizers
David J. Wooten, Sarah M. Groves, Darren R. Tyson, and 6 more authors
Adopting a systems approach, we devise a general workflow to define actionable subtypes in human cancers. Applied to small cell lung cancer (SCLC), the workflow identifies four subtypes based on global gene expression patterns and ontologies. Three correspond to known subtypes (SCLC-A, SCLC-N, and SCLC-Y), while the fourth is a previously undescribed ASCL1+ neuroendocrine variant (NEv2, or SCLC-A2). Tumor deconvolution with subtype gene signatures shows that all of the subtypes are detectable in varying proportions in human and mouse tumors. To understand how multiple stable subtypes can arise within a tumor, we infer a network of transcription factors and develop BooleaBayes, a minimally-constrained Boolean rule-fitting approach. In silico perturbations of the network identify master regulators and destabilizers of its attractors. Specific to NEv2, BooleaBayes predicts ELF3 and NR0B1 as master regulators of the subtype, and TCF3 as a master destabilizer. Since the four subtypes exhibit differential drug sensitivity, with NEv2 consistently least sensitive, these findings may lead to actionable therapeutic strategies that consider SCLC intratumoral heterogeneity. Our systems-level approach should generalize to other cancer types.
2016
Phys Rev C
Consistency of electron scattering data with a small proton radius
We determine the charge radius of the proton by analyzing the published low momentum transfer electron-proton scattering data from Mainz. We note that polynomial expansions of the form factor converge for momentum transfers squared below 4mπ2, where mπ is the pion mass. Expansions with enough terms to fit the data, but few enough not to overfit, yield proton radii smaller than the CODATA or Mainz values and in accord with the muonic atom results. We also comment on analyses using a wider range of data, and overall obtain a proton radius RE=0.840(16) fm.