A standalone humanitarian DNA identification database system to increase identification of human remains of foreign nationals

The identification of missing persons and human remains is a worldwide problem which has been exacerbated with increased migrations and rampant human trafficking and smuggling cases. DNA typing and DNA databases are primary tools and resources used to help identify human remains and missing persons. The foundation of most, if not all, national DNA database systems, e.g., CODIS, is law enforcement identification. With such database systems, compliance with statutory and operational requirements is necessary to ensure the integrity of the databases. However, because of conditions in their homelands, relatives of missing persons at times may not trust the government and may be reluctant to contact a law enforcement agency, making it difficult to satisfy the law enforcement nexus necessary for entry into a national DNA database. A potential solution to increase the identification of unidentified human remains found within the USA, such as those that may be of foreign nationals, the University of North Texas Center for Human Identification (UNTCHI) has created a Humanitarian DNA Identification DNA Database (HDID) that enables family reference sample DNA profiles from non-US citizens to be compared with the DNA profiles from unidentified human remains within its local database system. This short communication describes the needs, basis, policies, and practices to inform the scientific, investigative, and legal communities and the public so that various entities may become aware and consider submitting family reference sample (FRS) profiles from foreign nationals for the purpose of searching against UNTCHI's HDID. It is our hope that by creating this HDID, another vehicle is available to support identification of human remains within the USA and to bring much needed answers to the family members of missing persons. The HDID will merge high forensic quality and best practices with the broader accessibility for non-US families to voluntarily donate DNA profiles for searching for missing loved ones.

Numt identification and removal with RtN!

Motivation

Assays in mitochondrial genomics rely on accurate read mapping and variant calling. However, there are known and unknown nuclear paralogs that have fundamentally different genetic properties than that of the mitochondrial genome. Such paralogs complicate the interpretation of mitochondrial genome data and confound variant calling.

Results

Remove the Numts! (RtN!) was developed to categorize reads from massively parallel sequencing data not based on the expected properties and sequence identities of paralogous nuclear encoded mitochondrial sequences, but instead using sequence similarity to a large database of publicly available mitochondrial genomes. RtN! removes low-level sequencing noise and mitochondrial paralogs while not impacting variant calling, while competing methods were shown to remove true variants from mitochondrial mixtures.

Availability and implementation

https://github.com/Ahhgust/RtN

Supplementary information

Supplementary data are available at Bioinformatics online.

The lot-to-lot variability in the mitochondrial genome of controls

Current research in the biomedical field has illustrated how cell lines used as reference standards can change over time and, more importantly, can affect research and diagnostic results obtained from these cell lines. With the use of increasingly sensitive and highly resolving technologies (e.g., massively parallel sequencing), forensic scientists must be aware of and account for potential variability in the cell lines used as controls in their validation studies and day-to-day casework. In this study, multiple lot numbers from four commonly-used control cell line DNAs were sequenced with massively parallel sequencing on the Ion S5. The variability among these different lots was evaluated, and the effect on forensic laboratory work discussed.

Utility of the Ion S5™ and MiSeq FGx™ sequencing platforms to characterize challenging human remains

Often in missing persons' and mass disaster cases, the samples remaining for analysis are hard tissues such as bones, teeth, nails, and hair. These remains may have been exposed to harsh environmental conditions, which pose challenges for downstream genotyping. Short tandem repeat analysis (STR) via capillary electrophoresis (CE) is still the gold standard for DNA typing; however, a newer technology known as massively parallel sequencing (MPS) could improve upon our current techniques by typing different and more markers in a single analysis, and consequently improving the power of discrimination. In this study, bone and tooth samples exposed to a variety of DNA insults (cremation, embalming, decomposition, thermal degradation, and fire) were assessed and sequenced using the Precision ID chemistry and a custom AmpliSeq™ STR and iiSNP panel on the Ion S5™ System, and the ForenSeq DNA Signature Prep Kit on the MiSeq FGx™ system, as well as the GlobalFiler™ PCR Amplification Kit on the 3500™ Genetic Analyzer. The results demonstrated that using traditional CE-based genotyping performed as expected, producing a partial or full DNA profile for all samples, and that both sequencing chemistries and platforms were able to recover sufficient STR and SNP information from a majority of the same challenging samples. Run metrics including profile completeness and mean read depth produced good results with each system, considering the degree of damage of some samples. Most sample insults (except decomposed) produced similar numbers of alleles for both MPS systems. Comparable markers produced full concordance between the two platforms.

Evaluation of mitogenome sequence concordance, heteroplasmy detection, and haplogrouping in a worldwide lineage study using the Precision ID mtDNA Whole Genome Panel

The emergence of Massively Parallel Sequencing technologies enabled the analysis of full mitochondrial (mt)DNA sequences from forensically relevant samples that have, so far, only been typed in the control region or its hypervariable segments. In this study, we evaluated the performance of a commercially available multiplex-PCR-based assay, the Precision ID mtDNA Whole Genome Panel (Thermo Fisher Scientific), for the amplification and sequencing of the entire mitochondrial genome (mitogenome) from even degraded forensic specimens. For this purpose, more than 500 samples from 24 different populations were selected to cover the vast majority of established superhaplogroups. These are known to harbor different signature sequence motifs corresponding to their phylogenetic background that could have an effect on primer binding and, thus, could limit a broad application of this molecular genetic tool. The selected samples derived from various forensically relevant tissue sources and were DNA extracted using different methods. We evaluated sequence concordance and heteroplasmy detection and compared the findings to conventional Sanger sequencing as well as an orthogonal MPS platform. We discuss advantages and limitations of this approach with respect to forensic genetic workflow and analytical requirements.

Linkage, recombination, and mutation rate analyses of 19 X-chromosomal STR loci in Chinese Southern Han pedigrees

From Southern Han Chinese samples, we analyzed 19 X-STR markers for linkage, linkage disequilibrium (LD), and mutation rate. The data were collected from two- and three-generation Southern Han Chinese families. These data suggested that both linkage and linkage disequilibrium should be considered while calculating likelihood ratios with X-STR markers in relationship tests. The linkage disequilibrium of these 19 X-STR markers was calculated in our previous research study that was conducted on Southern Han Chinese population. In this study, the recombination fractions between pairs of markers and those obtained from the second-generation Rutgers combined linkage-physical map of the human genome were compared. The observed differences indicated that recombination was not homogeneous along the X chromosome. Therefore, we evaluated the effect on likelihood calculations by referring to haplotype frequencies obtained from allele distributions rather than haplotype counts of Southern Han Chinese population.

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.

Aircraft-Assisted Pilot Suicides in the General Aviation Increased for One-Year Period after 11 September 2001 Attack in the United States

Pilot aircraft-assisted suicides (AAS) are rare, and there is limited understanding of copycat phenomenon among aviators. The aim of this study was to evaluate the possible effect the 11 September 2001, terrorist attacks had on pilot AASs in the U.S. Fatal aviation accidents in the National Transportation Safety Board (NTSB) database were searched using the following search words: “suicide”, “murder-suicide” and “homicide-suicide”. The timeline between 11 September 1996, and 11 September 2004, was analyzed. Only those accidents in which NTSB judged that the cause of the accident was suicide were included in the final analysis. The relative risk (RR) of the pilot AASs in all fatal accidents in the U.S. was calculated in order to compare the one, two, and three-year periods after the September 11 terrorist attacks with five years preceding the event. The RR of a fatal general aviation aircraft accident being due to pilot suicide was 3.68-fold (95% confidence interval 1.04–12.98) during the first year after 11 September 2001, but there was not a statistically significant increase in the later years. This study showed an association, albeit not determinate causal effect, of a very specific series of simultaneous terrorist murder-suicides with subsequent pilot AASs.

Expanding beyond the current core STR loci: An exploration of 73 STR markers with increased diversity for enhanced DNA mixture deconvolution

Current approaches to mixture deconvolution of complex biological samples at times do not have the capability to resolve component contributors in DNA evidence. Additional short tandem repeat (STR) loci were sought that may improve the forensic genetic analysis of mixtures. This study presents exploratory data of a multiplex comprised of 73 highly polymorphic STRs (referred herein as the 73Plex) that were selected because of their high diversity due to sequence variation. These STRs (or a subset of them) may be considered as candidates that may augment current core markers capabilities for DNA mixture deconvolution. Population genetics analyses were performed for each locus using DNA samples from 451 individuals comprising three U.S. populations. Sequence-based heterozygosities ranged from 72% to 98%, where only two loci (D10A97 and D6A7) fell below 80%. Mixture deconvolution capabilities for two-person mixtures were assessed based on complete allele resolution per locus (i.e., four alleles observed) of pairwise mixtures using in silico methods. A subset of 20 highly informative loci (referred herein as the 20Plex) from the 73Plex was compared to the 20 CODIS core loci on all population samples with full DNA profiles for both panels (i.e., no locus dropout; n = 443). Based on proportion of loci displaying four alleles, the 20Plex outperformed the CODIS core loci with increases of 82.6% and 89.3% using length-based and sequence-based alleles, respectively. A combination of 17 STR from the 20Plex and 3 CODIS loci gave the highest capacity for resolving allelic components per locus. These data illustrate the increased value of utilizing sequenced-based alleles of additional STR loci. Furthermore, there are a number of candidate STR loci that could notably augment the current core STR loci and enhance mixture interpretation capabilities.

Assessment of impact of DNA extraction methods on analysis of human remain samples on massively parallel sequencing success

Skeletal remains recovered from missing persons' cases are often exposed to harsh environmental conditions resulting in the DNA being damaged, degraded, and/or the samples containing PCR inhibitors. In this study, the efficacy of common extraction methods was evaluated to remove high levels of PCR inhibitors commonly encountered with human remains, and their downstream compatibility with the two leading sequencing chemistries and platforms for human identification purposes. Blood, hair, and bone samples were spiked with high levels of inhibitors commonly identified in each particular substrate in order to test the efficiency of various DNA extraction methods prior to sequencing. Samples were extracted using three commercial extraction kits (DNA IQ™, DNA Investigator, and PrepFiler® BTA), organic (blood and hair only), and two total demineralization protocols (bone only)). Massively parallel sequencing (MPS) was performed using two different systems: Precision ID chemistry and a custom AmpliSeq™ STR and iiSNP panel on the Ion S5™ System and the ForenSeq DNA Signature Prep Kit on the MiSeq FGx™. The overall results showed that all DNA extraction methods were efficient and are fully compatible with both MPS systems. Key performance indicators such as STR and SNP reportable alleles, read depth, and heterozygote balance were comparable for each extraction method. In samples where CE-based STRs yielded partial profiles (bone), MPS-based STRs generated more complete or full profiles. Moreover, MPS panels contain more STR loci than current CE-based STR kits and also include SNPs, which can further increase the power of discrimination obtained from these samples, making MPS a desirable choice for the forensic analysis of such challenging samples.

Microbial enrichment and gene functional categories revealed on the walls of a spent fuel pool of a nuclear power plant

Microorganisms developing in the liner of the spent fuel pool (SFP) and the fuel transfer channel (FTC) of a Nuclear Power Plant (NPP) can form high radiation resistant biofilms and cause corrosion. Due to difficulties and limitations to obtain large samples from SFP and FTC, cotton swabs were used to collect the biofilm from the wall of these installations. Molecular characterization was performed using massively parallel sequencing to obtain a taxonomic and functional gene classification. Also, samples from the drainage system were evaluated because microorganisms may travel over the 12-meter column of the pool water of the Brazilian Nuclear Power Plant (Angra1), which has been functioning since 1985. Regardless of the treatment of the pool water, our data reveal the unexpected presence of Fungi (Basidiomycota and Ascomycota) as the main contaminators of the SFP and FTC. Ustilaginomycetes (Basidiomycota) was the major class contributor (70%) in the SFP and FTC reflecting the little diversity in these sites; nevertheless, Proteobacteria, Actinobacteria, Firmicutes (Bacilli) were present in small proportions. Mapping total reads against six fungal reference genomes indicate that there is, in fact, a high abundance of fungal sequences in samples collected from SFP and FTC. Analysis of the ribosomal internal transcribed spacer (ITS) 1 and 2 regions and the protein found in the mitochondria of eukaryotic cells, cytochrome b (cytb) grouped our sample fungi in the clade 7 as Ustilago and Pseudozyma. In contrast, in the drainage system, Alphaproteobacteria were present in high abundances (55%). The presence of Sphingopyxis, Mesorhizobium, Erythrobacter, Sphingomonas, Novosphingobium, Sphingobium, Chelativorans, Oceanicaulis, Acidovorax, and Cyanobacteria was observed. Based on genomic annotation data, the assessment of the biological function found a higher proportion of protein-coding sequences related to respiration and protein metabolism in SFP and FTC samples. The knowledge of this biological inventory present in the system may contribute to further studies of potential microorganisms that might be useful for bioremediation of nuclear waste.

The Probabilistic Genotyping Software STRmix: Utility and Evidence for its Validity

Forensic DNA interpretation is transitioning from manual interpretation based usually on binary decision-making toward computer-based systems that model the probability of the profile given different explanations for it, termed probabilistic genotyping (PG). Decision-making by laboratories to implement probability-based interpretation should be based on scientific principles for validity and information that supports its utility, such as criteria to support admissibility. The principles behind STRmix™ are outlined in this study and include standard mathematics and modeling of peak heights and variability in those heights. All PG methods generate a likelihood ratio (LR) and require the formulation of propositions. Principles underpinning formulations of propositions include the identification of reasonably assumed contributors. Substantial data have been produced that support precision, error rate, and reliability of PG, and in particular, STRmix™. A current issue is access to the code and quality processes used while coding. There are substantial data that describe the performance, strengths, and limitations of STRmix™, one of the available PG software.

Challenges in investigation of diabetes-related aviation fatalities-an analysis of 1491 subsequent aviation fatalities in USA during 2011-2016

Diabetes mellitus (DM) could cause pilot incapacitation and result in aviation fatalities. The mechanisms could be directly as a consequence of acute hypoglycemia/subacute diabetic ketoacidosis (DKA) or indirectly as an acute cardiovascular event by contributing to the development of atherosclerosis in coronary or carotid and cerebral arteries. In this study, DM-related fatal flight accidents in the US National Transport Bureau's database between years 2011-2016 were analyzed with special emphasis on postmortem (PM) glucose levels and correlation of toxicological reports with anamnestic information on DM. Additionally, autopsy results on coronary arteries were reviewed. In 43 out of 1491 (~ 3%) fatal accidents pilots had DM. Postmortem glucose or glycated hemoglobin percentage (Hb1Ac) was measured in 12 of the 43 cases; while antidiabetic medication was found in 14 of the cases (only two of the cases had both glucose measurements and medication). With the increasing prevalence of DM, a possibility of pilot incapacitation due to DM or complications of DM should be actively studied, even if no anamnestic information of DM was available. While PM hypoglycemia is difficult to assess, we propose a systematic investigation based on measurement of glucose, Hb1Ac%, and ketone bodies, and documentation of atherosclerotic lesions in major arteries to identify or rule out DM as a cause of pilot incapacitation.