The need to diversify genomic studies: Insights from Andean highlanders and Amazonians.

More globally diverse perspectives are needed in genomic studies and precision medicine practices on non-Europeans. Here, we illustrate this by discussing the distribution of clinically actionable genetic variants involved in drug response in Andean highlanders and Amazonians, considering their environment, history, genetic structure, and historical biases in the perception of biological diversity of Native Americans.

Native American Genomic Diversity through Ancient DNA.

Ancient DNA is a powerful tool to understand the evolutionary dynamics of both current and ancestral populations. Posth et al. use ancient DNA to elucidate important questions surrounding the peopling of Central and South America, giving us greater insights into the ancestry of genetically understudied populations.

A positively selected FBN1 missense variant reduces height in Peruvian individuals.

On average, Peruvian individuals are among the shortest in the world1. Here we show that Native American ancestry is associated with reduced height in an ethnically diverse group of Peruvian individuals, and identify a population-specific, missense variant in the FBN1 gene (E1297G) that is significantly associated with lower height. Each copy of the minor allele (frequency of 4.7%) reduces height by 2.2 cm (4.4 cm in homozygous individuals). To our knowledge, this is the largest effect size known for a common height-associated variant. FBN1 encodes the extracellular matrix protein fibrillin 1, which is a major structural component of microfibrils. We observed less densely packed fibrillin-1-rich microfibrils with irregular edges in the skin of individuals who were homozygous for G1297 compared with individuals who were homozygous for E1297. Moreover, we show that the E1297G locus is under positive selection in non-African populations, and that the E1297 variant shows subtle evidence of positive selection specifically within the Peruvian population. This variant is also significantly more frequent in coastal Peruvian populations than in populations from the Andes or the Amazon, which suggests that short stature might be the result of adaptation to factors that are associated with the coastal environment in Peru.

Racial, ethnic and country of origin disparities in aggressive endometrial cancer histologic subtypes.

OBJECTIVE: This study investigated the risk of an aggressive endometrial cancer (EC) diagnosis by race, ethnicity, and country of origin to further elucidate histologic disparities in non-Hispanic Black (NHB), Hispanic, Asian/Pacific Islander (API), American Indian/Alaskan Native (AIAN) vs. non-Hispanic White (NHW) patients, particularly in Hispanic or API subgroups. METHODS: Patient diagnosed between 2004 and 2020 with low grade (LG)-endometrioid endometrial cancer (ECC) or an aggressive EC including grade 3 EEC, serous carcinoma, clear cell carcinoma, mixed epithelial carcinoma, or carcinosarcoma in the National Cancer Database were studied. The odds ratio (OR) and 95% confidence interval (CI) for diagnosis of an aggressive EC histology was estimated using logistic modeling. RESULTS: There were 343,868 NHW, 48,897 NHB, 30,013 Hispanic, 15,015 API and 1646 AIAN patients. The OR (95% CI) for an aggressive EC diagnosis was 3.07 (3.01-3.13) for NHB, 1.08 (1.06-1.11) for Hispanic, 1.17 (1.13-1.21) for API and 1.07 (0.96-1.19) for AIAN, relative to NHW patients. Subset analyses by country of origin illustrated the diversity in the OR for an aggressive EC diagnosis among Hispanic (1.18 for Mexican to 1.87 for Dominican), Asian (1.14 Asian Indian-Pakistani to 1.48 Korean) and Pacific Islander (1.00 for Hawaiian to 1.33 for Samoan) descendants. Hispanic, API and AIAN patients were diagnosed 5-years younger that NHW patients, and the risk for an aggressive EC histology were all significantly higher than NHW patients after correcting for age. Insurance status was another independent risk factor for aggressive histology. CONCLUSIONS: Risk of an aggressive EC diagnosis varied by race, ethnicity, and country of origin. NHB patients had the highest risk, followed by Dominican, South/Central American, Cuban, Korean, Thai, Vietnamese, and Filipino descendants.

Survival disparities in non-Hispanic Black and White cervical cancer patients vary by histology and are largely explained by modifiable factors.

PURPOSE: We investigated racial disparities in survival by histology in cervical cancer and examined the factors contributing to these disparities. METHODS: Non-Hispanic Black and non-Hispanic White (hereafter known as Black and White) patients with stage I-IV cervical carcinoma diagnosed between 2004 and 2017 in the National Cancer Database were studied. Survival differences were compared using Cox modeling to estimate hazard ratio (HR) or adjusted HR (AHR) and 95% confidence interval (CI). The contribution of demographic, socioeconomic and clinical factors to the Black vs White differences in survival was estimated after applying propensity score weighting in patients with squamous cell carcinoma (SCC) or adenocarcinoma (AC). RESULTS: This study included 10,111 Black and 43,252 White patients with cervical cancer. Black patients had worse survival than White cervical cancer patients (HR = 1.40, 95% CI = 1.35-1.45). Survival disparities between Black and White patients varied significantly by histology (HR = 1.20, 95% CI = 1.15-1.24 for SCC; HR = 2.32, 95% CI = 2.12-2.54 for AC, interaction p < 0.0001). After balancing the selected demographic, socioeconomic and clinical factors, survival in Black vs. White patients was no longer different in those with SCC (AHR = 1.01, 95% CI 0.97-1.06) or AC (AHR = 1.09, 95% CI = 0.96-1.24). In SCC, the largest contributors to survival disparities were neighborhood income and insurance. In AC, age was the most significant contributor followed by neighborhood income, insurance, and stage. Diagnosis of AC (but not SCC) at ≥65 years old was more common in Black vs. White patients (26% vs. 13%, respectively). CONCLUSIONS: Histology matters in survival disparities and diagnosis at ≥65 years old between Black and White cervical cancer patients. These disparities were largely explained by modifiable factors.

X-Chromosome Association Study in Latin American Cohorts Identifies New Loci in Parkinson's Disease.

BACKGROUND: Sex differences in Parkinson's disease (PD) risk are well-known. However, the role of sex chromosomes in the development and progression of PD is still unclear. OBJECTIVE: The objective of this study was to perform the first X-chromosome-wide association study for PD risk in a Latin American cohort. METHODS: We used data from three admixed cohorts: (1) Latin American Research consortium on the Genetics of Parkinson's Disease (n = 1504) as discover cohort, and (2) Latino cohort from International Parkinson Disease Genomics Consortium (n = 155) and (3) Bambui Aging cohort (n = 1442) as replication cohorts. We also developed an X-chromosome framework specifically designed for admixed populations. RESULTS: We identified eight linkage disequilibrium regions associated with PD. We replicated one of these regions (top variant rs525496; discovery odds ratio [95% confidence interval]: 0.60 [0.478-0.77], P = 3.13 × 10-5 replication odds ratio: 0.60 [0.37-0.98], P = 0.04). rs5525496 is associated with multiple expression quantitative trait loci in brain and non-brain tissues, including RAB9B, H2BFM, TSMB15B, and GLRA4, but colocalization analysis suggests that rs5525496 may not mediate risk by expression of these genes. We also replicated a previous X-chromosome-wide association study finding (rs28602900), showing that this variant is associated with PD in non-European populations. CONCLUSIONS: Our results reinforce the importance of including X-chromosome and diverse populations in genetic studies. © 2023 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Evolutionary history of modern Samoans.

Archaeological studies estimate the initial settlement of Samoa at 2,750 to 2,880 y ago and identify only limited settlement and human modification to the landscape until about 1,000 to 1,500 y ago. At this point, a complex history of migration is thought to have begun with the arrival of people sharing ancestry with Near Oceanic groups (i.e., Austronesian-speaking and Papuan-speaking groups), and was then followed by the arrival of non-Oceanic groups during European colonialism. However, the specifics of this peopling are not entirely clear from the archaeological and anthropological records, and is therefore a focus of continued debate. To shed additional light on the Samoan population history that this peopling reflects, we employ a population genetic approach to analyze 1,197 Samoan high-coverage whole genomes. We identify population splits between the major Samoan islands and detect asymmetrical gene flow to the capital city. We also find an extreme bottleneck until about 1,000 y ago, which is followed by distinct expansions across the islands and subsequent bottlenecks consistent with European colonization. These results provide for an increased understanding of Samoan population history and the dynamics that inform it, and also demonstrate how rapid demographic processes can shape modern genomes.

The genetic structure and adaptation of Andean highlanders and Amazonians are influenced by the interplay between geography and culture.

Western South America was one of the worldwide cradles of civilization. The well-known Inca Empire was the tip of the iceberg of an evolutionary process that started 11,000 to 14,000 years ago. Genetic data from 18 Peruvian populations reveal the following: 1) The between-population homogenization of the central southern Andes and its differentiation with respect to Amazonian populations of similar latitudes do not extend northward. Instead, longitudinal gene flow between the northern coast of Peru, Andes, and Amazonia accompanied cultural and socioeconomic interactions revealed by archeology. This pattern recapitulates the environmental and cultural differentiation between the fertile north, where altitudes are lower, and the arid south, where the Andes are higher, acting as a genetic barrier between the sharply different environments of the Andes and Amazonia. 2) The genetic homogenization between the populations of the arid Andes is not only due to migrations during the Inca Empire or the subsequent colonial period. It started at least during the earlier expansion of the Wari Empire (600 to 1,000 years before present). 3) This demographic history allowed for cases of positive natural selection in the high and arid Andes vs. the low Amazon tropical forest: in the Andes, a putative enhancer in HAND2-AS1 (heart and neural crest derivatives expressed 2 antisense RNA1, a noncoding gene related to cardiovascular function) and rs269868-C/Ser1067 in DUOX2 (dual oxidase 2, related to thyroid function and innate immunity) genes and, in the Amazon, the gene encoding for the CD45 protein, essential for antigen recognition by T and B lymphocytes in viral-host interaction.

De novo mutations across 1,465 diverse genomes reveal mutational insights and reductions in the Amish founder population.

De novo mutations (DNMs), or mutations that appear in an individual despite not being seen in their parents, are an important source of genetic variation whose impact is relevant to studies of human evolution, genetics, and disease. Utilizing high-coverage whole-genome sequencing data as part of the Trans-Omics for Precision Medicine (TOPMed) Program, we called 93,325 single-nucleotide DNMs across 1,465 trios from an array of diverse human populations, and used them to directly estimate and analyze DNM counts, rates, and spectra. We find a significant positive correlation between local recombination rate and local DNM rate, and that DNM rate explains a substantial portion (8.98 to 34.92%, depending on the model) of the genome-wide variation in population-level genetic variation from 41K unrelated TOPMed samples. Genome-wide heterozygosity does correlate with DNM rate, but only explains <1% of variation. While we are underpowered to see small differences, we do not find significant differences in DNM rate between individuals of European, African, and Latino ancestry, nor across ancestrally distinct segments within admixed individuals. However, we did find significantly fewer DNMs in Amish individuals, even when compared with other Europeans, and even after accounting for parental age and sequencing center. Specifically, we found significant reductions in the number of C→A and T→C mutations in the Amish, which seem to underpin their overall reduction in DNMs. Finally, we calculated near-zero estimates of narrow sense heritability (h2), which suggest that variation in DNM rate is significantly shaped by nonadditive genetic effects and the environment.

STRIDE: a command-line HMM-based identifier and sub-classifier of Plasmodium falciparum RIFIN and STEVOR variant surface antigen families.

BACKGROUND: RIFINs and STEVORs are variant surface antigens expressed by P. falciparum that play roles in severe malaria pathogenesis and immune evasion. These two highly diverse multigene families feature multiple paralogs, making their classification challenging using traditional bioinformatic methods. RESULTS: STRIDE (STevor and RIfin iDEntifier) is an HMM-based, command-line program that automates the identification and classification of RIFIN and STEVOR protein sequences in the malaria parasite Plasmodium falciparum. STRIDE is more sensitive in detecting RIFINs and STEVORs than available PFAM and TIGRFAM tools and reports RIFIN subtypes and the number of sequences with a FHEYDER amino acid motif, which has been associated with severe malaria pathogenesis. CONCLUSIONS: STRIDE will be beneficial to malaria research groups analyzing genome sequences and transcripts of clinical field isolates, providing insight into parasite biology and virulence.

Characterizing the Genetic Architecture of Parkinson's Disease in Latinos.

OBJECTIVE: This work was undertaken in order to identify Parkinson's disease (PD) risk variants in a Latino cohort, to describe the overlap in the genetic architecture of PD in Latinos compared to European-ancestry subjects, and to increase the diversity in PD genome-wide association (GWAS) data. METHODS: We genotyped and imputed 1,497 PD cases and controls recruited from nine clinical sites across South America. We performed a GWAS using logistic mixed models; variants with a p-value <1 × 10-5 were tested in a replication cohort of 1,234 self-reported Latino PD cases and 439,522 Latino controls from 23andMe, Inc. We also performed an admixture mapping analysis where local ancestry blocks were tested for association with PD status. RESULTS: One locus, SNCA, achieved genome-wide significance (p-value <5 × 10-8 ); rs356182 achieved genome-wide significance in both the discovery and the replication cohorts (discovery, G allele: 1.58 OR, 95% CI 1.35-1.86, p-value 2.48 × 10-8 ; 23andMe, G allele: 1.26 OR, 95% CI 1.16-1.37, p-value 4.55 × 10-8 ). In our admixture mapping analysis, a locus on chromosome 14, containing the gene STXBP6, achieved significance in a joint test of ancestries and in the Native American single-ancestry test (p-value <5 × 10-5 ). A second locus on chromosome 6, containing the gene RPS6KA2, achieved significance in the African single-ancestry test (p-value <5 × 10-5 ). INTERPRETATION: This study demonstrated the importance of the SNCA locus for the etiology of PD in Latinos. By leveraging the demographic history of our cohort via admixture mapping, we identified two potential PD risk loci that merit further study. ANN NEUROL 2021;90:353-365.

Admixture/fine-mapping in Brazilians reveals a West African associated potential regulatory variant (rs114066381) with a strong female-specific effect on body mass and fat mass indexes.

BACKGROUND/OBJECTIVES: Admixed populations are a resource to study the global genetic architecture of complex phenotypes, which is critical, considering that non-European populations are severely underrepresented in genomic studies. Here, we study the genetic architecture of BMI in children, young adults, and elderly individuals from the admixed population of Brazil. SUBJECTS/METHODS: Leveraging admixture in Brazilians, whose chromosomes are mosaics of fragments of Native American, European, and African origins, we used genome-wide data to perform admixture mapping/fine-mapping of body mass index (BMI) in three Brazilian population-based cohorts from Northeast (Salvador), Southeast (Bambuí), and South (Pelotas). RESULTS: We found significant associations with African-associated alleles in children from Salvador (PALD1 and ZMIZ1 genes), and in young adults from Pelotas (NOD2 and MTUS2 genes). More importantly, in Pelotas, rs114066381, mapped in a potential regulatory region, is significantly associated only in females (p = 2.76e-06). This variant is rare in Europeans but with frequencies of ~3% in West Africa and has a strong female-specific effect (95% CI: 2.32-5.65 kg/m2 per each A allele). We confirmed this sex-specific association and replicated its strong effect for an adjusted fat mass index in the same Pelotas cohort, and for BMI in another Brazilian cohort from São Paulo (Southeast Brazil). A meta-analysis confirmed the significant association. Remarkably, we observed that while the frequency of rs114066381-A allele ranges from 0.8 to 2.1% in the studied populations, it attains ~9% among women with morbid obesity from Pelotas, São Paulo, and Bambuí. The effect size of rs114066381 is at least five times higher than the FTO SNPs rs9939609 and rs1558902, already emblematic for their high effects. CONCLUSIONS: We identified six candidate SNPs associated with BMI. rs114066381 stands out for its high effect that was replicated and its high frequency in women with morbid obesity. We demonstrate how admixed populations are a source of new relevant phenotype-associated genetic variants.

Genome-Wide Analysis of Copy Number Variation in Latin American Parkinson's Disease Patients.

BACKGROUND: Parkinson's disease is the second most common neurodegenerative disorder and affects people from all ethnic backgrounds, yet little is known about the genetics of Parkinson's disease in non-European populations. In addition, the overall identification of copy number variants at a genome-wide level has been understudied in Parkinson's patients. The objective of this study was to understand the genome-wide burden of copy number variants in Latinos and its association with Parkinson's disease. METHODS: We used genome-wide genotyping data from 747 Parkinson's disease patients and 632 controls from the Latin American Research Consortium on the Genetics of Parkinson's disease. RESULTS: Genome-wide copy number burden analysis showed that patients were significantly enriched for copy number variants overlapping known Parkinson's disease genes compared with controls (odds ratio, 3.97; 95%CI, 1.69-10.5; P = 0.018). PRKN showed the strongest copy number burden, with 20 copy number variant carriers. These patients presented an earlier age of disease onset compared with patients with other copy number variants (median age at onset, 31 vs 57 years, respectively; P = 7.46 × 10-7 ). CONCLUSIONS: We found that although overall genome-wide copy number variant burden was not significantly different, Parkinson's disease patients were significantly enriched with copy number variants affecting known Parkinson's disease genes. We also identified that of 250 patients with early-onset disease, 5.6% carried a copy number variant on PRKN in our cohort. Our study is the first to analyze genome-wide copy number variant association in Latino Parkinson's disease patients and provides insights about this complex disease in this understudied population. © 2020 International Parkinson and Movement Disorder Society.

Inferring the "Primordial Immune Complex": Origins of MHC Class I and Antigen Receptors Revealed by Comparative Genomics.

Comparative analyses suggest that the MHC was derived from a prevertebrate "primordial immune complex" (PIC). PIC duplicated twice in the well-studied two rounds of genome-wide duplications (2R) early in vertebrate evolution, generating four MHC paralogous regions (predominantly on human chromosomes [chr] 1, 6, 9, 19). Examining chiefly the amphibian Xenopus laevis, but also other vertebrates, we identified their MHC paralogues and mapped MHC class I, AgR, and "framework" genes. Most class I genes mapped to MHC paralogues, but a cluster of Xenopus MHC class Ib genes (xnc), which previously was mapped outside of the MHC paralogues, was surrounded by genes syntenic to mammalian CD1 genes, a region previously proposed as an MHC paralogue on human chr 1. Thus, this gene block is instead the result of a translocation that we call the translocated part of the MHC paralogous region (MHCtrans) Analyses of Xenopus class I genes, as well as MHCtrans, suggest that class I arose at 1R on the chr 6/19 ancestor. Of great interest are nonrearranging AgR-like genes mapping to three MHC paralogues; thus, PIC clearly contained several AgR precursor loci, predating MHC class I/II. However, all rearranging AgR genes were found on paralogues derived from the chr 19 precursor, suggesting that invasion of a variable (V) exon by the RAG transposon occurred after 2R. We propose models for the evolutionary history of MHC/TCR/Ig and speculate on the dichotomy between the jawless (lamprey and hagfish) and jawed vertebrate adaptive immune systems, as we found genes related to variable lymphocyte receptors also map to MHC paralogues.

Evolutionary genomic dynamics of Peruvians before, during, and after the Inca Empire.

Native Americans from the Amazon, Andes, and coastal geographic regions of South America have a rich cultural heritage but are genetically understudied, therefore leading to gaps in our knowledge of their genomic architecture and demographic history. In this study, we sequence 150 genomes to high coverage combined with an additional 130 genotype array samples from Native American and mestizo populations in Peru. The majority of our samples possess greater than 90% Native American ancestry, which makes this the most extensive Native American sequencing project to date. Demographic modeling reveals that the peopling of Peru began ∼12,000 y ago, consistent with the hypothesis of the rapid peopling of the Americas and Peruvian archeological data. We find that the Native American populations possess distinct ancestral divisions, whereas the mestizo groups were admixtures of multiple Native American communities that occurred before and during the Inca Empire and Spanish rule. In addition, the mestizo communities also show Spanish introgression largely following Peruvian Independence, nearly 300 y after Spain conquered Peru. Further, we estimate migration events between Peruvian populations from all three geographic regions with the majority of between-region migration moving from the high Andes to the low-altitude Amazon and coast. As such, we present a detailed model of the evolutionary dynamics which impacted the genomes of modern-day Peruvians and a Native American ancestry dataset that will serve as a beneficial resource to addressing the underrepresentation of Native American ancestry in sequencing studies.

Human-SARS-CoV-2 interactome and human genetic diversity: TMPRSS2-rs2070788, associated with severe influenza, and its population genetics caveats in Native Americans.

For human/SARS-CoV-2 interactome genes ACE2, TMPRSS2 and BSG, there is a convincing evidence of association in Asians with influenza-induced SARS for TMPRSS2-rs2070788, tag-SNP of the eQTL rs383510. This case illustrates the importance of population genetics and of sequencing data in the design of genetic association studies in different human populations: the high linkage disequilibrium (LD) between rs2070788 and rs383510 is Asian-specific. Leveraging on a combination of genotyping and sequencing data for Native Americans (neglected in genetic studies), we show that while their frequencies of the Asian tag-SNP rs2070788 is, surprisingly, the highest worldwide, it is not in LD with the eQTL rs383510, that therefore, should be directly genotyped in genetic association studies of SARS in populations with Native American ancestry.

Detecting geospatial patterns of Plasmodium falciparum parasite migration in Cambodia using optimized estimated effective migration surfaces.

BACKGROUND: Understanding the genetic structure of natural populations provides insight into the demographic and adaptive processes that have affected those populations. Such information, particularly when integrated with geospatial data, can have translational applications for a variety of fields, including public health. Estimated effective migration surfaces (EEMS) is an approach that allows visualization of the spatial patterns in genomic data to understand population structure and migration. In this study, we developed a workflow to optimize the resolution of spatial grids used to generate EEMS migration maps and applied this optimized workflow to estimate migration of Plasmodium falciparum in Cambodia and bordering regions of Thailand and Vietnam. METHODS: The optimal density of EEMS grids was determined based on a new workflow created using density clustering to define genomic clusters and the spatial distance between genomic clusters. Topological skeletons were used to capture the spatial distribution for each genomic cluster and to determine the EEMS grid density; i.e., both genomic and spatial clustering were used to guide the optimization of EEMS grids. Model accuracy for migration estimates using the optimized workflow was tested and compared to grid resolutions selected without the optimized workflow. As a test case, the optimized workflow was applied to genomic data generated from P. falciparum sampled in Cambodia and bordering regions, and migration maps were compared to estimates of malaria endemicity, as well as geographic properties of the study area, as a means of validating observed migration patterns. RESULTS: Optimized grids displayed both high model accuracy and reduced computing time compared to grid densities selected in an unguided manner. In addition, EEMS migration maps generated for P. falciparum using the optimized grid corresponded to estimates of malaria endemicity and geographic properties of the study region that might be expected to impact malaria parasite migration, supporting the validity of the observed migration patterns. CONCLUSIONS: Optimized grids reduce spatial uncertainty in the EEMS contours that can result from user-defined parameters, such as the resolution of the spatial grid used in the model. This workflow will be useful to a broad range of EEMS users as it can be applied to analyses involving other organisms of interest and geographic areas.

Ancestral characterization of 1018 cancer cell lines highlights disparities and reveals gene expression and mutational differences.

BACKGROUND: Although cell lines are an essential resource for studying cancer biology, many are of unknown ancestral origin, and their use may not be optimal for evaluating the biology of all patient populations. METHODS: An admixture analysis was performed using genome-wide chip data from the Catalogue of Somatic Mutations in Cancer (COSMIC) Cell Lines Project to calculate genetic ancestry estimates for 1018 cancer cell lines. After stratifying the analyses by tissue and histology types, linear models were used to evaluate the influence of ancestry on gene expression and somatic mutation frequency. RESULTS: For the 701 cell lines with unreported ancestry, 215 were of East Asian origin, 30 were of African or African American origin, and 453 were of European origin. Notable imbalances were observed in ancestral representation across tissue type, with the majority of analyzed tissue types having few cell lines of African American ancestral origin, and with Hispanic and South Asian ancestry being almost entirely absent across all cell lines. In evaluating gene expression across these cell lines, expression levels of the genes neurobeachin line 1 (NBEAL1), solute carrier family 6 member 19 (SLC6A19), HEAT repeat containing 6 (HEATR6), and epithelial cell transforming 2 like (ECT2L) were associated with ancestry. Significant differences were also observed in the proportions of somatic mutation types across cell lines with varying ancestral proportions. CONCLUSIONS: By estimating genetic ancestry for 1018 cancer cell lines, the authors have produced a resource that cancer researchers can use to ensure that their cell lines are ancestrally representative of the populations they intend to affect. Furthermore, the novel ancestry-specific signal identified underscores the importance of ancestral awareness when studying cancer.

Analysis of 6,515 exomes reveals the recent origin of most human protein-coding variants.

Establishing the age of each mutation segregating in contemporary human populations is important to fully understand our evolutionary history and will help to facilitate the development of new approaches for disease-gene discovery. Large-scale surveys of human genetic variation have reported signatures of recent explosive population growth, notable for an excess of rare genetic variants, suggesting that many mutations arose recently. To more quantitatively assess the distribution of mutation ages, we resequenced 15,336 genes in 6,515 individuals of European American and African American ancestry and inferred the age of 1,146,401 autosomal single nucleotide variants (SNVs). We estimate that approximately 73% of all protein-coding SNVs and approximately 86% of SNVs predicted to be deleterious arose in the past 5,000-10,000 years. The average age of deleterious SNVs varied significantly across molecular pathways, and disease genes contained a significantly higher proportion of recently arisen deleterious SNVs than other genes. Furthermore, European Americans had an excess of deleterious variants in essential and Mendelian disease genes compared to African Americans, consistent with weaker purifying selection due to the Out-of-Africa dispersal. Our results better delimit the historical details of human protein-coding variation, show the profound effect of recent human history on the burden of deleterious SNVs segregating in contemporary populations, and provide important practical information that can be used to prioritize variants in disease-gene discovery.