A pathogenic variant in RAB32 causes autosomal dominant Parkinson’s disease and activates LRRK2 kinase

Summary Background Parkinson’s disease (PD) is a progressive neurodegenerative disorder. Mendelian forms have revealed multiple genes, with a notable emphasis on membrane trafficking; RAB GTPases play an important role in PD as a subset are both regulators and substrates of LRRK2 protein kinase. To explore the role of RAB GTPases in PD, we undertook a comprehensive examination of their genetic variability in familial PD. Methods Affected probands from 130 multi-incident PD families underwent whole-exome sequencing and genotyping, Potential pathogenic variants in 61 RAB GTPases were genotyped in relatives to assess disease segregation. These variants were also genotyped in a larger case-control series, totaling 3,078 individuals (2,734 with PD). The single most significant finding was subsequently validated within genetic data (6,043 with PD). Clinical and pathologic findings were summarized for gene-identified patients, and haplotypes were constructed. In parallel, wild-type and mutant RAB GTPase structural variation, protein interactions, and resultant enzyme activities were assessed. Findings We found RAB32 c.213C>G (Ser71Arg) to co-segregate with autosomal dominant parkinsonism in three multi-incident families. RAB32 Ser71Arg was also significantly associated with PD in case-control samples: genotyping and database searches identified thirteen more patients with the same variant that was absent in unaffected controls. Notably, RAB32 Ser71Arg heterozygotes share a common haplotype. At autopsy, one patient had sparse neurofibrillary tangle pathology in the midbrain and thalamus, without Lewy body pathology. In transfected cells the RAB32 Arg71 was twice as potent as Ser71 wild type to activate LRRK2 kinase. Interpretation Our study provides unequivocal evidence to implicate RAB32 Ser71Arg in PD. Functional analysis demonstrates LRRK2 kinase activation. We provide a mechanistic explanation to expand and unify the etiopathogenesis of monogenic PD. Funding National Institutes of Health, the Canada Excellence Research Chairs program, Aligning Science Across Parkinson’s, the Michael J. Fox Foundation for Parkinson’s Research, and the UK Medical Research Council.


Pathologic diagnoses
The brain was removed within 24 hours of death, with one hemisphere frozen at 80 °C, the other fixed in formalin, and regions embedded in paraffin for dissection and histologic studies to provide a pathologic diagnosis of parkinsonism, as previously described 37 .

Whole exome sequencing (WES), variant selection and genotyping
Exonic regions were enriched using the Ion AmpliSeq exome kit (57•7Mb) and sequenced on the Ion Proton (Life Technologies, Carlsbad, CA, USA) with a minimum average coverage of 70 reads per base and an average read length of 150 bases.Reads were mapped to the NCBI Build 37•1 (hg19) reference genome using the Ion Torrent Suite 5•0.Sequences with a mapping Phred quality score under 20, fewer than 10 reads or over 95% strand bias were excluded from further analysis.Variants were annotated with ANNOVAR 1 and Combined Annotation Dependent Depletion (CADD) C-scores 2 to help rank functional, deleterious and disease causal variants.
Variants were selected for subsequent analysis if they were: (1) good quality calls, greater than 10 reads deep with balanced forward and reverse reads; (2) predicted to be SNVs or insertions/deletions (indels) in an exonic or splicing region; (3) had a minor allele frequency (MAF) of ≤ 0•01 from the genome Aggregation Database (gnomAD; http://gnomad.broadinstitute.org/),and; (4) were absent in our in-house control exomes.
Subsequent genotyping was performed by Sanger sequencing of specific exons, as previously described 3 , or using TaqMan probes (Life Technologies), following the manufacturer's instructions.

Haplotype analysis
Single nucleotide polymorphisms (SNP) were genotyped in the Tunisian families (TUN I and TUN II) using Affymetrix 500K NspI and StyI ships, and genotypes were extracted from .celintensity files using three algorithms, BBRML, JAPL and CHIAMO, and only included when there was consensus.One patient was genotyped on the HumanCyto-12_300K chip (FRA1).Otherwise, Illumina Multi-Ethnic Genome Arrays (MEGA) were used in all other familial probands and the French-Canadian family (CAN I) (Figure 1); GenomeStudio ® was used to provide genotypes for Illumina data.DNA from all available family members and samples identified with the putatively pathogenic RAB variant were genotyped, as specified above.Alternatively, genotypes for the same SNPs were retrieved from whole genome data using bcftools.Phase was established within pedigrees when possible.

Structural modelling
AlphaFold modelling 4 was used to predict the interaction between LRRK2 and RAB32.A local installation of ColabFold 5 using MMseqs2 sequence alignment and AlphaFold2-multimer-v3 model, with additional AMBER structure relaxation was employed to model interactions between the LRRK2 fragments (1-1000 or 350-550, Uniprot: Q5S007) and RAB32 (Ser71 WT and Arg71 mutant, Uniprot: Q13637).Resulting structures were visualized and analyzed using PyMOL 2•5•5.BIOVIA Discovery Studio Visualizer 2021 was used to determine possible intermolecular interactions (Non-bond Interaction Monitor) and to predict potential rotamers.

RAB32 mRNA and protein expression across tissues
For RAB32 mRNA expression we used RNA-seq data for 17,510 human samples originating from 54 different human tissues (GTEx, v8) that were downloaded using the R package recount (v 1•4•6)Click or tap here to enter text.. Cell lines, sex-specific tissues, and tissues with 10 samples or below were removed.Samples with large chromosomal deletions and duplications or large copy number variation previously associated with disease were filtered out.For RAB32 protein expression we downloaded normal tissue data from Human Protein Atlas https://www.proteinatlas.org/,accessed on 02/08/2023).In addition, we performed immunohistochemical staining in C57BL/6 mice to assess Rab32 expression in TH-expressing cells of the SNpc.Immunohistochemistry was performed, as previously describedClick or tap here to enter text., using primary rabbit anti-mouse RAB32 (ABC520, Millipore; 1:500) and chicken anti-Tyrosine Hydroxylase (ab76442, Abcam; 1:1000).Secondary antibodies used were Alexa Fluor® Goat anti-rabbit, chicken or mouse IgG (H+L) 488 and 568 secondary (1:1000).

Quantitative immunoblot analysis
A detailed description of the quantitative immunoblotting protocol has previously been described (dx.doi.org/10.17504/potocols.io.bsgrnby6).Briefly, cell lysates were mixed with a quarter of a volume of 4 x SDS-PAGE loading buffer (Invitrogen™ NuPAGE™ LDS Sample Buffer, cat# NP0007) and heated at 70 °C for 5 min.Samples were loaded onto NuPAGE 4-12% Bis-Tris Midi Gels (Thermo Fisher Scientific, Cat# WG1402BOX or Cat# WG1403BOX) and electrophoresed at 130 V for 2 hrs in NuPAGE MOPS SDS running buffer (Thermo Fisher Scientific, Cat# NP0001-02).Proteins were then electrophoretically transferred onto a nitrocellulose membrane (GE Healthcare, Amersham Protran Supported 0•45 μm NC) at 90 V for 100 min on ice in transfer buffer (48 mM Tris base and 39 mM glycine supplemented with 20% (v/v) methanol).The membranes were blocked with 5% (w/v) skim milk powder dissolved in TBS-T (50 mM Tris base, 150 mM sodium chloride (NaCl), 0•1% (v/v) Tween 20) at room temperature for 1 hr before overnight incubation at 4 °C in primary antibodies.Membranes were washed three times for 15 min each with TBS-T before being incubated with secondary antibodies for 1 hr at room temperature.Thereafter, membranes were washed with TBS-T three times with a 15-minute incubation for each wash, and protein bands were acquired via nearinfrared fluorescent detection using the LI-COR Odyssey CLx Western Blot imaging system, and intensities of bands were quantified using Image Studio Lite (version 5•2•5, RRID:SCR_013715).

Immunocytochemistry and co-localization
For experiments where GFP-RAB32 localization was examined, HeLa cells seeded on glass coverslips in 24-well plates were transfected with 1ug of cDNA using Lipofectamine 2000 (Invitrogen).The following day, coverslips were fixed, processed and imaged as previously described(16).For co-transfection of GFP-RAB32, HA-PINK1 and mCherry-Parkin, HEK293 cells were seeded on glass coverslips (Poly-D-Lysine coated) at a density of ~50,000 cells/well in 24-well plates in DMEM media without antibiotics supplemented with 10% FBS (heat inactivated).After 24 hours, cells were co-transfected using lipofectamine 3000 with plasmid constructs as follows: buffered solution (PBS) (pH 7•4) and fixed in 10% Formalin (buffered saline, pH 7•4) for 20 minutes at room temperature and permeabilized in 0•1% triton X-100 in 1x PBS for 15 minutes at room temperature with gentle shaking followed by blocking in 10% normal donkey serum in 1x PBS for 1 hour.Cells were incubated with primary mouse mAb anti-HA tag (Thermoscientific Pierce: PI26183) antibody (1:8000) overnight at 4°C.Cells were washed 3 times for 5 minutes each wash in 1x PBS.Secondary donkey anti-mouse Alexa Fluor 647 antibody (1:1000) was added to the cells and incubated for 2 hours at room temperature with gentle shaking.Primary and secondary antibodies were diluted in 0•02% triton X-100, 1% normal donkey serum in 1x PBS.Cells were washed 3 times for 5 minutes each wash in 1x PBS and mounted using Prolong Diamond anti-fade media (Life Technologies) and allowed to dry overnight before processing for confocal imaging.To assess co-localization between Rab32 and PINK1, a Leica TCS SP5 laser scanning confocal microscope (Leica Microsystem) was used.High resolution images were acquired using 63x oil-immersion (scan averaged, four times; 1024 × 1024 pixel resolution).Initial settings for pinhole, digital gain and noise reduction were optimized and all acquisition parameters were kept the same throughout the imaging.Co-localization was analyzed by Pearson's correlation co-efficient (ImageJ plug-in JACoP).12-15 cells with comparable cotransfection efficiency were selected for co-localization analysis.Statistical analysis was performed using two-way ANOVA (Graph pad Prism 8•1•1) and presented as mean ± SEM (p ≥0•0001).

Image processing and Analysis
Transfected cells were imaged under oil-immersion at 60x magnification on an Olympus FV-1000 confocal laser scanning microscope (9 x 0•33µm step size).Images were stacked using ImageJ software (NIH, USA), and masks were created to encompass the entirety of a cell, excluding the nucleus as labelled by DAPI.Masks and non-manipulated images were processed using a custom pipeline with Cell Profiler software (v•2•1•1) to identify number and density of GFP-RAB32 (-Ser71Arg)-positive structures between 1-3µm in diameter.Gametic phase for haplotypes in TUN1, TUN2 and CAN1 were determined from their pedigrees (Figure 1).Otherwise allele calls from genotyping is shown.*MAF=Minor allele frequency (dbSNP ALFA).Nevertheless, the minor allele by frequency is given for the Ref allele for variant rsIDs**.The minimal consensus haplotype in gray = 0.361 Mb