Immunotherapies that block inhibitory checkpoint receptors on T cells have transformed the clinical care of cancer patients 1 . However, whether the T cell response to checkpoint blockade relies on reinvigoration of pre-existing tumor infiltrating T cells (TILs) or on recruitment of novel T cells remains unclear 2 – 4 . Here, we performed paired single-cell RNA (scRNA) and T cell receptor (TCR)- sequencing on 79,046 cells from site-matched tumors from patients with basal cell carcinoma (BCC) or squamous cell carcinoma (SCC) pre- and post-anti-PD-1 therapy. Tracking TCR clones and transcriptional phenotypes revealed a coupling of tumor-recognition, clonal expansion, and T cell dysfunction marked by clonal expansions of CD8 + CD39 + T cells, which co-expressed markers of chronic T cell activation and exhaustion. However, this expansion did not derive from pre-existing TIL clones; rather, it was comprised of novel clonotypes not previously observed in the same tumor. Clonal replacement of T cells was preferentially observed in exhausted CD8 + T cells and evident in BCC and SCC patients. These results demonstrate that pre-existing tumor-specific T cells may have limited reinvigoration capacity, and that the T cell response to checkpoint blockade derives from a distinct repertoire of T cell clones that may have just recently entered the tumor.

In this multicenter, open-label, randomized phase II investigator-sponsored neoadjuvant trial with funding provided by Sanofi and GlaxoSmithKline (TRIO-US B07, Clinical Trials NCT00769470), participants with early-stage HER2-positive breast cancer (N = 128) were recruited from 13 United States oncology centers throughout the Translational Research in Oncology network. Participants were randomized to receive trastuzumab (T; N = 34), lapatinib (L; N = 36), or both (TL; N = 58) as HER2-targeted therapy, with each participant given one cycle of this designated anti-HER2 therapy alone followed by six cycles of standard combination chemotherapy with the same anti-HER2 therapy. The primary objective was to estimate the rate of pathologic complete response (pCR) at the time of surgery in each of the three arms. In the intent-to-treat population, we observed similar pCR rates between T (47%, 95% confidence interval [CI] 30–65%) and TL (52%, 95% CI 38–65%), and a lower pCR rate with L (25%, 95% CI 13–43%). In the T arm, 100% of participants completed all protocol-specified treatment prior to surgery, as compared to 69% in the L arm and 74% in the TL arm. Tumor or tumor bed tissue was collected whenever possible pre-treatment (N = 110), after one cycle of HER2-targeted therapy alone (N = 89), and at time of surgery (N = 59). Higher-level amplification of HER2 and hormone receptor (HR)-negative status were associated with a higher pCR rate. Large shifts in the tumor, immune, and stromal gene expression occurred after one cycle of HER2-targeted therapy. In contrast to pCR rates, the L-containing arms exhibited greater proliferation reduction than T at this timepoint. Immune expression signatures increased in all arms after one cycle of HER2-targeted therapy, decreasing again by the time of surgery. Our results inform approaches to early assessment of sensitivity to anti-HER2 therapy and shed light on the role of the immune microenvironment in response to HER2-targeted agents.

Copy number alterations (CNAs) play an important role in molding the genomes of breast cancers and have been shown to be clinically useful for prognostic and therapeutic purposes. However, our knowledge of intra-tumoral genetic heterogeneity of this important class of somatic alterations is limited. Here, using single-cell sequencing, we comprehensively map out the facets of copy number alteration heterogeneity in a cohort of breast cancer tumors. Ou/var/www/html/elife/12-05-2020/backup/r analyses reveal: genetic heterogeneity of non-tumor cells (i.e. stroma) within the tumor mass; the extent to which copy number heterogeneity impacts breast cancer genomes and the importance of both the genomic location and dosage of sub-clonal events; the pervasive nature of genetic heterogeneity of chromosomal amplifications; and the association of copy number heterogeneity with clinical and biological parameters such as polyploidy and estrogen receptor negative status. Our data highlight the power of single-cell genomics in dissecting, in its many forms, intra-tumoral genetic heterogeneity of CNAs, the magnitude with which CNA heterogeneity affects the genomes of breast cancers, and the potential importance of CNA heterogeneity in phenomena such as therapeutic resistance and disease relapse.

Genomic changes observed across treatment may result from either clonal evolution or geographically disparate sampling of heterogeneous tumors. Here we use computational modeling based on analysis of fifteen primary breast tumors and find that apparent clonal change between two tumor samples can frequently be explained by pre-treatment heterogeneity, such that at least two regions are necessary to detect treatment-induced clonal shifts. To assess for clonal replacement, we devise a summary statistic based on whole-exome sequencing of a pre-treatment biopsy and multi-region sampling of the post-treatment surgical specimen and apply this measure to five breast tumors treated with neoadjuvant HER2-targeted therapy. Two tumors underwent clonal replacement with treatment, and mathematical modeling indicates these two tumors had resistant subclones prior to treatment and rates of resistance-related genomic changes that were substantially larger than previous estimates. Our results provide a needed framework to incorporate primary tumor heterogeneity in investigating the evolution of resistance.

doi: medRxiv preprint NOTE: This preprint reports new research that has not been certified by peer review and should not be used to guide clinical practice.