Machine learning approaches to explore digenic inheritance.

Some rare genetic disorders, such as retinitis pigmentosa or Alport syndrome, are caused by the co-inheritance of DNA variants at two different genetic loci (digenic inheritance). To capture the effects of these disease-causing variants and their possible interactive effects, various statistical methods have been developed in human genetics. Analogous developments have taken place in the field of machine learning, particularly for the field that is now called Big Data. In the past, these two areas have grown independently and have started to converge only in recent years. We discuss an overview of each of the two fields, paying special attention to machine learning methods for uncovering the combined effects of pairs of variants on human disease.

Strategic validation of variants of uncertain significance in ECHS1 genetic testing.

BACKGROUND: Enoyl-CoA hydratase short-chain 1 (ECHS1) is an enzyme involved in the metabolism of branched chain amino acids and fatty acids. Mutations in the ECHS1 gene lead to mitochondrial short-chain enoyl-CoA hydratase 1 deficiency, resulting in the accumulation of intermediates of valine. This is one of the most common causative genes in mitochondrial diseases. While genetic analysis studies have diagnosed numerous cases with ECHS1 variants, the increasing number of variants of uncertain significance (VUS) in genetic diagnosis is a major problem. METHODS: Here, we constructed an assay system to verify VUS function for ECHS1 gene. A high-throughput assay using ECHS1 knockout cells was performed to index these phenotypes by expressing cDNAs containing VUS. In parallel with the VUS validation system, a genetic analysis of samples from patients with mitochondrial disease was performed. The effect on gene expression in cases was verified by RNA-seq and proteome analysis. RESULTS: The functional validation of VUS identified novel variants causing loss of ECHS1 function. The VUS validation system also revealed the effect of the VUS in the compound heterozygous state and provided a new methodology for variant interpretation. Moreover, we performed multiomics analysis and identified a synonymous substitution p.P163= that results in splicing abnormality. The multiomics analysis complemented the diagnosis of some cases that could not be diagnosed by the VUS validation system. CONCLUSIONS: In summary, this study uncovered new ECHS1 cases based on VUS validation and omics analysis; these analyses are applicable to the functional evaluation of other genes associated with mitochondrial disease.

Impact of measuring heteroplasmy of a pathogenic mitochondrial DNA variant at the single-cell level in individuals with mitochondrial disease.

Pathogenic mitochondrial DNA heteroplasmy has mainly been assessed with bulk sequencing in individuals with mitochondrial disease. However, the distribution of heteroplasmy at the single-cell level in skin fibroblasts obtained from individuals, together with detailed clinical and biochemical information, remains to be investigated. We used the mitochondrial DNA single-cell assay for the transposase-accessible chromatin sequencing method. Skin fibroblasts were obtained from six individuals with mitochondrial disease and pathogenic m.3243A>G variants of differing severity. Different distributions of heteroplasmy at the single-cell level were identified in skin fibroblasts from all six individuals. Four individuals with different outcomes showed similar averaged heteroplasmy rates with normal mitochondrial respiratory chain enzyme activity, while the distribution of single-cell heteroplasmy patterns differed among the individuals. This study showed different heteroplasmy distribution patterns at the single-cell level in individuals with the m.3243A>G variant, who had a similar averaged heteroplasmy rates with normal mitochondrial respiratory chain enzyme activity. Whether such different heteroplasmy distribution patterns explain the different clinical outcomes should be assessed further in future studies. Measuring heteroplasmy of pathogenic mitochondrial DNA variants at the single-cell level could be important in individuals with mitochondrial disease.

Population genetics: past, present, and future.

We present selected topics of population genetics and molecular phylogeny. As several excellent review articles have been published and generally focus on European and American scientists, here, we emphasize contributions by Japanese researchers. Our review may also be seen as a belated 50-year celebration of Motoo Kimura's early seminal paper on the molecular clock, published in 1968.

Case Report of Left Ventricular Noncompaction Cardiomyopathy Characterized by Undulating Phenotypes in Adult Patients.

Left ventricular noncompaction cardiomyopathy (LVNC) is a heart muscle disorder morphologically characterized by reticulated trabeculations and intertrabecular recesses in the left ventricular (LV) cavity. LVNC is a genetically and phenotypically heterogeneous condition, which has been increasingly recognized with the accumulation of evidence provided by genotype-phenotype correlation analyses. Here, we report 2 sporadic adult cases of LVNC; both developed acute heart failure as an initial clinical manifestation and harbored causal sarcomere gene mutations. One case was a 57-year-old male with digenic heterozygote mutations, p.R1344Q in myosin heavy chain 7 (MYH7) and p.R144W in troponin T2, cardiac type (TNNT2), who showed morphological characteristics of LVNC in the lateral to apical regions of the LV together with a comorbidity of non-transmural myocardial infarction, resulting from a coronary artery stenosis. After the removal of ischemic insult and standard heart failure treatment, LVNC became less clear, and LV function gradually improved. The other case was a 36-year-old male with a heterozygote mutation, p.E334K in myosin binding protein C3 (MYBPC3), who exhibited cardiogenic shock on admission with morphological characteristics of LVNC being most prominent in the apical segment of the LV. The dosage of beta-blocker was deliberately increased in an outpatient clinic over 6 months following hospitalization, which remarkably improved the LV ejection fraction from 21% to 54.3%. Via a combination of imaging and histopathological and genetic tests, we have found that these cases are not compatible with a persistent phenotype of primary cardiomyopathy, but their morphological features are changeable in response to treatment. Thus, we point out phenotypic plasticity or undulation as a noticeable element of LVNC in this case report.

Mortality of Japanese patients with Leigh syndrome: Effects of age at onset and genetic diagnosis.

Leigh syndrome is a major phenotype of mitochondrial diseases in children. With new therapeutic options being proposed, assessing the mortality and clinical condition of Leigh syndrome patients is crucial for evaluating therapeutics. As data are scarce in Japan, we analysed the mortality rate and clinical condition of Japanese Leigh syndrome patients that we diagnosed since 2007. Data from 166 Japanese patients diagnosed with Leigh syndrome from 2007 to 2017 were reviewed. Patients' present status, method of ventilation and feeding, and degree of disability as of April 2018 was analysed. Overall, 124 (74.7%) were living, 40 (24.1%) were deceased, and 2 (1.2%) were lost to follow-up. Median age of living patients was 8 years (1-39 years). Median length of disease course was 91 months for living patients and 23.5 months for deceased patients. Nearly 90% of deaths occurred by age 6. Mortality rate of patients with onset before 6 months of age was significantly higher than that of onset after 6 months. All patients with neonatal onset were either deceased or bedridden. MT-ATP6 deficiency caused by m.8993T>G mutation and MT-ND5 deficiency induced a severe form of Leigh syndrome. Patients with NDUFAF6, ECHS1, and SURF1 deficiency had relatively mild symptoms and better survival. The impact of onset age on prognosis varied across the genetic diagnoses. The clinical condition of many patients was poor; however, few did not require mechanical ventilation or tube-feeding and were not physically dependent. Early disease onset and genetic diagnosis may have prognostic value.

Heterozygosity mapping for human dominant trait variants.

Homozygosity mapping is a well-known technique to identify runs of homozygous variants that are likely to harbor genes responsible for autosomal recessive disease, but a comparable method for autosomal dominant traits has been lacking. We developed an approach to map dominant disease genes based on heterozygosity frequencies of sequence variants in the immediate vicinity of a dominant trait. We demonstrate through theoretical analysis that DNA variants surrounding an inherited dominant disease variant tend to have increased heterozygosity compared with variants elsewhere in the genome. We confirm existence of this phenomenon in sequence data with known dominant pathogenic variants obtained on family members and in unrelated population controls. A computer-based approach to estimating empirical significance levels associated with our test statistics shows genome-wide p-values smaller than 0.05 for many but not all of the individuals carrying a pathogenic variant.

Barth Syndrome: Different Approaches to Diagnosis.

The diagnosis of Barth syndrome is challenging owing to the wide phenotypic spectrum with allelic heterogeneity. Here we report 3 cases of Barth syndrome with phenotypic and allelic heterogeneity that were diagnosed by different approaches, including whole exome sequencing and final confirmation by reverse-transcription polymease chain reaction.

HDR-del: A tool based on Hamming distance for prioritizing pathogenic chromosomal deletions in exome sequencing.

High-density oligonucleotide arrays have widely been used to detect pathogenic chromosomal deletions. In addition to high-density oligonucleotide arrays, programs using whole-exome sequencing have become available for estimating copy-number variations using depth of coverage. Here, we propose a new statistical method, HDR-del, to prioritize pathogenic chromosomal deletions based on Hamming distance in exome sequencing. In vcf (variant call format) files generated from exome sequencing, hemizygous chromosomal deletion regions lack heterozygous variants and lead to apparent long runs of homozygosity (ROH). In our Hamming distance ratio (HDR)-del approach, we calculate the "difference" in heterozygous status between an affected individual and control individuals using the HDR over all candidate chromosomal deletion regions defined as ROH longer than 1Mbp. Using a suitable test statistic, which is expected to be large for a true pathogenic deletion region, we prioritize candidate chromosomal deletion regions based on this statistic. In our approach, we were able to considerably narrow down true pathogenic chromosomal deletion regions, which were confirmed by high-density oligonucleotide arrays in four mitochondrial disease patients. Our HDR-del approach represents an easy method for detecting chromosomal deletions.

Outcomes of labor epidural analgesia among women aged over 40: A single-institution retrospective study.

AIM: The purpose of this study was to investigate the effects of labor epidural analgesia (LEA) on maternal and neonatal outcomes among parturients aged 40 years or older. METHODS: We retrospectively reviewed medical records of all laboring, singleton and cephalic deliveries at ≥36 weeks' gestation at the Saitama Medical Center from April 2003 to September 2012. Women aged ≥40 years who received LEA (≥40 with LEA group) were compared with women aged ≥40 years who delivered without LEA (≥40 without LEA group) and women <40 years who received LEA (<40 with LEA group). Extracted outcomes included mode of delivery, oxytocin augmentation, duration of labor, amount of estimated blood loss, umbilical artery pH, Apgar scores, and neonatal intensive care unit admission. RESULTS: This study included 4441 women. There were 74 women in the ≥40 with LEA group, 369 in the ≥40 without LEA group, and 601 in the <40 with LEA group. The maternal outcomes of emergency cesarean delivery rate (9.5%, 12.5%, 9.0%), instrumental delivery rate (33.8%, 10.3%, 28.3%), duration of labor (521 min, 321 min, 565 min), and estimated blood loss (524 g, 351 g, 412 g) were reported for the ≥40 with LEA, ≥40 without LEA, and <40 with LEA groups, respectively. Neonatal outcomes were not different between these groups. LEA use was not associated with emergency cesarean delivery in the multivariable analysis. CONCLUSION: Our study showed that parturients aged ≥40 with LEA can expect similar LEA-associated labor outcomes to younger parturients with LEA.

Fetal heart rate pattern interpretation in the second stage of labor using the five-tier classification: impact of the degree and duration on severe fetal acidosis.

AIM: The aim of this study was to clarify the association between fetal heart rate (FHR) tracing interpretation levels in the second stage of labor and poor fetal acid-base balance. MATERIAL AND METHODS: The database at one tertiary hospital in Nagoya, Japan, was retrospectively reviewed for women with singleton fetuses in cephalic presentation and vaginal labor at ≥37 + 0 gestational weeks between 1 June 2011 and 30 April 2012. Continuous FHR tracings in the second stage of labor were subdivided into 15-min intervals, each of which we called a window, from the beginning of labor through delivery, and were assessed according to the five-tier classification proposed by the Japan Society of Obstetrics and Gynecology, in which level 1 is normal, level 2 is subnormal, and levels 3-5 are abnormal patterns. RESULTS: In total, 777 parturient women were eligible for the study protocol. The numbers of women with maximal levels of 1, 2, 3, 4, and 5 were 3, 77, 341, 349, and 7, respectively. No cases of severe fetal acidosis (pH < 7.0 or base excess <-12 mmol/L) were recorded when the maximal levels were below 3. Both the pH and base excess of the umbilical artery decreased with higher levels of FHR tracings interpretation (P < 0.001). Both the summations of level-4 windows and level-3 and level-4 windows were significantly higher in women with severe fetal acidosis than in women without (P < 0.001), indicating that the duration of abnormal levels is associated with severe fetal acidosis. CONCLUSIONS: Both the degree and duration of FHR tracing abnormalities correlate with severe fetal acidosis.

Microstripe pattern substrate consisting of alternating planar and nanoprotrusive regions improved hiPSC-derived cardiomyocytes' unidirectional alignment and functional properties.

The alignment of each cell in human myocardium is considered critical for the efficient movement of cardiac tissue. We investigated 96-well microstripe-patterned plates to align human induced pluripotent stem cell (hiPSC)-derived cardiomyocytes (CMs), which resemble fetal myocardium. The aligned CMs (ACMs) cultured on the microstripe-patterned plates exhibited pathology, motor function, gene expression, and drug response that more closely resembled those of adult cells than did unaligned CMs cultured on a flat plate (FCMs). We used these ACMs to evaluate drug side effects and efficacy, and to determine whether these were similar to adult-like responses. When CMs from patients with hypertrophic cardiomyopathy (HCMs) were seeded and cultured on the microstripe-patterned plates or layered on top of the ACMs, both sets of HCMs showed increased heart rate and synchronized contractions, indicating improved cardiac function. It is suggested that the ACMs could be used for drug screening as cells representative of adult-like CMs and be transplanted in the form of a cell sheet for regenerative treatment of heart failure.

Genetic Landscape of Masticatory Muscle Tendon-Aponeurosis Hyperplasia.

Limited mouth opening is a characteristic of masticatory muscle tendon-aponeurosis hyperplasia (MMTAH). Although genetic involvement is suspected where familial onset is frequently observed, the genetic background of MMTAH is yet to be elucidated. In this study, we conducted whole genome sequencing of 10 patients with MMTAH and their family members when available. We also conducted RNA sequencing of normal temporal tendon (as disease region) and Achilles tendon (as control region) from commercially available pig samples. We identified 51 genes that had rare variants in patients with MMTAH and were highly expressed in the temporal tendons of pigs. Among the 51 genes, 37 genes have not been reported to be causative for human genetic diseases so far. As an implication of genetic involvement in the pathogenesis of MMTAH, 21 of these 37 genes were identified in two independent families. In particular, PCDH1 and BAIAP3 were identified in one affected individual in a family and consistently segregated in unrelated family, indicating they could be candidate causative genes of MMTAH. Our findings will help elucidate the genetic landscape of MMTAH and provide insights into future possibilities for tendon regeneration treatment.

Association between maternally inherited deafness, epilepsy, and intellectual disability and the m.12207G > A MT-TS2 pathogenic variant in a Japanese family.

The identification of the m.12207G > A variant in MT-TS2, (NC_012920.1:m.12207G > A) was first reported in 2006. The affected individual presented with developmental delay, feeding difficulty, proximal muscle weakness, and lesions within her basal ganglia, with heteroplasmy levels of 92% in muscle and no evidence of maternal inheritance. Herein, we report a case involving a 16-year-old boy with the same pathogenic variation and different phenotype, including sensorineural deafness, epilepsy, and intellectual disability, without diabetes mellitus (DM). His mother and maternal grandmother had similar but milder symptoms with DM. Heteroplasmy levels of the proband in blood, saliva, and urinary sediments were 31.3%, 52.6%, and 73.9%, respectively, while those of his mother were 13.8%, 22.1%, and 29.4%, respectively. The differences in the symptoms might be explained by the different levels of heteroplasmy. To our knowledge, this is the first familial report of the m.12207G > A variant in MT-TS2 that causes DM. The present case showed milder neurological symptoms than did the former report, and suggests the presence of a good phenotype-genotype correlation within this family.

Mitochondrial DNA depletion syndrome with a mutation in SLC25A4 developing epileptic encephalopathy: A case report.

INTRODUCTION: Autosomal dominant mitochondrial DNA depletion syndrome (MTDPS-12A) is characterized by severe hypotonia from birth due to a mutation in the adenine nucleotide translocator 1 (ANT1). CASE REPORT: A 4-year-old female patient diagnosed with neonatal-onset mitochondrial disease, who had good cognitive function while receiving antiepileptic treatment, presented with sudden-onset status epilepticus with facial and limb myoclonus persisting for more than 30 min. Subsequently, she developed epileptic encephalopathy. Brain MRI showed progressive ventricular enlargement and marked white matter atrophy. She was unable to perform verbal communication or make eye contact and fingertip movements. She lacked any signs of cardiomyopathy. Sanger sequencing demonstrated a heterozygous de novo mutation of c.239G>A (p.Arg80His) in SLC25A4. Her right quadriceps muscle tissue showed lowered complexes I, III, and IV activities and mitochondria DNA depletion (mitochondria/nuclear DNA: 14.6 ± 2.2%) through the quantitative polymerase chain reaction. She was definitively diagnosed with MTDPS-12A. CONCLUSION: Status epilepticus causes encephalopathy in patients with MTDPS-12A. Reducing the energy requirement on the cardiac muscle and brain may be a treatment strategy for patients with MTDPS-12A. Therefore, seizure management and preventive treatment of status epilepticus are considered to be important for maintaining neurodevelopmental outcomes.

A heterozygous LAMA5 variant may contribute to slowly progressive, vinculin-enhanced familial FSGS and pulmonary defects.

The LAMA5 gene encodes laminin α5, an indispensable component of glomerular basement membrane and other types of basement membrane. A homozygous pathological variant in LAMA5 is known to cause a systemic developmental syndrome including glomerulopathy. However, the roles of heterozygous LAMA5 gene variants in human renal and systemic diseases have remained unclear. We performed whole-exome sequencing analyses of a family with slowly progressive nephropathy associated with hereditary focal segmental glomerulosclerosis, and we identified what we believe to be a novel probable pathogenic variant of LAMA5, NP_005551.3:p.Val3687Met. In vitro analyses revealed cell type-dependent changes in secretion of variant laminin α5 laminin globular 4-5 (LG4-5) domain. Heterozygous and homozygous knockin mice with a corresponding variant of human LAMA5, p.Val3687Met, developed focal segmental glomerulosclerosis-like pathology with reduced laminin α5 and increased glomerular vinculin levels, which suggested that impaired cell adhesion may underlie this glomerulopathy. We also identified pulmonary defects such as bronchial deformity and alveolar dilation. Reexaminations of the family revealed phenotypes compatible with reduced laminin α5 and increased vinculin levels in affected tissues. Thus, the heterozygous p.Val3687Met variant may cause a new syndromic nephropathy with focal segmental glomerulosclerosis through possibly defective secretion of laminin α5. Enhanced vinculin may be a useful disease marker.

Severe spinal cord hypoplasia due to a novel ATAD3A compound heterozygous deletion.

Biallelic deletions extending into the ATPase family AAA-domain containing protein 3A (ATAD3A) gene lead to infantile lethality with severe pontocerebellar hypoplasia (PCH). However, only 12 such cases have been reported worldwide to date, and the genotype-phenotype correlations are not well understood. We describe cases associated with the same novel biallelic deletions of the ATAD3A and ATAD3B/3A regions in Japanese siblings with severe spinal cord hypoplasia and multiple malformations, including PCH, leading to neonatal death. The ATAD3A protein is essential for normal interaction between mitochondria and endoplasmic reticulum and is important for mitochondrial biosynthesis. The cases were evaluated using whole-genome sequencing for genetic diagnosis of mitochondrial disease. Spinal cord lesions associated with biallelic compound heterozygous deletion extending into the ATAD3A gene have not been reported. In addition, the ATAD3A deletion was 19 base pairs long, which is short compared with those reported previously. This deletion introduced a frameshift, resulting in a premature termination codon, and was expected to be a null allele. The pathological findings of the atrophic spinal cord showed gliosis and tissue destruction of the gray and white matter. We describe spinal cord lesions as a new central nervous system phenotype associated with a biallelic compound heterozygous deletion extending into the ATAD3A gene. Biallelic ATAD3A deletions should be considered in cases of mitochondrial disease with spinal cord hypoplasia and PCH.

Development of Leigh syndrome with a high probability of cardiac manifestations in infantile-onset patients with m.14453G > A.

The m.14453G > A mutation in MT-ND6 has been described in a few patients with mitochondrial encephalomyopathy, lactic acidosis, and stroke-like episodes or Leigh syndrome.However, the clinical spectrum and molecular characteristics are unclear.Here, we present four infantile-onset patients with m.14453G > A-associated Leigh syndrome. All four patients had brainstem lesions with basal ganglia lesions, and two patients had cardiac manifestations. Decreased ND6 protein expression and immunoreactivity were observed in patient-derived samples. There was no clear correlation between heteroplasmy levels and onset age or between heteroplasmy levels and phenotype; however, infantile onset was associated with Leigh syndrome.

Neonatal-onset mitochondrial disease: clinical features, molecular diagnosis and prognosis.

OBJECTIVE: Neonatal-onset mitochondrial disease has not been fully characterised owing to its heterogeneity. We analysed neonatal-onset mitochondrial disease in Japan to clarify its clinical features, molecular diagnosis and prognosis. DESIGN: Retrospective observational study from January 2004 to March 2020. SETTING: Population based. PATIENTS: Patients (281) with neonatal-onset mitochondrial disease diagnosed by biochemical and genetic approaches. INTERVENTIONS: None. MAIN OUTCOME MEASURES: Disease types, initial symptoms, biochemical findings, molecular diagnosis and prognosis. RESULTS: Of the 281 patients, multisystem mitochondrial disease was found in 194, Leigh syndrome in 26, cardiomyopathy in 38 and hepatopathy in 23 patients. Of the 321 initial symptoms, 236 occurred within 2 days of birth. Using biochemical approaches, 182 patients were diagnosed by mitochondrial respiratory chain enzyme activity rate and 89 by oxygen consumption rate. The remaining 10 patients were diagnosed using a genetic approach. Genetic analysis revealed 69 patients had nuclear DNA variants in 36 genes, 11 of 15 patients had mitochondrial DNA variants in five genes and four patients had single large deletion. The Cox proportional hazards regression analysis showed the effects of Leigh syndrome (HR=0.15, 95% CI 0.04 to 0.63, p=0.010) and molecular diagnosis (HR=1.87, 95% CI 1.18 to 2.96, p=0.008) on survival. CONCLUSIONS: Neonatal-onset mitochondrial disease has a heterogenous aetiology. The number of diagnoses can be increased, and clarity regarding prognosis can be achieved by comprehensive biochemical and molecular analyses using appropriate tissue samples.