Improved computations for relationship inference using low-coverage sequencing data.

Pedigree inference, for example determining whether two persons are second cousins or unrelated, can be done by comparing their genotypes at a selection of genetic markers. When the data for one or more of the persons is from low-coverage next generation sequencing (lcNGS), currently available computational methods either ignore genetic linkage or do not take advantage of the probabilistic nature of lcNGS data, relying instead on first estimating the genotype. We provide a method and software (see familias.name/lcNGS) bridging the above gap. Simulations indicate how our results are considerably more accurate compared to some previously available alternatives. Our method, utilizing a version of the Lander-Green algorithm, uses a group of symmetries to speed up calculations. This group may be of further interest in other calculations involving linked loci.

Evaluation of microhaplotypes in forensic kinship analysis from a Swedish population perspective.

The development of massively parallel sequencing (MPS) technology has enabled the discovery of several new types of forensic markers where microhaplotypes are one of these promising novel genetic markers. Microhaplotypes are, commonly, less than 300 nucleotides in length and consist of two or more closely linked single-nucleotide polymorphisms (SNPs). In this study, we have examined a custom-made QIAseq Microhaplotype panel (Qiagen), including 45 different microhaplotype loci. DNA libraries were prepared according to the GeneRead DNAseq Targeted Panels V2 library preparation workflow (Qiagen) and sequenced on a MiSeq FGx instrument (Verogen). We evaluated the performance of the panel based on 75 samples of Swedish origin and haplotype frequencies were established. We performed sensitivity studies and could detect haplotypes at input amounts down to 0.8 ng. We also studied mixture samples with two contributors for which haplotypes, for the minor contributor, were detectable down to the level of 1:100. Furthermore, we executed kinship simulations to evaluate the usefulness of this panel in kinship analysis. The results showed that both paternity and full sibling cases can clearly be solved. When simulating a half sibling versus unrelated case scenario, there were, however, some overlap of the likelihood ratio distributions potentially resulting in inconclusiveness. To conclude, the results of this initial study are promising for further implementation of this microhaplotype assay into the forensic field, although we noticed some primer design issues that could be optimized, which possibly would increase the power of the assay.

The use of FTA cards to acquire DNA profiles from postmortem cases.

Filter papers have been used for many years in different applications of molecular biology and have been proven to be a stable way to store DNA waiting to be analyzed. Sampling of DNA on FTA (Flinders Technology Associates) cards is convenient and cost effective compared to alternative approaches involving DNA extractions and storage of DNA extracts. FTA cards are analyzed at many forensic laboratories, and the way to perform direct genetic profiling on buccal swab cards has developed into an almost industrial process. The possibility to include postmortem (PM) samples into an FTA-based workflow would facilitate and speed up the genetic identification process compared to conventional methods, both on a regular basis and in a mass casualty event. In this study, we investigated if FTA cards may be used to carry tissue DNA from deceased and present a high-quality DNA profile from the individual in order to be useful for the identification process. The study also aimed to investigate if a specific body tissue would be preferable, and if decomposed tissue is suitable at all to put on an FTA card in order to obtain a DNA profile. We have compared the quality of the DNA profiles acquired from postmortem tissue on FTA cards, with the results acquired with conventional methods from reference bone/muscle samples from the same individual. Several types of tissues have been tested from different identification cases and scenarios. We concluded that tissue cells from inner organs are suitable to put on FTA cards, and that the obtained DNA profiles have the potential to serve as PM data for identification purposes. In cases including compromised samples, however, it is recommended to keep the tissue sample as a backup if further DNA has to be extracted.

A SNP panel for identity and kinship testing using massive parallel sequencing.

Within forensic genetics, there is still a need for supplementary DNA marker typing in order to increase the power to solve cases for both identity testing and complex kinship issues. One major disadvantage with current capillary electrophoresis (CE) methods is the limitation in DNA marker multiplex capability. By utilizing massive parallel sequencing (MPS) technology, this capability can, however, be increased. We have designed a customized GeneRead DNASeq SNP panel (Qiagen) of 140 previously published autosomal forensically relevant identity SNPs for analysis using MPS. One single amplification step was followed by library preparation using the GeneRead Library Prep workflow (Qiagen). The sequencing was performed on a MiSeq System (Illumina), and the bioinformatic analyses were done using the software Biomedical Genomics Workbench (CLC Bio, Qiagen). Forty-nine individuals from a Swedish population were genotyped in order to establish genotype frequencies and to evaluate the performance of the assay. The analyses showed to have a balanced coverage among the included loci, and the heterozygous balance showed to have less than 0.5 % outliers. Analyses of dilution series of the 2800M Control DNA gave reproducible results down to 0.2 ng DNA input. In addition, typing of FTA samples and bone samples was performed with promising results. Further studies and optimizations are, however, required for a more detailed evaluation of the performance of degraded and PCR-inhibited forensic samples. In summary, the assay offers a straightforward sample-to-genotype workflow and could be useful to gain information in forensic casework, for both identity testing and in order to solve complex kinship issues.

Mixtures with relatives and linked markers.

Mixture DNA profiles commonly appear in forensic genetics, and a large number of statistical methods and software are available for such cases. However, most of the literature concerns mixtures where the contributors are assumed unrelated and the genetic markers are unlinked. In this paper, we consider mixtures of linked markers and related contributors. If no relationships are involved, linkage can be ignored. While unlinked markers can be treated independently, linkage introduces dependencies. The use of linked markers presents statistical and computational challenges, but may also lead to a considerable increase in power since the number of markers available is much larger if we do not require the markers to be unlinked. In addition, some cases that cannot be solved with an unlimited number of unlinked autosomal markers can be solved with linked markers. We focus on two special cases of linked markers: pairs of linked autosomal markers and X-chromosomal markers. A framework is presented for calculation of likelihood ratios for mixtures with general relationships and with linkage between any number of markers. Finally, we explore the effect of linkage disequilibrium, also called allelic association, on the likelihood ratio.

A general model for likelihood computations of genetic marker data accounting for linkage, linkage disequilibrium, and mutations.

Several applications necessitate an unbiased determination of relatedness, be it in linkage or association studies or in a forensic setting. An appropriate model to compute the joint probability of some genetic data for a set of persons given some hypothesis about the pedigree structure is then required. The increasing number of markers available through high-density SNP microarray typing and NGS technologies intensifies the demand, where using a large number of markers may lead to biased results due to strong dependencies between closely located loci, both within pedigrees (linkage) and in the population (allelic association or linkage disequilibrium (LD)). We present a new general model, based on a Markov chain for inheritance patterns and another Markov chain for founder allele patterns, the latter allowing us to account for LD. We also demonstrate a specific implementation for X chromosomal markers that allows for computation of likelihoods based on hypotheses of alleged relationships and genetic marker data. The algorithm can simultaneously account for linkage, LD, and mutations. We demonstrate its feasibility using simulated examples. The algorithm is implemented in the software FamLinkX, providing a user-friendly GUI for Windows systems (FamLinkX, as well as further usage instructions, is freely available at www.famlink.se ). Our software provides the necessary means to solve cases where no previous implementation exists. In addition, the software has the possibility to perform simulations in order to further study the impact of linkage and LD on computed likelihoods for an arbitrary set of markers.

Population genetics of 29 autosomal STRs and 17 Y-chromosomal STRs in a population sample from Afghanistan.

In this study, allele frequencies for 29 autosomal short tandem repeats (STRs) and haplotype frequencies for 17 Y-chromosomal STRs of an Afghan population have been generated. Samples from 348 men and women originating from Afghanistan were analysed for the autosomal STRs, and the combined match probability was estimated to be 7.5 × 10(-37). One hundred and sixty-nine men were analysed for the Y-chromosomal STRs, which resulted in 132 different haplotypes and a haplotype diversity of 0.995.

Ultrasensitive sequencing of STR markers utilizing unique molecular identifiers and the SiMSen-Seq method.

Massively parallel sequencing (MPS) is increasingly applied in forensic short tandem repeat (STR) analysis. The presence of stutter artefacts and other PCR or sequencing errors in the MPS-STR data partly limits the detection of low DNA amounts, e.g., in complex mixtures. Unique molecular identifiers (UMIs) have been applied in several scientific fields to reduce noise in sequencing. UMIs consist of a stretch of random nucleotides, a unique barcode for each starting DNA molecule, that is incorporated in the DNA template using either ligation or PCR. The barcode is used to generate consensus reads, thus removing errors. The SiMSen-Seq (Simple, multiplexed, PCR-based barcoding of DNA for sensitive mutation detection using sequencing) method relies on PCR-based introduction of UMIs and includes a sophisticated hairpin design to reduce unspecific primer binding as well as PCR protocol adjustments to further optimize the reaction. In this study, SiMSen-Seq is applied to develop a proof-of-concept seven STR multiplex for MPS library preparation and an associated bioinformatics pipeline. Additionally, machine learning (ML) models were evaluated to further improve UMI allele calling. Overall, the seven STR multiplex resulted in complete detection and concordant alleles for 47 single-source samples at 1 ng input DNA as well as for low-template samples at 62.5 pg input DNA. For twelve challenging mixtures with minor contributions of 10 pg to 150 pg and ratios of 1-15% relative to the major donor, 99.2% of the expected alleles were detected by applying the UMIs in combination with an ML filter. The main impact of UMIs was a substantially lowered number of artefacts as well as reduced stutter ratios, which were generally below 5% of the parental allele. In conclusion, UMI-based STR sequencing opens new means for improved analysis of challenging crime scene samples including complex mixtures.

Unearthing who and Y at Harewood Cemetery and inference of George Washington's Y-chromosomal haplotype.

An excavation conducted at Harewood Cemetery to identify the unmarked grave of Samuel Washington resulted in the discovery of burials presumably belonging to George Washington's paternal grandnephews and their mother, Lucy Payne. To confirm their identities this study examined Y-chromosomal, mitochondrial, and autosomal DNA from the burials and a living Washington descendant. The burial's Y-STR profile was compared to FamilyTreeDNA's database, which resulted in a one-step difference from the living descendant and an exact match to another Washington. A more complete Y-STR and Y-SNP profile from the descendant was inferred to be the Washington Y profile. Kinship comparisons performed in relation to the descendant, who is a 4th and 5th degree relative of the putative individuals, resulted in >37,000 overlapping autosomal SNPs and strong statistical support with likelihood ratios exceeding one billion. This study highlights the benefits of a multi-marker approach for kinship prediction and DNA-assisted identification of historical remains.

ABCB1 gene polymorphisms are associated with suicide in forensic autopsies.

BACKGROUND: Polymorphisms in ABCB1 have the ability to affect both the function and the expression of the transporter protein P-glycoprotein and may lead to an altered response for many drugs including some antidepressants and antipsychotics. OBJECTIVE: The aim of this study was to examine the impact of the ABCB1 polymorphisms 1199G>A, 1236C>T, 2677G>T/A, and 3435C>T in deaths by suicide. PATIENTS AND METHODS: A total of 998 consecutive Swedish forensic autopsies performed in 2008 in individuals 18 years of age or older, where femoral blood was available and a toxicological screening had been performed, were investigated. Genotypes were assessed with pyrosequencing and information on the cause and manner of each death was obtained from the forensic pathology and toxicology databases. RESULTS: There was a significantly higher frequency of the T allele at positions 1236, 2677, and 3435 among the suicide cases compared with the nonsuicide cases. CONCLUSION: Our result from forensic cases suggests that ABCB1 polymorphisms are associated with an increased risk for completed suicides. The biological mechanisms involved and the clinical implications for these findings are largely unknown and need to be examined further.

Applying Unique Molecular Indices with an Extensive All-in-One Forensic SNP Panel for Improved Genotype Accuracy and Sensitivity.

One of the major challenges in forensic genetics is being able to detect very small amounts of DNA. Massively parallel sequencing (MPS) enables sensitive detection; however, genotype errors may exist and could interfere with the interpretation. Common errors in MPS-based analysis are often induced during PCR or sequencing. Unique molecular indices (UMIs) are short random nucleotide sequences ligated to each template molecule prior to amplification. Applying UMIs can improve the limit of detection by enabling accurate counting of initial template molecules and removal of erroneous data. In this study, we applied the FORCE panel, which includes ~5500 SNPs, with a QIAseq Targeted DNA Custom Panel (Qiagen), including UMIs. Our main objective was to investigate whether UMIs can enhance the sensitivity and accuracy of forensic genotyping and to evaluate the overall assay performance. We analyzed the data both with and without the UMI information, and the results showed that both genotype accuracy and sensitivity were improved when applying UMIs. The results showed very high genotype accuracies (>99%) for both reference DNA and challenging samples, down to 125 pg. To conclude, we show successful assay performance for several forensic applications and improvements in forensic genotyping when applying UMIs.

Extended population genetic analysis of 12 X-STRs - Exemplified using a Norwegian population sample.

The use of X-chromosomal markers to resolve questions of relatedness has experienced a significant increase during the last years in forensic genetics. Perhaps primarily due to the emergence of commercial kits, but equally important due to an increased awareness of the utility of those markers. The X-chromosomal inheritance pattern entails that some cases, for instance paternal half-sisters, can potentially be resolved using a few X-chromosomal markers alone. For the statistical assessment in kinship cases it is of importance to have relevant population frequency data. In the present study 631 unrelated males from a Norwegian population sample are analyzed. The resulting haplotypes are compared to previously studied population samples and a deeper analysis of the linkage disequilibrium (LD) structure is conducted. We demonstrate that the power to detect LD will be low when few males, say below 300, are analyzed. We use entropy to describe the degree of LD between multiallelic loci and describe how this measure varies between different studied populations. Large population frequency databases have been recommended when using X-chromosomal markers, and we show that by combining reference databases from genetically similar populations, more precise haplotype frequency estimates can be obtained for rare haplotypes which improves the statistical assessment of the weight of evidence. In addition, we promote the use of simulations to assess the utility of STR markers in contrast to standard forensic parameters. Specifically we perform extensive simulations on cases where X-chromosomal markers are important and illustrate how the results can be used to infer the information gained from these markers.

Improved DNA Extraction and Illumina Sequencing of DNA Recovered from Aged Rootless Hair Shafts Found in Relics Associated with the Romanov Family.

This study describes an optimized DNA extraction protocol targeting ultrashort DNA molecules from single rootless hairs. It was applied to the oldest samples available to us: locks of hairs that were found in relics associated with the Romanov family. Published mitochondrial DNA genome sequences of Tsar Nicholas II and his wife, Tsarina Alexandra, made these samples ideal to assess this DNA extraction protocol and evaluate the types of genetic information that can be recovered by sequencing ultrashort fragments. Using this method, the mtGenome of the Tsarina's lineage was identified in hairs that were concealed in a pendant made by Karl Fabergé for Alexandra Feodorovna Romanov. In addition, to determine if the lock originated from more than one individual, two hairs from the locket were extracted independently and converted into Illumina libraries for shotgun sequencing on a NextSeq 500 platform. From these data, autosomal SNPs were analyzed to assess relatedness. The results indicated that the two hairs came from a single individual. Genetic testing of hairs that were found in the second artifact, a framed photograph of Louise of Hesse-Kassel, Queen of Denmark and maternal grandmother of Tsar Nicholas II, revealed that the hair belonged to a woman who shared Tsar Nicholas' maternal lineage, including the well-known point heteroplasmy at position 16169.

Investigative genetic genealogy: Current methods, knowledge and practice.

Investigative genetic genealogy (IGG) has emerged as a new, rapidly growing field of forensic science. We describe the process whereby dense SNP data, commonly comprising more than half a million markers, are employed to infer distant relationships. By distant we refer to degrees of relatedness exceeding that of first cousins. We review how methods of relationship matching and SNP analysis on an enlarged scale are used in a forensic setting to identify a suspect in a criminal investigation or a missing person. There is currently a strong need in forensic genetics not only to understand the underlying models to infer relatedness but also to fully explore the DNA technologies and data used in IGG. This review brings together many of the topics and examines their effectiveness and operational limits, while suggesting future directions for their forensic validation. We further investigated the methods used by the major direct-to-consumer (DTC) genetic ancestry testing companies as well as submitting a questionnaire where providers of forensic genetic genealogy summarized their operation/services. Although most of the DTC market, and genetic genealogy in general, has undisclosed, proprietary algorithms we review the current knowledge where information has been discussed and published more openly.

Getting the conclusive lead with investigative genetic genealogy - A successful case study of a 16 year old double murder in Sweden.

On the morning of October 19, 2004, an eight-year-old boy and a 56-year-old woman were stabbed to death on an open street in the city of Linköping, Sweden. The perpetrator left his DNA at the crime scene, and after 15 years of various investigation efforts, including more than 9000 interrogations and mass DNA screening of more than 6000 men, there were still no clues about the identity of the unknown murderer. The successful application of investigative genetic genealogy (IGG) in the US raised the interest for this tool within the Swedish Police Authority. After legal consultations it was decided that IGG could be applied in this double murder case as a pilot case study. From extensive DNA analysis, including whole-genome sequencing and genotype imputation, DNA data sets were established and searched within both GEDmatch and FamilyTree DNA genealogy databases. A number of fairly distant relatives were found from which family trees were created. The genealogy work resulted in two candidates, two brothers, one of whom matched the crime scene samples by routine STR profiling. The suspect confessed the murders at the initial police hearing and was later convicted of the murders. In this paper we describe the successful application of an emerging technology. We disclose details of the DNA analyses which, due to the poor quality and low quantity of the DNA, required reiterative sequencing and genotype imputation efforts. The successful application of IGG in this double murder case exemplifies its applicability not only in the US but also in Europe. The pressure is now high on the involved authorities to establish IGG as a tool for cold case criminal investigations and for missing person identifications. There is, however, a continuous need to accommodate legal, social and ethical aspects as well.

The FORCE Panel: An All-in-One SNP Marker Set for Confirming Investigative Genetic Genealogy Leads and for General Forensic Applications.

The FORensic Capture Enrichment (FORCE) panel is an all-in-one SNP panel for forensic applications. This panel of 5422 markers encompasses common, forensically relevant SNPs (identity, ancestry, phenotype, X- and Y-chromosomal SNPs), a novel set of 3931 autosomal SNPs for extended kinship analysis, and no clinically relevant/disease markers. The FORCE panel was developed as a custom hybridization capture assay utilizing ~20,000 baits to target the selected SNPs. Five non-probative, previously identified World War II (WWII) cases were used to assess the kinship panel. Each case included one bone sample and associated family reference DNA samples. Additionally, seven reference quality samples, two 200-year-old bone samples, and four control DNAs were processed for kit performance and concordance assessments. SNP recovery after capture resulted in a mean of ~99% SNPs exceeding 10X coverage for reference and control samples, and 44.4% SNPs for bone samples. The WWII case results showed that the FORCE panel could predict first to fifth degree relationships with strong statistical support (likelihood ratios over 10,000 and posterior probabilities over 99.99%). To conclude, SNPs will be important for further advances in forensic DNA analysis. The FORCE panel shows promising results and demonstrates the utility of a 5000 SNP panel for forensic applications.

Subdividing Y-chromosome haplogroup R1a1 reveals Norse Viking dispersal lineages in Britain.

The influence of Viking-Age migrants to the British Isles is obvious in archaeological and place-names evidence, but their demographic impact has been unclear. Autosomal genetic analyses support Norse Viking contributions to parts of Britain, but show no signal corresponding to the Danelaw, the region under Scandinavian administrative control from the ninth to eleventh centuries. Y-chromosome haplogroup R1a1 has been considered as a possible marker for Viking migrations because of its high frequency in peninsular Scandinavia (Norway and Sweden). Here we select ten Y-SNPs to discriminate informatively among hg R1a1 sub-haplogroups in Europe, analyse these in 619 hg R1a1 Y chromosomes including 163 from the British Isles, and also type 23 short-tandem repeats (Y-STRs) to assess internal diversity. We find three specifically Western-European sub-haplogroups, two of which predominate in Norway and Sweden, and are also found in Britain; star-like features in the STR networks of these lineages indicate histories of expansion. We ask whether geographical distributions of hg R1a1 overall, and of the two sub-lineages in particular, correlate with regions of Scandinavian influence within Britain. Neither shows any frequency difference between regions that have higher (≥10%) or lower autosomal contributions from Norway and Sweden, but both are significantly overrepresented in the region corresponding to the Danelaw. These differences between autosomal and Y-chromosomal histories suggest either male-specific contribution, or the influence of patrilocality. Comparison of modern DNA with recently available ancient DNA data supports the interpretation that two sub-lineages of hg R1a1 spread with the Vikings from peninsular Scandinavia.

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.

Whole-genome sequencing of human remains to enable genealogy DNA database searches - A case report.

Recently a number of high profile crime cases (e.g. the "Golden State Killer") have successfully been solved or given new leads with the use of genome wide DNA data in combination with pairwise matching from individuals present in genealogy DNA databases. Such databases will primarily involve distant relatives which in turn require a large amount of genetic information, in the range of several hundred thousand to millions of SNPs, to be genotyped. While it nowadays is fairly straightforward to obtain such as data from high quality and high quantity DNA, it is still a challenge for degraded DNA of low quantity such in the case of forensic samples. Here we present a successful effort in obtaining genome-wide genotype data from human remains. The goal was to get investigative leads in order to identify the remains of an unknown male ("the Ekeby man") that was found murdered in the south of Sweden in 2003. Whole-genome sequencing was performed on DNA originating from a bone sample. Three replicates of libraries were prepared using ThruPLEX DNA-seq Kit (Takara) which were sequenced on a HiSeq X instrument (Illumina). A mean coverage of 30X was obtained when the sequencing reads were mapped to a human reference genome. Following further bioinformatic processing, allele calling, quality checks and filtering to match the genealogy DNA database SNPs, genotypes for approximately one million SNPs were established. The resulting SNP genotypes were then used to search for relatives in the genealogy DNA database GEDmatch (www.gedmatch.com). A candidate list of relatives was obtained which was further processed using traditional genealogy methods in order to get leads about the identity of the unknown. In summary, this report shows how whole-genome sequencing successfully can be applied on forensic samples to create the SNP genotypes required for searches in genealogy DNA databases for the purpose of generating leads to identify missing or unknown persons, including perpetrators and victims.