Robert F. Oldt

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.

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.

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.

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.

SNP-based genetic characterization of the Tulane National Primate Research Center's conventional and specific pathogen-free rhesus macaque (Macaca mulatta) populations

The rhesus macaque is an important biomedical model organism, and the Tulane National Primate Research Center (TNPRC) has one of the largest rhesus macaque breeding colonies in the United States.

Determination of major histocompatibility class I and class II genetic composition of the Caribbean Primate Center specific pathogen-free rhesus macaque (Macaca mulatta) colony based on massively parallel sequencing

Knowledge of major histocompatibility complex (MHC) composition and distribution in rhesus macaque colonies is critical for management strategies that maximize the utility of this model for biomedical research.

Ancestry, Plasmodium cynomolgi prevalence and rhesus macaque admixture in cynomolgus macaques ( Macaca fascicularis ) bred for export in Chinese breeding farms

Most cynomolgus macaques (Macaca fascicularis) used in the United States as animal models are imported from Chinese breeding farms without documented ancestry. Cynomolgus macaques with varying rhesus macaque ancestry proportions may exhibit differences, such as susceptibility to malaria, that affect their suitability as a research model.

ABO Blood Group Phenotype Frequency Estimation Using Molecular Phenotyping in Rhesus and Cynomolgus Macaques

A much larger sample (N = 2369) was used to evaluate a previously reported distribution of the A, AB and B blood group phenotypes in rhesus and cynomolgus macaques from six different regional populations. These samples, acquired from 15 different breeding and research facilities in the United States, were analyzed using a real‐time quantitative polymerase chain reaction (qPCR) assay that targets single nucleotide polymorphisms (SNPs) responsible for the macaque A, B and AB phenotypes. The frequency distributions of blood group phenotypes of the two species differ significantly from each other and significant regional differentiation within the geographic ranges of each species was also observed. The B blood group phenotype was prevalent in rhesus macaques, especially those from India, while the frequencies of the A, B and AB phenotypes varied significantly among cynomolgus macaques from different geographic regions. The Mauritian cynomolgus macaques, despite having originated in Indonesia, showed significant (P ≪ .01) divergence from the Indonesian animals at the ABO blood group locus. Most Mauritian animals belonged to the B blood group while the Indonesian animals were mostly A. The close similarity in blood group frequency distributions between the Chinese rhesus and Indochinese cynomolgus macaques demonstrates that the introgression between these two species extends beyond the zone of intergradation in Indochina. This study underscores the importance of ABO blood group phenotyping of the domestic supply of macaques and their biospecimens.

Genetic analysis of samples from wild populations opens new perspectives on hybridization between long-tailed (Macaca fascicularis) and rhesus macaques (Macaca mulatta)

In the past decade, many researchers have published papers about hybridization between long‐tailed and rhesus macaques. These previous works have proposed unidirectional gene flow with the Isthmus of Kra as the zoogeographical barrier of hybridization. However, these reports analyzed specimens of unknown origin and/or did not include specimens from Thailand, the center of the proposed area of hybridization. Collected specimens of long‐tailed and rhesus macaques representing all suspected hybridization areas were examined. Blood samples from four populations each of long‐tailed and rhesus macaques inhabiting Thailand, Myanmar, and Laos were collected and analyzed with conspecific references from China (for rhesus macaques) and multiple countries from Sundaic regions (for long‐tailed macaques). Ninety‐six single nucleotide polymorphism (SNP) markers specifically designed to interrogate admixture and ancestry were used in genotyping. We found genetic admixture maximized at the hybrid zone (15–20°N), as well as admixture signals of varying strength in both directions outside of the hybrid zone. These findings show that the Isthmus of Kra is not a barrier to gene flow from rhesus to long‐tailed populations. However, to precisely identify a southernmost barrier, if in fact a boundary rather than simple isolation by distance exists, the samples from peninsular Malaysia must be included in the analysis. Additionally, a long‐tailed to rhesus gene flow boundary was found between northern Thailand and Myanmar. Our results suggest that selection of long‐tailed and rhesus macaques, the two most commonly used non‐human primates for biomedical research, should take into account not only the species identification but also the origin of and genetic admixture within and between the species.