Quantitative retrospective natural history modeling for orphan drug development.

The natural history of most rare diseases is incompletely understood and usually relies on studies with low level of evidence. Consistent with the goals for future research of rare disease research set by the International Rare Diseases Research Consortium in 2017, the purpose of this paper is to review the recently developed method of quantitative retrospective natural history modeling (QUARNAM) and to illustrate its usefulness through didactically selected analyses examples in an overall population of 849 patients worldwide with seven (ultra-) rare neurogenetic disorders. A quantitative understanding of the natural history of the disease is fundamental for the development of specific interventions and counseling afflicted families. QUARNAM has a similar relationship to a published case study as a meta-analysis has to an individual published study. QUARNAM relies on sophisticated statistical analyses of published case reports focusing on four research questions: How long does it take to make the diagnosis? How long do patients live? Which factors predict disease severity (eg, genotypes, signs/symptoms, biomarkers)? Where can patients be recruited for studies? Useful statistical techniques include Kaplan-Meier estimates, cluster analysis, regression techniques, binary decisions trees, word clouds, and geographic mapping. In comparison to other natural history study methods (prospective studies or retrospective studies such as chart reviews), QUARNAM can provide fast information on hard clinical endpoints (ie, survival, diagnostic delay) with a lower effort. The choice of method for a particular drug development program may be driven by the research question and may encompass combinatory approaches.

Biallelic mutations in nucleoporin NUP88 cause lethal fetal akinesia deformation sequence.

Nucleoporins build the nuclear pore complex (NPC), which, as sole gate for nuclear-cytoplasmic exchange, is of outmost importance for normal cell function. Defects in the process of nucleocytoplasmic transport or in its machinery have been frequently described in human diseases, such as cancer and neurodegenerative disorders, but only in a few cases of developmental disorders. Here we report biallelic mutations in the nucleoporin NUP88 as a novel cause of lethal fetal akinesia deformation sequence (FADS) in two families. FADS comprises a spectrum of clinically and genetically heterogeneous disorders with congenital malformations related to impaired fetal movement. We show that genetic disruption of nup88 in zebrafish results in pleiotropic developmental defects reminiscent of those seen in affected human fetuses, including locomotor defects as well as defects at neuromuscular junctions. Phenotypic alterations become visible at distinct developmental stages, both in affected human fetuses and in zebrafish, whereas early stages of development are apparently normal. The zebrafish phenotypes caused by nup88 deficiency are rescued by expressing wild-type Nup88 but not the disease-linked mutant forms of Nup88. Furthermore, using human and mouse cell lines as well as immunohistochemistry on fetal muscle tissue, we demonstrate that NUP88 depletion affects rapsyn, a key regulator of the muscle nicotinic acetylcholine receptor at the neuromuscular junction. Together, our studies provide the first characterization of NUP88 in vertebrate development, expand our understanding of the molecular events causing FADS, and suggest that variants in NUP88 should be investigated in cases of FADS.

Clinical characteristics of 248 patients with Krabbe disease: quantitative natural history modeling based on published cases.

PURPOSE: Krabbe disease (OMIM 245200) is an orphan neurometabolic disorder caused by a deficiency of the lysosomal enzyme galactocerebrosidase (GALC). Hard clinical endpoints and biomarker-phenotype correlations are useful for future clinical trials. METHODS: We performed a quantitative analysis of published cases (N = 248) with Krabbe disease, stratified by age at disease onset: early infantile (age 0-6 months), late infantile (age 7-36 months), juvenile/adolescent (age 37-180 months), and adult onset (>180 months). Main outcome measures were age of disease onset and survival. Cerebrospinal fluid (CSF) protein concentrations were explored as a potential predictor of survival. STROBE criteria were respected. RESULTS: Median age of onset was 4 months (early infantile), 14 months (late infantile), 48 months (juvenile), and 384 months (adult). Age of disease onset and therefore disease subtype determined survival rates. CSF protein concentrations predicted age at onset and survival rates in Krabbe disease. Patients with a CSF protein content ≤61.5 mg/dl survived significantly longer than patients with CSF protein values above this threshold. CONCLUSION: We define the estimated survival in published Krabbe disease cases and demonstrate an association of CSF protein concentration with disease severity. These data inform patient care and clinical trials.

Impairment of astrocytic glutaminolysis in glutaric aciduria type I.

Glutaric aciduria type I is a rare, autosomal recessive, inherited defect of glutaryl-CoA dehydrogenase. Deficiency of this protein in L-lysine degradation leads to the characteristic accumulation of nontoxic glutarylcarnitine and neurotoxic glutaric acid (GA), glutaryl-CoA, and 3-hydroxyglutaric acid. Untreated patients develop bilateral lesions of basal ganglia resulting in a complex movement disorder with predominant dystonia in infancy and early childhood. The current pathomechanistic concept strongly focuses on imbalanced neuronal energy metabolism due to accumulating metabolites, whereas little is known about the pathomechanistic role of astrocytes, which are thought to be in constant metabolic crosstalk with neurons. We found that glutaric acid (GA) causes astrocytic cell death under starvation cell culture conditions, i.e. low glucose, without glutamine and fetal calf serum. Glutamine completely abolished GA-induced toxicity, suggesting involvement of glutaminolysis. Increasing dependence on glutaminolysis by chemical induction of hypoxia signaling-potentiated GA-induced toxicity. We further show that GA disturbs glutamine degradation by specifically inhibiting glutamate dehydrogenase. Summarizing our study shows that pathologically relevant concentrations of GA block an important step in the metabolic crosstalk between neurons and astrocytes, ultimately leading to astrocytic cell death.

Multifactorial modulation of susceptibility to l-lysine in an animal model of glutaric aciduria type I.

Glutaric aciduria type I is an inherited defect in L-lysine, L-hydroxylysine and L-tryptophan degradation caused by deficiency of glutaryl-CoA dehydrogenase (GCDH). The majority of untreated patients presents with accumulation of neurotoxic metabolites - glutaric acid (GA) and 3-hydroxyglutaric acid (3-OHGA) - and striatal injury. Gcdh(-/-) mice display elevated levels of GA and 3-OH-GA but do not spontaneously develop striatal lesions. L-lysine-enriched diets (appr. 235 mg/d) were suggested to induce a neurological phenotype similar to affected patients. In our hands 93% of mice stressed according to the published protocol remained asymptomatic. To understand the underlying mechanism, we modified their genetic background (F1 C57BL6/Jx129/SvCrl) and increased the daily oral L-lysine supply (235-433 mg). We identified three modulating factors, (1) gender, (2) genetic background, and (3) amount of L-lysine. Male mice displayed higher vulnerability and inbreeding for more than two generations as well as elevating L-lysine supply increased the diet-induced mortality rate (up to 89%). Onset of first symptoms leads to strongly reduced intake of food and, thus, L-lysine suggesting a threshold for toxic metabolite production to induce neurological disease. GA and 3-OH-GA tissue concentrations did not correlate with dietary L-lysine supply but differed between symptomatic and asymptomatic mice. Cerebral activities of glyceraldehyde 3-phosphate dehydrogenase, 2-oxoglutarate dehydrogenase complex, and aconitase were decreased. Symptomatic mice did not develop striatal lesions or intracerebral hemorrhages. We found severe spongiosis in the hippocampus of Gcdh(-/-) mice which was independent of dietary L-lysine supply. In conclusion, the L-lysine-induced pathology in Gcdh(-/-) mice depends on genetic and dietary parameters.

First Japanese case of Zellweger syndrome with a mutation in PEX14.

Zellweger syndrome, one of the peroxisome biogenesis disorders, is an autosomal recessive disease caused by mutations in PEX genes. It is characterized by severe hypotonia, failure to thrive, psychomotor retardation, liver dysfunction, and sensorineural hearing impairment. Most of the patients with this disease die before the age of 1 year. PEX14 is the 13th PEX gene responsible for peroxisome biogenesis disorders. Thus far, only two patients with PEX14 deficiency have been reported. Here, we report the first case of a Japanese patient with a PEX14 mutation who showed severe hypotonia, psychomotor retardation, demyelination, and developed rickets at the age of 5 months. An increased excretion of 3,6-epoxydicarboxylic acids leads to the diagnosis of Zellweger syndrome and a mutation analysis of PEX14 revealed a homozygous mutation of c.538C>T (p.Q180X). The patient survived for a prolonged period of time but died of liver failure at the age of 46 months.

Mutations in genes encoding the glycine cleavage system predispose to neural tube defects in mice and humans.

Neural tube defects (NTDs), including spina bifida and anencephaly, are common birth defects of the central nervous system. The complex multigenic causation of human NTDs, together with the large number of possible candidate genes, has hampered efforts to delineate their molecular basis. Function of folate one-carbon metabolism (FOCM) has been implicated as a key determinant of susceptibility to NTDs. The glycine cleavage system (GCS) is a multi-enzyme component of mitochondrial folate metabolism, and GCS-encoding genes therefore represent candidates for involvement in NTDs. To investigate this possibility, we sequenced the coding regions of the GCS genes: AMT, GCSH and GLDC in NTD patients and controls. Two unique non-synonymous changes were identified in the AMT gene that were absent from controls. We also identified a splice acceptor site mutation and five different non-synonymous variants in GLDC, which were found to significantly impair enzymatic activity and represent putative causative mutations. In order to functionally test the requirement for GCS activity in neural tube closure, we generated mice that lack GCS activity, through mutation of AMT. Homozygous Amt(-/-) mice developed NTDs at high frequency. Although these NTDs were not preventable by supplemental folic acid, there was a partial rescue by methionine. Overall, our findings suggest that loss-of-function mutations in GCS genes predispose to NTDs in mice and humans. These data highlight the importance of adequate function of mitochondrial folate metabolism in neural tube closure.

Associations of common genetic risk variants of the muscarinic acetylcholine receptor M2 with cardiac autonomic dysfunction in patients with schizophrenia.

OBJECTIVES: Decreased vagal modulation, which has consistently been observed in schizophrenic patients, might contribute to increased cardiac mortality in schizophrenia. Previously, associations between CHRM2 (Cholinergic Receptor Muscarinic 2) and cardiac autonomic features have been reported. Here, we tested for possible associations between these polymorphisms and heart rate variability in patients with schizophrenia. METHODS: A total of three single nucleotide polymorphisms (SNPs) in CHRM2 (rs73158705 A>G, rs8191992 T>A and rs2350782 T>C) that achieved significance (p < 5 * 10-8) in genome-wide association studies for cardiac autonomic features were genotyped in 88 drug-naïve patients, 61 patients receiving antipsychotic medication and 144 healthy controls. Genotypes were analysed for associations with parameters of heart rate variability and complexity, in each diagnostic group. RESULTS: We observed a significantly altered heart rate variability in unmedicated patients with identified genetic risk status in rs73158705 A>G, rs8191992 T>A and rs2350782 T>C as compared to genotype non-risk status. In patients receiving antipsychotic medication and healthy controls, these associations were not observed. DISCUSSION: We report novel candidate genetic associations with cardiac autonomic dysfunction in schizophrenia, but larger cohorts are required for replication.

Analysis of CACNA1C and KCNH2 Risk Variants on Cardiac Autonomic Function in Patients with Schizophrenia.

BACKGROUND: Cardiac autonomic dysfunction (CADF) is a major contributor to increased cardiac mortality in schizophrenia patients. The aberrant function of voltage-gated ion channels, which are widely distributed in the brain and heart, may link schizophrenia and CADF. In search of channel-encoding genes that are associated with both CADF and schizophrenia, CACNA1C and KCNH2 are promising candidates. In this study, we tested for associations between genetic findings in both genes and CADF parameters in schizophrenia patients whose heart functions were not influenced by psychopharmaceuticals. METHODS: First, we searched the literature for single-nucleotide polymorphisms (SNPs) in CACNA1C and KCNH2 that showed genome-wide significant association with schizophrenia. Subsequently, we looked for such robust associations with CADF traits at these loci. A total of 5 CACNA1C SNPs and 9 KCNH2 SNPs were found and genotyped in 77 unmedicated schizophrenia patients and 144 healthy controls. Genotype-related impacts on heart rate (HR) dynamics and QT variability indices (QTvi) were analyzed separately in patients and healthy controls. RESULTS: We observed significantly increased QTvi in unmedicated patients with CADF-associated risk in CACNA1C rs2283274 C and schizophrenia-associated risk in rs2239061 G compared to the non-risk allele in these patients. Moreover, unmedicated patients with previously identified schizophrenia risk alleles in KCNH2 rs11763131 A, rs3807373 A, rs3800779 C, rs748693 G, and 1036145 T showed increased mean HR and QTvi as compared to non-risk alleles. CONCLUSIONS: We propose a potential pleiotropic role for common variation in CACNA1C and KCNH2 associated with CADF in schizophrenia patients, independent of antipsychotic medication, that predisposes them to cardiac arrhythmias and premature death.

A genome-wide association study identifies RNF213 as the first Moyamoya disease gene.

Moyamoya disease (MMD) shows progressive cerebral angiopathy characterized by bilateral internal carotid artery stenosis and abnormal collateral vessels. Although ∼ 15% of MMD cases are familial, the MMD gene(s) remain unknown. A genome-wide association study of 785,720 single-nucleotide polymorphisms (SNPs) was performed, comparing 72 Japanese MMD patients with 45 Japanese controls and resulting in a strong association of chromosome 17q25-ter with MMD risk. This result was further confirmed by a locus-specific association study using 335 SNPs in the 17q25-ter region. A single haplotype consisting of seven SNPs at the RNF213 locus was tightly associated with MMD (P = 5.3 × 10(-10)). RNF213 encodes a really interesting new gene finger protein with an AAA ATPase domain and is abundantly expressed in spleen and leukocytes. An RNA in situ hybridization analysis of mouse tissues indicated that mature lymphocytes express higher levels of Rnf213 mRNA than their immature counterparts. Mutational analysis of RNF213 revealed a founder mutation, p.R4859K, in 95% of MMD families, 73% of non-familial MMD cases and 1.4% of controls; this mutation greatly increases the risk of MMD (P = 1.2 × 10(-43), odds ratio = 190.8, 95% confidence interval = 71.7-507.9). Three additional missense mutations were identified in the p.R4859K-negative patients. These results indicate that RNF213 is the first identified susceptibility gene for MMD.

A common variation in HCN1 is associated with heart rate variability in schizophrenia.

BACKGROUND: There is growing evidence for a shared genetic basis between schizophrenia risk and cardiovascular disease. Reduced efferent vagal activity, indexed by reduced heart rate variability (HRV), has been consistently described in patients with schizophrenia and may potentially contribute to the increased cardiovascular risk in these patients. In this study, we tested the hypothesis whether the established schizophrenia risk variant HCN1 rs16902086 (A > G) is associated with reduced HRV. METHODS: We analyzed the risk status of HCN1 rs16902086 (AG/GG vs. AA genotype) in 83 unmedicated patients with schizophrenia and 96 healthy controls and investigated genotype-related impacts on various HRV parameters. RESULTS: We observed significantly increased resting heart rates and a marked decrease of vagal modulation in our patient cohort. Strikingly, HCN1 rs16902086 (A > G) was associated with reduced HRV parameters in patients only. A trend towards more pronounced HRV deviations was observed in homozygous (GG) compared to heterozygous patients (AG). CONCLUSION: We present first evidence for a genetic risk factor that is associated with decreased vagal modulation in unmedicated patients with schizophrenia. Moreover, our findings suggest that HCN1 might be involved in reduced vagal modulation and possibly in increased cardiac mortality in schizophrenia patients. Thus, our data indicate that reduced vagal modulation might be an endophenotype of schizophrenia.

Mutation analysis of the SHOC2 gene in Noonan-like syndrome and in hematologic malignancies.

Noonan syndrome is an autosomal dominant disease characterized by dysmorphic features, webbed neck, cardiac anomalies, short stature and cryptorchidism. It shows phenotypic overlap with Costello syndrome and cardio-facio-cutaneous (CFC) syndrome. Noonan syndrome and related disorders are caused by germline mutations in genes encoding molecules in the RAS/MAPK pathway. Recently, a gain-of-function mutation in SHOC2, p.S2G, has been identified as causative for a type of Noonan-like syndrome characterized by the presence of loose anagen hair. In order to understand the contribution of SHOC2 mutations to the clinical manifestations of Noonan syndrome and related disorders, we analyzed SHOC2 in 92 patients with Noonan syndrome and related disorders who did not exhibit PTPN11, KRAS, HRAS, BRAF, MAP2K1/2, SOS1 or RAF1 mutations. We found the previously identified p.S2G mutation in eight of our patients. We developed a rapid detection system to identify the p.S2G mutation using melting curve analysis, which will be a useful tool to screen for the apparently common mutation. All the patients with the p.S2G mutation showed short stature, sparse hair and atopic skin. Six of the mutation-positive patients showed severe mental retardation and easily pluckable hair, and one showed leukocytosis. No SHOC2 mutations were identified in leukemia cells from 82 leukemia patients. These results suggest that clinical manifestations in SHOC2 mutation-positive patients partially overlap with those in patients with typical Noonan or CFC syndrome and show that easily pluckable/loose anagen hair is distinctive in SHOC2 mutation-positive patients.

Livedoid vasculopathy: does hyperhomocysteinaemia play an aetiological role?

BACKGROUND: Livedoid vasculopathy (LV) has been shown to be associated with hypercoagulability. However, relevant genetic and exogenous thrombophilic factors are not fully determined. OBJECTIVES: To evaluate the frequency of hyperhomocysteinaemia (HHCE) and genotypes of hypercoagulative factors in LV patients. MATERIAL AND METHODS: Plasma homocysteine level was measured in 42 LV patients. Polymorphism of MTHFR (677C > T and 1298A > C), PAI1 (-675 5G/4G and -844A > G), and F2 (20210G > A), and the F5 Leiden mutation, as well as biochemical parameters for hypercoagulability, were analysed. RESULTS: Of the LV patients, 62% revealed mild HHCE. Polymorphisms of MTHFR were observed in 75% and 56% and the PAI1 -675 5G/4G polymorphism in 100% and 83% of patients with and without HHCE, respectively. All LV patients with renal failure had mild HHCE. A high level of comorbidity of hypertension (99%) and diabetes type 2 (44%) were noted. CONCLUSION: HHCE seems to play a major pathogenetic role in LV. A high prevalence of further procoagulative factors might support the view that LV is a "complex disease".

Model mice for mild-form glycine encephalopathy: behavioral and biochemical characterizations and efficacy of antagonists for the glycine binding site of N-methyl D-aspartate receptor.

Glycine encephalopathy (GE) is caused by an inherited deficiency of the glycine cleavage system (GCS) and characterized by accumulation of glycine in body fluids and various neurologic symptoms. Coma and convulsions develop in neonates in typical GE while psychomotor retardation and behavioral abnormalities in infancy and childhood are observed in mild GE. Recently, we have established a transgenic mouse line (low-GCS) with reduced GCS activity (29% of wild-type (WT) C57BL/6) and accumulation of glycine in the brain (Stroke, 2007; 38:2157). The purpose of the present study is to characterize behavioral features of the low-GCS mouse as a model of mild GE. Two other transgenic mouse lines were also analyzed: high-GCS mice with elevated GCS activity and low-GCS-2 mice with reduced GCS activity. As compared with controls, low-GCS mice manifested increased seizure susceptibility, aggressiveness and anxiety-like activity, which resembled abnormal behaviors reported in mild GE, whereas high-GCS mice were less sensitive to seizures, hypoactive and less anxious. Antagonists for the glycine-binding site of the N-methyl-D-aspartate receptor significantly ameliorated elevated locomotor activity and seizure susceptibility in the low-GCS mice. Our results suggest the usefulness of low-GCS mice as a mouse model for mild GE and a novel therapeutic strategy.