Estimation of the number of contributors to mixed samples of DNA by mitochondrial DNA analyses using massively parallel sequencing.

We evaluated whether the number of contributors to mixed DNA samples can be estimated by analyzing the D-loop of mitochondrial DNA using massively parallel sequencing. The A- (positions 16,209-16,400) and B- (positions 30-284) amplicons in hypervariable regions 1 and 2, respectively, were sequenced using MiSeq with 2 × 251 cycles. Sequence extraction and trimming were performed using CLC Genomics Workbench 11 and the number of observed haplotypes was counted for each amplicon type using Microsoft Excel. The haplotype ratios were calculated by dividing the number of counted reads of the corresponding haplotype by the total number of sequence reads. Haplotypes that were over the threshold (5%) were defined as positive haplotypes. The number of larger positive haplotypes in either of the two amplicon types was defined as the number of contributors. Samples were collected from seven individuals. Seventeen mixed samples were prepared by mixing DNA from two to five contributors at various ratios. The number of contributors was correctly estimated from almost all of the mixed samples containing equal amounts of DNA from two to five people. In mixed samples of two or three people, the minor components were detected down to a ratio of 20:1 or 8:2:1. However, heteroplasmy, base deletions, and sharing of the same haplotypes caused incorrect estimations of the number of contributors. Although this method still has room for improvement, it may be useful for estimating the number of contributors in a mixed sample, as it does not rely on forensic mathematics.

'Grasshopper sign': the novel imaging of post-COVID-19 myelopathy with delayed longitudinal white matter abnormalities.

Introduction: Recently, there have been a few reports of atypical post-coronavirus disease 2019 (COVID-19) myelopathy manifesting tract-specific lesions similar to those due to vitamin B12 deficiency. However, the precise characteristics of imaging or clinical course remain not well understood. Methods: A retrospective analysis of the clinical and imaging characteristics of four patients who were referred to our hospital with a unique post-COVID-19 myelopathy was performed. Results: Four-to-six weeks following COVID-19 infection in the summer of 2023, four middle-aged men developed paraparesis, hypo/dysesthesia and bladder/bowel disturbance, suggesting myelopathy. Although spinal MRI showed no abnormalities in the early stages, tract-specific longitudinal lesions along the dorsal and lateral columns became apparent as the symptoms progressed. Owing to the lack of MRI findings at the early stage, all cases were challenging to diagnose. However, the patients remained partially responsive to aggressive immunosuppressive therapies, even in the advanced stage. Discussion: We termed these tract-specific longitudinal lesions in the presented case series 'Grasshopper sign' because brain coronal and spine axial MRI findings looked like a grasshopper's antennae and face. Early identification of the characteristic MRI abnormality could allow for early intervention using intensive immunosuppressive therapy, which could improve patient outcomes.

A Case of Long-Term Exposure to Valproic Acid Mimicking Tremor-Dominant Parkinson's Disease.

Background: Valproic acid is associated with increased risks of tremor and parkinsonism. Case Report: A 67-year-old man with a diagnosis of epilepsy who had been treated with valproic acid (VPA) for 32 years noticed right-dominant upper-limb resting tremor accompanied by mild rigidity and bradykinesia. He was initially diagnosed with tremor-dominant Parkinson's disease (TDPD), but dopamine transporter single-photon emission computed tomography demonstrated no nigrostriatal degeneration. At 3 months after discontinuing VPA, his symptoms dramatically improved. Discussion: VPA-induced tremor usually consists of postural or kinetic tremor without asymmetry. Our case indicated that careful evaluation is needed, even in cases of asymmetrical resting tremor and mild parkinsonism resembling TDPD after long term exposure to VPA. Highlights: We report an atypical case of valproic acid-induced tremor and parkinsonism that mimics tremor-dominant Parkinson's disease. Physicians should not exclude the possible relation to valproic acid in patients presenting unilateral resting tremor and parkinsonism even in the absence of long-term side effects.

Poly_NumtS_430 or HSA_NumtS_587 observed in massively parallel sequencing of the mitochondrial HV1 and HV2 regions.

An increasing number of studies on massively parallel sequencing of mitochondrial DNA (mtDNA) have been reporting identification of various types of noise or off-target sequences. Herein, we report that an off-target haplotype (sequence length 192 bp) observed in MiSeq data of mtDNA at nucleotide position 16,209-16,400 was likely caused by polymorphic nuclear mitochondrial DNA sequences (NumtS). Buccal DNA samples from Volunteers #001-004 and Control DNA 007 were amplified with our multiplex system of the B (15,998-16,172), C (16,209-16,400), and E (30-289) regions using 2000 copies of mtDNA. A sample index was added using a Nextera XT index kit, and MiSeq Reagent Nano Kit v2 was used in 2 × 251 cycles on a MiSeq FGx. FASTQ files were analyzed by CLC Genomics Workbench 21.0.3. The extracted SAM files were analyzed using our original Excel macro to sum up the read counts as the phased variant calls for each region. An off-target haplotype differing at 19 sites against the revised Cambridge reference sequence was observed in Volunteer #001 (4 in 10 MiSeq data), Volunteer #002 (2 in 9), and Control DNA 007 (6 in 9). In a BLAST search, the sequence of the off-target haplotype matched perfectly to three polymorphic NumtS (Poly_NumtS_430 [KM281528.1], HSA_NumtS_587 [HE613849.1], and nuclear fossil [S80333.1]) and BAC clone of chromosome 11 (AC107937.2). The sequence also matched perfectly to a Filipino mtDNA (KC993973.1) which was inferred as a chimeric sequence of mtDNA and the HSA_NumtS_587. The sequence of the off-target haplotype was not contained in the latest human reference genome sequence (GRCh38.p13). In a phylogenetic tree, the off-target haplotype was genetically distant from modern human mtDNA and not directly connected to them. In conclusion, observed off-target haplotype amplified by our multiplex system was likely caused by Poly_NumtS_430 or HSA_NumtS_587.

Development and validation of a SYBR green-based mitochondrial DNA quantification method by following the MIQE and other guidelines.

In forensic mitochondrial DNA (mtDNA) analysis, quantitative PCR (qPCR) is usually performed to obtain high-quality sequence data for subsequent Sanger or massively parallel sequencing. Unlike methods for nuclear DNA quantification using qPCR, a calibrator is necessary to obtain mtDNA concentrations (i.e., copies/µL). Herein, we developed and validated a mtDNA quantification method based on a SYBR Green assay by following MIQE [Bustin et al., Clin. Chem. 55 (2009) 611-22] and other guidelines. Primers were designed to amplify nucleotide positions 16,190-16,420 in hypervariable region 1 for qPCR using PowerUp SYBR Green and QuantStudio 5. The optimized conditions were 0.3 µM each primer and an annealing temperature of 60 °C under a 2-step cycling protocol. K562 DNA at 100 pg/µL was converted into a mtDNA concentration of 16,400 copies/µL using linearized plasmid DNA. This mtDNA calibrator was obtained by cloning the synthesized DNA fragments of mtDNA (positions 16,140-16,470) containing a 100-bp inversion. The linear dynamic range of the K562 standard curve was 10,000-0.1 pg/µL (r2 ≥ 0.999). The accuracy was examined using NIST SRM 2372a, and its components A, B, and C were quantified with differences of -29.4%, -35.0%, and -22.0%, respectively, against the mtDNA concentrations calculated from published NIST data. We also examined the specificity of the primers, stability of the reaction mix, precision, tolerance against PCR inhibitors, and cross-reactivity against DNA from various animal taxa. Our newly developed mtDNA quantification method is expected to be useful for forensic mtDNA analysis.

Development and validation of Kongoh ver. 3.0.1: Open-source software for DNA mixture interpretation in the GlobalFiler system based on a quantitative continuous model.

Probabilistic genotyping software based on continuous models is effective for interpreting DNA profiles derived from DNA mixtures and small DNA samples. In this study, we updated our previously developed Kongoh software (to ver. 3.0.1) to interpret DNA profiles typed using the GlobalFiler™ PCR Amplification Kit. Recently, highly sensitive typing systems such as the GlobalFiler system have facilitated the detection of forward, double-back, and minus 2-nt stutters; therefore, we implemented statistical models for these stutters in Kongoh. In addition, we validated the new version of Kongoh using 2-4-person mixtures and DNA profiles with degradation in the GlobalFiler system. The likelihood ratios (LRs) for true contributors and non-contributors were well separated as the information increased (i.e., larger peak height and fewer contributors), and these LRs tended to neutrality as the information decreased. These trends were observed even in profiles with DNA degradation. The LR values were highly reproducible, and the accuracy of the calculation was also confirmed. Therefore, Kongoh ver. 3.0.1 is useful for interpreting DNA mixtures and degraded DNA samples in the GlobalFiler system.

Massively parallel sequencing data of 31 autosomal STR loci obtained using the Precision ID GlobalFiler NGS STR Panel v2 for 82 Japanese population samples.

Allele frequencies for 31 autosomal short tandem repeat (STR) loci (CSF1PO, D10S1248, D12ATA63, D12S391, D13S317, D14S1434, D16S539, D18S51, D19S433, D1S1656, D1S1677, D21S11, D22S1045, D2S1338, D2S1776, D2S441, D3S1358, D3S4529, D4S2408, D5S2800, D5S818, D6S1043, D6S474, D7S820, D8S1179, FGA, Penta D, Penta E, TH01, TPOX, and vWA) were obtained using Precision ID GlobalFiler NGS STR Panel v2 from 82 unrelated individuals sampled from the Japanese population. Autosomal STR alleles designated by NGS and conventional capillary electrophoresis were found to be concordant except at D2S441 allele 9.

Evaluation of probability distribution models for stutter ratios in the typing system of GlobalFiler and 3500xL Genetic Analyzer.

As DNA typing systems have become increasingly sensitive in recent years, probability distribution models for back, forward, double-back, and minus 2-nt stutter ratios have been desired to be considered in DNA evidence interpretation using specific software programs. However, experimental investigations have been insufficient, especially for forward, double-back, and minus 2-nt stutters. In this study, we experimentally reevaluated the probability distribution models for each stutter ratio in the typing systems of GlobalFiler™ PCR Amplification Kit and 3500xL Genetic Analyzer from Thermo Fisher Scientific. In addition, to enhance the reliability of longest uninterrupted stretch (LUS) values and corrected allele numbers used in previously developed models for stutter ratios using sequence information (i.e., LUS model and multi-seq model), we propose the weighted average of LUS values and corrected allele numbers based on the number of observations in sequence-based population data. Back stutter ratios demonstrated a positive correlation with allele numbers (allele model) in eight loci, LUS values (LUS model) in eight loci, and corrected allele numbers (multi-seq model) in five loci. The forward stutter ratios (FSRs) of D22S1045 followed the LUS model. FSRs other than D22S1045 and double-back stutter ratios followed the LUS model by considering multiple loci together. Minus 2-nt stutter ratios observed in SE33 and D1S1656 did not increase with the increase in the allele numbers. The adopted models for each stutter ratio can be implemented in software programs for DNA evidence interpretation and enable a reliable interpretation of crime stain profiles in forensic caseworks.

Evaluation of Rapid DNA system for buccal swab and disaster victim identification samples.

An evaluation of a Rapid DNA system was performed using buccal swab samples and mock Disaster Victim Identification (DVI) samples collected postmortem. The allelic ladder success rate was 90% and samples analyzed simultaneously with this allelic ladder were used for further analysis. Sample success rate of the Rapid DNA system for buccal swab samples, and blood and muscle DVI samples were calculated. Success rates of buccal swab samples were 100% and 75% using cassettes preloaded with all reagents suitable for high- and low-DNA content samples, respectively. Success rates of fresh DVI samples were 80% to 100%. Success rates of putrefied DVI samples varied widely between 0% and 20% and 50% to 80% depending on cassette and sample types. Conventional DNA analysis was performed for comparison with the results of the Rapid DNA system. DNA quantity and degradation of human DNA were measured using quantitative polymerase chain reaction. DVI samples that yielded more than 1 ng/µL of DNA when extracted with conventional protocols were suitable for analysis using cassettes for both high- and low-DNA content samples. DVI samples with less than 0.1 ng/µL of DNA were suitable only for analysis using cassettes for low-DNA content samples. All alleles called and exported by the Expert system software implemented in the Rapid DNA system were concordant with allele calls made by conventional capillary electrophoresis DNA analysis.

Heteroplasmies detected in an amplified mitochondrial DNA control region from a small amount of template.

When mitochondrial DNA (mtDNA) heteroplasmies are detected, they often confound forensic identification, especially if they are the result of poor biological sampling. In this study, we determined the ratio of heteroplasmy in samples that were amplified from a very small amount of template mtDNA or a few cells using a highly sensitive nested polymerase chain reaction (PCR) procedure and a direct sequencing analysis. As a result, more than half of the detected sequences (i.e., 17/20, 15/20, and 14/20) showed homoplasmy derived from a variation in the heteroplasmy proportion when only 10 copies of template mtDNA samples were amplified and analyzed. Additionally, with products amplified from one or several white blood cells (WBCs), several previously undetected heteroplasmies were detected. These results indicate the risks associated with using highly sensitive mtDNA techniques in forensic investigations because of the variable proportions of heteroplasmy or nucleotide substitutions that can possibly be detected from a very small biological sample.

Non-specific peaks generated by animal DNA during human STR analysis: Peak characteristics and a novel analysis method for mixed human/animal samples.

Forensic human identification (HID) laboratories occasionally encounter non-specific peaks generated by non-human DNA. Casework samples for human short tandem repeat (STR) profiling may be contaminated by animal DNA because of the specific environment or situation from which they were obtained. Validation studies for HID kits have reported that non-specific peaks generated from some animals are observed near the human amelogenin peak. In this study, we first revealed that DNA sequences associated with the non-specific peaks generated from animal DNA differ from one animal family to the other. However, non-specific peaks cannot be analyzed using the remainder of polymerase chain reaction (PCR) products left over from conventional HID kits when human and animal DNA are mixed. To overcome this issue, we have developed a novel analysis method of using non-specific peaks generated from animal DNA in human STR profiling to identify the source of contaminating animal DNA at the family level. The method applied here is termed as blocking PCR, which involves selective animal DNA re-amplification by blocking nontarget human amelogenin DNA amplification using an oligonucleotide probe that specifically binds to human amelogenin using the remaining PCR product from the HID kit. Our data demonstrated that HID and family discrimination among animals that are often encountered in forensic contexts could be performed simultaneously. This study enabled recovery of more information from limited quantities of casework samples contaminated with animal DNA, which would be useful for forensic HID scientists.

D5S818 Typing Discrepancy Between PowerPlex(®) Fusion and Other STR Kits Including GlobalFiler(®) Caused by a One-base Deletion in 31 Nucleotides Upstream of the Repeat Region.

Short tandem repeat (STR) typing is widely used in forensic investigation. When the same DNA sample is analyzed with different STR typing kits, a typing discrepancy is occasionally observed. In this study, we examined the cause of a typing discrepancy in a sample at D5S818 locus. This sample was designated as 10, 12 using Identifiler(®) , Identifiler(®) Plus, GlobalFiler(®) , PowerPlex(®) 16HS, and PowerPlex(®) 18D, but as 9.3, 12 using PowerPlex(®) Fusion. Sequencing results indicated that the shorter allele in the sample had a deletion (U31Tdel) at 31 nucleotides upstream of the repeat region (AGAT)10 . This deletion was located in the binding site of the published D5S818 forward primer in PowerPlex(®) 16 and was only 9 and 11 nucleotides downstream of our estimated 5' end position of D5S818 forward primer in GlobalFiler(®) and PowerPlex(®) 18D, respectively. We also examined the effect of primer length on the heterozygous peak balance in this sample.

Estimating allele dropout probabilities by logistic regression: Assessments using Applied Biosystems 3500xL and 3130xl Genetic Analyzers with various commercially available human identification kits.

Phenomena called allele dropouts are often observed in crime stain profiles. Allele dropouts are generated because one of a pair of heterozygous alleles is underrepresented by stochastic influences and is indicated by a low peak detection threshold. Therefore, it is important that such risks are statistically evaluated. In recent years, attempts to interpret allele dropout probabilities by logistic regression using the information on peak heights have been reported. However, these previous studies are limited to the use of a human identification kit and fragment analyzer. In the present study, we calculated allele dropout probabilities by logistic regression using contemporary capillary electrophoresis instruments, 3500xL Genetic Analyzer and 3130xl Genetic Analyzer with various commercially available human identification kits such as AmpFℓSTR® Identifiler® Plus PCR Amplification Kit. Furthermore, the differences in logistic curves between peak detection thresholds using analytical threshold (AT) and values recommended by the manufacturer were compared. The standard logistic curves for calculating allele dropout probabilities from the peak height of sister alleles were characterized. The present study confirmed that ATs were lower than the values recommended by the manufacturer in human identification kits; therefore, it is possible to reduce allele dropout probabilities and obtain more information using AT as the peak detection threshold.

Sibling assessment based on likelihood ratio and total number of shared alleles using 21 short tandem repeat loci included in the GlobalFiler™ kit.

Sibling assessment using the 15 autosomal short tandem repeat (STR) loci included in the Identifiler® kit can be difficult when comparing an unidentified party to an alleged sibling. Therefore, we investigated the likelihood ratio (LR) and the total number of shared alleles (TNSA) for sibship determination using the 21 autosomal STR loci included in the GlobalFiler™ kit. We computationally generated the genotypes of 10,000 sibling pairs and 10,000 unrelated pairs based on previously reported allele frequencies of the 15 Identifiler loci and the remaining 6 GlobalFiler loci. The LR and the TNSA were then calculated in each pair using the 15 and 21 loci. Next, these calculations were applied to 22 actual sibling pairs. LR values ⩾ 10,000 were observed in 48% of the sibling pairs using the 15 loci and in 80% of the sibling pairs using the 21 loci. The TNSA distribution between siblings and unrelated pairs was more divergent in GlobalFiler than in Identifiler. TNSA values ⩾ 20 were found only in true siblings in Identifiler, while TNSA values ⩾24 in GlobalFiler. In Identifiler, all pairs with TNSA ⩾ 24 had LR values ⩾ 10,000 and the same was true in GlobalFiler for TNSA ⩾29. Therefore, increasing the number of loci is very efficient for sibship determination. The LR is most reliable for determining sibship. However, TNSA values may be useful for the preliminary method of LR values because LR value demonstrated a significantly positive correlation with TNSA value in both Identifiler and GlobalFiler.

Next-generation sequencing analysis of off-ladder alleles due to migration shift caused by sequence variation at D12S391 locus.

In short tandem repeat (STR) analysis, length polymorphisms are detected by capillary electrophoresis (CE). At most STR loci, mobility shift due to sequence variation in the repeat region was thought not to affect the typing results. In our recent population studies of 1501 Japanese individuals, off-ladder calls were observed at the D12S391 locus using PowerPlex Fusion in nine samples for allele 22, one sample for allele 25, and one sample for allele 26. However, these samples were typed as ordinary alleles within the bins using GlobalFiler. In this study, next-generation sequencing analysis using MiSeq was performed for the D12S391 locus from the 11 off-ladder samples and 33 other samples, as well as the allelic ladders of PowerPlex Fusion and GlobalFiler. All off-ladder allele 22 in the nine samples had [AGAT]11[AGAC]11 as a repeat structure, while the corresponding allele was [AGAT]15[AGAC]6[AGAT] for the PowerPlex Fusion ladder, and [AGAT]13[AGAC]9 for the GlobalFiler ladder. Overall, as the number of [AGAT] in the repeat structure decreased at the D12S391 locus, the peak migrated more slowly using PowerPlex Fusion, the reverse strand of which was labeled, and it migrated more rapidly using GlobalFiler, the forward strand of which was labeled. The allelic ladders of both STR kits were reamplified with our small amplicon D12S391 primers and their mobility was also examined. In conclusion, off-ladder observations of allele 22 at the D12S391 locus using PowerPlex Fusion were mainly attributed to a relatively large difference of the repeat structure between its allelic ladder and off-ladder allele 22.

The common deletion found in patient reexamined after 33 years and comparison with complete mtDNA sequences of maternal relatives.

In 1966, a male (17 years old) was clinically examined at the National Institutes of Health (NIH) and diagnosed with Idiopathic Progressive External Ophthalmoplegia (IPEO). A muscle biopsy showing ragged-red fibers implicated mitochondrial involvement. Since the sequence of human mitochondrial DNA (mtDNA) was not determined until 1981, no genetic confirmation of the disease was possible at that time. In 1999, clinical reexamination and sequencing the entire mtDNA of the patient and living maternal relatives (mother and brother) indicated a progressive mitochondrial myopathy and the presence of the 4977 base pair (bp) deletion (the common deletion) in the patient.

DNA typing of bone specimens--the potential use of the profiler test as a tool for bone identification.

Twenty-six bone DNA identification cases are described. The postmortem periods of the studied remains ranged from three days to over 30 years, and the locations where the remains were found varied resulting in a variety of postmortem conditions. Nuclear DNA typing using an AmpFLSTR Profiler kit and mitochondrial DNA (mtDNA) typing of hypervariable regions 1 and 2 (HV1 and HV2) in a control region were performed both with decalcified and non-treated bone powder samples. Decalcification was shown to improve the success of DNA typing. The nucleotide sequences of the HV1 and HV2 regions were successfully determined in all cases examined. Nuclear DNA typing was very successful, more than half of the loci were typed during multiple amplifications (10 loci in one reaction) in 23 cases. Polymerase chain reaction (PCR) inhibition was observed in five cases including three samples that were found buried in soil. This inhibitory effect was identified as the result of unbalanced multiple PCR during the profiler test. These results revealed that DNA typing targeting nuclear DNA is a potentially powerful tool for bone identification.