The influence of genomic context on mutation patterns in the human genome inferred from rare variants

Understanding patterns of spontaneous mutations is of fundamental interest in studies of human genome evolution and genetic disease. Here, we used extremely rare variants in humans to model the molecular spectrum of single-nucleotide mutations. Compared to common variants in humans and human-chimpanzee fixed differences (substitutions), rare variants, on average, arose more recently in the human lineage and are less affected by the potentially confounding effects of natural selection, population demographic history, and biased gene conversion. We analyzed variants obtained from a population-based sequencing study of 202 genes in >14,000 individuals. We observed considerable variability in the per-gene mutation rate, which was correlated with local GC content, but not recombination rate. Using >20,000 variants with a derived allele frequency >= 10(-4), we examined the effect of local GC content and recombination rate on individual variant subtypes and performed comparisons with common variants and substitutions. The influence of local GC content on rare variants differed from that on common variants or substitutions, and the differences varied by variant subtype. Furthermore, recombination rate and recombination hotspots have little effect on rare variants of any subtype, yet both have a relatively strong impact on multiple variant subtypes in common variants and substitutions. This observation is consistent with the effect of biased gene conversion or selection-dependent processes. Our results highlight the distinct biases inherent in the initial mutation patterns and subsequent evolutionary processes that affect segregating variants

Replication Stress Induces Genome-wide Copy Number Changes in Human Cells that Resemble Polymorphic and Pathogenic Variants

Copy number variants (CNVs) are an important component of genomic variation in humans and other mammals. Similar de novo deletions and duplications, or copy number changes (CNCs), are now known to be a major cause of genetic and developmental disorders and to arise somatically in many cancers. A major mechanism leading to both CNVs and disease-associated CNCs is meiotic unequal crossing over, or nonallelic homologous recombination (NAHR), mediated by flanking repeated sequences or segmental duplications. Others appear to involve nonhomologous end joining (NHEJ) or aberrant replication suggesting a mitotic cell origin. Here we show that aphidicolin-induced replication stress in normal human cells leads to a high frequency of CNCs of tens to thousands of kilobases across the human genome that closely resemble CNVs and disease-associated CNCs. Most deletion and duplication breakpoint junctions were characterized by short (<6 bp) microhomologies, consistent with the hypothesis that these rearrangements were formed by NHEJ or a replication-coupled process, such as template switching. This is a previously unrecognized consequence of replication stress and suggests that replication fork stalling and subsequent error-prone repair are important mechanisms in the formation of CNVs and pathogenic CNCs in humans