Development and evaluations of the ancestry informative markers of the VISAGE Enhanced Tool for Appearance and Ancestry

The VISAGE Enhanced Tool for Appearance and Ancestry (ET) has been designed to combine markers for the prediction of bio-geographical ancestry plus a range of externally visible characteristics into a single massively parallel sequencing (MPS) assay. We describe the development of the ancestry panel markers used in ET, and the enhanced analyses they provide compared to previous MPS-based forensic ancestry assays. As well as established autosomal single nucleotide polymorphisms (SNPs) that differentiate sub-Saharan African, European, East Asian, South Asian, Native American, and Oceanian populations, ET includes autosomal SNPs able to efficiently differentiate populations from Middle East regions. The ability of the ET autosomal ancestry SNPs to distinguish Middle East populations from other continentally defined population groups is such that characteristic patterns for this region can be discerned in genetic cluster analysis using STRUCTURE. Joint cluster membership estimates showing individual co-ancestry that signals North African or East African origins were detected, or cluster patterns were seen that indicate origins from central and Eastern regions of the Middle East. In addition to an augmented panel of autosomal SNPs, ET includes panels of 85 Y-SNPs, 16 X-SNPs and 21 autosomal Microhaplotypes. The Y- and X-SNPs provide a distinct method for obtaining extra detail about co-ancestry patterns identified in males with admixed backgrounds. This study used the 1000 Genomes admixed African and admixed American sample sets to fully explore these enhancements to the analysis of individual co-ancestry. Samples from urban and rural Brazil with contrasting distributions of African, European, and Native American co-ancestry were also studied to gauge the efficiency of combining Y- and X-SNP data for this purpose. The small panel of Microhaplotypes incorporated in ET were selected because they showed the highest levels of haplotype diversity amongst the seven population groups we sought to differentiate. Microhaplotype data was not formally combined with single-site SNP genotypes to analyse ancestry. However, the haplotype sequence reads obtained with ET from these loci creates an effective system for de-convoluting two-contributor mixed DNA. We made simple mixture experiments to demonstrate that when the contributors have different ancestries and the mixture ratios are imbalanced (i.e., not 1:1 mixtures) the ET Microhaplotype panel is an informative system to infer ancestry when this differs between the contributors.     

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.     

Revision of the SNPforID 34-plex forensic ancestry test: Assay enhancements,standard reference sample genotypes and extended population studies

A revision of an established 34 SNP forensic ancestry test has been made by swapping the under-performing rs727811 component SNP with the highly informative rs3827760 that shows a near-fixed East Asian specific allele. We collated SNP variability data for the revised SNP set in 66 reference populations from 1000 Genomes and HGDP-CEPH panels and used this as reference data to analyse four U.S. populations showing a range of admixture patterns. The U.S. Hispanics sample in particular displayed heterogeneous values of co-ancestry between European, Native American and African contributors, likely to reflect in part, the way this disparate group is defined using cultural as well as population genetic parameters. The genotyping of over 700 U.S. population samples also provided the opportunity to thoroughly gauge peak mobility variation and peak height ratios observed from routine use of the single base extension chemistry of the 34-plex test. Finally, the genotyping of the widely used DNA profiling Standard Reference Material samples plus other control DNAs completes the audit of the 34-plex assay to allow forensic practitioners to apply this test more readily in their own laboratories.     

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.     

Comparison of the genetic background of different Colombian populations using the SNPforID 52plex identification panel

Various strategies for analysing SNP markers and genotyping have been published with the goal of obtaining informative profiles from biological samples that contain only small amounts of template and/or degraded DNA. In this study, a multiplex assay of 52 autosomal single-nucleotide polymorphisms (SNPs) was used to analyse 438 individuals from urban populations from different regions of Colombia, as well as a sample of 50 Native American individuals of the Pastos ethnic group from Nariño. To determine if significant differences in these 52 SNPs exist between the distinct regions of Colombia, genetic distance and admixture analyses were performed based on the available data for 17 different Colombian population groups and for population groups from Africa, Europe and America. The results demonstrate significant differences between the populations from the Southwest Andean, Central-West Andean, Central-East Andean, Orinoquian and northern Colombian Pacific Coast regions. Most of the regions exhibited a European and Native American admixture. One exception is the population from the region of Chocó (on the northern Pacific Coast), which exhibits a high proportion of African admixture (54 %). From the observed genetic backgrounds, it is possible to conclude that a single reference database for the entire country would not be suitable for forensic purposes. The allele frequencies and the forensically relevant parameters were calculated for all of the markers in each Colombian region with significant values for the combined matching probability (power of discrimination ≥0.99999999999999990) and the combined probability of exclusion (≥0.9990) in trios that were obtained from all of the population groups.     

Y chromosome comparative analysis of Rondônia with other Brazilian populations

In the present study, a Brazilian population, located in the Rond?nia state, was genetically characterized for a set of Y chromosome specific STRs included in the Applied Biosystems kit (AmpF?STR?Yfiler?), which allows the simultaneous amplification of 16 markers: DYS19, DYS389I, DYS389II, DYS390, DYS391, DYS392, DYS393, DYS385a/b, DYS438, DYS439, DYS437, DYS448, DYS456, DYS458, DYS635 and GATA H4. The studied population from Rond?nia state, in the North of Brazil, included individuals with admixed Native American, African and European ancestry. When comparing Rond?nia with other Brazilian populations no significant genetic distances were found. In the comparison with other worldwide populations, although a predominant male European influence could be detected, there were significant differences with some populations from Central and South America and Africa.     

Eurasiaplex: A forensic SNP assay for differentiating European and South Asian ancestries

We have selected a set of single nucleotide polymorphisms (SNPs) with the specific aim of differentiating European and South Asian ancestries. The SNPs were combined into a 23-plex SNaPshot primer extension assay: Eurasiaplex, designed to complement an existing 34-plex forensic ancestry test with both marker sets occupying well-spaced genomic positions, enabling their combination as single profile submissions to the Bayesian Snipper forensic ancestry inference system. We analyzed the ability of Eurasiaplex plus 34plex SNPs to assign ancestry to a total 1648 profiles from 16 European, 7 Middle East, 13 Central-South Asian and 21 East Asian populations. Ancestry assignment likelihoods were estimated from Snipper using training sets of five-group data (three Eurasian groups, East Asian and African genotypes) and four-group data (Middle East genotypes removed). Five-group differentiations gave assignment success of 91% for NW European populations, 72% for Middle East populations and 39% for Central-South Asian populations, indicating Middle East individuals are not reliably differentiated from either Europeans or Central-South Asians. Four-group differentiations provided markedly improved assignment success rates of 97% for most continental Europeans tested (excluding Turkish and Adygei at the far eastern edge of Europe) and 95% for Central-South Asians, despite applying a probability threshold for the highest likelihood ratio above ‘100 times more likely’. As part of the assessment of the sensitivity of Eurasiaplex to analyze challenging forensic material we detail Eurasiaplex and 34-plex SNP typing to infer ancestry of a cranium recovered from the sea, achieving 82% SNP genotype completeness. Therefore, Eurasiaplex provides an informative and forensically robust approach to the differentiation of European and South Asian ancestries amongst Eurasian populations.     

Exploring iris colour prediction and ancestry inference in admixed populations of South America

New DNA-based predictive tests for physical characteristics and inference of ancestry are highly informative tools that are being increasingly used in forensic genetic analysis. Two eye colour prediction models: a Bayesian classifier – Snipper and a multinomial logistic regression (MLR) system for the Irisplex assay, have been described for the analysis of unadmixed European populations. Since multiple SNPs in combination contribute in varying degrees to eye colour predictability in Europeans, it is likely that these predictive tests will perform in different ways amongst admixed populations that have European co-ancestry, compared to unadmixed Europeans. In this study we examined 99 individuals from two admixed South American populations comparing eye colour versus ancestry in order to reveal a direct correlation of light eye colour phenotypes with European co-ancestry in admixed individuals. Additionally, eye colour prediction following six prediction models, using varying numbers of SNPs and based on Snipper and MLR, were applied to the study populations. Furthermore, patterns of eye colour prediction have been inferred for a set of publicly available admixed and globally distributed populations from the HGDP-CEPH panel and 1000 Genomes databases with a special emphasis on admixed American populations similar to those of the study samples.     

Distribution of Y chromosomal STRs loci in Mayan and Mestizo populations from Guatemala

In this study, a sample of 225 Guatemalan males, comprising 115 Mestizo-Guatemalan and 110 Mayan-Guatemalan, was typed for 17 Y-short tandem repeats (STRs) loci (DYS19, DYS389I, DYS389II, DYS390, DYS391, DYS392, DYS393, DYS437, DYS438, DYS439, DYS448, DYS456, DYS458, DYS635, YGATA_H4.1 and DYS385a/b). The haplotype diversity (H=1) and discrimination capacity (96.86%) were calculated. Analysis of molecular variance (AMOVA) demonstrated a low but significant interpopulation differentiation when compared with the results obtained when we confront the Mestizo and Mayan populations with the European populations. Furthermore, the genetic variability and differences among the American, African, Asian, and European populations were analyzed with the software Statistica 9.1. In addition, the genetic distances were also calculated using other published data. Reynolds and Slatkin?s genetic distance was visualized using the multidimensional scaling (MDS) analysis. All the analysis performed locates the Mayan population next to the Native American population, while Guatemalan-Mestizo population was found to be between these populations and the European population, similar to other Mestizo one. The implementation of the estimation of individual ancestry proportions of the whole population sample showed the presence of two well-differentiated population groups.     

The Global AIMs Nano set: A 31-plex SNaPshot assay of ancestry-informative SNPs

A 31-plex SNaPshot assay, named ‘Global AIMs Nano’, has been developed by reassembling the most differentiated markers of the EUROFORGEN Global AIM-SNP set. The SNPs include three tri-allelic loci and were selected with the goal of maintaining a balanced differentiation of: Africans, Europeans, East Asians, Oceanians and Native Americans. The Global AIMs Nano SNP set provides higher divergence between each of the five continental population groups than previous small-scale AIM sets developed for forensic ancestry analysis with SNaPshot. Both of these characteristics minimise potential bias when estimating co-ancestry proportions in individuals with admixed ancestry; more likely to be observed when using markers disproportionately informative for only certain population group comparisons. The optimised multiplex is designed to be easily implemented using standard capillary electrophoresis regimes and has been used to successfully genotype challenging forensic samples from highly degraded material with low level DNA. The ancestry predictive performance of the Global AIMs Nano set has been evaluated by the analysis of samples previously characterised with larger AIM sets.     

Male lineage strata of Brazilian population disclosed by the simultaneous analysis of STRs and SNPs

Brazil has a large territory divided in five geographical regions harboring highly diverse populations that resulted from different degrees and modes of admixture between Native Americans, Europeans and Africans. In this study, a sample of 605 unrelated males was genotyped for 17 Y-STRs and 46 Y-SNPs aiming a deep characterization of the male gene pool of Rio de Janeiro and its comparison with other Brazilian populations. High values of Y-STR haplotype diversity (0.9999 ± 0.0001) and Y-SNP haplogroup diversity (0.7589 ± 0.0171) were observed. Population comparisons at both haplotype and haplogroup levels showed significant differences between Brazilian South Eastern and Northern populations that can be explained by differences in the proportion of African and Native American Y chromosomes. Statistical significant differences between admixed urban samples from the five regions of Brazil were not previously detected at haplotype level based on smaller size samples from South East, which emphasizes the importance of sample size to detected population stratification for an accurate interpretation of profile matches in kinship and forensic casework. Although not having an intra-population discrimination power as high as the Y-STRs, the Y-SNPs are more powerful to disclose differences in admixed populations. In this study, the combined analysis of these two types of markers proved to be a good strategy to predict population sub-structure, which should be taken into account when delineating forensic database strategies for Y chromosome haplotypes.     

A view of the maternal inheritance of Espírito Santo populations: The contrast between the admixed and Pomeranian descent groups

In addition to its valuable utility in forensic investigations, mitochondrial DNA (mtDNA) analysis is a reliable tool to uncover the origins of admixed populations, such as Brazilians. The state of Espírito Santo (ES), similar to other coastal Brazilian states, has a population shaped by 3 main ancestral roots: Amerindian, African and European. Among the latter, the Pomeranian descendants stand out due to the preservation of the traditional aspects of their culture, especially the Pomeranian language. Despite the genetic data already available, there is no mtDNA database that adequately reflects the diversity, the geographic distribution, and the origins of the maternal lineages from ES. To increase the knowledge of maternal ancestry and to investigate the population’s genetic stratification, a total of 291 samples were collected in the 4 macroregions (Metropolitan, South, Central and North) of ES and in the Pomeranian communities. Complete control region data were produced for the general (N=214) and Pomeranian (N=77) groups. Regarding the general population, the high values of haplotype diversity (H=99.9%) and pairwise differences (MNPD=16.9) found are in agreement with those reported for other populations in the southeast region of the country. Regarding maternal inheritance, the ES populations stood out due to the predominance of European haplogroups (49.5%), although the North macroregion had a higher African ancestry (47.1%). Among the Pomeranians, the lowest MNPD value (11.2) and the high percentage of shared haplotypes (15%) were indicative of founder events. The F_ST analysis showed that the Pomeranians (98.7% of European lineages) are genetically isolated from the other admixed populations in Brazil. This study demonstrated that the ES state contains singularities regarding the intrapopulational and interpopulational diversity of mtDNA. Even after 5 centuries of interethnic admixture, the present-day population of Espírito Santo harbors genetic marks that trace back to the historical aspects of its formation.     

Demand for larger Y-STR reference databases in ethnic melting-pot countries: Argentina as a test case

The Y chromosome behaves as a single locus. Its genetic information is useful in forensic casework, deficiency kinship testing, and population genetics studies. Continuous increases of loci number within commercial kits forced modification of worldwide reference databases. In Pan American countries, like Argentina, diverse parental ethnic groups contributed to the extant admixed urban populations. We report 509 additional haplotypes of 23 Y-STRs from donors inhabiting urban areas of six Argentinean provinces: Buenos Aires, Santiago del Estero, Santa Cruz, Rio Negro, Santa Fe, and Formosa. To better understand the demographic landscape of the admixed urban paternal lineages, structural analysis was performed using published data from other Argentinean provinces. AMOVA by Rst distance and inferred haplogroups by two predictive online software methods based on haplotypes yielded complementary results with respect to detected population structure, probably due to the different proportions of the Native American Q3-M3 haplogroup in the studied samples. This situation, which is common to most North, Meso, and South American countries, underscores the need for the additional step of typing specific SNPs for haplogroup diagnosis. We propose organizing Y-STR haplotype reference databases according to the most frequent haplogroups detected in a given admixed population.

     

Allelic frequencies and statistical data obtained from 48 AIM INDEL loci in an admixed population from the Brazilian Amazon

Allelic frequencies of 48 informative insert-delete (INDEL) loci were obtained from a sample set of 130 unrelated individuals living in Macapá, a city located in the northern Amazon region, in Brazil. The values of heterozygosity (H), polymorphic information content (PIC), power of discrimination (PD), power of exclusion (PE), matching probability (MP) and typical paternity index (TPI) were calculated and showed the forensic efficiency of these genetic markers. Based on the allele frequency obtained for the population of Macapá, we estimated an interethnic admixture for the three parental groups (European, Native American and African) of, respectively, 50%, 21% and 29%. Comparing these allele frequencies with those of other Brazilian populations and the parental populations, statistically significant distances were found. The interpopulation genetic distance (F(ST) coefficients) to the present database ranged from F(ST)=0.0431 (p<0.00001) between Macapá and Belém to F(ST)=0.266 (p<0.00001) between Macapá and the Native American group.     

Applicability of the SNPforID 52-plex panel for human identification and ancestry evaluation in a Brazilian population sample by next-generation sequencing

SNP analysis is of paramount importance in forensic genetics. The development of new technologies in next-generation sequencing allowed processing a large number of markers in various samples simultaneously. Although SNPs are less informative than STRs, they present lower mutation rates and perform better when using degraded samples. Some SNP systems were developed for forensic usage, such as the SNPforID 52-plex, from the SNPforID Consortium, containing 52 bi-allelic SNPs for human identification. In this paper we evaluated the informativeness of this system in a Brazilian population sample (n = 340). DNA libraries were prepared using a customized HaloPlex Target Enrichment System kit (Agilent Technologies, Inc.) and sequenced in the MiSeq Personal Sequencer platform (Illumina Inc.). The methodology presented here allowed the analysis of 51 out of 52 SNPforID markers. Allele frequencies and forensic parameters were estimated, revealing high informativeness: the combined match probability and power of exclusion were 6.48 × 10⁻²¹ and 0.9997, respectively. Population admixture analysis indicates high European contribution (more than 70%) and low Amerindian contribution (less than 10%) in our population, while individual admixture analyses were consistent with the majority of individuals presenting high European contribution. This study demonstrates that the 52-plex kit is suitable for forensic cases in a Brazilian population, presenting results comparable with those obtained using a 16 STR panel.     

Mitochondrial DNA diversity in a population from Santa Catarina (Brazil): predominance of the European input

The state of Santa Catarina (Brazil) is known to have represented a cultural crossroads in South America due to several historic migrations mainly from Europe and Africa. We set out to scrutinize whether the genetic imprint of these migrations could be traced through analysis of the matrilineal gene pool of the Catarinenses. The entire control region of the mitochondrial DNA was studied in 80 healthy and maternally unrelated individuals. The analysis of haplogroup distribution revealed that this population is extremely heterogeneous, showing the coexistence of matrilineal lineages with three different phylogeographic origins. European lineages are the most frequent due mainly to the impact of relatively recent migratory waves from Europe. In spite of this, Native American lineages and African lineages incorporated with the slave trade are also present in noticeable proportions. The strikingly high variability generated by intense gene flow is mirrored in a high sequence diversity (0.9930) and power of discrimination (0.9806). Thus, analysis of the entire mitochondrial DNA control region emerges as a valuable tool for forensic genetic purposes in this highly admixed population, an attribute common to several present-day Latin American populations.     

Population data for 12 Y-chromosome STR loci in a sample from El Salvador

Haplotype, allele frequencies and population data of 12 Y-chromosome STR loci DYS19, DYS385, DYS389I, DYS389II, DYS390, DYS391, DYS392, DYS393, DYS437, DYS438 and DYS439 were determined from a sample of 150 unrelated male individuals from El Salvador, Central America. A total of 131 haplotypes were identified by the 12 Y-STR loci of which 118 were unique. The haplotype diversity (99.08%) and the proportion of different haplotypes (87.33%) were estimated. R_ST genetic distances were calculated between El Salvador and other populations from Southern and Central America, Europe and Africa. The highest R_ST genetic distances were found when comparing El Salvador with African populations (0.334 ? R_ST ? 0.395). The lowest non-significant distance was found in the comparison with Honduras. The observed genetic distances between El Salvador and Southern and Central Native groups presented a wide range of values (from 0.024 to 0.210) that can be explained by the differences in the proportion of European versus Amerindian contributions in these population groups. The Multi Dimensional Scaling (MDS) plot analysis, based on pairwise R_ST values, showed that the general population of El Salvador is closer to the European cluster (composed by European and South American general population samples from Brazil, Argentina, Colombia and Venezuela) than to the Southern/Central American cluster of Native and Mestizo populations.     

Inferring ancestral origin using a single multiplex assay of ancestry-informative marker SNPs

Tests that infer the ancestral origin of a DNA sample have considerable potential in the development of forensic tools that can help to guide crime investigation. We have developed a single-tube 34-plex SNP assay for the assignment of ancestral origin by choosing ancestry-informative markers (AIMs) exhibiting highly contrasting allele frequency distributions between the three major population-groups. To predict ancestral origin from the profiles obtained, a classification algorithm was developed based on maximum likelihood. Sampling of two populations each from African, European and East Asian groups provided training sets for the algorithm and this was tested using the CEPH Human Genome Diversity Panel. We detected negligible theoretical and practical error for assignments to one of the three groups analyzed with consistently high classification probabilities, even when using reduced subsets of SNPs. This study shows that by choosing SNPs exhibiting marked allele frequency differences between population-groups a practical forensic test for assigning the most likely ancestry can be achieved from a single multiplexed assay.     

The maternal inheritance of Alto Paraná revealed by full mitogenome sequences

Most studies on maternal lineages of South America populations are restricted to control region (CR) markers and, for some geographical regions, the number of studied samples does not adequately represent the existing diversity. This is the case of mitochondrial DNA (mtDNA) studies on Paraguay that are limited to two Native ethnic groups. To overcome this deficiency, we analysed the mitogenomes from 105 individuals living in Alto Paraná, the second most populated department of the country. Using the Precision ID mtDNA Whole Genome Panel, the molecule was sequenced on Ion S5. The majority of the haplotypes belong to the Native American lineages A, B, C and D. Analyses of maximum parsimony using mitogenome data retrieved from publications and in The 1000 Genomes Project showed a high number of new native American subclades in Paraguay. Also, none of the haplotypes found in Alto Paraná match the remaining South American samples, which include admixed populations from Colombia, Peru and Ecuador, and natives from Colombia and Ecuador. F_ST genetic distance analysis showed that the native genetic background of Alto Paraná has an intermediate position between the Amazonian groups and the admixed populations from Peru and Ecuador, supporting the theory about the Amazonian origin of the Tupi-Guarani and, at the same time, showing the influence of other linguistic groups.     

Distribution and haplotype analysis of all the non-synonymous and autoimmunity-related single nucleotide polymorphisms in the human deoxyribonuclease II gene using worldwide populations

We have focused on the 14 SNPs including all the non-synonymous and autoimmunity-related ones in the DNase II gene (DNASE2). The distribution of each allele and haplotype in these SNPs was examined in eight Asian, three African, three Mexican and two Caucasian populations using the newly developed PCR–RFLP methods. Eight SNPs among nine non-synonymous ones were monomorphic, indicating that a specific allele generating the intact activity-harboring DNase II in these SNPs is well conserved in worldwide populations. On the other hand, five other SNPs (?1951G>A, ?1066G>C, ?390A>C, +2630T>C, and +6235G>C) related to autoimmunity exhibited polymorphism common in worldwide populations, and especially their distributions were ethnic-dependent in the same manner as those of haplotypes. Furthermore, a strong linkage between SNPs ?1951G>A and ?1066G>C was confirmed in most populations. This study was the first to report any worldwide population analysis regarding all the non-synonymous and autoimmunity-related SNPs in the DNASE2, providing genetic information on the DNASE2 as a genetic marker for personal identification and/or genetic factor for susceptibility to autoimmunity.