Evaluation of a custom QIAseq targeted DNA panel with 164 ancestry informative markers sequenced with the Illumina MiSeq.

Introduction of new methods requires meticulous evaluation before they can be applied to forensic genetic case work. Here, a custom QIAseq Targeted DNA panel with 164 ancestry informative markers was assessed using the MiSeq sequencing platform. Concordance, sensitivity, and the capability for analysis of mixtures were tested. The assay gave reproducible and nearly concordant results with an input of 10 and 2 ng DNA. Lower DNA input led to an increase in both locus and allele drop-outs, and a higher variation in heterozygote balance. Locus or allele drop-outs in the samples with less than 2 ng DNA input were not necessarily associated with the overall performance of a locus. Thus, the QIAseq assay will be difficult to implement in a forensic genetic setting where the sample material is often scarce and of poor quality. With equal or near equal mixture ratios, the mixture DNA profiles were easily identified by an increased number of imbalanced heterozygotes. For more skewed mixture ratios, the mixture DNA profiles were identified by an increased noise level. Lastly, individuals from Great Britain and the Middle East were investigated. The Middle Eastern individuals showed a greater affinity with South European populations compared to North European populations.

Assessment of the effectiveness of the EUROFORGEN NAME and Precision ID Ancestry panel markers for ancestry investigations.

The EUROFORGEN NAME panel is a regional ancestry panel designed to differentiate individuals from the Middle East, North Africa, and Europe. The first version of the panel was developed for the MassARRAY system and included 111 SNPs. Here, a custom AmpliSeq EUROFORGEN NAME panel with 102 of the original 111 loci was used to sequence 1098 individuals from 14 populations from Europe, the Middle East, North Africa, North-East Africa, and South-Central Asia. These samples were also sequenced with a global ancestry panel, the Precision ID Ancestry Panel. The GenoGeographer software was used to assign the AIM profiles to reference populations and calculate the weight of the evidence as likelihood ratios. The combination of the EUROFORGEN NAME and Precision ID Ancestry panels led to fewer ambiguous assignments, especially for individuals from the Middle East and South-Central Asia. The likelihood ratios showed that North African individuals could be separated from European and Middle Eastern individuals using the Precision ID Ancestry Panel. The separation improved with the addition of the EUROFORGEN NAME panel. The analyses also showed that the separation of Middle Eastern populations from European and South-Central Asian populations was challenging even when both panels were applied.

Body fluid identification and assignment to donors using a targeted mRNA massively parallel sequencing approach - results of a second EUROFORGEN / EDNAP collaborative exercise.

In a previous EUROFORGEN/EDNAP collaborative exercise, we tested two assays for targeted mRNA massively parallel sequencing for the identification of body fluids/tissues, optimized for the Illumina MiSeq/FGx and the Ion Torrent PGM/S5 platforms, respectively. The task of the second EUROFORGEN/EDNAP collaborative exercise was to analyze dried body fluid stains with two different multiplexes, the former Illumina 33plex mRNA panel for body fluid/tissue identification and a 35plex cSNP panel for assignment of body fluids/tissues to donors that was introduced in a proof-of-concept study recently. The coding region SNPs (cSNPs) are located within the body fluid specific mRNA transcripts and represent a direct link between the body fluid and the donor. We predicted the origin of the stains using a partial least squares discriminant analysis (PLS-DA) model, where most of the single source samples were correctly predicted. The mixed body fluid stains showed poorer results, however, at least one component was predicted correctly in most stains. The cSNP data demonstrated that coding region SNPs can give valuable information on linking body fluids/tissues with donors in mixed body fluid stains. However, due to the unfavorable performance of some cSNPs, the interpretation remains challenging. As a consequence, additional markers are needed to increase the discrimination power in each body fluid/tissue category.

Analysis of Uyghur and Kazakh populations using the Precision ID Ancestry Panel.

Autosomal ancestry informative markers (AIMs) are important markers for inferring ancestry of humans. In the present study, we typed 105 Uyghurs and 94 Kazakhs with the Precision ID Ancestry Panel that amplifies 165 autosomal AIMs. No statistically significant deviation from Hardy-Weinberg equilibrium and no linkage disequilibrium between loci was observed after Bonferroni correction. STRUCTURE and PCA analyses showed that Uyghurs and Kazakhs appeared as admixed individuals of primarily European and East Asian ancestry and were clearly differentiated from Europeans, Middle Easterners, South/Central Asians, and East Asians. However, it was not possible to differentiate the two populations from each other and they were also difficult to differentiate from Greenlanders, a population with European/Inuit admixture. GenoGeographer was used to evaluate the weight of the evidence. Initially, the results showed that the majority of AIM profiles from Uyghur and Kazakh individuals were not represented by any of the 36 reference populations of the GenoGeographer database. Consequently, it was not reasonable to infer the ancestry of these individuals. A randomly selected subset of the studied populations (75 Uyghur and 75 Kazakh individuals) was used to construct two new reference populations for GenoGeographer, and ancestry prediction was performed on the remaining test individuals. A total of 42 out of 49 test individuals were represented by at least one population after the introduction of Uyghur and Kazakh reference populations. Likelihood ratios ≥106 were obtained when the alternative hypothesis was that the individual belonged to the South/Central Asian, East Asian, Middle Eastern, European, or the admixed Greenlandic population.

Development and validation of the EUROFORGEN NAME (North African and Middle Eastern) ancestry panel.

Inference of biogeographic origin is an important factor in clinical, population and forensic genetics. The information provided by AIMs (Ancestry Informative Markers) can allow the differentiation of major continental population groups, and several AIM panels have been developed for this purpose. However, from these major population groups, Eurasia covers a wide area between two continents that is difficult to differentiate genetically. These populations display a gradual genetic cline from West Europe to South Asia in terms of allele frequency distribution. Although differences have been reported between Europe and South Asia, Middle East populations continue to be a target of further investigations due to the lack of genetic variability, therefore hampering their genetic differentiation from neighboring populations. In the present study, a custom-built ancestry panel was developed to analyze North African and Middle Eastern populations, designated the 'NAME' panel. The NAME panel contains 111 SNPs that have patterns of allele frequency differentiation that can distinguish individuals originating in North Africa and the Middle East when combined with a previous set of 126 Global AIM-SNPs.

The Danish STR sequence database: duplicate typing of 363 Danes with the ForenSeq™ DNA Signature Prep Kit.

Some STR loci have internal sequence variations, which are not revealed by the standard STR typing methods used in forensic genetics (PCR and fragment length analysis by capillary electrophoresis (CE)). Typing of STRs with next-generation sequencing (NGS) uncovers the sequence variation in the repeat region and in the flanking regions. In this study, 363 Danish individuals were typed for 56 STRs (26 autosomal STRs, 24 Y-STRs, and 6 X-STRs) using the ForenSeq™ DNA Signature Prep Kit to establish a Danish STR sequence database. Increased allelic diversity was observed in 34 STRs by the PCR-NGS assay. The largest increases were found in DYS389II and D12S391, where the numbers of sequenced alleles were around four times larger than the numbers of alleles determined by repeat length alone. Thirteen SNPs and one InDel were identified in the flanking regions of 12 STRs. Furthermore, 36 single positions and five longer stretches in the STR flanking regions were found to have dubious genotyping quality. The combined match probability of the 26 autosomal STRs was 10,000 times larger using the PCR-NGS assay than by using PCR-CE. The typical paternity indices for trios and duos were 500 and 100 times larger, respectively, than those obtained with PCR-CE. The assay also amplified 94 SNPs selected for human identification. Eleven of these loci were not in Hardy-Weinberg equilibrium in the Danish population, most likely because the minimum threshold for allele calling (30 reads) in the ForenSeq™ Universal Analysis Software was too low and frequent allele dropouts were not detected.

Body fluid identification using a targeted mRNA massively parallel sequencing approach - results of a EUROFORGEN/EDNAP collaborative exercise.

In a previous study we presented an assay for targeted mRNA sequencing for the identification of human body fluids, optimised for the Illumina MiSeq/FGx MPS platform. This assay, together with an additional in-house designed assay for the Ion Torrent PGM/S5 platform, was the basis for a collaborative exercise within 17 EUROFORGEN and EDNAP laboratories, in order to test the efficacy of targeted mRNA sequencing to identify body fluids. The task was to analyse the supplied dried body fluid stains and, optionally, participants' own bona fide or mock casework samples of human origin, according to specified protocols. The provided primer pools for the Illumina MiSeq/FGx and the Ion Torrent PGM/S5 platforms included 33 and 29 body fluid specific targets, respectively, to identify blood, saliva, semen, vaginal secretion, menstrual blood and skin. The results demonstrated moderate to high count values in the body fluid or tissue of interest with little to no counts in non-target body fluids. There was some inter-laboratory variability in read counts, but overall the results of the laboratories were comparable in that highly expressed markers showed high read counts and less expressed markers showed lower counts. We performed a partial least squares (PLS) analysis on the data, where blood, menstrual blood, saliva and semen markers and samples clustered well. The results of this collaborative mRNA massively parallel sequencing (MPS) exercise support targeted mRNA sequencing as a reliable body fluid identification method that could be added to the repertoire of forensic MPS panels.

Forensic SNP genotyping with SNaPshot: Technical considerations for the development and optimization of multiplexed SNP assays.

This review explores the key factors that influence the optimization, routine use, and profile interpretation of the SNaPshot single-base extension (SBE) system applied to forensic single-nucleotide polymorphism (SNP) genotyping. Despite being a mainly complimentary DNA genotyping technique to routine STR profiling, use of SNaPshot is an important part of the development of SNP sets for a wide range of forensic applications with these markers, from genotyping highly degraded DNA with very short amplicons to the introduction of SNPs to ascertain the ancestry and physical characteristics of an unidentified contact trace donor. However, this technology, as resourceful as it is, displays several features that depart from the usual STR genotyping far enough to demand a certain degree of expertise from the forensic analyst before tackling the complex casework on which SNaPshot application provides an advantage. In order to provide the basis for developing such expertise, we cover in this paper the most challenging aspects of the SNaPshot technology, focusing on the steps taken to design primer sets, optimize the PCR and single-base extension chemistries, and the important features of the peak patterns observed in typical forensic SNP profiles using SNaPshot. With that purpose in mind, we provide guidelines and troubleshooting for multiplex-SNaPshot-oriented primer design and the resulting capillary electrophoresis (CE) profile interpretation (covering the most commonly observed artifacts and expected departures from the ideal conditions).

A collaborative EDNAP exercise on SNaPshot™-based mtDNA control region typing.

A collaborative European DNA Profiling (EDNAP) Group exercise was undertaken to assess the performance of an earlier described SNaPshot™-based screening assay (denoted mini-mtSNaPshot) (Weiler et al., 2016) [1] that targets 18 single nucleotide polymorphism (SNP) positions in the mitochondrial (mt) DNA control region and allows for discrimination of major European mtDNA haplogroups. Besides the organising laboratory, 14 forensic genetics laboratories were involved in the analysis of 13 samples, which were centrally prepared and thoroughly tested prior to shipment. The samples had a variable complexity and comprised straightforward single-source samples, samples with dropout or altered peak sizing, a point heteroplasmy and two-component mixtures resulting in one to five bi-allelic calls. The overall success rate in obtaining useful results was high (97.6%) given that some of the participating laboratories had no previous experience with the typing technology and/or mtDNA analysis. The majority of the participants proceeded to haplotype inference to assess the feasibility of assigning a haplogroup and checking phylogenetic consistency when only 18 SNPs are typed. To mimic casework procedures, the participants compared the SNP typing data of all 13 samples to a set of eight mtDNA reference profiles that were described according to standard nomenclature (Parson et al., 2014) [2], and indicated whether these references matched each sample or not. Incorrect scorings were obtained for 2% of the comparisons and derived from a subset of the participants, indicating a need for training and guidelines regarding mini-mtSNaPshot data interpretation.

Thirty autosomal insertion-deletion polymorphisms analyzed using the Investigator® DIPplex Kit in populations from Iraq, Lithuania, Slovenia, and Turkey.

Thirty autosomal insertion-deletion (InDel) polymorphisms were analyzed in four populations from Iraq, Lithuania, Slovenia, and Turkey using the commercial kit Investigator® DIPplex. Genotyping issues were encountered for five of the 30 InDels. They were most probably caused by polymorphisms located in the primer binding sites. Population and forensic parameters were calculated. No significant deviations from Hardy-Weinberg equilibrium or significant linkage disequilibrium were detected. The observed heterozygosities ranged from 33% to 61% depending on the marker and the population. The combined probability of exclusion for the 30 markers was 99.7% in all four populations and the matching probabilities were 1 in 3-4×1012 individuals. The multidimensional scaling plot drawn from FST distances showed a good concordance between the relative position of the 15 populations included in the plot and their geographic locations.

DNA Commission of the International Society for Forensic Genetics: Recommendations on the validation of software programs performing biostatistical calculations for forensic genetics applications.

The use of biostatistical software programs to assist in data interpretation and calculate likelihood ratios is essential to forensic geneticists and part of the daily case work flow for both kinship and DNA identification laboratories. Previous recommendations issued by the DNA Commission of the International Society for Forensic Genetics (ISFG) covered the application of bio-statistical evaluations for STR typing results in identification and kinship cases, and this is now being expanded to provide best practices regarding validation and verification of the software required for these calculations. With larger multiplexes, more complex mixtures, and increasing requests for extended family testing, laboratories are relying more than ever on specific software solutions and sufficient validation, training and extensive documentation are of upmost importance. Here, we present recommendations for the minimum requirements to validate bio-statistical software to be used in forensic genetics. We distinguish between developmental validation and the responsibilities of the software developer or provider, and the internal validation studies to be performed by the end user. Recommendations for the software provider address, for example, the documentation of the underlying models used by the software, validation data expectations, version control, implementation and training support, as well as continuity and user notifications. For the internal validations the recommendations include: creating a validation plan, requirements for the range of samples to be tested, Standard Operating Procedure development, and internal laboratory training and education. To ensure that all laboratories have access to a wide range of samples for validation and training purposes the ISFG DNA commission encourages collaborative studies and public repositories of STR typing results.

Inter-laboratory evaluation of the EUROFORGEN Global ancestry-informative SNP panel by massively parallel sequencing using the Ion PGM™.

The EUROFORGEN Global ancestry-informative SNP (AIM-SNPs) panel is a forensic multiplex of 128 markers designed to differentiate an individual's ancestry from amongst the five continental population groups of Africa, Europe, East Asia, Native America, and Oceania. A custom multiplex of AmpliSeq™ PCR primers was designed for the Global AIM-SNPs to perform massively parallel sequencing using the Ion PGM™ system. This study assessed individual SNP genotyping precision using the Ion PGM™, the forensic sensitivity of the multiplex using dilution series, degraded DNA plus simple mixtures, and the ancestry differentiation power of the final panel design, which required substitution of three original ancestry-informative SNPs with alternatives. Fourteen populations that had not been previously analyzed were genotyped using the custom multiplex and these studies allowed assessment of genotyping performance by comparison of data across five laboratories. Results indicate a low level of genotyping error can still occur from sequence misalignment caused by homopolymeric tracts close to the target SNP, despite careful scrutiny of candidate SNPs at the design stage. Such sequence misalignment required the exclusion of component SNP rs2080161 from the Global AIM-SNPs panel. However, the overall genotyping precision and sensitivity of this custom multiplex indicates the Ion PGM™ assay for the Global AIM-SNPs is highly suitable for forensic ancestry analysis with massively parallel sequencing.

NGMSElect™ and Investigator(®) Argus X-12 analysis in population samples from Albania, Iraq, Lithuania, Slovenia, and Turkey.

The analysis of STRs is the main tool when studying genetic diversity in populations or when addressing individual identification in forensic casework. Population data are needed to establish reference databases that can be used in the forensic context. To that end, this work investigated five population samples from Albania, Iraq, Lithuania, Slovenia, and Turkey. Individuals were typed for 16 autosomal STRs and 12 X-chromosomal STRs using the NGMSElect™ and Investigator(®) Argus X-12 kits, respectively. The aim of the study was to characterize the diversity of both STR kits in these population samples and to expand our forensic database. The results showed that all markers were polymorphic in the five populations studied. No haplotype was shared between the males analysed for X-STRs. No statistically significant deviations from Hardy-Weinberg equilibrium were observed for any of the genetic markers included in both the kits. Pairwise LD was only detected in X-STRs between markers located in the same linkage group. Power of discrimination values for males and females and the probability of exclusion in duos and trios were high for the populations in this study.

Next-generation sequencing of 100 candidate genes in young victims of suspected sudden cardiac death with structural abnormalities of the heart.

BACKGROUND: In sudden, unexpected, non-traumatic death in young individuals, structural abnormalities of the heart are frequently identified at autopsy. However, the findings may be unspecific and cause of death may remain unclear. A significant proportion of these cases are most likely caused by inherited cardiac diseases, and the cases are categorized as sudden cardiac death (SCD). The purpose of this study was to explore the added diagnostic value of genetic testing by next-generation sequencing (NGS) of a broad gene panel, as a supplement to the traditional forensic investigation in cases with non-diagnostic structural abnormalities of the heart. METHODS AND RESULTS: We screened 72 suspected SCD cases (<50 years) using the HaloPlex Target Enrichment System (Agilent) and NGS (Illumina MiSeq) for 100 genes previously associated with inherited cardiomyopathies and channelopathies. Fifty-two cases had non-diagnostic structural cardiac abnormalities and 20 cases, diagnosed with a cardiomyopathy post-mortem (ARVC = 14, HCM = 6), served as comparators. Fifteen (29%) of the deceased individuals with non-diagnostic findings had variants with likely functional effects based on conservation, computational prediction, allele-frequency and supportive literature. The corresponding frequency in deceased individuals with cardiomyopathies was 35% (p = 0.8). CONCLUSION: The broad genetic screening revealed variants with likely functional effects at similar high rates, i.e. in 29 and 35% of the suspected SCD cases with non-diagnostic and diagnostic cardiac abnormalities, respectively. Although the interpretation of broad NGS screening is challenging, it can support the forensic investigation and help the cardiologist's decision to offer counselling and clinical evaluation to relatives of young SCD victims.

Forensic ancestry analysis with two capillary electrophoresis ancestry informative marker (AIM) panels: Results of a collaborative EDNAP exercise.

There is increasing interest in forensic ancestry tests, which are part of a growing number of DNA analyses that can enhance routine profiling by obtaining additional genetic information about unidentified DNA donors. Nearly all ancestry tests use single nucleotide polymorphisms (SNPs), but these currently rely on SNaPshot single base extension chemistry that can fail to detect mixed DNA. Insertion-deletion polymorphism (Indel) tests have been developed using dye-labeled primers that allow direct capillary electrophoresis detection of PCR products (PCR-to-CE). PCR-to-CE maintains the direct relationship between input DNA and signal strength as each marker is detected with a single dye, so mixed DNA is more reliably detected. We report the results of a collaborative inter-laboratory exercise of 19 participants (15 from the EDNAP European DNA Profiling group) that assessed a 34-plex SNP test using SNaPshot and a 46-plex Indel test using PCR-to-CE. Laboratories were asked to type five samples with different ancestries and detect an additional mixed DNA sample. Statistical inference of ancestry was made by participants using the Snipper online Bayes analysis portal plus an optional PCA module that analyzes the genotype data alongside calculation of Bayes likelihood ratios. Exercise results indicated consistent genotyping performance from both tests, reaching a particularly high level of reliability for the Indel test. SNP genotyping gave 93.5% concordance (compared to the organizing laboratory's data) that rose to 97.3% excluding one laboratory with a large number of miscalled genotypes. Indel genotyping gave a higher concordance rate of 99.8% and a reduced no-call rate compared to SNP analysis. All participants detected the mixture from their Indel peak height data and successfully assigned the correct ancestry to the other samples using Snipper, with the exception of one laboratory with SNP miscalls that incorrectly assigned ancestry of two samples and did not obtain informative likelihood ratios for a third. Therefore, successful ancestry assignments were achieved by participants in 92 of 95 Snipper analyses. This exercise demonstrates that ancestry inference tests based on binary marker sets can be readily adopted by laboratories that already have well-established CE regimes in place. The Indel test proved to be easy to use and allowed all exercise participants to detect the DNA mixture as well as achieving complete and concordant profiles in nearly all cases. Lastly, two participants successfully ran parallel next-generation sequencing analyses (each using different systems) and achieved high levels of genotyping concordance using the exercise PCR primer mixes unmodified.

Population and forensic data for three sets of forensic genetic markers in four ethnic groups from Iran: Persians, Lurs, Kurds and Azeris.

A total of 255 individuals (Persians, Lurs, Kurds and Azeris) from Iran were typed for three sets of forensic genetic markers with the NGM SElect™, DIPplex(®) and Argus X-12 kits. Statistically significant deviations (P≤0.002) from Hardy-Weinberg expectations were observed for the insertion-deletion markers HLD97 and HLD93 after Holm-Sidák correction. Statistically significant (P<0.05) levels of linkage disequilibrium were observed between markers within two of the four studied X-chromosomal linkage groups. AMOVA analyses of the three sets of markers did not show population structure when the individuals were grouped according to their ethnic group. The Iranian population grouped closely to populations living geographically near to Iran based on pairwise FST distances. The matching probabilities ranged from 1 in 3.2×10(7) males by using haplotype frequencies of four X-chromosomal haplogroups to 1 in 3.4×10(21) individuals for the 16 autosomal STRs.

Results for five sets of forensic genetic markers studied in a Greek population sample.

A population sample of 223 Greek individuals was typed for five sets of forensic genetic markers with the kits NGM SElect™, SNPforID 49plex, DIPplex®, Argus X-12 and PowerPlex® Y23. No significant deviation from Hardy-Weinberg expectations was observed for any of the studied markers after Holm-Sidák correction. Statistically significant (P<0.05) levels of linkage disequilibrium were observed between markers within two of the studied X-chromosome linkage groups. AMOVA analyses of the five sets of markers did not show population structure when the individuals were grouped according to their geographic origin. The Greek population grouped closely to the other European populations measured by F(ST)(*) distances. The match probability ranged from a value of 1 in 2×10(7) males by using haplotype frequencies of four X-chromosome haplogroups in males to 1 in 1.73×10(21) individuals for 16 autosomal STRs.