Deciphering intratumoral heterogeneity using integrated clonal tracking and single-cell transcriptome analyses

Cellular heterogeneity is a major cause of treatment resistance in cancer. Despite recent advances in single-cell genomic and transcriptomic sequencing, it remains difficult to relate measured molecular profiles to the cellular activities underlying cancer. Here, we present an integrated experimental system that connects single cell gene expression to heterogeneous cancer cell growth, metastasis, and treatment response. Our system integrates single cell transcriptome profiling with DNA barcode based clonal tracking in patient-derived xenograft models. We show that leukemia cells exhibiting unique gene expression respond to different chemotherapies in distinct but consistent manners across multiple mice. In addition, we uncover a form of leukemia expansion that is spatially confined to the bone marrow of single anatomical sites and driven by cells with distinct gene expression. Our integrated experimental system can interrogate the molecular and cellular basis of the intratumoral heterogeneity underlying disease progression and treatment resistance.

Intratumoral Heterogeneity and Immune Response Indicators to Predict Overall Survival in a Retrospective Study of HER2-Borderline (IHC 2+) Breast Cancer Patients

Breast cancer (BC) categorized as human epidermal growth factor receptor 2 (HER2) borderline [2+ by immunohistochemistry (IHC 2+)] presents challenges for the testing, frequently obscured by intratumoral heterogeneity (ITH). This leads to difficulties in therapy decisions. We aimed to establish prognostic models of overall survival (OS) of these patients, which take into account spatial aspects of ITH and tumor microenvironment by using hexagonal tiling analytics of digital image analysis (DIA). In particular, we assessed the prognostic value of Immunogradient indicators at the tumor–stroma interface zone (IZ) as a feature of antitumor immune response. Surgical excision samples stained for estrogen receptor (ER), progesterone receptor (PR), Ki67, HER2, and CD8 from 275 patients with HER2 IHC 2+ invasive ductal BC were used in the study. DIA outputs were subsampled by HexT for ITH quantification and tumor microenvironment extraction for Immunogradient indicators. Multiple Cox regression revealed HER2 membrane completeness (HER2 MC) (HR: 0.18, p = 0.0007), its spatial entropy (HR: 0.37, p = 0.0341), and ER contrast (HR: 0.21, p = 0.0449) as independent predictors of better OS, with worse OS predicted by pT status (HR: 6.04, p = 0.0014) in the HER2 non-amplified patients. In the HER2-amplified patients, HER2 MC contrast (HR: 0.35, p = 0.0367) and CEP17 copy number (HR: 0.19, p = 0.0035) were independent predictors of better OS along with worse OS predicted by pN status (HR: 4.75, p = 0.0018). In the non-amplified tumors, three Immunogradient indicators provided the independent prognostic value: CD8 density in the tumor aspect of the IZ and CD8 center of mass were associated with better OS (HR: 0.23, p = 0.0079 and 0.14, p = 0.0014, respectively), and CD8 density variance along the tumor edge predicted worse OS (HR: 9.45, p = 0.0002). Combining these three computational indicators of the CD8 cell spatial distribution within the tumor microenvironment augmented prognostic stratification of the patients. In the HER2-amplified group, CD8 cell density in the tumor aspect of the IZ was the only independent immune response feature to predict better OS (HR: 0.22, p = 0.0047). In conclusion, we present novel prognostic models, based on computational ITH and Immunogradient indicators of the IHC biomarkers, in HER2 IHC 2+ BC patients.

Intratumoral Heterogeneity Promotes Collective Cancer Invasion through NOTCH1 Variation

Cellular and molecular heterogeneity within tumors has long been associated with the progression of cancer to an aggressive phenotype and a poor prognosis. However, how such intratumoral heterogeneity contributes to the invasiveness of cancer is largely unknown. Here, using a tumor bioengineering approach, we investigate the interaction between molecular subtypes within bladder microtumors and the corresponding effects on their invasiveness. Our results reveal heterogeneous microtumors formed by multiple molecular subtypes possess enhanced invasiveness compared to individual cells, even when both cells are not invasive individually. To examine the molecular mechanism of intratumoral heterogeneity mediated invasiveness, live single cell biosensing, RNA interference, and CRISPR-Cas9 gene editing approaches were applied to investigate and control the composition of the microtumors. An agent-based computational model was also developed to evaluate the influence of NOTCH1 variation on DLL4 expression within a microtumor. The data indicate that intratumoral variation in NOTCH1 expression can lead to upregulation of DLL4 expression within the microtumor and enhancement of microtumor invasiveness. Overall, our results reveal a novel mechanism of heterogeneity mediated invasiveness through intratumoral variation of gene expression.

Multiscale heterogeneity in gastric adenocarcinoma evolution is an obstacle to precision medicine

Background Cancer is a somatic evolutionary disease and adenocarcinomas of the stomach and gastroesophageal junction (GC) may serve as a two-dimensional model of cancer expansion, in which tumor subclones are not evenly mixed during tumor progression but rather spatially separated and diversified. We hypothesize that precision medicine efforts are compromised when clinical decisions are based on a single-sample analysis, which ignores the mechanisms of cancer evolution and resulting intratumoral heterogeneity. Using multiregional whole-exome sequencing, we investigated the effect of somatic evolution on intratumoral heterogeneity aiming to shed light on the evolutionary biology of GC. Methods The study comprised a prospective discovery cohort of 9 and a validation cohort of 463 GCs. Multiregional whole-exome sequencing was performed using samples form 45 primary tumors and 3 lymph node metastases (range 3–10 tumor samples/patient) of the discovery cohort. Results In total, the discovery cohort harbored 16,537 non-synonymous mutations. Intratumoral heterogeneity of somatic mutations and copy number variants were present in all tumors of the discovery cohort. Of the non-synonymous mutations, 53–91% were not present in each patient’s sample; 399 genes harbored 2–4 different non-synonymous mutations in the same patient; 175 genes showed copy number variations, the majority being heterogeneous, including CD274 (PD-L1). Multi-sample tree-based analyses provided evidence for branched evolution being most complex in a microsatellite instable GC. The analysis of the mode of evolution showed a high degree of heterogeneity in deviation from neutrality within each tumor. We found evidence of parallel evolution and evolutionary trajectories: different mutations of SMAD4 aligned with different subclones and were found only in TP53 mutant GCs. Conclusions Neutral and non-neutral somatic evolution shape the mutational landscape in GC along its lateral expansions. It leads to complex spatial intratumoral heterogeneity, where lymph node metastases may stem from different areas of the primary tumor, synchronously. Our findings may have profound effects on future patient management. They illustrate the risk of mis-interpreting tumor genetics based on single-sample analysis and open new avenues for an evolutionary classification of GC, i.e., the discovery of distinct evolutionary trajectories which can be utilized for precision medicine.

Reviewing cancer’s biology: an eclectic approach

Background Cancer refers to a group of some of the worldwide most diagnosed and deadliest pathophysiological conditions that conquered researchers’ attention for decades and yet begs for more questions for a full comprehension of its complex cellular and molecular pathology. Main body The disease conditions are commonly characterized by unrestricted cell proliferation and dysfunctional replicative senescence pathways. In fact, the cell cycle operates under the rigorous control of complex signaling pathways involving cyclins and cyclin-dependent kinases assumed to be specific to each phase of the cycle. At each of these checkpoints, the cell is checked essentially for its DNA integrity. Genetic defects observed in these molecules (i.e., cyclins, cyclin-dependent kinases) are common features of cancer cells. Nevertheless, each cancer is different concerning its molecular and cellular etiology. These could range from the genetic defects mechanisms and/or the environmental conditions favoring epigenetically harbored homeostasis driving tumorigenesis alongside with the intratumoral heterogeneity with respect to the model that the tumor follows. Conclusions This review is not meant to be an exhaustive interpretation of carcinogenesis but to summarize some basic features of the molecular etiology of cancer and the intratumoral heterogeneity models that eventually bolster anticancer drug resistance for a more efficient design of drug targeting the pitfalls of the models.

A Novel Averaging Principle Provides Insights in the Impact of Intratumoral Heterogeneity on Tumor Progression

Typically stochastic differential equations (SDEs) involve an additive or multiplicative noise term. Here, we are interested in stochastic differential equations for which the white noise is nonlinearly integrated into the corresponding evolution term, typically termed as random ordinary differential equations (RODEs). The classical averaging methods fail to treat such RODEs. Therefore, we introduce a novel averaging method appropriate to be applied to a specific class of RODEs. To exemplify the importance of our method, we apply it to an important biomedical problem, in particular, we implement the method to the assessment of intratumoral heterogeneity impact on tumor dynamics. Precisely, we model gliomas according to a well-known Go or Grow (GoG) model, and tumor heterogeneity is modeled as a stochastic process. It has been shown that the corresponding deterministic GoG model exhibits an emerging Allee effect (bistability). In contrast, we analytically and computationally show that the introduction of white noise, as a model of intratumoral heterogeneity, leads to monostable tumor growth. This monostability behavior is also derived even when spatial cell diffusion is taken into account.

Multi-region exome sequencing reveals the intratumoral heterogeneity of surgically resected small cell lung cancer

Small cell lung cancer (SCLC) is a highly malignant tumor which is eventually refractory to any treatment. Intratumoral heterogeneity (ITH) may contribute to treatment failure. However, the extent of ITH in SCLC is still largely unknown. Here, we subject 120 tumor samples from 40 stage I-III SCLC patients to multi-regional whole-exome sequencing. The most common mutant genes are TP53 (88%) and RB1 (72%). We observe a medium level of mutational heterogeneity (0.30, range 0.0~0.98) and tumor mutational burden (TMB, 10.2 mutations/Mb, range 1.1~51.7). Our SCLC samples also exhibit somatic copy number variation (CNV) across all patients, with an average CNV ITH of 0.49 (range 0.02~0.99). In terms of mutation distribution, ITH, TMB, mutation clusters, and gene signatures, patients with combined SCLC behave roughly the same way as patients with pure SCLC. This condition also exists in smoking patients and patients with EGFR mutations. A higher TMB per cluster is associated with better disease-free survival while single-nucleotide variant ITH is linked to worse overall survival, and therefore these features may be used as prognostic biomarkers for SCLC. Together, these findings demonstrate the intratumoral genetic heterogeneity of surgically resected SCLC and provide insights into resistance to treatment.

Die Tumorheterogenität des uvealen Melanoms: Diagnostisch und prognostisch wichtig?

Tumor biopsies in uveal melanoma (UM) serve mainly the purpose of prognostication and assessment of individual metastatic risk, but can be used for diagnosis in selected cases. The importance of precise information is paramount for selecting adequate surveillance protocols, patient counseling, and optimization of treatment strategies. However, intratumoral heterogeneity and sample representativity are major concerns and can interfere with the correct prediction of the patient’s prognosis. We report a series of cases of UM with distinct morphologically identifiable areas, highlighting the differences in clinical behavior, as well as histopathological and genetic features.