ABSTRACT
Following introgression, Neanderthal DNA was initially purged from non-African genomes, but the evolutionary fate of remaining introgressed DNA has not been explored yet. To fill this gap, we analyzed 30,780 admixed genomes with African-like ancestry from the All of Us research program, in which Neanderthal alleles encountered novel genetic backgrounds during the last 15 generations. Observed amounts of Neanderthal DNA approximately match expectations based on ancestry proportions, suggesting neutral evolution. Nevertheless, we identified genomic regions that have significantly less or more Neanderthal ancestry than expected and are associated with spermatogenesis, innate immunity, and other biological processes. We also identified three novel introgression desert-like regions in recently admixed genomes, whose genetic features are compatible with hybrid incompatibilities and intrinsic negative selection. Overall, we find that much of the remaining Neanderthal DNA in human genomes is not under strong selection, and complex evolutionary dynamics have shaped introgression landscapes in our species.
ABSTRACT
Lengths of telomeres vary by an order of magnitude across mammalian species. Similarly, age- and sex-standardized telomere lengths differ by up to 1â kb (14%) across human populations. How to explain these differences? Telomeres play a central role in senescence and aging, and genes that affect telomere length are likely under weak selection (i.e. telomere length is a trait that is subject to nearly neutral evolution). Importantly, natural selection is more effective in large populations than in small populations. Here, we propose that observed differences in telomere length across species and populations are largely due to differences in effective population sizes. In this perspective, we present preliminary evolutionary genetic evidence supporting this hypothesis and highlight the need for more data.
Subject(s)
Population Density , Selection, Genetic , Telomere , Humans , Telomere/genetics , Animals , Evolution, Molecular , Telomere Homeostasis , Biological Evolution , Aging/geneticsABSTRACT
Androgenetic alopecia is a highly heritable trait. However, much of our understanding about the genetics of male pattern baldness comes from individuals of European descent. Here, we examined a novel dataset comprising 2,136 men from Ghana, Nigeria, Senegal, and South Africa that were genotyped using a custom array. We first tested how genetic predictions of baldness generalize from Europe to Africa, finding that polygenic scores from European GWAS yielded AUC statistics that ranged from 0.513 to 0.546, indicating that genetic predictions of baldness in African populations performed notably worse than in European populations. Subsequently, we conducted the first African GWAS of androgenetic alopecia, focusing on self-reported baldness patterns at age 45. After correcting for present age, population structure, and study site, we identified 266 moderately significant associations, 51 of which were independent (p-value < 10-5, r2 < 0.2). Most baldness associations were autosomal, and the X chromosomes does not appear to have a large impact on baldness in African men. Finally, we examined the evolutionary causes of continental differences in genetic architecture. Although Neanderthal alleles have previously been associated with skin and hair phenotypes, we did not find evidence that European-ascertained baldness hits were enriched for signatures of ancient introgression. Most loci that are associated with androgenetic alopecia are evolving neutrally. However, multiple baldness-associated SNPs near the EDA2R and AR genes have large allele frequency differences between continents. Collectively, our findings illustrate how evolutionary history contributes to the limited portability of genetic predictions across ancestries.
ABSTRACT
We conducted a multi-ancestry genome-wide association study of prostate-specific antigen (PSA) levels in 296,754 men (211,342 European ancestry; 58,236 African ancestry; 23,546 Hispanic/Latino; 3,630 Asian ancestry; 96.5% of participants were from the Million Veteran Program). We identified 318 independent genome-wide significant (p≤5e-8) variants, 184 of which were novel. Most demonstrated evidence of replication in an independent cohort (n=95,768). Meta-analyzing discovery and replication (n=392,522) identified 447 variants, of which a further 111 were novel. Out-of-sample variance in PSA explained by our new polygenic risk score reached 16.9% (95% CI=16.1%-17.8%) in European ancestry, 9.5% (95% CI=7.0%-12.2%) in African ancestry, 18.6% (95% CI=15.8%-21.4%) in Hispanic/Latino, and 15.3% (95% CI=12.7%-18.1%) in Asian ancestry, and lower for higher age. Our study highlights how including proportionally more participants from underrepresented populations improves genetic prediction of PSA levels, with potential to personalize prostate cancer screening.
ABSTRACT
BACKGROUND: Genetic factors play an important role in prostate cancer (PCa) susceptibility. OBJECTIVE: To discover common genetic variants contributing to the risk of PCa in men of African ancestry. DESIGN, SETTING, AND PARTICIPANTS: We conducted a meta-analysis of ten genome-wide association studies consisting of 19378 cases and 61620 controls of African ancestry. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: Common genotyped and imputed variants were tested for their association with PCa risk. Novel susceptibility loci were identified and incorporated into a multiancestry polygenic risk score (PRS). The PRS was evaluated for associations with PCa risk and disease aggressiveness. RESULTS AND LIMITATIONS: Nine novel susceptibility loci for PCa were identified, of which seven were only found or substantially more common in men of African ancestry, including an African-specific stop-gain variant in the prostate-specific gene anoctamin 7 (ANO7). A multiancestry PRS of 278 risk variants conferred strong associations with PCa risk in African ancestry studies (odds ratios [ORs] >3 and >5 for men in the top PRS decile and percentile, respectively). More importantly, compared with men in the 40-60% PRS category, men in the top PRS decile had a significantly higher risk of aggressive PCa (OR = 1.23, 95% confidence interval = 1.10-1.38, p = 4.4 × 10-4). CONCLUSIONS: This study demonstrates the importance of large-scale genetic studies in men of African ancestry for a better understanding of PCa susceptibility in this high-risk population and suggests a potential clinical utility of PRS in differentiating between the risks of developing aggressive and nonaggressive disease in men of African ancestry. PATIENT SUMMARY: In this large genetic study in men of African ancestry, we discovered nine novel prostate cancer (PCa) risk variants. We also showed that a multiancestry polygenic risk score was effective in stratifying PCa risk, and was able to differentiate risk of aggressive and nonaggressive disease.
Subject(s)
Genetic Predisposition to Disease , Prostatic Neoplasms , Male , Humans , Genome-Wide Association Study , Prostatic Neoplasms/genetics , Prostatic Neoplasms/epidemiology , Risk Factors , Black People/geneticsABSTRACT
As the ancestral homeland of our species, Africa contains elevated levels of genetic diversity and substantial population structure. Importantly, African genomes are heterogeneous: They contain mixtures of multiple ancestries, each of which have experienced different evolutionary histories. In this review, we view population genetics through the lens of admixture, highlighting how multiple demographic events have shaped African genomes. Each of these historical vignettes paints a recurring picture of population divergence followed by secondary contact. First, we give a brief overview of genetic variation in Africa and examine deep population structure within Africa, including the evidence of ancient introgression from archaic "ghost" populations. Second, we describe the genetic legacies of admixture events that have occurred during the past 10,000 years. This includes gene flow between different click-speaking Khoe-San populations, the stepwise spread of pastoralism from eastern to southern Africa, multiple migrations of Bantu speakers across the continent, as well as admixture from the Middle East and Europe into the Sahel region and North Africa. Furthermore, the genomic signatures of more recent admixture can be found in the Cape Peninsula and throughout the African diaspora. Third, we highlight how natural selection has shaped patterns of genetic variation across the continent, noting that gene flow provides a potent source of adaptive variation and that selective pressures vary across Africa. Finally, we explore the biomedical implications of population structure in Africa on health and disease and call for more ethically conducted studies of genetic variation in Africa.
Subject(s)
Genetic Variation , Genetics, Population , Africa, Southern , Biological Evolution , GenomeABSTRACT
Sex-biased admixture can be inferred from ancestry-specific proportions of X chromosome and autosomes. In a paper published in the American Journal of Human Genetics, Micheletti et al.1 used this approach to quantify male and female contributions following the transatlantic slave trade. Using a large dataset from 23andMe, they concluded that African and European contributions to gene pools in the Americas were much more sex biased than previously thought. We show that the reported extreme sex-specific contributions can be attributed to unassigned genetic ancestry as well as the limitations of simple models of sex-biased admixture. Unassigned ancestry proportions in the study by Micheletti et al. ranged from â¼1% to 21%, depending on the type of chromosome and geographic region. A sensitivity analysis illustrates how this unassigned ancestry can create false patterns of sex bias and that mathematical models are highly sensitive to slight sampling errors when inferring mean ancestry proportions, making confidence intervals necessary. Thus, unassigned ancestry and the sensitivity of the models effectively prohibit the interpretation of estimated sex biases for many geographic regions in Micheletti et al. Furthermore, Micheletti et al. assumed models of a single admixture event. Using simulations, we find that violations of demographic assumptions, such as subsequent gene flow and/or sex-specific assortative mating, may have confounded the analyses of Micheletti et al., but unassigned ancestry was likely the more important confounding factor. Our findings underscore the importance of using complete ancestry information, sufficiently large sample sizes, and appropriate models when inferring sex-biased patterns of demography. This Matters Arising paper is in response to Micheletti et al.,1 published in American Journal of Human Genetics. See also the response by Micheletti et al.,2 published in this issue.
Subject(s)
Genetics, Population , Sexism , Female , Humans , Male , Chromosomes , Gene Flow , Africa , EuropeABSTRACT
Many lizard species face extinction due to worldwide climate change. The Guatemalan Beaded Lizard, Heloderma charlesbogerti, is a member of the Family Helodermatidae that may be particularly imperiled; fewer than 600 mature individuals are believed to persist in the wild. In addition, H. charlesbogerti lizards are phenotypically remarkable. They are large in size, charismatically patterned, and possess a venomous bite. Here, we report the draft genome of the Guatemalan Beaded Lizard using DNA from a wild-caught individual. The assembled genome totals 2.31 Gb in length, similar in size to the genomes of related species. Single-copy orthologs were used to produce a novel molecular phylogeny, revealing that the Guatemalan Beaded Lizard falls into a clade with the Asian Glass Lizard (Anguidae) and in close association with the Komodo Dragon (Varanidae) and the Chinese Crocodile Lizard (Shinisauridae). In addition, we identified 31,411 protein-coding genes within the genome. Of the genes identified, we found 504 that evolved with a differential constraint on the branch leading to the Guatemalan Beaded Lizard. Lastly, we identified a decline in the effective population size of the Guatemalan Beaded Lizard approximately 400,000 years ago, followed by a stabilization before starting to dwindle again 60,000 years ago. The results presented here provide important information regarding a highly endangered, venomous reptile that can be used in future conservation, functional genetic, and phylogenetic analyses.
Subject(s)
Lizards , Humans , Animals , Lizards/genetics , Phylogeny , Population Density , Venoms/genetics , GenomeABSTRACT
BACKGROUND: Genome-wide association studies do not always replicate well across populations, limiting the generalizability of polygenic risk scores (PRS). Despite higher incidence and mortality rates of prostate cancer in men of African descent, much of what is known about cancer genetics comes from populations of European descent. To understand how well genetic predictions perform in different populations, we evaluated test characteristics of PRS from three previous studies using data from the UK Biobank and a novel dataset of 1298 prostate cancer cases and 1333 controls from Ghana, Nigeria, Senegal, and South Africa. RESULTS: Allele frequency differences cause predicted risks of prostate cancer to vary across populations. However, natural selection is not the primary driver of these differences. Comparing continental datasets, we find that polygenic predictions of case vs. control status are more effective for European individuals (AUC 0.608-0.707, OR 2.37-5.71) than for African individuals (AUC 0.502-0.585, OR 0.95-2.01). Furthermore, PRS that leverage information from African Americans yield modest AUC and odds ratio improvements for sub-Saharan African individuals. These improvements were larger for West Africans than for South Africans. Finally, we find that existing PRS are largely unable to predict whether African individuals develop aggressive forms of prostate cancer, as specified by higher tumor stages or Gleason scores. CONCLUSIONS: Genetic predictions of prostate cancer perform poorly if the study sample does not match the ancestry of the original GWAS. PRS built from European GWAS may be inadequate for application in non-European populations and perpetuate existing health disparities.
Subject(s)
Genome-Wide Association Study , Prostatic Neoplasms , Africa South of the Sahara/epidemiology , Genetic Predisposition to Disease , Humans , Male , Prostatic Neoplasms/genetics , Risk FactorsABSTRACT
Background: We recently developed a multi-ancestry polygenic risk score (PRS) that effectively stratifies prostate cancer risk across populations. In this study, we validated the performance of the PRS in the multi-ancestry Million Veteran Program and additional independent studies. Methods: Within each ancestry population, the association of PRS with prostate cancer risk was evaluated separately in each case-control study and then combined in a fixed-effects inverse-variance-weighted meta-analysis. We further assessed the effect modification by age and estimated the age-specific absolute risk of prostate cancer for each ancestry population. Results: The PRS was evaluated in 31,925 cases and 490,507 controls, including men from European (22,049 cases, 414,249 controls), African (8794 cases, 55,657 controls), and Hispanic (1082 cases, 20,601 controls) populations. Comparing men in the top decile (90-100% of the PRS) to the average 40-60% PRS category, the prostate cancer odds ratio (OR) was 3.8-fold in European ancestry men (95% CI = 3.62-3.96), 2.8-fold in African ancestry men (95% CI = 2.59-3.03), and 3.2-fold in Hispanic men (95% CI = 2.64-3.92). The PRS did not discriminate risk of aggressive versus nonaggressive prostate cancer. However, the OR diminished with advancing age (European ancestry men in the top decile: ≤55 years, OR = 7.11; 55-60 years, OR = 4.26; >70 years, OR = 2.79). Men in the top PRS decile reached 5% absolute prostate cancer risk ~10 years younger than men in the 40-60% PRS category. Conclusions: Our findings validate the multi-ancestry PRS as an effective prostate cancer risk stratification tool across populations. A clinical study of PRS is warranted to determine whether the PRS could be used for risk-stratified screening and early detection. Funding: This work was supported by the National Cancer Institute at the National Institutes of Health (grant numbers U19 CA214253 to C.A.H., U01 CA257328 to C.A.H., U19 CA148537 to C.A.H., R01 CA165862 to C.A.H., K99 CA246063 to B.F.D, and T32CA229110 to F.C), the Prostate Cancer Foundation (grants 21YOUN11 to B.F.D. and 20CHAS03 to C.A.H.), the Achievement Rewards for College Scientists Foundation Los Angeles Founder Chapter to B.F.D, and the Million Veteran Program-MVP017. This research has been conducted using the UK Biobank Resource under application number 42195. This research is based on data from the Million Veteran Program, Office of Research and Development, and the Veterans Health Administration. This publication does not represent the views of the Department of Veteran Affairs or the United States Government.
Subject(s)
Genome-Wide Association Study , Prostatic Neoplasms , Age Factors , Case-Control Studies , Genetic Predisposition to Disease , Humans , Male , Middle Aged , Multifactorial Inheritance , Prostatic Neoplasms/epidemiology , Prostatic Neoplasms/genetics , Risk Factors , United States/epidemiologyABSTRACT
Hybridization is a common occurrence in natural populations, and introgression is a major source of genetic variation. Despite the evolutionary importance of adaptive introgression, classical population genetics theory does not take into account hybrid fitness effects. Specifically, heterosis (i.e. hybrid vigor) and Dobzhansky-Muller incompatibilities influence the fates of introgressed alleles. Here, we explicitly account for polygenic, unlinked hybrid fitness effects when tracking a rare introgressed marker allele. These hybrid fitness effects quickly decay over time due to repeated backcrossing, enabling a separation-of-timescales approach. Using diffusion and branching process theory in combination with computer simulations, we formalize the intuition behind how hybrid fitness effects affect introgressed alleles. We find that hybrid fitness effects can significantly hinder or boost the fixation probability of introgressed alleles, depending on the relative strength of heterosis and Dobzhansky-Muller incompatibilities effects. We show that the inclusion of a correction factor (α, representing the compounded effects of hybrid fitness effects over time) into classic population genetics theory yields accurate fixation probabilities. Despite having a strong impact on the probability of fixation, hybrid fitness effects only subtly change the distribution of fitness effects of introgressed alleles that reach fixation. Although strong Dobzhansky-Muller incompatibility effects may expedite the loss of introgressed alleles, fixation times are largely unchanged by hybrid fitness effects.
Subject(s)
Genetics, Population , Models, Genetic , Alleles , Hybridization, Genetic , ProbabilityABSTRACT
Large numbers of expression quantitative trait loci (eQTLs) have recently been identified in humans, and many of these regulatory variants have large allele frequency differences between populations. Here, we conducted genome-wide scans of selection to identify adaptive eQTLs (i.e., eQTLs with large population branch statistics). We then tested if tissue pleiotropy affects whether eQTLs are more or less likely to be adaptive and identified tissues that have been key targets of positive selection during the last 100,000 years. Top adaptive eQTL outliers include rs1043809, rs66899053, and rs2814778 (a SNP that is associated with malaria resistance). We found that effect sizes of eQTLs were negatively correlated with population branch statistics and that adaptive eQTLs affect two-thirds as many tissues as do non-adaptive eQTLs. Because the tissue breadth of an eQTL can be viewed as a measure of pleiotropy, these results imply that pleiotropy inhibits adaptation. The proportion of eQTLs that are adaptive varies by tissue, and we found that eQTLs that regulate expression in testis, thyroid, blood, or sun-exposed skin are enriched for signatures of positive selection. By contrast, eQTLs that regulate expression in the cerebrum or female-specific tissues have a relative lack of adaptive outliers. Scans of selections also reveal that many adaptive eQTLs are closely linked to disease-associated loci. Taken together, our results indicate that eQTLs have played an important role in recent human evolution.
ABSTRACT
In this issue of Cancer Research, Emami and colleagues leveraged genetic data from over 200,000 men of European descent to implicate rare alleles that are associated with prostate cancer. However, this study went beyond a simple description of statistical associations between genetic variants and cancer risk. Polygenic risk scores were applied to large cohorts from Kaiser Permanente and the UK Biobank, demonstrating the clinical utility of genetic predictors of disease risk. Furthermore, by placing their results in an evolutionary framework and integrating genetic information with functional data, the authors of this major study were able to bridge the gap between genome-wide association studies and the biological mechanisms underlying prostate cancer risk.See related article by Emami et al., 1695.
Subject(s)
Genome-Wide Association Study , Prostatic Neoplasms , Genetic Predisposition to Disease , Genomics , Humans , Male , Polymorphism, Single Nucleotide , Prostatic Neoplasms/geneticsABSTRACT
Although prostate cancer is the leading cause of cancer mortality for African men, the vast majority of known disease associations have been detected in European study cohorts. Furthermore, most genome-wide association studies have used genotyping arrays that are hindered by SNP ascertainment bias. To overcome these disparities in genomic medicine, the Men of African Descent and Carcinoma of the Prostate (MADCaP) Network has developed a genotyping array that is optimized for African populations. The MADCaP Array contains more than 1.5 million markers and an imputation backbone that successfully tags over 94% of common genetic variants in African populations. This array also has a high density of markers in genomic regions associated with cancer susceptibility, including 8q24. We assessed the effectiveness of the MADCaP Array by genotyping 399 prostate cancer cases and 403 controls from seven urban study sites in sub-Saharan Africa. Samples from Ghana and Nigeria clustered together, whereas samples from Senegal and South Africa yielded distinct ancestry clusters. Using the MADCaP array, we identified cancer-associated loci that have large allele frequency differences across African populations. Polygenic risk scores for prostate cancer were higher in Nigeria than in Senegal. In summary, individual and population-level differences in prostate cancer risk were revealed using a novel genotyping array. SIGNIFICANCE: This study presents an Africa-specific genotyping array, which enables investigators to identify novel disease associations and to fine-map genetic loci that are associated with prostate and other cancers.
Subject(s)
Black People/genetics , Genetic Predisposition to Disease , Neoplasms/epidemiology , Neoplasms/genetics , Polymorphism, Single Nucleotide , Prostatic Neoplasms/epidemiology , Prostatic Neoplasms/genetics , Case-Control Studies , Cohort Studies , Genetic Loci , Genetics, Population , Genome-Wide Association Study , Humans , Male , Neoplasms/classification , Prostatic Neoplasms/classification , Risk Factors , South Africa/epidemiologyABSTRACT
To avoid an antibiotic resistance crisis, we need to develop antibiotics at a pace that matches the rate of evolution of resistance. However, the complex functions performed by antibiotics-combining, e.g., penetration of membranes, counteraction of resistance mechanisms, and interaction with molecular targets-have proven hard to achieve with current methods for drug development, including target-based screening and rational design. Here, we argue that we can meet the evolution of resistance in the clinic with evolution of antibiotics in the laboratory. On the basis of the results of experimental evolution studies of microbes in general and antibiotic production in Actinobacteria in particular, we propose methodology for evolving antibiotics to circumvent mechanisms of resistance. This exploits the ability of evolution to find solutions to complex problems without a need for design. We review evolutionary theory critical to this approach and argue that it is feasible and has important advantages over current methods for antibiotic discovery.
Subject(s)
Anti-Bacterial Agents/pharmacology , Bacteria/drug effects , Biological Products , Drug Development/methods , Drug Resistance, Microbial , Actinobacteria/chemistry , Actinobacteria/physiology , Evolution, MolecularABSTRACT
BACKGROUND: Accurate assessment of health disparities requires unbiased knowledge of genetic risks in different populations. Unfortunately, most genome-wide association studies use genotyping arrays and European samples. Here, we integrate whole genome sequence data from global populations, results from thousands of genome-wide association studies (GWAS), and extensive computer simulations to identify how genetic disease risks can be misestimated. RESULTS: In contrast to null expectations, we find that risk allele frequencies at known disease loci are significantly different for African populations compared to other continents. Strikingly, ancestral risk alleles are found at 9.51% higher frequency in Africa, and derived risk alleles are found at 5.40% lower frequency in Africa. By simulating GWAS with different study populations, we find that non-African cohorts yield disease associations that have biased allele frequencies and that African cohorts yield disease associations that are relatively free of bias. We also find empirical evidence that genotyping arrays and SNP ascertainment bias contribute to continental differences in risk allele frequencies. Because of these causes, polygenic risk scores can be grossly misestimated for individuals of African descent. Importantly, continental differences in risk allele frequencies are only moderately reduced if GWAS use whole genome sequences and hundreds of thousands of cases and controls. Finally, comparisons between uncorrected and corrected genetic risk scores reveal the benefits of considering whether risk alleles are ancestral or derived. CONCLUSIONS: Our results imply that caution must be taken when extrapolating GWAS results from one population to predict disease risks in another population.
Subject(s)
Genetic Predisposition to Disease , Genetics, Population , Genome, Human , Genome-Wide Association Study , Polymorphism, Single Nucleotide , Asian People , Black People , Gene Frequency , Global Health , Humans , Risk Factors , White PeopleABSTRACT
PURPOSE: Cancer of the prostate (CaP) is the leading cancer among men in sub-Saharan Africa (SSA). A substantial proportion of these men with CaP are diagnosed at late (usually incurable) stages, yet little is known about the etiology of CaP in SSA. METHODS: We established the Men of African Descent and Carcinoma of the Prostate Network, which includes seven SSA centers partnering with five US centers to study the genetics and epidemiology of CaP in SSA. We developed common data elements and instruments, regulatory infrastructure, and biosample collection, processing, and shipping protocols. We tested this infrastructure by collecting epidemiologic, medical record, and genomic data from a total of 311 patients with CaP and 218 matched controls recruited at the seven SSA centers. We extracted genomic DNA from whole blood, buffy coat, or buccal swabs from 265 participants and shipped it to the Center for Inherited Disease Research (Baltimore, MD) and the Centre for Proteomics and Genomics Research (Cape Town, South Africa), where genotypes were generated using the UK Biobank Axiom Array. RESULTS: We used common instruments for data collection and entered data into the shared database. Double-entered data from pilot participants showed a 95% to 98% concordance rate, suggesting that data can be collected, entered, and stored with a high degree of accuracy. Genotypes were obtained from 95% of tested DNA samples (100% from blood-derived DNA samples) with high concordance across laboratories. CONCLUSION: We provide approaches that can produce high-quality epidemiologic and genomic data in multicenter studies of cancer in SSA.
Subject(s)
Carcinoma/epidemiology , Carcinoma/genetics , Prostatic Neoplasms/epidemiology , Prostatic Neoplasms/genetics , Baltimore , Black People , Carcinoma/pathology , Genomics , Genotype , Humans , Male , Prostate/pathology , Prostatic Neoplasms/pathology , South Africa/epidemiologyABSTRACT
Phylogeny estimation is difficult for closely related populations and species, especially if they have been exchanging genes. We present a hierarchical Bayesian, Markov-chain Monte Carlo method with a state space that includes all possible phylogenies in a full Isolation-with-Migration model framework. The method is based on a new type of genealogy augmentation called a "hidden genealogy" that enables efficient updating of the phylogeny. This is the first likelihood-based method to fully incorporate directional gene flow and genetic drift for estimation of a species or population phylogeny. Application to human hunter-gatherer populations from Africa revealed a clear phylogenetic history, with strong support for gene exchange with an unsampled ghost population, and relatively ancient divergence between a ghost population and modern human populations, consistent with human/archaic divergence. In contrast, a study of five chimpanzee populations reveals a clear phylogeny with several pairs of populations having exchanged DNA, but does not support a history with an unsampled ghost population.
Subject(s)
Gene Flow , Genetic Techniques , Phylogeny , Animals , Bayes Theorem , Genetic Drift , Human Migration , Humans , Monte Carlo Method , Pan troglodytes/geneticsABSTRACT
Insertional activity of transposable elements (TEs) has had a major impact on the human genome; approximately one-half to two-thirds of the genome sequence is likely to be derived from TE insertions. Several families of human TEs - primarily Alu, L1 and SVA - continue to actively transpose, thereby generating insertional polymorphisms among individual genomes. The impact that TE insertions have on their human hosts' fitness, and accordingly the role that natural selection plays in shaping patterns of TE polymorphisms among populations, have yet to be systematically evaluated using whole genome sequence data. We present here a population genomic study of the effects of natural selection on human genetic variation that results from the recent activity of TEs. We developed a genome-wide scan for selection on human TE polymorphisms and applied it to a dataset of 14,384 locus-specific TE insertions characterized for 1511 individuals from 15 populations. Our TE selection scan looks for anomalously high population-specific TE insertion allele frequencies that are consistent with the action of positive (adaptive) selection. To control for the effects of demographic history, we compared the observed patterns of population-specific TE insertion allele frequencies to a neutral evolutionary model generated using time forward simulation of TE insertion allele frequencies among human population groups. This approach uncovered seven cases of polymorphic TE insertions that appear to have increased in frequency within specific human populations owing to the effects of positive selection. Five of the seven putatively selected TE insertions map to tissue-specific enhancers, and two cases correspond to expression quantitative trait loci that are associated with inter-individual gene regulatory differences. This study represents the first report of recent, local adaptation acting on polymorphic human TEs.