Temperature-sensitive mutation in PEX1 moderates the phenotypes of peroxisome deficiency disorders

The peroxisome biogenesis disorders (PBDs), including Zellweger syndrome (ZS), neonatal adrenoleukodystrophy (NALD) and infantile Refsum disease (IRD), are autosomal recessive diseases caused by deficiency of peroxisome assembly as well as malfunction of peroxisomes, where >10 genotypes have been reported. ZS patients manifest the most severe clinical and biochemical abnormalities, while those with NALD and IRD show the least severity and the mildest features, respectively. PEX1 is the causative gene for PBDs of complementation group I (CG1), the highest incidence PBD, and encodes the peroxin, Pex1p, a member of the AAA ATPase family. In the present work, we found that peroxisomes were morphologically and biochemically formed at 30 but not 37 degrees C, in the fibroblasts from all CG1 IRD patients examined, whereas almost no peroxisomes were seen in ZS and NALD cells, even at 30 degrees C. A point missense mutation, G843D, was identified in the PEX1 allele of most CG1 IRD patients. The mutant PEX1, termed HsPEX1G843D, gave rise to the same temperature-sensitive phenotype on CG1 CHO cell mutants upon transfection. Collectively, these results demonstrate temperature-sensitive peroxisome assembly to be responsible for the mildness of the clinical features of PEX1 - defective IRD of CG1.      

Nonsense and temperature-sensitive mutations in PEX13 are the cause of complementation group H of peroxisome biogenesis disorders.

Peroxisome biogenesis disorders, including Zellweger syndrome (ZS), neonatal adrenoleukodystrophy (NALD) and infantile Refsum disease, are lethal hereditary diseases caused by abnormalities in peroxisomal assembly. To date, 12 genotypes have been identified. We now have evidence that the complete human cDNA encoding Pex13p, an SH3 protein of a docking factor for the peroxisome targeting signal 1 receptor (Pex5p), rescues peroxisomal matrix protein import and its assembly in fibroblasts from PBD patients of complementation group H. In addition, we detected mutations on the human PEX13 cDNA in two patients of group H. A severe phenotype of a ZS patient (H-02) was homozygous for a nonsense mutation, W234ter, which results in the loss of not only the SH3 domain but also the putative transmembrane domain of Pex13p. A more mildly affected NALD patient (H-01), whose fibroblasts showed the temperature-sensitive (TS) phenotype, was homozygous for a missense mutation in the SH3 domain of Pex13p, I326T. This mutant PEX13 cDNA expression in a PEX13-defective CHO mutant showed I326T to be a TS mutation and thus suggested that Pex13p with the I326T mutation in the SH3 domain is stable at 30 degrees C but is somewhat unstable at 37 degrees C.     

Prader-Willi syndrome: causes of death in an international series of 27 cases

The peroxisome biogenesis disorders (PBDs) with generalized peroxisomal dysfunction include Zellweger syndrome (ZS), neonatal adrenoleukodystrophy (NALD), and infantile Refsum disease (IRD). There is clinical, biochemical, and genetic overlap among the three phenotypes, also known as Zellweger spectrum disorders. Clinical distinctions between the phenotypes are not sharply defined. Only limited sources are available to serve as a background for prognosis in PBD, especially in case of prolonged survival. We delineated the natural history of 31 PBD patients (age 1.2-24 years) through systematic clinical and biochemical investigations. We excluded classical ZS from our study, and included all patients with a biochemically confirmed generalized peroxisomal disorder over 1 year of age, irrespective of the previously diagnosed phenotype. The initial clinical suspicion, age at diagnosis, growth, development, neurological symptoms, organ involvements, and survival are summarized. Common to all patients were cognitive and motor dysfunction, retinopathy, sensorineural hearing impairment, and hepatic involvement. Many patients showed postnatal growth failure, 10 patients displayed hyperoxaluria of whom 4 had renal stones. Motor skills ranged from sitting with support to normal gait. Speech development ranged from non-verbal expression to grammatical speech and comprehensive reading. The neurodevelopmental course was variable with stable course, rapid decline with leukodystrophy, spinocerebellar syndrome, and slow decline over a wide range of faculties as outcome profiles. At the molecular level, 21 patients had mutations in the PEX1 gene. The two most common PEX1 mutations were the G843D (c.2528G-->A) missense and the c.2097insT frameshift mutation. Patients having the G843D/G843D or the G843D/c.2097insT genotypes were compared. Patients homozygous for G843D generally had a better developmental outcome. However, one patient who was homozygous for the "mild" G843D mutation had an early lethal disease, whereas two other patients had a phenotype overlapping with the G843D/c.2097insT group. This indicates that next to the PEX1 genotype other yet unknown factors determine the ultimate phenotype.     

Mutations in PEX10 is the cause of Zellweger peroxisome deficiency syndrome of complementation group B

Peroxisome biogenesis disorders (PBD), such as Zellweger syndrome, are autosomal recessive diseases caused by a deficiency in peroxisome assembly as well as a malfunction of the peroxisomes, where at least 10 genotypes have been reported. We have isolated a human PEX10 cDNA (HsPEX10) by an expressed sequence tag homology search on a human DNA database using yeast PEX10 from Hansenula polymorpha, followed by screening of a human liver cDNA library. This cDNA encodes a peroxisomal protein (a peroxin Pex10p) comprising 326 amino acids, with two putative transmembrane segments and a C3HC4zinc finger RING motif. Both the N- and C-terminal regions of Pex10p are exposed to the cytosol, as assessed by an expression study of epitope-tagged Pex10p. HsPEX10 expression morphologically and biochemically restored peroxisome biogenesis in fibroblasts from Zellweger patients of complementation group B in Japan (complementation group VII in the USA). One patient (PBDB-01) possessed a homozygous, inactivating mutation, a 2 bp deletion immediately upstream of the RING motif, which resulted in a frameshift, altering 65 amino acids from the normal. This implies that the C-terminal part, including the RING finger, is required for biological function of Pex10p. PEX10 cDNA derived from patient PBDB-01 was defective in peroxisome-restoring activity when expressed in patient fibroblasts. These results demonstrate that mutation in PEX10 is the genetic cause of complementation group B PBD.      

Isolation, characterization and mutation analysis of PEX13-defective Chinese hamster ovary cell mutants.

We isolated peroxisome biogenesis mutants ZP128 and ZP150 from rat PEX2 -transformed Chinese hamster ovary (CHO) cells, by the 9-(1'-pyrene)nonanol/ultraviolet method. The mutants lacked morphologically recognizable peroxisomes and showed a typical peroxisome assembly-defective phenotype such as a high sensitivity to 12-(1'-pyrene)dodecanoic acid/UV treatment. By means of PEX cDNA transfection and cell fusion, ZP128 and ZP150 were found to belong to a recently identified complementation group H. Expression of human PEX13 cDNA restored peroxisome assembly in ZP128 and ZP150. CHO cell PEX13 was isolated; its deduced sequence comprises 405 amino acids with 93% identity to human Pex13p. Mutation in PEX13 of mutant ZP150 was determined by RT-PCR: G to A transition resulted in one amino acid substitution, Ser319Asn, in one allele and truncation of a 42 amino acid sequence from Asp265 to Lys306 in another allele. Therefore, ZP128 and ZP150 are CHO cell lines with a phenotype of impaired PEX13.     

Novel PEX1 mutations and genotype-phenotype correlations in Australasian peroxisome biogenesis disorder patients

The peroxisome biogenesis disorders (PBDs) are a group of neuronal migration/neurodegenerative disorders that arise from defects in PEX genes. A major subgroup of the PBDs includes Zellweger syndrome (ZS), neonatal adrenoleukodystrophy (NALD), and infantile Refsum disease (IRD). These three disorders represent a clinical continuum with Zellweger syndrome the most severe. Mutations in the PEX1 gene, which encodes a protein of the AAA ATPase family involved in peroxisome matrix protein import, account for the genetic defect in more than half of the patients in this PBD subgroup. We report here on the results of PEX1 mutation detection in an Australasian cohort of PEX1-deficient PBD patients. This screen has identified five novel mutations, including nonsense mutations in exons 14 and 19 and single nucleotide deletions in exons 5 and 18. Significantly, the allele carrying the exon 18 frameshift mutation is present at moderately high frequency (approx. 10%) in this patient cohort. The fifth mutation is a missense mutation (R798G) that attenuates, but does not abolish PEX1 function. We have evaluated the cellular impact of these novel mutations, along with that of the two most common PEX1 mutations (c.2097-2098insT and G843D), in PBD patients by determining the levels of PEX1 mRNA, PEX1 protein, and peroxisome protein import. The findings are consistent with a close correlation between cellular phenotype, disease severity, and PEX1 genotype.     

Identification of PEX3 as the gene mutated in a Zellweger syndrome patient lacking peroxisomal remnant structures

Peroxisome biogenesis disorders, of which 13 complementation groups have been identified, are subdivided with regard to two major dysfunctions: peroxisomal matrix protein import and peroxisomal membrane synthesis. Detectable remnant membrane structures are evident only in the former. Molecular defects have been defined in 10 PEX genes, including eight related to protein import and two to membrane synthesis. We now have evidence that the human complete cDNA encoding Pex3p, a peroxisomal membrane protein (PMP) factor for the proper localization of PMPs, rescues the import of both PMP and the matrix protein in fibroblasts from a Zellweger syndrome patient of complementation group G. This patient was homozygous for a 1 base insertion in the codon for V182, which resulted in a change of codon (182-183) and introduced a termination codon (184), which inactivated PMP and matrix protein import by Pex3p. A PEX3-defective CHO mutant clone, ZPG208, was of the same complementation group as group G.     

A homozygous Gly470Ala variant in PEX6 causes severe Zellweger spectrum disorder

Zellweger spectrum disorder (ZSD) is a group of autosomal recessive disorders caused by biallelic pathogenic variants in any one of the 13 PEX genes essential for peroxisomal biogenesis. We report a cohort of nine infants who presented at birth with severe neonatal features suggestive of ZSD and found to be homozygous for a variant in PEX6 (NM_000287.4:c.1409G > C[p.Gly470Ala]). All were of Mixtec ancestry and identified by the California Newborn Screening (NBS) Program to have elevated C26:0-lysophosphatidylcholine but no reportable variants in ABCD1. The clinical and biochemical features of this cohort are described within. Gly470Ala may represent a founder variant in the Mixtec population of Central California. ZSD should be considered in patients who present at birth with severe hypotonia and enlarged fontanelles, especially in the setting of an abnormal NBS, Mixtec ancestry, or family history of infant death. There is a need to further characterize the natural history of ZSD, the Gly470Ala variant, and expand upon possible genotype–phenotype correlations.     

Autosomal dominant Zellweger spectrum disorder caused by de novo variants in PEX14 gene

PurposeZellweger spectrum disorders (ZSDs) are known as autosomal recessive disorders caused by defective peroxisome biogenesis due to bi-allelic pathogenic variants in any of at least 13 different PEX genes. Here, we report 2 unrelated patients who present with an autosomal dominant ZSD.MethodsWe performed biochemical and genetic studies in blood and skin fibroblasts of the patients and demonstrated the pathogenicity of the identified PEX14 variants by functional cell studies.ResultsWe identified 2 different single heterozygous de novo variants in the PEX14 genes of 2 patients diagnosed with ZSD. Both variants cause messenger RNA mis-splicing, leading to stable expression of similar C-terminally truncated PEX14 proteins. Functional studies indicated that the truncated PEX14 proteins lost their function in peroxisomal matrix protein import and cause increased degradation of peroxisomes, ie, pexophagy, thus exerting a dominant-negative effect on peroxisome functioning. Inhibition of pexophagy by different autophagy inhibitors or genetic knockdown of the peroxisomal autophagy receptor NBR1 resulted in restoration of peroxisomal functions in the patients’ fibroblasts.ConclusionOur finding of an autosomal dominant ZSD expands the genetic repertoire of ZSDs. Our study underscores that single heterozygous variants should not be ignored as possible genetic cause of diseases with an established autosomal recessive mode of inheritance.     

Impaired neuronal migration and endochondral ossification in Pex7 knockout mice: a model for rhizomelic chondrodysplasia punctata

Rhizomelic chondrodysplasia punctata is a human autosomal recessive disorder characterized by skeletal, eye and brain abnormalities. The disorder is caused by mutations in the PEX7 gene, which encodes the receptor for a class of peroxisomal matrix enzymes. We describe the generation and characterization of a Pex7 mouse knockout (Pex7(-/-)). Pex7(-/-) mice are born severely hypotonic and have a growth impairment. Mortality in Pex7(-/-) mice is highest in the perinatal period although some Pex7(-/-) mice survived beyond 18 months. Biochemically Pex7(-/-) mice display the abnormalities related to a Pex7 deficiency, i.e. a severe depletion of plasmalogens, impaired alpha-oxidation of phytanic acid and impaired beta-oxidation of very-long-chain fatty acids. In the intermediate zone of the developing cerebral cortex Pex7(-/-) mice have an increase in neuronal density. In vivo neuronal birthdating revealed that Pex7(-/-) mice have a delay in neuronal migration. Analysis of bone ossification in newborn Pex7(-/-) mice revealed a defect in ossification of distal bone elements of the limbs as well as parts of the skull and vertebrae. These findings demonstrate that Pex7 knockout mice provide an important model to study the role of peroxisomal functioning in the pathogenesis of the human disorder.     

Mutation analysis of PEX7 in 60 probands with rhizomelic chondrodysplasia punctata and functional correlations of genotype with phenotype

PEX7 encodes the cytosolic receptor for the set of peroxisomal matrix enzymes targeted to the organelle by the peroxisome targeting signal 2 (PTS2). Mutations in PEX7 cause rhizomelic chondrodysplasia punctata (RCDP), a distinct peroxisome biogenesis disorder. In previous work we described three novel PEX7 mutant alleles, including one, L292X, with a high frequency due to a founder effect. We have now extended our analysis to 60 RCDP probands and identified a total of 24 PEX7 alleles, accounting for 95% of the mutant PEX7 genes in our sample. Of these, 50% are L292X, 13% are IVS9+1G>C, and the remainder are mostly private. IVS9+1G>C occurs on at least three different haplotypes and thus appears to result from recurrent mutation. The phenotypic spectrum of RCDP is broader than commonly recognized and includes minimally affected individuals at the mild end of the spectrum. To relate PEX7 genotype and phenotype, we evaluated the consequence of the disease mutation on PEX7 RNA by Northern analysis and RT/PCR. We evaluated the function of the encoded Pex7 protein (Pex7p) by expressing selected alleles in fibroblasts from RCDP patients and assaying their ability to restore import of a PTS2 marker protein. We find that residual activity of mutant Pex7p and reduced amounts of normal Pex7p are associated with milder and variant phenotypes.     

PEX6 is Expressed in Photoreceptor Cilia and Mutated in Deafblindness with Enamel Dysplasia and Microcephaly

Deafblindness is part of several genetic disorders. We investigated a consanguineous Egyptian family with two siblings affected by congenital hearing loss and retinal degeneration, initially diagnosed as Usher syndrome type 1. At teenage, severe enamel dysplasia, developmental delay, and microcephaly became apparent. Genome‐wide homozygosity mapping and whole‐exome sequencing detected a homozygous missense mutation, c.1238G>T (p.Gly413Val), affecting a highly conserved residue of peroxisomal biogenesis factor 6, PEX6. Biochemical profiling of the siblings revealed abnormal and borderline plasma phytanic acid concentration, and cerebral imaging revealed white matter disease in both. We show that Pex6 localizes to the apical extensions of secretory ameloblasts and differentiated odontoblasts at early stages of dentin synthesis in mice, and to cilia of retinal photoreceptor cells. We propose PEX6, and possibly other peroxisomal genes, as candidate for the rare cooccurrence of deafblindness and enamel dysplasia. Our study for the first time links peroxisome biogenesis disorders to retinal ciliopathies.     

Genetic and clinical aspects of Zellweger spectrum patients with PEX1 mutations

Objective: To analyse the PEX1 gene, the most common cause for peroxisome biogenesis disorders (PBD), in a consecutive series of patients with Zellweger spectrum. Methods: Mutations were detected by different methods including SSCP analyses as a screening technique on the basis of genomic or cDNA, followed by direct sequencing of PCR fragments with an abnormal electrophoresis pattern. Results: 33 patients were studied. Two common mutations, c.2528G→A, G843D and c.2098_2098insT, I700YfsX42, accounted for over 80% of all abnormal PEX1 alleles, emphasising their diagnostic relevance. Most PEX1 mutations were distributed over the two AAA cassettes with the two functional protein domains, D1 and D2, and the highly conserved Walker motifs. Phenotypic severity of Zellweger spectrum in CG1 depended on the effect of the mutation on the PEX1 protein, peroxin 1. PEX1 mutations could be divided into two classes of genotype–phenotype correlation: class I mutations led to residual PEX1 protein levels and function and a milder phenotype; class II mutations almost abolished PEX1 protein levels and function, resulting in a severe phenotype. Compound heterozygote patients for a class I and class II mutation had an intermediate phenotype. Conclusions: Molecular confirmation of the clinical and biochemical diagnosis will allow the prediction of the clinical course of disease in individual PBD cases.     

Identification of a novel mutation in the coding region of the grey-lethal gene OSTM1 in human malignant infantile osteopetrosis

Autosomal recessive malignant infantile osteopetrosis (ARO) is characterized by severe osteosclerosis, pathologic fractures, hepatosplenomegaly, and pancytopenia. The pathophysiological basis is inadequate bone resorption due to osteoclast dysfunction. In the majority of cases, mutations in either of two human genes cause this fatal disorder: TCIRG1, encoding a subunit of the osteoclast H(+)-ATPase, and the voltage-gated chloride channel gene CLCN7. We excluded both genes in a small inbred family with malignant infantile osteopetrosis and undertook linkage analysis of several candidate loci that are involved in murine osteopetrosis. A region spanning more than 20 cM between the markers D6S1717 and D6S1608 on chromosome 6q21 was found to be homozygous in the affected child. This locus is syntenic to the genomic region harboring the gene for the osteopetrotic mutant mouse grey-lethal (gl). Recently, mutations in a novel gene of unknown function were described in the grey-lethal mouse and in one human patient. Mutation screening of the grey-lethal gene (OSTM1), revealed a homozygous 2-bp deletion in exon 2 (c.415_416delAG) in the affected child. No mutations could be found in six independent ARO patients who had tested negative for mutations in TCIRG1 and CLCN7. In summary, we describe the identification of a novel mutation in the coding sequence of the human grey-lethal gene, which is the second OSTM1 mutation found in human ARO, confirming the involvement of this gene in the pathogenesis of this severe bone disease.     

Anion exchanger 1 mutations associated with distal renal tubular acidosis in the Thai population

We have previously demonstrated that compound heterozygous (SAO/G701D) and homozygous (G701D/G701D) mutations of the anion exchanger 1 (AE1) gene, encoding erythroid and kidney AE1 proteins, cause autosomal recessive distal renal tubular acidosis (AR dRTA) in Thai patients. It is thus of interest to examine the prevalence of these mutations in the Thai population. The SAO and G701D mutations were examined in 844 individuals from north, northeast, central, and south Thailand. Other reported mutations including R602H, V850, and A858D were also examined in some groups of subjects. The SAO mutation was common in the southern Thai population; its heterozygote frequency was 7/206 and estimated allele frequency 1.70%. However, this mutation was not observed in populations of three other regions of Thailand. In contrast, the G701D mutation was not found in the southern population but was observed in the northern, northeastern, and central populations, with heterozygote frequencies of 1/216, 3/205, and 1/217, and estimated allele frequencies of 0.23%, 0.73%, and 0.23%, respectively. The higher allele frequency of the G701D mutation in the northeastern Thai population corresponds to our previous finding that all Thai patients with AR dRTA attributable to homozygous G701D mutation originate from this population. This suggests that the G701D allele that is observed in this region might arise in northeastern Thailand. The presence of patients with compound heterozygous SAO/G701D in southern Thailand and Malaysia and their apparently absence in northeastern Thailand indicate that the G701D allele may have migrated to the southern peninsular region where SAO is common, resulting in pathogenic allelic interaction.     

Digenic inheritance of mutations in HAMP and HFE results in different types of haemochromatosis

Haemochromatosis (HH) is a clinically and genetically heterogeneous disease caused by inappropriate iron absorption. Most HH patients are homozygous for the C282Y mutation in the HFE gene. However, penetrance of the C282Y mutation is incomplete, and other genetic factors may well affect the HH phenotype. Ferroportin and TFR2 mutations also cause HH, and two HAMP mutations have recently been reported that causes juvenile haemochromatosis (JH) in the homozygous state. Here, we report evidence for digenic inheritance of HH. We have detected two new HAMP mutations in two different families, in which there is concordance between severity of iron overload and heterozygosity for HAMP mutations when present with the HFE C282Y mutation. In family A, the proband has a JH phenotype and is heterozygous for C282Y and a novel HAMP mutation Met50del IVS2+1(-G). This is a four nucleotide ATGG deletion which causes a frameshift. The proband's unaffected mother is also heterozygous for Met50del IVS2+1(-G), but lacks the C282Y mutation and is heterozygous for the HFE H63D mutation. Met50del IVS2+1(-G) was absent from 642 control chromosomes. In family B, a second novel, less severe HAMP mutation, G71D, was identified. This was detected in the general population at an allele frequency of 0.3%. We propose that the phenotype of C282Y heterozygotes and homozygotes may be modified by heterozygosity for mutations which disrupt the function of hepcidin in iron homeostasis, with the severity of iron overload corresponding to the severity of the HAMP mutation.     

Production of MPS VII mouse (Gus(tm(hE540A x mE536A)Sly)) doubly tolerant to human and mouse beta-glucuronidase

Mucopolysaccharidosis VII (MPS VII, Sly syndrome) is an autosomal recessive lysosomal storage disease caused by beta-glucuronidase (GUS) deficiency. A naturally occurring mouse model of that disease has been very useful for studying experimental approaches to therapy. However, immune responses can complicate evaluation of the long-term benefits of enzyme replacement or gene therapy delivered to adult MPS VII mice. To make this model useful for studying the long-term effectiveness and side effects of experimental therapies delivered to adult mice, we developed a new MPS VII mouse model, which is tolerant to both human and murine GUS. To achieve this, we used homologous recombination to introduce simultaneously a human cDNA transgene expressing inactive human GUS into intron 9 of the murine Gus gene and a targeted active site mutation (E536A) into the adjacent exon 10. When the heterozygote products of germline transmission were bred to homozygosity, the homozygous mice expressed no GUS enzyme activity but expressed inactive human GUS protein highly and were tolerant to immune challenge with human enzyme. Expression of the mutant murine Gus gene was reduced to about 10% of normal levels, but the inactive murine GUS enzyme also conferred tolerance to murine GUS. This MPS VII mouse model should be useful to evaluate therapeutic responses in adult mice receiving repetitive doses of enzyme or mice receiving gene therapy as adults. Heterozygotes expressed only 9.5-26% of wild-type levels of murine GUS instead of the expected 50%, indicating a dominant-negative effect of the mutant enzyme monomers on the activity of GUS tetramers in different tissues. Corrective gene therapy in this model should provide high enough levels of expression of normal GUS monomers to overcome the dominant negative effect of mutant monomers on newly synthesized GUS tetramers in most tissues.     

Humanizing animal models: a key to autoimmune diabetes treatment

Preclinical evaluation of antibody-based immunotherapies for the treatment of type 1 diabetes (T1D) in animal models is often hampered by the fact that the human antibody drug does not cross-react with its mouse counterpart. In this issue of Science Translational Medicine, researchers describe a new mouse model that expresses the human isoform of a molecule targeted by T1D antibody therapies that are currently being tested in clinical trials--the human epsilon chain of the CD3 complex expressed on T cells. Anti-CD3 is capable of reducing insulin needs in individuals with recently diagnosed T1D; however, the precise underlying mechanisms of action and the minimal effective dose have been difficult to define. The new humanized mouse model will be instrumental in optimizing anti-CD3-based therapies and accelerating their clinical realization.     

Carrier frequency estimation of Zellweger spectrum disorder using ExAC database and bioinformatics tools

PurposeWe aimed to estimate the carrier frequency of Zellweger spectrum disorder (ZSD), a rare autosomal recessive disease, and the associated disease incidence based on data from the Exome Aggregation Consortium (ExAC) of approximately 60,000 individuals.MethodsWe obtained variants from ExAC in 13 PEX genes associated with ZSD. Potentially pathogenic missense variants were identified with computational variant analysis tools according to three stringency levels. Using variants classified as potentially pathogenic, we estimated the carrier frequency and the associated incidence for the entire ExAC population and its subpopulations. We also evaluated variants based on pathogenicity criteria for sequence variant interpretation outlined by the American College of Medical Genetics and Genomics (ACMG) and calculated the carrier frequency and incidence based on those variants.ResultsThe bioinformatically estimated incidence rate of ZSD in the ExAC population is 1 in 83,841 using our least stringent pathogenicity cutoff. Under clinical guidelines outlined by ACMG, the estimated incidence is 1 in 3,275,751 births.ConclusionWe outlined a process for estimating the ZSD disease carrier frequency and incidence in a large consortium using bioinformatics tools. Our results are close to current newborn screening estimates in New York of 1 in 90,000 births, estimated from 1.08 million screenings.