Single-cell decoding of drug induced transcriptomic reprogramming in triple negative breast cancers.

BACKGROUND: The encoding of cell intrinsic drug resistance states in breast cancer reflects the contributions of genomic and non-genomic variations and requires accurate estimation of clonal fitness from co-measurement of transcriptomic and genomic data. Somatic copy number (CN) variation is the dominant mutational mechanism leading to transcriptional variation and notably contributes to platinum chemotherapy resistance cell states. Here, we deploy time series measurements of triple negative breast cancer (TNBC) single-cell transcriptomes, along with co-measured single-cell CN fitness, identifying genomic and transcriptomic mechanisms in drug-associated transcriptional cell states. RESULTS: We present scRNA-seq data (53,641 filtered cells) from serial passaging TNBC patient-derived xenograft (PDX) experiments spanning 2.5 years, matched with genomic single-cell CN data from the same samples. Our findings reveal distinct clonal responses within TNBC tumors exposed to platinum. Clones with high drug fitness undergo clonal sweeps and show subtle transcriptional reversion, while those with weak fitness exhibit dynamic transcription upon drug withdrawal. Pathway analysis highlights convergence on epithelial-mesenchymal transition and cytokine signaling, associated with resistance. Furthermore, pseudotime analysis demonstrates hysteresis in transcriptional reversion, indicating generation of new intermediate transcriptional states upon platinum exposure. CONCLUSIONS: Within a polyclonal tumor, clones with strong genotype-associated fitness under platinum remained fixed, minimizing transcriptional reversion upon drug withdrawal. Conversely, clones with weaker fitness display non-genomic transcriptional plasticity. This suggests CN-associated and CN-independent transcriptional states could both contribute to platinum resistance. The dominance of genomic or non-genomic mechanisms within polyclonal tumors has implications for drug sensitivity, restoration, and re-treatment strategies.

Luminal breast epithelial cells from wildtype and BRCA mutation carriers harbor copy number alterations commonly associated with breast cancer.

Cancer-associated mutations have been documented in normal tissues, but the prevalence and nature of somatic copy number alterations and their role in tumor initiation and evolution is not well understood. Here, using single cell DNA sequencing, we describe the landscape of CNAs in >42,000 breast epithelial cells from women with normal or high risk of developing breast cancer. Accumulation of individual cells with one or two of a specific subset of CNAs (e.g. 1q gain and 16q, 22q, 7q, and 10q loss) is detectable in almost all breast tissues and, in those from BRCA1 or BRCA2 mutations carriers, occurs prior to loss of heterozygosity (LOH) of the wildtype alleles. These CNAs, which are among the most common associated with ductal carcinoma in situ (DCIS) and malignant breast tumors, are enriched almost exclusively in luminal cells not basal myoepithelial cells. Allele-specific analysis of the enriched CNAs reveals that each allele was independently altered, demonstrating convergent evolution of these CNAs in an individual breast. Tissues from BRCA1 or BRCA2 mutation carriers contain a small percentage of cells with extreme aneuploidy, featuring loss of TP53 , LOH of BRCA1 or BRCA2 , and multiple breast cancer-associated CNAs in addition to one or more of the common CNAs in 1q, 10q or 16q. Notably, cells with intermediate levels of CNAs are not detected, arguing against a stepwise gradual accumulation of CNAs. Overall, our findings demonstrate that chromosomal alterations in normal breast epithelium partially mirror those of established cancer genomes and are chromosome- and cell lineage-specific.

Single-cell genomic variation induced by mutational processes in cancer.

How cell-to-cell copy number alterations that underpin genomic instability1 in human cancers drive genomic and phenotypic variation, and consequently the evolution of cancer2, remains understudied. Here, by applying scaled single-cell whole-genome sequencing3 to wild-type, TP53-deficient and TP53-deficient;BRCA1-deficient or TP53-deficient;BRCA2-deficient mammary epithelial cells (13,818 genomes), and to primary triple-negative breast cancer (TNBC) and high-grade serous ovarian cancer (HGSC) cells (22,057 genomes), we identify three distinct 'foreground' mutational patterns that are defined by cell-to-cell structural variation. Cell- and clone-specific high-level amplifications, parallel haplotype-specific copy number alterations and copy number segment length variation (serrate structural variations) had measurable phenotypic and evolutionary consequences. In TNBC and HGSC, clone-specific high-level amplifications in known oncogenes were highly prevalent in tumours bearing fold-back inversions, relative to tumours with homologous recombination deficiency, and were associated with increased clone-to-clone phenotypic variation. Parallel haplotype-specific alterations were also commonly observed, leading to phylogenetic evolutionary diversity and clone-specific mono-allelic expression. Serrate variants were increased in tumours with fold-back inversions and were highly correlated with increased genomic diversity of cellular populations. Together, our findings show that cell-to-cell structural variation contributes to the origins of phenotypic and evolutionary diversity in TNBC and HGSC, and provide insight into the genomic and mutational states of individual cancer cells.

Accurate determination of CRISPR-mediated gene fitness in transplantable tumours.

Assessing tumour gene fitness in physiologically-relevant model systems is challenging due to biological features of in vivo tumour regeneration, including extreme variations in single cell lineage progeny. Here we develop a reproducible, quantitative approach to pooled genetic perturbation in patient-derived xenografts (PDXs), by encoding single cell output from transplanted CRISPR-transduced cells in combination with a Bayesian hierarchical model. We apply this to 181 PDX transplants from 21 breast cancer patients. We show that uncertainty in fitness estimates depends critically on the number of transplant cell clones and the variability in clone sizes. We use a pathway-directed allelic series to characterize Notch signaling, and quantify TP53 / MDM2 drug-gene conditional fitness in outlier patients. We show that fitness outlier identification can be mirrored by pharmacological perturbation. Overall, we demonstrate that the gene fitness landscape in breast PDXs is dominated by inter-patient differences.

Results of the phase I CCTG IND.231 trial of CX-5461 in patients with advanced solid tumors enriched for DNA-repair deficiencies.

CX-5461 is a G-quadruplex stabilizer that exhibits synthetic lethality in homologous recombination-deficient models. In this multicentre phase I trial in patients with solid tumors, 40 patients are treated across 10 dose levels (50-650 mg/m2) to determine the recommended phase II dose (primary outcome), and evaluate safety, tolerability, pharmacokinetics (secondary outcomes). Defective homologous recombination is explored as a predictive biomarker of response. CX-5461 is generally well tolerated, with a recommended phase II dose of 475 mg/m2 days 1, 8 and 15 every 4 weeks, and dose limiting phototoxicity. Responses are observed in 14% of patients, primarily in patients with defective homologous recombination. Reversion mutations in PALB2 and BRCA2 are detected on progression following initial response in germline carriers, confirming the underlying synthetic lethal mechanism. In vitro characterization of UV sensitization shows this toxicity is related to the CX-5461 chemotype, independent of G-quadruplex synthetic lethality. These results establish clinical proof-of-concept for this G-quadruplex stabilizer. Clinicaltrials.gov NCT02719977.

Iron- and Hepcidin-Independent Downregulation of the Iron Exporter Ferroportin in Macrophages during Salmonella Infection.

Retention of iron in tissue macrophages via upregulation of hepcidin (HAMP) and downregulation of the iron exporter ferroportin (FPN) is thought to participate in the establishment of anemia of inflammation after infection. However, an upregulation of FPN has been proposed to limit macrophages iron access to intracellular pathogens. Therefore, we studied the iron homeostasis and in particular the regulation of FPN after infection with Salmonella enterica serovar Typhimurium in mice presenting tissue macrophages with high iron (AcB61), basal iron (A/J and wild-type mice), or low iron (Hamp knock out, Hamp-/-) levels. The presence of iron in AcB61 macrophages due to extravascular hemolysis and strong erythrophagocytosis activity favored the proliferation of Salmonella in the spleen and liver with a concomitant decrease of FPN protein expression. Despite systemic iron overload, no or slight increase in Salmonella burden was observed in Hamp-/- mice compared to controls. Importantly, FPN expression at both mRNA and protein levels was strongly decreased during Salmonella infection in Hamp-/- mice. The repression of Fpn mRNA was also observed in Salmonella-infected cultured macrophages. In addition, the downregulation of FPN was associated with decreased iron stores in both the liver and spleen in infected mice. Our findings show that during Salmonella infection, FPN is repressed through an iron and hepcidin-independent mechanism. Such regulation likely provides the cellular iron indispensable for the growth of Salmonella inside the macrophages.