Power and limits of selection genome scans on temporal data from a selfing population

Tracking genetic changes of populations through time allows a more direct study of the evolutionary processes acting on the population than a single contemporary sample. Several statistical methods have been developed to characterize the demography and selection from temporal population genetic data. However, these methods are usually developed under the assumption of outcrossing reproduction and might not be applicable when there is substantial selfing in the population. Here, we focus on a method to detect loci under selection based on a genome scan of temporal differentiation, adapting it to the particularities of selfing populations. Selfing reduces the effective recombination rate and can extend hitch-hiking effects to the whole genome, erasing any local signal of selection on a genome scan. Therefore, selfing is expected to reduce the power of the test. By means of simulations, we evaluate the performance of the method under scenarios of adaptation from new mutations or standing variation at different rates of selfing. We find that the detection of loci under selection in predominantly selfing populations remains challenging even with the adapted method. Still, selective sweeps from standing variation on predominantly selfing populations can leave some signal of selection around the selected site thanks to historical recombination before the sweep. Under this scenario, ancestral advantageous alleles at low frequency leave the strongest local signal, while new advantageous mutations leave no local footprint of the sweep.

Historical recombination variability contributes to deciphering the genetic basis of phenotypic traits

Recombination is a main source of genetic variability. However, the potential role of the variation generated by recombination in phenotypic traits, including diseases, remains unexplored as there is currently no method to infer chromosomal subpopulations based on recombination patterns differences. We developed recombClust, a method that uses SNP-phased data to detect differences in historic recombination in a chromosome population. We validated our method by performing simulations and by using real data to accurately predict the alleles of well known recombination modifiers, including common inversions in Drosophila melanogaster and human, and the chromosomes under selective pressure at the lactase locus in humans. We then applied recombClust to the complex human 1q21.1 region, where nonallelic homologous recombination produces deleterious phenotypes. We discovered and validated the presence of two different recombination histories in these regions that significantly associated with the differential expression of ANKRD35 in whole blood and that were in high linkage with variants previously associated with hypertension. By detecting differences in historic recombination, our method opens a way to assess the influence of recombination variation in phenotypic traits.

LDJump: Estimating Variable Recombination Rates from Population Genetic Data

As recombination plays an important role in evolution, its estimation, as well as, the identification of hotspot positions is of considerable interest. We propose a novel approach for estimating historical recombination along a chromosome that involves a sequential multiscale change point estimator. Our method also permits to take demography into account. It uses a composite likelihood estimate and other summary statistics within a regression model fitted on suitable scenarios. Our proposed method is accurate, computationally fast, and provides a parsimonious solution by ensuring a type I error control against too many changes in the recombination rate. An application to human genome data suggests a good congruence between our estimated and experimentally identified hotspots. Our method is implemented in the R-package LDJump, which is freely available from https://github.com/PhHermann/LDJump.

A Coalescent Model of Recombination Hotspots

Recent experimental findings suggest that the assumption of a homogeneous recombination rate along the human genome is too naive. These findings point to block-structured recombination rates; certain regions (called hotspots) are more prone than other regions to recombination. In this report a coalescent model incorporating hotspot or block-structured recombination is developed and investigated analytically as well as by simulation. Our main results can be summarized as follows: (1) The expected number of recombination events is much lower in a model with pure hotspot recombination than in a model with pure homogeneous recombination, (2) hotspots give rise to large variation in recombination rates along the genome as well as in the number of historical recombination events, and (3) the size of a (nonrecombining) block in the hotspot model is likely to be overestimated grossly when estimated from SNP data. The results are discussed with reference to the current debate about block-structured recombination and, in addition, the results are compared to genome-wide variation in recombination rates. A number of new analytical results about the model are derived.

Simultaneous Detection and Fine Mapping of Quantitative Trait Loci in Mice Using Heterogeneous Stocks

We describe a method to simultaneously detect and fine map quantitative trait loci (QTL) that is especially suited to the mapping of modifier loci in mouse mutant models. The method exploits the high level of historical recombination present in a heterogeneous stock (HS), an outbred population of mice derived from known founder strains. The experimental design is an F2 cross between the HS and a genetically distinct line, such as one carrying a knockout or transgene. QTL detection is performed by a standard genome scan with ~100 markers and fine mapping by typing the same animals using densely spaced markers over those candidate regions detected by the scan. The analysis uses an extension of the dynamic-programming technique employed previously to fine map QTL in HS mice. We show by simulation that a QTL accounting for 5% of the total variance can be detected and fine mapped with >50% probability to within 3 cM by genotyping ~1500 animals.