Pivotal role of liver sinusoidal endothelial cells in NAFLD/NASH progression.

Liver sinusoidal endothelial cells (LSECs) are involved in the transport of nutrients, lipids, and lipoproteins, and LSEC injury occurs in various liver diseases including nonalcoholic fatty liver disease (NAFLD). However, the association between LSEC injury and NAFLD progression remains elusive. Accordingly, in this study, we aimed to elucidate the precise role of LSEC in the pathophysiology of NAFLD using two different mouse models, namely the choline-deficient, L-amino acid-defined and high-fat diet models. Administration of these diets resulted in liver metabolic dysregulation mimicking human NAFLD, such as steatosis, ballooning, lobular inflammation, and fibrosis, as well as central obesity, insulin resistance, and hyperlipidemia. LSEC injury appeared during the simple steatosis phase, and preceded the appearance of activated Kupffer cells and hepatic stellate cells (HSCs). These results indicate that LSEC injury may have a 'gatekeeper' role in the progression from simple steatosis to the early nonalcoholic steatohepatitis (NASH) stage, and LSEC injury may be necessary for the activation of Kupffer cells and HSCs, which in turn results in the development and perpetuation of chronic liver injuries. Taken together, our data provide new insights into the role of LSEC injury in NAFLD/NASH pathogenesis.

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.

Allele frequencies of 31 autosomal short tandem repeat (auSTR) loci obtained using the Precision ID GlobalFiler™ NGS STR Panel v2 in 322 individuals from the Japanese population.

In human identification methods that target short tandem repeats (STRs), massively parallel sequencing (MPS) technology has made it possible to genotype at the level of the specific sequence itself. This allows for the detection of repeat unit variants and single nucleotide polymorphisms (SNPs) adjacent to the STRs. Using the GlobalFiler™ NGS STR Panel v2, Ion S5, and Converge software, this study constructed a Japanese database of 31 autosomal STRs (auSTRs) and two sex markers from 322 individuals. After excluding some sequence errors and stutters, a total of 31 novel alleles were identified. Additionally, using the allele frequencies of 31 auSTR loci, the match probabilities for the length-based and sequence-based data were calculated to be 1.433 × 10-34 and 9.163 × 10-38, respectively. These values are at least nine orders of magnitude higher than that obtained from 21 auSTR loci in the Japanese population using the conventional capillary electrophoresis method. The database generated in this study is expected to be implemented in forensic practice and used to solve difficult casework.

Rapid semen identification from mixed body fluids using methylation-sensitive high-resolution melting analysis of the DACT1 gene.

In forensic genetics, a suspect is assigned to a component of a DNA mixture profile, and a probabilistic interpretation is then usually performed. However, it is difficult to determine what types of body fluid the component is from. Previous studies have reported that the fourth exon of the Dishevelled binding antagonist of beta catenin 1 (DACT1) gene is hypomethylated in a semen DNA-specific manner. In the present study, we evaluated whether the DACT1 gene could be effectively used to identify semen in body fluid mixtures and were able to semi-quantify the semen DNA content in mixed fluids. Our results showed that the DACT1 gene was useful in discriminating semen from venous blood and saliva. However, the amount of sperm in semen can affect semen identification. In addition, SI (the semen DNA content index), which we developed, was useful to determine whether the semen compromised majority, almost half, or was in the minority of the components in a mixed fluid. This technique is based on the methylation-sensitive high-resolution melting (MS-HRM) technology, which is time-, cost-, and labour-effective, and could be adopted in routine criminal investigations.

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.

Development of a software for kinship analysis considering linkage and mutation based on a Bayesian network.

In forensic DNA testing, the number of tested short tandem repeat loci has increased owing to new multiplex kits with additional loci. Although this advancement provides improved discrimination power, the effects of linkage and mutation must be considered during kinship analysis. However, no software currently includes both of these effects. In this study, we developed new freeware called KinBN for kinship analysis based on a Bayesian network. The software is graphical-user-interface-based and calculates the likelihood ratios (LRs) at multiple loci considering the effects of linkage and mutation. In addition, the software can simulate the LR distribution according to the specified relationship. We confirmed the accuracy of KinBN by comparing its LRs with those of other software and evaluated the effects of linkage and mutation on the LRs. Our results indicate that KinBN is a useful tool for kinship analysis, particularly if expanded locus sets are used for DNA testing.

Distinct spectrum of microRNA expression in forensically relevant body fluids and probabilistic discriminant approach.

MicroRNA is attracting worldwide attention as a new marker for the identification of forensically relevant body fluids. A probabilistic discriminant model was constructed to identify venous blood, saliva, semen, and vaginal secretion, based on microRNA expression assessed via RT-qPCR. We quantified 15 candidate microRNAs in four types of body fluids by RT-qPCR and found that miR-144-3p, miR-451a-5p, miR-888-5p, miR-891a-5p, miR-203a-3p, miR-223-3p and miR-1260b were helpful to discriminate body fluids. Using the relative expression of seven candidate microRNAs in each body fluid, we implemented a partial least squares-discriminant analysis (PLS-DA) as a probabilistic discriminant model and distinguished four types of body fluids. Of 14 testing samples, 13 samples were correctly identified with >90% posterior probability. We also investigated the effects of microRNA expression in skin, semen infertility, and vaginal secretion during different menstrual phases. Semen infertility and menstrual phases did not affect our body fluid identification system. Therefore, the selected microRNAs were effective in identifying the four types of body fluids, indicating that probabilistic evaluation may be practical in forensic casework.

Optimal small-molecular reference RNA for RT-qPCR-based body fluid identification.

MicroRNA (miRNA) -based body fluid identification (BFID) plays a prominent role in a forensic practice, and the selected reference RNA is indispensable for a robust normalization in BFID performed using reverse transcription-quantitative PCR. In this study, we first examined sample quality using RNA integrity number, then evaluated the consistency of expression of candidate reference RNAs in 4 forensically relevant body fluids using NormFinder and BestKeeper, and lastly used each rank and index output from these tools for selecting the optimal reference RNA and the combination of the multiple RNAs using the RankAggreg package of R. We found that RNA integrity number was small in our samples, despite the use of pristine body fluids; 5S-rRNA was the optimal reference RNA for the identification of forensically relevant body fluids; and the combination of 5S-rRNA and miR-92a-3p and/or miR-484 enhanced the normalization quality. Our findings enable us to perform stringent normalization of the expression of body fluid-specific RNAs, and thus, can contribute to the development of small RNA-based BFID systems.

Discrimination of relationships with the same degree of kinship using chromosomal sharing patterns estimated from high-density SNPs.

Distinguishing relationships with the same degree of kinship (e.g., uncle-nephew and grandfather-grandson) is generally difficult in forensic genetics by using the commonly employed short tandem repeat loci. In this study, we developed a new method for discerning such relationships between two individuals by examining the number of chromosomal shared segments estimated from high-density single nucleotide polymorphisms (SNPs). We computationally generated second-degree kinships (i.e., uncle-nephew and grandfather-grandson) and third-degree kinships (i.e., first cousins and great-grandfather-great-grandson) for 174,254 autosomal SNPs considering the effect of linkage disequilibrium and recombination for each SNP. We investigated shared chromosomal segments between two individuals that were estimated based on identity by state regions. We then counted the number of segments in each pair. Based on our results, the number of shared chromosomal segments in collateral relationships was larger than that in lineal relationships with both the second-degree and third-degree kinships. This was probably caused by differences involving chromosomal transitions and recombination between relationships. As we probabilistically evaluated the relationships between simulated pairs based on the number of shared segments using logistic regression, we could determine accurate relationships in >90% of second-degree relatives and >70% of third-degree relatives, using a probability criterion for the relationship ≥0.9. Furthermore, we could judge the true relationships of actual sample pairs from volunteers, as well as simulated data. Therefore, this method can be useful for discerning relationships between two individuals with the same degree of kinship.

Forensic age prediction for saliva samples using methylation-sensitive high resolution melting: exploratory application for cigarette butts.

There is high demand for forensic age prediction in actual crime investigations. In this study, a novel age prediction model for saliva samples using methylation-sensitive high resolution melting (MS-HRM) was developed. The methylation profiles of ELOVL2 and EDARADD showed high correlations with age and were used to predict age with support vector regression. ELOVL2 was first reported as an age predictive marker for saliva samples. The prediction model showed high accuracy with a mean absolute deviation (MAD) from chronological age of 5.96 years among 197 training samples. The model was further validated with an additional 50 test samples (MAD = 6.25). In addition, the age prediction model was applied to saliva extracted from seven cigarette butts, as in an actual crime scene. The MAD (7.65 years) for these samples was slightly higher than that of intact saliva samples. A smoking habit or the ingredients of cigarettes themselves did not significantly affect the prediction model and could be ignored. MS-HRM provides a quick (2 hours) and cost-effective (95% decreased compared to that of DNA chips) method of analysis. Thus, this study may provide a novel strategy for predicting the age of a person of interest in actual crime scene investigations.

Development and validation of open-source software for DNA mixture interpretation based on a quantitative continuous model.

In criminal investigations, forensic scientists need to evaluate DNA mixtures. The estimation of the number of contributors and evaluation of the contribution of a person of interest (POI) from these samples are challenging. In this study, we developed a new open-source software "Kongoh" for interpreting DNA mixture based on a quantitative continuous model. The model uses quantitative information of peak heights in the DNA profile and considers the effect of artifacts and allelic drop-out. By using this software, the likelihoods of 1-4 persons' contributions are calculated, and the most optimal number of contributors is automatically determined; this differs from other open-source software. Therefore, we can eliminate the need to manually determine the number of contributors before the analysis. Kongoh also considers allele- or locus-specific effects of biological parameters based on the experimental data. We then validated Kongoh by calculating the likelihood ratio (LR) of a POI's contribution in true contributors and non-contributors by using 2-4 person mixtures analyzed through a 15 short tandem repeat typing system. Most LR values obtained from Kongoh during true-contributor testing strongly supported the POI's contribution even for small amounts or degraded DNA samples. Kongoh correctly rejected a false hypothesis in the non-contributor testing, generated reproducible LR values, and demonstrated higher accuracy of the estimated number of contributors than another software based on the quantitative continuous model. Therefore, Kongoh is useful in accurately interpreting DNA evidence like mixtures and small amounts or degraded DNA samples.

New stutter ratio distribution for DNA mixture interpretation based on a continuous model.

In forensic science, DNA mixture interpretation is traditionally based on a binary model, which does not account for peak-height information in DNA profiles. In recent years, some countries have adopted a continuous model in which peak heights are used and stochastic effects are considered to enable rigorous calculation of likelihood ratios. However, this model requires certain biological parameters which affect the expected allelic and stutter peak heights. In this paper, we focused on estimating the distribution of the stutter ratio (SR) in 15 short tandem repeat loci in relation to the allele repeat number. We estimated the SR values of 234 single-source DNA samples by using a commercially available kit. In all loci except for D8S1179, D21S11, and D2S1338, a simple log-normal distribution model was fitted to the variability of SR. For D21S11, we developed a new distribution model in which distinct log-normal distributions between complete and incomplete repeat units are used (a separate log-normal distribution model). For D8S1179 and D2S1338, we developed another new distribution model that mixes two log-normal distributions to explain two types of repeat structures appearing within the same number of allele repeats. These two models were fitted to the observed SR values more accurately than the simple log-normal distribution model. We expected these new SR models to be applied to DNA mixture interpretation based on a continuous model.

Forensic age prediction for dead or living samples by use of methylation-sensitive high resolution melting.

Age prediction with epigenetic information is now edging closer to practical use in forensic community. Many age-related CpG (AR-CpG) sites have proven useful in predicting age in pyrosequencing or DNA chip analyses. In this study, a wide range methylation status in the ELOVL2 and FHL2 promoter regions were detected with methylation-sensitive high resolution melting (MS-HRM) in a labor-, time-, and cost-effective manner. Non-linear-distributions of methylation status and chronological age were newly fitted to the logistic curve. Notably, these distributions were revealed to be similar in 22 living blood samples and 52 dead blood samples. Therefore, the difference of methylation status between living and dead samples suggested to be ignorable by MS-HRM. Additionally, the information from ELOVL2 and FHL2 were integrated into a logistic curve fitting model to develop a final predictive model through the multivariate linear regression of logit-linked methylation rates and chronological age with adjusted R(2)=0.83. Mean absolute deviation (MAD) was 7.44 for 74 training set and 7.71 for 30 additional independent test set, indicating that the final predicting model is accurate. This suggests that our MS-HRM-based method has great potential in predicting actual forensic age.

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.

Pairwise Kinship Analysis by the Index of Chromosome Sharing Using High-Density Single Nucleotide Polymorphisms.

We developed a new approach for pairwise kinship analysis in forensic genetics based on chromosomal sharing between two individuals. Here, we defined "index of chromosome sharing" (ICS) calculated using 174,254 single nucleotide polymorphism (SNP) loci typed by SNP microarray and genetic length of the shared segments from the genotypes of two individuals. To investigate the expected ICS distributions from first- to fifth-degree relatives and unrelated pairs, we used computationally generated genotypes to consider the effect of linkage disequilibrium and recombination. The distributions were used for probabilistic evaluation of the pairwise kinship analysis, such as likelihood ratio (LR) or posterior probability, without allele frequencies and haplotype frequencies. Using our method, all actual sample pairs from volunteers showed significantly high LR values (i.e., ≥ 108); therefore, we can distinguish distant relationships (up to the fifth-degree) from unrelated pairs based on LR. Moreover, we can determine accurate degrees of kinship in up to third-degree relationships with a probability of > 80% using the criterion of posterior probability ≥ 0.90, even if the kinship of the pair is totally unpredictable. This approach greatly improves pairwise kinship analysis of distant relationships, specifically in cases involving identification of disaster victims or missing persons.

Ectopic cervical thymus: a clinicopathological study of consecutive, unselected infant autopsies.

OBJECTIVES: An ectopic cervical thymus (ECT) is regarded as a rare congenital anomaly; therefore, the optimal diagnostic and therapeutic strategy remains a debatable matter. We designed a study to elucidate the clinicopathological characteristics of ECTs in consecutive, unselected infant autopsies, to help guide case management. METHODS: We searched for ECTs in all of the 21 consecutive, unselected infant autopsy cases performed at our institution over a period of 3 years, and all ECT consensus diagnoses were confirmed by histological examination. The following clinical characteristics were evaluated in cases with ECTs: age, gender, birth week and weight, clinical symptoms due to the ECT(s), position on discovery of death, cause of death, ECT contribution to the cause of death, and concomitant congenital disorders. The anatomical features evaluated included the location, number, size, color, shape, and margins of the ECTs, and the presence of a mediastinal thymus. Histological findings of the ECT(s) and the mediastinal thymus were compared within each individual. Fusion of the parathyroid and the ECT was also investigated histologically. Spearman's rank correlation coefficient (ρ) and the corresponding P value were calculated to determine if there was an association between ECT diameter and age. RESULTS: We detected 10 ECT lesions in seven cases (33%) among the 21 infant autopsy cases. The ECT cases involved five boys and two girls, with ages ranging from 1 day to 4 months. There were no reports of a positive family history of sudden death or antemortem clinical symptoms due to ECT in any of the cases. The ECTs were considered incidental regarding the cause of death, with the exception of one case that was extremely rare. Most ECTs were localized to the inferior thyroid, ranging from 0.4 to 1.9 cm in size. Size demonstrated a significant negative correlation with age (ρ=-0.75 and P=0.034). CONCLUSIONS: This study revealed that ECT is an essentially benign anomaly that occurs frequently during the development of the thymus, and may disappear over the first few years of life. These results suggest a conservative approach to the management of ECTs would be appropriate.

Mixture interpretation: Experimental and simulated reevaluation of qualitative analysis.

We present here analytical data using the 15 STR typing (Identifiler) kit regarding heterozygote balance in experimental DNA samples including one or two persons. Surprisingly, the allelic imbalance was observed even in samples consisting of only one person but adequate DNA for the standard protocol. The variance of heterozygote balance was more expanded in two-person mixtures than in one-person samples. Therefore, it is not suitable to use allelic peak heights/areas for estimating the genotypes of the contributors such as the quantitative analysis. We also reevaluated the effectiveness of qualitative analysis by simulation, i.e. consideration of the probability of all possible genotype combinations from the typing results of a mixed DNA sample. As demonstrated, the qualitative analysis using 15 STR loci is still extremely effective even in a mixture from two or three individuals.

Bile canalicular abnormalities in the early phase of a mouse model of sclerosing cholangitis.

BACKGROUND: The bile canaliculus is the smallest and first biliary channel and is formed by two or three adjacent hepatocytes. Previous studies of chronic cholangiopathies such as primary sclerosing cholangitis have focused on the bile ductules. However, little is known about the pathological alterations in bile canaliculi in the early phase of cholangiopathies. AIM: To characterize the bile canalicular morphology in the early phase of sclerosing cholangitis we used 3,5-diethoxycarbonyl-1,4-dihydrocollidine-induced mouse model of sclerosing cholangitis. METHODS: Mice were fed a diet with 3,5-diethoxycarbonyl-1,4-dihydrocollidine (0.1%). Serum biochemical, histological, immunohistochemical, and electron microscopic analyses were performed 1, 2, 4, and 7 days after feeding. RESULTS: All experimental groups showed significantly increased serum aspartate aminotransferase, alanine aminotransferase, and alkaline phosphatase levels. From day 1, bile canalicular abnormalities such as dilatation and meandering and loss of microvilli were observed. After bile canalicular abnormalities had appeared, substantial infiltration of inflammatory cells was observed amongst the necrotic cells and periductal region. After these inflammatory changes, cholangiocytes proliferated in the portal area and formed ductular reactions. Finally, periductal fibrosis appeared. CONCLUSION: This study provides novel evidence of the occurrence of bile canalicular abnormalities during the early phase of sclerosing cholangitis.