Genomic evidence for shared common ancestry of East African hunting-gathering populations and insights into local adaptation.

Anatomically modern humans arose in Africa ∼300,000 years ago, but the demographic and adaptive histories of African populations are not well-characterized. Here, we have generated a genome-wide dataset from 840 Africans, residing in western, eastern, southern, and northern Africa, belonging to 50 ethnicities, and speaking languages belonging to four language families. In addition to agriculturalists and pastoralists, our study includes 16 populations that practice, or until recently have practiced, a hunting-gathering (HG) lifestyle. We observe that genetic structure in Africa is broadly correlated not only with geography, but to a lesser extent, with linguistic affiliation and subsistence strategy. Four East African HG (EHG) populations that are geographically distant from each other show evidence of common ancestry: the Hadza and Sandawe in Tanzania, who speak languages with clicks classified as Khoisan; the Dahalo in Kenya, whose language has remnant clicks; and the Sabue in Ethiopia, who speak an unclassified language. Additionally, we observed common ancestry between central African rainforest HGs and southern African San, the latter of whom speak languages with clicks classified as Khoisan. With the exception of the EHG, central African rainforest HGs, and San, other HG groups in Africa appear genetically similar to neighboring agriculturalist or pastoralist populations. We additionally demonstrate that infectious disease, immune response, and diet have played important roles in the adaptive landscape of African history. However, while the broad biological processes involved in recent human adaptation in Africa are often consistent across populations, the specific loci affected by selective pressures more often vary across populations.

Adaptive Landscape of Protein Variation in Human Exomes.

The human genome contains hundreds of thousands of missense mutations. However, only a handful of these variants are known to be adaptive, which implies that adaptation through protein sequence change is an extremely rare phenomenon in human evolution. Alternatively, existing methods may lack the power to pinpoint adaptive variation. We have developed and applied an Evolutionary Probability Approach (EPA) to discover candidate adaptive polymorphisms (CAPs) through the discordance between allelic evolutionary probabilities and their observed frequencies in human populations. EPA reveals thousands of missense CAPs, which suggest that a large number of previously optimal alleles experienced a reversal of fortune in the human lineage. We explored nonadaptive mechanisms to explain CAPs, including the effects of demography, mutation rate variability, and negative and positive selective pressures in modern humans. Many nonadaptive hypotheses were tested, but failed to explain the data, which suggests that a large proportion of CAP alleles have increased in frequency due to beneficial selection. This suggestion is supported by the fact that a vast majority of adaptive missense variants discovered previously in humans are CAPs, and hundreds of CAP alleles are protective in genotype-phenotype association data. Our integrated phylogenomic and population genetic EPA approach predicts the existence of thousands of nonneutral candidate variants in the human proteome. We expect this collection to be enriched in beneficial variation. The EPA approach can be applied to discover candidate adaptive variation in any protein, population, or species for which allele frequency data and reliable multispecies alignments are available.

A Locus at 5q33.3 Confers Resistance to Tuberculosis in Highly Susceptible Individuals.

Immunosuppression resulting from HIV infection increases the risk of progression to active tuberculosis (TB) both in individuals newly exposed to Mycobacterium tuberculosis (MTB) and in those with latent infections. We hypothesized that HIV-positive individuals who do not develop TB, despite living in areas where it is hyperendemic, provide a model of natural resistance. We performed a genome-wide association study of TB resistance by using 581 HIV-positive Ugandans and Tanzanians enrolled in prospective cohort studies of TB; 267 of these individuals developed active TB, and 314 did not. A common variant, rs4921437 at 5q33.3, was significantly associated with TB (odds ratio = 0.37, p = 2.11 × 10(-8)). This variant lies within a genomic region that includes IL12B and is embedded in an H3K27Ac histone mark. The locus also displays consistent patterns of linkage disequilibrium across African populations and has signals of strong selection in populations from equatorial Africa. Along with prior studies demonstrating that therapy with IL-12 (the cytokine encoded in part by IL12B, associated with longer survival following MTB infection in mice deficient in CD4 T cells), our results suggest that this pathway might be an excellent target for the development of new modalities for treating TB, especially for HIV-positive individuals. Our results also indicate that studying extreme disease resistance in the face of extensive exposure can increase the power to detect associations in complex infectious disease.

Recent human adaptation: genomic approaches, interpretation and insights.

The recent availability of genomic data has spurred many genome-wide studies of human adaptation in different populations worldwide. Such studies have provided insights into novel candidate genes and pathways that are putatively involved in adaptation to different environments, diets and disease prevalence. However, much work is needed to translate these results into candidate adaptive variants that are biologically interpretable. In this Review, we discuss methods that may help to identify true biological signals of selection and studies that incorporate complementary phenotypic and functional data. We conclude with recommendations for future studies that focus on opportunities to use integrative genomics methodologies in human adaptation studies.

e-GRASP: an integrated evolutionary and GRASP resource for exploring disease associations.

BACKGROUND: Genome-wide association studies (GWAS) have become a mainstay of biological research concerned with discovering genetic variation linked to phenotypic traits and diseases. Both discrete and continuous traits can be analyzed in GWAS to discover associations between single nucleotide polymorphisms (SNPs) and traits of interest. Associations are typically determined by estimating the significance of the statistical relationship between genetic loci and the given trait. However, the prioritization of bona fide, reproducible genetic associations from GWAS results remains a central challenge in identifying genomic loci underlying common complex diseases. Evolutionary-aware meta-analysis of the growing GWAS literature is one way to address this challenge and to advance from association to causation in the discovery of genotype-phenotype relationships. DESCRIPTION: We have created an evolutionary GWAS resource to enable in-depth query and exploration of published GWAS results. This resource uses the publically available GWAS results annotated in the GRASP2 database. The GRASP2 database includes results from 2082 studies, 177 broad phenotype categories, and ~8.87 million SNP-phenotype associations. For each SNP in e-GRASP, we present information from the GRASP2 database for convenience as well as evolutionary information (e.g., rate and timespan). Users can, therefore, identify not only SNPs with highly significant phenotype-association P-values, but also SNPs that are highly replicated and/or occur at evolutionarily conserved sites that are likely to be functionally important. Additionally, we provide an evolutionary-adjusted SNP association ranking (E-rank) that uses cross-species evolutionary conservation scores and population allele frequencies to transform P-values in an effort to enhance the discovery of SNPs with a greater probability of biologically meaningful disease associations. CONCLUSION: By adding an evolutionary dimension to the GWAS results available in the GRASP2 database, our e-GRASP resource will enable a more effective exploration of SNPs not only by the statistical significance of trait associations, but also by the number of studies in which associations have been replicated, and the evolutionary context of the associated mutations. Therefore, e-GRASP will be a valuable resource for aiding researchers in the identification of bona fide, reproducible genetic associations from GWAS results. This resource is freely available at http://www.mypeg.info/egrasp .

Challenges in Translating GWAS Results to Clinical Care.

Clinical genetic testing for Mendelian disorders is standard of care in many cases; however, it is less clear to what extent and in which situations clinical genetic testing may improve preventive efforts, diagnosis and/or prognosis of complex disease. One challenge is that much of the reported research relies on tag single nucleotide polymorphisms (SNPs) to act as proxies for assumed underlying functional variants that are not yet known. Here we use coronary artery disease and melanoma as case studies to evaluate how well reported genetic risk variants tag surrounding variants across population samples in the 1000 Genomes Project Phase 3 data. We performed a simulation study where we randomly assigned a "functional" variant and evaluated how often this simulated functional variant was correctly tagged in diverse population samples. Our results indicate a relatively large error rate when generalizing increased genetic risk of complex disease across diverse population samples, even when generalizing within geographic regions. Our results further highlight the importance of including diverse populations in genome-wide association studies. Future work focused on identifying functional variants will eliminate the need for tag SNPs; however, until functional variants are known, caution should be used in the interpretation of genetic risk for complex disease using tag SNPs.

Patterns of ancestry, signatures of natural selection, and genetic association with stature in Western African pygmies.

African Pygmy groups show a distinctive pattern of phenotypic variation, including short stature, which is thought to reflect past adaptation to a tropical environment. Here, we analyze Illumina 1M SNP array data in three Western Pygmy populations from Cameroon and three neighboring Bantu-speaking agricultural populations with whom they have admixed. We infer genome-wide ancestry, scan for signals of positive selection, and perform targeted genetic association with measured height variation. We identify multiple regions throughout the genome that may have played a role in adaptive evolution, many of which contain loci with roles in growth hormone, insulin, and insulin-like growth factor signaling pathways, as well as immunity and neuroendocrine signaling involved in reproduction and metabolism. The most striking results are found on chromosome 3, which harbors a cluster of selection and association signals between approximately 45 and 60 Mb. This region also includes the positional candidate genes DOCK3, which is known to be associated with height variation in Europeans, and CISH, a negative regulator of cytokine signaling known to inhibit growth hormone-stimulated STAT5 signaling. Finally, pathway analysis for genes near the strongest signals of association with height indicates enrichment for loci involved in insulin and insulin-like growth factor signaling.

Using the Coriell Personalized Medicine Collaborative Data to conduct a genome-wide association study of sleep duration.

Sleep is critical to health and functionality, and several studies have investigated the inherited component of insomnia and other sleep disorders using genome-wide association studies (GWAS). However, genome-wide studies focused on sleep duration are less common. Here, we used data from participants in the Coriell Personalized Medicine Collaborative (CPMC) (n = 4,401) to examine putative associations between self-reported sleep duration, demographic and lifestyle variables, and genome-wide single nucleotide polymorphism (SNP) data to better understand genetic contributions to variation in sleep duration. We employed stepwise ordered logistic regression to select our model and retained the following predictive variables: age, gender, weight, physical activity, physical activity at work, smoking status, alcohol consumption, ethnicity, and ancestry (as measured by principal components analysis) in our association testing. Several of our strongest candidate genes were previously identified in GWAS related to sleep duration (TSHZ2, ABCC9, FBXO15) and narcolepsy (NFATC2, SALL4). In addition, we have identified novel candidate genes for involvement in sleep duration including SORCS1 and ELOVL2. Our results demonstrate that the self-reported data collected through the CPMC are robust, and our genome-wide association analysis has identified novel candidate genes involved in sleep duration. More generally, this study contributes to a better understanding of the complexity of human sleep.

Colloquium paper: working toward a synthesis of archaeological, linguistic, and genetic data for inferring African population history.

Although Africa is the origin of modern humans, the pattern and distribution of genetic variation and correlations with cultural and linguistic diversity in Africa have been understudied. Recent advances in genomic technology, however, have led to genomewide studies of African samples. In this article, we discuss genetic variation in African populations contextualized with what is known about archaeological and linguistic variation. What emerges from this review is the importance of using independent lines of evidence in the interpretation of genetic and genomic data in the reconstruction of past population histories.

Postmortem toxicology findings from the Camden Opioid Research Initiative.

The United States continues to be impacted by decades of an opioid misuse epidemic, worsened by the COVID-19 pandemic and by the growing prevalence of highly potent synthetic opioids (HPSO) such as fentanyl. In instances of a toxicity event, first-response administration of reversal medications such as naloxone can be insufficient to fully counteract the effects of HPSO, particularly when there is co-occurring substance use. In an effort to characterize and study this multi-faceted problem, the Camden Opioid Research Initiative (CORI) has been formed. The CORI study has collected and analyzed post-mortem toxicology data from 42 cases of decedents who expired from opioid-related toxicity in the South New Jersey region to characterize substance use profiles. Co-occurring substance use, whether by intent or through possible contamination of the illicit opioid supply, is pervasive among deaths due to opioid toxicity, and evidence of medication-assisted treatment is scarce. Nearly all (98%) of the toxicology cases show the presence of the HPSO, fentanyl, and very few (7%) results detected evidence of medication-assisted treatment for opioid use disorder, such as buprenorphine or methadone, at the time of death. The opioid toxicity reversal drug, naloxone, was detected in 19% of cases, but 100% of cases expressed one or more stimulants, and sedatives including xylazine were detected in 48% of cases. These results showing complex substance use profiles indicate that efforts at mitigating the opioid misuse epidemic must address the complications presented by co-occurring stimulant and other substance use, and reduce barriers to and stigmas of seeking effective medication-assisted treatments.

Genome-wide insights into the patterns and determinants of fine-scale population structure in humans.

Studying genomic patterns of human population structure provides important insights into human evolutionary history and the relationship among populations, and it has significant practical implications for disease-gene mapping. Here we describe a principal component (PC)-based approach to studying intracontinental population structure in humans, identify the underlying markers mediating the observed patterns of fine-scale population structure, and infer the predominating evolutionary forces shaping local population structure. We applied this methodology to a data set of 650K SNPs genotyped in 944 unrelated individuals from 52 populations and demonstrate that, although typical PC analyses focus on the top axes of variation, substantial information about population structure is contained in lower-ranked PCs. We identified 18 significant PCs, some of which distinguish individual populations. In addition to visually representing sample clusters in PC biplots, we estimated the set of all SNPs significantly correlated with each of the most informative axes of variation. These polymorphisms, unlike ancestry-informative markers (AIMs), constitute a much larger set of loci that drive genomic signatures of population structure. The genome-wide distribution of these significantly correlated markers can largely be accounted for by the stochastic effects of genetic drift, although significant clustering does occur in genomic regions that have been previously implicated as targets of recent adaptive evolution.

Patient Perceptions and Potential Utility of Pharmacogenetic Testing in Chronic Pain Management and Opioid Use Disorder in the Camden Opioid Research Initiative.

Pharmacogenetics (PGx) has the potential to improve opioid medication management. Here, we present patient perception data, pharmacogenetic data and medication management trends in patients with chronic pain (arm 1) and opioid use disorder (arm 2) treated at Cooper University Health Care in Camden City, NJ. Our results demonstrate that the majority of patients in both arms of the study (55% and 65%, respectively) are open to pharmacogenetic testing, and most (66% and 69%, respectively) believe that genetic testing has the potential to improve their medical care. Our results further support the potential for CYP2D6 PGx testing to inform chronic pain medication management for poor metabolizers (PMs) and ultrarapid metabolizers (UMs). Future efforts to implement PGx testing in chronic pain management, however, must address patient concerns about genetic test result access and genetic discrimination.

Pharmacogenomics: From Basic Research to Clinical Implementation.

The established contribution of genetic variation to drug response has the potential to improve drug efficacy and reduce drug toxicity [...].

Common Treatment, Common Variant: Evolutionary Prediction of Functional Pharmacogenomic Variants.

Pharmacogenomics holds the promise of personalized drug efficacy optimization and drug toxicity minimization. Much of the research conducted to date, however, suffers from an ascertainment bias towards European participants. Here, we leverage publicly available, whole genome sequencing data collected from global populations, evolutionary characteristics, and annotated protein features to construct a new in silico machine learning pharmacogenetic identification method called XGB-PGX. When applied to pharmacogenetic data, XGB-PGX outperformed all existing prediction methods and identified over 2000 new pharmacogenetic variants. While there are modest pharmacogenetic allele frequency distribution differences across global population samples, the most striking distinction is between the relatively rare putatively neutral pharmacogene variants and the relatively common established and newly predicted functional pharamacogenetic variants. Our findings therefore support a focus on individual patient pharmacogenetic testing rather than on clinical presumptions about patient race, ethnicity, or ancestral geographic residence. We further encourage more attention be given to the impact of common variation on drug response and propose a new 'common treatment, common variant' perspective for pharmacogenetic prediction that is distinct from the types of variation that underlie complex and Mendelian disease. XGB-PGX has identified many new pharmacovariants that are present across all global communities; however, communities that have been underrepresented in genomic research are likely to benefit the most from XGB-PGX's in silico predictions.

Testing a best practices risk result format to communicate genetic risks.

OBJECTIVE: To investigate the effect of a genetic report format using risk communication "best-practices" on risk perceptions, in part to reduce risk overestimates. METHODS: Adults (N = 470) from the Coriell Personalized Medicine Collaborative (CPMC) were randomized to a 2 × 2 experimental design to receive a hypothetical "personalized" genetic risk result for leukemia (relative risk = 1.5 or 2.5) through either the standard CPMC report (N = 232) or an enriched report informed by best practices (N = 238). A one-time, online survey assessed numeracy and risk perceptions including "feelings of risk" and a numerical estimate. RESULTS: Regardless of numeracy, participants who received the enriched report had fewer overestimates of their lifetime risk estimate (LRE; odds ratio = 0.19, p < .001) and lower feelings of risk on two of three measures (p < .001). Participants with higher numeracy scores had fewer overestimates of LRE (OR = 0.66, p < .001) and lower feelings of risk on two out of three measures (p ≤ .01); the interaction between numeracy and report format was non-significant. CONCLUSION: The enriched report produced more accurate LRE and lower risk perceptions regardless of numeracy level, suggesting the enriched format was helpful to individuals irrespective of numeracy ability. PRACTICE IMPLICATIONS: Best practice elements in risk reports may help individuals form more accurate risk perceptions.

Population genomic analysis of ALMS1 in humans reveals a surprisingly complex evolutionary history.

Mutations in the human gene ALMS1 result in Alström Syndrome, which presents with early childhood obesity and insulin resistance leading to Type 2 diabetes. Previous genomewide scans for selection in the HapMap data based on linkage disequilibrium and population structure suggest that ALMS1 was subject to recent positive selection. Through a detailed population genomic analysis of existing genomewide data sets and new resequencing data obtained in geographically diverse populations, we find that the signature of selection at ALMS1 is considerably more complex than what would be expected for an idealized model of a selective sweep acting on a newly arisen advantageous mutation. Specifically, we observed three highly divergent and globally dispersed haplogroups, two of which carry a set of seven derived nonsynonymous single nucleotide polymorphisms that are nearly fixed in Asian populations. Our data suggest that the interaction of human demographic history and positive selection on standing variation in Eurasian populations approximately 15 thousand years ago parsimoniously explains the spectrum of extant ALMS1 variation. These results provide new insights into the evolutionary history of ALMS1 in humans and suggest that selective events identified in genomewide scans may be more complex than currently appreciated.

Pharmacogenomic (PGx) Counseling: Exploring Participant Questions about PGx Test Results.

As pharmacogenomic (PGx) use in healthcare increases, a better understanding of patient needs will be necessary to guide PGx result delivery. The Coriell Personalized Medicine Collaborative (CPMC) is a prospective study investigating the utility of personalized medicine. Participants received online genetic risk reports for 27 potentially actionable complex diseases and 7 drug-gene pairs and could request free, telephone-based genetic counseling (GC). To explore the needs of individuals receiving PGx results, we conducted a retrospective qualitative review of inquiries from CPMC participants who requested counseling from March 2009 to February 2017. Eighty out of 690 (12%) total GC inquiries were focused on the discussion of PGx results, and six salient themes emerged: "general help", "issues with drugs", "relevant disease experience", "what do I do now?", "sharing results", and "other drugs". The number of reported medications with a corresponding PGx result and participant engagement were significantly associated with PGx GC requests (p < 0.01 and p < 0.02, respectively). Our work illustrates a range of questions raised by study participants receiving PGx test results, most of which were addressed by a genetic counselor with few requiring referrals to prescribing providers or pharmacists. These results further support a role for genetic counselors in the team-based approach to optimal PGx result delivery.

rs11670527 Upstream of ZNF264 Associated with Body Mass Index in the Coriell Personalized Medicine Collaborative.

INTRODUCTION: the effects of obesity on health are a concern for the military as they affect the fitness to serve of active service members, increase costs to the Military Health System, and reduce quality of life for veterans and beneficiaries. Although obesity can be influenced by behavioral and environmental factors, it has also been shown to be associated with genetic risk factors that are not fully understood. MATERIALS AND METHODS: we performed a genome-wide association study of 5,251 participants in the Coriell Personalized Medicine Collaborative, which includes 2,111 Air Force participants. We applied a generalized linear model, using principal component analysis to account for population structure, and analyzed single-variant associations with body mass index (BMI) as a continuous variable, using a Bonferroni-corrected P-value threshold to account for multiplicity. RESULTS: we identified one genome-wide significant locus, rs11670527, upstream of the ZNF264 gene on chromosome 19, associated with BMI. CONCLUSIONS: the finding of an association between rs11670527 and BMI adds to the growing body of literature characterizing the complex genetics of obesity. These efforts may eventually inform personalized interventions aimed at achieving and maintaining healthy weight.

Habitat-Specific Clock Variation and Its Consequence on Reproductive Fitness.

The circadian clock controls daily activities at the cellular and organismic level, allowing an organism to anticipate incoming stresses and to use resources accordingly. The circadian clock has therefore been considered a fitness trait in multiple organisms. However, the mechanism of how circadian clock variation influences organismal reproductive fitness is still not well understood. Here we describe habitat-specific clock variation (HSCV) of asexual reproduction in Neurospora discreta, a species that is adapted to 2 different habitats, under or above tree bark. African (AF) N. discreta strains, whose habitat is above the tree bark in light-dark (LD) conditions, display a higher rhythmicity index compared with North American (NA) strains, whose habitat is under the tree bark in constant dark (DD). Although AF-type strains demonstrated an overall fitness advantage under LD and DD conditions, NA-type strains exhibit a habitat-specific fitness advantage in DD over the LD condition. In addition, we show that allelic variation of the clock-controlled gene, Ubiquinol cytochrome c oxidoreductase (NEUDI_158280), plays a role in HSCV by modulating cellular reactive oxygen species levels. Our results demonstrate a mechanism by which local adaptation involving circadian clock regulation influences reproductive fitness.