A Recurrent Missense Variant in AP2M1 Impairs Clathrin-Mediated Endocytosis and Causes Developmental and Epileptic Encephalopathy

The mechanism of mitochondrial autophagy regulating Clathrin-mediated endocytosis in epilepsy

OBJECTIVE

The objective of this study is to determine whether inhibition of mitophagy affects seizures through Clathrin-mediated endocytosis (CME).

METHODS

Pentylenetetrazol (PTZ) was intraperitoneally injected daily to establish a chronic PTZ-kindled seizure. The Western blot (WB) was used to compare the differences in Parkin protein expression between the epilepsy group and the control group. Immunofluorescence was used to detect the expression of MitoTracker and LysoTracker. Transferrin-Alexa488 (Tf-A488) was injected into the hippocampus of mice. We evaluated the effect of 3-methyladenine (3-MA) on epilepsy behavior through observation in PTZ-kindled models.

RESULTS

The methylated derivative of adenine, known as 3-MA, has been extensively utilized in the field of autophagy research. The transferrin protein is internalized from the extracellular environment into the intracellular space via the CME pathway. Tf-A488 uses a fluorescent marker to track CME. Western blot showed that the expression of Parkin was significantly increased in the PTZ-kindled model (p < 0.05), while 3-MA could reduce the expression (p < 0.05). The fluorescence uptake of MitoTracker and LysoTracker was increased in the primary cultured neurons induced by magnesium-free extracellular fluid (p < 0.05); the fluorescence uptake of Tf-A488 was significantly decreased in the 3-MA group compared with the control group (p < 0.05). Following hippocampal injection of Tf-A488, both the epilepsy group and the 3-MA group exhibited decreased fluorescence uptake, with a more pronounced effect observed in the 3-MA group. Inhibition of mitophagy by 3-MA from day 3 to day 9 progressively exacerbated seizure severity and shortened latency.

SIGNIFICANCE

It is speculated that the aggravation of seizures by 3-MA may be related to the failure to remove damaged mitochondria in time and effectively after inhibiting mitochondrial autophagy, affecting the vesicle endocytosis function of CME and increasing the susceptibility to epilepsy.

SUMMARY

Abnormal mitophagy was observed in a chronic pentylenetetrazol-induced seizure model and a Mg2+-free-induced spontaneous recurrent epileptiform discharge model. A fluorescent transferrin marker was utilized to track clathrin-mediated endocytosis. Using an autophagy inhibitor (3-methyladenine) on primary cultured neurons, we discovered that inhibition of autophagy led to a reduction in fluorescent transferrin uptake, while impairing clathrin-mediated endocytosis function mediated by mitophagy. Finally, we examined the effects of 3-methyladenine in an animal model of seizures showing that it exacerbated seizure severity. Ultimately, this study provides insights into potential mechanisms through which mitophagy regulates clathrin-mediated endocytosis in epilepsy.

AP2A2 mutation and defective endocytosis in a Malian family with hereditary spastic paraplegia

Hereditary spastic paraplegia (HSP) comprises a large group of neurogenetic disorders characterized by progressive lower extremity spasticity. Neurological evaluation and genetic testing were completed in a Malian family with early-onset HSP. Three children with unaffected consanguineous parents presented with symptoms consistent with childhood-onset complicated HSP. Neurological evaluation found lower limb weakness, spasticity, dysarthria, seizures, and intellectual disability. Brain MRI showed corpus callosum thinning with cortical and spinal cord atrophy, and an EEG detected slow background in the index patient. Whole exome sequencing identified a homozygous missense variant in the adaptor protein (AP) complex 2 alpha-2 subunit (AP2A2) gene. Western blot analysis showed reduced levels of AP2A2 in patient-iPSC derived neuronal cells. Endocytosis of transferrin receptor (TfR) was decreased in patient-derived neurons. In addition, we observed increased axon initial segment length in patient-derived neurons. Xenopus tropicalis tadpoles with ap2a2 knockout showed cerebral edema and progressive seizures. Immunoprecipitation of the mutant human AP-2-appendage alpha-C construct showed defective binding to accessory proteins. We report AP2A2 as a novel genetic entity associated with HSP and provide functional data in patient-derived neuron cells and a frog model. These findings expand our understanding of the mechanism of HSP and improve the genetic diagnosis of this condition.

Advances in big data and omics: Paving the way for discovery in childhood epilepsies

The insights gained from big data and omics approaches have transformed the field of childhood genetic epilepsy. With an increasing number of individuals receiving genetic testing for seizures, we are provided with an opportunity to identify clinically relevant subgroups and extract meaningful observations from this large-scale clinical data. However, the volume of data from electronic medical records and omics (e.g., genomics, transcriptomics) is so vast that standardized methods, such as the Human Phenotype Ontology, are necessary for reliable and comprehensive characterization. Here, we explore the integration of clinical and omics data, highlighting how these approaches pave the way for discovery in childhood epilepsies.

Generation of induced pluripotent stem cell line (UCSFi001-A-77) carrying a biallelic frameshift variant in exon 4 of SGIP1 through CRISPR/Cas9

SGIP1 encodes a protein Src homology 3-domain growth factor receptor-bound 2-like endophilin interacting protein 1. It is involved in the regulation of clathrin-mediated endocytosis along with having a role in energy homeostasis in neuronal systems. We generated an isogenic human induced pluripotent stem cell (iPSC) line with a biallelic frameshift variant in SGIP1. This exon has been shown to be subject to alternative splicing, leading to an isoform lacking 24 amino acids that are present in the longest SGIP isoform. The newly generated iPSC line will be helpful to dissect the differential properties of the two SGIP isoforms.

Multisystem disorder associated with a pathogenic variant in CLCN7 in the absence of osteopetrosis

BACKGROUND

We clinically and genetically evaluated a Taiwanese boy presenting with developmental delay, organomegaly, hypogammaglobulinemia and hypopigmentation without osteopetrosis. Whole-exome sequencing revealed a de novo gain-of-function variant, p.Tyr715Cys, in the C-terminal domain of ClC-7 encoded by CLCN7.

METHODS

Nicoli et al. (2019) assessed the functional impact of p.Tyr715Cys by heterologous expression in Xenopus oocytes and evaluating resulting currents.

RESULTS

The variant led to increased outward currents, indicating it underlies the patient's phenotype of lysosomal hyperacidity, storage defects and vacuolization. This demonstrates the crucial physiological role of ClC-7 antiporter activity in maintaining appropriate lysosomal pH.

CONCLUSION

Elucidating mechanisms by which CLCN7 variants lead to lysosomal dysfunction will advance understanding of genotype-phenotype correlations. Identifying modifier genes and compensatory pathways may reveal therapeutic targets. Ongoing functional characterization of variants along with longitudinal clinical evaluations will continue advancing knowledge of ClC-7's critical roles and disease mechanisms resulting from its dysfunction. Expanded cohort studies are warranted to delineate the full spectrum of associated phenotypes.

CCDC32 stabilizes clathrin-coated pits and drives their invagination

Clathrin-mediated endocytosis (CME) is essential for maintaining cellular homeostasis. Previous studies have reported more than 50 CME accessory proteins; however, the mechanism driving the invagination of clathrin-coated pits (CCPs) remains elusive. Quantitative live cell imaging reveals that CCDC32, a poorly characterized endocytic accessory protein, regulates CCP stabilization and is required for efficient CCP invagination. CCDC32 interacts with the α-appendage domain (AD) of AP2 via its coiled-coil domain to exert this function. Furthermore, we showed that the clinically observed nonsense mutations in CCDC32, which result in the development of cardio-facio-neuro-developmental syndrome (CFNDS), inhibit CME by abolishing CCDC32-AP2 interactions. Overall, our data demonstrates the function and molecular mechanism of a novel endocytic accessory protein, CCDC32, in CME regulation.

Significance Statement

Clathrin-mediated endocytosis (CME) happens via the initiation, stabilization, and invagination of clathrin-coated pits (CCPs). In this study, we used a combination of quantitative live cell imaging, ultrastructure electron microscopy and biochemical experiments to show that CCDC32, a poorly studied and functional ambiguous protein, acts as an important endocytic accessory protein that regulates CCP stabilization and invagination. Specifically, CCDC32 exerts this function via its interactions with AP2, and the coiled-coil domain of CCDC32 and the α-appendage domain (AD) of AP2 are essential in mediating CCDC32-AP2 interactions. Importantly, we demonstrate that clinically observed loss-of-function mutations in CCDC32 lose AP2 interaction capacity and inhibit CME, resulting in the development of cardio-facio-neuro-developmental syndrome (CFNDS).

A Longitudinal Exploration of CACNA1A -related Hemiplegic Migraine in Children

Introduction

Since the initial description of CACNA1A- related hemiplegic migraine (HM), the phenotypic spectrum has expanded from mild episodes in neurotypical individuals to potentially life-threatening events frequently seen in individuals with developmental and epileptic encephalopathies. However, the overall longitudinal course throughout childhood remains unknown.

Methods

We analyzed HM and seizure history in individuals with CACNA1A -related HM, delineating frequency and severity of events in monthly increments through a standardized approach. Combining these data with medication prescription information, we assessed the response of HM to different agents.

Results

Our cohort involved 15 individuals between 3 and 29 years (163 patient years) and included 11 unique and two recurrent variants (p.R1349Q and p.V1393M; both n= 2). The age of first confirmed HM ranged from 14 months to 13 years (average 3 years). 25% of all HM events were severe (lasting >3 days) and 73% of individuals had at least 1 severe occurrence. Spacing of HM events ranged from 1 month to 14 years and changes in HM severity over time of showed increases or decreases of >2 severity levels in 12/122 events. Eight individuals had epilepsy, but severity of epilepsy did not correlate with frequency and severity of HM events. While levetiracetam ( n= 6) and acetazolamide ( n= 5) were the most frequently used medications, they did not show efficacy in HM prevention or HM severity reduction. However, verapamil ( n= 3) showed efficacy in preventing HM episodes (OR 2.68, CI 1.39-5.67).

Significance

The longitudinal course of CACNA1A -related HM lacks recognizable patterns for timing and severity of HM events or correlation with seizure patterns. Our data underscores the unpredictability of CACNA1A -related HM, highlighting the need for close surveillance for reoccurring HM events even in individuals with symptom-free periods.

Key points

24% of hemiplegic migraines (HM) in CACNA1A- related disorders are severe, involving cerebral edema and greater than 4 days to recover Timing and severity of HM are unpredictable, with large changes in severity between events, and age of onset ranging from 1-13 yearsEpilepsy occurred in 53% of individuals, with neither the timing nor severity of seizures correlated with HM.

Neurodevelopmental and synaptic defects in DNAJC6 parkinsonism, amenable to gene therapy

DNAJC6 encodes auxilin, a co-chaperone protein involved in clathrin-mediated endocytosis (CME) at the presynaptic terminal. Biallelic mutations in DNAJC6 cause a complex, early-onset neurodegenerative disorder characterized by rapidly progressive parkinsonism-dystonia in childhood. The disease is commonly associated with additional neurodevelopmental, neurological and neuropsychiatric features. Currently, there are no disease-modifying treatments for this condition, resulting in significant morbidity and risk of premature mortality. To investigate the underlying disease mechanisms in childhood-onset DNAJC6 parkinsonism, we generated induced pluripotent stem cells (iPSC) from three patients harbouring pathogenic loss-of-function DNAJC6 mutations and subsequently developed a midbrain dopaminergic neuronal model of disease. When compared to age-matched and CRISPR-corrected isogenic controls, the neuronal cell model revealed disease-specific auxilin deficiency as well as disturbance of synaptic vesicle recycling and homeostasis. We also observed neurodevelopmental dysregulation affecting ventral midbrain patterning and neuronal maturation. To explore the feasibility of a viral vector-mediated gene therapy approach, iPSC-derived neuronal cultures were treated with lentiviral DNAJC6 gene transfer, which restored auxilin expression and rescued CME. Our patient-derived neuronal model provides deeper insights into the molecular mechanisms of auxilin deficiency as well as a robust platform for the development of targeted precision therapy approaches.

Generation of induced pluripotent stem cell lines from two unrelated patients affected by intellectual disability carrying homozygous variants in SGIP1

Intellectual disability (ID) is a diverse neurodevelopmental condition and almost half of the cases have a genetic etiology. SGIP1 acts as an endocytic protein that influences the signaling of receptors in neuronal systems related to energy homeostasis through its interaction with endophilins. This study focuses on the generation and characterization of induced pluripotent stem cells (iPSC) from two unrelated patients due to a frameshift variant (c.764dupA, NM_032291.4) and a splice donor site variant (c.74 + 1G > A, NM_032291.4) in the SGIP1 gene.

The clinical and genetic spectrum of paediatric speech and language disorders in 52,143 individuals

Speech and language disorders are known to have a substantial genetic contribution. Although frequently examined as components of other conditions, research on the genetic basis of linguistic differences as separate phenotypic subgroups has been limited so far. Here, we performed an in-depth characterization of speech and language disorders in 52,143 individuals, reconstructing clinical histories using a large-scale data mining approach of the Electronic Medical Records (EMR) from an entire large paediatric healthcare network. The reported frequency of these disorders was the highest between 2 and 5 years old and spanned a spectrum of twenty-six broad speech and language diagnoses. We used Natural Language Processing to assess to which degree clinical diagnosis in full-text notes were reflected in ICD-10 diagnosis codes. We found that aphasia and speech apraxia could be easily retrieved through ICD-10 diagnosis codes, while stuttering as a speech phenotype was only coded in 12% of individuals through appropriate ICD-10 codes. We found significant comorbidity of speech and language disorders in neurodevelopmental conditions (30.31%) and to a lesser degree with epilepsies (6.07%) and movement disorders (2.05%). The most common genetic disorders retrievable in our EMR analysis were STXBP1 (n=21), PTEN (n=20), and CACNA1A (n=18). When assessing associations of genetic diagnoses with specific linguistic phenotypes, we observed associations of STXBP1 and aphasia (P=8.57 × 10-7, CI=18.62-130.39) and MYO7A with speech and language development delay due to hearing loss (P=1.24 × 10-5, CI=17.46-Inf). Finally, in a sub-cohort of 726 individuals with whole exome sequencing data, we identified an enrichment of rare variants in synaptic protein and neuronal receptor pathways and associations of UQCRC1 with expressive aphasia and WASHC4 with abnormality of speech or vocalization. In summary, our study outlines the landscape of paediatric speech and language disorders, confirming the phenotypic complexity of linguistic traits and novel genotype-phenotype associations. Subgroups of paediatric speech and language disorders differ significantly with respect to the composition of monogenic aetiologies.

Spatiotemporal expression of AP-2/myosin Ⅵ in mouse cochlear IHCs and correlation with auditory function

BACKGROUND

Recycling of synaptic vesicles plays an important role in vesicle pool replenishment, neurotransmitter release and synaptic plasticity. Clathrin-mediated endocytosis (CME) is considered to be the main mechanism for synaptic vesicle replenishment. AP-2 (adaptor-related protein complex 2) and myosin Ⅵ are known as key proteins that regulate the structure and dynamics of CME.

OBJECTIVE

This study aims to reveal the spatiotemporal expression of AP-2/myosin Ⅵ in inner hair cells (IHCs) of the mouse cochlea and its correlation with auditory function.

MATERIAL AND METHODS

Immunofluorescence was used to detect the localization and expression of AP-2 and myosin Ⅵ in cochlear hair cells (HCs) of CBA/CaJ mice of various ages. qRT-PCR was used to verify the differential expression of AP-2 and myosin Ⅵ mRNA in the mouse cochlea, and ABR tests were administered to mice of various ages. A preliminary analysis of the correlation between AP-2/myosin Ⅵ levels and auditory function was conducted.

RESULTS

AP-2 was located in the cytoplasmic region of IHCs and was mainly expressed in the basal region of IHCs and the area near ribbon synapses, while myosin Ⅵ was expressed in the cytoplasmic region of IHCs and OHCs. Furthermore, AP-2 and myosin Ⅵ were not significant detected in the cochleae of P7 mice; the expression level reached a peak at P35 and then decreased significantly with age. The expression patterns and expression levels of AP-2 and myosin Ⅵ in the cochleae of the mice were consistent with the development of the auditory system.

CONCLUSIONS AND SIGNIFICANCE

AP-2 and myosin Ⅵ protein expression may differ in mice of different ages, and this variation probably leads to a difference in the efficiency in CME; it may also cause a defect in IHC function.

High-throughput screening of genetic and cellular drivers of syncytium formation induced by the spike protein of SARS-CoV-2

Mapping mutations and discovering cellular determinants that cause the spike protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) to induce infected cells to form syncytia would facilitate the development of strategies for blocking the formation of such cell-cell fusion. Here we describe high-throughput screening methods based on droplet microfluidics and the size-exclusion selection of syncytia, coupled with large-scale mutagenesis and genome-wide knockout screening via clustered regularly interspaced short palindromic repeats (CRISPR), for the large-scale identification of determinants of cell-cell fusion. We used the methods to perform deep mutational scans in spike-presenting cells to pinpoint mutable syncytium-enhancing substitutions in two regions of the spike protein (the fusion peptide proximal region and the furin-cleavage site). We also used a genome-wide CRISPR screen in cells expressing the receptor angiotensin-converting enzyme 2 to identify inhibitors of clathrin-mediated endocytosis that impede syncytium formation, which we validated in hamsters infected with SARS-CoV-2. Finding genetic and cellular determinants of the formation of syncytia may reveal insights into the physiological and pathological consequences of cell-cell fusion.

Hsc70 phosphorylation patterns and calmodulin regulate AP2 Clathrin-Coated-Vesicle life span for cell adhesion protein transport

AP2 forms AP2 CCV with clathrin and over 60 additional coat proteins. Due to this complexity, we have a limited understanding of CCV life cycle regulation. Synapses contain canonical AP2 CCV, canCCV, and more stable, thereby longer lived, AP2 CCV. The more stable AP2 CCV can be distinguished from canCCV due to the stable binding of Hsc70 to clathrin. The AP1/σ1B complex knockout leads to impaired synaptic vesicle recycling and altered endosomal protein sorting. This causes as a secondary phenotype the twofold upregulation of endocytosis by canCCV and by more stable AP2 CCV. These stable CCV are more stabilized than their wt counterpart, hence stCCV. They have less of the uncoating proteins synaptojanin1 and Hsc70, and more of the coat stabilizing AAK1. Hsc70 clathrin dissociation activity is regulated by complex phosphorylation patterns. Two major groups of hyper- and of hypo-phosphorylated Hsc70 proteins are formed. The latter are enriched in wt stable CCV and stabilized stCCV. Hsc70 T265 phosphorylation regulates binding of CaM/Ca2+. CaM/Ca2+ binding to the T265 domain blocks Hsc70 homodimerization and its concentration in stCCV required for clathrin disassembly. Kinases DYRK1A and CaMK-IIδ can phosphorylate T265 preventing CaM/Ca2+ binding. Their and the levels of STK38L and STK39/Cab39, which are able to phosphorylate additional Hsc70 residues are reduced in stCCV. The stCCV pathway sorts specifically the cell adhesion proteins CHL1 and Neurocan, supporting our model of that the stCCV pathway fulfills specific functions in synaptic plasticity.

Alterations in Synaptic Connectivity and Synaptic Transmission in Alzheimer's Disease with High Physical Activity

Background

Alzheimer's disease (AD) is a progressive neurodegeneration disease. Physical activity is one of the most promising modifiable lifestyles that can be effective in slowing down the progression of AD at an early stage.

Objective

Explore the molecular processes impaired in AD that were conversely preserved and enhanced by physical activity.

Methods

Integrated transcriptomic analyses were performed in datasets that contain AD patients and elders with different degrees of physical activity. The changes of the hub genes were validated through analyzing another two datasets. The expression of the hub genes was further detected in the hippocampus and cortexes of APP/PS1 transgenic mice with or without physical activity by Quantitative polymerase chain reaction (qPCR).

Results

Cross-comparison highlighted 195 DEGs displaying opposed regulation patterns between AD and high physical activity (HPA). The common DEGs were predominantly involved in synaptic vesicle recycling and synaptic transmission, largely downregulated in AD patients but upregulated in the elders with HPA. Two key modules and four hub genes that were related to synaptic vesicle turnover were obtained from the PPI network. The expression of these hub genes (SYT1, SYT4, SH3GL2, and AP2M1) was significantly decreased in AD transgenic mice and was reversed by HPA training.

Conclusions

HPA may reverse AD pathology by upregulating a range of synaptic vesicle transport related proteins which might improve the efficiency of synaptic vesicle turnover and facilitate inter-neuronal information transfer. The study provides novel insights into the mechanisms underlining the protective effects of HPA on AD.

Orchestrating vesicular and nonvesicular membrane dynamics by intrinsically disordered proteins

Compartmentalization by membranes is a common feature of eukaryotic cells and serves to spatiotemporally confine biochemical reactions to control physiology. Membrane-bound organelles such as the endoplasmic reticulum (ER), the Golgi complex, endosomes and lysosomes, and the plasma membrane, continuously exchange material via vesicular carriers. In addition to vesicular trafficking entailing budding, fission, and fusion processes, organelles can form membrane contact sites (MCSs) that enable the nonvesicular exchange of lipids, ions, and metabolites, or the secretion of neurotransmitters via subsequent membrane fusion. Recent data suggest that biomolecule and information transfer via vesicular carriers and via MCSs share common organizational principles and are often mediated by proteins with intrinsically disordered regions (IDRs). Intrinsically disordered proteins (IDPs) can assemble via low-affinity, multivalent interactions to facilitate membrane tethering, deformation, fission, or fusion. Here, we review our current understanding of how IDPs drive the formation of multivalent protein assemblies and protein condensates to orchestrate vesicular and nonvesicular transport with a special focus on presynaptic neurotransmission. We further discuss how dysfunction of IDPs causes disease and outline perspectives for future research.

α-Synuclein induces deficiency in clathrin-mediated endocytosis through inhibiting synaptojanin1 expression

Parkinson's disease (PD) is an age-related chronic neurological disorder, mainly characterized by the pathological feature of α-synuclein (α-syn) aggregation, with the exact disease pathogenesis unclear. During the onset and progression of PD, synaptic dysfunction, including dysregulation of axonal transport, impaired exocytosis, and endocytosis are identified as crucial events of PD pathogenesis. It has been reported that over-expression of α-syn impairs clathrin-mediated endocytosis (CME) in the synapses. However, the underlying mechanisms still needs to be explored. In this study, we investigated the molecular events underlying the synaptic dysfunction caused by over-expression of wild-type human α-syn and its mutant form, involving series of proteins participating in CME. We found that excessive human α-syn causes impaired fission and uncoating of clathrin-coated vesicles during synaptic vesicle recycling, leading to reduced clustering of synaptic vesicles near the active zone and increased size of plasma membrane and number of endocytic intermediates. Furthermore, over-expressed human α-syn induced changes of CME-associated proteins, among which synaptojanin1 (SYNJ1) showed significant reduction in various brain regions. Over-expression of SYNJ1 in primary hippocampal neurons from α-syn transgenic mice recovered the synaptic vesicle density, clustering and endocytosis. Using fluorescence-conjugated transferrin, we demonstrated that SYNJ1 re-boosted the CME activity by restoring the phosphatidylinositol-4,5-bisphosphate homeostasis. Our data suggested that over-expression of α-syn disrupts synaptic function through interfering with vesicle recycling, which could be alleviated by re-availing of SYNJ1. Our study unrevealed a molecular mechanism of the synaptic dysfunction in PD pathogenesis and provided a potential therapeutic target for treating PD.

Truncated DNM1 variant underlines developmental delay and epileptic encephalopathy

Background

Developmental and epileptic encephalopathies (DEEs) signify a group of heterogeneous neurodevelopmental disorder associated with early-onset seizures accompanied by developmental delay, hypotonia, mild to severe intellectual disability, and developmental regression. Variants in the DNM1 gene have been associated with autosomal dominant DEE type 31A and autosomal recessive DEE type 31B.

Methods

In the current study, a consanguineous Pakistani family consisting of a proband (IV-2) was clinically evaluated and genetically analyzed manifesting in severe neurodevelopmental phenotypes. WES followed by Sanger sequencing was performed to identify the disease-causing variant. Furthermore, 3D protein modeling and dynamic simulation of wild-type and mutant proteins along with reverse transcriptase (RT)-based mRNA expression were checked using standard methods.

Results

Data analysis of WES revealed a novel homozygous non-sense variant (c.1402G>T; p. Glu468*) in exon 11 of the DNM1 gene that was predicted as pathogenic class I. Variants in the DNM1 gene have been associated with DEE types 31A and B. Different bioinformatics prediction tools and American College of Medical Genetics guidelines were used to verify the identified variant. Sanger sequencing was used to validate the disease-causing variant. Our approach validated the pathogenesis of the variant as a cause of heterogeneous neurodevelopmental disorders. In addition, 3D protein modeling showed that the mutant protein would lose most of the amino acids and might not perform the proper function if the surveillance non-sense-mediated decay mechanism was skipped. Molecular dynamics analysis showed varied trajectories of wild-type and mutant DNM1 proteins in terms of root mean square deviation, root mean square fluctuation and radius of gyration. Similarly, RT-qPCR revealed a substantial reduction of the DNM1 gene in the index patient.

Conclusion

Our finding further confirms the association of homozygous, loss-of-function variants in DNM1 associated with DEE type 31B. The study expands the genotypic and phenotypic spectrum of pathogenic DNM1 variants related to DNM1-associated pathogenesis.

Reversal of cell, circuit and seizure phenotypes in a mouse model of DNM1 epileptic encephalopathy

Dynamin-1 is a large GTPase with an obligatory role in synaptic vesicle endocytosis at mammalian nerve terminals. Heterozygous missense mutations in the dynamin-1 gene (DNM1) cause a novel form of epileptic encephalopathy, with pathogenic mutations clustering within regions required for its essential GTPase activity. We reveal the most prevalent pathogenic DNM1 mutation, R237W, disrupts dynamin-1 enzyme activity and endocytosis when overexpressed in central neurons. To determine how this mutation impacted cell, circuit and behavioural function, we generated a mouse carrying the R237W mutation. Neurons from heterozygous mice display dysfunctional endocytosis, in addition to altered excitatory neurotransmission and seizure-like phenotypes. Importantly, these phenotypes are corrected at the cell, circuit and in vivo level by the drug, BMS-204352, which accelerates endocytosis. Here, we demonstrate a credible link between dysfunctional endocytosis and epileptic encephalopathy, and importantly reveal that synaptic vesicle recycling may be a viable therapeutic target for monogenic intractable epilepsies.

Genetic architecture of brain age and its causal relations with brain and mental disorders

The difference between chronological age and the apparent age of the brain estimated from brain imaging data-the brain age gap (BAG)-is widely considered a general indicator of brain health. Converging evidence supports that BAG is sensitive to an array of genetic and nongenetic traits and diseases, yet few studies have examined the genetic architecture and its corresponding causal relationships with common brain disorders. Here, we estimate BAG using state-of-the-art neural networks trained on brain scans from 53,542 individuals (age range 3-95 years). A genome-wide association analysis across 28,104 individuals (40-84 years) from the UK Biobank revealed eight independent genomic regions significantly associated with BAG (p < 5 × 10-8) implicating neurological, metabolic, and immunological pathways - among which seven are novel. No significant genetic correlations or causal relationships with BAG were found for Parkinson's disease, major depressive disorder, or schizophrenia, but two-sample Mendelian randomization indicated a causal influence of AD (p = 7.9 × 10-4) and bipolar disorder (p = 1.35 × 10-2) on BAG. These results emphasize the polygenic architecture of brain age and provide insights into the causal relationship between selected neurological and neuropsychiatric disorders and BAG.

Individualised human phenotype ontology gene panels improve clinical whole exome and genome sequencing analytical efficacy in a cohort of developmental and epileptic encephalopathies

BACKGROUND

The majority of genetic epilepsies remain unsolved in terms of specific genotype. Phenotype-based genomic analyses have shown potential to strengthen genomic analysis in various ways, including improving analytical efficacy.

METHODS

We have tested a standardised phenotyping method termed 'Phenomodels' for integrating deep-phenotyping information with our in-house developed clinical whole exome/genome sequencing analytical pipeline. Phenomodels includes a user-friendly epilepsy phenotyping template and an objective measure for selecting which template terms to include in individualised Human Phenotype Ontology (HPO) gene panels. In a pilot study of 38 previously solved cases of developmental and epileptic encephalopathies, we compared the sensitivity and specificity of the individualised HPO gene panels with the clinical epilepsy gene panel.

RESULTS

The Phenomodels template showed high sensitivity for capturing relevant phenotypic information, where 37/38 individuals' HPO gene panels included the causative gene. The HPO gene panels also had far fewer variants to assess than the epilepsy gene panel.

CONCLUSION

We have demonstrated a viable approach for incorporating standardised phenotype information into clinical genomic analyses, which may enable more efficient analysis.

The SphK1/S1P Axis Regulates Synaptic Vesicle Endocytosis via TRPC5 Channels

Sphingosine-1-phosphate (S1P), a bioactive sphingolipid concentrated in the brain, is essential for normal brain functions, such as learning and memory and feeding behaviors. Sphingosine kinase 1 (SphK1), the primary kinase responsible for S1P production in the brain, is abundant within presynaptic terminals, indicating a potential role of the SphK1/S1P axis in presynaptic physiology. Altered S1P levels have been highlighted in many neurologic diseases with endocytic malfunctions. However, it remains unknown whether the SphK1/S1P axis may regulate synaptic vesicle endocytosis in neurons. The present study evaluates potential functions of the SphK1/S1P axis in synaptic vesicle endocytosis by determining effects of a dominant negative catalytically inactive SphK1. Our data for the first time identify a critical role of the SphK1/S1P axis in endocytosis in both neuroendocrine chromaffin cells and neurons from mice of both sexes. Furthermore, our Ca2+ imaging data indicate that the SphK1/S1P axis may be important for presynaptic Ca2+ increases during prolonged stimulations by regulating the Ca2+ permeable TRPC5 channels, which per se regulate synaptic vesicle endocytosis. Collectively, our data point out a critical role of the regulation of TRPC5 by the SphK1/S1P axis in synaptic vesicle endocytosis.SIGNIFICANCE STATEMENT Sphingosine kinase 1 (SphK1), the primary kinase responsible for brain sphingosine-1-phosphate (S1P) production, is abundant within presynaptic terminals. Altered SphK1/S1P metabolisms has been highlighted in many neurologic disorders with defective synaptic vesicle endocytosis. However, whether the SphK1/S1P axis may regulate synaptic vesicle endocytosis is unknown. Here, we identify that the SphK1/S1P axis regulates the kinetics of synaptic vesicle endocytosis in neurons, in addition to controlling fission-pore duration during single vesicle endocytosis in neuroendocrine chromaffin cells. The regulation of the SphK1/S1P axis in synaptic vesicle endocytosis is specific since it has a distinguished signaling pathway, which involves regulation of Ca2+ influx via TRPC5 channels. This discovery may provide novel mechanistic implications for the SphK1/S1P axis in brain functions under physiological and pathologic conditions.

Enriching representation learning using 53 million patient notes through human phenotype ontology embedding

The Human Phenotype Ontology (HPO) is a dictionary of >15,000 clinical phenotypic terms with defined semantic relationships, developed to standardize phenotypic analysis. Over the last decade, the HPO has been used to accelerate the implementation of precision medicine into clinical practice. In addition, recent research in representation learning, specifically in graph embedding, has led to notable progress in automated prediction via learned features. Here, we present a novel approach to phenotype representation by incorporating phenotypic frequencies based on 53 million full-text health care notes from >1.5 million individuals. We demonstrate the efficacy of our proposed phenotype embedding technique by comparing our work to existing phenotypic similarity-measuring methods. Using phenotype frequencies in our embedding technique, we are able to identify phenotypic similarities that surpass current computational models. Furthermore, our embedding technique exhibits a high degree of agreement with domain experts' judgment. By transforming complex and multidimensional phenotypes from the HPO format into vectors, our proposed method enables efficient representation of these phenotypes for downstream tasks that require deep phenotyping. This is demonstrated in a patient similarity analysis and can further be applied to disease trajectory and risk prediction.

High-Resolution and Multidimensional Phenotypes Can Complement Genomics Data to Diagnose Diseases in the Neonatal Population

Advances in genomic medicine have greatly improved our understanding of human diseases. However, phenome is not well understood. High-resolution and multidimensional phenotypes have shed light on the mechanisms underlying neonatal diseases in greater details and have the potential to optimize clinical strategies. In this review, we first highlight the value of analyzing traditional phenotypes using a data science approach in the neonatal population. We then discuss recent research on high-resolution, multidimensional, and structured phenotypes in neonatal critical diseases. Finally, we briefly introduce current technologies available for the analysis of multidimensional data and the value that can be provided by integrating these data into clinical practice. In summary, a time series of multidimensional phenome can improve our understanding of disease mechanisms and diagnostic decision-making, stratify patients, and provide clinicians with optimized strategies for therapeutic intervention; however, the available technologies for collecting multidimensional data and the best platform for connecting multiple modalities should be considered.

The current landscape of epilepsy genetics: where are we, and where are we going?

PURPOSE OF REVIEW

In this review, we aim to analyse the progress in understanding the genetic basis of the epilepsies, as well as ongoing efforts to define the increasingly diverse and novel presentations, phenotypes and divergences from the expected that have continually characterized the field.

RECENT FINDINGS

A genetic workup is now considered to be standard of care for individuals with an unexplained epilepsy, due to mounting evidence that genetic diagnoses significantly influence treatment choices, prognostication, community support, and increasingly, access to clinical trials. As more individuals with epilepsy are tested, novel presentations of known epilepsy genes are being discovered, and more individuals with self-limited epilepsy are able to attain genetic diagnoses. In addition, new genes causative of epilepsy are being uncovered through both traditional and novel methods, including large international data-sharing collaborations and massive sequencing efforts as well as computational methods and analyses driven by the Human Phenotype Ontology (HPO).

SUMMARY

New approaches to gene discovery and characterization are advancing rapidly our understanding of the genetic and phenotypic architecture of the epilepsies. This review highlights relevant and groundbreaking studies published recently that have pushed forward the field of epilepsy genetics.

High expression of AP2M1 correlates with worse prognosis by regulating immune microenvironment and drug resistance to R-CHOP in diffuse large B cell lymphoma

BACKGROUND

First-line treatment with R-CHOP has cured 50%-60% patients of diffuse large B cell lymphoma (DLBCL), and more than one-third patients will eventually progressed to relapsed/refractory disease with dismal outcomes. Adaptor Related Protein Complex 2 Subunit Mu 1 (AP2M1) is required for the activity of a vacuolar ATPase and may also play an important role in regulating the intracellular trafficking and function of CTLA-4 protein. Herein, using both public databases and our own tumor samples, we aimed to demonstrate the prognostic role of AP2M1 and the potential tumor-promoting mechanisms in DLBCL.

METHOD

Using public datasets of DLBCL from both GEO and TCGA databases, we analyzed the role of AP2M1 in mediating chemoresistance to R-CHOP and its correlation with various clinical parameters and prognosis. By using various R packages, we evaluated the role of AP2M1 on regulating tumor immune microenvironment. Moreover, tumor samples of DLBCL from Beijing TongRen Hospital were used to validate our findings by immunohistochemistry staining.

RESULT

Expression of AP2M1 was significantly increased in DLBCL, which was correlated with poor prognosis and a variety of clinical indicators. On the basis of enrichment analysis, it was found that AP2M1 may be related to intracellular receptor signaling pathway. Through immune analysis and drug prediction, we found that the expression of AP2M1 affected the immune environment and drug response of DLBCL, which further revealed the important role of AP2M1 in DLBCL. By analyzing 61 patients treated uniformly with R-CHOP regimen in our center, we validated the above findings that high expression of AP2M1 correlated with inferior survival outcomes and affected sensitivity to R-CHOP treatment.

CONCLUSION

Expression of AP2M1 may affect the prognosis of DLBCL patients probably by affecting the immune environment and the responses to many drugs in treating DLBCL, indicating AP2M1 as a potential therapy target in DLBCL.

Routine Diagnostics Confirm Novel Neurodevelopmental Disorders

Routine diagnostics is biased towards genes and variants with satisfactory evidence, but rare disorders with only little confirmation of their pathogenicity might be missed. Many of these genes can, however, be considered relevant, although they may have less evidence because they lack OMIM entries or comprise only a small number of publicly available variants from one or a few studies. Here, we present 89 individuals harbouring variants in 77 genes for which only a small amount of public evidence on their clinical significance is available but which we still found to be relevant enough to be reported in routine diagnostics. For 21 genes, we present case reports that confirm the lack or provisionality of OMIM associations (ATP6V0A1, CNTN2, GABRD, NCKAP1, RHEB, TCF7L2), broaden the phenotypic spectrum (CC2D1A, KCTD17, YAP1) or substantially strengthen the confirmation of genes with limited evidence in the medical literature (ADARB1, AP2M1, BCKDK, BCORL1, CARS2, FBXO38, GABRB1, KAT8, PRKD1, RAB11B, RUSC2, ZNF142). Routine diagnostics can provide valuable information on disease associations and support for genes without requiring tremendous research efforts. Thus, our results validate and delineate gene-disorder associations with the aim of motivating clinicians and scientists in diagnostic departments to provide additional evidence via publicly available databases or by publishing short case reports.

A recurrent de novo splice site variant involving DNM1 exon 10a causes developmental and epileptic encephalopathy through a dominant-negative mechanism

Heterozygous pathogenic variants in DNM1 cause developmental and epileptic encephalopathy (DEE) as a result of a dominant-negative mechanism impeding vesicular fission. Thus far, pathogenic variants in DNM1 have been studied with a canonical transcript that includes the alternatively spliced exon 10b. However, after performing RNA sequencing in 39 pediatric brain samples, we find the primary transcript expressed in the brain includes the downstream exon 10a instead. Using this information, we evaluated genotype-phenotype correlations of variants affecting exon 10a and identified a cohort of eleven previously unreported individuals. Eight individuals harbor a recurrent de novo splice site variant, c.1197-8G>A (GenBank: NM_001288739.1), which affects exon 10a and leads to DEE consistent with the classical DNM1 phenotype. We find this splice site variant leads to disease through an unexpected dominant-negative mechanism. Functional testing reveals an in-frame upstream splice acceptor causing insertion of two amino acids predicted to impair oligomerization-dependent activity. This is supported by neuropathological samples showing accumulation of enlarged synaptic vesicles adherent to the plasma membrane consistent with impaired vesicular fission. Two additional individuals with missense variants affecting exon 10a, p.Arg399Trp and p.Gly401Asp, had a similar DEE phenotype. In contrast, one individual with a missense variant affecting exon 10b, p.Pro405Leu, which is less expressed in the brain, had a correspondingly less severe presentation. Thus, we implicate variants affecting exon 10a as causing the severe DEE typically associated with DNM1-related disorders. We highlight the importance of considering relevant isoforms for disease-causing variants as well as the possibility of splice site variants acting through a dominant-negative mechanism.

Computational analysis of neurodevelopmental phenotypes: Harmonization empowers clinical discovery

Making a specific diagnosis in neurodevelopmental disorders is traditionally based on recognizing clinical features of a distinct syndrome, which guides testing of its possible genetic etiologies. Scalable frameworks for genomic diagnostics, however, have struggled to integrate meaningful measurements of clinical phenotypic features. While standardization has enabled generation and interpretation of genomic data for clinical diagnostics at unprecedented scale, making the equivalent breakthrough for clinical data has proven challenging. However, increasingly clinical features are being recorded using controlled dictionaries with machine readable formats such as the Human Phenotype Ontology (HPO), which greatly facilitates their use in the diagnostic space. Improving the tractability of large-scale clinical information will present new opportunities to inform genomic research and diagnostics from a clinical perspective. Here, we describe novel approaches for computational phenotyping to harmonize clinical features, improve data translation through revising domain-specific dictionaries, quantify phenotypic features, and determine clinical relatedness. We demonstrate how these concepts can be applied to longitudinal phenotypic information, which represents a critical element of developmental disorders and pediatric conditions. Finally, we expand our discussion to clinical data derived from electronic medical records, a largely untapped resource of deep clinical information with distinct strengths and weaknesses.

Genetic Testing for Malformations of Cortical Development

Background and Objectives

Malformations of cortical development (MCD), though individually rare, constitute a significant burden of disease. The diagnostic yield of next-generation sequencing (NGS) in these patients varies across studies and methods, and novel genes and variants continue to emerge.

Methods

Patients (n = 123) with a definite radiologic or histopathologic diagnosis of MCD, with or without epilepsy were included in this study. They underwent NGS-based targeted gene panel (TGP) testing, whole-exome sequencing (WES), or WES-based virtual panel testing. Selected patients who underwent epilepsy surgery (n = 69) also had somatic gene testing of brain tissue–derived DNA. We analyzed predictors of positive germline genetic finding and diagnostic yield of respective methods.

Results

Pathogenic or likely pathogenic germline genetic variants were detected in 21% of patients (26/123). In the surgical subgroup (69/123), we performed somatic sequencing in 40% of cases (28/69) and detected causal variants in 18% (5/28). Diagnostic yield did not differ between TGP, WES-based virtual gene panel, and open WES (p = 0.69). Diagnosis of focal cortical dysplasia type 2A, epilepsy, and intellectual disability were associated with positive results of germline testing. We report previously unpublished variants in 16/26 patients and 4 cases of MCD with likely pathogenic variants in non-MCD genes.

Discussion

In this study, we are reporting genetic findings of a large cohort of MCD patients with epilepsy or potentially epileptogenic MCD. We determine predictors of successful ascertainment of a genetic diagnosis in real-life setting and report novel, likely pathogenic variants in MCD and non-MCD genes alike.

PheNominal: an EHR-integrated web application for structured deep phenotyping at the point of care

BACKGROUND

Clinical phenotype information greatly facilitates genetic diagnostic interpretations pipelines in disease. While post-hoc extraction using natural language processing on unstructured clinical notes continues to improve, there is a need to improve point-of-care collection of patient phenotypes. Therefore, we developed "PheNominal", a point-of-care web application, embedded within Epic electronic health record (EHR) workflows, to permit capture of standardized phenotype data.

METHODS

Using bi-directional web services available within commercial EHRs, we developed a lightweight web application that allows users to rapidly browse and identify relevant terms from the Human Phenotype Ontology (HPO). Selected terms are saved discretely within the patient's EHR, permitting reuse both in clinical notes as well as in downstream diagnostic and research pipelines.

RESULTS

In the 16 months since implementation, PheNominal was used to capture discrete phenotype data for over 1500 individuals and 11,000 HPO terms during clinic and inpatient encounters for a genetic diagnostic consultation service within a quaternary-care pediatric academic medical center. An average of 7 HPO terms were captured per patient. Compared to a manual workflow, the average time to enter terms for a patient was reduced from 15 to 5 min per patient, and there were fewer annotation errors.

CONCLUSIONS

Modern EHRs support integration of external applications using application programming interfaces. We describe a practical application of these interfaces to facilitate deep phenotype capture in a discrete, structured format within a busy clinical workflow. Future versions will include a vendor-agnostic implementation using FHIR. We describe pilot efforts to integrate structured phenotyping through controlled dictionaries into diagnostic and research pipelines, reducing manual effort for phenotype documentation and reducing errors in data entry.

Computational insights into missense mutations in HTT gene causing Huntington's disease and its interactome networks

BACKGROUND

Huntington's disease is a rare neurodegenerative illness of the central nervous system that is inherited in an autosomal dominant pattern. Mutant huntingtin protein is produced as a result of enlargement of CAG repeat in the N-terminal of the polyglutamine tract.

AIM OF THE STUDY

Herein, we aim to investigate the mutations and their effects on the HTT gene and its genetic variants. Additionally, the protein-protein interaction of HTT with other proteins and receptor-ligand interaction with the three-dimensional structure of huntingtin protein were identified.

METHODS

A comprehensive analysis of the HTT interactome and protein-ligand interaction has been carried out to provide a global picture of structure-function analysis of huntingtin protein. Mutations were analyzed and mutation verification tools were used to check the effect of mutation on protein function.

RESULTS

The results showed, mutations in a single gene are not only responsible for causing a particular disease but may also cause other hereditary disorders as well. Moreover, the modification at the nucleotide level also cause the change in the specific amino acid which may disrupt the function of HTT and its interacting proteins contributing in disease pathogenesis. Furthermore, the interaction between MECP2 and BDNF lowers the rate of transcriptional activity. Molecular docking further confirmed the strong interaction between MECP2 and BDNF with highest affinity. Amino acid residues of the HTT protein, involved in the interaction with tetrabenazine were N912, Y890, G2385, and V2320. These findings proved, tetrabenazine as one of the potential therapeutic agent for treatment of Huntington's disease.

CONCLUSION

These results give further insights into the genetics of Huntington's disease for a better understanding of disease models which will be beneficial for the future therapeutic studies.

Creation and evaluation of full-text literature-derived, feature-weighted disease models of genetically determined developmental disorders

There are >2500 different genetically determined developmental disorders (DD), which, as a group, show very high levels of both locus and allelic heterogeneity. This has led to the wide-spread use of evidence-based filtering of genome-wide sequence data as a diagnostic tool in DD. Determining whether the association of a filtered variant at a specific locus is a plausible explanation of the phenotype in the proband is crucial and commonly requires extensive manual literature review by both clinical scientists and clinicians. Access to a database of weighted clinical features extracted from rigorously curated literature would increase the efficiency of this process and facilitate the development of robust phenotypic similarity metrics. However, given the large and rapidly increasing volume of published information, conventional biocuration approaches are becoming impractical. Here, we present a scalable, automated method for the extraction of categorical phenotypic descriptors from the full-text literature. Papers identified through literature review were downloaded and parsed using the Cadmus custom retrieval package. Human Phenotype Ontology terms were extracted using MetaMap, with 76-84% precision and 65-73% recall. Mean terms per paper increased from 9 in title + abstract, to 68 using full text. We demonstrate that these literature-derived disease models plausibly reflect true disease expressivity more accurately than widely used manually curated models, through comparison with prospectively gathered data from the Deciphering Developmental Disorders study. The area under the curve for receiver operating characteristic (ROC) curves increased by 5-10% through the use of literature-derived models. This work shows that scalable automated literature curation increases performance and adds weight to the need for this strategy to be integrated into informatic variant analysis pipelines. Database URL: https://doi.org/10.1093/database/baac038.

Assessing the landscape of STXBP1-related disorders in 534 individuals

Disease-causing variants in STXBP1 are among the most common genetic causes of neurodevelopmental disorders. However, the phenotypic spectrum in STXBP1-related disorders is wide and clear correlations between variant type and clinical features have not been observed so far. Here, we harmonized clinical data across 534 individuals with STXBP1-related disorders and analysed 19 973 derived phenotypic terms, including phenotypes of 253 individuals previously unreported in the scientific literature. The overall phenotypic landscape in STXBP1-related disorders is characterized by neurodevelopmental abnormalities in 95% and seizures in 89% of individuals, including focal-onset seizures as the most common seizure type (47%). More than 88% of individuals with STXBP1-related disorders have seizure onset in the first year of life, including neonatal seizure onset in 47%. Individuals with protein-truncating variants and deletions in STXBP1 (n = 261) were almost twice as likely to present with West syndrome and were more phenotypically similar than expected by chance. Five genetic hotspots with recurrent variants were identified in more than 10 individuals, including p.Arg406Cys/His (n = 40), p.Arg292Cys/His/Leu/Pro (n = 30), p.Arg551Cys/Gly/His/Leu (n = 24), p.Pro139Leu (n = 12), and p.Arg190Trp (n = 11). None of the recurrent variants were significantly associated with distinct electroclinical syndromes, single phenotypic features, or showed overall clinical similarity, indicating that the baseline variability in STXBP1-related disorders is too high for discrete phenotypic subgroups to emerge. We then reconstructed the seizure history in 62 individuals with STXBP1-related disorders in detail, retrospectively assigning seizure type and seizure frequency monthly across 4433 time intervals, and retrieved 251 anti-seizure medication prescriptions from the electronic medical records. We demonstrate a dynamic pattern of seizure control and complex interplay with response to specific medications particularly in the first year of life when seizures in STXBP1-related disorders are the most prominent. Adrenocorticotropic hormone and phenobarbital were more likely to initially reduce seizure frequency in infantile spasms and focal seizures compared to other treatment options, while the ketogenic diet was most effective in maintaining seizure freedom. In summary, we demonstrate how the multidimensional spectrum of phenotypic features in STXBP1-related disorders can be assessed using a computational phenotype framework to facilitate the development of future precision-medicine approaches.

Cellular abundance of sodium phosphate cotransporter SLC20A1/PiT1 and phosphate uptake are controlled post-transcriptionally by ESCRT

Inorganic phosphate is essential for human life. The widely expressed mammalian sodium/phosphate cotransporter SLC20A1/PiT1 mediates phosphate uptake into most cell types; however, while SLC20A1 is required for development, and elevated SLC20A1 expression is associated with vascular calcification and aggressive tumor growth, the mechanisms regulating SLC20A1 protein abundance are unknown. Here, we found that SLC20A1 protein expression is low in phosphate-replete cultured cells but is strikingly induced following phosphate starvation, whereas mRNA expression is high in phosphate-replete cells and only mildly increased by phosphate starvation. To identify regulators of SLC20A1 protein levels, we performed a genome-wide CRISPR-based loss-of-function genetic screen in phosphate-replete cells using SLC20A1 protein induction as readout. Our screen revealed that endosomal sorting complexes required for transport (ESCRT) machinery was essential for proper SLC20A1 protein downregulation in phosphate-replete cells. We show that SLC20A1 colocalizes with ESCRT and that ESCRT deficiency increases SLC20A1 protein and phosphate uptake into cells. We also found numerous additional candidate regulators of mammalian phosphate homeostasis, including genes modifying protein ubiquitination and the Krebs cycle and oxidative phosphorylation pathways. Many of these targets have not been previously implicated in this process. We present here a model in which SLC20A1 protein abundance and phosphate uptake are tonically negatively regulated post-transcriptionally in phosphate-replete cells through direct ESCRT-mediated SLC20A1 degradation. Moreover, our screening results provide a comprehensive resource for future studies to elucidate the mechanisms governing cellular phosphate homeostasis. We conclude that genome-wide CRISPR-based genetic screening is a powerful tool to discover proteins and pathways relevant to physiological processes.

Cardiofacioneurodevelopmental syndrome: Report of a novel patient and expansion of the phenotype

The cardiofacioneurodevelopmental syndrome (CFNDS) is characterized by craniofacial anomalies including bilateral cleft lip and palate, cardiac, skeletal, and neurodevelopmental features and additional variable manifestations. Whole-exome sequencing revealed homozygous loss-of-function variants in CCDC32 (alternative name: C15orf57) in both previously described patients. ccdc32 deletion in zebrafish suggests a ciliary contribution to the pathomechanism. We report a 9-year-old female patient with CFNDS caused by a homozygous 32,583-bp deletion affecting CCDC32. Independent of the affected CCDC32 transcript variant this deletion likely leads to loss of the encoded protein. The patient had intellectual disability, marked hypertelorism, bilateral cleft lip and palate, and short stature. She had bilateral conductive hearing loss, small hands and feet, and finger abnormalities. Brain imaging disclosed hypoplastic corpus callosum. We describe a core phenotype comprising developmental delay and bilateral cleft lip and palate in the three individuals with CFNDS. Variable abnormalities of the face, brain, heart, fingers, and toes and postnatal growth retardation or microcephaly can be present. Possible involvement of the uncharacterized CCDC32 protein in the adapter protein 2 (AP2) complex regulating clathrin-mediated endocytosis has been reported. Cleft palate and cardiac defects observed in mice deficient of different AP2 subunits support a CCDC32 function in the AP2 complex.

Rab13 sustains breast cancer stem cells by supporting tumor-stroma crosstalk

Cancer stem cells (CSC) are supported by the tumor microenvironment, and non-CSCs can regain CSC phenotypes in certain niches, leading to limited clinical benefits of CSC-targeted therapy. A better understanding of the mechanisms governing the orchestration of the CSC niche could help improve the therapeutic targeting of CSCs. Here, we report that Rab13, a small GTPase, is highly expressed in breast CSCs (BCSCs). Rab13 depletion suppressed breast cancer cell stemness, tumorigenesis, and chemoresistance by reducing tumor-stroma crosstalk. Accordingly, Rab13 controlled the membrane translocation of CXCR1/2, allowing tumor cells to interact with tumor-associated macrophages and cancer-associated fibroblasts to establish a supportive BCSC niche. Targeting the Rab13-mediated BCSC niche with bardoxolone-methyl (CDDO-Me) prevented BCSC stemness in vitro and in vivo. These findings highlight the novel regulatory mechanism of Rab13 in BCSC, with important implications for the development of therapeutic strategies for disrupting the BCSC niche.

Imaging Endocytosis Dynamics in Health and Disease

Endocytosis is a critical process for cell growth and viability. It mediates nutrient uptake, guarantees plasma membrane homeostasis, and generates intracellular signaling cascades. Moreover, it plays an important role in dead cell clearance and defense against external microbes. Finally, endocytosis is an important cellular route for the delivery of nanomedicines for therapeutic treatments. Thus, it is not surprising that both environmental and genetic perturbation of endocytosis have been associated with several human conditions such as cancer, neurological disorders, and virus infections, among others. Over the last decades, a lot of research has been focused on developing advanced imaging methods to monitor endocytosis events with high resolution in living cells and tissues. These include fluorescence imaging, electron microscopy, and correlative and super-resolution microscopy. In this review, we outline the major endocytic pathways and briefly discuss how defects in the molecular machinery of these pathways lead to disease. We then discuss the current imaging methodologies used to study endocytosis in different contexts, highlighting strengths and weaknesses.

Computational normal mode analysis accurately replicates the activity and specificity profiles of CRISPR-Cas9 and high-fidelity variants

The CRISPR-Cas system has transformed the field of gene-editing and created opportunities for novel genome engineering therapeutics. The field has significantly progressed, and recently, CRISPR-Cas9 was utilized in clinical trials to target disease-causing mutations. Existing tools aim to predict the on-target efficacy and potential genome-wide off-targets by scoring a particular gRNA according to an array of gRNA design principles or machine learning algorithms based on empirical results of large numbers of gRNAs. However, such tools are unable to predict the editing outcome by variant Cas enzymes and can only assess potential off-targets related to reference genomes. Here, we employ normal mode analysis (NMA) to investigate the structure of the Cas9 protein complexed with its gRNA and target DNA and explore the function of the protein. Our results demonstrate the feasibility and validity of NMA to predict the activity and specificity of SpyCas9 in the presence of mismatches by comparison to empirical data. Furthermore, despite the absence of their exact structures, this method accurately predicts the enzymatic activity of known high-fidelity engineered Cas9 variants.

3q27.1 microdeletion causes prenatal and postnatal growth restriction and neurodevelopmental abnormalities

BACKGROUND

Overlapping microdeletions of chromosome 3q26-3q28 have been reported in eight individuals. The common phenotype observed in these individuals include intrauterine growth restriction, short stature, microcephaly, feeding difficulties, facial dysmorphisms, limb abnormalities and developmental delay. The most striking clinical features shared among all reported cases is prenatal and postnatal growth restriction and neurodevelopmental abnormalities.

CASE PRESENTATION

We identified two additional individuals with overlapping deletions and shared clinical features by high-resolution SNP oligonucleotide microarray, and refined the smallest region of overlap (SRO) to a 1.2 Mb genomic location in chromosome 3q27.1 by reviewing and comparing all published cases. We evaluated the SRO using ACMG/ClinGen current recommendations for classifying copy number variants (CNVs), and discussed the contribution of the genes deleted in the SRO to the abnormal phenotype observed in these individuals.

CONCLUSIONS

This study provides further evidence supporting the existence of a novel 3q27.1 microdeletion syndrome and suggests that haploinsufficiency of potential candidate genes, DVL3, AP2M1, and PARL in the SRO in 3q27.1 is responsible for the phenotype.

Cardiovascular Phenotypes Profiling for L-Transposition of the Great Arteries and Prognosis Analysis

Objectives

Congenitally corrected transposition of the great arteries (ccTGA) is a rare and complex congenital heart disease with the characteristics of double discordance. Enormous co-existed anomalies are the culprit of prognosis evaluation and clinical decision. We aim at delineating a novel ccTGA clustering modality under human phenotype ontology (HPO) instruction and elucidating the relationship between phenotypes and prognosis in patients with ccTGA.

Methods

A retrospective review of 270 patients diagnosed with ccTGA in Fuwai hospital from 2009 to 2020 and cross-sectional follow-up were performed. HPO-instructed clustering method was administered in ccTGA risk stratification. Kaplan-Meier survival, Landmark analysis, and cox regression analysis were used to investigate the difference of outcomes among clusters.

Results

The median follow-up time was 4.29 (2.07–7.37) years. A total of three distinct phenotypic clusters were obtained after HPO-instructed clustering with 21 in cluster 1, 136 in cluster 2, and 113 in cluster 3. Landmark analysis revealed significantly worse mid-term outcomes in all-cause mortality (p = 0.021) and composite endpoints (p = 0.004) of cluster 3 in comparison with cluster 1 and cluster 2. Multivariate analysis indicated that pulmonary arterial hypertension (PAH), atrioventricular septal defect (AVSD), and arrhythmia were risk factors for composite endpoints. Moreover, the surgical treatment was significantly different among the three groups (p < 0.001) and surgical strategies had different effects on the prognosis of the different phenotypic clusters.

Conclusions

Human phenotype ontology-instructed clustering can be a potentially powerful tool for phenotypic risk stratification in patients with complex congenital heart diseases, which may improve prognosis prediction and clinical decision.

Clathrin-independent endocytic retrieval of SV proteins mediated by the clathrin adaptor AP-2 at mammalian central synapses

Neurotransmission is based on the exocytic fusion of synaptic vesicles (SVs) followed by endocytic membrane retrieval and the reformation of SVs. Conflicting models have been proposed regarding the mechanisms of SV endocytosis, most notably clathrin/adaptor protein complex 2 (AP-2)-mediated endocytosis and clathrin-independent ultrafast endocytosis. Partitioning between these pathways has been suggested to be controlled by temperature and stimulus paradigm. We report on the comprehensive survey of six major SV proteins to show that SV endocytosis in mouse hippocampal neurons at physiological temperature occurs independent of clathrin while the endocytic retrieval of a subset of SV proteins including the vesicular transporters for glutamate and GABA depend on sorting by the clathrin adaptor AP-2. Our findings highlight a clathrin-independent role of the clathrin adaptor AP-2 in the endocytic retrieval of select SV cargos from the presynaptic cell surface and suggest a revised model for the endocytosis of SV membranes at mammalian central synapses.

In-silico phenotype prediction by normal mode variant analysis in TUBB4A-related disease

TUBB4A-associated disorder is a rare condition affecting the central nervous system. It displays a wide phenotypic spectrum, ranging from isolated late-onset torsion dystonia to a severe early-onset disease with developmental delay, neurological deficits, and atrophy of the basal ganglia and cerebellum, therefore complicating variant interpretation and phenotype prediction in patients carrying TUBB4A variants. We applied entropy-based normal mode analysis (NMA) to investigate genotype–phenotype correlations in TUBB4A-releated disease and to develop an in-silico approach to assist in variant interpretation and phenotype prediction in this disorder. Variants included in our analysis were those reported prior to the conclusion of data collection for this study in October 2019. All TUBB4A pathogenic missense variants reported in ClinVar and Pubmed, for which associated clinical information was available, and all benign/likely benign TUBB4A missense variants reported in ClinVar, were included in the analysis. Pathogenic variants were divided into five phenotypic subgroups. In-silico point mutagenesis in the wild-type modeled protein structure was performed for each variant. Wild-type and mutated structures were analyzed by coarse-grained NMA to quantify protein stability as entropy difference value (ΔG) for each variant. Pairwise ΔG differences between all variant pairs in each structural cluster were calculated and clustered into dendrograms. Our search yielded 41 TUBB4A pathogenic variants in 126 patients, divided into 11 partially overlapping structural clusters across the TUBB4A protein. ΔG-based cluster analysis of the NMA results revealed a continuum of genotype–phenotype correlation across each structural cluster, as well as in transition areas of partially overlapping structural clusters. Benign/likely benign variants were integrated into the genotype–phenotype continuum as expected and were clearly separated from pathogenic variants. We conclude that our results support the incorporation of the NMA-based approach used in this study in the interpretation of variant pathogenicity and phenotype prediction in TUBB4A-related disease. Moreover, our results suggest that NMA may be of value in variant interpretation in additional monogenic conditions.

Bi-allelic variants in OGDHL cause a neurodevelopmental spectrum disease featuring epilepsy, hearing loss, visual impairment, and ataxia

The 2-oxoglutarate dehydrogenase-like (OGDHL) protein is a rate-limiting enzyme in the Krebs cycle that plays a pivotal role in mitochondrial metabolism. OGDHL expression is restricted mainly to the brain in humans. Here, we report nine individuals from eight unrelated families carrying bi-allelic variants in OGDHL with a range of neurological and neurodevelopmental phenotypes including epilepsy, hearing loss, visual impairment, gait ataxia, microcephaly, and hypoplastic corpus callosum. The variants include three homozygous missense variants (p.Pro852Ala, p.Arg244Trp, and p.Arg299Gly), three compound heterozygous single-nucleotide variants (p.Arg673Gln/p.Val488Val, p.Phe734Ser/p.Ala327Val, and p.Trp220Cys/p.Asp491Val), one homozygous frameshift variant (p.Cys553Leufs^∗16), and one homozygous stop-gain variant (p.Arg440Ter). To support the pathogenicity of the variants, we developed a novel CRISPR-Cas9-mediated tissue-specific knockout with cDNA rescue system for dOgdh, the Drosophila ortholog of human OGDHL. Pan-neuronal knockout of dOgdh led to developmental lethality as well as defects in Krebs cycle metabolism, which was fully rescued by expression of wild-type dOgdh. Studies using the Drosophila system indicate that p.Arg673Gln, p.Phe734Ser, and p.Arg299Gly are severe loss-of-function alleles, leading to developmental lethality, whereas p.Pro852Ala, p.Ala327Val, p.Trp220Cys, p.Asp491Val, and p.Arg244Trp are hypomorphic alleles, causing behavioral defects. Transcript analysis from fibroblasts obtained from the individual carrying the synonymous variant (c.1464T>C [p.Val488Val]) in family 2 showed that the synonymous variant affects splicing of exon 11 in OGDHL. Human neuronal cells with OGDHL knockout exhibited defects in mitochondrial respiration, indicating the essential role of OGDHL in mitochondrial metabolism in humans. Together, our data establish that the bi-allelic variants in OGDHL are pathogenic, leading to a Mendelian neurodevelopmental disease in humans.

A modified density gradient proteomic-based method to analyze endolysosomal proteins in cardiac tissue

The importance of lysosomes in cardiac physiology and pathology is well established, and evidence for roles in calcium signaling is emerging. We describe a label-free proteomics method suitable for small cardiac tissue biopsies based on density-separated fractionation, which allows study of endolysosomal (EL) proteins. Density gradient fractions corresponding to tissue lysate; sarcoplasmic reticulum (SR), mitochondria (Mito) (1.3 g/mL); and EL with negligible contamination from SR or Mito (1.04 g/mL) were analyzed using Western blot, enzyme activity assay, and liquid chromatography with tandem mass spectrometry (LC-MS/MS) analysis (adapted discontinuous Percoll and sucrose differential density gradient). Kyoto Encyclopedia of Genes and Genomes, Reactome, Panther, and Gene Ontology pathway analysis showed good coverage of RAB proteins and lysosomal cathepsins (including cardiac-specific cathepsin D) in the purified EL fraction. Significant EL proteins recovered included catalytic activity proteins. We thus present a comprehensive protocol and data set of guinea pig atrial EL organelle proteomics using techniques also applicable for non-cardiac tissue.

Diagnostic yield of chromosomal microarray and trio whole exome sequencing in cryptogenic cerebral palsy

OBJECTIVE

To determine the yield of genetic diagnoses using chromosomal microarray (CMA) and trio whole exome sequencing (WES), separately and combined, among patients with cryptogenic cerebral palsy (CP).

METHODS

Trio WES of patients with prior CMA analysis for cryptogenic CP, defined as disabling, non-progressive motor symptoms beginning before the age of 3 years without known cause.

RESULTS

Given both CMA analysis and trio WES, clinically significant genetic findings were identified for 58% of patients (26 of 45). Diagnoses were eight large CNVs detected by CMA and 18 point mutations detected by trio WES. None had more than one severe mutation. Approximately half of events (14 of 26) were de novo. Yield was significantly higher in patients with CP with comorbidities (69%, 22 of 32) than in those with pure motor CP (31%, 4 of 13; p=0.02). Among patients with genetic diagnoses, CNVs were more frequent than point mutations among patients with congenital anomalies (OR 7.8, 95% CI 1.2 to 52.4) or major dysmorphic features (OR 10.5, 95% CI 1.4 to 73.7). Clinically significant mutations were identified in 18 different genes: 14 with known involvement in CP-related disorders and 4 responsible for other neurodevelopmental conditions. Three possible new candidate genes for CP were ARGEF10, RTF1 and TAOK3.

CONCLUSIONS

Cryptogenic CP is genetically highly heterogeneous. Genomic analysis has a high yield and is warranted in all these patients. Trio WES has higher yield than CMA, except in patients with congenital anomalies or major dysmorphic features, but these methods are complementary. Patients with negative results with one approach should also be tested by the other.

Phenotypic homogeneity in childhood epilepsies evolves in gene-specific patterns across 3251 patient-years of clinical data

While genetic studies of epilepsies can be performed in thousands of individuals, phenotyping remains a manual, non-scalable task. A particular challenge is capturing the evolution of complex phenotypes with age. Here, we present a novel approach, applying phenotypic similarity analysis to a total of 3251 patient-years of longitudinal electronic medical record data from a previously reported cohort of 658 individuals with genetic epilepsies. After mapping clinical data to the Human Phenotype Ontology, we determined the phenotypic similarity of individuals sharing each genetic etiology within each 3-month age interval from birth up to a maximum age of 25 years. 140 of 600 (23%) of all 27 genes and 3-month age intervals with sufficient data for calculation of phenotypic similarity were significantly higher than expect by chance. 11 of 27 genetic etiologies had significant overall phenotypic similarity trajectories. These do not simply reflect strong statistical associations with single phenotypic features but appear to emerge from complex clinical constellations of features that may not be strongly associated individually. As an attempt to reconstruct the cognitive framework of syndrome recognition in clinical practice, longitudinal phenotypic similarity analysis extends the traditional phenotyping approach by utilizing data from electronic medical records at a scale that is far beyond the capabilities of manual phenotyping. Delineation of how the phenotypic homogeneity of genetic epilepsies varies with age could improve the phenotypic classification of these disorders, the accuracy of prognostic counseling, and by providing historical control data, the design and interpretation of precision clinical trials in rare diseases.

Modeling seizures in the Human Phenotype Ontology according to contemporary ILAE concepts makes big phenotypic data tractable

Objective:

The clinical features of epilepsy determine how it is defined, which in turn guides management. Therefore, consideration of the fundamental clinical entities that comprise an epilepsy is essential in the study of causes, trajectories, and treatment responses. The Human Phenotype Ontology (HPO) is used widely in clinical and research genetics for concise communication and modeling of clinical features, allowing extracted data to be harmonized using logical inference. We sought to redesign the HPO seizure subontology to improve its consistency with current epileptological concepts, supporting the use of large clinical data sets in high-throughput clinical and research genomics.

Methods:

We created a new HPO seizure subontology based on the 2017 International League Against Epilepsy (ILAE) Operational Classification of Seizure Types, and integrated concepts of status epilepticus, febrile, reflex, and neonatal seizures at different levels of detail. We compared the HPO seizure subontology prior to, and following, our revision, according to the information that could be inferred about the seizures of 791 individuals from three independent cohorts: 2 previously published and 150 newly recruited individuals. Each cohort’s data were provided in a different format and harmonized using the two versions of the HPO.

Results:

The new seizure subontology increased the number of descriptive concepts for seizures 5-fold. The number of seizure descriptors that could be annotated to the cohort increased by 40% and the total amount of information about individuals’ seizures increased by 38%. The most important qualitative difference was the relationship of focal to bilateral tonic-clonic seizure to generalized-onset and focal-onset seizures.

Mild Phenotype of Wolfram Syndrome Associated With a Common Pathogenic Variant Is Predicted by a Structural Model of Wolframin

Objective

To describe the WFS1 c.1672C>T; p.R558C missense variant, found in 1.34% of Ashkenazi Jews, that has a relatively mild phenotype and to use computational normal mode analysis (NMA) to explain the genotype-phenotype relationship.

Methods

The clinical, laboratory, and genetic features of 8 homozygotes were collected. A model of the wolframin protein was constructed, and NMA was used to simulate the effect of the variant on protein thermodynamics.

Results

Mean age at Wolfram syndrome (WS) diagnosis among homozygotes was 30 years; diabetes (7/8) was diagnosed at mean age 19 years (15–21 years), and bilateral optic atrophy (with MRI evidence of optic/chiasm atrophy) (6/8) at mean age 29 years (15–48 years). The oldest patient (62 years) also had gait difficulties, memory problems, parietal and cerebellar atrophy, and white matter hyperintense lesions. All retained functional vision with independent ambulation and self-care; none had diabetes insipidus or hearing loss. The p.R558C variant caused less impairment of protein entropy than WFS1 variants associated with a more severe phenotype.

Conclusions

The p.R558C variant causes a milder, late-onset phenotype of WS. We report a structural model of wolframin protein based on empirical functional studies and use NMA modeling to show a genotype-phenotype correlation across all homozygotes. Clinicians should be alert to this condition in patients with juvenile diabetes and patients of any age with a combination of diabetes and optic atrophy. Computational NMA has potential benefit for prediction of the genotype-phenotype relationship.

Ap2s1 mutation causes hypercalcaemia in mice and impairs interaction between calcium-sensing receptor and adaptor protein-2

Adaptor protein 2 (AP2), a heterotetrameric complex comprising AP2α, AP2β2, AP2μ2 and AP2σ2 subunits, is ubiquitously expressed and involved in endocytosis and trafficking of membrane proteins, such as the calcium-sensing receptor (CaSR), a G-protein coupled receptor that signals via Gα₁₁. Mutations of CaSR, Gα₁₁ and AP2σ2, encoded by AP2S1, cause familial hypocalciuric hypercalcaemia types 1–3 (FHH1–3), respectively. FHH3 patients have heterozygous AP2S1 missense Arg15 mutations (p.Arg15Cys, p.Arg15His or p.Arg15Leu) with hypercalcaemia, which may be marked and symptomatic, and occasional hypophosphataemia and osteomalacia. To further characterize the phenotypic spectrum and calcitropic pathophysiology of FHH3, we used CRISPR/Cas9 genome editing to generate mice harboring the AP2S1 p.Arg15Leu mutation, which causes the most severe FHH3 phenotype. Heterozygous (Ap2s1^+/L15) mice were viable, and had marked hypercalcaemia, hypermagnesaemia, hypophosphataemia, and increases in alkaline phosphatase activity and fibroblast growth factor-23. Plasma 1,25-dihydroxyvitamin D was normal, and no alterations in bone mineral density or bone turnover were noted. Homozygous (Ap2s1^L15/L15) mice invariably died perinatally. Co-immunoprecipitation studies showed that the AP2S1 p.Arg15Leu mutation impaired protein–protein interactions between AP2σ2 and the other AP2 subunits, and also with the CaSR. Cinacalcet, a CaSR positive allosteric modulator, decreased plasma calcium and parathyroid hormone concentrations in Ap2s1^+/L15 mice, but had no effect on the diminished AP2σ2-CaSR interaction in vitro. Thus, our studies have established a mouse model that is representative for FHH3 in humans, and demonstrated that the AP2S1 p.Arg15Leu mutation causes a predominantly calcitropic phenotype, which can be ameliorated by treatment with cinacalcet.