A method for determining the origin of crude drugs derived from animals using MinION, a compact next-generation sequencer.

We evaluated whether MinION, an inexpensive portable sequencer, can be used to identify the origin of crude drugs derived from animals. Standard and nonstandard crude drugs with different species of origin were examined. In addition, standards mixed with nonstandard samples were used. As a target gene, cytochrome c oxidase I was amplified and sequenced. The fast mode results had a slightly lower match ratio than high-accuracy mode, but the animals of origin were correctly determined by BLAST for all samples. For antler velvet derived from Rangifer tarandus, even when the sequences were aligned based on Cervus elaphus, the animal of origin was determined correctly. Minor contents could be detected from mixtures of two animals, if the mixtures contained at least 19:1 mtDNA when the coverage allele-fraction threshold was 0.05. By contrast, in fast mode, two sequences could not be separated due to the low accuracy of the base-calling for each read. For fieldwork, the species of origin of crude drugs could be identified with only simple DNA extraction and library preparation. Therefore, MinION appears to be a convenient tool for identifying the origins of crude drugs derived from animals.

Development of an Okinawa panel for biogeographic inference of Okinawans.

BACKGROUND: The Precision ID Ancestry Panel with 165 SNP markers was unable to differentiate between mainland Japanese and Okinawa Japanese or to distinguish either of them from other East Asian populations. AIM: An Okinawa panel was developed with the aim of further separating Okinawa Japanese individuals from mainland Japanese and other Asian groups. Seventy-five SNPs were selected using the most informative markers from the literature. Further, 22 SNPs were selected to separate Okinawa Japanese at minimum SNPs. SUBJECTS AND METHODS: Samples were collected from 48 unrelated individuals from mainland Japan and 46 unrelated residents of the Okinawa prefecture. Data were evaluated by STRUCTURE, principal component, and GenoGeographer analyses. RESULTS: The 22 SNP set had similar levels of differentiation in STRUCTURE and PCA analyses as the 75 SNP set. GenoGeographer analysis showed that, out of the 46 Okinawa Japanese individuals, the 75 SNP and 22 SNP sets correctly assigned the Okinawan population as the most likely population of origin for 32 and 31 individuals, respectively. CONCLUSION: Neither SNP set could completely differentiate between Okinawa Japanese and other Asian groups, however, these sets should be useful for crime investigation, when the sample, cost and time are limited.

Estimating individual mtDNA haplotypes in mixed DNA samples by combining MinION and MiSeq.

We tried to estimate individual mtDNA haplotypes in mixed DNA samples by combining MinION and MiSeq. The BAM files produced by MiSeq were viewed using Integrative Genomics Viewer (IGV) to verify mixed bases. By sorting the reads according to base type for each mixed base, partial haplotypes were determined. Then, the BAM files produced by MinKNOW were viewed using IGV. To determine haplotypes with IGV, only mixed bases determined by MiSeq were used as target bases. By sorting the reads according to base type for each target base, each contributor's haplotype was estimated. In mixed samples from two contributors, even a haplotype with a minor contribution of 5% could be distinguished from the haplotype of the major contributor. In mixed samples of three contributors (mixture ratios of 1:1:1 and 4:2:1), each haplotype could also be distinguished. Sequences of C-stretches were determined very inaccurately in the MinION analysis. Although the analysis method was simple, each haplotype was correctly detected in all mixed samples with two or three contributors in various mixture ratios by combining MinION and MiSeq. This should be useful for identifying contributors to mixed samples.

Dissecting lesions in a culprit artery of a hemorrhagic focus in the basal ganglia: Histopathological analysis by serial sectioning.

In a hypertensive hemorrhagic focus of the basal ganglia, the culprit arteries have been reported to be associated with dissecting lesions, whose topographical relationship to the rupture sites remains to be clarified. Herein we describe multiple dissecting lesions in the culprit artery of hypertensive hemorrhage of the basal ganglia. A 1.0 × 0.8 × 0.8 cm-sized bleeding globe was confirmed at a left lenticulostriate artery and histologically analyzed by serial sectioning. Three independent dissecting lesions were identified in the culprit artery. They were situated near the bifurcations, ranging from 240 to 3200 µm in length. The dissections mainly occurred between the intima and media with disruption of the internal elastic lamina (IEL), forming a fresh thrombus within the false lumen. Two rupture sites causing the cerebral hematoma were confirmed away from the dissecting lesions. One was situated close but not adjacent to the longest dissecting lesion; the other, measuring approximately 150 µm in diameter, was adjacent to the bifurcation of an artery. The histopathological findings suggest that the dissecting lesion resulted from medial detachment following IEL disruption in the process of arterial rupture of the culprit artery. We conclude that this was a secondary manifestation during the rupture rather than a cause of the arterial rupture.

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.

Deaths Associated With Brotizolam Poisoning From a Single Drug Overdose: Four Reported Cases.

Brotizolam is a short-acting hypnotic in the benzodiazepine family, and fatal poisonings by an overdose of brotizolam are rare. This report describes 4 cases of deaths associated with brotizolam poisoning from a single drug overdose. The ages ranged from 51 to 90 years, and the postmortem interval between death and tissue sampling was 1.5 to 2.5 days. These deaths were classified as 1 homicide and 3 suicides. The concentration of the brotizolam ranged from 0.05 to 0.21 mg/L in the blood samples. Ethanol, which could cause mild alcohol intoxication, was detected in the blood samples from 2 cases. Postmortem examinations did not find any significant pathologic conditions, except for a case of death by drowning in a bathtub due to brotizolam poisoning. These 4 cases suggest that a brotizolam overdose should not be underestimated in terms of its fatal effects, particularly when situations involve alcohol intoxication, injury subsequent to the poisoning, or underlying medical conditions including aging.

Fatal Intracranial Aneurysms and Dissections Causing Subarachnoid Hemorrhage: An Epidemiological and Pathological Analysis of 607 Legal Autopsy Cases.

BACKGROUND: There are no detailed reports, in terms of epidemiology and pathology, on intracranial aneurysms and on dissections that were found in unexpected fatal subarachnoid hemorrhage (SAH) cases. In this report we analyzed, based on large-sized medicolegal autopsy cases, the detailed epidemiology and pathological aspects of both lesions. METHODS: We analyzed 607 autopsy cases of unexpected fatal SAHs including 496 aneurysms and 111 dissections. RESULTS: The following results were obtained: (1) Patients who died of dissections were younger than those who died of aneurysms; (2) symptom prevalence rates of aneurysms were 31.9%, appearing to be lower than those in previous studies; (3) a significantly higher prevalence of clinical symptoms was found in patients with dissections (60.5%) than patients with aneurysms; (4) hypertensive cardiomegaly was indicated in more than 80%, while no obvious difference in incidence in hypertensive cardiomegaly was noted between aneurysms and dissections; (5) aneurysms were found to occur much more frequently in the anterior communicating artery (31.9%) and vertebral arteries (7.5%), while dissections were found much more commonly in vertebral arteries (93.7%); and (6) the size of aneurysms was much smaller in general than that previously regarded as a risk factor of rupturing. CONCLUSIONS: These data might help in the prompt intervention in SAH and also in the prevention of lethal SAH in clinical settings.

Identification of canine saliva using mRNA-based assay.

A dog saliva analysis in addition to a bite-mark analysis may be important for evidence when a crime involves a dog bite. In this study, the utility of detecting canine saliva-specific mRNAs to identify canine saliva was evaluated. Canine saliva swabs (n = 20), urine swabs (n = 20), body surface swabs (n = 20), whole blood samples (n = 10), human saliva (n = 20), human skin surface swabs (n = 20), and human whole blood (n = 20) were tested. The saliva-specific genes encoding statherin (STATH), carbonic anhydrase VI (CA-VI), and dog allergens (Canf1 and Canf2) were analyzed as candidate genes. Moreover, glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was used as confirmation of canine mRNA extraction. STATH, CA-VI, Canf1, Canf2, and GAPDH mRNAs were detected in 19/20, 1/20, 11/20, 4/20, and 20/20 saliva samples, respectively. The STATH, CA-VI, Canf1, Canf2, and GAPDH mRNAs did not exhibit cross-reactivity with samples of human origin. This mRNA-based assay was also able to detect canine saliva in mock forensic samples. The results of this study indicated that the detection of STATH mRNA is useful for the identification of canine saliva, and GAPDH is a suitable marker for canine mRNA extraction.

A novel method for sex determination by detecting the number of X chromosomes.

A novel method for sex determination, based on the detection of the number of X chromosomes, was established. Current methods, based on the detection of the Y chromosome, can directly identify an unknown sample as male, but female gender is determined indirectly, by not detecting the Y chromosome. Thus, a direct determination of female gender is important because the quality (e.g., fragmentation and amelogenin-Y null allele) of the Y chromosome DNA may lead to a false result. Thus, we developed a novel sex determination method by analyzing the number of X chromosomes using a copy number variation (CNV) detection technique (the comparative Ct method). In this study, we designed a primer set using the amelogenin-X gene without the CNV region as the target to determine the X chromosome copy number, to exclude the influence of the CNV region from the comparative Ct value. The number of X chromosomes was determined statistically using the CopyCaller software with real-time PCR. All DNA samples from participants (20 males, 20 females) were evaluated correctly using this method with 1-ng template DNA. A minimum of 0.2-ng template DNA was found to be necessary for accurate sex determination with this method. When using ultraviolet-irradiated template DNA, as mock forensic samples, the sex of the samples could not be determined by short tandem repeat (STR) analysis but was correctly determined using our method. Thus, we successfully developed a method of sex determination based on the number of X chromosomes. Our novel method will be useful in forensic practice for sex determination.

Uncertainty in the number of contributor estimation methods applied to a Y-STR profile.

Maximum allele count (MAC) and total allele count (TAC) methods are widely used for estimating the number of contributors (NoC) of autosomal short tandem repeat (STR) profile in many forensic laboratories. In this study, we applied NoC estimation methods to mixed Y-STR profiles and evaluated its uncertainty and performance. For the MAC method, as recent Y-STR typing kits involve single- and multi-copy loci, we defined "MAC-single" for use across only single-copy loci and "MAC-multi" for use across only multi-copy loci. We generated a dataset containing 120,000 Y-STR profiles for a one to six-person mixture in silico based on previously reported haplotype frequencies of 27 Y-STR loci in Yfiler Plus for the U.S. population (reported by NIST) and the Henan Han population. The dataset was randomly split into a training set and a test set. The training set was used to construct a TAC distribution (TAC curve), whereas the test set was used to calculate the performance metrics (accuracy, precision, recall, and F1-score). In addition, the effect of the upper limit of NoC considered for estimation on overall accuracy was evaluated. The overall accuracies of MAC-single, MAC-multi, and TAC methods when the upper limit of NoC was set to six-person were 0.7920, 0.4329, and 0.7877 for the U.S. population and 0.8207, 0.4609, and 0.8385 for the Henan Han population. Our results suggest that the MAC-single and TAC methods can estimate the NoC for mixed Y-STR profiles with high levels of accuracy.

Estimating bloodstain age in the short term based on DNA fragment length using nanopore sequencer.

We used a nanopore sequencer to quantify DNA fragments > 10,000 bp in size and then evaluated their relationship with short-term bloodstain age. Moreover, DNA degradation was investigated after bloodstains were wetted once with water. Bloodstain samples on cotton gauze were stored at room temperature and low humidity for up to 6 months. Bloodstains stored for 1 day were wetted with nuclease-free water, allowed to dry, and stored at room temperature and low humidity for up to 1 week. The proportion of fragments > 20,000 bp in dry bloodstains tended to decrease over time, particularly for fragments > 50,000 bp in size. This trend was modeled using a power approximation curve, with the highest R2 value (0.6475) noted for fragments > 50,000 bp in size; lower values were recorded for shorter fragments. The proportion of longer fragments was significantly reduced in bloodstains that were dried after being wetted once, and there was significant difference in fragments > 50,000 bp between dry conditions and once-wetted. This result suggests that even temporary exposure to water causes significant DNA fragmentation, but not extensive degradation. Thus, bloodstains that appear fresh but have a low proportion of long DNA fragments may have been wetted previously. Our results indicate that evaluating the proportion of long DNA fragments yields information on both bloodstain age and the environment in which they were stored.

Effect existence of aging on stutter ratio evaluated via Bayesian inference.

The stochastic behavior of the stutter ratio (SR) in capillary electrophoresis-based DNA typing is currently described and predicted using statistical models in forensic genetics. Clarifying this behavior can help obtain more objective and robust evidence to the court in terms of mixture interpretation. This study aimed to investigate the effect existence of aging on SR via a Bayesian framework. Nail scrapings and clippings were collected from 68 healthy individuals with informed consent. Samples were classified by age-class: young group (0-16 years; n = 36) and older-adult group (>61 years; n = 32). Then, they were compared in terms of their SRs for each simple repeat locus included in GlobalFiler Kit. Bayesian modeling was performed with lognormal distribution model, which implemented multiple linear regression, allele and age-class as explanatory variables. For all simple repeat loci, the median of the posterior distribution of the age-class parameter was a positive value. For CSF1PO and D7S820, the 95% credible interval of the posterior distribution did not include 0. Our data suggested that aging slightly increases the SR. These findings might help elucidate the stochastic behavior of SR.

Simplified detection of the species of origin of antler velvets using single-stranded tag hybridization chromatographic printed-array strip.

In this study, we developed a convenient and easy-to-use origin identification method for antler velvets based on a simple DNA extraction technique and single-stranded tag hybridization chromatographic printed-array strip (STH-PAS). The primer sets used to detect Cervus elaphus, Rangifer tarandus, and 12S rRNA did not engage in non-specific reactions such as primer dimer formation. In both the triplex and singleplex assays, the sensitivity was < 1 ng DNA. Moreover, Cervus elaphus DNA could be detected in OTC crude drug products. Although the detection sensitivity resulting from the simplified extraction was slightly lower than that obtained with extraction by conventional methods, the amount of DNA was sufficient even from a small sample. The choice of a triplex or singleplex assay will depend on the purpose of the test. For example, if it is important to determine whether the antler velvet is derived from Cervus elaphus or Rangifer tarandus, a triplex assay is appropriate. If it is necessary to explore whether antler velvet from Cervus elaphus is included in an OTC crude drug product, a singleplex assay using the Cervus elaphus primer set is informative. If it is necessary to explore whether powdered antler velvet includes counterfeit products (from Rangifer tarandus), a singleplex assay employing the Rangifer tarandus primer is appropriate. The singleplex assay detects minor components even at a 1,000:1 ratio. Our study thus demonstrated the utility of a method combining simple DNA extraction with STH-PAS for efficient identification of the origin of antler velvets.

Autopsy findings of post-COVID-19 vaccination deaths in Tokyo Metropolis, Japan, 2021.

BACKGROUND: COVID-19 vaccines have been used across Japan since 17 February 2021, and as of 17 April 2022, 1690 deaths potentially caused by vaccine-related adverse effects have been reported to the Ministry of Health, Labour and Welfare. However, the causal relationship between vaccination and death could not be fully evaluated because of a lack of sufficient information. METHODS: Autopsy cases in which deaths occurred within seven days after COVID-19 vaccination in Tokyo Metropolis and were handled by medical examiners were selected (n = 54). Age, sex, vaccine-related information, cause of death, and possible causal relationship between vaccination and death were examined. RESULTS: The mean age of the deceased individuals was 68.1 years, and the study sample consisted of 34 males (63.9%) and 20 females (37.0%). Thirty-seven and six individuals received Comirnaty and Spikevax, respectively (68.5% and 11.1% respectively). The manner of death included natural (n = 43), non-natural (n = 8), and undetermined (n = 3). The most frequent cause of death was ischemic heart disease (n = 16). Regarding causal relationships, 46 cases (85.2%) did not show a causal relationship to vaccination, except for myocarditis (n = 3), thrombosis-related death (n = 4), and others (n = 1). CONCLUSION: Although many cases of deaths after COVID-19 vaccination in this study showed no definite causal relationship between the vaccination and deaths, some cases showed possible adverse events such as myocarditis. Autopsies are essential for detecting vaccine-related deaths, and the Japanese death investigation system needs to be reinforced from this viewpoint.

A sudden infant death associated with atrial septal aneurysm and abnormality of the connecting site between the inferior vena cava and right atrium.

An atrial septal aneurysm (ASA) is a rare cardiac anomaly characterized by varicose bulgingofthe atrial septum (oval fossa) into the left or right atrium. Pathogenesis and clinical significance of ASA are controversial. We report an autopsy case of a huge undiagnosed ASA with abnormality of the connecting site between the inferior vena cava and the right atrial ostium in a 2-month-oldJapanesefemale who died suddenly and unexpectedly. She was born at 36 weeks 4 days (body weight 3,110 g). No abnormality was detected during pregnancy or delivery. The postnatal growth was normal with no cardiac problem detected at the 1-month checkup. The ASA bulged off in a mass to the left atrium (width, 0.8 cm; excursion ratio, 53%), reaching close to the inflow site of the right pulmonary vein, with dilation of the pulmonary vein. The connecting site between the inferior vena cava and the right atrium was atypically located 1.6 cm away from the atrioventricular groove. Although most cases of ASA in an infant resolve physiologically as the infant grows, the infant in the present case is thought to have had an exceptional pathological ASA, possibly causing supraventricular arrhythmia. The abnormality of the connecting site between the inferior vena cava and the right atrium might have affected the development and continuation of the ASA.

Diagnostic significance of the histopathology of bone marrow macrophages in forensic autopsies.

Forensic pathologists often encounter autopsies that require an assessment of antemortem general conditions (e.g., infection, metabolic disorders). To establish evaluation clues for such cases, we quantitatively examined macrophages and the general pathology of bone marrow in samples from 180 forensic autopsy cases of decedents with various conditions. Hematoxylin-eosin staining, Berlin blue staining, and immunostainings for CD163, CD138, and CD61 were performed. We determined the numbers per field (density) of total macrophages, swollen macrophages, macrophages with hemophagocytosis, and hemosiderin-laden macrophages. Each density was standardized by identifying its ratio to the total number of macrophages. The decedents' background data (cause of death, other pathological findings, postmortem interval, antemortem symptoms, and presence of resuscitation) were extracted. No correlations were found between the postmortem interval and the other decedent data, indicating that these data are not affected by postmortem changes. In the group in which inflammatory disease was the cause of death, there were significant elevations in the ratio of the swollen macrophage density to total macrophages. Significantly higher ratios of the density of swollen and hemophagocytic macrophages were observed in the group in which conditions with a prolonged agonal period were the cause of death. The group with a return of spontaneous circulation to resuscitation showed a significantly higher ratio of macrophage density with hemophagocytosis. This study provides the first statistical analysis focused on bone marrow histopathology in forensic autopsies. The results will be useful for elucidating causes of death and agonal-period conditions.

Highly sensitive sex determination method using the exon 1 region of the amelogenin gene.

Sex determination is a crucial factor in the identification of unidentified human remains. Sex determination by DNA analysis is particularly useful because it can be applied to samples for which morphological characteristics are unavailable. Because samples handled in forensic DNA typing are easily degraded by environmental factors and microorganisms, there is a need for a method that can accurately determine sex even in highly decayed samples. Previous studies mainly used sex differences in an intron of the amelogenin gene. However, this region is highly polymorphic, and there are cases where accurate sexing cannot be performed because of genetic mutations in the target region. Thus, for reliable sex determination, it is desirable to select loci with as few non-sexual polymorphisms as possible. In this study, we focused on the exon 1 region of the amelogenin gene, which has very little polymorphism other than sex differences. We developed a primer set for sex determination and compared it with the GlobalFiler™ PCR Amplification Kit (GF), which is widely used for forensic DNA typing. The results showed that the amount of DNA required for accurate sex determination was 25 pg for both methods, achieving equivalent sensitivity. Next, we compared the two methods using ancient human skeletons and found that the present method with its shorter amplicon was considerably superior to GF. The present method is simple, rapid, inexpensive, and suitable for analyzing highly degraded samples. Therefore, this method is expected to contribute to forensic sciences and physical anthropology.

Unexpected death in a young child associated with anomalous aortic origin of the left main coronary artery without physical exertion: A case of an anomalous coronary artery with highly abundant elastic fibers.

Sudden death due to anomalous aortic origin of a coronary artery is far less common among young children in the absence of exercise stress. This report describes the case of a 2-year-old boy with a lower respiratory tract infection who suffered sudden cardiac arrest in his bed at home. The autopsy revealed that the left coronary artery (LCA) originated from the right sinus of Valsalva with an acute angle takeoff and traveled between the aorta and the pulmonary trunk (an interarterial course). Upon histological examination, the LCA, before reaching its major branches, was located adjacent to the outside of the aortic wall without an intramural passage, and the arterial wall was composed almost exclusively of elastic fibers without media containing smooth muscle cells throughout the entire length of the abnormal running. Screening tests for respiratory virus infection detected enterovirus in the lung tissue. In association with an acute angle takeoff and interarterial course, the wall structure with highly abundant elastic fibers that are more flexible tissues among blood vessel components might suggest their vulnerability to compression during the great vessels' systolic expansion in the sympathetic activation induced by the viral infection, leading to fatal myocardial ischemia without physical exertion.

Estimating included animal species in mixed crude drugs derived from animals using massively parallel sequencing.

We developed a method that can detect each animal species of origin for crude drugs derived from multiple animal species based on massively parallel sequencing analysis of mitochondrial genes. The crude drugs derived from animals investigated in this study were Cervi Parvum Cornu and Trogopterorum feces, which are derived from a mix of different animal species, two chopped cicada sloughs, and two commercial Kampo drugs. The mitochondrial 12S rRNA, 16S rRNA, and cytochrome oxidase subunit I gene regions were amplified and sequenced using MiSeq. The ratios of haplotype to total number of sequences reads were calculated after sequence extraction and trimming. Haplotypes that exceeded the threshold were defined as positive haplotypes, which were compared with all available sequences using BLAST. In the Cervi Parvum Cornu and Trogopterorum feces samples, the haplotype ratios corresponded roughly to the mixture ratios, although there was a slight difference from mixture ratios depending on the gene examined. This method could also roughly estimate the compositions of chopped cicada sloughs and Kampo drugs. This analysis, whereby the sequences of several genes are elucidated, is better for identifying the included animal species. This method should be useful for quality control of crude drugs and Kampo drugs.

Modeling the minus two base pair stutter ratio of the D1S1656 locus: A sequence-based mixture distribution model.

In this study, we propose a stutter ratio for a minus two base pair stutter (-2bpSR) model of the D1S1656 locus in capillary electrophoresis (CE)-based short tandem repeat (STR) typing. DNA from a total of 108 Japanese individuals was analyzed via massively parallel sequencing to investigate the length of the longest uninterrupted stretch of two base repeat motif (2bpLUS value) within repetitive structures involving the flanking region. Additionally, -2bpSR data was collected using the GlobalFiler Kit on a 3500xL Genetic Analyzer. As a result of sequencing analysis, all alleles were classified into two types by their 2bpLUS values. The -2bpSR differed significantly between the types. Then, we modeled the -2bpSR with a mixture log-normal distribution using the classification of alleles based on the 2bpLUS values. Furthermore, probabilities of the sequence type within each repeat number in the mixture log-normal distribution model were estimated using logistic regression for each of the five major detected populations. This study is expected to enable interpretation of STR typing while considering minus two base pair stutter at the D1S1656 locus.