Abstract

The processes by which tumors evolve are essential to the efficacy of treatment, but quantitative understanding of intratumoral dynamics has been limited. Although intratumoral heterogeneity is common, quantification of evolution is difficult from clinical samples because treatment replicates cannot be performed and because matched serial samples are infrequently available. To circumvent these problems we derived and assayed large sets of human triple-negative breast cancer xenografts and cell cultures from two patients, including 86 xenografts from cyclophosphamide, doxorubicin, cisplatin, docetaxel, or vehicle treatment cohorts as well as 45 related cell cultures. We assayed these samples via exome-seq and/or high-resolution droplet digital PCR, allowing us to distinguish complex therapy-induced selection and drift processes among endogenous cancer subclones with cellularity uncertainty <3%. For one patient, we discovered two predominant subclones that were granularly intermixed in all 48 co-derived xenograft samples. These two subclones exhibited differential chemotherapy sensitivity -- when xenografts were treated with cisplatin for 3 weeks, the post-treatment volume change was proportional to the post-treatment ratio of subclones on a xenograft-to-xenograft basis. A subsequent cohort in which xenografts were treated with cisplatin, allowed a drug holiday, then treated a second time continued to exhibit this proportionality. In contrast, xenografts from other treatment cohorts, spatially dissected xenograft fragments, and cell cultures evolved unsystematically but with substantial population bottlenecks. These results show that ecologies susceptible to successive retreatment can arise spontaneously in breast cancer in spite of a background of irregular subclonal bottlenecks, and our work provides to our knowledge the first quantification of the population genetics of such a system. Intriguingly, in such an ecology the ratio of common subclones is predictive of the state of treatment susceptibility, suggesting that this ratio can be measured to optimize dynamic treatment protocols in patients.\n\nAUTHOR SUMMARYAn overarching challenge of cancer is that patients develop resistance to treatment -- an essentially evolutionary process. However, there is currently very little understanding of how tumor evolution can be exploited to improve treatment. One reason for this is that usually only 1-2 samples can be obtained per patient, so cancer evolutionary processes are still poorly understood. To solve this problem, we created many dozens of copies of the tumors from two breast cancer patients using xenografting and cell culture methods. We then compared the evolution in these tumor copies in response to different treatments, including four of the most common breast cancer chemotherapies. These studies present the most exhaustive comparisons of treatment-induced evolution that have yet been performed for individual cancer patients. Unexpectedly, high-resolution sequencing of these samples revealed a special dynamically treatable ecology in one tumor, in which tumor growth during platinum therapy was determined by the ecological balance of two tumor cell populations. Our work shows that ecologies that can be targeted by dynamic treatment strategies arise spontaneously in breast cancers. Population heterogeneity is common within cancers, and our work suggests how tracking of intratumoral evolution can be used to optimize treatment.