Non-invasive prenatal paternity testing using a standard forensic genetic massively parallel sequencing assay for amplification of human identification SNPs.

Prenatal paternity testing often relies on invasive procedures that cause risk to both the mother and the foetus. Non-invasive, prenatal paternity testing by investigating paternally inherited single nucleotide polymorphisms (SNPs) in cell-free foetal DNA (cffDNA) in maternal plasma was performed at consecutive time points during early gestation. Plasma from 15 pregnant women was investigated at consecutive time points from gestational weeks (GWs) 4-20. The Precision ID Identity Panel and an Ion S5 Sequencer was used to analyse the cffDNA. Paternally inherited foetal SNP alleles were detected from GW7. The median foetal fractions were 0%, 3.9%, 5.1%, 5.2%, and 4.7% at GWs 4, 7, 12, 16, and 20, respectively. The corresponding median numbers of detected paternally inherited foetal autosomal SNP alleles were 0, 3, 9, 10, and 12, respectively. The typical (i.e. geometric mean) paternity indices at GW12 and GW20 were 24 (range 0.0035-8389) and 199 (range 5.1-30,137), respectively. The method is very promising. However, the method can be improved by shortening the lengths of the PCR amplicons and increasing the number of SNPs. To our knowledge, this is the first study to successfully identify paternally inherited foetal SNP alleles at consecutive time points in early gestation independently of the foetal gender.

High Y-chromosomal Differentiation Among Ethnic Groups of Dir and Swat Districts, Pakistan.

The ethnic groups that inhabit the mountainous Dir and Swat districts of northern Pakistan are marked by high levels of cultural and phenotypic diversity. To obtain knowledge of the extent of genetic diversity in this region, we investigated Y-chromosomal diversity in five population samples representing the three main ethnic groups residing within these districts, including Gujars, Pashtuns and Kohistanis. A total of 27 Y-chromosomal short tandem repeats (Y-STRs) and 331 Y-chromosomal single nucleotide polymorphisms (Y-SNPs) were investigated. In the Y-STRs, we observed very high and significant levels of genetic differentiation in nine of the 10 pairwise between-group comparisons (RST 0.179-0.746), and the differences were mirrored in the Y-SNP haplogroup frequency distribution. No genetic differences were found between the two Pashtun subethnic groups Tarklanis and Yusafzais (RST = 0.000). Utmankhels, also considered Pashtuns culturally, were not closely related to any of the other population samples (RST 0.451-0.746). Thus, our findings provide examples of both associations and dissociations between cultural and genetic legacies. When analyzed within a larger continental-scale context, these five ethnic groups fall mostly outside the previously characterized Y-chromosomal gene pools of the Indo-Pakistani subcontinent. Male founder effects, coupled with culturally and topographically based constraints upon marriage and movement, are likely responsible for the high degree of genetic structure in this region.

ISO 17025 validation of a next-generation sequencing assay for relationship testing.

The HID-Ion AmpliSeq™ Identity Panel is a next-generation sequencing assay with 90 autosomal and 34 Y-chromosome SNPs that are amplified in one PCR step and subsequently sequenced using the Ion Personal Genome Machine (Ion PGM™) System. This assay was validated for relationship testing in our ISO 17025 accredited laboratory in 2015. Here, the essential parts of the validation report submitted to the Danish Accreditation Fund are presented. A total of 100 unrelated Danes were typed in duplicates and the locus balance, heterozygote balance (Hb) and noise levels were analysed in detail. Two loci were disregarded for casework because genotyping was uncertain. Hb for rs7520386 was skewed and high levels of noise were observed in rs576261. Three general acceptance criteria for analysis of single-source samples were defined: (i) sequencing depth > 200 reads, (ii) noise level < 3% and (iii) Hb > 0.3. A Python script named SNPonPGM was developed to assist the analyst by highlighting loci that do not fulfil the general acceptance criteria. Furthermore, SNPonPGM has functions that reduce the hands-on time of the reporting officer to a few minutes per case. Mixtures with DNA from two individuals in a 1:24 ratio were readily identified using the three criteria and the SNPonPGM script.

High-throughput sequencing of forensic genetic samples using punches of FTA cards with buccal swabs.

Here, we demonstrate that punches from buccal swab samples preserved on FTA cards can be used for high-throughput DNA sequencing, also known as massively parallel sequencing (MPS). We typed 44 reference samples with the HID-Ion AmpliSeq Identity Panel using washed 1.2 mm punches from FTA cards with buccal swabs and compared the results with those obtained with DNA extracted using the EZ1 DNA Investigator Kit. Concordant profiles were obtained for all samples. Our protocol includes simple punch, wash, and PCR steps, reducing cost and hands-on time in the laboratory. Furthermore, it facilitates automation of DNA sequencing.

Introduction of the Python script STRinNGS for analysis of STR regions in FASTQ or BAM files and expansion of the Danish STR sequence database to 11 STRs.

This work introduces the in-house developed Python application STRinNGS for analysis of STR sequence elements in BAM or FASTQ files. STRinNGS identifies sequence reads with STR loci by their flanking sequences, it analyses the STR sequence and the flanking regions, and generates a report with the assigned SNP-STR alleles. The main output file from STRinNGS contains all sequences with read counts above 1% of the total number of reads per locus. STR sequences are automatically named according to the nomenclature used previously and according to the repeat unit definitions in STRBase (http://www.cstl.nist.gov/strbase/). The sequences are named with (1) the locus name, (2) the length of the repeat region divided by the length of the repeat unit, (3) the sequence(s) of the repeat unit(s) followed by the number of repeats and (4) variations in the flanking regions. Lower case letters in the main output file are used to flag sequences with previously unknown variations in the STRs. SNPs in the flanking regions are named by their "rs" numbers and the nucleotides in the SNP position. Data from 207 Danes sequenced with the Ion Torrent™ HID STR 10-plex that amplified nine STRs (CSF1PO, D3S1358, D5S818, D7S820, D8S1179, D16S539, TH01, TPOX, vWA), and Amelogenin was analysed with STRinNGS. Sequencing uncovered five common SNPs near four STRs and revealed 20 new alleles in the 207 Danes. Three short homopolymers in the D8S1179 flanking regions caused frequent sequencing errors. In 29 of 3726 allele calls (0.8%), sequences with homopolymer errors were falsely assigned as true alleles. An in-house developed script in R compensated for these errors by compiling sequence reads that had identical STR sequences and identical nucleotides in the five common SNPs. In the output file from the R script, all SNP-STR haplotype calls were correct. The 207 samples and six additional samples were sequenced for D3S1358, D12S391, and D21S11 using the 454 GS Junior platform in this and a previous work. Overall, next generation sequencing (NGS) of the 11 STRs lowered the mean match probability 386 times and increased the typical paternity indexes (i.e. the geometric mean) for trios and duos 47 and 23 times, respectively, compared to the traditional PCR-CE typing of the same population.

Forensic and population genetic analyses of Danes, Greenlanders and Somalis typed with the Yfiler® Plus PCR amplification kit.

Recently, the Yfiler® Plus PCR Amplification Kit (Yfiler® Plus, Thermo Fisher Scientific, Waltham, MA, USA) was introduced. Yfiler® Plus amplifies 27 Y-chromosomal short tandem repeat loci (Y-STRs) and adds ten new Y-STRs to those analysed with the commonly used AmpFlSTR® Yfiler® PCR Amplification Kit (Yfiler®, Thermo Fisher Scientific, Waltham, MA, USA). Seven of the new Y-STRs are rapidly mutating Y-STRs (RM Y-STRs). In this study, 551 male individuals from Denmark, Greenland and Somalia were typed with Yfiler® Plus. The results were compared to those obtained with Yfiler® in the same individuals. Forensic and population genetic parameters were estimated for Yfiler® Plus. Yfiler® Plus had a higher power of discrimination than Yfiler® in all three populations. Compared to Yfiler®, Yfiler® Plus offers increased power of discrimination, which is obviously an advantage in crime case investigations. However, the inclusion of seven RM Y-STRs in Yfiler® Plus makes it less attractive for relationship testing because of the relatively high combined mutation rate, approximately 15%.

Gene dose effects of GSTM1, GSTT1 and GSTP1 polymorphisms on outcome in childhood acute lymphoblastic leukemia.

Children with acute lymphoblastic leukemia (ALL) react very differently to chemotherapy. One explanation for this is inherited genetic variation. The glutathione S-transferase (GST) enzymes inactivate a number of chemotherapeutic drugs administered in childhood ALL therapy. Two multiplexing methods were applied for genotyping the GSTM1 and GSTT1 genes (distinguishing between 0, 1, or 2 gene copies) and the GSTP1 313 A>G polymorphism, simultaneously. A total of 263 childhood ALL patients were genotyped. No gene dose effect on outcome was demonstrated with either GST polymorphisms. Grouping of GSTM1 and GSTT1 into poor (0 or 1 gene copy)-and good metabolizers (at least 2 gene copies)-showed that the poor metabolizers had a trend toward a better outcome (event-free survival =91.8%) compared with the good metabolizers (event-free survival =83.2%). Similarly, in the adjusted analysis the good metabolizers demonstrated a 2.2-fold higher risk trend of experiencing an event (resistant disease or relapse) compared with the poor metabolizers (P=0.066; hazard ratio =2.248; 95% confidence interval, 0.948-5.327). In conclusion, our results suggest that the combined gene dose of GSTM1 and GSTT1 may influence outcome in childhood ALL.

The role of the glutathione S-transferase genes GSTT1, GSTM1, and GSTP1 in acetaminophen-poisoned patients.

The aim of this study was to assess if genetic variants in the glutathione-S-transferase genes GST-T1, M1, and P1 reflect risk factors in acetaminophen (APAP)-poisoned patients assessed by investigation of the relation to prothrombin time (PT), which is a sensitive marker of survival in these patients. A total of 104 APAP-poisoned patients were genotyped for deletion polymorphisms in the GSTT1 and GSTM1 genes and for the GSTP1 Ile105Val polymorphism. We found a borderline association (p = 0.05) between the GSTT1 homozygous deletion genotype and high trough PT (a marker of prognosis in APAP poisoning) compared to carrying two functioning copies of the gene. No significant association was found between any of the GSTM1 and GSTP1 genotypes and PT. The frequency of GSTP1 Val/Val genotypes was significantly lower in the patients than in the background population (p = 0.047). The results suggest that the GSTT1 homozygous deletion genotype may be associated with a better prognosis after APAP poisoning and that carriers of the GSTP1 homozygous variant genotype may have a decreased risk of being APAP poisoned.

Postmortem blood concentrations of R- and S-enantiomers of methadone and EDDP in drug users: influence of co-medication and p-glycoprotein genotype.

We investigated toxicological and pharmacogenetic factors that could influence methadone toxicity using postmortem samples. R- and S-methadone were measured in femoral blood from 90 postmortem cases, mainly drug users. The R-enantiomer concentrations significantly exceeded that of the S-enantiomers (Wilcoxon's test, p < 0.001). The samples were divided into four groups according to other drugs detected (methadone only, methadone and strong analgesics, methadone and benzodiazepines, or methadone and other drugs). There was no significant difference in any of the R-methadone/total methadone ratios among the four groups. The median R/S ratio was 1.38, which tends to be higher than that reported for the plasma of living subjects. In addition, we investigated whether small nucleotide polymorphisms in the MDR1 gene that encode the drug transporter P-glycoprotein were associated with the concentrations of R- and S-methadone and its metabolite 2-ethylidene-1,5-dimethyl-3,3-diphenylpyrrolidine. No significant association was detected.

Dominance of pre-analytical over analytical variation for measurement of methadone and its main metabolite in postmortem femoral blood.

On the basis of simultaneously sampled postmortem blood specimens from the left and right femoral veins the pre-analytical variation of methadone measurements was evaluated and compared to the analytical variation. The material consisted of a series of 27 duplicate samples from routine autopsy cases comprising mainly drug addicts. A chiral LC-MS/MS method was used for measurement of the R- and S-enantiomers of methadone and its main metabolite 2-ethyl-1,5-dimethyl-3,3-diphenylpyrrolinium (EDDP). The analytical CV% was determined to be in the range 3-4% for methadone enantiomers and 4-6% for EDDP enantiomers. The total measurement uncertainty (CV(T)) was estimated from the pre-analytical variation (CV(PA)), analytical variation proper (CV(A)), and variation related to calibration (traceability) (CV(Cal)) according to the relationship CV(T) = [CV(2)(PA) + CV(2)(A) + CV(2)(cal)](0.5). Uncertainty related to calibration concerned a component related to the purity of drug reference compound and a contribution from the production of calibrator solutions (CV(Cal)<1%). Pre-analytical sampling variation was estimated from the duplicate measurements of blood samples after subtraction of the analytical component. The pre-analytical variation amounted to a CV% of 19-21% for R- and S-methadone and 30-38% for R- and S-EDDP, i.e. considerably larger than the other components. Due to the squared addition principle, the resulting total uncertainty (CV(T)) became largely identical to the CV(PA), i.e. 19-21% for R- and S-methadone and 31-38% for R- and S-EDDP enantiomers. Accordingly, CV(T) exceeded CV(A) by a factor 5 or more. Dominance of the pre-analytical component of variation may also be likely for other compounds measured in postmortem blood samples. Thus, the width of the 95%-uncertainty interval (+/-2CV(T)) for a postmortem measurement is largely determined by the pre-analytical component of variation. This should be kept in mind when judging on the uncertainty of postmortem measurement results.

Multiplex PCR detection of GSTM1, GSTT1, and GSTP1 gene variants: simultaneously detecting GSTM1 and GSTT1 gene copy number and the allelic status of the GSTP1 Ile105Val genetic variant.

The glutathione S-transferase (GST) genes GSTM1, GSTT1, and GSTP1 are involved in the detoxification of a broad range of toxic substances. Genetic polymorphisms in these genes have been studied intensively for their potential role in cancer susceptibility and drug response. In Caucasians, the enzyme activity of GSTM1 and GSTT1 is absent in approximately 50 and 15% of the population, respectively, due to deletions of both chromosomal copies of the genes. A trimodal phenotype pattern exists in which individuals with two, one, or no functional genes are fast, intermediate, or slow "conjugators," respectively. Most studies investigating the effect of the GSTM1 and GSTT1 deletions do not distinguish between fast and intermediate conjugators because the applied genotyping assays only detect if at least one copy of either gene is present. We present a multiplex PCR assay that detects if an individual has none, one, or two copies of the GSTM1 and GSTT1 genes and simultaneously detects the allelic status of the GSTP1 Ile105Val genetic variant. A total of 200 Danes, 100 Somalis, and 100 Greenlanders were genotyped. This multiplex PCR assay enables future large-scale studies to investigate the role of GSTs.

Multiplex PCR and minisequencing of SNPs--a model with 35 Y chromosome SNPs.

We have developed a robust single nucleotide polymorphism (SNPs) typing assay with co-amplification of 25 DNA-fragments and the detection of 35 human Y chromosome SNPs. The sizes of the PCR products ranged from 79 to 186 base pairs. PCR primers were designed to have a theoretical Tm of 60 +/- 5 degrees C at a salt concentration of 180 mM. The sizes of the primers ranged from 19 to 34 nucleotides. The concentration of amplification primers was adjusted to obtain balanced amounts of PCR products in 8mM MgCl2. For routine purposes, 1 ng of genomic DNA was amplified and the lower limit was approximately 100 pg DNA. The minisequencing reactions were performed simultaneously for all 35 SNPs with fluorescently labelled dideoxynucleotides. The size of the minisequencing primers ranged from 19 to 106 nucleotides. The minisequencing reactions were analysed by capillary electrophoresis and multicolour fluorescence detection. Female DNA did not influence the results of Y chromosome SNP typing when added in concentrations more than 300 times the concentrations of male DNA. The frequencies of the 35 SNPs were determined in 194 male Danes. The gene diversity of the SNPs ranged from 0.01 to 0.5.