Single-cell Analysis Reveals Inter- and Intratumour Heterogeneity in Metastatic Breast Cancer.

Metastasis is the leading cause of cancer-related deaths of breast cancer patients. Some cancer cells in a tumour go through successive steps, referred to as the metastatic cascade, and give rise to metastases at a distant site. We know that the plasticity and heterogeneity of cancer cells play critical roles in metastasis but the precise underlying molecular mechanisms remain elusive. Here we aimed to identify molecular mechanisms of metastasis during colonization, one of the most important yet poorly understood steps of the cascade. We performed single-cell RNA-Seq (scRNA-Seq) on tumours and matched lung macrometastases of patient-derived xenografts of breast cancer. After correcting for confounding factors such as the cell cycle and the percentage of detected genes (PDG), we identified cells in three states in both tumours and metastases. Gene-set enrichment analysis revealed biological processes specific to proliferation and invasion in two states. Our findings suggest that these states are a balance between epithelial-to-mesenchymal (EMT) and mesenchymal-to-epithelial transitions (MET) traits that results in so-called partial EMT phenotypes. Analysis of the top differentially expressed genes (DEGs) between these cell states revealed a common set of partial EMT transcription factors (TFs) controlling gene expression, including ZNF750, OVOL2, TP63, TFAP2C and HEY2. Our data suggest that the TFs related to EMT delineate different cell states in tumours and metastases. The results highlight the marked interpatient heterogeneity of breast cancer but identify common features of single cells from five models of metastatic breast cancer.

Detection of isoforms and genomic alterations by high-throughput full-length single-cell RNA sequencing in ovarian cancer.

Understanding the complex background of cancer requires genotype-phenotype information in single-cell resolution. Here, we perform long-read single-cell RNA sequencing (scRNA-seq) on clinical samples from three ovarian cancer patients presenting with omental metastasis and increase the PacBio sequencing depth to 12,000 reads per cell. Our approach captures 152,000 isoforms, of which over 52,000 were not previously reported. Isoform-level analysis accounting for non-coding isoforms reveals 20% overestimation of protein-coding gene expression on average. We also detect cell type-specific isoform and poly-adenylation site usage in tumor and mesothelial cells, and find that mesothelial cells transition into cancer-associated fibroblasts in the metastasis, partly through the TGF-ß/miR-29/Collagen axis. Furthermore, we identify gene fusions, including an experimentally validated IGF2BP2::TESPA1 fusion, which is misclassified as high TESPA1 expression in matched short-read data, and call mutations confirmed by targeted NGS cancer gene panel results. With these findings, we envision long-read scRNA-seq to become increasingly relevant in oncology and personalized medicine.

Dynamic thresholding and tissue dissociation optimization for CITE-seq identifies differential surface protein abundance in metastatic melanoma.

Multi-omics profiling by CITE-seq bridges the RNA-protein gap in single-cell analysis but has been largely applied to liquid biopsies. Applying CITE-seq to clinically relevant solid biopsies to characterize healthy tissue and the tumor microenvironment is an essential next step in single-cell translational studies. In this study, gating of cell populations based on their transcriptome signatures for use in cell type-specific ridge plots allowed identification of positive antibody signals and setting of manual thresholds. Next, we compare five skin dissociation protocols by taking into account dissociation efficiency, captured cell type heterogeneity and recovered surface proteome. To assess the effect of enzymatic digestion on transcriptome and epitope expression in immune cell populations, we analyze peripheral blood mononuclear cells (PBMCs) with and without dissociation. To further assess the RNA-protein gap, RNA-protein we perform codetection and correlation analyses on thresholded protein values. Finally, in a proof-of-concept study, using protein abundance analysis on selected surface markers in a cohort of healthy skin, primary, and metastatic melanoma we identify CD56 surface marker expression on metastatic melanoma cells, which was further confirmed by multiplex immunohistochemistry. This work provides practical guidelines for processing and analysis of clinically relevant solid tissue biopsies for biomarker discovery.

Genome assemblies of the simultaneously hermaphroditic flatworms Macrostomum cliftonense and Macrostomum hystrix.

The free-living, simultaneously hermaphroditic flatworms of the genus Macrostomum are increasingly used as model systems in various contexts. In particular, Macrostomum lignano, the only species of this group with a published genome assembly, has emerged as a model for the study of regeneration, reproduction, and stem-cell function. However, challenges have emerged due to M. lignano being a hidden polyploid, having recently undergone whole-genome duplication and chromosome fusion events. This complex genome architecture presents a significant roadblock to the application of many modern genetic tools. Hence, additional genomic resources for this genus are needed. Here, we present such resources for Macrostomum cliftonense and Macrostomum hystrix, which represent the contrasting mating behaviors of reciprocal copulation and hypodermic insemination found in the genus. We use a combination of PacBio long-read sequencing and Illumina shot-gun sequencing, along with several RNA-Seq data sets, to assemble and annotate highly contiguous genomes for both species. The assemblies span ∼227 and ∼220 Mb and are represented by 399 and 42 contigs for M. cliftonense and M. hystrix, respectively. Furthermore, high BUSCO completeness (∼84-85%), low BUSCO duplication rates (8.3-6.2%), and low k-mer multiplicity indicate that these assemblies do not suffer from the same assembly ambiguities of the M. lignano genome assembly, which can be attributed to the complex karyology of this species. We also show that these resources, in combination with the prior resources from M. lignano, offer an excellent foundation for comparative genomic research in this group of organisms.

Live slow-frozen human tumor tissues viable for 2D, 3D, ex vivo cultures and single-cell RNAseq.

Biobanking of surplus human healthy and disease-derived tissues is essential for diagnostics and translational research. An enormous amount of formalin-fixed and paraffin-embedded (FFPE), Tissue-Tek OCT embedded or snap-frozen tissues are preserved in many biobanks worldwide and have been the basis of translational studies. However, their usage is limited to assays that do not require viable cells. The access to intact and viable human material is a prerequisite for translational validation of basic research, for novel therapeutic target discovery, and functional testing. Here we show that surplus tissues from multiple solid human cancers directly slow-frozen after resection can subsequently be used for different types of methods including the establishment of 2D, 3D, and ex vivo cultures as well as single-cell RNA sequencing with similar results when compared to freshly analyzed material.

Tracking and characterization of partial and full epithelial-mesenchymal transition cells in a mouse model of metastatic breast cancer.

The various stages of epithelial-mesenchymal transition (EMT) generate phenotypically heterogeneous populations of cells. Here, we detail a dual recombinase lineage tracing system using a transgenic mouse model of metastatic breast cancer to trace and characterize breast cancer cells at different EMT stages. We describe analytical steps to label cancer cells at an early partial or a late full EMT state, followed by tracking their behavior in tumor slice cultures. We then characterize their transcriptome by five-cell RNA sequencing. For complete details on the use and execution of this protocol, please refer to Luond et al. (2021).

Statistical tests for intra-tumour clonal co-occurrence and exclusivity.

Tumour progression is an evolutionary process in which different clones evolve over time, leading to intra-tumour heterogeneity. Interactions between clones can affect tumour evolution and hence disease progression and treatment outcome. Intra-tumoural pairs of mutations that are overrepresented in a co-occurring or clonally exclusive fashion over a cohort of patient samples may be suggestive of a synergistic effect between the different clones carrying these mutations. We therefore developed a novel statistical testing framework, called GeneAccord, to identify such gene pairs that are altered in distinct subclones of the same tumour. We analysed our framework for calibration and power. By comparing its performance to baseline methods, we demonstrate that to control type I errors, it is essential to account for the evolutionary dependencies among clones. In applying GeneAccord to the single-cell sequencing of a cohort of 123 acute myeloid leukaemia patients, we find 1 clonally co-occurring and 8 clonally exclusive gene pairs. The clonally exclusive pairs mostly involve genes of the key signalling pathways.

Distinct contributions of partial and full EMT to breast cancer malignancy.

Epithelial-mesenchymal transition (EMT) is a transient, reversible process of cell de-differentiation where cancer cells transit between various stages of an EMT continuum, including epithelial, partial EMT, and mesenchymal cell states. We have employed Tamoxifen-inducible dual recombinase lineage tracing systems combined with live imaging and 5-cell RNA sequencing to track cancer cells undergoing partial or full EMT in the MMTV-PyMT mouse model of metastatic breast cancer. In primary tumors, cancer cells infrequently undergo EMT and mostly transition between epithelial and partial EMT states but rarely reach full EMT. Cells undergoing partial EMT contribute to lung metastasis and chemoresistance, whereas full EMT cells mostly retain a mesenchymal phenotype and fail to colonize the lungs. However, full EMT cancer cells are enriched in recurrent tumors upon chemotherapy. Hence, cancer cells in various stages of the EMT continuum differentially contribute to hallmarks of breast cancer malignancy, such as tumor invasion, metastasis, and chemoresistance.

Collection and preprocessing of fine needle aspirate patient samples for single cell profiling and data analysis.

High cell viability and recovered cell concentration are typical quality control requirements for single-cell processing and quality data. This protocol describes procedures for sampling, live-cell biobanking, preprocessing for single-cell RNA sequencing, and analysis of fine-needle aspiration (FNA) samples of the skin. The minimally invasive nature of FNA collection is more accepted by patients and allows for frequent longitudinal sampling, resulting in high-quality single-cell sequencing data that capture cellular heterogeneity in clinical samples.

Spatial Distribution of Private Gene Mutations in Clear Cell Renal Cell Carcinoma.

Intra-tumour heterogeneity is the molecular hallmark of renal cancer, and the molecular tumour composition determines the treatment outcome of renal cancer patients. In renal cancer tumourigenesis, in general, different tumour clones evolve over time. We analysed intra-tumour heterogeneity and subclonal mutation patterns in 178 tumour samples obtained from 89 clear cell renal cell carcinoma patients. In an initial discovery phase, whole-exome and transcriptome sequencing data from paired tumour biopsies from 16 ccRCC patients were used to design a gene panel for follow-up analysis. In this second phase, 826 selected genes were targeted at deep coverage in an extended cohort of 89 patients for a detailed analysis of tumour heterogeneity. On average, we found 22 mutations per patient. Pairwise comparison of the two biopsies from the same tumour revealed that on average, 62% of the mutations in a patient were detected in one of the two samples. In addition to commonly mutated genes (VHL, PBRM1, SETD2 and BAP1), frequent subclonal mutations with low variant allele frequency (<10%) were observed in TP53 and in mucin coding genes MUC6, MUC16, and MUC3A. Of the 89 ccRCC tumours, 87 (~98%) harboured private mutations, occurring in only one of the paired tumour samples. Clonally exclusive pathway pairs were identified using the WES data set from 16 ccRCC patients. Our findings imply that shared and private mutations significantly contribute to the complexity of differential gene expression and pathway interaction and might explain the clonal evolution of different molecular renal cancer subgroups. Multi-regional sequencing is central for the identification of subclones within ccRCC.

JAK2-V617F and interferon-α induce megakaryocyte-biased stem cells characterized by decreased long-term functionality.

We studied a subset of hematopoietic stem cells (HSCs) that are defined by elevated expression of CD41 (CD41hi) and showed bias for differentiation toward megakaryocytes (Mks). Mouse models of myeloproliferative neoplasms (MPNs) expressing JAK2-V617F (VF) displayed increased frequencies and percentages of the CD41hi vs CD41lo HSCs compared with wild-type controls. An increase in CD41hi HSCs that correlated with JAK2-V617F mutant allele burden was also found in bone marrow from patients with MPN. CD41hi HSCs produced a higher number of Mk-colonies of HSCs in single-cell cultures in vitro, but showed reduced long-term reconstitution potential compared with CD41lo HSCs in competitive transplantations in vivo. RNA expression profiling showed an upregulated cell cycle, Myc, and oxidative phosphorylation gene signatures in CD41hi HSCs, whereas CD41lo HSCs showed higher gene expression of interferon and the JAK/STAT and TNFα/NFκB signaling pathways. Higher cell cycle activity and elevated levels of reactive oxygen species were confirmed in CD41hi HSCs by flow cytometry. Expression of Epcr, a marker for quiescent HSCs inversely correlated with expression of CD41 in mice, but did not show such reciprocal expression pattern in patients with MPN. Treatment with interferon-α further increased the frequency and percentage of CD41hi HSCs and reduced the number of JAK2-V617F+ HSCs in mice and patients with MPN. The shift toward the CD41hi subset of HSCs by interferon-α provides a possible mechanism of how interferon-α preferentially targets the JAK2 mutant clone.

Oncolytic virotherapy-mediated anti-tumor response: a single-cell perspective.

Talimogene laherparepvec (T-VEC) is a genetically modified herpes simplex 1 virus (HSV-1) approved for cancer therapy. We investigate its effect on the clinical, histological, single-cell transcriptomic, and immune repertoire level using repeated fine-needle aspirates (FNAs) of injected and noninjected lesions in primary cutaneous B cell lymphoma (pCBCL). Thirteen patients received intralesional T-VEC, 11 of which demonstrate tumor response in the injected lesions. Using single-cell sequencing of the FNAs, we identify the malignant population and separate three pCBCL subtypes. Twenty-four hours after the injection, we detect HSV-1T-VEC transcripts in malignant and nonmalignant cells of the injected lesion but not of the noninjected lesion. Oncolytic virotherapy results in a rapid eradication of malignant cells. It also leads to interferon pathway activation and early influx of natural killer cells, monocytes, and dendritic cells. These events are followed by enrichment in cytotoxic T cells and a decrease of regulatory T cells in injected and noninjected lesions.

The Tumor Profiler Study: integrated, multi-omic, functional tumor profiling for clinical decision support.

The application and integration of molecular profiling technologies create novel opportunities for personalized medicine. Here, we introduce the Tumor Profiler Study, an observational trial combining a prospective diagnostic approach to assess the relevance of in-depth tumor profiling to support clinical decision-making with an exploratory approach to improve the biological understanding of the disease.

Convergent network effects along the axis of gene expression during prostate cancer progression.

BACKGROUND: Tumor-specific genomic aberrations are routinely determined by high-throughput genomic measurements. It remains unclear how complex genome alterations affect molecular networks through changing protein levels and consequently biochemical states of tumor tissues. RESULTS: Here, we investigate the propagation of genomic effects along the axis of gene expression during prostate cancer progression. We quantify genomic, transcriptomic, and proteomic alterations based on 105 prostate samples, consisting of benign prostatic hyperplasia regions and malignant tumors, from 39 prostate cancer patients. Our analysis reveals the convergent effects of distinct copy number alterations impacting on common downstream proteins, which are important for establishing the tumor phenotype. We devise a network-based approach that integrates perturbations across different molecular layers, which identifies a sub-network consisting of nine genes whose joint activity positively correlates with increasingly aggressive tumor phenotypes and is predictive of recurrence-free survival. Further, our data reveal a wide spectrum of intra-patient network effects, ranging from similar to very distinct alterations on different molecular layers. CONCLUSIONS: This study uncovers molecular networks with considerable convergent alterations across tumor sites and patients. It also exposes a diversity of network effects: we could not identify a single sub-network that is perturbed in all high-grade tumor regions.

Hypoxia Triggers the Intravasation of Clustered Circulating Tumor Cells.

Circulating tumor cells (CTCs) are shed from solid cancers in the form of single or clustered cells, and the latter display an extraordinary ability to initiate metastasis. Yet, the biological phenomena that trigger the shedding of CTC clusters from a primary cancerous lesion are poorly understood. Here, when dynamically labeling breast cancer cells along cancer progression, we observe that the majority of CTC clusters are undergoing hypoxia, while single CTCs are largely normoxic. Strikingly, we find that vascular endothelial growth factor (VEGF) targeting leads to primary tumor shrinkage, but it increases intra-tumor hypoxia, resulting in a higher CTC cluster shedding rate and metastasis formation. Conversely, pro-angiogenic treatment increases primary tumor size, yet it dramatically suppresses the formation of CTC clusters and metastasis. Thus, intra-tumor hypoxia leads to the formation of clustered CTCs with high metastatic ability, and a pro-angiogenic therapy suppresses metastasis formation through prevention of CTC cluster generation.

Clonal hematopoiesis in donors and long-term survivors of related allogeneic hematopoietic stem cell transplantation.

Clonal hematopoiesis (CH) is associated with age and an increased risk of myeloid malignancies, cardiovascular risk, and all-cause mortality. We tested for CH in a setting where hematopoietic stem cells (HSCs) of the same individual are exposed to different degrees of proliferative stress and environments, ie, in long-term survivors of allogeneic hematopoietic stem cell transplantation (allo-HSCT) and their respective related donors (n = 42 donor-recipient pairs). With a median follow-up time since allo-HSCT of 16 years (range, 10-32 years), we found a total of 35 mutations in 23 out of 84 (27.4%) study participants. Ten out of 42 donors (23.8%) and 13 out of 42 recipients (31%) had CH. CH was associated with older donor and recipient age. We identified 5 cases of donor-engrafted CH, with 1 case progressing into myelodysplastic syndrome in both donor and recipient. Four out of 5 cases showed increased clone size in recipients compared with donors. We further characterized the hematopoietic system in individuals with CH as follows: (1) CH was consistently present in myeloid cells but varied in penetrance in B and T cells; (2) colony-forming units (CFUs) revealed clonal evolution or multiple independent clones in individuals with multiple CH mutations; and (3) telomere shortening determined in granulocytes suggested ∼20 years of added proliferative history of HSCs in recipients compared with their donors, with telomere length in CH vs non-CH CFUs showing varying patterns. This study provides insight into the long-term behavior of the same human HSCs and respective CH development under different proliferative conditions.

JAK2-mutant hematopoietic cells display metabolic alterations that can be targeted to treat myeloproliferative neoplasms.

Increased energy requirement and metabolic reprogramming are hallmarks of cancer cells. We show that metabolic alterations in hematopoietic cells are fundamental to the pathogenesis of mutant JAK2-driven myeloproliferative neoplasms (MPNs). We found that expression of mutant JAK2 augmented and subverted metabolic activity of MPN cells, resulting in systemic metabolic changes in vivo, including hypoglycemia, adipose tissue atrophy, and early mortality. Hypoglycemia in MPN mouse models correlated with hyperactive erythropoiesis and was due to a combination of elevated glycolysis and increased oxidative phosphorylation. Modulating nutrient supply through high-fat diet improved survival, whereas high-glucose diet augmented the MPN phenotype. Transcriptomic and metabolomic analyses identified numerous metabolic nodes in JAK2-mutant hematopoietic stem and progenitor cells that were altered in comparison with wild-type controls. We studied the consequences of elevated levels of Pfkfb3, a key regulatory enzyme of glycolysis, and found that pharmacological inhibition of Pfkfb3 with the small molecule 3PO reversed hypoglycemia and reduced hematopoietic manifestations of MPNs. These effects were additive with the JAK1/2 inhibitor ruxolitinib in vivo and in vitro. Inhibition of glycolysis by 3PO altered the redox homeostasis, leading to accumulation of reactive oxygen species and augmented apoptosis rate. Our findings reveal the contribution of metabolic alterations to the pathogenesis of MPNs and suggest that metabolic dependencies of mutant cells represent vulnerabilities that can be targeted for treating MPNs.

Neutrophils escort circulating tumour cells to enable cell cycle progression.

A better understanding of the features that define the interaction between cancer cells and immune cells is important for the development of new cancer therapies1. However, focus is often given to interactions that occur within the primary tumour and its microenvironment, whereas the role of immune cells during cancer dissemination in patients remains largely uncharacterized2,3. Circulating tumour cells (CTCs) are precursors of metastasis in several types of cancer4-6, and are occasionally found within the bloodstream in association with non-malignant cells such as white blood cells (WBCs)7,8. The identity and function of these CTC-associated WBCs, as well as the molecular features that define the interaction between WBCs and CTCs, are unknown. Here we isolate and characterize individual CTC-associated WBCs, as well as corresponding cancer cells within each CTC-WBC cluster, from patients with breast cancer and from mouse models. We use single-cell RNA sequencing to show that in the majority of these cases, CTCs were associated with neutrophils. When comparing the transcriptome profiles of CTCs associated with neutrophils against those of CTCs alone, we detect a number of differentially expressed genes that outline cell cycle progression, leading to more efficient metastasis formation. Further, we identify cell-cell junction and cytokine-receptor pairs that define CTC-neutrophil clusters, representing key vulnerabilities of the metastatic process. Thus, the association between neutrophils and CTCs drives cell cycle progression within the bloodstream and expands the metastatic potential of CTCs, providing a rationale for targeting this interaction in treatment of breast cancer.

Expression and Mutation Patterns of PBRM1, BAP1 and SETD2 Mirror Specific Evolutionary Subtypes in Clear Cell Renal Cell Carcinoma.

Bi-allelic inactivation of the VHL gene on chromosome 3p is the characteristic feature in most clear cell renal cell carcinomas (ccRCC). Frequent gene alterations were also identified in SETD2, BAP1 and PBRM1, all of which are situated on chromosome 3p and encode histone/chromatin regulators. The relationship between gene mutation, loss of protein expression and the correlations with clinicopathological parameters is important for the understanding of renal cancer progression. We analyzed PBRM1 and BAP1 protein expression as well as the tri-methylation state of H3K36 as a surrogate marker for SETD2 activity in more than 700 RCC samples. In ccRCC loss of nuclear PBRM1 (68%), BAP1 (40%) and H3K36me3 (47%) expression was significantly correlated with each other, advanced tumor stage, poor tumor differentiation (P < .0001 each), and necrosis (P < .005) Targeted next generation sequencing of 83 ccRCC samples demonstrated a significant association of genetic mutations in PBRM1, BAP1, and SETD2 with absence of PBRM1, BAP1, and HEK36me3 protein expression (P < .05, each). By assigning the protein expression patterns to evolutionary subtypes, we revealed similar clinical phenotypes as suggested by TRACERx Renal. Given their important contribution to tumor suppression, we conclude that combined functional inactivation of PBRM1, BAP1, SETD2 and pVHL is critical for ccRCC progression.

Circulating Tumor Cell Clustering Shapes DNA Methylation to Enable Metastasis Seeding.

The ability of circulating tumor cells (CTCs) to form clusters has been linked to increased metastatic potential. Yet biological features and vulnerabilities of CTC clusters remain largely unknown. Here, we profile the DNA methylation landscape of single CTCs and CTC clusters from breast cancer patients and mouse models on a genome-wide scale. We find that binding sites for stemness- and proliferation-associated transcription factors are specifically hypomethylated in CTC clusters, including binding sites for OCT4, NANOG, SOX2, and SIN3A, paralleling embryonic stem cell biology. Among 2,486 FDA-approved compounds, we identify Na+/K+ ATPase inhibitors that enable the dissociation of CTC clusters into single cells, leading to DNA methylation remodeling at critical sites and metastasis suppression. Thus, our results link CTC clustering to specific changes in DNA methylation that promote stemness and metastasis and point to cluster-targeting compounds to suppress the spread of cancer.