A novel autism-associated UBLCP1 mutation impacts proteasome regulation/activity – Translational Psychiatry

Abstract The landscape of autism spectrum disorder (ASD) in Lebanon is unique because of high rates of consanguinity, shared ancestry, and increased remote consanguinity. ASD prevalence in Lebanon is 1 in 68 with a male-to-female ratio of 2:1. This study aims to investigate the impact of an inherited deletion in UBLCP1 (Ubiquitin-Like Domain-Containing CTD Phosphatase 1) on the ubiquitin-proteasome system (UPS) and proteolysis. Whole exome sequencing in a Lebanese family with ASD without pathogenic copy number variations (CNVs) uncovered a deletion in UBLCP1 . Functional evaluation of the identified variant is described in fibroblasts from the affected. The deletion in UBLCP1 exon 10 (g.158,710,261CAAAG > C) generates a premature stop codon interrupting the phosphatase domain and is predicted as pathogenic. It is absent from databases of normal variation worldwide and in Lebanon. Wild-type UBLCP1 is widely expressed in mouse brains. The mutation results in decreased UBLCP1 protein expression in patient-derived fibroblasts from the autistic patient compared to controls. The truncated UBLCP1 protein results in increased proteasome activity decreased ubiquitinated protein levels, and downregulation in expression of other proteasome subunits in samples from the affected compared to controls. Inhibition of the proteasome by using MG132 in proband cells reverses alterations in gene expression due to the restoration of protein levels of the common transcription factor, NRF1. Finally, treatment with gentamicin, which promotes premature termination codon read-through, restores UBLCP1 expression and function. Discovery of an ASD-linked mutation in UBLCP1 leading to overactivation of cell proteolysis is reported. This, in turn, leads to dysregulation of proteasome subunit transcript levels as a compensatory response. Introduction Autism spectrum disorder (ASD) is a neurodevelopmental disorder diagnosed with: (1) deficits in social communication and interaction and (2) restricted, repetitive patterns of behavior [1]. ASD has strong genetic components [2, 3]. The concordance rate is 60–70% in monozygotic twins and 5–30% in siblings [4], with recurrence rates of 18% in infant siblings and 33% in multiplex families [4, 5]. ASD is diagnosed at a rate of 1 in 44 children in the USA, with males three to fourfold [6] more likely to be diagnosed than females. In Lebanon, an incidence of 1 in 68 is reported, with a male-to-female ratio of 2:1 [7, 8]. Next-generation sequencing has led to the identification of hundreds of genes linked to ASD [9,10,11]. The Lebanese population is genetically distinguished by a distinctive blend of DNA recombination events occurring over millennia [12, 13]. Significant shared ancestry and high rates of consanguinity make it ideally suited for genetic analysis of inherited causes of autism [14]. In this study, a novel homozygous pathogenic variant in UBLCP1 in a Lebanese consanguineous family with ASD is described. The ubiquitin-like domain-containing C-terminal domain phosphatase 1 (UBLCP1) is a phosphatase that regulates proteasome activity. The ubiquitin-proteasome system (UPS) is responsible for the degradation of thousands of short-lived and damaged proteins. Increasing focus on the UPS role in neurodevelopmental disorders, specifically ASD, is emerging. Mutations in several genes encoding members of the UPS are associated with ASD. These genes encode ubiquitin E3 ligases ( UBE3A/B/C , CUL3 , PRKN , FBXO40 , and HUWE1 [15,16,17,18]) and proteasomal subunits ( PSMD10 [19] and PSMD12 [20]). The 26 S proteasome consists of a 20 S core particle (CP) and a 19 S regulatory particle (RP) [21]. CP is formed of subunits with peptidase activities. AAA + -type ATPases (Rpt1-6) together with non-ATPase subunits Rpn1, 2, and 13 form the “base” of RP, which directly binds to the CP. The “lid” of RP is composed of Rpn proteins (Rpn3, 5-12, and 15). Polyubiquitinated protein substrates are captured by RP and ubiquitin chains are cleaved off by deubiquitinases [22]. Unfolded substrates are then translocated into the CP catalytic chamber for degradation. To maintain cell viability, UPS activity must be regulated via translational and post-translational mechanisms. Phosphatases reversibly control phosphorylation of proteasome subunits by kinases [23, 24]. UBLCP1 was the first proteasome phosphatase identified [25]. It harbors the N-terminal ubiquitin-like (UBL) domain [26] and the CTD phosphatase domain at the C-terminus. UBLCP1 interacts directly with the 19 S subunit Rpn1 through the UBL domain and dephosphorylates Rpn1-Ser361, a phosphosite required for proper assembly and activity of the 26 S proteasome [25, 27, 28]. UBLCP1 also dephosphorylates Rpt1 at multiple phosphorylation sites impairing ATPase activity [29]. Overall, these multiple dephosphorylations result in inhibition of proteasome activity by interfering with proteasome assembly [25, 29, 30]. A novel UBLCP1 mutation is identified in a Lebanese family with ASD, resulting in a premature stop codon that disrupts UBLCP1 protein phosphatase domain, disturbing proteasome regulation and activity. Material and methods Subjects and specimens Blood was collected from a 20-year-old ASD patient, who fulfilled DSM-5 criteria for autism, his parents, and two healthy siblings according to a protocol and consent forms reviewed and approved by the Institutional Review Board (IRB) of the American University of Beirut (IRB protocol number: BIOCH.RB.06, protocol name: Autism susceptibility genes in the Lebanese population). All experimental protocols were carried out in accordance with institutional guidelines. Informed consent was obtained by CITI-certified investigators from all family members, including both parents, siblings, and by parents on behalf of the affected individual. Blood samples from a cohort of 100 unaffected controls were collected. DNA samples were extracted from peripheral blood using QIAamp® blood midi kit (Qiagen, Inc., Valencia, CA). Whole exome sequencing, data processing, and variant identification Whole exome sequencing (WES) was performed on all family members, including the individual with ASD, both parents and healthy siblings. DNA was sequenced on the Illumina HiSeq 2000 platform at the Broad Institute (Massachusetts Institute of Technology) with a mean target coverage of 90% at ≥20× and a mean read depth of 102×. WES data was processed following best practices recommended by the Broad Institute [31]. Reads were aligned to the human reference genome version GRCh37/hg19 using the Burrows-Wheeler Aligner (BWA, version 0.7.10). Duplicate reads were removed using Picard (version 1.117). Local realignment and quality recalibration were performed using the Genome Analysis Toolkit (GATK; version 3.3). Variants (single nucleotide variants (SNVs) and insertions or deletions (indels) were detected using GATK with HaplotypeCaller. Quality control checks for (i) duplicate samples, (ii) samples per platform, (iii) genome call rate, (iv) missingness rate, (v) singleton rate, (vi) heterozygosity rate, (vii) homozygosity rate, (viii) Ti/Tv ratio, (ix) inbreeding coefficient, and (x) sex inference were performed. Variant call format (VCF) files for SNVs and indels were annotated with ANNOVAR using allele frequencies from the 1000 Genomes project (2015; 1000 G), the Genome Aggregation Database (gnomAD), and the Greater Middle East Variome Project (GME). Annotated VCF files were uploaded into an SQL database for working storage and analysis. Genome data was stored and analyses were performed on the Texas Advanced Computing Center (TACC) high-performance computing server, a resource of the University of Texas (Austin, TX). Sanger sequencing was performed on DNA samples from the affected individual and family members to validate the WES finding of the UBLCP1 variant (g. 158,710,261CAAAG > C) by amplification using primers UBLCP1_Foward 5′-GCTGGTGTTGTCAGGTTCAG-3′ and UBLCP1_Reverse 5′-GGGGTCACTGATGAAAAATGA-3′). Dermal fibroblast derivation and expansion Skin biopsy from the proband and from control age-matched Lebanese males were transferred into sterile glass Petri dishes containing 4 ml freshly prepared fibroblast media (Dulbecco’s modified Eagle’s medium high glucose with GlutaMAX (Life Technologies, Carlsbad, CA, USA), 10% certified fetal bovine serum (Life Technologies, Carlsbad, CA, USA), 1% antibiotics (Sigma, St. Louis, MO, USA)). Four scratches were made in each well of a 12-well plate using a 19-gauge needle to produce a checkerboard pattern. Biopsy was minced into small 2-mm pieces by sharply cutting edges using a new disposable stainless-steel blade scalpel. The biopsy pieces were transferred onto the 12-well culture plate (3 pieces/well) to the intersection of two scratches so that cells could adhere and migrate from the sample onto the culture plate. Tissue pieces were allowed to dry for 3 minutes (min) under the hood, and a few drops of media were added. After incubating for 30 min at 37 °C, an appropriate volume of media was added to submerge the tissue pieces. Explants were subjected weekly to a complete medium change. As fibroblasts grew out from the explants and expanded to 50% confluency, tissue pieces were removed. Fibroblast cultures were established from skin biopsies after 3–4 weeks. Fibroblasts were detached from the flask using prewarmed 37 °C Trypsin (Lonza, Houston, TX, USA) and split into larger culture flasks for expansion. Cells were frozen at a density of 1.0 × 10 6 cells/milliliter/cryovial. All cell lines used/frozen tested negative for mycoplasma contamination. Fibroblast read-through experiments and proteasome inhibition Early passage primary skin fibroblasts were plated in 25 cm 2 flasks and grown to 80% confluency. For read-through experiments, cells were treated with gentamicin (Sigma-Aldrich, St. Louis, MO, USA) at 100 μM to promote premature termination codon (PTC) read-through, then harvested after 24 hours (hr) of treatment. For proteasome inhibition, cells were treated with MG132 (Calbiochem, San Diego, CA, USA) at 10 μM, or just with DMSO vehicle, then harvested after 10 hr of treatment. Control fibroblasts without treatment, control fibroblasts treated with gentamicin or MG132, patient-derived fibroblasts without treatment, and patient-derived fibroblasts treated with gentamicin or MG132 were assayed using biological triplicates, and technical triplicates for each biological replicate. Proteasome activity Protein samples (2 μg) were incubated with 100 μM succinyl-Leu-Leu-Val-Tyr-7-amido-4-methylcoumarin (Suc-LLVY-AMC) in a proteasome activity buffer (50 mM Tris–HCl (pH 7.5), 40 mM KCl, 5 mM MgCl 2 , 1 mM DTT, 2 mM ATP) for a final volume of 100 μl. The release of fluorescent free-AMC was quantified at an excitation of 380 nm and emission of 460 nm. After 1 hr incubation in the dark at 37 °C, fluorescence was measured using a spectrophotometer. All experiments were done using biological triplicates, as well as technical triplicates for each biological replicate. Western blot and immunohistochemistry For western blot analysis, experiments were performed using biological triplicates, as well as technical triplicates for each biological replicate. Cells were harvested and lysed in RIPA lysis buffer with a protease inhibitor cocktail (Roche, Basel, Switzerland). Protein concentrations of whole cell lysates were determined by Bradford protein assay using Coomassie Plus Protein Assay Reagent (Fisher Scientific, Ottawa, ON, Canada). Proteins were separated after denaturation for 5 min at 95 °C on a 10% acrylamide gel in a Tris-glycine running buffer system and transferred onto a PVDF membrane with a Tris-glycine buffer system using a wet transfer unit (BioRad, Hercules, CA, USA) at 100 volts for 40 min. Membranes were blocked using 5% bovine serum albumin and incubated with antibody in TBST overnight at 4 °C with antibodies to Ubiquitin (1:1000; Cell Signaling, Massachusetts, USA) and GAPDH (1:1000; Cell Signaling, Massachusetts, USA) or Vinculin (1:1000; Cell Signaling, Massachusetts, USA). Membranes were then incubated in HRP-conjugated secondary antibodies. Signals were detected by chemiluminescence using the SuperSignal West Femto Maximum Sensitivity Substrate (ThermoFisher Scientific, MA, USA) and visualized using ChemiDoc MP Imaging System (BioRad, Hercules, CA, USA). Quantification was performed using ImageJ. Bands corresponding to the protein of interest, or their respective loading control (GAPDH or Vinculin), were selected, as well as an additional “blank” background band which was used to subtract background signal/noise from the measurements. The resulting values were then normalized to their respective reference band (GAPDH or Vinculin). For immunofluorescence experiments, cells previously plated on coverslips in 24-well plates were fixed with 4% paraformaldehyde/PBS for 20 min, washed three times with PBS, and incubated for 45 min in blocking buffer (PBS, 10% fetal calf serum, and 0.1% Triton X-100). Cells were incubated overnight with primary antibody (anti-UBLCP1, 1:200, Sigma, St. Louis, MO, USA), then secondary antibody coupled with fluorescence. The nucleus was counterstained with Hoechst dye. Fluorescence was visualized on an inverted confocal microscope (Leica, Germany). Parameters (such as brightness, contrast, and fluorescence threshold) remained unchanged and were standardized for all quantified images. When comparing mutated versus control fibroblasts, parameters were set for low-brightness signals within the dynamic range, and maintained equal for all images. Quantification of fluorescence intensity was performed using the ImageJ software. For brain tissue collection and processing, C57BL-6 adult mice ( n  = 3, 2 males and 1 female) were perfused using 4% paraformaldehyde (4% PFA), and brains harvested. Post-natal day 4 C57BL-6 mice ( n  = 4, 2 males, and 2 females) were sacrificed by decapitation, dissected, and fixed in 4% PFA for 4 hours. Fixed brains were subsequently immersed in 20% sucrose for 24 hours before freezing, then embedded in a cryo-embedding matrix and frozen. Frozen brains were sectioned using a cryotome at a thickness of 20 microns per section. Antibody concentrations used for brain tissue immunofluorescence are as follows: anti-Map2, 1:500 (Invitrogen, MA, USA); anti-UBLCP1, 1:200 (Sigma, St. Louis, MO, USA); anti-CD11b, 1:50 (Cell signaling, Massachusetts, USA); anti-TH, 1:500 (Chemicon, Massachusetts, USA); anti-Olig2, 1:100 (Millipore, Massachusetts, USA). Fluorescence was visualized on an inverted confocal microscope (Leica, Germany). No randomization or blinding was used for these mice. All animal experiments were performed with approval from the Animal Care Committee of the American University of Beirut, Lebanon. RNA extraction and qRT-PCR Total RNA was extracted from fibroblasts and reverse transcribed using RevertAid Reverse Transcriptase (ThermoFisher Scientific, MA, USA) with 100-1000 ng of input RNA and random primers (ThermoFisher Scientific, MA, USA). Quantitative real-time PCR reactions (qRT-PCR) were performed in 96- or 384-well plates using specific primers (TIB MOLBIOL, Germany), and the iQTM SYBR® Green Supermix (BioRad, CA, USA) was used as a fluorescent detection dye in a final volume of 10 μl (CFX96TM Real-Time PCR (BioRad, CA, USA)). Each reaction was performed using biological triplicates, as well as technical triplicates for each biological replicate. All results are normalized to β-actin or GAPDH mRNA levels and calculated using the 2 ∆Ct method. Melt curve analysis was applied to characterize generated amplicons and to control contamination by unspecific by-products. Statistical analysis Statistical analysis was conducted using GraphPad Prism 6 software (GraphPad software, CA, USA). Continuous data was expressed as mean ± SEM. Statistical analysis was performed using two-tailed Student’s t test and one-way analysis of variance (ANOVA) followed by Tukey post hoc test for multiple group comparisons. Values of P  < 0.05 were considered statistically significant. Results Clinical presentation, whole exome sequencing, and variant identification F2-05, a boy diagnosed with ASD, is the third child of a Lebanese family from a first-cousin marriage (Fig. 1a). Previous microarray analysis in the affected individual showed the absence of known or predicted ASD-associated CNVs. To identify putative causative variants, whole exome sequencing (WES) was performed on all family members, including the proband, both parents, and both unaffected siblings. Assuming a homozygous recessive model of inheritance, led to the identification of variants in 5 candidate genes (Fig. 1b, Table 1). Of these variants, three were intronic identified in candidate genes GABRB2 , HDAC9, and TAS1R2 . These variants were all predicted to be probably benign by MutationTaster. Additional transcriptional analysis of GABRB2 , which encodes the γ-aminobutyric acid (GABA) A receptor beta 2 subunit, ruled out splice site involvement and showed no functional impact for the identified variant (data not shown). The remaining 2 identified variants were non-synonymous exonic: a missense variant in WNK1 (g.1,005,622 C > T, p.P2488L) and a frameshift variant in UBLCP1 (g.158,710,261CAAAG > C, p.K282Nfs*3). The WNK1 variant was predicted to be disease causing by MutationTaster, PolyPhen2, and SIFT, but was predicted as a polymorphism by UMD predictor. In addition, the mutated residue was not conserved across species (Table 1). This homozygous frameshift deletion in UBLCP1 , predicted to be deleterious by MutationTaster (Table 1), leads to a downstream premature termination codon (PTC) resulting in truncation of the phosphatase domain (Fig. 2a–c). The affected residue was partially conserved across species (Fig. 2d). To confirm WES result, we performed Sanger sequencing on samples from all family members. The homozygous g.158,710,261CAAAG > C variant was confirmed in the proband. In addition, it was confirmed that both parents and an unaffected sister are heterozygous carriers, and the unaffected brother is homozygous for the reference allele. The variant was absent from 100 unaffected and unrelated Lebanese subjects, and it has not been reported in the medical literature or in disease-related variant databases (Human Gene Mutation Database (HGMD) or ClinVar). Characterization of wild-type UBLCP1 expression and localization in the brain To assess the relevance of UBLCP1 to ASD, its expression in the mammalian brain was examined. UBLCP1 expression was detected in early post-natal brain (post-natal day 4) (Fig. 3a) and in adult mouse brain (Fig. 3b) by immunohistochemical staining. UBLCP1 was present in all regions of the brain including cortex, diencephalon, rhombencephalon, hippocampus, cerebellum, and ventral tegmental area. This pattern of expression is consistent with the distribution of Ublcp1 mRNA as determined by in situ hybridization data from the Allen Brain Atlas project [32]. There was no difference noted between post-natal day 4 and adult mice without the impact of male and female gender on mouse brains at the two ages. Double immunolabeling was used to assess UBLCP1 cellular and subcellular localization in the adult. The presence of UBLCP1 was observed in neurons (NeuN and MAP2 positive), dopaminergic neurons [tyrosine hydroxylase (TH) positive], and Purkinje cells (calbindin positive). UBLCP1 was also expressed in the glial population, in astrocytes (GFAP positive), microglia (CD11b positive), and oligodendrocytes (Olig2 positive) (Fig. 3c). Interestingly, UBLCP1 expression in glial cells is mostly localized to the nucleolus. Neuronal cells display additional expression of UBLCP1 in the nucleoplasm (Fig. 3d). UBLCP1 colocalized with the nucleolar protein fibrillarin in neurons, a component of the small nuclear ribonucleoprotein complex required for processing of pre-rRNA molecules [33] (Fig. 3e). A weak signal was detected in the cytoplasm of some but not all cells on post-natal day 4 and in adult mice. The UBLCP1 ASD mutation (del g.158710261CAAAG > C) alters UBLCP1 levels in fibroblasts UBLCP1 ASD mutation (del g.158710261CAAAG > C) alters UBLCP1 levels in fibroblasts The described UBLCP1 frameshift mutation leads to a premature stop codon. The impact on mRNA levels was assessed by qRT-PCR. Analysis revealed significantly reduced UBLCP1 mRNA levels in fibroblasts from the affected individual compared to unaffected controls (Fig. 4a). Immunostaining of fibroblasts showed that UBLCP1 localized mostly to the nucleus in control fibroblasts and exhibited a diffuse distribution in the cytoplasm (Fig. 4d, upper panel), whereas mutant UBLCP1 had a marked reduction in expression in fibroblasts from the ASD patient (Fig. 4d, lower panel). Quantification showed an 80% and 40% decrease in fluorescence of mutant protein within the nucleus and cytoplasm, respectively, in ASD patient fibroblasts (Fig. 4b, c). Proteasome activity is enhanced in fibroblasts with the UBLCP1 (del g.158710261CAAAG > C) mutation UBLCP1 (del g.158710261CAAAG > C) mutation The ASD mutation identified in family F2 is a deletion within exon 10 of UBLCP1 (g. 158,710,261CAAAG > C) and is predicted to generate a stop codon that interrupts the phosphatase domain of the protein. To study the effect of this mutation on the UPS, proteasome activity was measured in fibroblasts using a fluorogenic synthetic peptide substrate, Suc-LLVY-AMC, which is converted to a fluorescent degradation product, 7-amino-4-methylcoumarin (AMC) upon cleavage by the proteasome. Fibroblasts from the proband displayed a statistically significant increase in the rate of Suc-LLVY-AMC cleavage compared to fibroblasts from unaffected controls (Fig. 4e). In addition, Western blot analysis uncovered a statistically significant decrease in the level of ubiquitinated proteins in fibroblasts from the proband compared to those from unaffected controls (Fig. 4f). Differential expression of 26 S proteasome subunits In addition to proteolysis in cells, the UPS impacts the regulation of the cell cycle, signal transduction, apoptosis, and antigen presentation, as well as other cellular functions. Proteasome homeostasis is tightly regulated. Mammalian cells respond to impairment of cellular proteasomal enzyme activity through a coordinated upregulation of proteasome subunits at both the transcriptional and translational levels. The hypothesis is that defective UBLCP1 protein leads to decreased inhibition of proteasome activity. The effect of the UBLCP1 ASD mutation on RNA expression of different proteasome subunits ( PSMA1 , PSMB5 , PSMC2 , PSMC4 , and PSMD2 ) in patient-derived fibroblasts was investigated by qRT-PCR. A general tendency for reduced RNA expression of all the subunits in the proband’s cells compared to expression in cells from unaffected controls was observed. This reduction was statistically significant for PSMA1 ( P  = 0.0318), PSMC4 ( P  = 0.0291), and PSMD2 ( P  = 0.0165) (Fig. 4h). When proband cells were treated with MG132, a proteasome inhibitor, the mRNA levels of the 5 subunits increased significantly compared to DMSO treated cells ( PSMA1 ( P  = 0.0053), PSMB5 ( P  = 0.0128), PSMC2 ( P  = 0.0084), PSMC4 ( P  = 0.0124), and PSMD2 ( P  = 0.0369) (Fig. 4i). Although mutant UBLCP1 mRNA levels were slightly increased, this difference did not reach statistical significance (P = 0.1096). The effect of proteasome inhibition with MG132 on wild-type UBLCP1 gene expression was determined in control fibroblasts. UBLCP1 mRNA levels were significantly increased in treated cells as compared to DMSO-treated cells (P = 0.0384) (Fig. 4g). Gentamicin treatment restored expression and function of mutant UBLCP1 in ASD patient-derived fibroblasts Gentamicin, an aminoglycoside antibiotic, can promote premature termination codons (PTC) read-through, restoring full-length protein without inducing significant misreading at normal stop codons [34]. Untreated fibroblasts generated from the ASD patient revealed significantly reduced UBLCP1 protein expression in the nucleus and in cytoplasm compared to control fibroblasts. Treatment of patient fibroblasts with gentamicin restored UBLCP1 expression in the nucleus and cytoplasm (Fig. 5a). Quantification revealed that UBLCP1 expression was significantly reduced in the nucleus and the cytoplasm to 22.4% ± 1.4% and 83.8% ± 1.2% of levels observed in control fibroblasts, respectively. Treatment with gentamicin upregulated UBLCP1 expression in nucleus and cytoplasm to 79.5 ± 3.8% and 125.6 ± 4.7% of levels observed in control fibroblasts, respectively (Fig. 5b, c). To verify if UBLCP1 proteasome inhibitor function was restored, proteasome activity was investigated after treatment with gentamicin. Fibroblasts generated from the ASD patient displayed a statistically significant decrease in the rate of Suc-LLVY-AMC cleavage after treatment with gentamicin, confirming restored function (Fig. 5d). Discussion In this study, whole exome sequencing of a consanguineous Lebanese family with a boy affected with ASD, led to the identification of a novel candidate gene, UBLCP1 . The gene encodes a protein with a phosphatase domain implicated in the regulation of the proteasome. Under normal conditions, UBLCP1 functions to inhibit or downregulate proteasome activity. Histological analysis in the mouse revealed that UBLCP1 is ubiquitously expressed early during postnatal development and in the adult brain. This consolidates the idea that UBLCP1 plays an essential role in cellular functioning. A defective UBLCP1 may cause dysregulation in cellular homeostasis leading to altered development and function of the central nervous system. Ubiquitous expression in neuronal and glial cell types cannot but have a role to play in the pathogenesis of ASD. Functional analysis of fibroblasts bearing truncated UBLCP1 revealed an increase in proteasome activity, and subsequent decrease in polyubiquitinated protein levels. This data confirms the hypothesis that mutant UBLCP1, in this instance fails to phosphorylate the proteasome, leading to decreased inhibition, hence, increased proteasomal activation. Dysfunction of the UPS is observed in many diseases, including disorders of the brain [30]. Numerous mutations in genes encoding proteins within the UPS are linked to neurodegenerative and neurodevelopmental disorders, including ASD [35, 36]. Guo et al. showed that UBLCP1 localizes to the nucleus due to its interaction with nuclear proteasomes via the Lys44 residue within the UBL domain [25]. In this study, UBLCP1 was detected primarily in the nucleus, but also at lower levels in the cytoplasm. These results are in line with observations of proteasome location in the nucleolus and the cytoplasm [37]. Despite the opinion that proteasomes degrade proteins inside the nucleus, recent studies diverge and report that catalytically active proteasomes exist almost exclusively in the cytosol [38]. Here, there is a 40% decrease in expression of mutant UBLCP1 in the nucleus. The hypothesis here is that wild-type UBLCP1, by binding and dephosphorylating the proteasome decreases its assembly, and contributes to its nuclear localization in a dissociated state. In the case of mutant UBLCP1 truncated within the phosphatase domain, it binds to the proteasome via the UBL domain but does not dephosphorylate it, hence, resulting in increased assembly and trafficking of active proteasome to the cytoplasm. In yeast and plants, the coordinated expression of proteasome subunit genes is controlled by unrelated but functionally analogous transcription factors, Rpn4 [39, 40] and NAC53/78 [41], respectively. In mammals, there are different analogous transcription factors. NRF1 (Nuclear Factor Erythroid-derived 2-Related Factor 1) is a central regulator of proteasome expression in the presence of proteasomal dysfunction [39, 40]. Under normal conditions, NRF1 is a type II integral ER membrane protein continuously retro-translocated from the ER back to the cytosol via the ER-associated protein degradation (ERAD) pathway, where it is rapidly ubiquitinated and degraded by the proteasome [39, 40, 42]. With proteotoxic stress, NRF1 is proteolytically released from the ER, enters the nucleus, and binds to antioxidant response element (ARE) core sequences (5’-TGA[C/T]NNNGC-3’) [43, 44] found upstream of most proteasome subunit genes. In this study, overactivation of the proteasome induced by mutant UBLCP1 is countered by downregulation of different subunits of the proteasome and UBLCP1. Here, it was determined that ARE sequences exist in the promotor of all subunits genes investigated as previously described in the literature [45] and also within the promotor of UBLCP1 (5’-TGACACAGC-3’ at ~360 bp from the start site). Overactivation of the proteasome in cells due to defective de-phosphorylation by mutant UBLCP1 leads to degradation of a common transcription factor in the cytosol, thus leading to a coordinated downregulation of levels of core proteasome subunits and associated regulators, such as UBLCP1, as a protective mechanism (Fig. 6). This was demonstrated by the upregulation of these genes observed upon inhibition of the proteasome in proband cells. The significant increase in UBLCP1 mRNA levels observed in treated control cells attests to the fact that UBLCP1 expression is regulated by proteasomal activity. This, however, did not translate to a statistically significant result in proband cells. A potential reason for this observation could be nonsense-mediated decay of the truncated UBLCP1 transcript. This truncation could be leading to a reduced half-life of the transcript, which masks the increase in gene expression after treatment in proband cells. Western blot analysis revealed an increase in NRF1 protein levels after treatment with MG132, as previously reported by Redhakrishnan et al. in 2010 (see Supplementary Fig. 1). While inhibition of the proteasome is widely investigated, very few studies report overactivation of proteasome function, and little is known about the ensuing impact on cells. During the acute phase of intracerebral hemorrhage, overactivation of the proteasome causes early degradation of the endoplasmic reticulum (ER) chaperone GRP78 and IκB proteins, which in turn activate NF-κB increasing neuroinflammation via expression of pro-inflammatory cytokines. According to Liew et al., overactivation of proteasome activity exacerbates ER stress [46]. A number of studies suggest an important role of ER stress in the pathophysiology of autism [47, 48]. The expression of ER stress-related genes ATF4 , ATF6 , PERK , XBP1 , CHOP , and IRE1 is significantly increased in the middle frontal gyrus of ASD subjects [49]. Fujita et al. [50] state that low levels of ER stress might alter membrane trafficking of synaptic functional molecules possibly leading to ASD pathophysiology. A missense mutation in SLC6A1 , encoding GABA transporter 1 (GAT-1) in a patient with epilepsy and ASD was recently described with resultant relocation of the mutated protein to the ER. The increased ER retention, along with excessive ER-associated degradation, led to reduced GABA uptake function in this case [51]. Emerging evidence suggests that ASD pathophysiology may be linked to dysfunction of GABAergic inhibitory signaling in the brain. One of the GABAA receptor subunits, GABAAα1, exhibits a significant reduction at the protein level in the frontal cortex of ASD subjects [52]. Using postmortem brain samples and mouse primary cortical neurons, Crider et al. demonstrate a critical role for UPS-mediated mechanisms in the post-translational regulation of GABAAα1. Conversely, expression of the E3 ligase SYVN1 associated with GABAAα1 was higher in tissue samples from ASD compared to controls [49]. Overall, excessive UPS-mediated GABAAα1 turnover in ASD might lead to dysfunction of GABA signaling in cortical neurons. These findings, while promising, remain distanced from the context of ASD due to tissue specificity differences between human fibroblasts and neurons. Nonetheless, it establishes the framework for future work in which iPSCs would be generated from control and proband fibroblasts prior to differentiation to neuronal and glial cellular subtypes. Only then can one assess the ensuing effects on neuronal and glial pathways. Emerging evidence suggests that proteasomes are relatively enriched in the brain and may play critical roles in neuronal physiology [53, 54]. A defective regulator such as ubiquitously expressed UBLCP1 may impact various aspects of neuronal development, including neurogenesis, differentiation, and synapse formation. It could also affect numerous post-developmental functions, such as synaptic plasticity and signaling. Other limitations of the current study include the use of male subjects only. This hinders the assessment of sex-specific differences in proteasome activity, which have been previously reported [55, 56]. This is the first discovery of a mutation leading to overactivation of the proteasome in ASD. Overall, the present work suggests that under or overactivation of proteostasis may cause ASD. Compensatory changes necessary to maintain a viable rate of proteolysis probably contribute to the neurological symptoms of ASD. Neutralizing the effect of this mutation by using gentamicin to promote read-through of the premature stop codon in UBLCP1 mRNA corrected the deficient phenotype in cells. This result is a prelude to using the less toxic drug, Ataluren, a small non-aminoglycoside molecule, as a safe and effective treatment for this patient, and perhaps additional cases of autism with premature stop codons in other genes. Data availability Raw whole exome sequencing data from this study are available through the NCBI Sequence Read Archive (BioProject accession number: PRJNA1047123). VCFs are available upon request from Dr. Maria Chahrour (maria.chahrour@utsouthwestern.edu). References Lai MC, Lombardo MV, Baron-Cohen S. Autism. Lancet. 2014;383:896–910. De Rubeis S, Buxbaum JD. Genetics and genomics of autism spectrum disorder: embracing complexity. Hum Mol Genet. 2015;24:R24–31. Sandin S, Lichtenstein P, Kuja-Halkola R, Hultman C, Larsson H, Reichenberg A. The heritability of autism spectrum disorder. JAMA. 2017;318:1182–4. Bailey A, Le Couteur A, Gottesman I, Bolton P, Simonoff E, Yuzda E, et al. Autism as a strongly genetic disorder: evidence from a British twin study. Psychol Med. 1995;25:63–77. Ozonoff S, Young GS, Carter A, Messinger D, Yirmiya N, Zwaigenbaum L, et al. Recurrence risk for autism spectrum disorders: a Baby Siblings Research Consortium study. Pediatrics. 2011;128:e488–495. Loomes R, Hull L, Mandy WPL. What is the male-to-female ratio in autism spectrum disorder? A systematic review and meta-analysis. J Am Acad Child Adolesc Psychiatry. 2017;56:466–74. Chaaya M, Saab D, Maalouf FT, Boustany RM. Prevalence of autism spectrum disorder in nurseries in Lebanon: a cross sectional study. J Autism Dev Disord. 2016;46:514–22. Saab DC, M; Boustany, RM. Prevalence & correlates of autism spectrum disorder: Lebanese cross-sectional study. Autism 2018;8. Gerges P, Bitar T, Laumonnier F, Marouillat S, Nemer G, Andres CR, et al. Identification of novel gene variants for autism spectrum disorders in the lebanese population using whole-exome sequencing. Genes. 2022;13:186. Iossifov I, O’Roak BJ, Sanders SJ, Ronemus M, Krumm N, Levy D, et al. The contribution of de novo coding mutations to autism spectrum disorder. Nature. 2014;515:216–21. Jimenez-Barron LT, O’Rawe JA, Wu Y, Yoon M, Fang H, Iossifov I, et al. Genome-wide variant analysis of simplex autism families with an integrative clinical-bioinformatics pipeline. Cold Spring Harb Mol Case Stud. 2015;1:a000422. Barbour B, Salameh P. Consanguinity in Lebanon: prevalence, distribution and determinants. J Biosoc Sci. 2009;41:505–17. El-Sibai M, Platt DE, Haber M, Xue Y, Youhanna SC, Wells RS, et al. Geographical structure of the Y-chromosomal genetic landscape of the Levant: a coastal-inland contrast. Ann Hum Genet. 2009;73:568–81. Yu TW, Chahrour MH, Coulter ME, Jiralerspong S, Okamura-Ikeda K, Ataman B, et al. Using whole-exome sequencing to identify inherited causes of autism. Neuron. 2013;77:259–73. Chahrour MH, Yu TW, Lim ET, Ataman B, Coulter ME, Hill RS, et al. Whole-exome sequencing and homozygosity analysis implicate depolarization-regulated neuronal genes in autism. PLoS Genet. 2012;8:e1002635. Glessner JT, Wang K, Cai G, Korvatska O, Kim CE, Wood S, et al. Autism genome-wide copy number variation reveals ubiquitin and neuronal genes. Nature. 2009;459:569–73. Kishino T, Lalande M, Wagstaff J. UBE3A/E6-AP mutations cause Angelman syndrome. Nat Genet. 1997;15:70–73. O’Roak BJ, Deriziotis P, Lee C, Vives L, Schwartz JJ, Girirajan S, et al. Exome sequencing in sporadic autism spectrum disorders identifies severe de novo mutations. Nat Genet. 2011;43:585–9. Piton A, Gauthier J, Hamdan FF, Lafreniere RG, Yang Y, Henrion E, et al. Systematic resequencing of X-chromosome synaptic genes in autism spectrum disorder and schizophrenia. Mol Psychiatry. 2011;16:867–80. Khalil R, Kenny C, Hill RS, Mochida GH, Nasir R, Partlow JN, et al. PSMD12 haploinsufficiency in a neurodevelopmental disorder with autistic features. Am J Med Genet Part B, Neuropsychiatr Genet. 2018;177:736–45. Finley D. Recognition and processing of ubiquitin-protein conjugates by the proteasome. Annu Rev Biochem. 2009;78:477–513. Lee MJ, Lee BH, Hanna J, King RW, Finley D. Trimming of ubiquitin chains by proteasome-associated deubiquitinating enzymes. Mol Cell Proteomics. 2011;10:R110 003871. Zhang F, Hu Y, Huang P, Toleman CA, Paterson AJ, Kudlow JE. Proteasome function is regulated by cyclic AMP-dependent protein kinase through phosphorylation of Rpt6. J Biol Chem. 2007;282:22460–71. Bingol B, Wang CF, Arnott D, Cheng D, Peng J, Sheng M. Autophosphorylated CaMKIIalpha acts as a scaffold to recruit proteasomes to dendritic spines. Cell. 2010;140:567–78. Guo X, Engel JL, Xiao J, Tagliabracci VS, Wang X, Huang L, et al. UBLCP1 is a 26S proteasome phosphatase that regulates nuclear proteasome activity. Proc Natl Acad Sci USA. 2011;108:18649–54. Yun JH, Ko S, Lee CK, Cheong HK, Cheong C, Yoon JB, et al. Solution structure and Rpn1 interaction of the UBL domain of human RNA polymerase II C-terminal domain phosphatase. PloS One. 2013;8:e62981. Wang X, Huang L. Identifying dynamic interactors of protein complexes by quantitative mass spectrometry. Mol Cell Proteomics. 2008;7:46–57. Liu X, Xiao W, Zhang Y, Wiley SE, Zuo T, Zheng Y, et al. Reversible phosphorylation of Rpn1 regulates 26S proteasome assembly and function. Proc Natl Acad Sci USA. 2020;117:328–36. Sun S, Liu S, Zhang Z, Zeng W, Sun C, Tao T, et al. Phosphatase UBLCP1 controls proteasome assembly. Open Biol. 2017;7:170042. Zheng Q, Huang T, Zhang L, Zhou Y, Luo H, Xu H, et al. Dysregulation of ubiquitin-proteasome system in neurodegenerative diseases. Front Aging Neurosci. 2016;8:303. Van der Auwera GA, Carneiro MO, Hartl C, Poplin R, Del Angel G, Levy-Moonshine A, et al. From FastQ data to high confidence variant calls: the Genome Analysis Toolkit best practices pipeline. Curr Protoc Bioinformatics. 2013;43:11 10 11–11 10 33. Allen Institute for Brain Science. Allen Mouse Brain Atlas [Dataset]. Available from mouse.brain-map.org (2004). Ochs RL, Lischwe MA, Spohn WH, Busch H. Fibrillarin: a new protein of the nucleolus identified by autoimmune sera. Biol Cell. 1985;54:123–33. Bidou L, Allamand V, Rousset JP, Namy O. Sense from nonsense: therapies for premature stop codon diseases. Trends Mol Med. 2012;18:679–88. Ciechanover A, Brundin P. The ubiquitin proteasome system in neurodegenerative diseases: sometimes the chicken, sometimes the egg. Neuron. 2003;40:427–46. Louros SR, Osterweil EK. Perturbed proteostasis in autism spectrum disorders. J Neurochem. 2016;139:1081–92. Wojcik C, DeMartino GN. Intracellular localization of proteasomes. Int J Biochem Cell Biol. 2003;35:579–89. Kulichkova VA, Artamonova TO, Zaykova JJ, Ermolaeva JB, Khodorkovskii MA, Barlev NA, et al. Simultaneous EGFP and tag labeling of the beta7 subunit for live imaging and affinity purification of functional human proteasomes. Mol Biotechnol. 2015;57:36–44. Radhakrishnan SK, den Besten W, Deshaies RJ. p97-dependent retrotranslocation and proteolytic processing govern formation of active Nrf1 upon proteasome inhibition. eLife. 2014;3:e01856. Sha Z, Goldberg AL. Proteasome-mediated processing of Nrf1 is essential for coordinate induction of all proteasome subunits and p97. Curr Biol. 2014;24:1573–83. Gladman NP, Marshall RS, Lee KH, Vierstra RD. The proteasome stress regulon is controlled by a pair of NAC transcription factors in arabidopsis. Plant Cell. 2016;28:1279–96. Steffen J, Seeger M, Koch A, Kruger E. Proteasomal degradation is transcriptionally controlled by TCF11 via an ERAD-dependent feedback loop. Mol Cell. 2010;40:147–58. Baird L, Tsujita T, Kobayashi EH, Funayama R, Nagashima T, Nakayama K, et al. A homeostatic shift facilitates endoplasmic reticulum proteostasis through transcriptional integration of proteostatic stress response pathways. Mol Cell Biol. 2017;37:e00439–16. Chepelev NL, Bennitz JD, Huang T, McBride S, Willmore WG. The Nrf1 CNC-bZIP protein is regulated by the proteasome and activated by hypoxia. PloS One. 2011;6:e29167. Koizumi S, Hamazaki J, Murata S. Transcriptional regulation of the 26S proteasome by Nrf1. Proc Jpn Acad Ser B Phys Biol Sci. 2018;94:325–36. Liew HK, Hu WF, Lin PB, Wang PK, Tsai AP, Pang CY, et al. Over-activated proteasome mediates neuroinflammation on acute intracerebral hemorrhage in rats. Cells. 2019;8:1326. Etherton M, Foldy C, Sharma M, Tabuchi K, Liu X, Shamloo M, et al. Autism-linked neuroligin-3 R451C mutation differentially alters hippocampal and cortical synaptic function. Proc Natl Acad Sci USA. 2011;108:13764–9. Momoi T, Fujita E, Senoo H, Momoi M. Genetic factors and epigenetic factors for autism: endoplasmic reticulum stress and impaired synaptic function. Cell Biol Int. 2009;34:13–19. Crider A, Ahmed AO, Pillai A. Altered expression of endoplasmic reticulum stress-related genes in the middle frontal cortex of subjects with autism spectrum disorder. Mol Neuropsychiatry. 2017;3:85–91. Fujita E, Dai H, Tanabe Y, Zhiling Y, Yamagata T, Miyakawa T, et al. Autism spectrum disorder is related to endoplasmic reticulum stress induced by mutations in the synaptic cell adhesion molecule, CADM1. Cell Death Dis. 2010;1:e47. Wang J, Poliquin S, Mermer F, Eissman J, Delpire E, Wang J, et al. Endoplasmic reticulum retention and degradation of a mutation in SLC6A1 associated with epilepsy and autism. Mol Brain. 2020;13:76. Fatemi SH, Folsom TD, Reutiman TJ, Thuras PD. Expression of GABA(B) receptors is altered in brains of subjects with autism. Cerebellum. 2009;8:64–9. He HY, Ahsan A, Bera R, McLain N, Faulkner R, Ramachandran KV, et al. Neuronal membrane proteasomes regulate neuronal circuit activity in vivo and are required for learning-induced behavioral plasticity. Proc Natl Acad Sci USA. 2023;120:e2216537120. Turker F, Cook EK, Margolis SS. The proteasome and its role in the nervous system. Cell Chem Biol. 2021;28:903–17. Beamish SB, Gross KS, Anderson MM, Helmstetter FJ, Frick KM. Sex differences in training-induced activity of the ubiquitin proteasome system in the dorsal hippocampus and medial prefrontal cortex of male and female mice. Learn Mem. 2022;29:302–11. Jenkins EC, Shah N, Gomez M, Casalena G, Zhao D, Kenny TC, et al. Proteasome mapping reveals sexual dimorphism in tissue-specific sensitivity to protein aggregations. EMBO Rep. 2020;21:e48978. Ashkenazy H, Abadi S, Martz E, Chay O, Mayrose I, Pupko T, et al. ConSurf 2016: an improved methodology to estimate and visualize evolutionary conservation in macromolecules. Nucleic Acids Res. 2016;44:W344–350. Acknowledgements We thank families and staff from the Lebanese Autism Society (LAS), SESOBEL, and the AUBMC Special Kids Clinic for their help on this project. The work is supported by generous grants to R-MB from both OpenMinds (grant number 620229) and the AUBMC Medical Practice Plan (grant number 320046). Ethics declarations Competing interests The authors declare no competing interests. Additional information Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Supplementary information Rights and permissions Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/. About this article Cite this article Soueid, J., Hamze, Z., Bedran, J. et al. A novel autism-associated UBLCP1 mutation impacts proteasome regulation/activity.
Transl Psychiatry 13 , 404 (2023). https://doi.org/10.1038/s41398-023-02702-0 Received: Revised: Accepted: Published: DOI: https://doi.org/10.1038/s41398-023-02702-0