Pitfalls and challenges with population assignments of individuals from admixed populations: Applying Genogeographer on Brazilian individuals.

The assignment of individuals to a population can be of importance for the identification of mass disaster victims or criminal offenders in the field of forensic genetics. This assignment is based on biostatistical methods that process data of ancestry informative markers (AIMs), which are selected based on large allele frequency differences between the populations of interest. However, population assignments of individuals with an admixed genetic background are challenging. Admixed individuals are genetic mosaics of chromosomal segments from the parental populations, which may lead to ambiguous or no population assignment. This is problematic since admixture events are a substantial part of human history. In this study, we present challenges of interpreting the evidential weight of population assignments. We used Genogeographer for likelihood ratio (LR) calculations and Brazilians as examples of admixed individuals. Brazilians are a very heterogenous population representing a three-way admixture between Native Americans, Europeans, and Africans. Ancestry informative markers were typed in a total of 589 individuals from Brazil using the Precision ID Ancestry Panel. The Brazilians were assigned to six metapopulations (East Asia, Europe, Middle East, North Africa, South-Central Asia, Sub-Saharan Africa) defined in the Genogeographer software and LRs were calculated if the AIM profile was not an outlier in all metapopulations and simulated two-way (1:1) admixtures of the six metapopulations. Population assignments failed for 55% of the samples. These samples had significantly higher genetic contributions from East Asia, South-Central Asia and Sub-Saharan Africa, and significantly lower genetic contributions from Europe. Most of the individuals with population assignments were assigned to the metapopulations of Middle East (58%) or North Africa (36%), followed by Europe (4%), South-Central Asia (1%), and Sub-Saharan Africa (1%). For 8% of the samples, population assignments were only possible when assignments to simulated two-way (1:1) admixtures of the six metapopulations were considered. Most of these individuals were assigned to two-way admixtures of North Africa, South-Central Asia, or Sub-Saharan Africa. Relatively low median likelihood ratios (LRs<1000) were observed when comparing population likelihoods for Europe, Middle East, North Africa, South-Central Asia, or simulated 1:1 admixtures of these metapopulations. Comparisons including East Asian or Sub-Saharan African populations resulted in larger median LRs (LR>1010). The results suggested that the Precision ID Ancestry Panel provided too little information and that additional markers specifically selected for sub-continental differentiation may be required for accurate population assignment of admixed individuals. Furthermore, a Genogeographer database with additional populations including admixed populations would be advantageous for interpretation of admixed AIM profiles. It would likely increase the number of population assignments and illustrate alternatives to the most likely population, which would be valuable information for the case officer when writing the case report.

Skin pigmentation and genetic variants in an admixed Brazilian population of primarily European ancestry.

Although many genes have been shown to be associated with human pigmentary traits and forensic prediction assays exist (e.g. HIrisPlex-S), the genetic knowledge about skin colour remains incomplete. The highly admixed Brazilian population is an interesting study population for investigation of the complex genotype-phenotype architecture of human skin colour because of its large variation. Here, we compared variants in 22 pigmentary genes with quantitative skin pigmentation levels on the buttock, arm, and forehead areas of 266 genetically admixed Brazilian individuals. The genetic ancestry of each individual was estimated by typing 46 AIM-InDels. The mean proportion of genetic ancestry was 68.8% European, 20.8% Sub-Saharan African, and 10.4% Native American. A high correlation (adjusted R2 = 0.65, p < 0.05) was observed between nine SNPs and quantitative skin pigmentation using multiple linear regression analysis. The correlations were notably smaller between skin pigmentation and biogeographic ancestry (adjusted R2 = 0.45, p < 0.05), or markers in the leading forensic skin colour prediction system, the HIrisPlex-S (adjusted R2 = 0.54, p < 0.05). Four of the nine SNPs, OCA2 rs1448484 (rank 2), APBA2 rs4424881 (rank 4), MFSD12 rs10424065 (rank 8), and TYRP1 1408799 (rank 9) were not investigated as part of the HIrisPlex-S selection process, and therefore not included in the HIrisPlex-S model. Our results indicate that these SNPs account for a substantial part of the skin colour variation in individuals of admixed ancestry. Hence, we suggest that these SNPs are considered when developing future skin colour prediction models.

Analysis of ancestry informative markers in three main ethnic groups from Ecuador supports a trihybrid origin of Ecuadorians.

Ancestry inference is traditionally done using autosomal SNPs that present great allele frequency differences among populations from different geographic regions. These ancestry informative markers (AIMs) are useful for determining the most likely biogeographic ancestry or population of origin of an individual. Due to the growing interest in AIMs and their applicability in different fields, commercial companies have started to develop AIM multiplexes targeted for Massive Parallel Sequencing platforms. This project focused on the study of three main ethnic groups from Ecuador (Kichwa, Mestizo, and Afro-Ecuadorian) using the Precision ID Ancestry panel (Thermo Fisher Scientific). In total, 162 Ecuadorian individuals were investigated. The Afro-Ecuadorian and Mestizo showed higher average genetic diversities compared to the Kichwa. These results are consistent with the highly admixed nature of the first two groups. The Kichwa showed the highest proportion of Native Amerindian (NAM) ancestry relative to the other two groups. The Mestizo had an admixed ancestry of NAM and European with a larger European component, whereas the Afro-Ecuadorian were highly admixed presenting proportions of African, Native Amerindian, and European ancestries. The comparison of our results with previous studies based on uniparental markers (i.e. Y chromosome and mtDNA) highlighted the sex-biased admixture process in the Ecuadorian Mestizo. Overall, the data generated in this work represent one important step to assess the application of ancestry inference in admixed populations in a forensic context.

Typing of Amerindian Kichwas and Mestizos from Ecuador with the SNPforID multiplex.

A total of 119 unrelated individuals from two of the major ethnic groups in Ecuador were typed for 49 of the autosomal single nucleotide polymorphisms (SNPs) in the SNPforID 52plex using the SNapShot(®) assay. Of the above, 42 samples originated from Mestizos (an admixed population) and the remaining 77 were from Native Amerindian Kichwas. We obtained full SNP profiles in all individuals and concordance of duplicated analyses. No deviation from Hardy-Weinberg equilibrium (HWE) was observed for any SNP in the Mestizo and Kichwa populations and only one and four pairs of loci, respectively showed significant linkage disequilibrium. A relatively low genetic diversity and global positive F(IS) value was observed in Kichwas. A statistically significant global F(ST) value was obtained when the two Ecuadorian populations were compared with populations in Spain, Portugal, Argentina, Denmark, Greenland, China, Somalia and Mozambique. All pairwise F(ST) values were statistically significant. A multi-dimensional scaling based on pairwise F(ST) values showed that the Kichwa population differed from all other populations investigated and that the Mestizos had an intermediate position between Kichwas and Europeans. An admixture analysis indicated that the greater contributor to the Mestizo population was the Kichwas (71.2%) compared to the European contribution. The combined mean match probability and mean paternity exclusion probability were 3.3 × 10(-17) and 0.998, respectively, for the Mestizo population and 3.3 × 10(-14) and 0.993, respectively, for the Kichwa population.

Introduction of an single nucleodite polymorphism-based "Major Y-chromosome haplogroup typing kit" suitable for predicting the geographical origin of male lineages.

The European Consortium "High-throughput analysis of single nucleotide polymorphisms for the forensic identification of persons--SNPforID", has performed a selection of candidate Y-chromosome single nucleotide polymorphisms (SNPs) for making inferences on the geographic origin of an unknown sample. From more than 200 SNPs compiled in the phylogenetic tree published by the Y-Chromosome Consortium, and looking at the population studies previously published, a package of 29 SNPs has been selected for the identification of major population haplogroups. A "Major Y-chromosome haplogroup typing kit" has been developed, which allows the multiplex amplification of all 29 SNPs in a single reaction. Allele genotyping was performed with a single base extension reaction (minisequencing) detected by CE. The validation of the multiplex was performed in a total of 1126 unrelated males distributed among 12 worldwide populations. The approach takes advantage of the specific geographic distribution of the Y-chromosome haplogroups and demonstrates the utility of binary polymorphisms to infer the origin of a male lineage.