Abstract Background Autism spectrum disorder (ASD), a developmental disorder of early childhood onset, affects males four times more frequently than females, suggesting a role for the sex chromosomes. In this study, we describe a family with ASD in which a predicted pathogenic nonsense mutation in the X-chromosome gene RAB39B segregates with ASD phenotype. Methods Clinical phenotyping, microarray, and whole genome sequencing (WGS) were performed on the five members of this family. Maternal and female sibling X inactivation ratio was calculated, and phase was investigated. Mutant-induced pluripotent stem cells engineered for an exon 2 nonsense mutation were generated and differentiated into cortical neurons for expression and pathway analyses. Results Two males with an inherited RAB39B mutation both presented with macrocephaly, intellectual disability (ID), and ASD. Their female sibling with the same mutation presented with ID and a broad autism phenotype. In contrast, their transmitting mother has no neurodevelopmental diagnosis. Our investigation of phase indicated maternal preferential inactivation of the mutated allele, with no such bias observed in the female sibling. We offer the explanation that this bias in X inactivation may explain the absence of a neurocognitive phenotype in the mother. Our cellular knockout model of RAB39B revealed an impact on expression in differentiated neurons for several genes implicated in brain development and function, supported by our pathway enrichment analysis. Conclusions Penetrance for ASD is high among males but more variable among females with RAB39B mutations. A critical role for this gene in brain development and function is demonstrated. Electronic supplementary material The online version of this article (10.1186/s13229-017-0175-3) contains supplementary material, which is available to authorized users. Keywords: RAB39B, Intellectual disability (ID), RNAseq, Whole genome sequencing (WGS) Background Genomic microarray technology and genome sequencing are changing our understanding of the genetic architecture of autism spectrum disorder (ASD). This lifelong neurodevelopmental disorder of early childhood onset [[51]1] is now known to be associated with rare inherited and de novo genetic variation, comprising single nucleotide variants (SNV), smaller and larger insertions and deletions (indels and copy number variants, CNVs, respectively), and other complex structural variation [[52]2–[53]4]. ASD affects males on average four times more frequently than females, which strongly suggests an etiological role for the sex chromosomes. However, the majority of genes implicated in ASD, so far, are autosomal, with notable exceptions including MeCP2, the NLGN genes, and FMR1. With an estimated 400 or more genes involved in ASD’s etiology [[54]5], there may remain several undiscovered ASD susceptibility genes on the sex chromosome. For example, some X-linked genes involved in intellectual disability (ID) may be contributory factors in ASD, particularly as more than 100 genes on the X chromosome have now been described in association with ID [[55]6]. Our own whole genome sequencing (WGS) study ([56]www.mss.ng), currently comprising 2620 ASD genomes from 2066 unique families, has identified a number of X chromosome genes in which rare predicted damaging mutations are enriched, including MECP2, NLGN3, NLGN4, and PCDH11X [[57]7]. One such gene, RAB39B, has been described in the literature in association with X-linked ID along with variable phenotypic manifestations including ASD, seizures, macrocephaly, delayed psychomotor development, and early-onset Parkinson’s disease or Parkinsonism [[58]8–[59]11]. This two-exon gene, which comprises two protein domains, is a member of the RAS oncogene family with known expression in neurons and neuronal precursors, particularly in the hippocampus [[60]8]. Along with other RAS proteins, it plays a key role in intracellular vesicular trafficking. Specifically, through its trafficking role, RAB39B is now known to mediate the surface expression of GluA2, a subunit of the glutamate AMPA receptor [[61]12]. Here, we present a detailed description of a family in which two brothers with ASD have a maternally inherited loss of function (LoF) mutation in exon 2 of RAB39B. Both males are now in their early 30s and have ID with no functional language; they also have fine motor difficulties and macrocephaly. Additionally, one presents with a childhood onset unilateral hand tremor and has marked slowing in the execution of simple, routinized movements. In addition, their female sibling, who shares the RAB39B mutation, presents with mild ID and a broad autism phenotype (BAP), but not ASD. Finally, the transmitting mother has a unilateral fine upper intention tremor of unknown etiology. We examined the impact of an exon 2 mutation on gene expression in differentiated neurons using CRISPR/Cas. Methods The family described was recruited as part of ongoing studies investigating the genetic etiology of ASD [[62]7]. ASD diagnoses are made by expert clinicians using the Autism Diagnostic Interview (ADI) [[63]13] and the Autism Diagnostic Observation Schedule (ADOS) [[64]14] combined with clinical judgment. All data were collected following informed consent from participants or substitute decision-makers, and the study was conducted with approval from the local research ethics board. The mother has provided specific written consent for this case report. Phenotypes Both affected males were recruited into the study at the same time, one aged 9 years (III-3, hereafter “proband”) and the other 8 years (III-4) (Fig. [65]1). They had both been diagnosed with ASD at age 3 years, and on recruitment into the study, diagnostic clarification was sought with ADI-R and ADOS-G. Additionally, nonverbal IQ was measured using the Leiter-R, language using the Oral and Written Language Scales (OWLS-II), and adaptive function using the Vineland Adaptive Behavior Scales-II (VABS). Their female sibling (III-5) was first assessed in a specialist pediatric clinic because of developmental concerns age 5 years. At the time of recruitment (age 6 years), she completed a measure of IQ (Wechsler Intelligence Scale for Children) and her adaptive function was measured using the VABS. Assessment for ASD was not completed during her childhood, but age 29 years the presence of autism traits of childhood onset were ascertained by a maternally completed Social Responsiveness Scale (SRS) [[66]15]. Head circumference was also collected, and craniofacial features were evaluated. Medical history was collected by way of a detailed questionnaire. The two males with ASD were re-evaluated age 33 years (III-3) and 32 years (III-4). The Leiter-R and VABS were both repeated, and interim medical history was collected. Both parents completed the Toronto Alexithymia Scale [[67]16]. This scale measures alexithymia, characterized by an inability to self-reflect on and label emotions. Although neuropsychiatric and developmental history was not specifically sought for either parent, both are professionally employed and neither endorsed an ASD diagnosis, or broader phenotype, during the collection of family history. Fig. 1. Fig. 1 [68]Open in a new tab Annotated pedigree of family described in the text Genotypes Members of the family were genotyped using either an Affymetrix 500k (II-3, II-4, and III-3) [[69]17] and/or Illumina 1M (II-3, II-4, and III-4) [[70]18, [71]19] microarray, and CNVs were detected as previously published [[72]18–[73]20]. Briefly, multiple different CNV calling algorithms were used to generate high-confidence CNV calls, with support from a minimum of two algorithms required, and CNVs with a frequency of < 0.1% in controls (N = 9611) prioritized as “rare.” CNVs with > 50% reciprocal overlap were deemed identical. We also removed all CNVs that had > 50% overlap with a known segmental duplication. No CNVs of clinical significance, or of uncertain clinical significance, were identified in this family. The proband, his male and female siblings, and both parents also underwent whole genome sequencing (WGS) by Complete Genomics (Mountain View, CA) as previously published [[74]21]. Identified variants were annotated and likely deleterious mutations prioritized to capture those that were rare (MAF ≤ 1%), LoF (nonsense, splice site, and frameshift), and damaging de novo missense mutations (damaging as evidenced by SIFT ≤ 0.05, Polyphen2 ≥ 0.95, CADD ≥ 15, Mutation Assessor score ≥ 2 and PhyloP ≥ 2.4) [[75]21]. For filtering, MAF was based on the following control samples: 1000 Genomes Project, NHLBI-ESP, ExAC, and Complete Genomics whole genome sequencing control data (CGI diversity panel and Wellderly samples). Variants reported in this paper have been validated by Sanger sequencing. Extended family Several members of the extended family (Fig. [76]1) were genotyped for the RAB39B mutation. Although these individuals were unavailable for neuropsychiatric evaluation, phenotypic data was collected through a family history review provided by the mother (II-3). X inactivation The X inactivation (XI) ratio was calculated by characterizing genotyped alleles according to the number of CAG repeats on the highly polymorphic androgen receptor at Xq11-q12. The methylation-sensitive restriction enzymes Hpa II and Hha I were used to cut the repeat region on the active (but not inactive) strand before and after quantitative PCR to determine the X inactivation ratio. Non-methylated (active X) DNA segments digest with the enzymes and are thereby unavailable for PCR amplification. Methylated (inactive X) Hpa II/Hha I sites do not digest with enzyme and remain intact for amplification. Post-digestion PCR products therefore represent methylated (inactive X) DNA sequences only. The androgen receptor was also genotyped on paternal and male offspring X-chromosomes to facilitate phasing of the mutated allele. This thereby allowed us to determine whether the mutated or non-mutated allele was preferentially inactivated in mother (II-3) and daughter (III-5). Investigation of familial segregation We estimated likely causality for the family’s RAB39B mutation, and for all variants in RAB39B described in the literature, using the method of Thompson et al. [[77]22, [78]23]. This considers the segregation of the variant across the full pedigree from which a Bayes factor (BF) is calculated. As articulated in the new ACMG-AMP guidelines for variant classification [[79]24], a BF of 8 or more provides “supportive” evidence of pathogenicity, of 16 or more “moderate” evidence of pathogenicity, and of 32 or more “strong” evidence of pathogenicity. iPS culture and differentiation Skin fibroblasts from a healthy man were obtained at the Hospital for Sick Children under the approval of the SickKids Research Ethics Board, and induced pluripotent stem cells (iPSCs) were generated and characterized as previously published [[80]25]. Cultured iPSCs were maintained on Matrigel® (Corning). One million iPSCs were nucleofected using Nucleofector II (Amaxa) device with program A-023 in Solution I with 5 μg Cas9D10A plasmid (Addgene, #44720), 3 μg each of two gRNA cloning plasmids (Addgene, #41824) containing paired gRNAs A (AAGAGGTTGTCAAATCAGAGAGG) and B (TCTGAAGAGTGAACCACATTTGG) targeting exon 2 of RAB39B, and 1 μl of 100 μM single-stranded oligonucleotide (ssODN) containing two 60 nucleotide-long arms of homology flanking a 3xSTOP, i.e., stop codons in all three reading frames. This introduced mutation is not the same as that involved in the family described. Cells were evenly distributed in a 96-well plate. At confluency, genomic DNA was extracted and droplet digital PCR (ddPCR) was performed to determine the absolute quantification of mutant alleles in each well. The mutant stem cell clone was then differentiated into cortical excitatory neurons. A comparison cell line was generated from the isogenic control line, i.e., the same line from which the knockout described above was made. RNAseq RNeasy® mini kit (Qiagen) was used to extract total RNA from mature neurons. RNA samples were submitted to an Agilent Bioanalyzer 2100 RNA Nano chip for quality control. Concentration was determined using Qubit RNA HS Assay on a Qubit fluorometer (ThermoFisher). RNA libraries were prepared using NEBNext Ultra RNA Library Preparation kit for Illumina. Briefly, 500 ng of total RNA was used for poly-A mRNA enrichment before being split into 200–300 bp fragments for 4 min at 94 °C. Fragments were converted to double-stranded cDNA, end-repaired, and adenylated in 3′ to create an overhang A, allowing ligation of Illumina adapters with an overhang T. Fragments were amplified and RNA library integrity was verified on a Bioanalyzer 2100 DNA high-sensitivity chip (Agilent Technologies) and quantified using Kapa Library Quantification Illumina/ABI Prism Kit protocol (KAPA Biosystems). Stranded libraries were pooled in equimolar amounts and sequenced on an Illumina HiSeq 2500 platform using a high-throughput run mode flowcell and the V4 sequencing chemistry to generate paired-end reads of 126-bases in length, following Illumina’s recommended protocol. Data quality was assessed using FastQC v.0.11.2. Trimmed reads were screened for presence of rRNA and mtRNA sequences using FastQ-Screen v.0.4.3. RSeQC package v.2.3.7 was used to assess read distribution and positional read duplication and confirm strandedness of alignments. Raw trimmed reads were aligned to the reference genome hg19 using Tophat v.2.0.11. Tophat alignments were processed to extract raw read counts for genes using htseq-count v.0.6.1p2. Raw gene counts were loaded and sample-normalized using DESeq v.1.18.0. Principal component analysis (PCA) was performed to assess relation among samples. Two-condition differential expression was done with the edgeR R package, v.3.8.6 on four (RAB39B knockout) and six (control) independent experiments. Cpm filter was adjusted to the size of each library set to obtain at least 8–10 reads in at least two samples. EdgeR results were evaluated using FDR with a cutoff of 0.05. Pathway enrichment Pathway enrichment analysis was then performed using the R package goseq version 1.24.0 and R version 3.3.1 (2016-06-21) using a custom gene set collection including 6008 from gene ontology (GO, obtained from the R package GO.db version 3.3.0) and 1723 from pathways (KEGG, Reactome, NCI collections downloaded from the respective websites on 20170120). GO terms were obtained from the R package [81]org.Hs.eg.db v3.4.0. Gene sets were filtered by size to remove both very generic and very specific terms. To aid in reproducibility, the filtered gene set definitions are in Additional file [82]1: human_GO_15_800_PATH_10_500.gmt. Results Phenotypes The proband (III-3, Fig. [83]1) was evaluated age 3 years because of parental concerns regarding speech delay and poor socialization and diagnosed with ASD. Parents reported that his speech was developing well up until 18 months with the acquisition of a small number of well-articulated words, but then plateaued. At 36 months, he presented as a loner who had no interest in communicating or interacting with others, and he failed to use eye-to-eye contact. Gross and fine motor milestones developed within normal limits, and there were no other medical concerns. Morphological abnormalities included a head circumference on the 98th percentile and finger syndactyly. By 8 years, his language had improved, although was comprised predominantly of echolalia. Assessment of intellectual ability and adaptive functioning age 9 years were consistent with a diagnosis of ID, and ADI-R [[84]13] and ADOS [[85]14] assessments completed at this time were both consistent with a diagnosis of ASD (Table [86]1 and Additional file [87]2: Table S1). He was further assessed age 33 years. Interim medical history was negative for major medical or psychiatric illness, and re-evaluation at this time revealed very little change in symptoms. A number of movement-related symptoms were noted, including impaired fine motor control and episodic difficulties initiating movements. He was also noted to have a fine bilateral tremor of his hands, worse when agitated or anxious. Table 1. Genetic and clinical data from family segregating a RAB39B mutation II-3 II-4 III-3 III-4 III-5 Age at first assessment (months) NA NA 36 36 60 Sex F M M M F Microarray^a NS NS NS NS NS WGS^b Xq28, RAB39B, g.154,490,171C>A 3p26.3, CNTN6, g.1,443,180G>A WT RAB39B WT CNTN6 Xq28, RAB39B, g.154,490,171C>A 3p26.3, CNTN6, g.1,443,180G>A Xq28, RAB39B, g.154,490,171C>A 3p26.3, CNTN6 g.1,443,180G>A Xq28, RAB39B, g.154,490,171C>A WT CNTN6 Growth  Head circumference (%ile) NA NA 98th 98th 50th Neurodevelopment  Full-scale IQ 110 109 39 32 69  Speech delay NA NA + + +  ASD − − + + +  Other neurodevelopmental NA NA − − − Neurological  Epilepsy − − − − −  Parkinson − − − − −  Other neurological Tremor NA Tremor, bradykinesia, poor fine motor Poor fine motor −  Congenital NA NA Syndactyly Syndactyly, epicanthus, high-arched palate Left epicanthus  Other medical NA NA − − − [88]Open in a new tab NA information not available, +/− positive/negative for attribute ^a NS implies no variants of pathological significance, or variants of unknown significance contributing to ASD identified ^b WT wild type The proband’s male sibling (III-4) was evaluated age 3 years because of parental concerns with language delay and diagnosed with ASD. At 36 months, he used no functional language and was socially aloof. His motor milestones were normal. He was on the 98th percentile for head circumference and had epicanthic folds, a high-arched palate, and a finger syndactyly. By 6 years, he had developed some language skills, although failed to use this to facilitate interaction with others. Assessment of intellectual ability and adaptive functioning age 8 years were consistent with a diagnosis of ID, and ADI-R and ADOS at the time were both consistent with ASD (Table [89]1). When re-evaluated age 32 years, he had some limited language skills although his nonverbal IQ and adaptive skills had changed little from childhood. An interim medical history did not identify any major medical or psychiatric illness. No movement-related symptoms were identified, although he, too, had major difficulties with fine motor control, impacting feeding, dressing, and other instrumental activities of daily living. Their younger female sibling (III-5) was first evaluated when she was 5 years because of concerns with speech delay. She was described as responsive and curious as an infant and became securely attached to her primary caregivers. However, her speech developed late, and articulation and intonation abnormalities were noted. She was described as aloof from her peers and content to play by herself. Assessment of intellectual ability and adaptive functioning age 6 years were consistent with a diagnosis of ID (Table [90]1). Morphologically, she was on the 50th percentile for head circumference at the time of assessment and had a left-sided epicanthus. Although she was not assessed for ASD during childhood, a maternally completed SRS [[91]15] when she was 29 years indicated a moderate degree of impairment, consistent with BAP. She was unavailable for further diagnostic testing. Both parents completed a range of assessments that indicated intellectual ability and language functioning in the average range. Both parents’ scores on the Toronto Alexithymia Scale [[92]16] were consistent with alexithymia (father = 70, mother = 60, scores 61 and above being consistent with alexithymia). At the time of assessment, parents were in full-time professional employment. Maternal self-report indicated a diagnosis of a persistent left-sided intention tremor of unknown etiology, but no other neurological history. Members of the extended family were unavailable for phenotypic evaluation (Fig. [93]1). However, review of the family history according to the mother was negative for ASD, other developmental disorders, neuropsychiatric diagnoses, or major medical illness in extended family members. Genetics A nonsense variant in exon 2 of RAB39B (Xq28, g.154,490,171C>A, c.559G>T, [p.E187X], [94]NM_171998, CADD=12.46) was identified in both males with ASD, inherited from their mother. The female sibling was found to be heterozygous for this variant by WGS and Sanger sequencing. The variant was absent from all members of the extended family who were tested (Fig. [95]1). Analysis of X inactivation in the mother and daughter showed the mother had preferential inactivation of the mutated allele (67%:33%, mutated allele: non-mutated allele), whereas her daughter showed no such bias (48%:52%, mutated allele: non-mutated alleles). RAB39B is given a pLI score of 0.72 in the Exome Aggregation Consortium (ExAC), consistent with intolerance to mutation [[96]26]. No individuals with this mutation are listed in ExAC. A Bayes factor calculated from the pattern of segregation of this mutation in the family is consistent with moderate to strong evidence of pathogenicity according to the ACMG-AMP guidelines (BF = 26.1) [[97]24]. The two affected males also inherited from their mother a nonsense variant in CNTN6 (3p26.3, g.1,443,180G>A [p.W923X], [98]NM_001289080, CADD=13.85), a gene implicated in ASD. This variant is observed in one non-Finnish European individual in ExAC, which assigns this gene a pLI score of 0, suggesting greater tolerance to mutation than RAB39B. The frequency of this mutation in ExAC is 8.44e-06. This mutation was not shared by their female sibling. However, one maternal sibling reported to be normal from a neurodevelopmental perspective also had this mutation (II-7). No clinically significant CNVs were identified. We investigated our own data for evidence of RAB39B mutations in other individuals with ASD (N = 2620) ([99]www.mss.ng). One male was identified with a maternally inherited 13 kb deletion (chrX:154,477,501-154,490,661 [hg19]). He is described as exhibiting classic symptoms of ASD with speech delay (age of single words = 78 months). No further information is available. We also investigated for the presence of CNVs and SNVs in CNTN6. Among cases with ASD, five individuals have exon impacting SNVs and 15 exon impacting CNVs (3 deletions and 12 duplications). No individuals were identified with variants impacting both CNTN6 and RAB39B. iPSC and expression studies We used CRISPR/Cas to create a mutation in exon 2 of RAB39B. Although the mutation generated in this way is not identical to the family’s, it does provide insight into the impact of an exon 2 loss of function mutation on differential patterns of gene expression. A number of genes were differentially expressed between mutant and control lines (Fig. [100]2a, b; Additional file [101]3: RAB39B4_vs_Controls6_RPKM_Filt.txt). Figure [102]2 specifically highlights those differentially expressed genes that were previously found associated with ASD and are present in the Simons Foundation Autism Research Initiative (SFARI) database. Among these top hits are known SFARI syndromic genes (HCN1, C12orf57) and one high confidence ASD-risk gene (GRIN2B). There are also two genes, GABRA4 and GABRB2, that code for subunits of the GABA receptor, the major inhibitory neurotransmission unit in the mammalian brain. Moreover, several genes encode elements of the potassium channel (KCNIP2, KCNK13, KCNIP2). Finally, NXPH1 encodes a protein that forms complexes with alpha-neurexins, implicating its role in cellular adhesion in the brain. We next conducted gene set enrichment analysis to identify the associated GO terms with the 260 differentially expressed genes with FDR < 0.05 (Additional file [103]4: enrichment_RAB39B4_vs_Controls6.txt). Only downregulated genes were significant, with their GO terms provided in the Additional file [104]5: Table S2. The enriched terms point to the involvement of these downregulated genes in neuron fate, neurogenesis, and regulation of nervous system development. Fig. 2. Fig. 2 [105]Open in a new tab Levels of transcript of selected genes in control and RAB39B-/γ-glutaminergic neurons as determined by RNAseq experiment. Values are presented as mean ± SD of four (RAB39B) and six (control) independent experiments We curated all of the descriptions of RAB39B patients from the literature (Fig. [106]3). Five families have been described, each with a unique variant [[107]8, [108]9, [109]11], along with one further individual identified from a sample of 1348 patients with Parkinson’s disease (PD) [[110]10]. In one family (Fig. [111]3, Reference 1), six males with ID and macrocephaly are described with a RAB39B splice donor variant (hg19: 154,493,358, c.215+1G>A) [[112]8]. Two males are also diagnosed with ASD. In another family (Fig. [113]3, Reference 2), four males with ID, three of whom were also diagnosed with epilepsy, inherited a RAB39B nonsense mutation in exon 1 (hg19: 154,493,553, c.21C>A [p.Tyr7X]) [[114]8]. In two other families (Fig. [115]3, References 3 and 4), macrocephaly, ID, and Parkinsonian symptoms are