Theoretical framework of population genetics with somatic mutations taken into account: application to copy number variations in humans

Traditionally, population genetics focuses on the dynamics of frequencies of alleles acquired by mutations on germ-lines, because only such mutations are heritable. Typical genotyping experiments, however, use DNA from some somatic tissues such as blood, which harbors somatic mutations at the current generation in addition to germ-line mutations accumulated since the most recent common ancestor of the sample. This common practice may sometimes cause erroneous interpretations of polymorphism data, unless we properly understand the role of somatic mutations in population genetics. We here introduce a very basic theoretical framework of population genetics with somatic mutations taken into account. It is easy to imagine that somatic mutations at the current generation simply add individual-specific variations, as errors in mutation detection do. Our theory quantifies this increment under various conditions. We find that the major contribution of somatic mutations plus errors is to very rare variants, particularly to singletons. The relative contribution is markedly large when mutations are deleterious. Because negative selection also increases rare variants, it is important to distinguish the roles of these mutually confounding factors when we interpret the data, even after correcting for demography. We apply this theory to human copy number variations (CNVs), for which the composite effect of somatic mutations and errors may not be negligible. Using genome-wide CNV data, we demonstrate how the joint action of the two factors, selection and somatic mutations plus errors, shapes the observed pattern of polymorphism