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

Història natural de les malalties genètiques mendelianes i complexes

Las enfermedades genéticas se clasifican típicamente en dos grandes grupos: las enfermedades mendelianas y las enfermedades complejas. Mientras que las enfermedades mendelianas se caracterizan por ser de baja frecuencia en la población y estar causadas por mutaciones en un gen particular, las enfermedades complejas son el principal problema sanitario en los países desarrollados y se encuentran producidas por la interacción de factores ambientales y factores genéticos. En este caso no se puede hablar de mutación en un determinado gen, sino de polimorfismo que incrementa en una pequeña fracción el riesgo a padecer la enfermedad. En la presente tesis se ha estudiado la distribución espacial de la variabilidad genética tanto en enfermedades mendelianas (en concreto la fibrosis quística, la fenilcetonuria y la b-talasemia) como en una enfermedad compleja (la enfermedad coronaria) en poblaciones europeas y de todo el mundo. Los resultados obtenidos sugieren que la distribución geográfica de la variabilidad genética de las enfermedades mendelianas depende principalmente de factores demográficos y de la historia de las poblaciones. Ahora bien, este efecto no es independiente de factores selectivos. En particular, fenómenos de selección equilibradora pueden incrementar o disminuir la variabilidad genética en una población dependiendo de el momento en el que se dio el evento selectivo. En el caso de la enfermedad compleja estudiada, la enfermedad coronaria, nuestros resultados indican que la distribución espacial de los polimorfismos de riesgo en poblaciones europeas depende, al igual que sucede con otros marcadores genéticos, principalmente de la historia de poblaciones, especialmente del poblamiento del continente europeo, la posterior reexpansión después del último periodo glacial y de las gran expansión poblacional de los agricultores durante el neolítico

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

Proportioning whole-genome single-nucleotide-polymorphism diversity for the identification of geographic population structure and genetic ancestry

The identification of geographic population structure and genetic ancestry on the basis of a minimal set of genetic markers is desirable for a wide range of applications in medical and forensic sciences. However, the absence of sharp discontinuities in the neutral genetic diversity among human populations implies that, in practice, a large number of neutral markers will be required to identify the genetic ancestry of one individual. We showed that it is possible to reduce the amount of markers required for detecting continental population structure to only 10 single-nucleotide polymorphisms (SNPs), by applying a newly developed ascertainment algorithm to Affymetrix GeneChip Mapping 10K SNP array data that we obtained from samples of globally dispersed human individuals (the Y Chromosome Consortium panel). Furthermore, this set of SNPs was able to recover the genetic ancestry of individuals from all four continents represented in the original data set when applied to an independent, much larger, worldwide population data set (Centre d'Etude du Polymorphisme Humain-Human Genome Diversity Project Cell Line Panel). Finally, we provide evidence that the unusual patterns of genetic variation we observed at the respective genomic regions surrounding the five most informative SNPs is in agreement with local positive selection being the explanation for the striking SNP allele-frequency differences we found between continental groups of human populations

PHOX2B polyalanine repeat length is associated with sudden infant death syndrome and unclassified sudden infant death in the Dutch population

Unclassified sudden infant death (USID) is the sudden and unexpected death of an infant that remains unexplained after thorough case investigation including performance of a complete autopsy and review of the circumstances of death and the clinical history. When the infant is below 1 year of age and with onset of the fatal episode apparently occurring during sleep, this is referred to as sudden infant death syndrome (SIDS). USID and SIDS remain poorly understood despite the identification of several environmental and some genetic risk factors. In this study, we investigated genetic risk factors involved in the autonomous nervous system in 195 Dutch USID/SIDS cases and 846 Dutch, age-matched healthy controls. Twenty-five DNA variants from 11 genes previously implicated in the serotonin household or in the congenital central hypoventilation syndrome, of which some have been associated with SIDS before, were tested. Of all DNA variants considered, only the length variation of the polyalanine repeat in exon 3 of the PHOX2B gene was found to be statistically significantly associated with USID/SIDS in the Dutch population after multiple test correction. Interestingly, our data suggest that contraction of the PHOX2B exon 3 polyalanine repeat that we found in six of 160 SIDS and USID cases and in six of 814 controls serves as a probable genetic risk factor for USID/SIDS at least in the Dutch population. Future studies are needed to confirm this finding and to understand the functional effect of the polyalanine repeat length variation, in particular contraction, in exon 3 of the PHOX2B gene

Clinal distribution of human genomic diversity across the Netherlands despite archaeological evidence for genetic discontinuities in Dutch population history

Background: The presence of a southeast to northwest gradient across Europe in human genetic diversity is a well-established observation and has recently been confirmed by genome-wide single nucleotide polymorphism (SNP) data. This pattern is traditionally explained by major prehistoric human migration events in Palaeolithic and Neolithic times. Here, we investigate whether (similar) spatial patterns in human genomic diversity also occur on a micro-geographic scale within Europe, such as in the Netherlands, and if so, whether these patterns could also be explained by more recent demographic events, such as those that occurred in Dutch population history.Methods: We newly collected data on a total of 999 Dutch individuals sampled at 54 sites across the country at 443,816 autosomal SNPs using the Genome-Wide Human SNP Array 5.0 (Affymetrix). We studied the individual genetic relationships by means of classical multidimensional scaling (MDS) using different genetic distance matrices, spatial ancestry analysis (SPA), and ADMIXTURE software. We further performed dedicated analyses to search for spatial patterns in the genomic variation and conducted simulations (SPLATCHE2) to provide a historical interpretation of the observed spatial patterns.Results: We detected a subtle but clearly noticeable genomic population substructure in the Dutch population, allowing differentiation of a north-eastern, central-western, central-northern and a southern group. Furthermore, we observed a statistically significant southeast to northwest cline in the distribution of genomic diversity across the Netherlands, similar to earlier findings from across Europe. Simulation analyses indicate that this genomic gradient could similarly be caused by ancient as well as by the more recent events in Dutch history.Conclusions: Considering the strong archaeological evidence for genetic discontinuity in the Netherlands, we interpret the observed clinal pattern of genomic diversity as being caused by recent rather than ancient events in Dutch population history. We therefore suggest that future human population genetic studies pay more attention to recent demographic history in interpreting genetic clines. Furthermore, our study demonstrates that genetic population substructure is detectable on a small geographic scale in Europe despite recent demographic events, a finding we consider potentially relevant for future epidemiological and forensic studies

Melanesian and Asian origins of Polynesians: mtDNA and Y chromosome gradients across the Pacific

The human settlement of the Pacific Islands represents one of the most recent major migration events of mankind. Polynesians originated in Asia according to linguistic evidence or in Melanesia according to archaeological evidence. To shed light on the genetic origins of Polynesians, we investigated over 400 Polynesians from 8 island groups, in comparison with over 900 individuals from potential parental populations of Melanesia, Southeast and East Asia, and Australia, by means of Y chromosome (NRY) and mitochondrial DNA (mtDNA) markers. Overall, we classified 94.1% of Polynesian Y chromosomes and 99.8% of Polynesian mtDNAs as of either Melanesian (NRY-DNA: 65.8%, mtDNA: 6%) or Asian (NRY-DNA: 28.3%, mtDNA: 93.8%) origin, suggesting a dual genetic origin of Polynesians in agreement with the "Slow Boat" hypothesis. Our data suggest a pronounced admixture bias in Polynesians toward more Melanesian men than women, perhaps as a result of matrilocal residence in the ancestral Polynesian society. Although dating methods are consistent with somewhat similar entries of NRY/mtDNA haplogroups into Polynesia, haplotype sharing suggests an earlier appearance of Melanesian haplogroups than those from Asia. Surprisingly, we identified gradients in the frequency distribution of some NRY/mtDNA haplogroups across Polynesia and a gradual west-to-east decrease of overall NRY/mtDNA diversity, not only providing evidence for a west-to-east direction of Polynesian settlements but also suggesting that Pacific voyaging was regular rather than haphazard. We also demonstrate that Fiji played a pivotal role in the history of Polynesia: humans probably first migrated to Fiji, and subsequent settlement of Polynesia probably came from Fiji

BMD loci contribute to ethnic and developmental differences in skeletal fragility across populations: Assessment of evolutionary selection pressures

Bone mineral density (BMD) is a highly heritable trait used both for the diagnosis of osteoporosis in adults and to assess bone health in children. Ethnic differences in BMD have been documented, with markedly higher levels in individuals of African descent, which partially explain disparity in osteoporosis risk across populations. To date, 63 independent genetic variants have been associated with BMD in adults of Northern-European ancestry. Here, we demonstrate that at least 61 of these variants are predictive of BMD early in life by studying their compound effect within two multiethnic pediatric cohorts. Furthermore, we show that within these cohorts and across populations worldwide the frequency of those alleles associated with increased BMD is systematically elevated in individuals of Sub-Saharan African ancestry. The amount of differentiation in the BMD genetic scores among Sub-Saharan and non-Sub-Saharan populations together with neutrality tests, suggest that these allelic differences are compatible with the hypothesis of selective pressures acting on the genetic determinants of BMD. These findings constitute an explorative contribution to the role of selection on ethnic BMD differences and likely a new example of polygenic adaptation acting on a human trait