The peopling of South America and the trans-Andean gene flow of the first settlers.

Genetic and archaeological data indicate that the initial Paleoindian settlers of South America followed two entry routes separated by the Andes and the Amazon rainforest. The interactions between these paths and their impact on the peopling of South America remain unclear. Analysis of genetic variation in the Peruvian Andes and regions located south of the Amazon River might provide clues on this issue. We analyzed mitochondrial DNA variation at different Andean locations and >360,000 autosomal SNPs from 28 Native American ethnic groups to evaluate different trans-Andean demographic scenarios. Our data reveal that the Peruvian Altiplano was an important enclave for early Paleoindian expansions and point to a genetic continuity in the Andes until recent times, which was only marginally affected by gene flow from the Amazonian lowlands. Genomic variation shows a good fit with the archaeological evidence, indicating that the genetic interactions between the descendants of the settlers that followed the Pacific and Atlantic routes were extremely limited.

Evaluation of forensic DNA mixture evidence: protocol for evaluation, interpretation, and statistical calculations using the combined probability of inclusion.

BACKGROUND: The evaluation and interpretation of forensic DNA mixture evidence faces greater interpretational challenges due to increasingly complex mixture evidence. Such challenges include: casework involving low quantity or degraded evidence leading to allele and locus dropout; allele sharing of contributors leading to allele stacking; and differentiation of PCR stutter artifacts from true alleles. There is variation in statistical approaches used to evaluate the strength of the evidence when inclusion of a specific known individual(s) is determined, and the approaches used must be supportable. There are concerns that methods utilized for interpretation of complex forensic DNA mixtures may not be implemented properly in some casework. Similar questions are being raised in a number of U.S. jurisdictions, leading to some confusion about mixture interpretation for current and previous casework. RESULTS: Key elements necessary for the interpretation and statistical evaluation of forensic DNA mixtures are described. Given the most common method for statistical evaluation of DNA mixtures in many parts of the world, including the USA, is the Combined Probability of Inclusion/Exclusion (CPI/CPE). Exposition and elucidation of this method and a protocol for use is the focus of this article. Formulae and other supporting materials are provided. CONCLUSIONS: Guidance and details of a DNA mixture interpretation protocol is provided for application of the CPI/CPE method in the analysis of more complex forensic DNA mixtures. This description, in turn, should help reduce the variability of interpretation with application of this methodology and thereby improve the quality of DNA mixture interpretation throughout the forensic community.

The mitochondrial DNA history of a former native American village in northern Uruguay.

OBJECTIVES: In 1828, between 8,000 and 15,000 Indians from the Jesuit Missions were brought to Uruguay. There, they were settled in a village, presently named Bella Unión, in the northwest corner of the country. According to historic sources, the Indians abandoned the settlement shortly thereafter, with the village subsequently repopulated by "criollos" and immigrants from abroad. As a first approach to reconstruct the genetic history of the population, data about the living population genetic structure will be used. Based on the analysis of the maternal lineages of the inhabitants of Bella Unión, and of those from two nearby villages, we expect to partially answer what happened with the first and subsequent inhabitants. METHODS: We analyzed the maternal lineages of the present inhabitants of Bella Unión and neighboring localities through the sequencing of the mitochondrial DNA control region. RESULTS: A total of 64.3%, 5.7%, and 30% of the mtDNAs were of Native, African, and West Eurasian origin, respectively. These figures are quite similar to that of the population of Tacuarembó, which is located in northeastern Uruguay. The four main Native American founding haplogroups were detected, with B2 being the most frequent, while some rare subhaplogroups (B2h, C1b2, D1f1) were also found. When compared with other Native American sequences, near- matches most consistently pointed to an Amazonian Indian origin which, when considered with historical evidence, suggested a probable Guaraní-Missionary-related origin. CONCLUSIONS: The data support the existence of a relationship between the historic and present inhabitants of the extreme northwest Uruguay, with a strong contribution of Native Americans to the mitochondrial DNA diversity observed there.

Toward male individualization with rapidly mutating y-chromosomal short tandem repeats.

Relevant for various areas of human genetics, Y-chromosomal short tandem repeats (Y-STRs) are commonly used for testing close paternal relationships among individuals and populations, and for male lineage identification. However, even the widely used 17-loci Yfiler set cannot resolve individuals and populations completely. Here, 52 centers generated quality-controlled data of 13 rapidly mutating (RM) Y-STRs in 14,644 related and unrelated males from 111 worldwide populations. Strikingly, >99% of the 12,272 unrelated males were completely individualized. Haplotype diversity was extremely high (global: 0.9999985, regional: 0.99836-0.9999988). Haplotype sharing between populations was almost absent except for six (0.05%) of the 12,156 haplotypes. Haplotype sharing within populations was generally rare (0.8% nonunique haplotypes), significantly lower in urban (0.9%) than rural (2.1%) and highest in endogamous groups (14.3%). Analysis of molecular variance revealed 99.98% of variation within populations, 0.018% among populations within groups, and 0.002% among groups. Of the 2,372 newly and 156 previously typed male relative pairs, 29% were differentiated including 27% of the 2,378 father-son pairs. Relative to Yfiler, haplotype diversity was increased in 86% of the populations tested and overall male relative differentiation was raised by 23.5%. Our study demonstrates the value of RM Y-STRs in identifying and separating unrelated and related males and provides a reference database.

Comparison of the performance of different models for the interpretation of low level mixed DNA profiles.

DNA analyses from forensic casework samples commonly result in complex DNA profiles. Often, these profiles consist of multiple contributors and display multiple stochastic events such as peak height imbalance, allelic or locus drop-out, allelic drop-in, and excessive or indistinguishable stutter. This increased complexity has established a need for more sophisticated methods of DNA mixture interpretation. This study compares the effectiveness of statistical models in the interpretation of artificially created low template two person mixed DNA profiles at varying proportions and template quantities. Two binary models (combined probability of inclusion and random match probability), a semicontinuous (Lab retriever), and continuous model (STRmix™) were compared. Generally, as the sophistication of the models increases, the power of discrimination increases. Differences in discrimination often correlate to each model's ability to use observed data effectively. Binary models require static thresholds resulting in unused data and outliers that may lead to difficult or incorrect interpretation. Semicontinuous and continuous models eliminate the stochastic threshold, however Lab Retriever does not account for stochastic events beyond drop-out and drop-in leading to possible less effective use of the data. STRmix™ incorporates all stochastic events listed above into the calculation making the most effective use of the observed data.

Mutation rates of 15 X chromosomal short tandem repeat markers.

Though allele frequency data for a variety of X chromosomal short tandem repeat (STR) markers in a range of populations have been reported, fewer studies of mutation rates in these same markers or populations are available. In order to address possible mismatches during kinship analysis due to mutation, a robust estimate of the rate of mutation must be established. Here, mutation rates in three US populations have been determined for a total of 15 markers (DXS6789, DXS9902, DXS7132, DXS7130, DXS6795, DXS10147, DXS8378, DXS7423, HPRTB, DXS101, DXS7424, GATA31E08, GATA172D05, GATA165B12, and DXS6803). Eighteen mutations over 20,625 meioses were observed, and the overall X STR mutation rate in this study was found to be 8.73 × 10(-4) (95 % CI, 5.2-13.8 × 10(-4)). A review of published mutation rate studies revealed similar findings in other global populations, and allowed the compilation of a combined dataset of 81,310 meioses which can be employed by the forensic community.

Mitochondrial control region variation in a Korean population sample.

Mitochondrial DNA (mtDNA) control region (16024-576) sequences were generated from 281 individuals from South Korea. Robotic liquid handling, a redundant sequencing strategy, and a series of quality control checks were implemented to ensure the high quality of the dataset. This population sample showed a low random match probability (0.25 %) and high genetic diversity (0.9933). The haplogroup breakdown was consistent with previous studies describing Korean mtDNA variation. The 224 unique haplotypes (33 shared) presented will supplement the data already publically available.

Mixture detection with Demixtify.

The de facto genetic markers of forensics are short tandem repeats (STRs). There are many analytical tools designed to work with STRs, including techniques for analyzing and assessing DNA mixtures. In contrast, the nascent field of forensic genetic genealogy often relies on biallelic single nucleotide polymorphisms (SNPs). Tools designed for the forensic assessment of SNPs are somewhat lacking, especially for DNA mixtures. In this paper we introduce Demixtify, a program that detects DNA mixtures using biallelic SNPs. Demixtify is quite powerful; highly imbalanced mixtures can be detected (≤1:99, considering in silico and in vitro mixtures) when coverage is ample. Demixtify can also detect mixtures in low coverage (∼1×) samples (when the mixture is relatively balanced). Demixtify includes an empirical estimator of sequence error that is specific to the markers assayed, making it especially relevant to the forensic community. Orthogonal techniques are also developed to characterize in vitro mixtures, as well as samples thought to be single source, and the results of these approaches serve to validate the techniques presented.

Identifying distant relatives using benchtop-scale sequencing.

The genetic component of forensic genetic genealogy (FGG) is an estimate of kinship, often conducted at genome scales between a great number of individuals. The promise of FGG is substantial: in concert with genealogical records and other nongenetic information, it can indirectly identify a person of interest. A downside of FGG is cost, as it is currently expensive and requires chemistries uncommon to forensic genetic laboratories (microarrays and high throughput sequencing). The more common benchtop sequencers can be coupled with a targeted PCR assay to conduct FGG, though such approaches have limited resolution for kinship. This study evaluates low-pass sequencing, an alternative strategy that is accessible to benchtop sequencers and can produce resolutions comparable to high-pass sequencing. Samples from a three-generation pedigree were augmented to include up to 7th degree relatives (using whole genome pedigree simulations) and the ability to recover the true kinship coefficient was assessed using algorithms qualitatively similar to those found in GEDmatch. We show that up to 7th degree relatives can be reliably inferred from 1 × whole genome sequencing obtainable from desktop sequencers.

Demographic expansions in South America: enlightening a complex scenario with genetic and linguistic data.

Native Americans are characterized by specific and unique patterns of genetic and cultural/linguistic diversities, and this information has been used to understand patterns of geographic dispersion, and the relationship between these peoples. Particularly interesting are the Tupi and Je speaker dispersions. At present, a large number of individuals speak languages of these two stocks; for instance, Tupi-Guarani is one of the official languages in Paraguay, Bolivia, and the Mercosul economic block. Although the Tupi expansion can be compared in importance to the Bantu migration in Africa, little is known about this event relative to others. Equal and even deeper gaps exist concerning the Je-speakers' expansion. This study aims to elucidate some aspects of these successful expansions. To meet this purpose, we analyzed Native American mtDNA complete control region from nine different populations and included HVS-I sequences available in the literature, resulting in a total of 1,176 samples investigated. Evolutionary relationships were explored through median-joining networks and genetic/geographic/linguistic correlations with Mantel tests and spatial autocorrelation analyses. Both Tupi and Je showed general traces of ancient or more recent fission-fusion processes, but a very different pattern of demographic expansion. Tupi populations displayed a classical isolation-by-distance pattern, while Je groups presented an intricate and nonlinear mode of dispersion. We suggest that the collective memory and other cultural processes could be important factors influencing the fission-fusion events, which likely contributed to the genetic structure, evolution, and dispersion of Native American populations.

Study of CTS DNA Proficiency Tests with Regard to DNA Mixture Interpretation: A NIST Scientific Foundation Review.

The National Institute of Standards and Technology has released a document entitled DNA Mixture Interpretation: A NIST Scientific Foundation Review for public comment. This has become known as the Draft NIST Foundation Review. It contains the statement: "Across these 69 data sets, there were 80 false negatives and 18 false positives reported from 110,408 possible responses (27,602 participants × two evidence items × two reference items). In the past five years, the number of participants using PGS has grown." We examine a set of proficiency test results to determine if these NIST statements could be justified. The summary reports for each relevant forensic biology test (Forensic Biology, Semen, and Mixture) in the years 2018-2021 were reviewed. Data were also provided to us by CTS upon our request. None of the false positives or negatives could be attributed to the mixture interpretation strategy and certainly not to the use of PGS.

Assessment of human nuclear and mitochondrial DNA qPCR assays for quantification accuracy utilizing NIST SRM 2372a.

In forensic DNA casework, a highly accurate real-time quantitative polymerase chain reaction (qPCR) assay is recommended per the Scientific Working Group on DNA Analysis Methods (SWGDAM) (SWGDAM Validation Guidelines for DNA Analysis Methods [1]) to determine whether a DNA sample is of sufficient quantity and robust quality to move forward with downstream short tandem repeats (STR) or sequencing analyses. Most of these assays rely on a standard curve, referred to herein and traditionally as absolute qPCR, in which an unknown is compared, relative to that curve. However, one fundamental issue with absolute qPCR is the quantifiable concentration of commercial assay standards can vary depending on (1) origin, i.e., whether from a cell line or a human subject, (2) supplier, (3) lot number, (4) shipping method, etc. In 2018, the National Institute for Standards and Technology (NIST) released a human DNA standard reference material for evaluating qPCR quantification standards, Standard Reference Material (SRM) 2372a, Romsos et al. (2018) [2] which contains three well-characterized human genomic DNA samples: Component A) a single male1 donor, Component B) a single female1 donor, and Component C) a 1:3 male2:female2 donor, each with certification data for nDNA and informational mitochondrial DNA(mtDNA)/nuclear DNA (nDNA) ratio data. The SRM 2372a was used to assess four qPCR assays: (1) Quantifiler Trio (Thermo Fisher Scientific, Waltham, MA) for nDNA quantification, (2) NovaQUANT (EMD Millipore Corporation, San Diego, CA) for nDNA and mtDNA quantification, (3) a custom duplex mtDNA assay, and (4) a custom triplex mtDNA assay. Additionally, extracts from eighteen (18) skeletal remains were tested with the latter three assays for concordance of DNA concentration and with assays (2) and (3), for the degradation state. Our assessment revealed that an accurate, efficient, and reproducible qPCR assay is dependent on (1) the quality and reliability of the DNA standard, (2) the qPCR chemistry, and (3) the specific primers, and probes (if applicable), used in an assay. Our findings indicate qPCR assays may not always quantify as expected and that performance of each lot should be verified using a well-characterized DNA standard such as the NIST SRM 2372a and adjusted if warranted.

A new implementation of a semi-continuous method for DNA mixture interpretation.

A new calculation module within the PopStats module of the CODIS software package, based on the underlying mathematics presented in the MixKin software package, has been developed for assigning the Likelihood Ratio (LR) of DNA mixture profiles. This module uses a semi-continuous model that allows for population structure and allelic drop-out and drop-in but does not require allelic peak heights or other laboratory-specific parameters. This new implementation (named SC Mixture), like MixKin, does not specify or estimate a probability of drop-out. Instead, each contributor to a mixture has an independent drop-out rate, and the probability of the mixture profile for a specified proposition concerning the contributors is integrated over the range of possible drop-out rates. The allelic drop-in rate and the population structure parameter, theta, used by the software are specified by the user. The user can examine up to five contributors to a mixture, however, conditioning on assumed contributors and limiting the number of unknowns in both numerator and denominator hypotheses greatly improves performance. We report results from an extensive validation study performed for ten mixtures with each of one (single source), two, three, four, or five contributors, with four combinations of drop-in rate and a population structure parameter. Each mixture was run as a complete profile or with the random removal of alleles to simulate drop-out. All 1620 combinations were evaluated with PopStats, MixKin, and LRmix and considerable consistency was found among the results with all three packages.

Evaluating self-declared ancestry of U.S. Americans with autosomal, Y-chromosomal and mitochondrial DNA.

The current U.S. population represents an amalgam of individuals originating mainly from four continental regions (Africa, Europe, Asia and America). To study the genetic ancestry and compare with self-declared ancestry we have analyzed paternally, maternally and bi-parentally inherited DNA markers sensitive for indicating continental genetic ancestry in all four major U.S. American groups. We found that self-declared U.S. Hispanics and U.S. African Americans tend to show variable degrees of continental genetic admixture among the three genetic systems, with evidence for a marked sex-biased admixture history. Moreover, for these two groups we observed significant regional variation across the country in genetic admixture. In contrast, self-declared U.S. European and U.S. Asian Americans were genetically more homogeneous at the continental ancestry level. Two autosomal ancestry-sensitive markers located in skin pigmentation candidate genes showed significant differences in self-declared U.S. African Americans or U.S. European Americans, relative to their assumed parental populations from Africa or Europe. This provides genetic support for the importance of skin color in the complex process of ancestry identification.

Homogeneity in mitochondrial DNA control region sequences in Swedish subpopulations.

In order to promote mitochondrial DNA (mtDNA) testing in Sweden we have typed 296 Swedish males, which will serve as a Swedish mtDNA frequency database. The tested males were taken from seven geographically different regions representing the contemporary Swedish population. The complete mtDNA control region was typed and the Swedish population was shown to have high haplotype diversity with a random match probability of 0.5%. Almost 47% of the tested samples belonged to haplogroup H and further haplogroup comparison with worldwide populations clustered the Swedish mtDNA data together with other European populations. AMOVA analysis of the seven Swedish subregions displayed no significant maternal substructure in Sweden (F (ST) = 0.002). Our conclusion from this study is that the typed Swedish individuals serve as good representatives for a Swedish forensic mtDNA database. Some caution should, however, be taken for individuals from the northernmost part of Sweden (provinces of Norrbotten and Lapland) due to specific demographic conditions. Furthermore, our analysis of a small sample set of a Swedish Saami population confirmed earlier findings that the Swedish Saami population is an outlier among European populations.

Evaluation of modified Yfiler amplification strategy for compromised samples.

AIM: To characterize the data produced using a modified amplification protocol for the AmpFlSTR Yfiler PCR Amplification Kit (Applied Biosystems) and explore the potential of Y-chromosomal short tandem repeat (Y-STR) recovery from severely degraded skeletal remains encountered at the Armed Forces DNA Identification Laboratory. METHODS: Experiments were performed using two sets of Yfiler amplification parameters. One set of parameters reflected the manufacturer's recommendations. The second set of parameters included twice the recommended Taq concentration and 6 additional cycles. Recovery of authentic alleles and the incidence of drop-in alleles were assessed for 3 data sets: 8 different quantities of pristine DNA, 8 artificially-degraded samples, and 31 non-probative case samples. RESULTS: Samples tested with both protocols from all 3 data sets yielded twice as many authentic alleles under the modified parameters than under the standard parameters (62% vs 31%), with only a nominal associated increase in the occurrence of non-authentic alleles (1.36% of all alleles detected). When applied to a range of representative casework samples, the modified protocol leveraged 9 or more reproducible alleles from over half of the specimens tested. CONCLUSION: Reproducible and informative Y-STR profiles can be recovered from a broad range of degraded and inhibited skeletal remains extracts when a commercially available kit is employed under modified amplification parameters.

Probabilistic genotyping software: An overview.

The interpretation of mixed profiles from DNA evidentiary material is one of the more challenging duties of the forensic scientist. Traditionally, analysts have used a "binary" approach to interpretation where inferred genotypes are either included or excluded from the mixture using a stochastic threshold and other biological parameters such as heterozygote balance, mixture ratio, and stutter ratios. As the sensitivity of STR multiplexes and capillary electrophoresis instrumentation improved over the past 25 years, coupled with the change in the type of evidence being submitted for analysis (from high quality and quantity (often single-source) stains to low quality and quantity (often mixed) "touch" samples), the complexity of DNA profile interpretation has equally increased. This review provides a historical perspective on the movement from binary methods of interpretation to probabilistic methods of interpretation. We describe the two approaches to probabilistic genotyping (semi-continuous and fully continuous) and address issues such as validation and court acceptance. Areas of future needs for probabilistic software are discussed.

STRmix™ collaborative exercise on DNA mixture interpretation.

An intra and inter-laboratory study using the probabilistic genotyping (PG) software STRmix™ is reported. Two complex mixtures from the PROVEDIt set, analysed on an Applied Biosystems™ 3500 Series Genetic Analyzer, were selected. 174 participants responded. For Sample 1 (low template, in the order of 200 rfu for major contributors) five participants described the comparison as inconclusive with respect to the POI or excluded him. Where LRs were assigned, the point estimates ranging from 2 × 104 to 8 × 106. For Sample 2 (in the order of 2000 rfu for major contributors), LRs ranged from 2 × 1028 to 2 × 1029. Where LRs were calculated, the differences between participants can be attributed to (from largest to smallest impact): This study demonstrates a high level of repeatability and reproducibility among the participants. For those results that differed from the mode, the differences in LR were almost always minor or conservative.

Characterization of 26 miniSTR loci for improved analysis of degraded DNA samples.

An additional 20 novel mini-short tandem repeat (miniSTR) loci have been developed and characterized beyond the six previously developed by our laboratory for a total of 26 non-CODIS miniSTR markers. These new markers produce short PCR products in the target range of 50-150 base pairs (bp) by moving the primer sequences as close as possible-often directly next to the identified repeat region. These candidate loci were initially screened based on their small amplicon sizes and locations on chromosomes currently unoccupied by the 13 CODIS STR loci or at least 50 Mb away from them on the same chromosome. They were sequenced and evaluated across more than 600 samples, and their population statistics were determined. The heterozygosities of the new loci were compared with those of the 13 CODIS loci and all were found to be comparable. Only five of the new loci had lower values than the CODIS loci; however, all of these were much smaller in size. This data suggests that these 26 miniSTR loci will serve as useful complements to the CODIS loci to aid in the forensic analysis of degraded DNA, as well as missing persons work and parentage testing with limited next-of-kin reference samples.