High-Resolution Recombination Patterns in a Region of Human Chromosome 21 Measured by Sperm Typing

For decades, classical crossover studies and linkage disequilibrium (LD) analysis of genomic regions suggested that human meiotic crossovers may not be randomly distributed along chromosomes but are focused instead in “hot spots.” Recent sperm typing studies provided data at very high resolution and accuracy that defined the physical limits of a number of hot spots. The data were also used to test whether patterns of LD can predict hot spot locations. These sperm typing studies focused on several small regions of the genome already known or suspected of containing a hot spot based on the presence of LD breakdown or previous experimental evidence of hot spot activity. Comparable data on target regions not specifically chosen using these two criteria is lacking but is needed to make an unbiased test of whether LD data alone can accurately predict active hot spots. We used sperm typing to estimate recombination in 17 almost contiguous ~5 kb intervals spanning 103 kb of human Chromosome 21. We found two intervals that contained new hot spots. The comparison of our data with recombination rates predicted by statistical analyses of LD showed that, overall, the two datasets corresponded well, except for one predicted hot spot that showed little crossing over. This study doubles the experimental data on recombination in men at the highest resolution and accuracy and supports the emerging genome-wide picture that recombination is localized in small regions separated by cold areas. Detailed study of one of the new hot spots revealed a sperm donor with a decrease in recombination intensity at the canonical recombination site but an increase in crossover activity nearby. This unique finding suggests that the position and intensity of hot spots may evolve by means of a concerted mechanism that maintains the overall recombination intensity in the region

Family reunion via error correction: an efficient analysis of duplex sequencing data

Background Duplex sequencing is the most accurate approach for identification of sequence variants present at very low frequencies. Its power comes from pooling together multiple descendants of both strands of original DNA molecules, which allows distinguishing true nucleotide substitutions from PCR amplification and sequencing artifacts. This strategy comes at a cost—sequencing the same molecule multiple times increases dynamic range but significantly diminishes coverage, making whole genome duplex sequencing prohibitively expensive. Furthermore, every duplex experiment produces a substantial proportion of singleton reads that cannot be used in the analysis and are thrown away. Results In this paper we demonstrate that a significant fraction of these reads contains PCR or sequencing errors within duplex tags. Correction of such errors allows “reuniting” these reads with their respective families increasing the output of the method and making it more cost effective. Conclusions We combine an error correction strategy with a number of algorithmic improvements in a new version of the duplex analysis software, Du Novo 2.0. It is written in Python, C, AWK, and Bash. It is open source and readily available through Galaxy, Bioconda, and Github: https://github.com/galaxyproject/dunovo

The Molecular Anatomy of Spontaneous Germline Mutations in Human Testes

The frequency of the most common sporadic Apert syndrome mutation (C755G) in the human fibroblast growth factor receptor 2 gene (FGFR2) is 100–1,000 times higher than expected from average nucleotide substitution rates based on evolutionary studies and the incidence of human genetic diseases. To determine if this increased frequency was due to the nucleotide site having the properties of a mutation hot spot, or some other explanation, we developed a new experimental approach. We examined the spatial distribution of the frequency of the C755G mutation in the germline by dividing four testes from two normal individuals each into several hundred pieces, and, using a highly sensitive PCR assay, we measured the mutation frequency of each piece. We discovered that each testis was characterized by rare foci with mutation frequencies 103 to >104 times higher than the rest of the testis regions. Using a model based on what is known about human germline development forced us to reject (p < 10−6) the idea that the C755G mutation arises more frequently because this nucleotide simply has a higher than average mutation rate (hot spot model). This is true regardless of whether mutation is dependent or independent of cell division. An alternate model was examined where positive selection acts on adult self-renewing Ap spermatogonial cells (SrAp) carrying this mutation such that, instead of only replacing themselves, they occasionally produce two SrAp cells. This model could not be rejected given our observed data. Unlike the disease site, similar analysis of C-to-G mutations at a control nucleotide site in one testis pair failed to find any foci with high mutation frequencies. The rejection of the hot spot model and lack of rejection of a selection model for the C755G mutation, along with other data, provides strong support for the proposal that positive selection in the testis can act to increase the frequency of premeiotic germ cells carrying a mutation deleterious to an offspring, thereby unfavorably altering the mutational load in humans. Studying the anatomical distribution of germline mutations can provide new insights into genetic disease and evolutionary change

Evaluating the Evidence for Transmission Distortion in Human Pedigrees

Children of a heterozygous parent are expected to carry either allele with equal probability. Exceptions can occur, however, due to meiotic drive, competition among gametes, or viability selection, which we collectively term “transmission distortion” (TD). Although there are several well-characterized examples of these phenomena, their existence in humans remains unknown. We therefore performed a genome-wide scan for TD by applying the transmission disequilibrium test (TDT) genome-wide to three large sets of human pedigrees of European descent: the Framingham Heart Study (FHS), a founder population of European origin (HUTT), and a subset of the Autism Genetic Resource Exchange (AGRE). Genotyping error is an important confounder in this type of analysis. In FHS and HUTT, despite extensive quality control, we did not find sufficient evidence to exclude genotyping error in the strongest signals. In AGRE, however, many signals extended across multiple SNPs, a pattern highly unlikely to arise from genotyping error. We identified several candidate regions in this data set, notably a locus in 10q26.13 displaying a genome-wide significant TDT in combined female and male transmissions and a signature of recent positive selection, as well as a paternal TD signal in 6p21.1, the same region in which a significant TD signal was previously observed in 30 European males. Neither region replicated in FHS, however, and the paternal signal was not visible in sperm competition assays or as allelic imbalance in sperm. In maternal transmissions, we detected no strong signals near centromeres or telomeres, the regions predicted to be most susceptible to female-specific meiotic drive, but we found a significant enrichment of top signals among genes involved in cell junctions. These results illustrate both the potential benefits and the challenges of using the TDT to study transmission distortion and provide candidates for investigation in future studies

Massively Parallel Haplotyping on Microscopic Beads for the High-Throughput Phase Analysis of Single Molecules

In spite of the many advances in haplotyping methods, it is still very difficult to characterize rare haplotypes in tissues and different environmental samples or to accurately assess the haplotype diversity in large mixtures. This would require a haplotyping method capable of analyzing the phase of single molecules with an unprecedented throughput. Here we describe such a haplotyping method capable of analyzing in parallel hundreds of thousands single molecules in one experiment. In this method, multiple PCR reactions amplify different polymorphic regions of a single DNA molecule on a magnetic bead compartmentalized in an emulsion drop. The allelic states of the amplified polymorphisms are identified with fluorescently labeled probes that are then decoded from images taken of the arrayed beads by a microscope. This method can evaluate the phase of up to 3 polymorphisms separated by up to 5 kilobases in hundreds of thousands single molecules. We tested the sensitivity of the method by measuring the number of mutant haplotypes synthesized by four different commercially available enzymes: Phusion, Platinum Taq, Titanium Taq, and Phire. The digital nature of the method makes it highly sensitive to detecting haplotype ratios of less than 1∶10,000. We also accurately quantified chimera formation during the exponential phase of PCR by different DNA polymerases

Crossovers are associated with mutation and biased gene conversion at recombination hotspots

Meiosis is a potentially important source of germline mutations, as sites of meiotic recombination experience recurrent double-strand breaks (DSBs). However, evidence for a local mutagenic effect of recombination from population sequence data has been equivocal, likely because mutation is only one of several forces shaping sequence variation. By sequencing large numbers of single crossover molecules obtained from human sperm for two recombination hotspots, we find direct evidence that recombination is mutagenic: Crossovers carry more de novo mutations than nonrecombinant DNA molecules analyzed for the same donors and hotspots. The observed mutations were primarily CG to TA transitions, with a higher frequency of transitions at CpG than non-CpGs sites. This enrichment of mutations at CpG sites at hotspots could predominate in methylated regions involving frequent single-stranded DNA processing as part of DSB repair. In addition, our data set provides evidence that GC alleles are preferentially transmitted during crossing over, opposing mutation, and shows that GC-biased gene conversion (gBGC) predominates over mutation in the sequence evolution of hotspots. These findings are consistent with the idea that gBGC could be an adaptation to counteract the mutational load of recombination

BEH on different template lengths and SNP positions.

<p>The left panel is a schematic of the location of the amplified SNP region relative to the template. The right panel represents the sum of the fluorescence intensity obtained for both alleles for the possible haplotypes. The error bars are the standard deviation of the fluorescence intensity. Shown above each pair of haplotypes is the observed allelic ratio of the captured haplotypes. A. Assay of the same two SNPs (422 bp apart) on templates of different lengths. B. Assay of different SNPs pairs 400, 1200, and 2200 base pairs apart on the same 2733 bp template.</p