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.

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.

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 × 10⁴ to 8 × 10⁶. For Sample 2 (in the order of 2000 rfu for major contributors), LRs ranged from 2 × 10²⁸ to 2 × 10²⁹. 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.

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.

A response to "Likelihood ratio as weight of evidence: A closer look" by Lund and Iyer

Recently, Lund and Iyer (L&I) raised an argument regarding the use of likelihood ratios in court. In our view, their argument is based on a lack of understanding of the paradigm. L&I argue that the decision maker should not accept the expert's likelihood ratio without further consideration. This is agreed by all parties. In normal practice, there is often considerable and proper exploration in court of the basis for any probabilistic statement. We conclude that L&I argue against a practice that does not exist and which no one advocates. Further we conclude that the most informative summary of evidential weight is the likelihood ratio. We state that this is the summary that should be presented to a court in every scientific assessment of evidential weight with supporting information about how it was constructed and on what it was based.

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-016-0429-7) contains supplementary material, which is available to authorized users.

Genetic mapping of 15 human X chromosomal forensic short tandem repeat (STR) loci by means of multi-core parallelization

Typing of X chromosomal short tandem repeat (X STR) markers has become a standard element of human forensic genetic analysis. Joint consideration of many X STR markers at a time increases their discriminatory power but, owing to physical linkage, requires inter-marker recombination rates to be accurately known. We estimated the recombination rates between 15 well established X STR markers using genotype data from 158 families (1041 individuals) and following a previously proposed likelihood-based approach that allows for single-step mutations. To meet the computational requirements of this family-based type of analysis, we modified a previous implementation so as to allow multi-core parallelization on a high-performance computing system. While we obtained recombination rate estimates larger than zero for all but one pair of adjacent markers within the four previously proposed linkage groups, none of the three X STR pairs defining the junctions of these groups yielded a recombination rate estimate of 0.50. Corroborating previous studies, our results therefore argue against a simple model of independent X chromosomal linkage groups. Moreover, the refined recombination fraction estimates obtained in our study will facilitate the appropriate joint consideration of all 15 investigated markers in forensic analysis.

Low-template DNA: A single DNA analysis or two replicates?

This study investigates the following two questions: (1) Should the DNA analyst concentrate the DNA extract into a single amplification or should he/she split it up to do two replicates? (2) Given the electropherogram obtained from a first analysis, is it worthwhile for the DNA analyst to invest in obtaining a second replicate? A decision-theoretic approach addresses these questions by quantitatively expressing the expected net gain (ENG) of each DNA analysis of interest. The results indicate that two replicates generally have a greater ENG than a single DNA analysis for DNA quantities capable of producing two replicates having an average allelic peak height as low as 43rfu. This supports the position that two replicates increase the information content with regard to a single analysis.

D5S2500 is an ambiguously characterized STR: Identification and description of forensic microsatellites in the genomics age

In the process of establishing short tandem repeat (STR) sequence variant nomenclature guidelines in anticipation of expanded forensic multiplexes for massively parallel sequencing (MPS), it was discovered that the STR D5S2500 has multiple positions and genomic characteristics reported. This ambiguity is because the marker named D5S2500 consists of two different microsatellites forming separate components in the capillary electrophoresis multiplexes of Qiagen's HDplex (Hilden, Germany) and AGCU ScienTech's non-CODIS STR 21plex (Wuxi, Jiangsu, China). This study outlines the genomic details used to identify each microsatellite and reveals the D5S2500 marker in HDplex has the correctly assigned STR name, while the D5S2500 marker in the AGCU 21plex, closely positioned a further 1643 nucleotides in the human reference sequence, is an unnamed microsatellite. The fact that the D5S2500 marker has existed as two distinct STR loci undetected for almost ten years, even with reported discordant genotypes for the standard control DNA, underlines the need for careful scrutiny of the genomic properties of forensic STRs, as they become adapted for sequence analysis with MPS systems. We make the recommendation that precise chromosome location data must be reported for any forensic marker under development but not in common use, so that the genomic characteristics of the locus are validated to the same level of accuracy as its allelic variation and forensic performance. To clearly differentiate each microsatellite, we propose the name D5S2800 be used to identify the Chromosome-5 STR in the AGCU 21plex.

Uncertainty in the number of contributors in the proposed new CODIS set

The probability that multiple contributors are detected within a forensic DNA profile improves as more highly polymorphic loci are analysed. The assignment of the correct number of contributors to a profile is important when interpreting the DNA profiles. In this work we investigate the probability of a mixed DNA profile appearing as having originated from a fewer number of contributors for the African American, Asian, Caucasian and Hispanic US populations. We investigate a range of locus configurations from the proposed new CODIS set. These theoretical calculations are based on allele frequencies only and ignore peak heights. We show that the probability of a higher order mixture (five or six contributors) appearing as having originated from one less individual is high. This probability decreases as the number of loci tested increases.

A global analysis of Y-chromosomal haplotype diversity for 23 STR loci

In a worldwide collaborative effort, 19,630 Y-chromosomes were sampled from 129 different populations in 51 countries. These chromosomes were typed for 23 short-tandem repeat (STR) loci (DYS19, DYS389I, DYS389II, DYS390, DYS391, DYS392, DYS393, DYS385ab, DYS437, DYS438, DYS439, DYS448, DYS456, DYS458, DYS635, GATAH4, DYS481, DYS533, DYS549, DYS570, DYS576, and DYS643) and using the PowerPlex Y23 System (PPY23, Promega Corporation, Madison, WI). Locus-specific allelic spectra of these markers were determined and a consistently high level of allelic diversity was observed. A considerable number of null, duplicate and off-ladder alleles were revealed. Standard single-locus and haplotype-based parameters were calculated and compared between subsets of Y-STR markers established for forensic casework. The PPY23 marker set provides substantially stronger discriminatory power than other available kits but at the same time reveals the same general patterns of population structure as other marker sets. A strong correlation was observed between the number of Y-STRs included in a marker set and some of the forensic parameters under study. Interestingly a weak but consistent trend toward smaller genetic distances resulting from larger numbers of markers became apparent.

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 base composition analysis and Sanger sequencing of mitochondrial DNA for four U.S. population groups

A set of 711 samples from four U.S. population groups was analyzed using a novel mass spectrometry based method for mitochondrial DNA (mtDNA) base composition profiling. Comparison of the mass spectrometry results with Sanger sequencing derived data yielded a concordance rate of 99.97%. Length heteroplasmy was identified in 46% of samples and point heteroplasmy was observed in 6.6% of samples in the combined mass spectral and Sanger data set. Using discrimination capacity as a metric, Sanger sequencing of the full control region had the highest discriminatory power, followed by the mass spectrometry base composition method, which was more discriminating than Sanger sequencing of just the hypervariable regions. This trend is in agreement with the number of nucleotides covered by each of the three assays.

Characterising the STR locus D6S1043 and examination of its effect on stutter rates

The forensic analysis of DNA is most often undertaken by the amplification of short tandem repeats (STR) using the polymerase chain reaction (PCR). DNA amplification can result in production of the target allele amplicon and a by-product called stutter. Stutter is the result of the miscopy of the target allele and is typically one repeat smaller. Stutter is traditionally described as a ratio of stutter and allele height; stutter ratio (SR). The challenge to DNA profile interpretation is most serious whenever stutter products are of a similar height to the minor allelic peaks in a mixed DNA profile. An accurate assignment of peaks and the prediction of their height is important when objectively interpreting forensic DNA profiles. The longest uninterrupted stretch (LUS) of tandem repeats within the allele has previously been shown to be a good predictor of stutter ratio. LUS is determined by sequencing a range of observed alleles at a locus. The locus D6S1043 is a relatively new locus to appear in commercial forensic DNA testing kits. To date however, there has been no comprehensive report of sequencing of this locus. In this work, we sequence a sample of D6S1043 alleles to determine LUS values and investigate allele repeat number and LUS as explanatory variables for SR.

An optimized protocol for forensic application of the PreCR™ Repair Mix to multiplex STR amplification of UV-damaged DNA

Damage to the DNA molecule can occur through exposure to environmental conditions such as ultraviolet light, heat and humidity. Forensic samples are particularly prone to such damage due to their prolonged exposure after deposition at crime scenes or mass disasters. Current methods for typing such samples rely heavily on the intact DNA template, and can be adversely affected by damage that is present. Proposed solutions center around increased access to the smaller remaining fragments and/or increased sensitivity. However, all rely on the polymerase chain reaction to copy the starting material; the required polymerase can be impeded by certain types of damage such as dimer-formation after ultraviolet light exposure, resulting in stochastic effects that can complicate interpretation. In vitro repair of such damage offers the ability to generate high quality profiles using traditional methods without changes to the current amplification reagents or conditions. Typically, repair reactions required large quantities of starting material and a separate repair reaction. Forensic samples, however, usually consist of small quantities, and quality control measures necessitate laboratory procedures that minimize sample manipulation. Here, an optimized protocol for forensic application of the PreCR™ Repair Mix to current typing methods is demonstrated for samples damaged by ultraviolet light exposure.

The identification of the Romanovs: Can we (finally) put the controversies to rest?

For much of the 20^th century the fate of the last Imperial family of Russia, the Romanovs, was a mystery after their execution in 1918. In the mid 1970s the mass grave of the Romanov family (minus two of the children) was discovered and officially exhumed after the fall of the Soviet Union. Forensic DNA testing of the remains in the early 1990s was used to identify the family. Despite the overwhelming evidence for establishing the identity of the Romanov family, a small but vocal number of scientists have tried to raise doubt about the DNA testing during the late 1990s and early 2000s. With the discovery of the two missing Romanov children in 2007, there was an opportunity to re-analyze all of the evidence associated with the case which confirmed the initial DNA testing and brought finality to the mystery. This article will discuss the controversies associated with the Romanov identification and reflect upon the importance of the case to the field of forensic DNA typing over the last 20 years.

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.

Inspecting close maternal relatedness: Towards better mtDNA population samples in forensic databases

Reliable data are crucial for all research fields applying mitochondrial DNA (mtDNA) as a genetic marker. Quality control measures have been introduced to ensure the highest standards in sequence data generation, validation and a posteriori inspection. A phylogenetic alignment strategy has been widely accepted as a prerequisite for data comparability and database searches, for forensic applications, for reconstructions of human migrations and for correct interpretation of mtDNA mutations in medical genetics.

There is continuing effort to enhance the number of worldwide population samples in order to contribute to a better understanding of human mtDNA variation. This has often lead to the analysis of convenience samples collected for other purposes, which might not meet the quality requirement of random sampling for mtDNA data sets. Here, we introduce an additional quality control means that deals with one aspect of this limitation: by combining autosomal short tandem repeat (STR) marker with mtDNA information, it helps to avoid the bias introduced by related individuals included in the same (small) sample. By STR analysis of individuals sharing their mitochondrial haplotype, pedigree construction and subsequent software-assisted calculation of likelihood ratios based on the allele frequencies found in the population, closely maternally related individuals can be identified and excluded. We also discuss scenarios that allow related individuals in the same set.

An ideal population sample would be representative for its population: this new approach represents another contribution towards this goal.

Development and characterization of two mini-X chromosomal short tandem repeat multiplexes

This study presents the development and characterization of two X chromosomal short tandem repeat (STR) multiplexes utilizing reduced-size amplicons (less than 200 base pairs) for identity and kinship testing with degraded DNA. Approximately 1360 samples across 4 U.S. population groups were typed for 15 X chromosomal STR markers: DXS6789, DXS7130, DXS9902, GATA31E08, DXS7424, GATA165B12, DXS101, DXS6795, GATA172D05, DXS10147, DXS8378, DXS7132, DXS6803, HPRTB, and DXS7423. A high degree of polymorphism was observed for each marker and both multiplexes were sensitive down to 200pg of pristine DNA. The two proposed multiplexes are suitable for forensic use, and show potential for improved analysis of compromised bone samples.

Population genetic data for 17 STR markers from Lebanon

Seventeen autosomal STRs were analyzed (D2S1338, D3S1358, D5S818, D7S820, D8S1179, D13S317, D16S539, D18S51, D19S433, D21S11, CSF1PO, FGA, TH01, TPOX, vWA, Penta D, and Penta E) in the Lebanese population. A total of 192 unrelated individuals were genotyped for the 15 autosomal STRs in the Promega PowerPlex 16 STR kit. An additional 275 unrelated individuals were genotyped for the Applied Biosystems AmpFlSTR Identifiler and SGM+STR kits. Allele frequencies for the shared CODIS 13 loci among the three STR kits tested were not significantly different among individuals within the Lebanese population. Forensic and population genetic parameters for the 17 loci were calculated. We also compared the allele frequencies from this population with other populations in the same geographic vicinity.

Autosomal SNP typing of forensic samples with the GenPlex™ HID System: results of a collaborative study

The GenPlex™ HID System (Applied Biosystems - AB) offers typing of 48 of the 52 SNPforID SNPs and amelogenin. Previous studies have shown a high reproducibility of the GenPlex™ HID System using 250-500pg DNA of good quality. An international exercise was performed by 14 laboratories (9 in Europe and 5 in the US) in order to test the robustness and reliability of the GenPlex™ HID System on forensic samples. Three samples with partly degraded DNA and 10 samples with low amounts of DNA were analyzed in duplicates using various amounts of DNA. In order to compare the performance of the GenPlex™ HID System with the most commonly used STR kits, 500pg of partly degraded DNA from three samples was typed by the laboratories using one or more STR kits. The median SNP typing success rate was 92.3% with 500pg of partly degraded DNA. Three of the fourteen laboratories counted for more than two thirds of the locus dropouts. The median percentage of discrepant results was 0.2% with 500pg degraded DNA. An increasing percentage of locus dropouts and discrepant results were observed when lower amounts of DNA were used. Different success rates were observed for the various SNPs. The rs763869 SNP was the least successful. With the exception of the MiniFiler™ kit (AB), GenPlex™ HID performed better than five other tested STR kits. When partly degraded DNA was analyzed, GenPlex™ HID showed a very low mean mach probability, while all STR kits except MiniFiler™ had very limited discriminatory power.

Titanic's unknown child: the critical role of the mitochondrial DNA coding region in a re-identification effort

This report describes a re-examination of the remains of a young male child recovered in the Northwest Atlantic following the loss of the Royal Mail Ship Titanic in 1912 and buried as an unknown in Halifax, Nova Scotia shortly thereafter. Following exhumation of the grave in 2001, mitochondrial DNA (mtDNA) hypervariable region 1 sequencing and odontological examination of the extremely limited skeletal remains resulted in the identification of the child as Eino Viljami Panula, a 13-month-old Finnish boy. This paper details recent and more extensive mitochondrial genome analyses that indicate the remains are instead most likely those of an English child, Sidney Leslie Goodwin. The case demonstrates the benefit of targeted mtDNA coding region typing in difficult forensic cases, and highlights the need for entire mtDNA sequence databases appropriate for forensic use.

WITHDRAWN: Development and Characterization of Two Mini-X Chromosomal Short Tandem Repeat Multiplexes

This article has been withdrawn at the request of the author(s) and/or editor. The Publisher apologizes for any inconvenience this may cause. The full Elsevier Policy on Article Withdrawal can be found at http://www.elsevier.com/locate/withdrawalpolicy.

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.

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.

Mystery solved: the identification of the two missing Romanov children using DNA analysis

One of the greatest mysteries for most of the twentieth century was the fate of the Romanov family, the last Russian monarchy. Following the abdication of Tsar Nicholas II, he and his wife, Alexandra, and their five children were eventually exiled to the city of Yekaterinburg. The family, along with four loyal members of their staff, was held captive by members of the Ural Soviet. According to historical reports, in the early morning hours of July 17, 1918 the entire family along with four loyal members of their staff was executed by a firing squad. After a failed attempt to dispose of the remains in an abandoned mine shaft, the bodies were transported to an open field only a few kilometers from the mine shaft. Nine members of the group were buried in one mass grave while two of the children were buried in a separate grave. With the official discovery of the larger mass grave in 1991, and subsequent DNA testing to confirm the identities of the Tsar, the Tsarina, and three of their daughters--doubt persisted that these remains were in fact those of the Romanov family. In the summer of 2007, a group of amateur archeologists discovered a collection of remains from the second grave approximately 70 meters from the larger grave. We report forensic DNA testing on the remains discovered in 2007 using mitochondrial DNA (mtDNA), autosomal STR, and Y-STR testing. Combined with additional DNA testing of material from the 1991 grave, we have virtually irrefutable evidence that the two individuals recovered from the 2007 grave are the two missing children of the Romanov family: the Tsarevich Alexei and one of his sisters.

The phylogeny of the four pan-American MtDNA haplogroups: implications for evolutionary and disease studies

Only a limited number of complete mitochondrial genome sequences belonging to Native American haplogroups were available until recently, which left America as the continent with the least amount of information about sequence variation of entire mitochondrial DNAs. In this study, a comprehensive overview of all available complete mitochondrial DNA (mtDNA) genomes of the four pan-American haplogroups A2, B2, C1, and D1 is provided by revising the information scattered throughout GenBank and the literature, and adding 14 novel mtDNA sequences. The phylogenies of haplogroups A2, B2, C1, and D1 reveal a large number of sub-haplogroups but suggest that the ancestral Beringian population(s) contributed only six (successful) founder haplotypes to these haplogroups. The derived clades are overall starlike with coalescence times ranging from 18,000 to 21,000 years (with one exception) using the conventional calibration. The average of about 19,000 years somewhat contrasts with the corresponding lower age of about 13,500 years that was recently proposed by employing a different calibration and estimation approach. Our estimate indicates a human entry and spread of the pan-American haplogroups into the Americas right after the peak of the Last Glacial Maximum and comfortably agrees with the undisputed ages of the earliest Paleoindians in South America. In addition, the phylogenetic approach also indicates that the pathogenic status proposed for various mtDNA mutations, which actually define branches of Native American haplogroups, was based on insufficient grounds.

Development and expansion of high-quality control region databases to improve forensic mtDNA evidence interpretation

In an effort to increase the quantity, breadth and availability of mtDNA databases suitable for forensic comparisons, we have developed a high-throughput process to generate approximately 5000 control region sequences per year from regional US populations, global populations from which the current US population is derived and global populations currently under-represented in available forensic databases. The system utilizes robotic instrumentation for all laboratory steps from pre-extraction through sequence detection, and a rigorous eight-step, multi-laboratory data review process with entirely electronic data transfer. Over the past 3 years, nearly 10,000 control region sequences have been generated using this approach. These data are being made publicly available and should further address the need for consistent, high-quality mtDNA databases for forensic testing.