Revisiting the genetic background and phylogenetic structure of five Sino-Tibetan-speaking populations: insights from autosomal InDels.

The Tibetan-Yi Corridor, located on the eastern edge of the Tibetan Plateau, is the main route of the people of the plateau. Human settlement and diffusion along the corridor have played a pivotal role in shaping the genetic architecture of Sino-Tibetan-speaking (STs) populations in China. In this study, five STs groups (Chengdu Tibetan, Chengdu Han, Muli Tibetan, Lugu Lake Mosuo and Xichang Yi) settling in the Tibetan-Yi Corridor were genotyped via AGCU InDel 50 kit on the capillary electrophoresis platform to decrypt the genetic landscape and phylogenetic relationship of STs populations and investigate the forensic characteristics. Allele frequency distributions of all autosomal insertion/deletion polymorphisms (InDels) in studied groups comply with Hardy-Weinberg equilibrium. The combined power of discrimination values are 0.9999999999999999998, 0.9999999999999999995, 0.9999999999999999999, 0.999999999999999993 and 0.99999999999999999994, respectively, and all the combined probability of exclusion values exceed 0.9990. Forensically relevant statistics implied that these InDels could be used for individual identification and as a promising alternative to STR profiling in paternity testing. Typical population comparisons showed strikingly high homogeneity among studied STs people, indicating complicated genetic admixture among populations in the Tibetan-Yi Corridor. The STs groups in the Tibetan-Yi Corridor keep close genetic affinity with geographically or linguistically close populations, and the genetic components of investigated populations arose from a mixture of multiple ancestral gene pools (resulting from the admixture from the ancestral Highland Tibetans and ancestral Lowland indigenous populations).

The yield of postmortem genetic testing in sudden death cases with structural findings at autopsy.

Sudden cardiac death (SCD) is often associated with structural abnormalities of the heart during autopsy. This study sought to compare the diagnostic yield of postmortem genetic testing in (1) cases with structural findings of uncertain significance at autopsy to (2) cases with autopsy findings diagnostic of cardiomyopathy. We evaluated 57 SCD cases with structural findings at cardiac autopsy. Next-generation sequencing using a panel of 77 primary electrical disorder and cardiomyopathy genes was performed. Pathogenic and likely pathogenic variants were classified using American College of Medical Genetics (ACMG) consensus guidelines. In 29 cases (51%) autopsy findings of uncertain significance were identified whereas in 28 cases (49%) a diagnosis of cardiomyopathy was established. We identified a pathogenic or likely pathogenic variant in 10 cases (18%); in 1 (3%) case with non-specific autopsy findings compared with 9 (32%) cases with autopsy findings diagnostic of cardiomyopathy (p = 0.0054). The yield of genetic testing in SCD cases with autopsy findings consistent with cardiomyopathy is comparable with the yield in cardiomyopathy patients that are alive. Genetic testing in cases with findings of uncertain significance offers lower clinical utility than in cardiomyopathy, with lower yields than detected previously. This highlights the need for stringent evaluation of variant pathogenicity.

Assessing the forensic efficiency of the GlobalFiler STR loci among the genetically isolated Chechen subpopulation in Jordan.

Short tandem repeats (STRs) are a widely utilized tool in forensic applications, the latter of which range from human identification and paternity testing to population analysis. The GlobalFiler STR loci, which includes 21 autosomal STRS, were analyzed in the Chechen subpopulation of Jordan. Whole blood samples were withdrawn from 159 Jordanian Chechen individuals, and genomic DNA was extracted from each sample. The GlobalFiler™ kit PCR Amplification Kit amplified and analyzed the STR loci on the 3130xl Genetic Analyzer using GeneMapper ID-X software. The combined match probability for the 21 autosomal STR loci was calculated to be 1.06 × 10-24, a number that is highly discriminatory and informative. The SE33 (0.983) and TPOX (0.806) loci exhibited the highest and lowest powers of discrimination, respectively. Conclusively, the current study indicates that the GlobalFiler loci have a high utility in the Jordanian Chechen population, possibly paving the way for the future establishment of a reference population database in Jordan.

How many individuals share a mitochondrial genome?

Mitochondrial DNA (mtDNA) is useful to assist with identification of the source of a biological sample, or to confirm matrilineal relatedness. Although the autosomal genome is much larger, mtDNA has an advantage for forensic applications of multiple copy number per cell, allowing better recovery of sequence information from degraded samples. In addition, biological samples such as fingernails, old bones, teeth and hair have mtDNA but little or no autosomal DNA. The relatively low mutation rate of the mitochondrial genome (mitogenome) means that there can be large sets of matrilineal-related individuals sharing a common mitogenome. Here we present the mitolina simulation software that we use to describe the distribution of the number of mitogenomes in a population that match a given mitogenome, and investigate its dependence on population size and growth rate, and on a database count of the mitogenome. Further, we report on the distribution of the number of meioses separating pairs of individuals with matching mitogenome. Our results have important implications for assessing the weight of mtDNA profile evidence in forensic science, but mtDNA analysis has many non-human applications, for example in tracking the source of ivory. Our methods and software can also be used for simulations to help validate models of population history in human or non-human populations.

NIST interlaboratory studies involving DNA mixtures (MIX05 and MIX13): Variation observed and lessons learned.

Interlaboratory studies are a type of collaborative exercise in which many laboratories are presented with the same set of data to interpret, and the results they produce are examined to get a "big picture" view of the effectiveness and accuracy of analytical protocols used across participating laboratories. In 2005 and again in 2013, the Applied Genetics Group of the National Institute of Standards and Technology (NIST) conducted interlaboratory studies involving DNA mixture interpretation. In the 2005 NIST MIX05 study, 69 laboratories interpreted data in the form of electropherograms of two-person DNA mixtures representing four different mock sexual assault cases with different contributor ratios. In the 2013 NIST MIX13 study,108 laboratories interpreted electropherogram data for five different case scenarios involving two, three, or four contributors, with some of the contributors potentially related. This paper describes the design of these studies, the variations observed among laboratory results, and lessons learned.

An assessment of the information content of likelihood ratios derived from complex mixtures.

With the increasing sensitivity of DNA typing methodologies, as well as increasing awareness by law enforcement of the perceived capabilities of DNA typing, complex mixtures consisting of DNA from two or more contributors are increasingly being encountered. However, insufficient research has been conducted to characterize the ability to distinguish a true contributor (TC) from a known non-contributor (KNC) in these complex samples, and under what specific conditions. In order to investigate this question, sets of six 15-locus Caucasian genotype profiles were simulated and used to create mixtures containing 2-5 contributors. Likelihood ratios were computed for various situations, including varying numbers of contributors and unknowns in the evidence profile, as well as comparisons of the evidence profile to TCs and KNCs. This work was intended to illustrate the best-case scenario, in which all alleles from the TC were detected in the simulated evidence samples. Therefore the possibility of drop-out was not modeled in this study. The computer program DNAMIX was then used to compute LRs comparing the evidence profile to TCs and KNCs. This resulted in 140,000 LRs for each of the two scenarios. These complex mixture simulations show that, even when all alleles are detected (i.e. no drop-out), TCs can generate LRs less than 1 across a 15-locus profile. However, this outcome was rare, 7 of 140,000 replicates (0.005%), and associated only with mixtures comprising 5 contributors in which the numerator hypothesis includes one or more unknown contributors. For KNCs, LRs were found to be greater than 1 in a small number of replicates (75 of 140,000 replicates, or 0.05%). These replicates were limited to 4 and 5 person mixtures with 1 or more unknowns in the numerator. Only 5 of these 75 replicates (0.004%) yielded an LR greater than 1,000. Thus, overall, these results imply that the weight of evidence that can be derived from complex mixtures containing up to 5 contributors, under a scenario in which no drop-out is required to explain any of the contributors, is remarkably high. This is a useful benchmark result on top of which to layer the effects of additional factors, such as drop-out, peak height, and other variables.

Lab Retriever: a software tool for calculating likelihood ratios incorporating a probability of drop-out for forensic DNA profiles.

BACKGROUND: Technological advances have enabled the analysis of very small amounts of DNA in forensic cases. However, the DNA profiles from such evidence are frequently incomplete and can contain contributions from multiple individuals. The complexity of such samples confounds the assessment of the statistical weight of such evidence. One approach to account for this uncertainty is to use a likelihood ratio framework to compare the probability of the evidence profile under different scenarios. While researchers favor the likelihood ratio framework, few open-source software solutions with a graphical user interface implementing these calculations are available for practicing forensic scientists. RESULTS: To address this need, we developed Lab Retriever, an open-source, freely available program that forensic scientists can use to calculate likelihood ratios for complex DNA profiles. Lab Retriever adds a graphical user interface, written primarily in JavaScript, on top of a C++ implementation of the previously published R code of Balding. We redesigned parts of the original Balding algorithm to improve computational speed. In addition to incorporating a probability of allelic drop-out and other critical parameters, Lab Retriever computes likelihood ratios for hypotheses that can include up to four unknown contributors to a mixed sample. These computations are completed nearly instantaneously on a modern PC or Mac computer. CONCLUSIONS: Lab Retriever provides a practical software solution to forensic scientists who wish to assess the statistical weight of evidence for complex DNA profiles. Executable versions of the program are freely available for Mac OSX and Windows operating systems.

Evaluation of the iPLEX® Sample ID Plus Panel designed for the Sequenom MassARRAY® system. A SNP typing assay developed for human identification and sample tracking based on the SNPforID panel.

Sequenom launched the first commercial SNP typing kit for human identification, named the iPLEX(®) Sample ID Plus Panel. The kit amplifies 47 of the 52 SNPs in the SNPforID panel, amelogenin and two Y-chromosome SNPs in one multiplex PCR. The SNPs were analyzed by single base extension (SBE) and Matrix Assisted Laser Desorption/Ionization-Time of Flight Mass Spectrometry (MALDI-TOF MS). In this study, we evaluated the accuracy and sensitivity of the iPLEX(®) Sample ID Plus Panel by comparing the typing results of the iPLEX(®) Sample ID Plus Panel with those obtained with our ISO 17025 accredited SNPforID assay. The average call rate for duplicate typing of any one SNPs in the panel was 90.0% when the mass spectra were analyzed automatically with the MassARRAY(®) TYPER 4.0 genotyping software in real time. Two reproducible inconsistencies were observed (error rate: 0.05%) at two different SNP loci. In addition, four inconsistencies were observed once. The optimal amount of template DNA in the PCR was ≥10ng. There was a relatively high risk of allele and locus drop-outs when ≤1ng template DNA was used. We developed an R script with a stringent set of "forensic analysis parameters" based on the peak height and the signal to noise data exported from the TYPER 4.0 software. With the forensic analysis parameters, all inconsistencies were eliminated in reactions with ≥10ng DNA. However, the average call rate decreased to 69.9%. The iPLEX(®) Sample ID Plus Panel was tested on 10 degraded samples from forensic case-work. Two samples could not be typed, presumably because the samples contained PCR and SBE inhibitors. The average call rate was generally lower for degraded DNA samples and the number of inconsistencies higher than for pristine DNA. However, none of the inconsistencies were reproduced and the highest match probability for the degraded samples typed with the panel was 1.7E-9 using the stringent forensic analysis parameters. Although the relatively low sensitivity of the iPLEX(®) Sample ID Plus Panel makes it inappropriate for typing of trace samples from crime scenes, the panel may be interesting for relationship testing and for identification of e.g. samples in biobanks because of the low reagent costs, the limited hands-on time of the iPLEX(®) assay and the automatic analysis of the mass spectra.

Genetic analyses of the human eye colours using a novel objective method for eye colour classification.

In this study, we present a new objective method for measuring the eye colour on a continuous scale that allows researchers to associate genetic markers with different shades of eye colour. With the use of the custom designed software Digital Iris Analysis Tool (DIAT), the iris was automatically identified and extracted from high resolution digital images. DIAT was made user friendly with a graphical user interface. The software counted the number of blue and brown pixels in the iris image and calculated a Pixel Index of the Eye (PIE-score) that described the eye colour quantitatively. The PIE-score ranged from -1 to 1 (brown to blue). The software eliminated the need for user based interpretation and qualitative eye colour categories. In 94% (570) of 605 analyzed eye images, the iris region was successfully extracted and a PIE-score was calculated. A very high correlation between the PIE-score and the human perception of eye colour was observed. The correlations between the PIE-scores and the six IrisPlex SNPs (HERC2 rs12913832, OCA2 rs1800407, SLC24A4 rs12896399, TYR rs1393350, SLC45A2 rs16891982 and IRF4 rs12203592) were analyzed in 570 individuals. Significant differences (p<10(-6)) in the PIE-scores of the individuals typed as HERC2 rs12913832 G (PIE=0.99) and rs12913832 GA (PIE=-0.71) or A (PIE=-0.87) were observed. We adjusted for the effect of HERC2 rs12913832 and showed that the quantitative PIE-scores were significantly associated with SNPs with minor effects (OCA2 rs1800407, SLC24A4 rs12896399 and TYR rs1393350) on the eye colour. We evaluated the two published prediction models for eye colour (IrisPlex [1] and Snipper[2]) and compared the predictions with the PIE-scores. We found good concordance with the prediction from individuals typed as HERC2 rs12913832 G. However, both methods had difficulties in categorizing individuals typed as HERC2 rs12913832 GA because of the large variation in eye colour in HERC2 rs12913832 GA individuals. With the use of the DIAT software and the PIE-score, it will be possible to automatically compare the iris colour of large numbers of iris images obtained by different studies and to perform large meta-studies that may reveal loci with small effects on the eye colour.

The interpretation of single source and mixed DNA profiles.

A method for interpreting autosomal mixed DNA profiles based on continuous modelling of peak heights is described. MCMC is applied with a model for allelic and stutter heights to produce a probability for the data given a specified genotype combination. The theory extends to handle any number of contributors and replicates, although practical implementation limits analyses to four contributors. The probability of the peak data given a genotype combination has proven to be a highly intuitive probability that may be assessed subjectively by experienced caseworkers. Whilst caseworkers will not assess the probabilities per se, they can broadly judge genotypes that fit the observed data well, and those that fit relatively less well. These probabilities are used when calculating a subsequent likelihood ratio. The method has been trialled on a number of mixed DNA profiles constructed from known contributors. The results have been assessed against a binary approach and also compared with the subjective judgement of an analyst.

Ancestry analysis reveals a predominant Native American component with moderate European admixture in Bolivians.

We have genotyped 46 Ancestry Informative Markers (AIMs) in two of the most populated areas in Bolivia, namely, La Paz (Andean region; n=105), and Chuquisaca (Sub-Andean region; n=73). Using different analytical tools, we inferred admixture proportions of these two American communities by comparing the genetic profiles with those publicly available from the CEPH (Centre d'Etude du Polymorphisme Humain) panel representing three main continental groups (Africa, Europe, and America). By way of simulations, we first evaluated the minimum sample size needed in order to obtain accurate estimates of ancestry proportions. The results indicated that sample sizes above 30 individuals could be large enough to estimate main continental ancestry proportions using the 46 AIMs panel. With the exception of a few individuals, the results also indicated that Bolivians showed a predominantly Native American ancestry with variable levels of European admixture. The proportions of ancestry were statistically different in La Paz and Chuquisaca: the Native American component was 86% and 77% (Mann-Whitney U-test: un-adjusted P-value=2.1×10(-5)), while the European ancestry was 13% and 21% (Mann-Whitney U-test: un-adjusted P-value=3.6×10(-5)), respectively. The African ancestry in Bolivians captured by the AIMs analyzed in the present study was below 2%. The inferred ancestry of Bolivians fits well with previous studies undertaken on haplotype data, indicating a major proportion of Native American lineages. The genetic differences observed in these two groups suggest that forensic genetic analysis should be better performed based on local databases built in the main Bolivian areas.

Evaluation of next generation mtGenome sequencing using the Ion Torrent Personal Genome Machine (PGM).

Insights into the human mitochondrial phylogeny have been primarily achieved by sequencing full mitochondrial genomes (mtGenomes). In forensic genetics (partial) mtGenome information can be used to assign haplotypes to their phylogenetic backgrounds, which may, in turn, have characteristic geographic distributions that would offer useful information in a forensic case. In addition and perhaps even more relevant in the forensic context, haplogroup-specific patterns of mutations form the basis for quality control of mtDNA sequences. The current method for establishing (partial) mtDNA haplotypes is Sanger-type sequencing (STS), which is laborious, time-consuming, and expensive. With the emergence of Next Generation Sequencing (NGS) technologies, the body of available mtDNA data can potentially be extended much more quickly and cost-efficiently. Customized chemistries, laboratory workflows and data analysis packages could support the community and increase the utility of mtDNA analysis in forensics. We have evaluated the performance of mtGenome sequencing using the Personal Genome Machine (PGM) and compared the resulting haplotypes directly with conventional Sanger-type sequencing. A total of 64mtGenomes (>1 million bases) were established that yielded high concordance with the corresponding STS haplotypes (<0.02% differences). About two-thirds of the differences were observed in or around homopolymeric sequence stretches. In addition, the sequence alignment algorithm employed to align NGS reads played a significant role in the analysis of the data and the resulting mtDNA haplotypes. Further development of alignment software would be desirable to facilitate the application of NGS in mtDNA forensic genetics.

Evaluating forensic DNA profiles using peak heights, allowing for multiple donors, allelic dropout and stutters.

Increases in the sensitivity of DNA profiling technology now allow profiles to be obtained from smaller and more degraded DNA samples than was previously possible. The resulting profiles can be highly informative, but the subjective elements in the interpretation make it problematic to achieve the valid and efficient evaluation of evidential strength required in criminal cases. The problems arise from stochastic phenomena such as "dropout" (absence of an allele in the profile that is present in the underlying DNA) and experimental artefacts such as "stutter" that can generate peaks of ambiguous allelic status. Currently in the UK, evidential strength evaluation uses an approach in which the complex signals in the DNA profiles are interpreted in a semi-manual fashion by trained experts aided by a set of guidelines, but also relying substantially on professional judgment. We introduce a statistical model to calculate likelihood ratios for evaluating DNA evidence arising from multiple known and unknown contributors that allows for such stochastic phenomena by incorporating peak heights. Efficient use of peak heights allows for more crime scene profiles to be reported to courts than is currently possible. The model parameters are estimated from experimental data incorporating multiple sources of variability in the profiling system. We report and analyse experimental results from the SGMPlus system, run at 28 amplification cycles with no enhancements, currently used in the UK. Our methods are readily adapted to other DNA profiling systems provided that the experimental data for the parameter estimation is available.

A study of East Timor variability using the SNPforID 52-plex SNP panel.

A set of 52 autosomal single nucleotide polymorphism (SNP) loci was analyzed in 46 unrelated individuals from the East Timor population using the forensic assay previously described by Sanchez et al. (2006) [J.J. Sanchez, C. Phillips, C. Børsting, K. Balogh, M. Bogus, M. Fondevila, C.D. Harrison, E. Musgrave-Brown, A. Salas, D. Syndercombe Court, PM. Schneider, A. Carracedo, N. Morling, A multiplex assay with 52 single nucleotide polymorphisms for human identification, Electrophoresis 27 (2006) 1713-1724]. Allele frequencies are presented for the 52 SNPs with all loci in Hardy-Weinberg equilibrium for the study population. Comparison with African, European, East Asian and Oceanian populations of the CEPH human genome diversity panel (CEPH-HGDP) revealed significant differences in allele frequency distributions between East Timor and each of the above population groups. Statistical parameters measuring forensic informativeness were also calculated and the values obtained reached comparable levels to those previously described for the other global population groups. This is the first study of variability in these SNPs in an Oceanian population outside of the CEPH-HGDP.

Comparison of Bayesian and maximum-likelihood phylogenetic approaches in two legal cases involving accusations of transmission of HIV.

We investigated two cases of alleged criminal transmission of HIV-1 using Bayesian and maximum-likelihood phylogenetic approaches to determine whether the inference method used influenced the outcome in these cases. In the first case, Bayesian methods were used to reexamine gag and env sequences from an earlier investigation in which the HIV-1 strains infecting one of several contacts could not be linked phylogenetically to that of the accused despite strongly suggestive epidemiological evidence. In the second case, maximum-likelihood and Bayesian inference methods were used to investigate the relatedness of gag and env sequences from HIV-1 strains infecting a man accused of intentionally transmitting the virus to several contacts. Bayesian analysis of HIV-1 strains from the first case confirmed earlier results obtained by maximum-likelihood analysis. A monophyletic cluster linking viruses from the accused and three of his direct and indirect contacts was supported, but a linkage between these viruses and a fourth epidemiologically linked contact could not be demonstrated. In the second case, a strong virological link between the accused and two of his contacts, and the absence of links with four other contacts, was confirmed by both maximum-likelihood and Bayesian inference methods. It is important that phylogenetic programs applied in a legal setting are conservative in their outcome. Although Bayesian methods offer computational tractability for large data sets and complex evolutionary models, this study demonstrates they do not assist when clear linkages between viruses are demonstrated using maximum-likelihood methods.

[Organizational problems of disaster victim identification in mass casualties as exemplified by Tu 154-M and Airbus A310 passenger plane crashes].

Managerial experience is described that was gained during the large-scale work on victim identification following mass casualties in the Tu 154-M and Airbus A310 passenger plane crashes. The authors emphasize the necessity to set up a specialized agency of constant readiness meeting modern requirements for the implementation of a system of measures for personality identification. This agency must incorporate relevant departments of the Ministries of Health, Defense, and Emergency Situations as well as investigative authorities and other organizations.

Locked nucleic acids in PCR primers increase sensitivity and performance.

The incorporation of locked nucleic acids (LNAs) into oligonucleotide primers has been shown to increase template binding strength and specificity for DNA amplification. Real-time PCR and DNA sequencing have been shown to be significantly enhanced by the use of LNAs. Theoretically, increasing primers' binding strength may also increase the sensitivity of conventional PCR, reducing minimum template requirements. We compared LNA-modified PCR primers with their standard DNA counterparts for amplification sensitivity with template amounts as low as 5 pg. Although the results are highly dependent on the design of the LNA primers, large increases in peak height can be achieved from as little as 75 pg, as well as clearer and more complete profiles. Increased amplification success with lower template amounts may also be seen. Additionally, the use of LNAs can enhance multiplexing. Thus, incorporating LNAs into PCR primers can increase amplification success, sensitivity, and performance under a wide range of conditions.

DNA Commission of the International Society for Forensic Genetics (ISFG): recommendations regarding the role of forensic genetics for disaster victim identification (DVI).

The ISFG membership consists of scientists and medical professionals specialized in using genetic testing for kinship analysis and the individualization of biological material. This expertise makes the forensic geneticist a resource of advice to international and national organizations dealing with human identifications and causes many DNA laboratories to get involved in DVI tasks. The present recommendations are meant to educate more forensic geneticists about their potential involvement in mass fatality preparedness and possible DVI efforts, as well as to provide practical guidance for each of the laboratories' individual tasks. The idea to work on DNA-specific recommendations was born after a round table discussion dealing with the 2004 Tsunami disaster in south east Asia during the 21st congress of the International Society for Forensic Genetics on the Azores, Portugal, in September 2005. The ensuing discussion between scientists and pathologists that had been involved in the International Center in Khao Lak, Thailand, revealed the need for the scientific community to be better prepared to answer the local authorities' questions by formulating generally acceptable scientific standards for the most efficient use of DNA-based victim identification methods. These recommendations, as well as the many cited references, are intended to provide guidance on establishing preparedness for the forensic genetics laboratory, on collecting and storing ante-mortem and post-mortem samples suitable for DNA analysis, on DNA extraction and genetic typing strategies, on data management, and on issues related to the biostatistical interpretation and reporting of results.