RNA sequencing-based cell proliferation analysis across 19 cancers identifies a subset of proliferation-informative cancers with a common survival signature

The evolution of lung adenocarcinoma precursors is associated with chromosomal instability and transition from innate to adaptive immune response/evasion

Studying lung adenocarcinoma (LUAD) early carcinogenesis is challenging, primarily due to the lack of LUAD precursors specimens. We amassed multi-omics data from 213 LUAD and LUAD precursors to identify molecular features underlying LUAD precancer evolution. We observed progressively increasing mutations, chromosomal aberrations, whole genome doubling and genomic instability from precancer to invasive LUAD, indicating aggravating chromosomal instability (CIN). Telomere shortening, a crucial genomic alteration linked to CIN, emerged at precancer stage. Moreover, later-stage lesions demonstrated increasing cancer stemness and decreasing alveolar identity, suggesting epithelial de-differentiation during early LUAD carcinogenesis. The innate immune cells progressively diminished from precancer to invasive LUAD, concomitant with a gradual recruitment of adaptive immune cells (except CD8+ and gamma-delta T cells that decreased in later stages) and upregulation of numerous immune checkpoints, suggesting LUAD precancer evolution is associated with a shift from innate to adaptive immune response and immune evasion mediated by various mechanisms.

Evidence of a pan-tissue decline in stemness during human aging

Despite their biological importance, the role of stem cells in human aging remains to be elucidated. In this work, we applied a machine learning methodology to GTEx transcriptome data and assigned stemness scores to 17,382 healthy samples from 30 human tissues aged between 20 and 79 years. We found that ~60% of the studied tissues exhibit a significant negative correlation between the subject's age and stemness score. The only significant exception was the uterus, where we observed an increased stemness with age. Moreover, we observed that stemness is positively correlated with cell proliferation and negatively correlated with cellular senescence. Finally, we also observed a trend that hematopoietic stem cells derived from older individuals might have higher stemness scores. In conclusion, we assigned stemness scores to human samples and show evidence of a pan-tissue loss of stemness during human aging, which adds weight to the idea that stem cell deterioration may contribute to human aging.

A signature of enhanced proliferation associated with response and survival to anti-PD-L1 therapy in early-stage non-small cell lung cancer

In early-stage non-small cell lung cancer, the combination of neoadjuvant anti-PD-L1 and subablative stereotactic body radiation therapy (SBRT) is associated with higher rates of major pathologic response compared to anti-PD-L1 alone. Here, we identify a 140-gene set, enriched in genes characteristic of highly proliferating cells, associated with response to the dual therapy. Analysis of on-treatment transcriptome data indicate roles for T and B cells in response. The 140-gene set is associated with disease-free survival when applied to the combined trial arms. This 140-gene set identifies a subclass of tumors in all 7 of The Cancer Genome Atlas tumor types examined. Worse survival is associated with the 140-gene signature in 5 of these tumor types. Collectively, our data support that this 140-gene set, discovered in association with response to combined anti-PD-L1 and SBRT, identifies a clinically aggressive subclass of solid tumors that may be more likely to respond to immunotherapies.

Pathway level subtyping identifies a slow-cycling biological phenotype associated with poor clinical outcomes in colorectal cancer

Molecular stratification using gene-level transcriptional data has identified subtypes with distinctive genotypic and phenotypic traits, as exemplified by the consensus molecular subtypes (CMS) in colorectal cancer (CRC). Here, rather than gene-level data, we make use of gene ontology and biological activation state information for initial molecular class discovery. In doing so, we defined three pathway-derived subtypes (PDS) in CRC: PDS1 tumors, which are canonical/LGR5+ stem-rich, highly proliferative and display good prognosis; PDS2 tumors, which are regenerative/ANXA1+ stem-rich, with elevated stromal and immune tumor microenvironmental lineages; and PDS3 tumors, which represent a previously overlooked slow-cycling subset of tumors within CMS2 with reduced stem populations and increased differentiated lineages, particularly enterocytes and enteroendocrine cells, yet display the worst prognosis in locally advanced disease. These PDS3 phenotypic traits are evident across numerous bulk and single-cell datasets, and demark a series of subtle biological states that are currently under-represented in pre-clinical models and are not identified using existing subtyping classifiers.

High-resolution transcriptomics analysis of CXCL13+ EPSTI1+ CDK1+ cells with a specific focus on lung adenocarcinoma

Background

Programmed cell death ligand 1 (PD-L1) blocking therapy has transformed the treatment of lung adenocarcinoma (LUAD), which has significantly changed the landscape of immunotherapy. We aimed to explore specific cell subpopulations to understand tumor progression and identify markers of response to PD-L1 blocking therapy.

Methods

Bulk, fluorescence-activated cell sorting (FACS), and single-cell RNA (scRNA) sequencing were used to profile CXCL13, EPSTI1, and CDK1. The gene set variation analysis (GSVA) R package was utilized for score calculation, and prognostic analyses included receiver operating characteristic (ROC) curves, Cox proportional hazard models, and meta-analysis. Additionally, we analyzed tumor microenvironment (TME), genomics, compound perturbations, and clinical indicators. The high-dimensional analysis captured the intrinsic characteristics of the subpopulation. Furthermore, subpopulation differential genes were used for enrichment analysis of transcription factors and compounds.

Results

Literature and website analyses supported the essential role of CXCL13, CDK1, and EPSTI1 in immunotherapy. This led us to focus specifically on LUAD by representing a pan-cancer profile of immune-sensitive genes. Logically, the high-characteristic population may consist of samples positive for CXCL13, EPSTI1, and CDK1. The three-gene signature was a favorable indicator of immunotherapy response in the Stand Up to Cancer-Mark Foundation (SU2C-MARK) LUAD cohort but showed a poor prognosis before treatment in the Lung Cancer Explorer (LCE) database. Further mechanistic exploration revealed specific mutations associated with the three-gene signature in SU2C-MARK LUAD, such as STK11. In The Cancer Genome Atlas (TCGA)-LUAD cohort, the high-scoring group exhibited a higher tumor mutational burden (TMB) and global methylation but a lower fraction genome altered (FGA) and estimated tumor purity. Moreover, dasatinib demonstrated sensitivity in the high-scoring group. The co-localization of the CXCL13, EPSTI1, and CDK1 subpopulation was validated through spatial transcriptome and immunohistochemical databases. Assessment of the subpopulation depicted high-resolution intercellular communication. Maintenance of specific pathways, such as TNF, CD74, and CD44, contributed to immunotherapy sensitivity. Finally, the subpopulation-enriched targets and drugs were confirmed through ConnectivityMap (CMAP) analysis and multi-omics, respectively.

Conclusions

In this study, positive samples for CXCL13, EPSTI1, and CDK1 exhibited poor prognostic significance in treatment-naïve LUAD cases but demonstrated benefits from PD-L1 blockade and dasatinib therapies.

A Novel Mathematical Approach for Analysis of Integrated Cell-Patient Data Uncovers a 6-Gene Signature Linked to Endocrine Therapy Resistance

A significant number of breast cancers develop resistance to hormone therapy. This progression, while posing a major clinical challenge, is difficult to predict. Despite important contributions made by cell models and clinical studies to tackle this problem, both present limitations when taken individually. Experiments with cell models are highly reproducible but do not reflect the indubitable heterogenous landscape of breast cancer. On the other hand, clinical studies account for this complexity but introduce uncontrolled noise due to external factors. Here, we propose a new approach for biomarker discovery that is based on a combined analysis of sequencing data from controlled MCF7 cell experiments and heterogenous clinical samples that include clinical and sequencing information from The Cancer Genome Atlas. Using data from differential gene expression analysis and a Bayesian logistic regression model coupled with an original simulated annealing-type algorithm, we discovered a novel 6-gene signature for stratifying patient response to hormone therapy. The experimental observations and computational analysis built on independent cohorts indicated the superior predictive performance of this gene set over previously known signatures of similar scope. Together, these findings revealed a new gene signature to identify patients with breast cancer with an increased risk of developing resistance to endocrine therapy.

Histone methyltransferase activity affects metabolism in human cells independently of transcriptional regulation

The N-terminal tails of eukaryotic histones are frequently posttranslationally modified. The role of these modifications in transcriptional regulation is well-documented. However, the extent to which the enzymatic processes of histone posttranslational modification might affect metabolic regulation is less clear. Here, we investigated how histone methylation might affect metabolism using metabolomics, proteomics, and RNA-seq data from cancer cell lines, primary tumour samples and healthy tissue samples. In cancer, the expression of histone methyltransferases (HMTs) was inversely correlated to the activity of NNMT, an enzyme previously characterised as a methyl sink that disposes of excess methyl groups carried by the universal methyl donor S-adenosyl methionine (SAM or AdoMet). In healthy tissues, histone methylation was inversely correlated to the levels of an alternative methyl sink, PEMT. These associations affected the levels of multiple histone marks on chromatin genome-wide but had no detectable impact on transcriptional regulation. We show that HMTs with a variety of different associations to transcription are co-regulated by the Retinoblastoma (Rb) tumour suppressor in human cells. Rb-mutant cancers show increased total HMT activity and down-regulation of NNMT. Together, our results suggest that the total activity of HMTs affects SAM metabolism, independent of transcriptional regulation.

Defining the landscape of circular RNAs in neuroblastoma unveils a global suppressive function of MYCN

Circular RNAs (circRNAs) are a regulatory RNA class. While cancer-driving functions have been identified for single circRNAs, how they modulate gene expression in cancer is not well understood. We investigate circRNA expression in the pediatric malignancy, neuroblastoma, through deep whole-transcriptome sequencing in 104 primary neuroblastomas covering all risk groups. We demonstrate that MYCN amplification, which defines a subset of high-risk cases, causes globally suppressed circRNA biogenesis directly dependent on the DHX9 RNA helicase. We detect similar mechanisms in shaping circRNA expression in the pediatric cancer medulloblastoma implying a general MYCN effect. Comparisons to other cancers identify 25 circRNAs that are specifically upregulated in neuroblastoma, including circARID1A. Transcribed from the ARID1A tumor suppressor gene, circARID1A promotes cell growth and survival, mediated by direct interaction with the KHSRP RNA-binding protein. Our study highlights the importance of MYCN regulating circRNAs in cancer and identifies molecular mechanisms, which explain their contribution to neuroblastoma pathogenesis.

Improved Delineation of Colorectal Cancer Molecular Subtypes and Functional Profiles with a 62-Gene Panel

Since its establishment in 2015, the transcriptomics-based consensus molecular subtype (CMS) classification has unified our understanding of colorectal cancer. Each of the four CMS exhibited distinctive high-level molecular signatures that correlated well with prognosis and treatment response. Nonetheless, many key aspects of colorectal cancer progression and intra-subtype heterogeneity remain unresolved. This is partly because the bulk transcriptomic data used to define CMS contain substantial interference from non-tumor cells. Here, we propose a concise panel of 62 genes that not only accurately recapitulates all key characteristics of the four original CMS but also identifies three additional subpopulations with unique molecular signatures. Validation on independent cohorts confirms that the new CMS4 intra-subtypes coincide with single-cell-derived intrinsic subtypes and that the panel consists of many immune cell-type markers that can capture the status of tumor microenvironment. Furthermore, a 2D embedding of CMS structure based on the proposed gene panel provides a high-resolution view of the functional pathways and cell-type markers that underlie each CMS intra-subtype and the continuous progression from CMS2 to CMS4 subtypes. Our gene panel and 2D visualization refined the delineation of colorectal cancer subtypes and could aid further discovery of molecular mechanisms in colorectal cancer.

IMPLICATIONS

: Well-selected gene panel and representation can capture both the continuum of cancer cell states and tumor microenvironment status.

Molecular phenotyping of malignant canine mammary tumours: Detection of high-risk group and its relationship with clinicomolecular characteristics

Canine mammary gland tumours (CMTs) constitute the most common cancer in female dogs and comprise approximately 50% of all canine cancers. With the advent of high-throughput technologies such as microarray and next-generation sequencing, the molecular phenotyping (classification) of various cancers has been extensively developed. The present study used a canine RNA-sequencing dataset, namely GSE119810, to classify 113 malignant CMTs and 64 matched normal samples via an unsupervised hierarchical algorithm with a view to evaluating the association between the resulting subtypes (clusters) (n = 4) and clinical and molecular characteristics. Finally, a molecular classifier was developed, and it detected 1 high-risk molecular subtype in the training dataset (GSE119810) and 2 independent validation datasets (GSE20718 and GSE22516). Our results revealed four molecular subtypes (C2-C5) in malignant CMTs. Furthermore, the normal samples constituted a distinct group in the clustering analysis. Marked significant associations were observed between the molecular subtypes (especially C5) and clinical/molecular features, including positive lymphatic invasion, high tumour grades, histopathology diagnoses, short survival and high TP53 mutation rates (ps <.05). The high-risk subtype (C5) was further characterized through the development of a cell cycle-based gene signature, which comprised 37 proliferation-related genes according to the support vector machine algorithm. This signature identified the high-risk group in both training and validation datasets (ps <.001). In the validation analysis, our potential classifier robustly predicted patients with positive lymphatic invasion, metastases and short survival.

Deciphering radiological stable disease to immune checkpoint inhibitors

BACKGROUND

'Stable disease (SD)' as per RECIST is a common but ambiguous outcome in patients receiving immune checkpoint inhibitors (ICIs). This study aimed to characterize SD and identify the subset of patients with SD who are benefiting from treatment. Understanding SD would facilitate drug development and improve precision in correlative research.

PATIENTS AND METHODS

A systematic review was carried out to characterize SD in ICI trials. SD and objective response were compared to proliferation index using The Cancer Genome Atlas gene expression data. To identify a subgroup of SD with outcomes mirroring responders, we examined a discovery cohort of non-small-cell lung cancer (NSCLC). Serial cutpoints of two variables, % best overall response and progression-free survival (PFS), were tested to define a subgroup of patients with SD with similar survival as responders. Results were then tested in external validation cohorts.

RESULTS

Among trials of ICIs (59 studies, 14 280 patients), SD ranged from 16% to 42% in different tumor types and was associated with disease-specific proliferation index (ρ = -0.75, P = 0.03), a proxy of tumor kinetics, rather than relative response to ICIs. In a discovery cohort of NSCLC [1220 patients, 313 (26%) with SD to ICIs], PFS ranged widely in SD (0.2-49 months, median 4.9 months). The subset with PFS >6 months and no tumor growth mirrored partial response (PR) minor (overall survival hazard ratio 1.0) and was proposed as the definition of SD responder. This definition was confirmed in two validation cohorts from trials of NSCLC treated with durvalumab and found to apply in tumor types treated with immunotherapy in which depth and duration of benefit were correlated.

CONCLUSIONS

RECIST-defined SD to immunotherapy is common, heterogeneous, and may largely reflect tumor growth rate rather than ICI response. In patients with NSCLC and SD to ICIs, PFS >6 months and no tumor growth may be considered 'SD responders'. This definition may improve the efficiency of and insight derivable from clinical and translational research.

ENDORSE: a prognostic model for endocrine therapy in estrogen-receptor-positive breast cancers

Advanced and metastatic estrogen receptor-positive (ER+ ) breast cancers are often endocrine resistant. However, endocrine therapy remains the primary treatment for all advanced ER+ breast cancers. Treatment options that may benefit resistant cancers, such as add-on drugs that target resistance pathways or switching to chemotherapy, are only available after progression on endocrine therapy. Here we developed an endocrine therapy prognostic model for early and advanced ER+ breast cancers. The endocrine resistance (ENDORSE) model is composed of two components, each based on the empirical cumulative distribution function of ranked expression of gene signatures. These signatures include a feature set associated with long-term survival outcomes on endocrine therapy selected using lasso-regularized Cox regression and a pathway-based curated set of genes expressed in response to estrogen. We extensively validated ENDORSE in multiple ER+ clinical trial datasets and demonstrated superior and consistent performance of the model over clinical covariates, proliferation markers, and multiple published signatures. Finally, genomic and pathway analyses in patient data revealed possible mechanisms that may help develop rational stratification strategies for endocrine-resistant ER+ breast cancer patients.

A pan-cancer landscape of telomeric content shows that RAD21 and HGF alterations are associated with longer telomeres

Background

Cancer cells can proliferate indefinitely through telomere maintenance mechanisms. These mechanisms include telomerase-dependent elongation, mediated by TERT activation, and alternative lengthening of telomeres (ALT), linked to loss of ATRX or DAXX.

Methods

We analyzed the telomeric content of 89,959 tumor samples within the Foundation Medicine dataset and investigated the genomic determinants of high telomeric content, linking them to clinical outcomes, when available.

Results

Telomeric content varied widely by disease type with leiomyosarcoma having the highest and Merkel cell carcinoma having the lowest telomeric content. In agreement with previous studies, telomeric content was significantly higher in samples with alterations in TERC, ATRX, and DAXX. We further identified that amplifications in two genes, RAD21 and HGF, were enriched in samples with high telomeric content, which was confirmed using the PCAWG/ICGC dataset. We identified the minimal amplified region associated with high telomeric content for RAD21 (8q23.1–8q24.12), which excludes MYC, and for HGF (7q21.11). Our results demonstrated that RAD21 and HGF exerted an additive telomere lengthening effect on samples with existing alterations in canonical genes previously associated with telomere elongation. Furthermore, patients with breast cancer who harbor RAD21 alterations had poor median overall survival and trended towards higher levels of Ki-67 staining.

Conclusions

This study highlights the importance of the role played by RAD21 (8q23.1–8q24.12) and HGF (7q21.11) in the lengthening of telomeres, supporting unlimited replication in tumors. These findings open avenues for work aimed at targeting this crucial pathway in tumorigenesis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-022-01029-7.

Malignancy and NF-κB signalling strengthen coordination between expression of mitochondrial and nuclear-encoded oxidative phosphorylation genes

BACKGROUND

Mitochondria are ancient endosymbiotic organelles crucial to eukaryotic growth and metabolism. The mammalian mitochondrial genome encodes for 13 mitochondrial proteins, and the remaining mitochondrial proteins are encoded by the nuclear genome. Little is known about how coordination between the expression of the two sets of genes is achieved.

RESULTS

Correlation analysis of RNA-seq expression data from large publicly available datasets is a common method to leverage genetic diversity to infer gene co-expression modules. Here we use this method to investigate nuclear-mitochondrial gene expression coordination. We identify a pitfall in correlation analysis that results from the large variation in the proportion of transcripts from the mitochondrial genome in RNA-seq data. Commonly used normalisation techniques based on total read counts, such as FPKM or TPM, produce artefactual negative correlations between mitochondrial- and nuclear-encoded transcripts. This also results in artefactual correlations between pairs of nuclear-encoded genes, with important consequences for inferring co-expression modules beyond mitochondria. We show that these effects can be overcome by normalizing using the median-ratio normalisation (MRN) or trimmed mean of M values (TMM) methods. Using these normalisations, we find only weak and inconsistent correlations between mitochondrial and nuclear-encoded mitochondrial genes in the majority of healthy human tissues from the GTEx database.

CONCLUSIONS

We show that a subset of healthy tissues with high expression of NF-κB show significant coordination, suggesting a role for NF-κB in ensuring balanced expression between mitochondrial and nuclear genes. Contrastingly, most cancer types show robust coordination of nuclear and mitochondrial OXPHOS gene expression, identifying this as a feature of gene regulation in cancer.

Datopotamab Deruxtecan, a Novel TROP2-directed Antibody-drug Conjugate, Demonstrates Potent Antitumor Activity by Efficient Drug Delivery to Tumor Cells

Trophoblast cell surface antigen 2 (TROP2) is highly expressed on various epithelial tumors and correlates with poor prognosis. We developed the novel TROP2-directed antibody-drug conjugate (ADC), datopotamab deruxtecan (Dato-DXd, DS-1062a), with a potent DNA topoisomerase I inhibitor (DXd), and evaluated its antitumor activity and safety profiles in preclinical models.The pharmacologic activity and mechanism of action of Dato-DXd were investigated in several human cancer cell lines and xenograft mouse models including patient-derived xenograft (PDX) models. Safety profiles were also assessed in rats and cynomolgus monkeys.Dato-DXd bound specifically to TROP2 and was internalized into tumor cells followed by intracellular trafficking to lysosome and DXd release, which induced DNA damage and apoptosis in TROP2-expressing tumor cells in vitro. Dato-DXd exhibited in vivo antitumor activity with DNA damage induced by the accumulated DXd in TROP2-expressing xenograft tumors, but neither isotype control IgG-ADC nor anti-TROP2 antibody had this effect. Dato-DXd also showed potent antitumor activity with tumor regression in several TROP2-expressing xenograft tumors including NSCLC PDX models. Safety profiles of Dato-DXd in rats and cynomolgus monkeys were acceptable.Dato-DXd demonstrated potent antitumor activity against TROP2-expressing tumors by efficient payload delivery into tumors and acceptable safety profiles in preclinical models. These results suggest Dato-DXd could be a valuable treatment option for patients with TROP2-expressing tumors in the clinical setting.

Pan-Cancer Molecular Patterns and Biological Implications Associated with a Tumor-Specific Molecular Signature

Studying tissue-independent components of cancer and defining pan-cancer subtypes could be addressed using tissue-specific molecular signatures if classification errors are controlled. Since PAM50 is a well-known, United States Food and Drug Administration (FDA)-approved and commercially available breast cancer signature, we applied it with uncertainty assessment to classify tumor samples from over 33 cancer types, discarded unassigned samples, and studied the emerging tumor-agnostic molecular patterns. The percentage of unassigned samples ranged between 55.5% and 86.9% in non-breast tissues, and gene set analysis suggested that the remaining samples could be grouped into two classes (named C1 and C2) regardless of the tissue. The C2 class was more dedifferentiated, more proliferative, with higher centrosome amplification, and potentially more TP53 and RB1 mutations. We identified 28 gene sets and 95 genes mainly associated with cell-cycle progression, cell-cycle checkpoints, and DNA damage that were consistently exacerbated in the C2 class. In some cancer types, the C1/C2 classification was associated with survival and drug sensitivity, and modulated the prognostic meaning of the immune infiltrate. Our results suggest that PAM50 could be repurposed for a pan-cancer context when paired with uncertainty assessment, resulting in two classes with molecular, biological, and clinical implications.

MYC Activity Inference Captures Diverse Mechanisms of Aberrant MYC Pathway Activation in Human Cancers

c-MYC (MYC) is deregulated in more than 50% of all cancers. While MYC amplification is the most common MYC-deregulating event, many other alterations can increase MYC activity. We thus systematically investigated MYC pathway activity across different tumor types. Using a logistic regression framework, we established tumor type-specific, transcriptomic-based MYC activity scores that can accurately capture MYC activity. We show that MYC activity scores reflect a variety of MYC-regulating mechanisms, including MYCL and/or MYCN amplification, MYC promoter methylation, MYC mRNA expression, lncRNA PVT1 expression, MYC mutations, and viral integrations near the MYC locus. Our MYC activity score incorporates all of these mechanisms, resulting in better prognostic predictions compared with MYC amplification status, MYC promoter methylation, and MYC mRNA expression in several cancer types. In addition, we show that tumor proliferation and immune evasion are likely contributors to this reduction in survival. Finally, we developed a MYC activity signature for liquid tumors in which MYC translocation is commonly observed, suggesting that our approach can be applied to different types of genomic alterations. In conclusion, we developed a MYC activity score that captures MYC pathway activity and is clinically relevant. IMPLICATIONS: By using cancer type-specific MYC activity profiles, we were able to assess MYC activity across many more tumor types than previously investigated. The range of different MYC-related alterations captured by our MYC activity score can be used to facilitate the application of future MYC inhibitors and aid physicians to preselect patients for targeted therapy.

Predictive biomarkers for sacituzumab govitecan efficacy in Trop-2-expressing triple-negative breast cancer

Sacituzumab govitecan (SG) is an antibody-drug conjugate composed of a humanized anti-Trop-2 IgG antibody conjugated via a hydrolysable linker to SN-38, the topoisomerase I-inhibitory active component of irinotecan. We investigated whether Trop-2-expression and homologous recombination repair (HRR) of SN-38-mediated double-strand DNA (dsDNA) breaks play a role in the sensitivity of triple-negative breast cancer (TNBC) to SG. Activation of HRR pathways, as evidenced by Rad51 expression, was assessed in SG-sensitive cell lines with low and moderate Trop-2-expression (SK-MES-1 squamous cell lung carcinoma and HCC1806 TNBC, respectively), compared to a low Trop-2-expressing, less SG-sensitive TNBC cell line (MDA-MB-231). Further, two Trop-2-transfectants of MDA-MB-231, C13 and C39 (4- and 25-fold higher Trop-2, respectively), were treated in mice with SG to determine whether increasing Trop-2 expression improves SG efficacy. SG mediated >2-fold increase in Rad51 in MDA-MB-231 but had no effect in SK-MES-1 or HCC1806, resulting in lower levels of dsDNA breaks in MDA-MB-231. SG and saline produced similar effects in parental MDA-MB-231 tumor-bearing mice (median survival time (MST) = 21d and 19.5d, respectively). However, in mice bearing higher Trop-2-expressing C13 and C39 tumors after Trop-2 transfection, SG provided a significant survival benefit, even compared to irinotecan (MST = 97d vs. 35d for C13, and 81d vs. 28d for C39, respectively; P < 0.0007). These results suggest that SG could provide better clinical benefit than irinotecan in patients with HRR-proficient tumors expressing high levels of Trop-2, as well as to patients with HRR-deficient tumors expressing low/moderate levels of Trop-2.

Non-Genetic Intra-Tumor Heterogeneity Is a Major Predictor of Phenotypic Heterogeneity and Ongoing Evolutionary Dynamics in Lung Tumors

Impacts of genetic and non-genetic intra-tumor heterogeneity (ITH) on tumor phenotypes and evolvability remain debated. We analyze ITH in lung squamous cell carcinoma at the levels of genome, transcriptome, and tumor-immune interactions and histopathological characteristics by multi-region bulk and single-cell sequencing. Genomic heterogeneity alone is a weak indicator of intra-tumor non-genetic heterogeneity at immune and transcriptomic levels that impact multiple cancer-related pathways, including those related to proliferation and inflammation, which in turn contribute to intra-tumor regional differences in histopathology and subtype classification. Tumor subclones have substantial differences in proliferation score, suggestive of non-neutral clonal dynamics. Proliferation and other cancer-related pathways also show intra-tumor regional differences, sometimes even within the same subclones. Neo-epitope burden negatively correlates with immune infiltration, indicating immune-mediated purifying selection on somatic mutations. Taken together, our observations suggest that non-genetic heterogeneity is a major determinant of heterogeneity in histopathological characteristics and impacts evolutionary dynamics in lung cancer.

A kinase-deficient NTRK2 splice variant predominates in glioma and amplifies several oncogenic signaling pathways

Independent scientific achievements have led to the discovery of aberrant splicing patterns in oncogenesis, while more recent advances have uncovered novel gene fusions involving neurotrophic tyrosine receptor kinases (NTRKs) in gliomas. The exploration of NTRK splice variants in normal and neoplastic brain provides an intersection of these two rapidly evolving fields. Tropomyosin receptor kinase B (TrkB), encoded NTRK2, is known for critical roles in neuronal survival, differentiation, molecular properties associated with memory, and exhibits intricate splicing patterns and post-translational modifications. Here, we show a role for a truncated NTRK2 splice variant, TrkB.T1, in human glioma. TrkB.T1 enhances PDGF-driven gliomas in vivo, augments PDGF-induced Akt and STAT3 signaling in vitro, while next generation sequencing broadly implicates TrkB.T1 in the PI3K signaling cascades in a ligand-independent fashion. These TrkB.T1 findings highlight the importance of expanding upon whole gene and gene fusion analyses to include splice variants in basic and translational neuro-oncology research.

Molecular natural history of breast cancer: Leveraging transcriptomics to predict breast cancer progression and aggressiveness

BACKGROUND

Characterizing breast cancer progression and aggressiveness relies on categorical descriptions of tumor stage and grade. Interpreting these categorical descriptions is challenging because stage convolutes the size and spread of the tumor and no consensus exists to define high/low grade tumors.

METHODS

We address this challenge of heterogeneity in patient-specific cancer samples by adapting and applying several tools originally created for understanding heterogeneity and phenotype development in single cells (specifically, single-cell topological data analysis and Wanderlust) to create a continuous metric describing breast cancer progression using bulk RNA-seq samples from individual patient tumors. We also created a linear regression-based method to predict tumor aggressiveness in vivo from bulk RNA-seq data.

RESULTS

We found that breast cancer proceeds along three convergent phenotype trajectories: luminal, HER2-enriched, and basal-like. Furthermore, 31 296 genes (for luminal cancers), 17 827 genes (for HER2-enriched), and 18 505 genes (for basal-like) are dynamically differentially expressed during breast cancer progression. Across progression trajectories, our results show that expression of genes related to ADP-ribosylation decreased as tumors progressed (while PARP1 and PARP2 increased or remained stable), suggesting the potential for a differential response to PARP inhibitors based on cancer progression. Additionally, we developed a 132-gene expression regression equation to predict mitotic index and a 23-gene expression regression equation to predict growth rate from a single breast cancer biopsy.

CONCLUSION

Our results suggest that breast cancer dynamically changes during disease progression, and growth rate of the cancer cells is associated with distinct transcriptional profiles.

Cardelino: computational integration of somatic clonal substructure and single-cell transcriptomes

Bulk and single-cell DNA sequencing has enabled reconstructing clonal substructures of somatic tissues from frequency and cooccurrence patterns of somatic variants. However, approaches to characterize phenotypic variations between clones are not established. Here we present cardelino (https://github.com/single-cell-genetics/cardelino), a computational method for inferring the clonal tree configuration and the clone of origin of individual cells assayed using single-cell RNA-seq (scRNA-seq). Cardelino flexibly integrates information from imperfect clonal trees inferred based on bulk exome-seq data, and sparse variant alleles expressed in scRNA-seq data. We apply cardelino to a published cancer dataset and to newly generated matched scRNA-seq and exome-seq data from 32 human dermal fibroblast lines, identifying hundreds of differentially expressed genes between cells from different somatic clones. These genes are frequently enriched for cell cycle and proliferation pathways, indicating a role for cell division genes in somatic evolution in healthy skin.

A phenotypic and genomics approach in a multi-ethnic cohort to subtype systemic lupus erythematosus

Systemic lupus erythematous (SLE) is a heterogeneous autoimmune disease in which outcomes vary among different racial groups. Here, we aim to identify SLE subgroups within a multiethnic cohort using an unsupervised clustering approach based on the American College of Rheumatology (ACR) classification criteria. We identify three patient clusters that vary according to disease severity. Methylation association analysis identifies a set of 256 differentially methylated CpGs across clusters, including 101 CpGs in genes in the Type I Interferon pathway, and we validate these associations in an external cohort. A cis-methylation quantitative trait loci analysis identifies 744 significant CpG-SNP pairs. The methylation signature is enriched for ethnic-associated CpGs suggesting that genetic and non-genetic factors may drive outcomes and ethnic-associated methylation differences. Our computational approach highlights molecular differences associated with clusters rather than single outcome measures. This work demonstrates the utility of applying integrative methods to address clinical heterogeneity in multifactorial multi-ethnic disease settings.

Proliferative potential and resistance to immune checkpoint blockade in lung cancer patients

Background

Resistance to immune checkpoint inhibitors (ICIs) has been linked to local immunosuppression independent of major ICI targets (e.g., PD-1). Clinical experience with response prediction based on PD-L1 expression suggests that other factors influence sensitivity to ICIs in non-small cell lung cancer (NSCLC) patients.

Methods

Tumor specimens from 120 NSCLC patients from 10 institutions were evaluated for PD-L1 expression by immunohistochemistry, and global proliferative profile by targeted RNA-seq.

Results

Cell proliferation, derived from the mean expression of 10 proliferation-associated genes (namely BUB1, CCNB2, CDK1, CDKN3, FOXM1, KIAA0101, MAD2L1, MELK, MKI67, and TOP2A), was identified as a marker of response to ICIs in NSCLC. Poorly, moderately, and highly proliferative tumors were somewhat equally represented in NSCLC, with tumors with the highest PD-L1 expression being more frequently moderately proliferative as compared to lesser levels of PD-L1 expression. Proliferation status had an impact on survival in patients with both PD-L1 positive and negative tumors. There was a significant survival advantage for moderately proliferative tumors compared to their combined highly/poorly counterparts (p = 0.021). Moderately proliferative PD-L1 positive tumors had a median survival of 14.6 months that was almost twice that of PD-L1 negative highly/poorly proliferative at 7.6 months (p = 0.028). Median survival in moderately proliferative PD-L1 negative tumors at 12.6 months was comparable to that of highly/poorly proliferative PD-L1 positive tumors at 11.5 months, but in both instances less than that of moderately proliferative PD-L1 positive tumors. Similar to survival, proliferation status has impact on disease control (DC) in patients with both PD-L1 positive and negative tumors. Patients with moderately versus those with poorly or highly proliferative tumors have a superior DC rate when combined with any classification schema used to score PD-L1 as a positive result (i.e., TPS ≥ 50% or ≥ 1%), and best displayed by a DC rate for moderately proliferative tumors of no less than 40% for any classification of PD-L1 as a negative result. While there is an over representation of moderately proliferative tumors as PD-L1 expression increases this does not account for the improved survival or higher disease control rates seen in PD-L1 negative tumors.

Conclusions

Cell proliferation is potentially a new biomarker of response to ICIs in NSCLC and is applicable to PD-L1 negative tumors.

Electronic supplementary material

The online version of this article (10.1186/s40425-019-0506-3) contains supplementary material, which is available to authorized users.

Genome doubling shapes the evolution and prognosis of advanced cancers

Ploidy abnormalities are a hallmark of cancer, but their impact on the evolution and outcomes of cancers is unknown. Here, we identified whole-genome doubling (WGD) in the tumors of nearly 30% of 9,692 prospectively sequenced advanced cancer patients. WGD varied by tumor lineage and molecular subtype, and arose early in carcinogenesis after an antecedent transforming driver mutation. While associated with TP53 mutations, 46% of all WGD arose in TP53-wild-type tumors and in such cases was associated with an E2F-mediated G1 arrest defect, although neither aberration was obligate in WGD tumors. The variability of WGD across cancer types can be explained in part by cancer cell proliferation rates. WGD predicted for increased morbidity across cancer types, including KRAS-mutant colorectal cancers and estrogen receptor-positive breast cancers, independently of established clinical prognostic factors. We conclude that WGD is highly common in cancer and is a macro-evolutionary event associated with poor prognosis across cancer types.

RELN signaling modulates glioblastoma growth and substrate-dependent migration

Glioblastoma (GBM) represents the most common and most malignant type of primary brain tumor and significantly contributes to cancer morbidity and mortality. Invasion into the healthy brain parenchyma is a major feature of glioblastoma aggressiveness. Reelin (RELN) is a large secreted extracellular matrix glycoprotein that regulates neuronal migration and positioning in the developing brain and sustains functionality in the adult brain. We here show that both RELN and its main downstream effector DAB1 are silenced in glioblastoma as compared to non-neoplastic tissue and mRNA expression is inversely correlated with malignancy grade. Furthermore, RELN expression is positively correlated with patient survival in two large, independent clinically annotated datasets. RELN silencing occurs via promoter hypermethylation as shown by both database mining and bisulfite sequencing of the RELN promoter. Consequently, treatment with 5'-Azacytidine and trichostatin A induced RELN expression in vitro. On the functional level, we found RELN to regulate glioblastoma cell migration both in a DAB1 (tyrosine phosphorylation)-dependent and -independent fashion, depending on the substrate provided. Moreover, stimulation of RELN signaling strongly reduced proliferation in glioblastoma cells. This phenotype depends on DAB1 stimulation by RELN, as a mutant that lacks all RELN induced tyrosine phosphorylation sites (DAB1-5F) failed to induce a growth arrest. Proteomic analyzes revealed that these effects are mediated by a reduction in E2F targets and dephosphorylation of ERK1/2. Taken together, our data establish a relevance of RELN signaling in glioblastoma pathology and thereby might unearth novel, yet unrecognized treatment options.