Shifting patterns of genomic variation in the somatic evolution of papillary thyroid carcinoma

Jill C. Rubinstein, Taylor C. Brown, Emily R. Christison-Lagay, Yawei Zhang, John W. Kunstman, C. Christofer Juhlin, Carol Nelson-Williams, Gerald Goh, Courtney E. Quinn, Glenda G. Callender, Robert Udelsman, Richard P. Lifton, Reju Korah, Tobias Carling

Cancer is increasingly understood to arise in the context of dynamically evolving genomes with continuously generated variants subject to selective pressures. Diverse mutations have been identified in papillary thyroid carcinoma (PTC), but unifying theories underlying genomic change are lacking. Applying a framework of somatic evolution, we sought to broaden understanding of the PTC genome through identification of global trends that help explain risk of tumorigenesis. Exome sequencing was performed on 53 PTC and matched adjacent non-tumor thyroid tissues (ANT). Single nucleotide substitution (SNS) signatures from each sample pair were divided into three subsets based on their presence in tumor, non-tumor thyroid, or both. Nine matched blood samples were sequenced and SNS signatures intersected with these three subsets. The intersected genomic signatures were used to define branch-points in the evolution of the tumor genome, distinguishing variants present in the tissues' common ancestor cells from those unique to each tissue type and therefore acquired after genomic divergence of the tumor, non-tumor, and blood samples. Single nucleotide substitutions shared by the tumor and the non-tumor thyroid were dominated by C-to-T transitions, whereas those unique to either tissue type were enriched for C-to-A transversions encoding non-synonymous, predicted-deleterious variants. On average, SNSs of matched blood samples were 81 % identical to those shared by tumor and non-tumor thyroid, but only 12.5 % identical to those unique to either tissue. Older age and BRAF mutation were associated with increased SNS burden. The current study demonstrates novel patterns of genomic change in PTC, supporting a theory of somatic evolution in which the zygote's germline genome undergoes continuous remodeling to produce progressively differentiated, tissue-specific signatures. Late somatic events in thyroid tissue demonstrate shifted mutational spectra compared to earlier polymorphisms. These late events are enriched for predicted-deleterious variants, suggesting a mechanism of genomic instability in PTC tumorigenesis.