Human genomic regions with exceptionally high levels of population differentiation identified from 911 whole-genome sequences

Background

Population differentiation has proved to be effective for identifying loci under geographically localized positive selection, and has the potential to identify loci subject to balancing selection. We have previously investigated the pattern of genetic differentiation among human populations at 36.8 million genomic variants to identify sites in the genome showing high frequency differences. Here, we extend this dataset to include additional variants, survey sites with low levels of differentiation, and evaluate the extent to which highly differentiated sites are likely to result from selective or other processes.

Results

We demonstrate that while sites with low differentiation represent sampling effects rather than balancing selection, sites showing extremely high population differentiation are enriched for positive selection events and that one half may be the result of classic selective sweeps. Among these, we rediscover known examples, where we actually identify the established functional SNP, and discover novel examples including the genes ABCA12, CALD1 and ZNF804, which we speculate may be linked to adaptations in skin, calcium metabolism and defense, respectively.

Conclusions

We identify known and many novel candidate regions for geographically restricted positive selection, and suggest several directions for further research.     

A Selective Sweep on a Deleterious Mutation in CPT1A in Arctic Populations

Arctic populations live in an environment characterized by extreme cold and the absence of plant foods for much of the year and are likely to have undergone genetic adaptations to these environmental conditions in the time they have been living there. Genome-wide selection scans based on genotype data from native Siberians have previously highlighted a 3 Mb chromosome 11 region containing 79 protein-coding genes as the strongest candidates for positive selection in Northeast Siberians. However, it was not possible to determine which of the genes might be driving the selection signal. Here, using whole-genome high-coverage sequence data, we identified the most likely causative variant as a nonsynonymous G>A transition (rs80356779; c.1436C>T [p.Pro479Leu] on the reverse strand) in CPT1A, a key regulator of mitochondrial long-chain fatty-acid oxidation. Remarkably, the derived allele is associated with hypoketotic hypoglycemia and high infant mortality yet occurs at high frequency in Canadian and Greenland Inuits and was also found at 68% frequency in our Northeast Siberian sample. We provide evidence of one of the strongest selective sweeps reported in humans; this sweep has driven this variant to high frequency in circum-Arctic populations within the last 6–23 ka despite associated deleterious consequences, possibly as a result of the selective advantage it originally provided to either a high-fat diet or a cold environment.     

Where west meets east: the complex mtDNA landscape of the southwest and Central Asian corridor

The southwestern and Central Asian corridor has played a pivotal role in the history of humankind, witnessing numerous waves of migration of different peoples at different times. To evaluate the effects of these population movements on the current genetic landscape of the Iranian plateau, the Indus Valley, and Central Asia, we have analyzed 910 mitochondrial DNAs (mtDNAs) from 23 populations of the region. This study has allowed a refinement of the phylogenetic relationships of some lineages and the identification of new haplogroups in the southwestern and Central Asian mtDNA tree. Both lineage geographical distribution and spatial analysis of molecular variance showed that populations located west of the Indus Valley mainly harbor mtDNAs of western Eurasian origin, whereas those inhabiting the Indo-Gangetic region and Central Asia present substantial proportions of lineages that can be allocated to three different genetic components of western Eurasian, eastern Eurasian, and south Asian origin. In addition to the overall composite picture of lineage clusters of different origin, we observed a number of deep-rooting lineages, whose relative clustering and coalescent ages suggest an autochthonous origin in the southwestern Asian corridor during the Pleistocene. The comparison with Y-chromosome data revealed a highly complex genetic and demographic history of the region, which includes sexually asymmetrical mating patterns, founder effects, and female-specific traces of the East African slave trade.     

Ethiopian genetic diversity reveals linguistic stratification and complex influences on the Ethiopian gene pool

Humans and their ancestors have traversed the Ethiopian landscape for millions of years, and present-day Ethiopians show great cultural, linguistic, and historical diversity, which makes them essential for understanding African variability and human origins. We genotyped 235 individuals from ten Ethiopian and two neighboring (South Sudanese and Somali) populations on an Illumina Omni 1M chip. Genotypes were compared with published data from several African and non-African populations. Principal-component and STRUCTURE-like analyses confirmed substantial genetic diversity both within and between populations, and revealed a match between genetic data and linguistic affiliation. Using comparisons with African and non-African reference samples in 40-SNP genomic windows, we identified "African" and "non-African" haplotypic components for each Ethiopian individual. The non-African component, which includes the SLC24A5 allele associated with light skin pigmentation in Europeans, may represent gene flow into Africa, which we estimate to have occurred ~3 thousand years ago (kya). The non-African component was found to be more similar to populations inhabiting the Levant rather than the Arabian Peninsula, but the principal route for the expansion out of Africa ~60 kya remains unresolved. Linkage-disequilibrium decay with genomic distance was less rapid in both the whole genome and the African component than in southern African samples, suggesting a less ancient history for Ethiopian populations.     

Identification of neutral alleles at pollen sterility gene loci of cultivated rice (<Emphasis Type="Italic">Oryza sativa</Emphasis> L.) from wild rice (<Emphasis Type="Italic">O. rufipogon</Emphasis> Griff.)

Oryza rufipogon Griff. is the ancestor of Asian cultivated rice (O. sativa L.) and possesses valuable genes for rice breeding. Pollen abortion is one of the major causes of indica–japonica hybrid sterility in rice and it happens due to allelic interaction at the pollen sterility gene loci. A total of six loci (Sa, Sb, Sc, Sd, Se, and Sf) have been found to be associated with F₁ pollen sterility between indica and japonica rice, and five of them (all except Sf) have been mapped. Neutral alleles (S ⁿ ) at each locus have the potential to overcome the pollen sterility associated with the respective locus. Therefore, exploitation and utilization of neutral alleles have significant importance in overcoming indica–japonica hybrid sterility. In this study, an accession (IRW28) of O. rufipogon, native to China, was selected as paternal to cross with typical japonica (Taichung 65) and indica (Guanglu’ai 4) tester lines, and two F₂ populations were developed. The simple sequence repeat markers tightly linked to five pollen sterility loci were applied for genotyping the F₂ populations. Chi-squared tests were applied to examine the normal segregation/distortion at each locus. The expected and observed pollen sterility for each locus were estimated. As a result, the genotypes at five pollen sterility gene loci for IRW28 were identified as: Sa ⁱ⁻¹/Sa ⁱ⁻¹, Sb ⁿ /Sb ⁿ , Sc ⁱ⁻²/Sc ⁱ⁻², Sd ⁿ /Sd ⁿ and Se ⁿ /Se ⁿ . Our results suggest that IRW28 (O. rufipogon) has the neutral alleles for pollen fertility at the Sb, Sd and Se loci, and these alleles have a good affinity with indica and japonica rice. These neutral alleles provide valuable gene resources to overcome the inter-subspecific hybrid pollen sterility in rice.     

Investigating the Link Between Imipenem Resistance and Biofilm Formation by Pseudomonas aeruginosa

Pseudomonas aeruginosa, a ubiquitous environmental organism, is a difficult-to-treat opportunistic pathogen due to its broad-spectrum antibiotic resistance and its ability to form biofilms. In this study, we investigate the link between resistance to a clinically important antibiotic, imipenem, and biofilm formation. First, we observed that the laboratory strain P. aeruginosa PAO1 carrying a mutation in the oprD gene, which confers resistance to imipenem, showed a modest reduction in biofilm formation. We also observed an inverse relationship between imipenem resistance and biofilm formation for imipenem-resistant strains selected in vitro, as well as for clinical isolates. We identified two clinical isolates of P. aeruginosa from the sputum of cystic fibrosis patients that formed robust biofilms, but were sensitive to imipenem (MIC?≤?2 μg/ml). To test the hypothesis that there is a general link between imipenem resistance and biofilm formation, we performed transposon mutagenesis of these two clinical strains to identify mutants defective in biofilm formation, and then tested these mutants for imipenem resistance. Analysis of the transposon mutants revealed a role for previously described biofilm factors in these clinical isolates of P. aeruginosa, including mutations in the pilY1, pilX, pilW, algC, and pslI genes, but none of the biofilm-deficient mutants became imipenem resistant (MIC?≥?8 μg/ml), arguing against a general link between biofilm formation and resistance to imipenem. Thus, assessing biofilm formation capabilities of environmental isolates is unlikely to serve as a good predictor of imipenem resistance. We also discuss our findings in light of the limited literature addressing planktonic antibiotic resistance factors that impact biofilm formation.     

Independent introduction of two lactase-persistence alleles into human populations reflects different history of adaptation to milk culture

The T₋₁₃₉₁₀ variant located in the enhancer element of the lactase (LCT) gene correlates perfectly with lactase persistence (LP) in Eurasian populations whereas the variant is almost nonexistent among Sub-Saharan African populations, showing high prevalence of LP. Here, we report identification of two new mutations among Saudis, also known for the high prevalence of LP. We confirmed the absence of the European T₋₁₃₉₁₀ and established two new mutations found as a compound allele: T/G₋₁₃₉₁₅ within the −13910 enhancer region and a synonymous SNP in the exon 17 of the MCM6 gene T/C₋₃₇₁₂, −3712 bp from the LCT gene. The compound allele is driven to a high prevalence among Middle East population(s). Our functional analyses in vitro showed that both SNPs of the compound allele, located 10 kb apart, are required for the enhancer effect, most probably mediated through the binding of the hepatic nuclear factor 1 α (HNF1α). High selection coefficient (s) ~0.04 for LP phenotype was found for both T₋₁₃₉₁₀ and the compound allele. The European T₋₁₃₉₁₀ and the earlier identified East African G₋₁₃₉₀₇ LP allele share the same ancestral background and most likely the same history, probably related to the same cattle domestication event. In contrast, the compound Arab allele shows a different, highly divergent ancestral haplotype, suggesting that these two major global LP alleles have arisen independently, the latter perhaps in response to camel milk consumption. These results support the convergent evolution of the LP in diverse populations, most probably reflecting different histories of adaptation to milk culture.