Evaluating a subset of ancestry informative SNPs for discriminating among Southwest Asian and circum-Mediterranean populations

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.

Insights into AIM-InDel diversities in Yunnan Miao and Hani ethnic groups of China for forensic and population genetic purposes

Background

Ancestry informative markers are regarded as useful tools for inferring the ancestral information of an individual, which have been widely used in the criminal investigations and population genetic studies. Previously, a multiplex amplification panel containing 39 AIM-InDel loci was constructed. This study aims to investigate the genetic polymorphisms of these 39 AIM-InDel loci in Yunnan Hani and Miao ethnic groups, and to uncover their genetic affinities with reference populations based on the AIM-InDel markers.

Materials and methods

In this research, 39 AIM-InDel profiles of 203 unrelated Miao individuals and 203 unrelated Hani individuals in Yunnan province of China were acquired. Additionally, we evaluated the genetic polymorphisms of 39 InDel loci in Yunnan Miao and Hani groups. Moreover, the genetic relationships among Yunnan Miao, Hani and reference populations were also clarified based on Nei’s genetic distances, pairwise fixation indexes, principal component analyses, phylogenetic analyses, and STRUCTURE analyses.

Results

Genetic diversity analyses demonstrated that these InDel loci showed varying degrees of genetic polymorphisms, and could be utilized in forensic identifications in Yunnan Miao and Hani groups. The results of principal component analyses, phylogenetic analyses and Structure analyses revealed that Yunnan Miao and Hani groups had closer genetic relationships with East Asian populations, especially with the populations from Southern China. This research enriched the genetic data of Chinese ethnic minority, and provided ancestral information of Yunnan Miao and Hani groups from the perspective of population genetics.

Supplementary Information

The online version contains supplementary material available at 10.1186/s41065-022-00238-9.

Evaluation of a custom QIAseq targeted DNA panel with 164 ancestry informative markers sequenced with the Illumina MiSeq

Introduction of new methods requires meticulous evaluation before they can be applied to forensic genetic case work. Here, a custom QIAseq Targeted DNA panel with 164 ancestry informative markers was assessed using the MiSeq sequencing platform. Concordance, sensitivity, and the capability for analysis of mixtures were tested. The assay gave reproducible and nearly concordant results with an input of 10 and 2 ng DNA. Lower DNA input led to an increase in both locus and allele drop-outs, and a higher variation in heterozygote balance. Locus or allele drop-outs in the samples with less than 2 ng DNA input were not necessarily associated with the overall performance of a locus. Thus, the QIAseq assay will be difficult to implement in a forensic genetic setting where the sample material is often scarce and of poor quality. With equal or near equal mixture ratios, the mixture DNA profiles were easily identified by an increased number of imbalanced heterozygotes. For more skewed mixture ratios, the mixture DNA profiles were identified by an increased noise level. Lastly, individuals from Great Britain and the Middle East were investigated. The Middle Eastern individuals showed a greater affinity with South European populations compared to North European populations.

Assessment of the effectiveness of the EUROFORGEN NAME and Precision ID Ancestry panel markers for ancestry investigations

The EUROFORGEN NAME panel is a regional ancestry panel designed to differentiate individuals from the Middle East, North Africa, and Europe. The first version of the panel was developed for the MassARRAY system and included 111 SNPs. Here, a custom AmpliSeq EUROFORGEN NAME panel with 102 of the original 111 loci was used to sequence 1098 individuals from 14 populations from Europe, the Middle East, North Africa, North-East Africa, and South-Central Asia. These samples were also sequenced with a global ancestry panel, the Precision ID Ancestry Panel. The GenoGeographer software was used to assign the AIM profiles to reference populations and calculate the weight of the evidence as likelihood ratios. The combination of the EUROFORGEN NAME and Precision ID Ancestry panels led to fewer ambiguous assignments, especially for individuals from the Middle East and South-Central Asia. The likelihood ratios showed that North African individuals could be separated from European and Middle Eastern individuals using the Precision ID Ancestry Panel. The separation improved with the addition of the EUROFORGEN NAME panel. The analyses also showed that the separation of Middle Eastern populations from European and South-Central Asian populations was challenging even when both panels were applied.

Reliability of phenotype estimation and extended classification of ancestry using decedent samples

The Illumina® MiSeq FGx™, in conjunction with the ForenSeq™ DNA Signature Prep kit, produces genotypes of the CODIS-required short tandem repeats and provides phenotype and biogeographical ancestry estimations via phenotype-informative and ancestry-informative markers, respectively. Although both markers have been validated for use in forensic biology, there is little data to determine the practical utility of these estimations to assist in identifying missing persons using decedent casework samples. The accuracy and utility of phenotypic and ancestral estimations were investigated for 300 samples received by the Los Angeles County Department of Medical Examiner-Coroner. piSNP genotypes were translated into hair and eye colors using the Forenseq™ Universal Analysis Software (UAS) on the MiSeq FGx™ and the HIrisPlex System, and statistical accuracy was evaluated in context with the reported decedent characteristics. Similarly, estimates of each decedent's biogeographical ancestry were compared to assess the efficacy of these markers to predict ancestry correctly. The average UAS and the HIrisPlex system prediction accuracy for brown and blue eyes were 95.3% and 96.2%, respectively. Intermediate eye color could not be predicted with high accuracy using either system. Other than the black hair phenotype reporting an accuracy that exceeded 90% using either system, hair color was also too variable to be predicted with high accuracy. The FROG-kb database distinguishes decedents adequately beyond the Asian, African, European, and Admixed American global ancestries provided by the MiSeq FGx™ UAS PCA plots. FROG-kb correctly identified Middle Eastern, Pacific Islander, Latin American, or Jewish ancestries with accuracies of 70.0%, 81.8%, 73.8%, and 86.7%, respectively.

Genetic relationships of Southwest Asian and Mediterranean populations

The Southwest Asian, circum-Mediterranean, and Southern European populations (collectively, SWAMSE) together with Northern European populations form one of five "continental" groups of global populations in many analyses of population relationships. This region is of great anthropologic and forensic interest but relationships of large numbers of populations within the region have not been able to be cleanly resolved with autosomal genetic markers. To examine the genetic boundaries to the SWAMSE region and whether internal structure can be detected we have assembled data for a total of 151 separate autosomal genetic markers on populations in this region and other parts of the world for a global set of 95 populations. The markers include 83 ancestry informative SNPs as singletons and 68 microhaplotype loci defined by 204 SNPs. The 151 loci are ancestry informative on a global scale, identifying at least five biogeographic clusters. One of those clusters is a clear grouping of 37 populations containing the SWAMSE plus northern European populations to the exclusion of populations in South Central Asia and populations from farther East. A refined analysis of the 37 populations shows the northern European populations clustering separately from the SWAMSE populations. Within Southwest Asia the Samaritans and Shabaks are distinct outliers. The Yemenite Jews, Saudi, Kuwaiti, Palestinian Arabs, and Southern Tunisians cluster together loosely while the remaining populations from Northern Iraq, Mediterranean Europe, the Caucasus region, and Iran cluster in a more complex graded fashion. The majority of the SWAMSE populations from the mainland of Southwest Asia form a cluster with little internal structure reflecting a very complex history of endogamy and migrations. The set of 151 DNA polymorphisms not only distinguishes major geographical regions globally but can distinguish ancestry to a small degree within geographical regions such as SWAMSE. We discuss forensic characteristics of the polymorphisms and also identify those that rank highest by Rosenberg's I_n measure for the SWAMSE region populations and for the global set of populations analyzed. DATA AVAILABILITY: Genotypes on all 151 markers for all 3790 individuals typed in the Kidd Lab on the 72 Kidd lab populations have been deposited in the Zenodo archive and can be freely accessed at https://doi.org/10.5281/zenodo.4658892. Some of the data has been made public previously as supplemental files appended to publications. Data for the additional individuals included in the analyses was taken from already public datasets as indicated in the text.

Evaluation of the Ion AmpliSeq™ PhenoTrivium Panel: MPS-Based Assay for Ancestry and Phenotype Predictions Challenged by Casework Samples

As the field of forensic DNA analysis has started to transition from genetics to genomics, new methods to aid in crime scene investigations have arisen. The development of informative single nucleotide polymorphism (SNP) markers has led the forensic community to question if DNA can be a reliable "eye-witness" and whether the data it provides can shed light on unknown perpetrators. We have developed an assay called the Ion AmpliSeq™ PhenoTrivium Panel, which combines three groups of markers: 41 phenotype- and 163 ancestry-informative autosomal SNPs together with 120 lineage-specific Y-SNPs. Here, we report the results of testing the assay's sensitivity and the predictions obtained for known reference samples. Moreover, we present the outcome of a blind study performed on real casework samples in order to understand the value and reliability of the information that would be provided to police investigators. Furthermore, we evaluated the accuracy of admixture prediction in Converge™ Software. The results show the panel to be a robust and sensitive assay which can be used to analyze casework samples. We conclude that the combination of the obtained predictions of phenotype, biogeographical ancestry, and male lineage can serve as a potential lead in challenging police investigations such as cold cases or cases with no suspect.

Evaluation of the Precision of Ancestry Inferences in South American Admixed Populations

Ancestry informative markers (AIMs) are used in forensic genetics to infer biogeographical ancestry (BGA) of individuals and may also have a prominent role in future police and identification investigations. In the last few years, many studies have been published reporting new AIM sets. These sets include markers (usually around 100 or less) selected with different purposes and different population resolutions. Regardless of the ability of these sets to separate populations from different continents or regions, the uncertainty associated with the estimates provided by these panels and their capacity to accurately report the different ancestral contributions in individuals of admixed populations has rarely been investigated. This issue is addressed in this study by evaluating different AIM sets. Ancestry inference was carried out in admixed South American populations, both at population and individual levels. The results of ancestry inferences using AIM sets with different numbers of markers among admixed reference populations were compared. To evaluate the performance of the different ancestry panels at the individual level, expected and observed estimates among families and their offspring were compared, considering that (1) the apportionment of ancestry in the offspring should be closer to the average ancestry of the parents, and (2) full siblings should present similar ancestry values. The results obtained illustrate the importance of having a good balance/compromise between not only the number of markers and their ability to differentiate ancestral populations, but also a balanced differentiation among reference groups, to obtain more precise values of genetic ancestry. This work also highlights the importance of estimating errors associated with the use of a limited number of markers. We demonstrate that although these errors have a moderate effect at the population level, they may have an important impact at the individual level. Considering that many AIM-sets are being described for inferences at the individual level and not at the population level, e.g., in association studies or the determination of a suspect's BGA, the results of this work point to the need of a more careful evaluation of the uncertainty associated with the ancestry estimates in admixed populations, when small AIM-sets are used.

Insights into the genetic characteristics and population structures of Chinese two Tibetan groups using 35 insertion/deletion polymorphic loci

Studying the genetic structure of each ethnic group is helpful to clarify the genetic background and trace back to the ethnic origin. Tibetan people have lived in the Qinghai-Tibet Plateau (mean elevation over 4500 m) for generations, and have well adapted to the high-altitude environment. Due to the relatively closed geographical environment, Tibetans have preserved their representative physical characteristics and genetic information, thereby become an important research group in human genetics. In this study, genetic characteristics and population structures of two Tibetan groups (Qinghai Tibetans and Tibet Tibetans) were revealed by 35 insertion/deletion polymorphism (DIP) loci, aiming to provide valuable genetic information for population genetic differentiation analyses and forensic identifications. The combined discrimination power, cumulative exclusion probability and combined match probability of the 35 DIP loci in Qinghai Tibetan and Tibet Tibetan groups were 0.9999999999999945, 0.9988, 5.56623 × 10^-15; and 0.9999999999999904, 0.9990, 9.69071 × 10^-15, respectively, indicating that the panel possessed a strong capability for Tibetan personal identifications. Population differentiations and genetic relationship analyses among the two studied Tibetan groups and other 27 comparison populations were carried out using the Nei's D_A genetic distances, population pairwise genetic distances F-statistics (F_ST), analysis of molecular variance (AMOVA), phylogenetic tree reconstruction, principal component analysis and STRUCTURE methods. Results demonstrated that the most intimate genetic relationships existed in these two Tibetan groups; and genetic similarities between two Tibetan groups and the populations from East Asia were much stronger than that between the Tibetan groups and other geographical populations. Furthermore, forensic ancestral informativeness assessments suggested that several loci could be regarded as ancestry informative markers inferring individual biogeographic origins as well as contributing to forensic anthropology and population genetic researches.

Development and validation of the EUROFORGEN NAME (North African and Middle Eastern) ancestry panel

Inference of biogeographic origin is an important factor in clinical, population and forensic genetics. The information provided by AIMs (Ancestry Informative Markers) can allow the differentiation of major continental population groups, and several AIM panels have been developed for this purpose. However, from these major population groups, Eurasia covers a wide area between two continents that is difficult to differentiate genetically. These populations display a gradual genetic cline from West Europe to South Asia in terms of allele frequency distribution. Although differences have been reported between Europe and South Asia, Middle East populations continue to be a target of further investigations due to the lack of genetic variability, therefore hampering their genetic differentiation from neighboring populations. In the present study, a custom-built ancestry panel was developed to analyze North African and Middle Eastern populations, designated the 'NAME' panel. The NAME panel contains 111 SNPs that have patterns of allele frequency differentiation that can distinguish individuals originating in North Africa and the Middle East when combined with a previous set of 126 Global AIM-SNPs.

Microhaplotypes in forensic genetics

Microhaplotype loci (microhaps, MHs) are a novel type of molecular marker of less than 300 nucleotides, defined by two or more closely linked SNPs associated in multiple allelic combinations. The value of these markers is enhanced by massively parallel sequencing (MPS), which allows the sequencing of both parental haplotypes at each of the many multiplexed loci. This review describes the features of these multi-SNP markers and documents their value in forensic genetics, focusing on individualization, biogeographic ancestry inference, and mixture deconvolution. Foreseeable applications also include missing person identification, relationship testing, and medical diagnostic applications. The technique is not restricted to humans.

Ancestry inference of 96 population samples using microhaplotypes

Microhaplotypes have become a new type of forensic marker with a great ability to identify and deconvolute mixtures because massively parallel sequencing (MPS) allows the alleles (haplotypes) of the multi-SNP loci to be determined directly for an individual. As originally defined, a microhaplotype locus is a short segment of DNA with two or more SNPs defining three or more haplotypes. The length is short enough, less than about 300 bp, that the read length of current MPS technology can produce a phase-known sequence of each chromosome of an individual. As part of the discovery phase of our studies, data on 130 microhaplotype loci with estimates of haplotype frequency data on 83 populations have been published. To provide a better picture of global allele frequency variation, we have now tested 13 more populations for 65 of the microhaplotype loci from among those with higher levels of inter-population gene frequency variation, including 8 loci not previously published. These loci provide clear distinctions among 6 biogeographic regions and provide some information distinguishing up to 10 clusters of populations.