Revisiting the false positive rate in detecting recent positive selection

Power of neutrality tests for detecting natural selection

Detection of natural selection is one of the main interests in population genetics. Thus, many tests have been developed for detecting natural selection using genomic data. Although it is recognized that the utility of tests depends on several evolutionary factors, such as the timing of selection, strength of selection, frequency of selected alleles, demographic events, and initial frequency of selected allele when selection started acting (softness of selection), the relationships between such evolutionary factors and the power of tests are not yet entirely clear. In this study, we investigated the power of 4 tests: Tajiama's D, Fay and Wu's H, relative extended haplotype homozygosity (rEHH), and integrated haplotype score (iHS), under ranges of evolutionary parameters and demographic models to quantitatively expand the understanding of approaches for detecting selection. The results show that each test detects selection within a limited parameter range, and there are still wide ranges of parameters for which none of these tests work effectively. In addition, the parameter space in which each test shows the highest power overlaps the empirical results of previous research. These results indicate that our present perspective of adaptation is limited to only a part of actual adaptation.

Bray-Curtis (AFD) differentiation in molecular ecology: Forecasting, an adjustment ( A A), and comparative performance in selection detection

Geographic genetic differentiation measures are used for purposes such as assessing genetic diversity and connectivity, and searching for signals of selection. Confirmation by unrelated measures can minimize false positives. A popular differentiation measure, Bray‐Curtis, has been used increasingly in molecular ecology, renamed AFD (hereafter called BCAFD). Critically, BCAFD is expected to be partially independent of the commonly used Hill “Q‐profile” measures. BCAFD needs scrutiny for potential biases, by examining limits on its value, and comparing simulations against expectations. BCAFD has two dependencies on within‐population (alpha) variation, undesirable for a between‐population (beta) measure. The first dependency is derived from similarity to GST and FST. The second dependency is that BCAFD cannot be larger than the highest allele proportion in either location (alpha variation), which can be overcome by data‐filtering or by a modified statistic ^AA or “Adjusted AFD”. The first dependency does not forestall applications such as assessing connectivity or selection, if we know the measure's null behavior under selective neutrality with specified conditions—which is shown in this article for ^AA, for equilibrium, and nonequilibrium, for the commonly used data type of single‐nucleotide‐polymorphisms (SNPs) in two locations. Thus, ^AA can be used in tandem with mathematically contrasting differentiation measures, with the aim of reducing false inferences. For detecting adaptive loci, the relative performance of ^AA and other measures was evaluated, showing that it is best to use two mathematically different measures simultaneously, and that ^AA is in one of the best such pairwise criteria. For any application, using ^AA, rather than BCAFD, avoids the counterintuitive limitation by maximum allele proportion within localities.

Human-modified canids in human-modified landscapes: The evolutionary consequences of hybridization for grey wolves and free-ranging domestic dogs

Introgressive hybridization between domestic animals and their wild relatives is an indirect form of human-induced evolution, altering gene pools and phenotypic traits of wild and domestic populations. Although this process is well documented in many taxa, its evolutionary consequences are poorly understood. In this study, we assess introgression patterns in admixed populations of Eurasian wolves and free-ranging domestic dogs (FRDs), identifying chromosomal regions with significantly overrepresented hybrid ancestry and assessing whether genes located within these regions show signatures of selection. Although the dog admixture proportion in West Eurasian wolves (2.7%) was greater than the wolf admixture proportion in FRDs (0.75%), the number and average length of chromosomal blocks showing significant overrepresentation of hybrid ancestry were smaller in wolves than FRDs. In wolves, 6% of genes located within these blocks showed signatures of positive selection compared to 23% in FRDs. We found that introgression from wolves may provide a considerable adaptive advantage to FRDs, counterbalancing some of the negative effects of domestication, which can include reduced genetic diversity and excessive tameness. In wolves, introgression from FRDs is mostly driven by drift, with a small number of positively selected genes associated with brain function and behaviour. The predominance of drift may be the consequence of small effective size of wolf populations, which reduces efficiency of selection for weakly advantageous or against weakly disadvantageous introgressed variants. Small wolf population sizes result largely from human-induced habitat loss and hunting, thus linking introgression rates to anthropogenic processes. Our results imply that maintenance of large population sizes should be an important element of wolf management strategies aimed at reducing introgression rates of dog-derived variants.

Is there adaptation in the human genome for taste perception and phase I biotransformation?

Background

During the modern human expansion, new environmental pressures may have driven adaptation, especially in genes related to the perception of ingested substances and their detoxification. Consequently, positive (adaptive) selection may have occurred in genes related to taste, and in those related to the CYP450 system due to its role in biotransformation of potentially toxic compounds. A total of 91 genes (taste receptors and CYP450 superfamily) have been studied using Hierarchical Boosting, a powerful combination of different selection tests, to detect signatures of recent positive selection in three continental human populations: Northern Europeans (CEU), East Asians (CHB) and Africans (YRI). Analyses have been refined with selection analyses of the 26 populations of 1000 Genomes Project Phase 3.

Results

Genes related to taste perception have not been positively selected in the three continental human populations. This finding suggests that, contrary to results of previous studies, different allele frequencies among populations in genes such as TAS2R38 and TAS2R16 are not due to positive selection but to genetic drift. CYP1 and CYP2 genes, also previously considered to be under positive selection, did not show signatures of selective sweeps. However, three genes belonging to the CYP450 system have been identified by the Hierarchical Boosting as positively selected: CYP3A4 and CYP3A43 in CEU, and CYP27A1 in CHB.

Conclusions

No main adaptive differences are found in known taste receptor genes among the three continental human populations studied. However, there are important genetic adaptations in the cytochrome P450 system related to the Out of Africa expansion of modern humans. We confirmed that CYP3A4 and CYP3A43 are under selection in CEU, and we report for the first time CYP27A1 to be under positive selection in CHB.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1366-7) contains supplementary material, which is available to authorized users.

Genome-Wide Detection of Genes Under Positive Selection in Worldwide Populations of the Barley Scald Pathogen

Coevolution between hosts and pathogens generates strong selection pressures to maintain resistance and infectivity, respectively. Genomes of plant pathogens often encode major effect loci for the ability to successfully infect specific host genotypes. Hence, spatial heterogeneity in host genotypes coupled with abiotic factors could lead to locally adapted pathogen populations. However, the genetic basis of local adaptation is poorly understood. Rhynchosporium commune, the pathogen causing barley scald disease, interacts at least partially in a gene-for-gene manner with its host. We analyzed global field populations of 125 R. commune isolates to identify candidate genes for local adaptation. Whole genome sequencing data showed that the pathogen is subdivided into three genetic clusters associated with distinct geographic and climatic regions. Using haplotype-based selection scans applied independently to each genetic cluster, we found strong evidence for selective sweeps throughout the genome. Comparisons of loci under selection among clusters revealed little overlap, suggesting that ecological differences associated with each cluster led to variable selection regimes. The strongest signals of selection were found predominantly in the two clusters composed of isolates from Central Europe and Ethiopia. The strongest selective sweep regions encoded protein functions related to biotic and abiotic stress responses. Selective sweep regions were enriched in genes encoding functions in cellular localization, protein transport activity, and DNA damage responses. In contrast to the prevailing view that a small number of gene-for-gene interactions govern plant pathogen evolution, our analyses suggest that the evolutionary trajectory is largely determined by spatially heterogeneous biotic and abiotic selection pressures.

Detecting Recent Positive Selection with a Single Locus Test Bipartitioning the Coalescent Tree

Many population genomic studies have been conducted in the past to search for traces of recent events of positive selection. These traces, however, can be obscured by temporal variation of population size or other demographic factors. To reduce the confounding impact of demography, the coalescent tree topology has been used as an additional source of information for detecting recent positive selection in a population or a species. Based on the branching pattern at the root, we partition the hypothetical coalescent tree, inferred from a sequence sample, into two subtrees. The reasoning is that positive selection could impose a strong impact on branch length in one of the two subtrees while demography has the same effect on average on both subtrees. Thus, positive selection should be detectable by comparing statistics calculated for the two subtrees. Simulations demonstrate that the proposed test based on these principles has high power to detect recent positive selection even when DNA polymorphism data from only one locus is available, and that it is robust to the confounding effect of demography. One feature is that all components in the summary statistics ([Formula: see text]) can be computed analytically. Moreover, misinference of derived and ancestral alleles is seen to have only a limited effect on the test, and it therefore avoids a notorious problem when searching for traces of recent positive selection.