Purpose: Recent studies have highlighted the existence of subclones in tumors. Lymph nodes are generally the first location of metastasis for most solid epithelial tumors, including colorectal cancer. We sought to understand the genetic origin of lymph node metastasis in colorectal cancer by evaluating the relationship between colorectal cancer subclones present in primary tumors and lymph nodes.Experimental Design: A total of 33 samples from seven colorectal cancers, including two or three spatially disparate regions from each primary tumor and one to four matched lymph nodes for each tumor, underwent next-generation whole-exome DNA sequencing, Affymetrix OncoScan SNP arrays, and targeted deep confirmatory sequencing. We performed mapping between SNPs and copy number events from the primary tumor and matched lymph node samples, allowing us to profile heterogeneity and the mutational origin of lymph node metastases. The computational method PyClone was used to define subclones within each tumor. The method Clonality Inference in Tumors Using Phylogeny (CITUP) was subsequently used to infer phylogenetic relationships among subclones.Results: We found that there was substantial heterogeneity in mutations and copy number changes among all samples from any given patient. For each patient, the primary tumor regions and matched lymph node metastases were each polyclonal, and the clonal populations differed from one lymph node to another. In some patients, the cancer cell populations in a given lymph node originated from multiple distinct regions of a tumor.Conclusions: Our data support a model of lymph node metastatic spread in colorectal cancer whereby metastases originate from multiple waves of seeding from the primary tumor over time. Clin Cancer Res; 24(9); 2214-24. ©2017 AACRSee related commentary by Gerlinger, p. 2032.