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.

Novel ITPA variants identified by whole genome sequencing and RNA sequencing.

Approximately 80% of rare diseases have a genetic cause, and an accurate genetic diagnosis is necessary for disease management, prognosis prediction, and genetic counseling. Whole-exome sequencing (WES) is a cost-effective approach for exploring the genetic cause, but several cases often remain undiagnosed. We combined whole genome sequencing (WGS) and RNA sequencing (RNA-seq) to identify the pathogenic variants in an unsolved case using WES. RNA-seq revealed aberrant exon 4 and exon 6 splicing of ITPA. WGS showed a previously unreported splicing donor variant, c.263+1G>A, and a novel heterozygous deletion, including exon 6. Detailed examination of the breakpoint indicated the deletion caused by recombination between Alu elements in different introns. The proband was found to have developmental and epileptic encephalopathies caused by variants in the ITPA gene. The combination of WGS and RNA-seq may be effective in diagnosing conditions in proband who could not be diagnosed using WES.

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.

Genome sequencing and RNA-seq analyses of mitochondrial complex I deficiency revealed Alu insertion-mediated deletion in NDUFV2.

Isolated complex I deficiency is the most common cause of pediatric mitochondrial disease. Exome sequencing (ES) has revealed many complex I causative genes. However, there are limitations associated with identifying causative genes by ES analysis. In this study, we performed multiomics analysis to reveal the causal variants. We here report two cases with mitochondrial complex I deficiency. In both cases, ES identified a novel c.580G>A (p.Glu194Lys) variant in NDUFV2. One case additionally harbored c.427C>T (p.Arg143*), but no other variants were observed in the other case. RNA sequencing showed aberrant exon splicing of NDUFV2 in the unsolved case. Genome sequencing revealed a novel heterozygous deletion in NDUFV2, which included one exon and resulted in exon skipping. Detailed examination of the breakpoint revealed that an Alu insertion-mediated rearrangement caused the deletion. Our report reveals that combined use of transcriptome sequencing and GS was effective for diagnosing cases that were unresolved by ES.

A homozygous variant in NDUFA8 is associated with developmental delay, microcephaly, and epilepsy due to mitochondrial complex I deficiency.

Mitochondrial complex I deficiency is caused by pathogenic variants in mitochondrial and nuclear genes associated with complex I structure and assembly. We report the case of a patient with　NDUFA8-related mitochondrial disease. The patient presented with developmental delay, microcephaly, and epilepsy. His fibroblasts showed apparent biochemical defects in mitochondrial complex I. Whole-exome sequencing revealed that the patient carried a homozygous variant in NDUFA8. His fibroblasts showed a reduction in the protein expression level of not only NDUFA8, but also the other complex I subunits, consistent with assembly defects. The enzyme activity of complex I and oxygen consumption rate were restored by reintroducing wild-typeNDUFA8 cDNA into patient fibroblasts. The functional properties of the variant in NDUFA8 were also investigated using NDUFA8 knockout cells expressing wild-type or mutated NDUFA8 cDNA. These experiments further supported the pathogenicity of the variant in complex I assembly. This is the first report describing that the loss of NDUFA8, which has not previously been associated with mitochondrial disease, causes severe defect in the assembly of mitochondrial complex I, leading to progressive neurological and developmental abnormalities.

Mitochondrial ribosomal protein PTCD3 mutations cause oxidative phosphorylation defects with Leigh syndrome.

Pentatricopeptide repeat domain proteins are a large family of RNA-binding proteins involved in mitochondrial RNA editing, stability, and translation. Mitochondrial translation machinery defects are an expanding group of genetic diseases in humans. We describe a patient who presented with low birth weight, mental retardation, and optic atrophy. Brain MRI showed abnormal bilateral signals at the basal ganglia and brainstem, and the patient was diagnosed as Leigh syndrome. Exome sequencing revealed two potentially loss-of-function variants [c.415-2A>G, and c.1747_1748insCT (p.Phe583Serfs*3)] in PTCD3 (also known as MRPS39). PTCD3, a member of the pentatricopeptide repeat domain protein family, is a component of the small mitoribosomal subunit. The patient had marked decreases in mitochondrial complex I and IV levels and activities, oxygen consumption and ATP biosynthesis, and generalized mitochondrial translation defects in fibroblasts. Quantitative proteomic analysis revealed decreased levels of the small mitoribosomal subunits. Complementation experiments rescued oxidative phosphorylation complex I and IV levels and activities, ATP biosynthesis, and MT-RNR1 rRNA transcript level, providing functional validation of the pathogenicity of identified variants. This is the first report of an association of PTCD3 mutations with Leigh syndrome along with combined oxidative phosphorylation deficiencies caused by defects in the mitochondrial translation machinery.

A Comprehensive Genomic Analysis Reveals the Genetic Landscape of Mitochondrial Respiratory Chain Complex Deficiencies.

Mitochondrial disorders have the highest incidence among congenital metabolic disorders characterized by biochemical respiratory chain complex deficiencies. It occurs at a rate of 1 in 5,000 births, and has phenotypic and genetic heterogeneity. Mutations in about 1,500 nuclear encoded mitochondrial proteins may cause mitochondrial dysfunction of energy production and mitochondrial disorders. More than 250 genes that cause mitochondrial disorders have been reported to date. However exact genetic diagnosis for patients still remained largely unknown. To reveal this heterogeneity, we performed comprehensive genomic analyses for 142 patients with childhood-onset mitochondrial respiratory chain complex deficiencies. The approach includes whole mtDNA and exome analyses using high-throughput sequencing, and chromosomal aberration analyses using high-density oligonucleotide arrays. We identified 37 novel mutations in known mitochondrial disease genes and 3 mitochondria-related genes (MRPS23, QRSL1, and PNPLA4) as novel causative genes. We also identified 2 genes known to cause monogenic diseases (MECP2 and TNNI3) and 3 chromosomal aberrations (6q24.3-q25.1, 17p12, and 22q11.21) as causes in this cohort. Our approaches enhance the ability to identify pathogenic gene mutations in patients with biochemically defined mitochondrial respiratory chain complex deficiencies in clinical settings. They also underscore clinical and genetic heterogeneity and will improve patient care of this complex disorder.

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.

Fam57b (family with sequence similarity 57, member B), a novel peroxisome proliferator-activated receptor γ target gene that regulates adipogenesis through ceramide synthesis.

This report identifies a novel gene encoding Fam57b (family with sequence similarity 57, member B) as a novel peroxisome proliferator-activated receptor γ (PPARγ)-responsive transmembrane gene that is related to obesity. The gene was identified based on an integrated bioinformatics analysis of the following three expression profiling data sets: adipocyte differentiation of mouse stromal cells (ST2 cells), adipose tissues from obesity mice, and siRNA-mediated knockdown of Pparγ using ST2 cells. Fam57b consists of three variants expressed from different promoters and contains a Tram-Lag1-CLN8 domain that is related to ceramide synthase. Reporter and ChIP assays showed that Fam57b variant 2 is a bona fide PPARγ target gene in ST2 cells. Fam57b was up-regulated during adipocyte differentiation, suggesting that FAM57B is involved in this process. Surprisingly, FAM57B overexpression inhibited adipogenesis, and siRNA-mediated knockdown promoted adipocyte differentiation. Analysis of the ceramide content by lipid assay found that ceramides were in fact augmented in FAM57B-overexpressing ST2 cells. We also confirmed that ceramide inhibits adipogenesis. Therefore, the aforementioned results of FAM57B overexpression and siRNA experiments are reconciled by ceramide synthesis. In summary, we present in vitro evidence showing that PPARγ regulates Fam57b transcription during the adipogenesis of ST2 cells. In addition, our results suggest that PPARγ activation contributes to the regulation of ceramide metabolism during adipogenesis via FAM57B.

The paired-box homeodomain transcription factor Pax6 binds to the upstream region of the TRAP gene promoter and suppresses receptor activator of NF-κB ligand (RANKL)-induced osteoclast differentiation.

Osteoclast formation is regulated by balancing between the receptor activator of nuclear factor-κB ligand (RANKL) expressed in osteoblasts and extracellular negative regulatory cytokines such as interferon-γ (IFN-γ) and interferon-ß (IFN-ß), which can suppress excessive bone destruction. However, relatively little is known about intrinsic negative regulatory factors in RANKL-mediated osteoclast differentiation. Here, we show the paired-box homeodomain transcription factor Pax6 acts as a negative regulator of RANKL-mediated osteoclast differentiation. Electrophoretic mobility shift and reporter assays found that Pax6 binds endogenously to the proximal region of the tartrate acid phosphatase (TRAP) gene promoter and suppresses nuclear factor of activated T cells c1 (NFATc1)-induced TRAP gene expression. Introduction of Pax6 retrovirally into bone marrow macrophages attenuates RANKL-induced osteoclast formation. Moreover, we found that the Groucho family member co-repressor Grg6 contributes to Pax6-mediated suppression of the TRAP gene expression induced by NFATc1. These results suggest that Pax6 interferes with RANKL-mediated osteoclast differentiation together with Grg6. Our results demonstrate that the Pax6 pathway constitutes a new aspect of the negative regulatory circuit of RANKL-RANK signaling in osteoclastogenesis and that the augmentation of Pax6 might therefore represent a novel target to block pathological bone resorption.

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.

Id4, a new candidate gene for senile osteoporosis, acts as a molecular switch promoting osteoblast differentiation.

Excessive accumulation of bone marrow adipocytes observed in senile osteoporosis or age-related osteopenia is caused by the unbalanced differentiation of MSCs into bone marrow adipocytes or osteoblasts. Several transcription factors are known to regulate the balance between adipocyte and osteoblast differentiation. However, the molecular mechanisms that regulate the balance between adipocyte and osteoblast differentiation in the bone marrow have yet to be elucidated. To identify candidate genes associated with senile osteoporosis, we performed genome-wide expression analyses of differentiating osteoblasts and adipocytes. Among transcription factors that were enriched in the early phase of differentiation, Id4 was identified as a key molecule affecting the differentiation of both cell types. Experiments using bone marrow-derived stromal cell line ST2 and Id4-deficient mice showed that lack of Id4 drastically reduces osteoblast differentiation and drives differentiation toward adipocytes. On the other hand knockdown of Id4 in adipogenic-induced ST2 cells increased the expression of Ppargamma2, a master regulator of adipocyte differentiation. Similar results were observed in bone marrow cells of femur and tibia of Id4-deficient mice. However the effect of Id4 on Ppargamma2 and adipocyte differentiation is unlikely to be of direct nature. The mechanism of Id4 promoting osteoblast differentiation is associated with the Id4-mediated release of Hes1 from Hes1-Hey2 complexes. Hes1 increases the stability and transcriptional activity of Runx2, a key molecule of osteoblast differentiation, which results in an enhanced osteoblast-specific gene expression. The new role of Id4 in promoting osteoblast differentiation renders it a target for preventing the onset of senile osteoporosis.

Changes in histopathology and heteroplasmy rates over 8 years and effectiveness of taurine supplementation in a patient with mitochondrial nephropathy caused by MT-TL1 mutation: A case report.

The m.3243A > G mutation in the mitochondrially encoded tRNA leucine 1 (MT-TL1) gene is known to cause mitochondrial nephropathy. However, its long-term effects of the m.3243A > G mutation on renal histopathology or heteroplasmy rates remain unknown. Here we present the case of a female patient who underwent renal biopsy at 34 years of age to investigate the reason for a low estimated glomerular filtration rate (eGFR) of 47.9 mL/min/1.73 m2. Light microscopy revealed nephrosclerosis with granular swollen epithelial cells (GSECs) in the renal tubules. Genetic testing revealed the m.3243A > G mutation in the MT-TL1 gene. Over a follow-up period of 8 years, the eGFR declined at a rate of 1.50 mL/min/1.73 m2/year. A second renal biopsy was performed at the age of 42 years; the patient's glomerular sclerosis rate had increased from 45.5% to 63.2%, and the frequency of GSECs in the collecting ducts had increased from 5.8% to 20.8%. Furthermore, the heteroplasmy rate in blood cells and urinary sediment cells increased from 9% to 20% and 20% to 53%, respectively. Taurine therapy was initiated just after the second kidney biopsy. To date, after approximately 3 years of taurine administration, the rate of eGFR decline has markedly decreased to 0.26 mL/min/1.73 m2/year. This experience suggests that an increased heteroplasmy rate may be associated with the progression of mitochondrial nephropathy caused by MT-TL1 mutation. Furthermore, our case is the first to suggest the effectiveness of taurine for mitochondrial nephropathy caused by the m.3243A > G mutation in the MT-TL1 gene.

Focal segmental glomerulosclerosis with a mutation in the mitochondrially encoded NADH dehydrogenase 5 gene: A case report.

NADH dehydrogenase 5 (ND5) is one of 44 subunits composed of Complex I in mitochondrial respiratory chain. Therefore, a mitochondrially encoded ND5 (MT-ND5) gene mutation causes mitochondrial oxidative phosphorylation (OXPHOS) disorder, resulting in the development of mitochondrial diseases. Focal segmental glomerulosclerosis (FSGS) which had podocytes filled with abnormal mitochondria is induced by mitochondrial diseases. An MT-ND5 mutation also causes FSGS. We herein report a Japanese woman who was found to have proteinuria and renal dysfunction in an annual health check-up at 29 years old. Because her proteinuria and renal dysfunction were persistent, she had a kidney biopsy at 33 years of age. The renal histology showed FSGS with podocytes filled with abnormal mitochondria. The podocytes also had foot process effacement and cytoplasmic vacuolization. In addition, the renal pathological findings showed granular swollen epithelial cells (GSECs) in tubular cells, age-inappropriately disarranged and irregularly sized vascular smooth muscle cells (AiDIVs), and red-coloured podocytes (ReCPos) by acidic dye. A genetic analysis using peripheral mononuclear blood cells and urine sediment cells detected the m.13513 G > A variant in the MT-ND5 gene. Therefore, this patient was diagnosed with FSGS due to an MT-ND5 gene mutation. Although this is not the first case report to show that an MT-ND5 gene mutation causes FSGS, this is the first to demonstrate podocyte injuries accompanied with accumulation of abnormal mitochondria in the cytoplasm.

Adult-onset Leigh Syndrome with a m.9176T>C Mutation Manifested As Reversible Cerebral Vasoconstriction Syndrome.

A 26-year-old woman developed a sudden headache, ptosis, and diplopia. Magnetic resonance imaging and angiography demonstrated a symmetrical lesion from the midbrain to the brainstem, involving the solitary nucleus and multifocal cerebral artery narrowing. Reversible cerebral vasoconstriction syndrome (RCVS) was suspected, and the patient improved after vasodilatation. Leigh syndrome was suspected due to the elevated serum pyruvate level, so mitochondrial DNA was analyzed, and an m.9176T>C mutation was detected. The final diagnosis was adult-onset Leigh syndrome manifesting as RCVS. An uncontrolled baroreflex due to a solitary nuclear lesion or endothelial dysfunction may have contributed to her unique presentation.

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.

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.

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.