Genetic differentiation between and within Northern Native American language groups: an argument for the expansion of the Native American CODIS database.

The National Research Council recommends that genetic differentiation among subgroups of ethnic samples be lower than 3% of the total genetic differentiation within the ethnic sample to be used for estimating reliable random match probabilities for forensic use. Native American samples in the United States' Combined DNA Index System (CODIS) database represent four language families: Algonquian, Na-Dene, Eskimo-Aleut, and Salishan. However, a minimum of 27 Native American language families exists in the US, not including language isolates. Our goal was to ascertain whether genetic differences are correlated with language groupings and, if so, whether additional language families would provide a more accurate representation of current genetic diversity among tribal populations. The 21 short tandem repeat (STR) loci included in the Globalfiler® PCR Amplification Kit were used to characterize six indigenous language families, including three of the four represented in the CODIS database (i.e. Algonquian, Na-Dene, and Eskimo-Aleut), and two language isolates (Miwok and Seri) using major population genetic diversity metrics such as F statistics and Bayesian clustering analysis of genotype frequencies. Most of the genetic variation (97%) was found to be within language families instead of among them (3%). In contrast, when only the three of the four language families represented in both the CODIS database and the present study were considered, 4% of the genetic variation occurred among the language groups. Bayesian clustering resulted in a maximum posterior probability indicating three genetically distinct groups among the eight language families and isolates: (1) Eskimo, (2) Seri, and (3) all other language groups and isolates, thus confirming genetic subdivision among subgroups of the CODIS Native American database. This genetic structure indicates the need for an increased number of Native American populations based on language affiliation in the CODIS database as well as more robust sample sets for those language families. Supplemental data for this article is available online at https://doi.org/10.1080/20961790.2021.1963088 .

Flanking region variation of ForenSeq™ DNA Signature Prep Kit STR and SNP loci in Yavapai Native Americans.

Massively parallel sequencing (MPS) offers advantages over current capillary electrophoresis-based analysis of short tandem repeat (STR) loci for human identification testing. In particular STR repeat motif sequence information can be obtained, thereby increasing the discrimination power of some loci. While sequence variation within the repeat region is observed relatively frequently in some of the commonly used STRs, there is an additional degree of variation found in the flanking regions adjacent to the repeat motif. Repeat motif and flanking region sequence variation have been described for major population groups, however, not for more isolated populations. Flanking region sequence variation in STR and single nucleotide polymorphism (SNP) loci in the Yavapai population was analyzed using the ForenSeq™ DNA Signature Prep Kit and STRait Razor v2s. Seven and 14 autosomal STRs and identity-informative single nucleotide polymorphisms (iiSNPs), respectively, had some degree of flanking region variation. Three and four of these identity-informative loci, respectively, showed ≥5% increase in expected heterozygosity. The combined length- and sequence-based random match probabilities (RMPs) for 27 autosomal STRs were 6.11×10-26 and 2.79×10-29, respectively. When combined with 94 iiSNPs (a subset of which became microhaplotypes) the combined RMP was 5.49×10-63. Analysis of length-based and sequence-based autosomal STRs in STRUCTURE indicated that the Yavapai are most similar to the Hispanic population. While producing minimal increase in X- and Y-STR discrimination potential, access to flanking region data enabled identification of one novel X-STR and three Y-STR alleles relative to previous reports. Five ancestry-informative SNPs (aiSNPs) and two phenotype-informative SNPs (piSNPs) exhibited notable flanking region variation.

The genetic structure of native Americans in North America based on the Globalfiler® STRs.

Current forensic STR databases, such as CODIS, lack population genetic data on Native American populations. Information from a geographically diverse array of tribes is necessary to provide improved statistical estimates of the strength of associations with DNA evidence. The Globalfiler® STR markers were used to characterize the genetic structure of ten tribal populations from seven geographic regions in North America, including those not presently represented in forensic databases. Samples from the Arctic region, Baja California, California/Great Basin, the Southeast, Mexico, the Midwest, and the Southwest were analyzed for allele frequencies, observed and expected heterozygosities, and F-statistics. The tribal samples exhibited an FST or θ value above the conservative 0.03 estimate recommended by the National Research Council (NRC) for calculating random match probabilities among Native Americans. The greater differentiation among tribal populations computed here (θ=0.04) warrants the inclusion of additional regional Native American samples into STR databases.

Native American population data based on the Globalfiler(®) autosomal STR loci.

Native American population data are limited and thus impact computing accurate statistical parameters for forensic investigations. Thus, additional information should be generated from geographically representative tribes in North America, particularly from those that are not included in existing population databases for forensic use. The Globafiler(®) PCR Amplification kit was used to produce STR genotypic data for 533 individuals who represent 31 Native American tribal populations derived from eight geographically diverse regions in North America. Population genetic estimates from 21 autosomal STRs are reported.

Genetic analysis of the Yavapai Native Americans from West-Central Arizona using the Illumina MiSeq FGx™ forensic genomics system.

Forensically-relevant genetic markers were typed for sixty-two Yavapai Native Americans using the ForenSeq™ DNA Signature Prep Kit.These data are invaluable to the human identity community due to the greater genetic differentiation among Native American tribes than among other subdivisions within major populations of the United States. Autosomal, X-chromosomal, and Y-chromosomal short tandem repeat (STR) and identity-informative (iSNPs), ancestry-informative (aSNPs), and phenotype-informative (pSNPs) single nucleotide polymorphism (SNP) allele frequencies are reported. Sequence-based allelic variants were observed in 13 autosomal, 3 X, and 3 Y STRs. These observations increased observed and expected heterozygosities for autosomal STRs by 0.081±0.068 and 0.073±0.063, respectively, and decreased single-locus random match probabilities by 0.051±0.043 for 13 autosomal STRs. The autosomal random match probabilities (RMPs) were 2.37×10-26 and 2.81×10-29 for length-based and sequence-based alleles, respectively. There were 22 and 25 unique Y-STR haplotypes among 26 males, generating haplotype diversities of 0.95 and 0.96, for length-based and sequencebased alleles, respectively. Of the 26 haplotypes generated, 17 were assigned to haplogroup Q, three to haplogroup R1b, two each to haplogroups E1b1b and L, and one each to haplogroups R1a and I1. Male and female sequence-based X-STR random match probabilities were 3.28×10-7 and 1.22×10-6, respectively. The average observed and expected heterozygosities for 94 iSNPs were 0.39±0.12 and 0.39±0.13, respectively, and the combined iSNP RMP was 1.08×10-32. The combined STR and iSNP RMPs were 2.55×10-58 and 3.02×10-61 for length-based and sequence-based STR alleles, respectively. Ancestry and phenotypic SNP information, performed using the ForenSeq™ Universal Analysis Software, predicted black hair, brown eyes, and some probability of East Asian ancestry for all but one sample that clustered between European and Admixed American ancestry on a principal components analysis. These data serve as the first population assessment using the ForenSeq™ panel and highlight the value of employing sequence-based alleles for forensic DNA typing to increase heterozygosity, which is beneficial for identity testing in populations with reduced genetic diversity.