Copy number variants (CNVs) are a pervasive source of genetic variation and evolutionary potential, but the dynamics and diversity of CNVs within evolving populations remain unclear. Long-term evolution experiments in chemostats provide an ideal system for studying the molecular processes underlying CNV formation and the temporal dynamics with which they are generated, selected, and maintained. Here, we developed a fluorescent CNV reporter to detect de novo gene amplifications and deletions in individual cells. We used the CNV reporter in Saccharomyces cerevisiae to study CNV formation at the GAP1 locus, which encodes the general amino acid permease, in different nutrient-limited chemostat conditions. We find that under strong selection, GAP1 CNVs are repeatedly generated and selected during the early stages of adaptive evolution, resulting in predictable dynamics. Molecular characterization of CNV-containing lineages shows that the CNV reporter detects different classes of CNVs, including aneuploidies, nonreciprocal translocations, tandem duplications, and complex CNVs. Despite GAP1’s proximity to repeat sequences that facilitate intrachromosomal recombination, breakpoint analysis revealed that short inverted repeat sequences mediate formation of at least 50% of GAP1 CNVs. Inverted repeat sequences are also found at breakpoints at the DUR3 locus, where CNVs are selected in urea-limited chemostats. Analysis of 28 CNV breakpoints indicates that inverted repeats are typically 8 nucleotides in length and separated by 40 bases. The features of these CNVs are consistent with origin-dependent inverted-repeat amplification (ODIRA), suggesting that replication-based mechanisms of CNV formation may be a common source of gene amplification. We combined the CNV reporter with barcode lineage tracking and found that 102–104 independent CNV-containing lineages initially compete within populations, resulting in extreme clonal interference. However, only a small number (18–21) of CNV lineages ever constitute more than 1% of the CNV subpopulation, and as selection progresses, the diversity of CNV lineages declines. Our study introduces a novel means of studying CNVs in heterogeneous cell populations and provides insight into their dynamics, diversity, and formation mechanisms in the context of adaptive evolution.

CTCF and cohesin play a key role in organizing chromatin into topologically associating domain (TAD) structures. Disruption of a single CTCF binding site is sufficient to change chromosomal interactions leading to alterations in chromatin modifications and gene regulation. However, the extent to which alterations in chromatin modifications can disrupt 3D chromosome organization leading to transcriptional changes is unknown. In multiple myeloma, a 4;14 translocation induces overexpression of the histone methyltransferase, NSD2, resulting in expansion of H3K36me2 and shrinkage of antagonistic H3K27me3 domains. Using isogenic cell lines producing high and low levels of NSD2, here we find oncogene activation is linked to alterations in H3K27ac and CTCF within H3K36me2 enriched chromatin. A logistic regression model reveals that differentially expressed genes are significantly enriched within the same insulated domain as altered H3K27ac and CTCF peaks. These results identify a bidirectional relationship between 2D chromatin and 3D genome organization in gene regulation.

CTCF and cohesin play a key role in organizing chromatin into TAD structures. Disruption of a single CTCF binding site is sufficient to change chromosomal interactions leading to alterations in chromatin modifications and gene regulation. However, the extent to which alterations in chromatin modifications can disrupt 3D chromosome organization leading to transcriptional changes is unknown. In multiple myeloma a 4;14 translocation induces overexpression of the histone methyltransferase, NSD2 resulting in expansion of H3K36me2 and shrinkage of antagonistic H3K27me3 domains. Using isogenic cell lines producing high and low levels of NSD2, we find oncogene activation is linked to alterations in H3K27ac and CTCF within H3K36me2 enriched chromatin. A linear regression model reveals that changes in both CTCF and/or H3K27ac significantly increase the probability that a gene sharing the same insulated domain will be differentially expressed. These results identify a bidirectional relationship between 2D chromatin and 3D genome organization in gene regulation.

The NSD2 p.E1099K (EK) mutation is observed in 10% of acute lymphoblastic leukemia (ALL) samples with enrichment at relapse indicating a role in clonal evolution and drug resistance. To discover mechanisms that mediate clonal expansion, we engineered B-ALL cell lines (Reh, 697) to overexpress wildtype (WT) and EK NSD2, but observed no differences in proliferation, clonal growth, or chemosensitivity. To address whether NSD2 EK acts collaboratively with other pathways, we used shRNAs to knockdown expression of NSD2 in B-ALL cell lines heterozygous for NSD2 EK (RS4;11, RCH-ACV, SEM). Knockdown resulted in decreased proliferation in all lines, decreased clonal growth in RCH-ACV, and increased sensitivity to cytotoxic chemotherapeutic agents, although the pattern of drug sensitivity varied among cell lines implying that the oncogenic properties of NSD2 mutations are likely cell context specific and rely on cooperative pathways. Knockdown of both Type II and REIIBP EK isoforms had a greater impact than knockdown of Type II alone, suggesting that both SET containing EK isoforms contribute to phenotypic changes driving relapse. Furthermore, in vivo models using both cell lines and patient samples revealed dramatically enhanced proliferation of NSD2 EK compared to WT and reduced sensitivity to 6-mercaptopurine in the relapse sample relative to diagnosis. Finally, EK-mediated changes in chromatin state and transcriptional output differed dramatically among cell lines further supporting a cell context specific role of NSD2 EK. These results demonstrate a unique role of NSD2 EK in mediating clonal fitness through pleiotropic mechanisms dependent on the genetic and epigenetic landscape. IMPLICATIONSNSD2 p.E1099K mutation leads to drug resistance and a clonal advantage in childhood B-ALL.

Activation of regulatory elements is thought to be inversely correlated with DNA methylation levels. However, it is difficult to determine whether DNA methylation is compatible with chromatin accessibility or transcription factor (TF) binding if assays are performed separately. We developed a fast, low-input, low sequencing depth method, EpiMethylTag, that combines ATAC-seq or ChIP-seq (M-ATAC or M-ChIP) with bisulfite conversion, to simultaneously examine accessibility/TF binding and methylation on the same DNA. Here we demonstrate that EpiMethylTag can be used to study the functional interplay between chromatin accessibility and TF binding (CTCF and KLF4) at methylated sites.

Copy number variants (CNVs) are a pervasive, but understudied source of genetic variation and evolutionary potential. Long-term evolution experiments in chemostats provide an ideal system for studying the molecular processes underlying CNV formation and the temporal dynamics of de novo CNVs. Here, we developed a fluorescent reporter to monitor gene amplifications and deletions at a specific locus with single-cell resolution. Using a CNV reporter in nitrogen-limited chemostats, we find that GAP1 CNVs are repeatedly generated and selected during the early stages of adaptive evolution resulting in predictable dynamics of CNV selection. However, subsequent diversification of populations defines a second phase of evolutionary dynamics that cannot be predicted. Using whole genome sequencing, we identified a variety of GAP1 CNVs that vary in size and copy number. Despite GAP1 's proximity to tandem repeats that facilitate intrachromosomal recombination, we find that non-allelic homologous recombination (NAHR) between flanking tandem repeats occurs infrequently. Rather, breakpoint characterization revealed that for at least 50% of GAP1 CNVs, origin-dependent inverted-repeat amplification (ODIRA), a DNA replication mediated process, is the likely mechanism. We also find evidence that ODIRA generates DUR3 CNVs, indicating that it may be a common mechanism of gene amplification. We combined the CNV reporter with barcode lineage tracking and found that 10 3 -10 4 independent CNV-containing lineages initially compete within populations, which results in extreme clonal interference. Our study introduces a novel means of studying CNVs in heterogeneous cell populations and provides insight into the underlying dynamics of CNVs in evolution.condition [32,33] . A high rate of CNV formation suggests that multiple, independent CNV-containing lineages may compete during adaptive evolution resulting in clonal interference, which is characteristic of large, evolving populations [29,[34][35][36] . However, the extent of clonal interference among CNV-containing lineages is unknown.The general amino acid permease, GAP1 , is ideally suited to studying the role of CNVs in adaptive evolution. GAP1 encodes a high-affinity transporter for all naturally occurring amino acids and analogues, and it is highly expressed in nitrogen-poor conditions [37,38] . We have previously shown that two classes of CNVs are selected at the GAP1 locus in S. cerevisiae : amplification alleles in glutamine and glutamate-limited chemostats and deletion alleles in ureaand allantoin-limited chemostats [24,25] . GAP1 CNVs are also found in natural populations.Multiple, tandem copies of GAP1 have been identified in wild populations of the nectar yeast, Metschnikowia reukaufii , which result in a competitive advantage over other microbes when amino acids are scarce [39] . As a frequent target of selection in adverse environments in both experimental and natural populations, GAP1 is a model locus for studying the dynamics and mechanisms underlying both gene amplification a...

Although B-cell acute lymphoblastic leukemia (B-ALL) is the most common malignancy in children and while highly curable, it remains a leading cause of cancer-related mortality. The outgrowth of tumor subclones carrying mutations in genes responsible for resistance to therapy has led to a Darwinian model of clonal selection. Previous work has indicated that alterations in the epigenome might contribute to clonal selection, yet the extent to which the chromatin state is altered under the selective pressures of therapy is unknown. To address this, we performed chromatin immunoprecipitation, gene expression analysis, and enhanced reduced representation bisulfite sequencing on a cohort of paired diagnosis and relapse samples from individual patients who all but one relapsed within 36 months of initial diagnosis. The chromatin state at diagnosis varied widely among patients, while the majority of peaks remained stable between diagnosis and relapse. Yet a significant fraction was either lost or newly gained, with some patients showing few differences and others showing massive changes of the epigenetic state. Evolution of the epigenome was associated with pathways previously linked to therapy resistance as well as novel candidate pathways through alterations in pyrimidine biosynthesis and downregulation of polycomb repressive complex 2 targets. Three novel, relapse-specific superenhancers were shared by a majority of patients including one associated with S100A8, the top upregulated gene seen at relapse in childhood B-ALL. Overall, our results support a role of the epigenome in clonal evolution and uncover new candidate pathways associated with relapse.Significance: This study suggests a major role for epigenetic mechanisms in driving clonal evolution in B-ALL and identifies novel pathways associated with drug resistance.

Activation of regulatory elements is thought to be inversely correlated with DNA methylation levels. However, it is difficult to determine whether DNA methylation is compatible with chromatin accessibility or transcription factor (TF) binding if assays are performed separately. We developed a low input, low sequencing depth method, EpiMethylTag that combines ATAC-seq or ChIP-seq (M-ATAC or M-ChIP) with bisulfite conversion, to simultaneously examine accessibility/TF binding and methylation on the same DNA.