Human genome diversity: frequently asked questions

Genetic heritage of the Baphuthi highlights an over-ethnicized notion of "Bushman" in the Maloti-Drakensberg, southern Africa

Using contemporary people as proxies for ancient communities is a contentious but necessary practice in anthropology. In southern Africa, the distinction between the Cape KhoeSan and eastern KhoeSan remains unclear, as ethnicity labels have been changed through time and most communities were decimated if not extirpated. The eastern KhoeSan may have had genetic distinctions from neighboring communities who speak Bantu languages and KhoeSan further away; alternatively, the identity may not have been tied to any notion of biology, instead denoting communities with a nomadic "lifeway" distinct from African agro-pastoralism. The Baphuthi of the 1800s in the Maloti-Drakensberg, southern Africa had a substantial KhoeSan constituency and a lifeway of nomadism, cattle raiding, and horticulture. Baphuthi heritage could provide insights into the history of the eastern KhoeSan. We examine genetic affinities of 23 Baphuthi to discern whether the narrative of KhoeSan descent reflects distinct genetic ancestry. Genome-wide SNP data (Illumina GSA) were merged with 52 global populations, for 160,000 SNPs. Genetic analyses show no support for a unique eastern KhoeSan ancestry distinct from other KhoeSan or southern Bantu speakers. The Baphuthi have strong affinities with early-arriving southern Bantu-speaking (Nguni) communities, as the later-arriving non-Nguni show strong evidence of recent African admixture possibly related to late-Iron Age migrations. The references to communities as "San" and "Bushman" in historic literature has often been misconstrued as notions of ethnic/biological distinctions. The terms may have reflected ambiguous references to non-sedentary polities instead, as seems to be the case for the eastern "Bushman" heritage of the Baphuthi.

The 'ethic of knowledge' and responsible science: Responses to genetically motivated racism

This study takes off from the ethical problem that racism grounded in population genetics raises. It is an analysis of four standard scientific responses to the problem of genetically motivated racism, seen in connection with the Human Genome Diversity Project (HGDP): (1) Discriminatory uses of scientific facts and arguments are in principle 'misuses' of scientific data that the researcher cannot be further responsible for. (2) In a strict scientific sense, genomic facts 'disclaim racism', which means that an epistemically correct grasp of genomics should be ethically justified. (3) Ethical difficulties are issues to be 'resolved' by an ethics institution or committee, which will guarantee the ethical quality of the research scrutinized. (4) Although population genetics occasionally may lead to racism, its overall 'value' for humankind justifies its cause as a desirable pursuit. I argue that these typical responses to genetically motivated racism supervene on a principle called the 'ethic of knowledge', which implies that an epistemically correct account has intrinsic ethical value. This principle, and its logically related ideas concerning the ethic of science, effectively avoids a deeper ethical question of responsibility in science from being raised.

Race-Specific Genetic Profiles of Homologous Recombination Deficiency in Multiple Cancers

Homologous recombination deficiency (HRD) has been used to predict both cancer prognosis and the response to DNA-damaging therapies in many cancer types. HRD has diverse manifestations in different cancers and even in different populations. Many screening strategies have been designed for detecting the sensitivity of a patient's HRD status to targeted therapies. However, these approaches suffer from low sensitivity, and are not specific to each cancer type and population group. Therefore, identifying race-specific and targetable HRD-related genes is of clinical importance. Here, we conducted analyses using genomic sequencing data that was generated by the Pan-Cancer Atlas. Collapsing non-synonymous variants with functional damage to HRD-related genes, we analyzed the association between these genes and race within cancer types using the optimal sequencing kernel association test (SKAT-O). We have identified race-specific mutational patterns of curated HRD-related genes across cancers. Overall, more significant mutation sites were found in ATM, BRCA2, POLE, and TOP2B in both the 'White' and 'Asian' populations, whereas PTEN, EGFG, and RIF1 mutations were observed in both the 'White' and 'African American/Black' populations. Furthermore, supported by pathogenic tendency databases and previous reports, in the 'African American/Black' population, several associations, including BLM with breast invasive carcinoma, ERCC5 with ovarian serous cystadenocarcinoma, as well as PTEN with stomach adenocarcinoma, were newly described here. Although several HRD-related genes are common across cancers, many of them were found to be specific to race. Further studies, using a larger cohort of diverse populations, are necessary to identify HRD-related genes that are specific to race, for guiding gene testing methods.

The vocal tract as a time machine: inferences about past speech and language from the anatomy of the speech organs

While speech and language do not fossilize, they still leave traces that can be extracted and interpreted. Here, we suggest that the shape of the hard structures of the vocal tract may also allow inferences about the speech of long-gone humans. These build on recent experimental and modelling studies, showing that there is extensive variation between individuals in the precise shape of the vocal tract, and that this variation affects speech and language. In particular, we show that detailed anatomical information concerning two components of the vocal tract (the lower jaw and the hard palate) can be extracted and digitized from the osteological remains of three historical populations from The Netherlands, and can be used to conduct three-dimensional biomechanical simulations of vowel production. We could recover the signatures of inter-individual variation between these vowels, in acoustics and articulation. While 'proof-of-concept', this study suggests that older and less well-preserved remains could be used to draw inferences about historic and prehistoric languages. Moreover, it forces us to clarify the meaning and use of the uniformitarian principle in linguistics, and to consider the wider context of language use, including the anatomy, physiology and cognition of the speakers. This article is part of the theme issue 'Reconstructing prehistoric languages'.

Genotype imputation and reference panel: a systematic evaluation on haplotype size and diversity

Here, 622 imputations were conducted with 394 customized reference panels for Han Chinese and European populations. Besides validating the fact that imputation accuracy could always benefit from the increased panel size when the reference panel was population specific, the results brought two new thoughts. First, when the haplotype size of the reference panel was fixed, the imputation accuracy of common and low-frequency variants (Minor Allele Frequency (MAF) > 0.5%) decreased while the population diversity of the reference panel increased, but for rare variants (MAF < 0.5%), a small fraction of diversity in panel could improve imputation accuracy. Second, when the haplotype size of the reference panel was increased with extra population-diverse samples, the imputation accuracy of common variants (MAF > 5%) for the European population could always benefit from the expanding sample size. However, for the Han Chinese population, the accuracy of all imputed variants reached the highest when reference panel contained a fraction of an extra diverse sample (8-21%). In addition, we evaluated the imputation performances in the existing reference panels, such as the Haplotype Reference Consortium (HRC), 1000 Genomes Project Phase 3 and the China, Oxford and Virginia Commonwealth University Experimental Research on Genetic Epidemiology (CONVERGE). For the European population, the HRC panel showed the best performance in our analysis. For the Han Chinese population, we proposed an optimum imputation reference panel constituent ratio if researchers would like to customize their own sequenced reference panel, but a high-quality and large-scale Chinese reference panel was still needed. Our findings could be generalized to the other populations with conservative genome; a tool was provided to investigate other populations of interest (https://github.com/Abyss-bai/reference-panel-reconstruction).

Tracing the genetic legacy in the French Caribbean islands: A study of mitochondrial and Y-chromosome lineages in the Guadeloupe archipelago

OBJECTIVES

The history of the Caribbean region is marked by numerous and various successive migration waves that resulted in a global blending of African, European, and Amerindian lineages. As the origin and genetic composition of the current population of French Caribbean islands has not been studied to date, we used both mitochondrial DNA and Y-chromosome markers to complete the characterization of the dynamics of admixture in the Guadeloupe archipelago.

MATERIALS AND METHODS

We sequenced the mitochondrial hypervariable regions and genotyped mitochondrial and Y-chromosomal single nucleotide polymorphisms (SNPs) of 198 individuals from five localities of the Guadeloupe archipelago.

RESULTS

The maternal haplogroups revealed a blend of 85% African lineages (mainly traced to Western, West-Central, and South-Eastern Africa), 12.5% Eurasian lineages, and 0.5% Amerindian lineages. We highlighted disequilibria between European paternal contribution (44%) and European maternal contribution (7%), pointing out an important sexual asymmetry. Finally, the estimated Native American component was strikingly low and supported the near-extinction of native lineages in the region.

DISCUSSION

We confirmed that all historically known migratory events indeed left a visible genetic imprint in the contemporary Caribbean populations. The data gathered clearly demonstrated the significant impact of the transatlantic slave trade on the Guadeloupean population's constitution. Altogether, the data in our study confirm that in the Caribbean region, human population variation is correlated with colonial and postcolonial policies and unique island histories.

Genomic Ancestry, CYP2D6, CYP2C9, and CYP2C19 Among Latin Americans

We present the distribution of CYP2D6, CYP2C9, and CYP2C19 variants and predicted phenotypes in 33 native and admixed populations from Ibero-America (n > 6,000) in the context of genetic ancestry (n = 3,387). Continental ancestries are the major determinants of frequencies of the increased-activity allele CYP2C19*17 and CYP2C19 gUMs (negatively associated with Native American ancestry), decreased-activity alleles CYP2D6*41 and CYP2C9*2 (positively associated with European ancestry), and decreased-activity alleles CYP2D6*17 and CYP2D6*29 (positively associated with African ancestry). For the rare alleles, CYP2C9*2 and CYPC19*17, European admixture accounts for their presence in Native American populations, but rare alleles CYP2D6*5 (null-activity), CYP2D6-multiplication alleles (increased activity), and CYP2C9*3 (decreased-activity) were present in the pre-Columbian Americas. The study of a broad spectrum of Native American populations from different ethno-linguistic groups show how autochthonous diversity shaped the distribution of pharmaco-alleles and give insights on the prevalence of clinically relevant phenotypes associated with drugs, such as paroxetine, tamoxifen, warfarin, and clopidogrel.

Counting Oceanians of Non-European, Non-Asian Descent (ONENA) in the South Pacific to Make Them Count in Global Health

Several diseases and vulnerabilities associated with genetic or microbial factors are more frequent among populations of Oceanian, Non-European, Non-Asian descent (ONENA). ONENA are specific and have long been isolated geographically. To our knowledge, there are no published official, quantitative, aggregated data on the populations impacted by these excess vulnerabilities in Oceania. We searched official census reports for updated estimates of the total population for each of the Pacific Island Countries and Territories (including Australia) and the US State of Hawaii, privileging local official statistical or censual sources. We multiplied the most recent total population estimate by the cumulative percentage of the ONENA population as determined in official reports. Including Australia and the US State of Hawaii, Oceania counts 27 countries and territories, populated in 2016 by approximately 41 M inhabitants (17 M not counting Australia) among which approximately 12.5 M (11.6 M not counting Australia) consider themselves of entire or partial ONENA ancestry. Specific genetic and microbiome traits of ONENA may be unique and need further investigation to adjust risk estimates, risk prevention, diagnostic and therapeutic strategies, to the benefit of populations in the Pacific and beyond.

Humans and Chimpanzees Display Opposite Patterns of Diversity in Arylamine N-Acetyltransferase Genes

Among the many genes involved in the metabolism of therapeutic drugs, human arylamine N-acetyltransferases (NATs) genes have been extensively studied, due to their medical importance both in pharmacogenetics and disease epidemiology. One member of this small gene family, NAT2, is established as the locus of the classic human acetylation polymorphism in drug metabolism. Current hypotheses hold that selective processes favoring haplotypes conferring lower NAT2 activity have been operating in modern humans’ recent history as an adaptation to local chemical and dietary environments. To shed new light on such hypotheses, we investigated the genetic diversity of the three members of the NAT gene family in seven hominid species, including modern humans, Neanderthals and Denisovans. Little polymorphism sharing was found among hominids, yet all species displayed high NAT diversity, but distributed in an opposite fashion in chimpanzees and bonobos (Pan genus) compared to modern humans, with higher diversity in Pan species at NAT1 and lower at NAT2, while the reverse is observed in humans. This pattern was also reflected in the results returned by selective neutrality tests, which suggest, in agreement with the predicted functional impact of mutations detected in non-human primates, stronger directional selection, presumably purifying selection, at NAT1 in modern humans, and at NAT2 in chimpanzees. Overall, the results point to the evolution of divergent functions of these highly homologous genes in the different primate species, possibly related to their specific chemical/dietary environment (exposome) and we hypothesize that this is likely linked to the emergence of controlled fire use in the human lineage.

The Impact of Population Variation in the Analysis of microRNA Target Sites

The impact of population variation in the analysis of regulatory interactions is an underdeveloped area. MicroRNA target recognition occurs via pairwise complementarity. Consequently, a number of computational prediction tools have been developed to identify potential target sites that can be further validated experimentally. However, as microRNA target predictions are done mostly considering a reference genome sequence, target sites showing variation among populations are neglected. Here, we studied the variation at microRNA target sites in human populations and quantified their impact in microRNA target prediction. We found that African populations carry a significant number of potential microRNA target sites that are not detectable in the current human reference genome sequence. Some of these targets are conserved in primates and only lost in Out-of-Africa populations. Indeed, we identified experimentally validated microRNA/transcript interactions that are not detected in standard microRNA target prediction programs, yet they have segregating target alleles abundant in non-European populations. In conclusion, we show that ignoring population diversity may leave out regulatory elements essential to understand disease and gene expression, particularly neglecting populations of African origin.

Effect of prothrombotic genotypes on the risk of venous thromboembolism in patients with and without ischemic stroke. The Tromsø Study

Essentials Prothrombotic genotypes may agument the risk of venous thromboembolism (VTE) after ischemic stroke. We studied this effect in a case-cohort study using a genetic risk score. In stroke patients, a one-category increase in the genetic risk score was associated with a 50% higher relative risk of VTE. The risk of VTE in stroke patients increased with an increasing number of risk alleles. SUMMARY: Background Patients with ischemic stroke have a transiently increased risk of subsequent venous thromboembolism (VTE). Prothrombotic genotypes may augment VTE risk under conditions of high thrombosis risk related to stroke. Aims To investigate the effect of prothrombotic genotypes in patients with ischemic stroke on the risk of VTE in a population-based case-cohort study. Methods Cases with incident VTE (n = 664) and a randomly selected age-weighted subcohort (n = 1817) were sampled from three surveys of the Tromsø Study (1994-2008). Participants were genotyped for ABO (rs8176719), F5 (rs6025), F2 (rs1799963), FGG (rs2066865) and F11 (rs2036914) single-nucleotide polymorphisms (SNPs). Cox regression models were used to calculate hazard ratios (HRs) for incident VTE according to individual SNPs and categories of risk alleles (5-SNP score; 0-1, 2, 3-4 and ≥ 5) in participants with and without ischemic stroke. Results There were 192 patients with incident stroke, of whom 43 developed VTE during a median of 15.2 years of follow-up. The risk alleles of individual SNPs augmented the elevated VTE risk brought about by ischemic stroke. In stroke patients, a one-category increase in the genetic risk score was associated with a 50% higher relative risk of overall VTE (HR 1.5, 95% confidence interval [CI] 1.3-1.8) and an 80% higher relative risk of provoked VTE (HR 1.8, 95% CI 1.5-2.1). Stroke patients with ≥ 5 risk alleles had a 12-fold (HR 11.7, 95% CI 4.1-33.3) higher relative risk of VTE than stroke-free participants with 0-1 risk alleles. Conclusions Prothrombotic genotypes increased the risk of VTE in stroke patients, and the risk increased with an increasing number of risk alleles.

Genes, Race, and Causation: US Public Perspectives About Racial Difference

Concerns have been raised that the increase in popular interest in genetics may herald a new era within which racial inequities are seen as 'natural' or immutable. In the following study, we provide data from a nationally representative survey on how the US population perceives general ability, athleticism, and intellect being determined by race and/or genetics and whether they believe racial health inequities to be primarily the product of genetic or social factors. We find that self-described race is of primary importance in attributing general ability to race, increasing age is a significant factor in attributing athleticism and intellect to genes and race, and education is a significant factor in decreasing such racially and genetically deterministic views . Beliefs about the meaning of race are statistically significantly associated with respect to the perception of athletic abilities and marginally associated with the perception of racial health inequalities being either socially or genetically derived. Race, education, socioeconomic status, and concepts of race were frequently found to be multiplicative in their statistical effects. The persistent acceptance of a genetically and racially deterministic view of athleticism among the White and older population group is discussed in respect to its social impact, as is the high level of agreement that general abilities are determined by race among non-White respondents and those of lower socioeconomic status. We argue that these findings highlight that both biological and non-biological forms of understanding race continue to play a role into the politics of race and social difference within contemporary US society.

Modelling human hard palate shape with Bézier curves

People vary at most levels, from the molecular to the cognitive, and the shape of the hard palate (the bony roof of the mouth) is no exception. The patterns of variation in the hard palate are important for the forensic sciences and (palaeo)anthropology, and might also play a role in speech production, both in pathological cases and normal variation. Here we describe a method based on Bézier curves, whose main aim is to generate possible shapes of the hard palate in humans for use in computer simulations of speech production and language evolution. Moreover, our method can also capture existing patterns of variation using few and easy-to-interpret parameters, and fits actual data obtained from MRI traces very well with as little as two or three free parameters. When compared to the widely-used Principal Component Analysis (PCA), our method fits actual data slightly worse for the same number of degrees of freedom. However, it is much better at generating new shapes without requiring a calibration sample, its parameters have clearer interpretations, and their ranges are grounded in geometrical considerations.

Validation of genotype imputation in Southeast Asian populations and the effect of single nucleotide polymorphism annotation on imputation outcome

Background

Imputation involves the inference of untyped single nucleotide polymorphisms (SNPs) in genome-wide association studies. The haplotypic reference of choice for imputation in Southeast Asian populations is unclear. Moreover, the influence of SNP annotation on imputation results has not been examined.

Methods

This study was divided into two parts. In the first part, we applied imputation to genotyped SNPs from Southeast Asian populations from the Pan-Asian SNP database. Five percent of the total SNPs were removed. The remaining SNPs were applied to imputation with IMPUTE2. The imputed outcomes were verified with the removed SNPs. We compared imputation references from Chinese and Japanese haplotypes from the HapMap phase II (HMII) and the complete set of haplotypes from the 1000 Genomes Project (1000G). The second part was imputation accuracy and yield in Thai patient dataset. Half of the autosomal SNPs was removed to create Set 1. Another dataset, Set 2, was then created where we switched which half of the SNPs were removed. Both Set 1 and Set 2 were imputed with HMII to create a complete imputed SNPs dataset. The dataset was used to validate association testing, SNPs annotation and imputation outcome.

Results

The accuracy was highest for all populations when using the HMII reference, but at the cost of a lower yield. Thai genotypes showed the highest accuracy over other populations in both HMII and 1000G panels, although accuracy and yield varied across chromosomes. Imputation was tested in a clinical dataset to compare accuracy in gene-related regions, and coding regions were found to have a higher accuracy and yield.

Conclusions

This work provides the first evidence of imputation reference selection for Southeast Asian studies and highlights the effects of SNP locations respective to genes on imputation outcome. Researchers will need to consider the trade-off between accuracy and yield in future imputation studies.

Electronic supplementary material

The online version of this article (10.1186/s12881-018-0534-8) contains supplementary material, which is available to authorized users.

Infant mortality in the United States

The infant mortality rate (IMR) of 6.0 per 1000 live births in the United States in 2013 is nearly the highest among developed countries. Moreover, the IMR among blacks is >twice that among whites-11.11 versus 5.06 deaths per 1000 live births.This higher IMR and racial disparity in IMR is due to a higher preterm birth rate (11.4% of live births in 2013) and higher IMR among term infants. The United States also ranks near the bottom for maternal mortality and life expectancy among the developed nations-despite ranking highest in the proportion of gross national product spent on health care. This suggests that factors other than health care contribute to the higher IMR and racial disparity in IMR. One factor is disadvantaged socioeconomic status. All of the actionable determinates that negatively impact health-personal behavior, social factors, heath-care access and quality and the environment-disproportionately affect the poor. Addressing disadvantaged socioeconomic status by improving access to quality health care and increasing social expenditures would have the greatest impact on the USA's IMR and racial disparity in IMR.

Negative control of CSL gene transcription by stress/DNA damage response and p53

CSL is a key transcriptional repressor and mediator of Notch signaling. Despite wide interest in CSL, mechanisms responsible for its own regulation are little studied. CSL down-modulation in human dermal fibroblasts (HDFs) leads to conversion into cancer associated fibroblasts (CAF), promoting keratinocyte tumors. We show here that CSL transcript levels differ among HDF strains from different individuals, with negative correlation with genes involved in DNA damage/repair. CSL expression is negatively regulated by stress/DNA damage caused by UVA, Reactive Oxygen Species (ROS), smoke extract, and doxorubicin treatment. P53, a key effector of the DNA damage response, negatively controls CSL gene transcription, through suppression of CSL promoter activity and, indirectly, by increased p21 expression. CSL was previously shown to bind p53 suppressing its activity. The present findings indicate that p53, in turn, decreases CSL expression, which can serve to enhance p53 activity in acute DNA damage response of cells.

African Ancestry, Social Factors, and Hypertension Among Non-Hispanic Blacks in the Health and Retirement Study

The biomedical literature contains much speculation about possible genetic explanations for the large and persistent black-white disparities in hypertension, but profound social inequalities are also hypothesized to contribute to this outcome. Our goal is to evaluate whether socioeconomic status (SES) differences provide a plausible mechanism for associations between African ancestry and hypertension in a U.S. cohort of older non-Hispanic blacks. We included only non-Hispanic black participants (N = 998) from the Health and Retirement Study who provided genetic data. We estimated percent African ancestry based on 84,075 independent single nucleotide polymorphisms using ADMIXTURE V1.23, imposing K = 4 ancestral populations, and categorized into quartiles. Hypertension status was self-reported in the year 2000. We used linear probability models (adjusted for age, sex, and southern birth) to predict prevalent hypertension with African ancestry quartile, before and after accounting for a small set of SES measures. Respondents with the highest quartile of African ancestry had 8 percentage points' (RD = 0.081; 95% CI: -0.001, 0.164) higher prevalence of hypertension compared to the lowest quartile. Adjustment for childhood disadvantage, education, income, and wealth explained over one-third (RD = 0.050; 95% CI: -0.034, 0.135) of the disparity. Explanations for the residual disparity remain unspecified and may include other indicators of SES or diet, lifestyle, and psychosocial mechanisms.

Variation in NAT2 acetylation phenotypes is associated with differences in food-producing subsistence modes and ecoregions in Africa

Background

Dietary changes associated to shifts in subsistence strategies during human evolution may have induced new selective pressures on phenotypes, as currently held for lactase persistence. Similar hypotheses exist for arylamine N-acetyltransferase 2 (NAT2) mediated acetylation capacity, a well-known pharmacogenetic trait with wide inter-individual variation explained by polymorphisms in the NAT2 gene. The environmental causative factor (if any) driving its evolution is as yet unknown, but significant differences in prevalence of acetylation phenotypes are found between hunter-gatherer and food-producing populations, both in sub-Saharan Africa and worldwide, and between agriculturalists and pastoralists in Central Asia. These two subsistence strategies also prevail among sympatric populations of the African Sahel, but knowledge on NAT2 variation among African pastoral nomads was up to now very scarce. Here we addressed the hypothesis of different selective pressures associated to the agriculturalist or pastoralist lifestyles having acted on the evolution of NAT2 by sequencing the gene in 287 individuals from five pastoralist and one agriculturalist Sahelian populations.

Results

We show that the significant NAT2 genetic structure of African populations is mainly due to frequency differences of three major haplotypes, two of which are categorized as decreased function alleles (NAT2*5B and NAT2*6A), particularly common in populations living in arid environments, and one fast allele (NAT2*12A), more frequently detected in populations living in tropical humid environments. This genetic structure does associate more strongly with a classification of populations according to ecoregions than to subsistence strategies, mainly because most Sahelian and East African populations display little to no genetic differentiation between them, although both regions hold nomadic or semi-nomadic pastoralist and sedentary agriculturalist communities. Furthermore, we found significantly higher predicted proportions of slow acetylators in pastoralists than in agriculturalists, but also among food-producing populations living in the Sahelian and dry savanna zones than in those living in humid environments, irrespective of their mode of subsistence.

Conclusion

Our results suggest a possible independent influence of both the dietary habits associated with subsistence modes and the chemical environment associated with climatic zones and biomes on the evolution of NAT2 diversity in sub-Saharan African populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0543-6) contains supplementary material, which is available to authorized users.

Evaluating the Photoprotective Effects of Ochre on Human Skin by In Vivo SPF Assessment: Implications for Human Evolution, Adaptation and Dispersal

Archaeological indicators of cognitively modern behaviour become increasingly prevalent during the African Middle Stone Age (MSA). Although the exploitation of ochre is viewed as a key feature of the emergence of modern human behaviour, the uses to which ochre and ochre-based mixtures were put remain ambiguous. Here we present the results of an experimental study exploring the efficacy of ochre as a topical photoprotective compound. This is achieved through the in vivo calculation of the sun protection factor (SPF) values of ochre samples obtained from Ovahimba women (Kunene Region, Northern Namibia) and the Palaeozoic Bokkeveld Group deposits of the Cape Supergroup (Western Cape Province, South Africa). We employ visible spectroscopy, energy-dispersive X-ray fluorescence (ED-XRF), X-ray diffraction (XRD) and granulometric analyses to characterise ochre samples. The capacity of ochre to inhibit the susceptibility of humans to the harmful effects of exposure to ultraviolet radiation (UVR) is confirmed and the mechanisms implicated in the efficacy of ochre as a sunscreen identified. It is posited that the habitual application of ochre may have represented a crucial innovation for MSA humans by limiting the adverse effects of ultraviolet exposure. This may have facilitated the colonisation of geographic regions largely unfavourable to the constitutive skin colour of newly arriving populations.

Across language families: Genome diversity mirrors linguistic variation within Europe

Objectives: The notion that patterns of linguistic and biological variation may cast light on each other and on population histories dates back to Darwin's times; yet, turning this intuition into a proper research program has met with serious methodological difficulties, especially affecting language comparisons. This article takes advantage of two new tools of comparative linguistics: a refined list of Indo‐European cognate words, and a novel method of language comparison estimating linguistic diversity from a universal inventory of grammatical polymorphisms, and hence enabling comparison even across different families. We corroborated the method and used it to compare patterns of linguistic and genomic variation in Europe. Materials and Methods: Two sets of linguistic distances, lexical and syntactic, were inferred from these data and compared with measures of geographic and genomic distance through a series of matrix correlation tests. Linguistic and genomic trees were also estimated and compared. A method (Treemix) was used to infer migration episodes after the main population splits. Results: We observed significant correlations between genomic and linguistic diversity, the latter inferred from data on both Indo‐European and non‐Indo‐European languages. Contrary to previous observations, on the European scale, language proved a better predictor of genomic differences than geography. Inferred episodes of genetic admixture following the main population splits found convincing correlates also in the linguistic realm. Discussion: These results pave the ground for previously unfeasible cross‐disciplinary analyses at the worldwide scale, encompassing populations of distant language families. Am J Phys Anthropol 157:630–640, 2015. © 2015 Wiley Periodicals, Inc.

Detecting individual ancestry in the human genome

Detecting and quantifying the population substructure present in a sample of individuals are of main interest in the fields of genetic epidemiology, population genetics, and forensics among others. To date, several algorithms have been proposed for estimating the amount of genetic ancestry within an individual. In the present review, we introduce the most widely used methods in population genetics for detecting individual genetic ancestry. We further show, by means of simulations, the performance of popular algorithms for detecting individual ancestry in various controlled demographic scenarios. Finally, we provide some hints on how to interpret the results from these algorithms.

Women who deliver twins are more likely to smoke and have high frequencies of specific SNPs: Results from a sample of African-American women who delivered preterm, low birth weight babies

OBJECTIVES

We examine if there are genetic and environmental differences between mothers of singleton and multiple pregnancies in a sample of African-American mothers.

METHODS

We focus on genomic areas suggested to increase or decrease the odds of multiple pregnancies. We computed the odds ratio (OR) and the 95% confidence interval (CI) for each SNP unadjusted or adjusted with smoking. SNPs' allelic differences between mothers of multiple pregnancies and singletons were also tested using Fisher's exact test. We considered additive terms for the SNPs' genotypes, smoking, and a multiplicative interaction term of two selected SNPs' genotypes.

RESULTS

We found significant interactions between smoking and SNPs of the CYP19A, MDM4, MTHFR and TP53 genes which correlated with higher odds of twinning. We also found a significant interaction between SNPs at the TP53 (rs8079544) and MTHFR gene (rs4846049), where the interaction between the homozygotes (TT for rs8079544, GG for rs4846049) correlated with lowered odds of multiple pregnancy.

CONCLUSIONS

We provide a mechanistic explanation and preliminary evidence for previous reports that mothers of twins are more likely to have smoked, despite seemingly conflicting evidence for the fertility-reducing effects of nicotine. Nicotine, as an aromatase inhibitor, inhibits estrogen synthesis and may allow for greater production of gonadotropins. While smoking may have deleterious effects on fertility across many genotypes, in women of specific genotypes it may raise their odds of producing twins. TP53 involvement suggests the necessity of future work examining relationships between women who bear multiples and cancer risk.

Population-specific common SNPs reflect demographic histories and highlight regions of genomic plasticity with functional relevance

BACKGROUND

Population differentiation is the result of demographic and evolutionary forces. Whole genome datasets from the 1000 Genomes Project (October 2012) provide an unbiased view of genetic variation across populations from Europe, Asia, Africa and the Americas. Common population-specific SNPs (MAF > 0.05) reflect a deep history and may have important consequences for health and wellbeing. Their interpretation is contextualised by currently available genome data.

RESULTS

The identification of common population-specific (CPS) variants (SNPs and SSV) is influenced by admixture and the sample size under investigation. Nine of the populations in the 1000 Genomes Project (2 African, 2 Asian (including a merged Chinese group) and 5 European) revealed that the African populations (LWK and YRI), followed by the Japanese (JPT) have the highest number of CPS SNPs, in concordance with their histories and given the populations studied. Using two methods, sliding 50-SNP and 5-kb windows, the CPS SNPs showed distinct clustering across large genome segments and little overlap of clusters between populations. iHS enrichment score and the population branch statistic (PBS) analyses suggest that selective sweeps are unlikely to account for the clustering and population specificity. Of interest is the association of clusters close to recombination hotspots. Functional analysis of genes associated with the CPS SNPs revealed over-representation of genes in pathways associated with neuronal development, including axonal guidance signalling and CREB signalling in neurones.

CONCLUSIONS

Common population-specific SNPs are non-randomly distributed throughout the genome and are significantly associated with recombination hotspots. Since the variant alleles of most CPS SNPs are the derived allele, they likely arose in the specific population after a split from a common ancestor. Their proximity to genes involved in specific pathways, including neuronal development, suggests evolutionary plasticity of selected genomic regions. Contrary to expectation, selective sweeps did not play a large role in the persistence of population-specific variation. This suggests a stochastic process towards population-specific variation which reflects demographic histories and may have some interesting implications for health and susceptibility to disease.

Nine things to remember about human genome diversity

Understanding how and why humans are biologically different is indispensable to get oriented in the ever-growing body of genomic data. Here we discuss the evidence based on which we can confidently state that humans are the least genetically variable primate, both when individuals and when populations are compared, and that each individual genome can be regarded as a mosaic of fragments of different origins. Each population is somewhat different from any other population, and there are geographical patterns in that variation. These patterns clearly indicate an African origin for our species, and keep a record of the main demographic changes accompanying the peopling of the whole planet. However, only a minimal fraction of alleles, and a small fraction of combinations of alleles along the chromosome, is restricted to a single geographical region (and even less so to a single population), and diversity between members of the same population is very large. The small genomic differences between populations and the extensive allele sharing across continents explain why historical attempts to identify, once and for good, major biological groups in humans have always failed. Nevertheless, racial categorization is all but gone, especially in clinical studies. We argue that racial labels may not only obscure important differences between patients but also that they have become positively useless now that cheap and reliable methods for genotyping are making it possible to pursue the development of truly personalized medicine.

GAGA: a new algorithm for genomic inference of geographic ancestry reveals fine level population substructure in Europeans

Attempts to detect genetic population substructure in humans are troubled by the fact that the vast majority of the total amount of observed genetic variation is present within populations rather than between populations. Here we introduce a new algorithm for transforming a genetic distance matrix that reduces the within-population variation considerably. Extensive computer simulations revealed that the transformed matrix captured the genetic population differentiation better than the original one which was based on the T1 statistic. In an empirical genomic data set comprising 2,457 individuals from 23 different European subpopulations, the proportion of individuals that were determined as a genetic neighbour to another individual from the same sampling location increased from 25% with the original matrix to 52% with the transformed matrix. Similarly, the percentage of genetic variation explained between populations by means of Analysis of Molecular Variance (AMOVA) increased from 1.62% to 7.98%. Furthermore, the first two dimensions of a classical multidimensional scaling (MDS) using the transformed matrix explained 15% of the variance, compared to 0.7% obtained with the original matrix. Application of MDS with Mclust, SPA with Mclust, and GemTools algorithms to the same dataset also showed that the transformed matrix gave a better association of the genetic clusters with the sampling locations, and particularly so when it was used in the AMOVA framework with a genetic algorithm. Overall, the new matrix transformation introduced here substantially reduces the within population genetic differentiation, and can be broadly applied to methods such as AMOVA to enhance their sensitivity to reveal population substructure. We herewith provide a publically available (http://www.erasmusmc.nl/fmb/resources/GAGA) model-free method for improved genetic population substructure detection that can be applied to human as well as any other species data in future studies relevant to evolutionary biology, behavioural ecology, medicine, and forensics.

Evaluation of large-scale genetic structure in complex demographic and historical scenarios: the mitochondrial DNA and Y-chromosome pools of the Iberian Atlantic façade

Genetic structural patterns of human populations are usually a combination of long-term evolutionary forces and short-term social, cultural, and demographic processes. Recently, using mitochondrial DNA and Y-chromosome loci, various studies in northern Spain have found evidence that the geographical distribution of Iron Age tribal peoples might have influenced current patterns of genetic structuring in several autochthonous populations. Using the wealth of data that are currently available from the whole territory of the Iberian Peninsula, we have evaluated its genetic structuring in the spatial scale of the Atlantic façade. Hierarchical tree modeling procedures, combined with a classic analysis of molecular variance (AMOVA), were used to model known sociocultural divisions from the third century BCE to the eighth century CE, contrasting them with uniparental marker data. Our results show that, while mountainous and abrupt areas of the Iberian North bear the signals of long-term isolation in their maternal and paternal gene pools, the makeup of the Atlantic façade as a whole can be related to tribal population groups that predate the Roman conquest of the Peninsula. The maintenance through time of such a structure can be related to the numerous geographic barriers of the Iberian mainland, which have historically conditioned its settlement patterns and the occurrence of genetic drift processes.

The influence of clan structure on the genetic variation in a single Ghanaian village

Socioeconomic and cultural factors are thought to have an important role in influencing human population genetic structure. To explain such population structure differences, most studies analyse genetic differences among widely dispersed human populations. In contrast, we have studied the genetic structure of an ethnic group occupying a single village in north-eastern Ghana. We found a markedly skewed male population substructure because of an almost complete lack of male gene flow among Bimoba clans in this village. We also observed a deep male substructure within one of the clans in this village. Among all males, we observed only three Y-single-nucleotide polymorphism (SNP) haplogroups: E1b1a*-M2, E1b1a7a*-U174 and E1b1a8a*-U209, P277, P278. In contrast to the marked Y-chromosomal substructure, mitochondrial DNA HVS-1 sequence variation and autosomal short-tandem repeats variation patterns indicate high genetic diversities and a virtually random female-mediated gene flow among clans. On the extreme micro-geographical scale of this single Bimoba village, correspondence between the Y-chromosome lineages and clan membership could be due to the combined effects of the strict patrilocal and patrilineal structure. If translated to larger geographic scales, our results would imply that the extent of variation in uniparentally inherited genetic markers, which are typically associated with historical migration on a continental scale, could equally likely be the result of many small and different cumulative effects of social factors such as clan membership that act at a local scale. Such local scale effects should therefore be considered in genetic studies, especially those that use uniparental markers, before making inferences about human history at large.

Detecting genetic isolation in human populations: a study of European language minorities

The identification of isolation signatures is fundamental to better understand the genetic structure of human populations and to test the relations between cultural factors and genetic variation. However, with current approaches, it is not possible to distinguish between the consequences of long-term isolation and the effects of reduced sample size, selection and differential gene flow. To overcome these limitations, we have integrated the analysis of classical genetic diversity measures with a Bayesian method to estimate gene flow and have carried out simulations based on the coalescent. Combining these approaches, we first tested whether the relatively short history of cultural and geographical isolation of four “linguistic islands” of the Eastern Alps (Lessinia, Sauris, Sappada and Timau) had left detectable signatures in their genetic structure. We then compared our findings to previous studies of European population isolates. Finally, we explored the importance of demographic and cultural factors in shaping genetic diversity among the groups under study. A combination of small initial effective size and continued genetic isolation from surrounding populations seems to provide a coherent explanation for the diversity observed among Sauris, Sappada and Timau, which was found to be substantially greater than in other groups of European isolated populations. Simulations of micro-evolutionary scenarios indicate that ethnicity might have been important in increasing genetic diversity among these culturally related and spatially close populations.

Human loci involved in drug biotransformation: worldwide genetic variation, population structure, and pharmacogenetic implications

Understanding the role of inheritance in individual variation in drug response is the focus of pharmacogenetics (PGx). A key part of this understanding is quantifying the role of genetic ancestry in this phenotypic outcome. To provide insight into the relationship between ethnicity and drug response, this study first infers the global distribution of PGx variation and defines its structure. Second, the study evaluates if geographic population structure stems from all PGx loci in general, or if structure is caused by specific genes. Lastly, we identify the genetic variants contributing the greatest proportion of such structure. Our study describes the global genetic structure of PGx loci across the 52 populations of the Human Genome Diversity Cell-Line Panel, the most inclusive set of human populations freely available for studies on human genetic variation. By analysing genetic variation at 1,001 single nucleotide polymorphisms (SNPs) involved in biotransformation of exogenous substances, we describe the between-populations PGx variation, as well geographical groupings of diversity. In addition, with discriminant analysis of principal component (DAPC), we infer how many and which groups of populations are supported by PGx variation, and identify which SNPs actually contribute to the PGx structure between such groups. Our results show that intergenic, synonymous and non-synonymous SNPs show similar levels of genetic variation across the globe. Conversely, loci coding for Cytochrome P450s (mainly metabolizing exogenous substances) show significantly higher levels of genetic diversity between populations than the other gene categories. Overall, genetic variation at PGx loci correlates with geographic distances between populations, and the apportionment of genetic variation is similar to that observed for the rest of the genome. In other words, the pattern of PGx variation has been mainly shaped by the demographic history of our species, as in the case of most of our genes. The population structure defined by PGx loci supports the presence of six genetic clusters reflecting geographic location of samples. In particular, the results of the DAPC analyses show that 27 SNPs substantially contribute to the first three discriminant functions. Among these SNPs, some, such as the intronic rs1403527 of NR1I2 and the non-synonymous rs699 of AGT, are known to be associated with specific drug responses. Their substantial variation between different groups of populations may have important implications for PGx practical applications.

Human neutral genetic variation and forensic STR data

The forensic genetics field is generating extensive population data on polymorphism of short tandem repeats (STR) markers in globally distributed samples. In this study we explored and quantified the informative power of these datasets to address issues related to human evolution and diversity, by using two online resources: an allele frequency dataset representing 141 populations summing up to almost 26 thousand individuals; a genotype dataset consisting of 42 populations and more than 11 thousand individuals. We show that the genetic relationships between populations based on forensic STRs are best explained by geography, as observed when analysing other worldwide datasets generated specifically to study human diversity. However, the global level of genetic differentiation between populations (as measured by a fixation index) is about half the value estimated with those other datasets, which contain a much higher number of markers but much less individuals. We suggest that the main factor explaining this difference is an ascertainment bias in forensics data resulting from the choice of markers for individual identification. We show that this choice results in average low variance of heterozygosity across world regions, and hence in low differentiation among populations. Thus, the forensic genetic markers currently produced for the purpose of individual assignment and identification allow the detection of the patterns of neutral genetic structure that characterize the human population but they do underestimate the levels of this genetic structure compared to the datasets of STRs (or other kinds of markers) generated specifically to study the diversity of human populations.

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.

A population-genetic perspective on the similarities and differences among worldwide human populations

Recent studies have produced a variety of advances in the investigation of genetic similarities and differences among human populations. Here, I pose a series of questions about human population-genetic similarities and differences, and I then answer these questions by numerical computation with a single shared population-genetic data set. The collection of answers obtained provides an introductory perspective for understanding key results on the features of worldwide human genetic variation.

The impact of recent events on human genetic diversity

The historical record tells us stories of migrations, population expansions and colonization events in the last few thousand years, but what was their demographic impact? Genetics can throw light on this issue, and has mostly done so through the maternally inherited mitochondrial DNA (mtDNA) and the male-specific Y chromosome. However, there are a number of problems, including marker ascertainment bias, possible influences of natural selection, and the obscuring layers of the palimpsest of historical and prehistorical events. Y-chromosomal lineages are particularly affected by genetic drift, which can be accentuated by recent social selection. A diversity of approaches to expansions in Europe is yielding insights into the histories of Phoenicians, Roma, Anglo-Saxons and Vikings, and new methods for producing and analysing genome-wide data hold much promise. The field would benefit from more consensus on appropriate methods, and better communication between geneticists and experts in other disciplines, such as history, archaeology and linguistics.

Nothing in medicine makes sense, except in the light of evolution

The practice of medicine is a fruitful marriage of classic diagnostic and healing arts with modern advancements in many relevant sciences. The scientific aspects of medicine are rooted in understanding the biology of our species and those of other organisms that interact with us in health and disease. Thus, it is reasonable to paraphrase Dobzhansky, stating that, "nothing in the biological aspects of medicine makes sense except in the light of evolution." However, the art and science of medicine are also rooted in the unusual cognitive abilities of humans and the cultural evolutionary processes arising. This explains the rather bold and inclusive title of this essay. The near complete absence of evolution in medical school curricula is a historical anomaly that needs correction. Otherwise, we will continue to train generations of physicians who lack understanding of some fundamental principles that should guide both medical practice and research. I here recount my attempts to correct this deficiency at my own medical school and the lessons learned. I also attempt to summarize what I teach in the limited amount of time allowed for the purpose. Particular attention is given to the value of comparing human physiology and disease with those of other closely related species. There is a long way to go before the teaching of evolution can be placed in its rightful context within the medical curriculum. However, the trend is in the right direction. Let us aim for a day when an essay like this will no longer be relevant.

KIR genotypic diversity can track ancestries in heterogeneous populations: a potential confounder for disease association studies

Killer cell immunoglobulin-like receptors (KIR) are encoded by highly polymorphic genes that regulate the activation of natural killer (NK) cells and other lymphocyte subsets and likely play key roles in innate and adaptive immunity. Association studies increasingly implicate KIR in disease predisposition and outcome but could be confounded by unknown KIR genetic structure in heterogeneous populations. To examine this, we characterized the diversity of 16 KIR genes in 712 Northern Californians (NC) stratified by self-assigned ethnicities and compared the profiles of KIR polymorphism with other US and global populations using a reference database. Sixty-eight distinct KIR genotypes were characterized: 58 in 457 Caucasians (NCC), 17 in 47 African Americans (NCAA), 21 in 80 Asians (NCA), 20 in 74 Hispanics (NCH), and 18 in 54 "other" ethnicities (NCO). KIR genotype patterns and frequencies in the 4 defined ethnicities were compared with each other and with 34 global populations by phylogenetic analysis. Although there were no population-specific genotypes, the KIR genotype frequency patterns faithfully traced the ancestry of NCC, NCAA, and NCA but not of NCH whose ancestries are known to be more heterogeneous. KIR genotype frequencies can therefore track ethnic ancestries in modern urban populations. Our data emphasize the importance of selecting ethnically matched controls in KIR-based studies to avert spurious associations.

Arylamine N-acetyltransferase 2 (NAT2) genetic diversity and traditional subsistence: a worldwide population survey

Arylamine N-acetyltransferase 2 (NAT2) is involved in human physiological responses to a variety of xenobiotic compounds, including common therapeutic drugs and exogenous chemicals present in the diet and the environment. Many questions remain about the evolutionary mechanisms that have led to the high prevalence of slow acetylators in the human species. Evidence from recent surveys of NAT2 gene variation suggests that NAT2 slow-causing variants might have become targets of positive selection as a consequence of the shift in modes of subsistence and lifestyle in human populations in the last 10,000 years. We aimed to test more extensively the hypothesis that slow acetylation prevalence in humans is related to the subsistence strategy adopted by the past populations. To this end, published frequency data on the most relevant genetic variants of NAT2 were collected from 128 population samples (14,679 individuals) representing different subsistence modes and dietary habits, allowing a thorough analysis at both a worldwide and continent scale. A significantly higher prevalence of the slow acetylation phenotype was observed in populations practicing farming (45.4%) and herding (48.2%) as compared to populations mostly relying on hunting and gathering (22.4%) (P = 0.0007). This was closely mirrored by the frequency of the slow 590A variant that was found to occur at a three-fold higher frequency in food producers (25%) as compared to hunter-gatherers (8%). These findings are consistent with the hypothesis that the Neolithic transition to subsistence economies based on agricultural and pastoral resources modified the selective regime affecting the NAT2 acetylation pathway. Furthermore, the vast amount of data collected enabled us to provide a comprehensive and up-to-date description of NAT2 worldwide genetic diversity, thus building up a useful resource of frequency data for further studies interested in epidemiological or anthropological research questions involving NAT2.

Helicobacter pylori from Peruvian amerindians: traces of human migrations in strains from remote Amazon, and genome sequence of an Amerind strain

Background

The gastric pathogen Helicobacter pylori is extraordinary in its genetic diversity, the differences between strains from well-separated human populations, and the range of diseases that infection promotes.

Principal Findings

Housekeeping gene sequences from H. pylori from residents of an Amerindian village in the Peruvian Amazon, Shimaa, were related to, but not intermingled with, those from Asia. This suggests descent of Shimaa strains from H. pylori that had infected the people who migrated from Asia into The Americas some 15,000+ years ago. In contrast, European type sequences predominated in strains from Amerindian Lima shantytown residents, but with some 12% Amerindian or East Asian-like admixture, which indicates displacement of ancestral purely Amerindian strains by those of hybrid or European ancestry. The genome of one Shimaa village strain, Shi470, was sequenced completely. Its SNP pattern was more Asian- than European-like genome-wide, indicating a purely Amerind ancestry. Among its unusual features were two cagA virulence genes, each distinct from those known from elsewhere; and a novel allele of gene hp0519, whose encoded protein is postulated to interact with host tissue. More generally, however, the Shi470 genome is similar in gene content and organization to those of strains from industrialized countries.

Conclusions

Our data indicate that Shimaa village H. pylori descend from Asian strains brought to The Americas many millennia ago; and that Amerind strains are less fit than, and were substantially displaced by, hybrid or European strains in less isolated communities. Genome comparisons of H. pylori from Amerindian and other communities should help elucidate evolutionary forces that have shaped pathogen populations in The Americas and worldwide.