PopAffiliator: online calculator for individual affiliation to a major population group based on 17 autosomal short tandem repeat genotype profile

Open-Access Worldwide Population STR Database Constructed Using High-Coverage Massively Parallel Sequencing Data Obtained from the 1000 Genomes Project

Achieving accurate STR genotyping by using next-generation sequencing data has been challenging. To provide the forensic genetics community with a reliable open-access STR database, we conducted a comprehensive genotyping analysis of a set of STRs of broad forensic interest obtained from 1000 Genome populations. We analyzed 22 STR markers using files of the high-coverage dataset of Phase 3 of the 1000 Genomes Project. We used HipSTR to call genotypes from 2504 samples obtained from 26 populations. We were not able to detect the D21S11 marker. The Hardy-Weinberg equilibrium analysis coupled with a comprehensive analysis of allele frequencies revealed that HipSTR was not able to identify longer alleles, which resulted in heterozygote deficiency. Nevertheless, AMOVA, a clustering analysis that uses STRUCTURE, and a Principal Coordinates Analysis showed a clear-cut separation between the four major ancestries sampled by the 1000 Genomes Consortium. Except for larger Penta D and Penta E alleles, and two very small Penta D alleles (2.2 and 3.2) usually observed in African populations, our analyses revealed that allele frequencies and genotypes offered as an open-access database are consistent and reliable.

Forensic genetics through the lens of Lewontin: population structure, ancestry and race

In his famous 1972 paper, Richard Lewontin used ‘classical’ protein-based markers to show that greater than 85% of human genetic diversity was contained within, rather than between, populations. At that time, these same markers also formed the basis of forensic technology aiming to identify individuals. This review describes the evolution of forensic genetic methods into DNA profiling, and how the field has accounted for the apportionment of genetic diversity in considering the weight of forensic evidence. When investigative databases fail to provide a match to a crime-scene profile, specific markers can be used to seek intelligence about a suspect: these include inferences on population of origin (biogeographic ancestry) and externally visible characteristics, chiefly pigmentation of skin, hair and eyes. In this endeavour, ancestry and phenotypic variation are closely entangled. The markers used show patterns of inter- and intrapopulation diversity that are very atypical compared to the genome as a whole, and reinforce an apparent link between ancestry and racial divergence that is not systematically present otherwise. Despite the legacy of Lewontin's result, therefore, in a major area in which genetics coincides with issues of public interest, methods tend to exaggerate human differences and could thereby contribute to the reification of biological race.

This article is part of the theme issue ‘Celebrating 50 years since Lewontin's apportionment of human diversity’.

Could routine forensic STR genotyping data leak personal phenotypic information?

The application of forensic genetic markers must comply with privacy rights and legal policies on a premise that the markers do not expose phenotypic information. The most widely-used short tandem repeats (STRs) are generally viewed as 'junk' DNA because most STRs are located in non-coding regions and therefore refrain from leaking phenotypic traits. But with a deepening understanding of phenotypes and underlying genetic structure, whether STRs could potentially reflect any phenotypic information may need re-examining. Therefore, we performed the following analyses. First, we analyzed the association between 15 STRs and three facial characteristics (single or double eyelid, with or without epicanthus, unattached or attached earlobe) on 721 unrelated Han Chinese individuals. Then, we collected 27199 individuals' STRs and geographic data from the literature to investigate the association between STRs and bio-geographic information, and predict geographic information by STRs on additional 1993 unrelated individuals. We found that there was scarcely any association between STRs with studied facial characteristics. Although allele19 in D2S1338 and allele 18 in FGA (P = 0.0032, P = 0.0030, respectively after Bonferroni correction) showed statistical significance, the prediction effectiveness was very low. For the STRs and bio-geographic information, the principal component analysis showed the first three components could explain 87.7% of the variance, but the prediction accuracy only reached 25.2%. We demonstrated that the forensic phenotypes are usually complex traits, it is hardly possible to uncover phenotypic information by testing only dozens of STR loci.

Comprehensive Insights Into Forensic Features and Genetic Background of Chinese Northwest Hui Group Using Six Distinct Categories of 231 Molecular Markers

The Hui minority is predominantly composed of Chinese-speaking Islamic adherents distributed throughout China, of which the individuals are mainly concentrated in Northwest China. In the present study, we employed the length and sequence polymorphisms-based typing system of 231 molecular markers, i.e., amelogenin, 22 phenotypic-informative single nucleotide polymorphisms (PISNPs), 94 identity-informative single nucleotide polymorphisms (IISNPs), 24 Y-chromosomal short tandem repeats (Y-STRs), 56 ancestry-informative single nucleotide polymorphisms (AISNPs), 7 X-chromosomal short tandem repeats (X-STRs), and 27 autosomal short tandem repeats (A-STRs), into 90 unrelated male individuals from the Chinese Northwest Hui group to comprehensively explore its forensic characteristics and genetic background. Total of 451 length-based and 652 sequence-based distinct alleles were identified from 58 short tandem repeats (STRs) in 90 unrelated Northwest Hui individuals, denoting that the sequence-based genetic markers could pronouncedly provide more genetic information than length-based markers. The forensic characteristics and efficiencies of STRs and IISNPs were estimated, both of which externalized high polymorphisms in the Northwest Hui group and could be further utilized in forensic investigations. No significant departure from the Hardy-Weinberg equilibrium (HWE) expectation was observed after the Bonferroni correction. Additionally, four group sets of reference population data were exploited to dissect the genetic background of the Northwest Hui group separately from different perspectives, which contained 26 populations for 93 IISNPs, 58 populations for 17 Y-STRs, 26 populations for 55 AISNPs (raw data), and 109 populations for 55 AISNPs (allele frequencies). As a result, the analyses based on the Y-STRs indicated that the Northwest Hui group primarily exhibited intimate genetic relationships with reference Hui groups from Chinese different regions except for the Sichuan Hui group and secondarily displayed close genetic relationships with populations from Central and West Asia, as well as several Chinese groups. However, the AISNP analyses demonstrated that the Northwest Hui group shared more intimate relationships with current East Asian populations apart from reference Hui group, harboring the large proportion of ancestral component contributed by East Asia.

Comparative analysis of allele variation using allele frequencies according to sample size in Korean population

Background

Allele frequency using short tandem repeats (STRs) is used to calculate likelihood ratio for database match, to interpret DNA mixture and to estimate ethnic groups in forensic genetics. In Korea, three population studies for 23 STR loci have been conducted with different sample size for forensic purposes.

Objective

We performed comparative analysis to determine how the difference of sample size affects the allele frequency and allele variation within same ethnic population (i.e. Korean). Furthermore, this study was conducted to check how the sampling group and multiplex kit also affect allele variation such as rare alleles and population specific alleles.

Methods

To compare allele variation, we used allele frequencies of three population data published from three Korean forensic research groups. Allele frequencies were calculated using different sample sizes and multiplex kits: 526, 1000, and 2000 individuals, respectively.

Results

The results showed the different distribution of allele frequencies in some loci. There was also a difference in the number of rare alleles observed by the sample size and sampling bias. In particular, an allele of 9.1 in the D2S441 locus was not observed in population study with 526 individuals due to multiplex kits.

Conclusion

Because the allele frequencies play an important role in forensic genetics, even if the samples are derived from the same population, it is important to consider the effects of sample size, sampling bias, and selection of multiplex kits in population studies.

The analysis of ancestry with small-scale forensic panels of genetic markers

In the last 10 years, forensic genetic analysis has been extended beyond identification tests that link a suspect to crime scene evidence using standard DNA profiling, to new supplementary tests that can provide information to investigators about a suspect in the absence of a database hit or eyewitness testimony. These tests now encompass the prediction of physical appearance, ancestry and age. In this review, we give a comprehensive overview of the full range of DNA-based ancestry inference tests designed to work with forensic contact traces, when the level of DNA is often very low or highly degraded. We outline recent developments in the design of ancestry-informative marker sets, forensic assays that use capillary electrophoresis or massively parallel sequencing, and the statistical analysis frameworks that examine the test profile and compares it to reference population variation. Three casework ancestry analysis examples are described which were successfully accomplished in the authors' laboratory, where the ancestry information obtained was critical to the outcome of the DNA analyses made.

Next generation sequencing of a set of ancestry-informative SNPs: ancestry assignment of three continental populations and estimating ancestry composition for Mongolians

When traditional short tandem repeat profiling fails to provide valuable information to arrest the criminal, forensic ancestry inference of the biological samples left at the crime scene will probably offer investigative leads and facilitate the investigation process of the case. That is why there are consistent efforts in developing panels for ancestry inference in forensic science. Presently, a 30-plex next generation sequencing-based assay was exploited in this study by assembling well-differentiated single nucleotide polymorphisms for ancestry assignment of unknown individuals from three continental populations (African, European and East Asian). And meanwhile, relatively balanced population-specific differentiation values were maintained to avoid the over-estimation or under-estimation of co-ancestry proportions in individuals with admixed ancestry. The principal component analysis and STRUCTURE analysis of reference populations, test populations and the studied Mongolian group indicated that the novel assay was efficient enough to determine the ancestry origin of an unknown individual from the three continental populations. Besides, ancestry membership proportion estimations for the Mongolian group revealed that a large fraction of the ancestry was contributed by East Asian genetic component (approximately 83.9%), followed by European (approximately 12.6%) and African genetic components (approximately 3.5%), respectively. And next generation sequencing technology applied in this study offers possibility to incorporate more single nucleotide polymorphisms for individual identification and phenotype prediction into the same assay to provide as many as possible investigative clues in the future.

A multivariate statistical approach for the estimation of the ethnic origin of unknown genetic profiles in forensic genetics

DNA typing and genetic profile data interpretation are among the most relevant topics in forensic science; among other applications, genetic profile's capability to distinguish biogeographic information about population groups, subgroups and affiliations have been largely explored in the last decade. In fact, for investigative and intelligence purposes, it is extremely useful to identify subjects and estimate their biogeographic origins by examining the recovered DNA profiles from evidence on a crime scene. Current approaches for BiogeoGraphic Ancestry (BGA) estimation using STRs profiles are usually based on Bayesian methods, which quantify the evidence in terms of likelihood ratio, supporting or not the hypothesis that a certain profile belongs to a specific ethnic group. The present study provides an alternative approach to the likelihood ratio method that involves multivariate data analysis strategies for the estimation of multiple populations. Starting from the well-known NIST US autosomal STRs dataset involving African-American, Asian, and Caucasian individuals, and moving towards further and more geographically restricted populations (such as Northern Africans vs sub-Saharan Africans, Afghans vs Iraqis and Italians vs Romanians), powerful multivariate techniques such as Sparse and Logistic Principal Component Analysis (SL-PCA), Sparse Partial Least Squares-Discriminant Analysis (sPLS-DA) and Support Vector Machines (SVM) were employed and their discriminating power was also compared. Both sPLS-DA and SVM techniques provided robust classifications, yielding high sensitivity and specificity models capable of discriminating populations on ethnic basis. This application may represent a powerful and dynamic tool for law enforcement agencies whenever a standard autosomal STR profile is obtained from the biological evidence collected at a crime scene or recovered during mass-disaster and missing person investigations.

Population genetics data for 22 autosomal STR loci in European, South Asian and African populations using SureID® 23comp Human DNA Identification Kit

Allele frequency data for 22 short tandem repeat loci; D18S1364, D1S1656, D13S325, D5S2800, D9S1122, D4S2366, D3S1744, D12S391, D11S2368, D21S2055, D20S482, D8S1132, D7S3048, D2S441, D19S253, D10S1248, D17S1301, D22-GATA198B05, D16S539, D6S474, D14S1434 and D15S659 from the SureID^® 23comp Human DNA Identification Kit have been determined for unrelated individuals in European, South Asian and African populations. Deviations from Hardy-Weinberg equilibrium were observed in loci D1S1656 and D19S253 in European; D18S1364, D6S474 and D14S1434 in South Asian; and D9S1122 and D8S1132 in African populations (p-value <0.05). However, after Bonferroni correction no significant deviations were observed (p-value <0.002). The most discriminating loci were D1S1656 and D12S391 for European (PD=0.977), D21S2055 for South Asian (PD=0.980), and D21S2055 and D7S3048 for African (PD=0.972) populations. The match probabilities were 1 in 6.7×10²⁵ for European, 1 in 1.4×10²⁶ for South Asian and 1 in 1.6×10²⁶ for African populations. These findings established the high discriminatory capacity and robustness of the tested STR loci for forensic identification and kinship testing.

A genetic portrait of the South-Eastern Carpathians based on autosomal short tandem repeats loci used in forensics

OBJECTIVES

This work aimed to describe the genetic landscape of the Balkan Peninsula, as revealed by STR markers commonly used in forensics and spatial methods specifically developed for genetic data.

METHODS

We generated and analyzed 16 short tandem repeats (STRs) autosomal genotypes in 287 subjects from ten administrative/geographical regions of Eastern Europe (Romania and the Republic of Moldova). We report estimates of the allele frequencies in these sub-populations, their fixation indexes, and use these results to complement previous spatial analyses of Southern Europe.

RESULTS

In seven out of ten analyzed regional samples the heterozygosity, averaged across loci, was lower than expected. The average Fis was 0.011. Among the 16 loci, five returned a significant fixation index Fst. The composite Fst across the 16 loci, among the 10 regional samples, was 0.00417, a figure twice as large as that obtained with the same markers across the entire Northern Mediterranean. The first spatial principal component (sPC1) returned the picture of a Central-European pattern of frequencies for the Carpathians, which extended to the Southern boundary of the Balkan Peninsula. However, the 8 alleles extracted by sPC1 returned a picture of a strong reduction of the migration rate in the Carpathian region, mostly between the inner locations.

CONCLUSIONS

Our results revealed an unexpected heterogeneity in the area. We believe that populations from some regions will require treatment as distinct entities when considered in forensic applications.

Enlarging the gene-geography of Europe and the Mediterranean area to STR loci of common forensic use: longitudinal and latitudinal frequency gradients

BACKGROUND

Tetranucleotide Short Tandem Repeats (STRs) for human identification and common use in forensic cases have recently been used to address the population genetics of the North-Eastern Mediterranean area. However, to gain confidence in the inferences made using STRs, this kind of analysis should be challenged with changes in three main aspects of the data, i.e. the sizes of the samples, their distance across space and the genetic background from which they are drawn.

AIM

To test the resilience of the gradients previously detected in the North-Eastern Mediterranean to the enlargement of the surveyed area and population set, using revised data.

SUBJECTS AND METHODS

STR genotype profiles were obtained from a publicly available database (PopAffilietor databank) and a dataset was assembled including >7000 subjects from the Arabian Peninsula to Scandinavia, genotyped at eight loci. Spatial principal component analysis (sPCA) was applied and the frequency maps of the nine alleles which contributed most strongly to sPC1 were examined in detail.

RESULTS

By far the greatest part of diversity was summarised by a single spatial principal component (sPC1), oriented along a SouthEast-to-NorthWest axis. The alleles with the top 5% squared loadings were TH01(9.3), D19S433(14), TH01(6), D19S433(15.2), FGA(20), FGA(24), D3S1358(14), FGA(21) and D2S1338(19). These results confirm a clinal pattern over the whole range for at least four loci (TH01, D19S433, FGA, D3S1358).

CONCLUSIONS

Four of the eight STR loci (or even alleles) considered here can reproducibly capture continental arrangements of diversity. This would, in principle, allow for the exploitation of forensic data to clarify important aspects in the formation of local gene pools.

Age Estimation with DNA: From Forensic DNA Fingerprinting to Forensic (Epi)Genomics: A Mini-Review

Forensic genetics developed from protein-based techniques a quarter of a century ago and became famous as “DNA fingerprinting,” this being based on restriction fragment length polymorphisms (RFLPs) of high-molecular-weight DNA. The amplification of much smaller short tandem repeat (STR) sequences using the polymerase chain reaction soon replaced RFLP analysis and advanced to become the gold standard in genetic identification. Meanwhile, STR multiplexes have been developed and made commercially available which simultaneously amplify up to 30 STR loci from as little as 15 cells or fewer. The enormous information content that comes with the large variety of observed STR genotypes allows for genetic individualisation (with the exception of identical twins). Carefully selected core STR loci form the basis of intelligence-led DNA databases that provide investigative leads by linking unsolved crime scenes and criminals through their matched STR profiles. Nevertheless, the success of modern DNA fingerprinting depends on the availability of reference material from suspects. In order to provide new investigative leads in cases where such reference samples are absent, forensic scientists started to explore the prediction of phenotypic traits from the DNA of the evidentiary sample. This paradigm change now uses DNA and epigenetic markers to forecast characteristics that are useful to triage further investigative work. So far, the best investigated externally visible characteristics are eye, hair and skin colour, as well as geographic ancestry and age. Information on the chronological age of a stain donor (or any sample donor) is elemental for forensic investigations in a number of aspects and has, therefore, been explored by researchers in some detail. Among different methodological approaches tested to date, the methylation-sensitive analysis of carefully selected DNA markers (CpG sites) has brought the most promising results by providing prediction accuracies of ±3–4 years, which can be comparable to, or even surpass those from, eyewitness reports. This mini-review puts recent developments in age estimation via (epi)genetic methods in the context of the requirements and goals of forensic genetics and highlights paths to follow in the future of forensic genomics.

Genetic variation and population structure of Botswana populations as identified with AmpFLSTR Identifiler short tandem repeat (STR) loci

Population structure was investigated in 990 Botswana individuals according to ethno-linguistics, Bantu and Khoisan, and geography (the nine administrative districts) using the Identifiler autosomal microsatellite markers. Genetic diversity and forensic parameters were calculated for the overall population, and according to ethno-linguistics and geography. The overall combined power of exclusion (CPE) was 0.9999965412 and the combined match probability 6,28 × 10^-19. CPE was highest for the Khoisan Tuu ethnolinguistic group and the Northeast District at 0.9999582029 and 0.9999922652 respectively. CMP ranged from 6.28 × 10^-19 (Khoisan Tuu) to 1,02 × 10^-18 (Northwest district). Using pairwise genetic distances (F_ST), analysis of molecular variance (AMOVA), factorial correspondence analysis (FCA), and the unsupervised Bayesian clustering method found in STRUCTURE and TESS, ethno-linguistics were found to have a greater influence on population structure than geography. FCA showed clustering between Bantu and Khoisan, and within the Bantu. This Bantu sub-structuring was not seen with STRUCTURE and TESS, which detected clustering only between Bantu and Khoisan. The patterns of population structure revealed highlight the need for regional reference databases that include ethno-linguistic and geographic location information. These markers have important potential for bio-anthropological studies as well as for forensic applications.

Genetic structure of Tibetan populations in Gansu revealed by forensic STR loci

The origin and diversification of Sino-Tibetan speaking populations have been long-standing hot debates. However, the limited genetic information of Tibetan populations keeps this topic far from clear. In the present study, we genotyped 15 forensic autosomal short tandem repeats (STRs) from 803 unrelated Tibetan individuals from Gansu Province (635 from Gannan and 168 from Tianzhu) in northwest China. We combined these data with published dataset to infer a detailed population affinities and genetic substructure of Sino-Tibetan populations. Our results revealed Tibetan populations in Gannan and Tianzhu are genetically very similar with Tibetans from other regions. The Tibetans in Tianzhu have received more genetic influence from surrounding lowland populations. The genetic structure of Sino-Tibetan populations was strongly correlated with linguistic affiliations. Although the among-population variances are relatively small, the genetic components for Tibetan, Lolo-Burmese, and Han Chinese were quite distinctive, especially for the Deng, Nu, and Derung of Lolo-Burmese. Han Chinese but not Tibetans are suggested to share substantial genetic component with southern natives, such as Tai-Kadai and Hmong-Mien speaking populations, and with other lowland East Asian populations, which implies there might be extensive gene flow between those lowland groups and Han Chinese after Han Chinese were separated from Tibetans. The dataset generated in present study is also valuable for forensic identification and paternity tests in China.

Genetic evidence for an East Asian origin of Chinese Muslim populations Dongxiang and Hui

There is a long-going debate on the genetic origin of Chinese Muslim populations, such as Uygur, Dongxiang, and Hui. However, genetic information for those Muslim populations except Uygur is extremely limited. In this study, we investigated the genetic structure and ancestry of Chinese Muslims by analyzing 15 autosomal short tandem repeats in 652 individuals from Dongxiang, Hui, and Han Chinese populations in Gansu province. Both genetic distance and Bayesian-clustering methods showed significant genetic homogeneity between the two Muslim populations and East Asian populations, suggesting a common genetic ancestry. Our analysis found no evidence of substantial gene flow from Middle East or Europe into Dongxiang and Hui people during their Islamization. The dataset generated in present study are also valuable for forensic identification and paternity tests in China.

Inference of Ancestry in Forensic Analysis I: Autosomal Ancestry-Informative Marker Sets

An expanding choice of ancestry-informative marker single nucleotide polymorphisms (AIM-SNPs) is becoming available for the forensic user in the form of sensitive SNaPshot-based tests or in alternative single-base extension genotyping systems (e.g., Sequenom iPLEX) that can be adapted for analysis with SNaPshot. In addition, alternative ancestry-informative variation: Indels and STRs can be analyzed using direct PCR-to-CE techniques that offer the possibility to detect mixed profiles. We review the current forensically viable AIM panels, their optimized PCR multiplexes, and the population differentiation power they offer. We also describe how improved population divergence balance can be achieved with the enlarged multiplex scales of next-generation sequencing approaches to enable analysis of admixed individuals without biased estimation of co-ancestry proportions.

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.

Introducing the Forensic Research/Reference on Genetics knowledge base, FROG-kb

BACKGROUND

Online tools and databases based on multi-allelic short tandem repeat polymorphisms (STRPs) are actively used in forensic teaching, research, and investigations. The Fst value of each CODIS marker tends to be low across the populations of the world and most populations typically have all the common STRP alleles present diminishing the ability of these systems to discriminate ethnicity. Recently, considerable research is being conducted on single nucleotide polymorphisms (SNPs) to be considered for human identification and description. However, online tools and databases that can be used for forensic research and investigation are limited.

METHODS

The back end DBMS (Database Management System) for FROG-kb is Oracle version 10. The front end is implemented with specific code using technologies such as Java, Java Servlet, JSP, JQuery, and GoogleCharts.

RESULTS

We present an open access web application, FROG-kb (Forensic Research/Reference on Genetics-knowledge base, http://frog.med.yale.edu), that is useful for teaching and research relevant to forensics and can serve as a tool facilitating forensic practice. The underlying data for FROG-kb are provided by the already extensively used and referenced ALlele FREquency Database, ALFRED (http://alfred.med.yale.edu). In addition to displaying data in an organized manner, computational tools that use the underlying allele frequencies with user-provided data are implemented in FROG-kb. These tools are organized by the different published SNP/marker panels available. This web tool currently has implemented general functions possible for two types of SNP panels, individual identification and ancestry inference, and a prediction function specific to a phenotype informative panel for eye color.

CONCLUSION

The current online version of FROG-kb already provides new and useful functionality. We expect FROG-kb to grow and expand in capabilities and welcome input from the forensic community in identifying datasets and functionalities that will be most helpful and useful. Thus, the structure and functionality of FROG-kb will be revised in an ongoing process of improvement. This paper describes the state as of early June 2012.