key: cord-0824482-ukl2145l
authors: Sanford, Erica; Jones, Marilyn C.; Brigger, Matthew; Hammer, Monia; Giudugli, Lucia; Kingsmore, Stephen F.; Dimmock, David; Bainbridge, Matthew N.
title: Postmortem diagnosis of PPA2-associated sudden cardiac death from dried blood spot in a neonate presenting with vocal cord paralysis
date: 2020-10-03
journal: Cold Spring Harb Mol Case Stud
DOI: 10.1101/mcs.a005611
sha: 944382c565ff37c67c00a32d387f8087464d3e3a
doc_id: 824482
cord_uid: ukl2145l

Biallelic variants in inorganic pyrophosphatase 2 (PPA2) are known to cause infantile sudden cardiac failure (OMIM #617222), but relatively little is known about phenotypic variability of these patients prior to their death. We report a 5-wk-old male with bilateral vocal cord paralysis and hypertension who had a sudden unexpected cardiac death. Subsequently, molecular autopsy via whole-genome sequencing from newborn dried blood spot identified compound heterozygous mutations in PPA2, with a paternally inherited, pathogenic missense variant (c.514G > A; p.Glu172Lys) and a novel, maternally inherited missense variant of uncertain significance (c.442A > T; p.Thr148Ser). This report expands the presenting phenotype of patients with PPA2 variants. It also highlights the utility of dried blood spots for postmortem molecular diagnosis.

Sudden unexpected death in infancy (SUDI) often remains unexplained despite comprehensive postmortem investigations. More recently, molecular autopsy has yielded answers for some of these patients (Torkamani et al. 2016) . Variants that cause SUDI are typically very rare and may be dominant or recessive (Bainbridge 2020) . Damaging variants in cardiacassociated genes have been found in 20%-34% of SUDI cases (Klaver et al. 2011; Hertz et al. 2016) , including genes encoding cardiac ion channels and genes implicated in lethal cardiomyopathies (Klaver et al. 2011; Hertz et al. 2016 ). Next-generation sequencing (NGS) facilitates simultaneous analysis of a much larger number of genes, and this technologic advancement is likely responsible for the higher variant detection rate in cardiac genes in more recent publications (Hertz et al. 2016) .

In 2016, Kennedy et al. and Guimier et al. separately identified biallelic missense variants in the gene PPA2 in several infants with sudden cardiac arrest (Guimier et al. 2016; Kennedy et al. 2016) . PPA2 is a nuclear-encoded enzyme that has a critical role in mitochondrial maintenance and function (Guimier et al. 2016) . The gene encodes an inorganic pyrophosphatase 2 that is responsible for hydrolyzing inorganic pyrophosphate (also known as diphosphate, PPi) into two molecules of inorganic phosphate (also known as orthophosphate, Pi) (Mansurova 1989) . This enzymatic reaction is essential for energy metabolism as it provides substrate for ATP synthesis (Mansurova 1989; Kennedy et al. 2016) . It is also important for synthesis of DNA, RNA, protein, polysaccharides, and lipids (Mansurova 1989; Guimier et al. 2016) .

To date there have been 19 patients described in the literature who suffered from sudden cardiac death attributed to disease-causing variants in PPA2 (https://omim.org/entry .617222, Guimier et al. 2016; Kennedy et al. 2016; Vasilescu et al. 2018; Phoon et al. 2020) . In addition to the first patients identified by Kennedy et al. and Guimier et al., a Finnish cardiomyopathy cohort (KidCMP) of 66 children with severe cardiomyopathy found two siblings with compound heterozygous PPA2 variants who both died from rapidly progressive cardiomyopathy in infancy (Vasilescu et al. 2018) . Additionally, in 2020, Phoon et al. reported two siblings who also died unexpectedly in infancy, and on postmortem analysis a molecular diagnosis of biallelic PPA2 variants was made (Phoon et al. 2020) .

The majority of published patients were <1 yr of age at the time of cardiac arrest, with no comorbidities at their time of death (Table 1 ; Guimier et al. 2016; Kennedy et al. 2016; Vasilescu et al. 2018; Phoon et al. 2020) . This report describes a 5-wk-old who was admitted and discharged from the Pediatric Intensive Care Unit (PICU) twice for symptoms of respiratory distress due to vocal cord paralysis, who then re-presented to the emergency room and had a cardiac arrest. Molecular autopsy via newborn dried blood spot identified compound heterozygous variants in PPA2. The novel clinical presentation in this infant potentially expands the spectrum of clinical phenotype associated with PPA2 variants.

The male proband initially presented to his pediatrician at 3 wk of age for a weight check. He was noted to have fussiness, decreased oral intake, and emesis in the prior 3 d. He had previously been diagnosed at 1 wk of age with laryngomalacia by his pediatrician. He was sent to the emergency room for respiratory distress and then transferred to the PICU at Rady Children's Hospital, San Diego, for continuous positive airway pressure with 50% oxygen Because of persistent stridor, a bedside fiberoptic bronchoscopy was performed and revealed mild to moderate laryngomalacia as well as bilateral vocal cord hypomobility. The patient then went to the operating room for a more extensive airway evaluation, which confirmed bilateral vocal cord immobility. A dysphagia study showed aspiration of thin liquids, and his feeds were subsequently thickened. Brain magnetic resonance imaging was unremarkable. His oxygen requirement resolved and he was discharged home after 6 d in the hospital.

The patient returned the following day for recurrence of respiratory distress and was readmitted to the PICU. Oxygen support was initiated and a nasogastric tube was placed for feeding. Systemic hypertension was noted. Computed tomography angiogram of the chest and abdomen was remarkable for the finding of two left renal arteries without evidence of renal artery stenosis, and renal ultrasonography with Doppler waveform showed diminished diastolic flow consistent with elevated renal vascular resistance. The infant was started on amlodipine at the recommendation of the consulting nephrologist. Suture lateralization of the right vocal cord was performed to address persistent stridor. He was again weaned off of oxygen supplementation and discharged home on nasogastric tube feeds after 12 d in the hospital.

One day after his second discharge from the hospital, the patient returned to the emergency room following four episodes of vomiting initially thought to be due to malposition of the enteric tube. A radiograph was obtained that showed the tube terminating in the duodenum. During the period of observation in the emergency room, the patient had an unprovoked episode of unresponsiveness and cyanosis and subsequently had a bradycardic arrest. Chest compressions were initiated and the patient was intubated with an endotracheal tube. The underlying rhythm was primarily pulseless electrical activity (PEA), with one incidence of pulseless ventricular tachycardia, for which the patient was unsuccessfully defibrillated. Despite 1 h of cardiopulmonary resuscitation, including multiple doses of epinephrine, he was unable to be resuscitated.

At autopsy, the lungs were heavy (115 gm vs. expected 75 gm) with lymphocytic infiltration in the larynx and trachea, pulmonary congestion, and increased alveolar macrophages, consistent with a known history of silent aspiration. The liver weighed 227 gm (expected 140 gm) and evidenced mild to moderate panacinar, anisomorphic, mostly macrovesicular steatosis. The heart was large (34 gm vs. expected 23 gm; Pryce et al. 2014) but was grossly and microscopically normal. There was no dilation or hypertrophy of the cardiac chambers. Skeletal muscle was microscopically normal. Hepatic steatosis was noted. There was no renal pathology. The brain was edematous but structurally normal.

Proband sequencing (whole-genome sequencing [WGS] ) was undertaken to determine a possible cause of death. This patient was found to be compound heterozygous for a paternally inherited, known pathogenic missense variant (c.514G > A; p.Glu172Lys; Chr 4: 106359121) and a maternally inherited, novel missense variant of uncertain significance (c.442A > T; p.Thr148Ser; Chr 4:106359193) in the PPA2 gene (Table 2) . No additional sequence variants or copy-number variants (CNVs) were identified by WGS that could explain or be related to the patient's clinical presentation.

The paternally inherited p.Glu172Lys variant in PPA2 has been previously reported as a compound heterozygous change in three individuals from two families with sudden unexpected cardiac arrest in infancy (Guimier et al. 2016 ) and in two individuals from one family who developed a rapidly progressive dilated cardiomyopathy and cardiac failure, with only a few days from disease onset to death (Vasilescu et al. 2018) . This glutamine to lysine substitution is at a highly conserved residue and is predicted to disrupt at least three hydrogen bonds between interacting protein chains near the surface of the enzyme's active site and subsequently impair the enzymatic function of PPA2 (Guimier et al. 2016 ). This variant is present in the heterozygous state in the gnomAD population database at a frequency of 0.053% (147/275080) overall and 0.095% in non-Finnish Europeans. Based on the available evidence, the c.514G > A (p.Glu172Lys) variant was classified as Pathogenic according to American College of Medical Genetics and Genomics (ACMG) and Association for Molecular Pathology (AMP) standards and guidelines (Richards et al. 2015) . The specific scoring criteria used for classification are listed in Supplemental Table S1 .

The maternally inherited c.442A > T is a novel missense variant. This variant has not been previously reported or functionally characterized in the literature to our knowledge. It is absent from the gnomAD population database and thus is presumed to be rare. The c.442A > T (p.Thr148Ser) variant affects an amino acid that is conserved down to Arabidopsis thaliana and is predicted by multiple in silico tools to have a deleterious effect on protein function. The specific scoring criteria used for classification are listed in Supplemental Table S1 . Results of protein sequence alignment from human to yeast shows remarkable sequence conservation of the Thr and its adjoining amino acids. This suggests an important evolutionary pressure to prevent a change of this Thr. We used three different algorithms to assess the mutations in the corresponding threonines (T100 in pdb code: 6C45 and T97 in pdb code: 4QLZ). These three methods of prediction are mCSM, SDM, and DUET (Topham et al. 1997; Worth et al. 2011; Pires et al. 2014 ). All three prediction tools point to a destabilizing effect of the missense change in both structures (Supplemental Table S2 ). Notably, the c.514G > A (p.Glu172Lys) was characterized as disruptive in a functional complementation assay in yeast where the human wild-type construct rescued a yeast growth defect, whereas the mutant p. Glu172Lys could not. Furthermore, the alteration was shown to affect the yeast mitochondrial function (Guimier et al. 2016) . Based on the available evidence, the c.442A > T (p.Thr148Ser) variant was classified as a variant of uncertain significance.

NGS is increasingly being utilized for postmortem diagnosis of sudden unexplained death in infancy (Hertz et al. 2016) . The genetic landscape for this condition is still being defined, in part because of the rare nature of this event and the novelty of the methodologic tools available (Vasilescu et al. 2018) . Genes that were not previously analyzed in cases of sudden unexpected death in infancy have more recently been implicated as a probable cause of death (Hertz et al. 2016) . The use of an NGS approach has shown diagnostic superiority compared with standard genetic testing and thus also makes for an ideal molecular autopsy tool, especially because of the ability to analyze recently discovered genes and genes not found in standard gene panels and to interrogate the genome for disorders that were not included in the premortem differential diagnosis because of an atypical or previously unrecognized presentation (Miller et al. 2015; Clark et al. 2018 Clark et al. , 2019 Lionel et al. 2018; Vasilescu et al. 2018; Sanford et al. 2019) . WGS is an ideal tool for the molecular autopsy as analysis is not limited by premortem differential diagnosis. However, it is worth noting that WGS does have its own limitations, including difficulty identifying trinucleotide repeat disorders as well as pathogenic deep intronic mutations (Smith et al. 2016) .

We identified compound heterozygous mutations in PPA2 in a 5-wk-old infant who was admitted and discharged from the PICU twice for symptoms related to an unexplained vocal cord paralysis and then sustained a sudden unexpected cardiac arrest. The p.Glu172Lys variant identified in this patient has been reported as pathogenic for a sudden cardiac arrest in infancy and progressive lethal cardiomyopathy in multiple families (Guimier et al. 2016; Vasilescu et al. 2018) , whereas this is the first report of the p.Thr148Ser variant to our knowledge and with the current available evidence it was classified as a variant of uncertain significance. The molecular diagnosis in this case was an unexpected finding, as the patient did not display any evidence of cardiac pathology (by telemetry or echocardiography) during two PICU admissions in the weeks prior to his death. This case underscores the knowledge that disease-causing PPA2 variants can cause sudden death in infancy without preceding symptoms (Phoon et al. 2020 ).

The p.Glu172Lys variant found in our patient has been previously reported as a compound heterozygous change in three individuals from two families with sudden unexpected cardiac arrest in infancy (Guimier et al. 2016 ) and in two brothers who developed a rapidly progressive dilated cardiomyopathy (Vasilescu et al. 2018) . The siblings in the Vasilescu et al. cohort displayed a different phenotype than patients described by Guimier et al. (2016) , Kennedy et al. (2016) , and Phoon et al. (2020) as they had a few days between disease presentation and death (Vasilescu et al. 2018 ). The second brother was actually followed with serial echocardiograms, which were normal until the sudden manifestation of disease and rapid progression to cardiac death at 5 mo of age (Vasilescu et al. 2018) . The patient in our case also had a documented normal echocardiogram ∼3 wk prior to his death, although he experienced sudden death and no dilation was noted on cardiac autopsy.

As described previously, PPA2 variants occupy a range of the clinical spectrum, with sudden cardiac death documented from the neonatal period through adolescence, although the majority of patients died at <1 yr of age (Guimier et al. 2016; Kennedy et al. 2016; Vasilescu et al. 2018; Phoon et al. 2020) . Two siblings of probands who suffered from sudden cardiac death were actually identified by genetic testing to have the same PPA2 variants, but received implantable defibrillators and were reportedly alive in their 30s with no events (Kennedy et al. 2016) . Some cases of sudden cardiac arrest appeared to be associated with viral illness or alcohol ingestion, whereas for others there was no evidence of provocation (Kennedy et al. 2016; Phoon et al. 2020) . Findings on necroscopy have ranged from nonspecific to cardiac inflammation, fibrosis, acute myocyte loss, and/or myocarditis (Kennedy et al. 2016; Vasilescu et al. 2018; Phoon et al. 2020 ). Similar to the two patients described by Phoon et al. the autopsy in our case was nonspecific (Phoon et al. 2020) .

PPA2 was initially identified in yeast and was not associated with human disease until recently (Guimier et al. 2016; Kennedy et al. 2016; Vasilescu et al. 2018 ). The nuclear gene encodes an inorganic pyrophosphatase that functions within the mitochondrial matrix and hydrolyzes inorganic pyrophosphate (diphosphate) to two molecules of inorganic phosphate (orthophosphate), which is necessary for energy metabolism and numerous cellular reactions, including synthesis of DNA, RNA, proteins, polysaccharides, and lipids (Mansurova 1989; Guimier et al. 2016) . Functional characterization in yeast has demonstrated its essential role in multiple cellular functions: In vitro PPA2 mutants, including the p.Glu172Lys variant, have demonstrated decreased oxidative respiration, reduced ATP synthesis, inability to maintain electrical potential across the mitochondrial inner membrane, and reduction in complex IV of the mitochondrial respiratory chain (Guimier et al. 2016; Kennedy et al. 2016; Vasilescu et al. 2018) . Likewise, Eschericia coli expressing recombinant protein with the p.Glu172Lys variant had only 5%-10% residual activity (Kennedy et al. 2016) . Studies of human cardiac and skeletal tissue have shown decreased mitochondrial complex I and IV in affected patients (Phoon et al. 2020) . Examination of cardiac tissue on autopsy showed decreased levels of mitochondrial complex I subunit (two patients), decreased amount of PPA2 protein in heart tissue, and evidence indicating fewer mitochondria overall (Kennedy et al. 2016) . Unfortunately, the mitochondria were not visualized on autopsy in the case presented here.

Further studies are needed to truly understand how the in vitro observations translate to cardiac dysfunction and ultimately cardiac death. Guimier hypothesized that the high energy needs of cardiac myocytes may preclude return of spontaneous circulation (ROSC) despite early and optimal resuscitation in patients with pathologic PPA2 variants, potentially because of an inability of the mitochondria to meet the energy demands of the cell (Guimier et al. 2016) . This theory may be especially applicable to our patient as the majority of pediatric patients with in-hospital cardiac arrests do achieve ROSC (Sutton et al. 2018) .

Our patient displayed a novel clinical presentation of bilateral vocal cord paresis. We suggest that this be considered as a potential phenotype expansion for PPA2 mutations. As noted, PPA2 variants result in a disorder affecting mitochondria. Several patients described had features consistent with classical descriptions of mitochondrial diseases, although nonspecific: hypotonia, myopathy, seizures, lactic acidosis, vomiting, and failure to thrive (Guimier et al. 2016; Kennedy et al. 2016) . Vocal cord paralysis has been previously described in patients with a number of mitochondrial disorders (Supplemental Table S3 ), including Kearns-Sayre syndrome, Leigh disease, Brown-Vialetto-van Laere syndrome, Charcot-Marie-Tooth disease, and flavin adenine dinucleotide synthase (FADS) deficiency (Lin et al. 1999; Diamantopoulou et al. 2001; Sevilla et al. 2008; Chen et al. 2010; Ando et al. 2017; Thulasi et al. 2017; Yamada et al. 2019) . Phenotypic expansion could have important implications for which patients receive more comprehensive genetic testing, such as WGS, and potentially a change in their clinical management-for example, implantation of a defibrillator, if appropriate, as in the cases of two affected siblings (Kennedy et al. 2016) . However, the autopsy in this particular case was nonspecific and did not include an examination of the mitochondria. Additionally, the patient died at the age of 5 wk, potentially prior to the manifestation of additional disease features. Thus, although we present this case as an interesting possibility, we are unable to say with certainty that bilateral vocal cord paresis is a phenotypic expansion of PPA2-related cardiac disease. More patients will need to be identified and more in vitro studies conducted to follow up on this potential association.

Molecular autopsy has the potential to influence family planning and familial screening and to more accurately define the genetic landscape of unexpected sudden death in infancy, regardless of premortem assumptions. This case illustrates that example, but more research is required to further elucidate the pathophysiology at play, and whether more comprehensive genetic screening of infants should be pursued earlier on.

Parental consent was obtained for genome sequencing. Proband WGS was performed at Rady Children's Institute for Genomic Medicine as described (Torkamani et al. 2016; Kingsmore et al. 2019) . The newborn dried blood spot was released from the state laboratory for this purpose. Following DNA isolation from newborn dried blood spot (patient) and whole blood (parents), sequencing libraries were generated using the KAPA HyperPlus Library Prep kit (Roche) according to the manufacturer's instructions. Paired-end sequencing was performed on a NovaSeq 6000 instrument and S1 flow cell (Illumina). Read alignment to the reference human genome assembly GRCh37/hg19 and single-nucleotide variant (SNV)/ insertion-deletion (indel) calling was performed using the DRAGEN Bio-IT Platform (Illumina). CNV calling was performed using CNVnator and Manta (Kingsmore et al. 2019) .

SNVs/indels and CNVs were annotated and analyzed using Fabric Enterprise version 6.5.2 (Fabric Genomics). The median genomic coverage was 35× with 93.3% of OMIM genes fully covered at ≥10× (Table 3) .

Human Phenotype Ontology (HPO) terms used during variant interpretation included laryngomalacia (HP:0001601), bilateral vocal cord paresis (HP:0012822), aspiration (HP:0002835), hypoglycemia (HP:0001943), gastroesophageal reflux (HP:0002020), respiratory failure (HP:0002878), and cardiac arrest (HP:0001695). SNVs/indels were filtered to retain variants with allelic balance between 0.3 and 0.7 and allele frequency <0.5% in the Genome Aggregation Database (gnomAD V2.1.1), prioritized by Phevor Gene Rank (Kingsmore et al. 2019) , and classified according to ACMG and AMP standards and guidelines (Richards et al. 2015) . CNVs were filtered to retain intragenic and intergenic variants shown to have a disease association or possibly associated with disease and classified according to ACMG and ClinGen recommendations (Riggs et al. 2020) . Sequence variants of interest were orthogonally confirmed by Sanger sequencing, whereas CNVs of interest are typically confirmed by MLPA (none in this case).

The causative variants have been submitted to ClinVar (https://www.ncbi.nlm.nih.gov/ clinvar/) with accession numbers SCV001250701.1 and SCV001250702.1. Raw sequencing data could not be deposited because patient consent was not obtained.

Informed and signed consent forms were obtained for all sequenced individuals in this study. The project is approved by the Institutional Review Board of the University of California at San Diego under protocol #20172719 and has received nonsignificant risk status in a pre-Investigational Device Exemption submission to the Food and Drug Administration. 

Clinical and genetic diversities of Charcot-Marie-Tooth disease with MFN2 mutations in a large case study

Determining the incidence of rare diseases

CMT2C with vocal cord paresis associated with short stature and mutations in the TRPV4 gene

Meta-analysis of the diagnostic and clinical utility of genome and exome sequencing and chromosomal microarray in children with suspected genetic diseases

Diagnosis of genetic diseases in seriously ill children by rapid whole-genome sequencing and automated phenotyping and interpretation

Kearns-Sayre syndrome: presenting with vocal fold palsy

Biallelic PPA2 mutations cause sudden unexpected cardiac arrest in infancy

Genetic investigations of sudden unexpected deaths in infancy using next-generation sequencing of 100 genes associated with cardiac diseases

Sudden cardiac death due to deficiency of the mitochondrial inorganic pyrophosphatase PPA2

A randomized, controlled trial of the analytic and diagnostic performance of singleton and trio, rapid genome and exome sequencing in ill infants

Cardiac ion channel mutations in the sudden infant death syndrome

Vocal cord paralysis and hypoventilation in a patient with suspected Leigh disease

Improved diagnostic yield compared with targeted gene sequencing panels suggests a role for whole-genome sequencing as a first-tier genetic test

Inorganic pyrophosphate in mitochondrial metabolism

A 26-hour system of highly sensitive whole genome sequencing for emergency management of genetic diseases

Sudden unexpected death in asymptomatic infants due to PPA2 variants

DUET: a server for predicting effects of mutations on protein stability using an integrated computational approach

Reference ranges for organ weights of infants at autopsy: results of >1,000 consecutive cases from a single centre

Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology

Technical standards for the interpretation and reporting of constitutional copy-number variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics (ACMG) and the Clinical Genome Resource (ClinGen)

Rapid whole genome sequencing has clinical utility in children in the PICU

Vocal cord paresis and diaphragmatic dysfunction are severe and frequent symptoms of GDAP1-associated neuropathy

Whole-exome sequencing and whole-genome sequencing in critically ill neonates suspected to have single-gene disorders

Chest compression rates and pediatric in-hospital cardiac arrest survival outcomes

A case of Brown-Vialetto-van Laere syndrome due to a novel mutation in SLC52A3 gene: clinical course and response to riboflavin

Prediction of the stability of protein mutants based on structural environment-dependent amino acid substitution and propensity tables

Molecular autopsy for sudden unexpected death

Genetic basis of severe childhood-onset cardiomyopathies

SDM-a server for predicting effects of mutations on protein stability and malfunction

Flavin adenine dinucleotide synthase deficiency due to FLAD1 mutation presenting as multiple acyl-CoA dehydrogenation deficiency-like disease: a case report

We are extremely thankful for the family's willingness to participate in this study.

M.N.B. conceived of the study. E.S., L.G., and M.H. wrote the manuscript. All authors reviewed and edited the manuscript.

Funding for this study was provided in part by the National Heart, Lung, and Blood Institute (NHLBI) (R01HL145175-02) and by the Scripps Research Institute (5 UL1 TR002550).