Comprehensive whole-genome sequence analyses provide insights into the genomic architecture of cerebral palsy.

We performed whole-genome sequencing (WGS) in 327 children with cerebral palsy (CP) and their biological parents. We classified 37 of 327 (11.3%) children as having pathogenic/likely pathogenic (P/LP) variants and 58 of 327 (17.7%) as having variants of uncertain significance. Multiple classes of P/LP variants included single-nucleotide variants (SNVs)/indels (6.7%), copy number variations (3.4%) and mitochondrial mutations (1.5%). The COL4A1 gene had the most P/LP SNVs. We also analyzed two pediatric control cohorts (n = 203 trios and n = 89 sib-pair families) to provide a baseline for de novo mutation rates and genetic burden analyses, the latter of which demonstrated associations between de novo deleterious variants and genes related to the nervous system. An enrichment analysis revealed previously undescribed plausible candidate CP genes (SMOC1, KDM5B, BCL11A and CYP51A1). A multifactorial CP risk profile and substantial presence of P/LP variants combine to support WGS in the diagnostic work-up across all CP and related phenotypes.

ANK2 loss-of-function variants are associated with epilepsy, and lead to impaired axon initial segment plasticity and hyperactive network activity in hiPSC-derived neuronal networks.

PURPOSE: To characterize a novel neurodevelopmental syndrome due to loss-of-function (LoF) variants in Ankyrin 2 (ANK2), and to explore the effects on neuronal network dynamics and homeostatic plasticity in human-induced pluripotent stem cell-derived neurons. METHODS: We collected clinical and molecular data of 12 individuals with heterozygous de novo LoF variants in ANK2. We generated a heterozygous LoF allele of ANK2 using CRISPR/Cas9 in human-induced pluripotent stem cells (hiPSCs). HiPSCs were differentiated into excitatory neurons, and we measured their spontaneous electrophysiological responses using micro-electrode arrays (MEAs). We also characterized their somatodendritic morphology and axon initial segment (AIS) structure and plasticity. RESULTS: We found a broad neurodevelopmental disorder (NDD), comprising intellectual disability, autism spectrum disorders and early onset epilepsy. Using MEAs, we found that hiPSC-derived neurons with heterozygous LoF of ANK2 show a hyperactive and desynchronized neuronal network. ANK2-deficient neurons also showed increased somatodendritic structures and altered AIS structure of which its plasticity is impaired upon activity-dependent modulation. CONCLUSIONS: Phenotypic characterization of patients with de novo ANK2 LoF variants defines a novel NDD with early onset epilepsy. Our functional in vitro data of ANK2-deficient human neurons show a specific neuronal phenotype in which reduced ANKB expression leads to hyperactive and desynchronized neuronal network activity, increased somatodendritic complexity and AIS structure and impaired activity-dependent plasticity of the AIS.

Disruption of DDX53 coding sequence has limited impact on iPSC-derived human NGN2 neurons.

BACKGROUND: The X-linked PTCHD1 locus is strongly associated with autism spectrum disorder (ASD). Males who carry chromosome microdeletions of PTCHD1 antisense long non-coding RNA (PTCHD1-AS)/DEAD-box helicase 53 (DDX53) have ASD, or a sub-clinical form called Broader Autism Phenotype. If the deletion extends beyond PTCHD1-AS/DDX53 to the next gene, PTCHD1, which is protein-coding, the individuals typically have ASD and intellectual disability (ID). Three male siblings with a 90 kb deletion that affects only PTCHD1-AS (and not including DDX53) have ASD. We performed a functional analysis of DDX53 to examine its role in NGN2 neurons. METHODS: We used the clustered regularly interspaced short palindromic repeats (CRISPR) gene editing strategy to knock out DDX53 protein by inserting 3 termination codons (3TCs) into two different induced pluripotent stem cell (iPSC) lines. DDX53 CRISPR-edited iPSCs were differentiated into cortical excitatory neurons by Neurogenin 2 (NGN-2) directed differentiation. The functional differences of DDX53-3TC neurons compared to isogenic control neurons with molecular and electrophysiological approaches were assessed. RESULTS: Isogenic iPSC-derived control neurons exhibited low levels of DDX53 transcripts. Transcriptional analysis revealed the generation of excitatory cortical neurons and DDX53 protein was not detected in iPSC-derived control neurons by western blot. Control lines and DDX53-3TC neurons were active in the multi-electrode array, but no overt electrophysiological phenotype in either isogenic line was observed. CONCLUSION: DDX53-3TC mutation does not alter NGN2 neuronal function in these experiments, suggesting that synaptic deficits causing ASD are unlikely in this cell type.

CAPRIN1 haploinsufficiency causes a neurodevelopmental disorder with language impairment, ADHD and ASD.

We describe an autosomal dominant disorder associated with loss-of-function variants in the Cell cycle associated protein 1 (CAPRIN1; MIM*601178). CAPRIN1 encodes a ubiquitous protein that regulates the transport and translation of neuronal mRNAs critical for synaptic plasticity, as well as mRNAs encoding proteins important for cell proliferation and migration in multiple cell types. We identified 12 cases with loss-of-function CAPRIN1 variants, and a neurodevelopmental phenotype characterized by language impairment/speech delay (100%), intellectual disability (83%), attention deficit hyperactivity disorder (82%) and autism spectrum disorder (67%). Affected individuals also had respiratory problems (50%), limb/skeletal anomalies (50%), developmental delay (42%) feeding difficulties (33%), seizures (33%) and ophthalmologic problems (33%). In patient-derived lymphoblasts and fibroblasts, we showed a monoallelic expression of the wild-type allele, and a reduction of the transcript and protein compatible with a half dose. To further study pathogenic mechanisms, we generated sCAPRIN1+/- human induced pluripotent stem cells via CRISPR-Cas9 mutagenesis and differentiated them into neuronal progenitor cells and cortical neurons. CAPRIN1 loss caused reduced neuronal processes, overall disruption of the neuronal organization and an increased neuronal degeneration. We also observed an alteration of mRNA translation in CAPRIN1+/- neurons, compatible with its suggested function as translational inhibitor. CAPRIN1+/- neurons also showed an impaired calcium signalling and increased oxidative stress, two mechanisms that may directly affect neuronal networks development, maintenance and function. According to what was previously observed in the mouse model, measurements of activity in CAPRIN1+/- neurons via micro-electrode arrays indicated lower spike rates and bursts, with an overall reduced activity. In conclusion, we demonstrate that CAPRIN1 haploinsufficiency causes a novel autosomal dominant neurodevelopmental disorder and identify morphological and functional alterations associated with this disorder in human neuronal models.

Genetic variants in DDX53 contribute to Autism Spectrum Disorder associated with the Xp22.11 locus.

Autism Spectrum Disorder (ASD) exhibits an ~4:1 male-to-female sex bias and is characterized by early-onset impairment of social/communication skills, restricted interests, and stereotyped behaviors. Disruption of the Xp22.11 locus has been associated with ASD in males. This locus includes the three-exon PTCHD1 gene, an adjacent multi-isoform long noncoding RNA (lncRNA) named PTCHD1-AS (spanning ~1Mb), and a poorly characterized single-exon RNA helicase named DDX53 that is intronic to PTCHD1-AS. While the relationship between PTCHD1/PTCHD1-AS and ASD is being studied, the role of DDX53 has not been examined, in part because there is no apparent functional murine orthologue. Through clinical testing, here, we identified 6 males and 1 female with ASD from 6 unrelated families carrying rare, predicted-damaging or loss-of-function variants in DDX53. Then, we examined databases, including the Autism Speaks MSSNG and Simons Foundation Autism Research Initiative, as well as population controls. We identified 24 additional individuals with ASD harboring rare, damaging DDX53 variations, including the same variants detected in two families from the original clinical analysis. In this extended cohort of 31 participants with ASD (28 male, 3 female), we identified 25 mostly maternally-inherited variations in DDX53, including 18 missense changes, 2 truncating variants, 2 in-frame variants, 2 deletions in the 3' UTR and 1 copy number deletion. Our findings in humans support a direct link between DDX53 and ASD, which will be important in clinical genetic testing. These same autism-related findings, coupled with the observation that a functional orthologous gene is not found in mouse, may also influence the design and interpretation of murine-modelling of ASD.

Genomic architecture of autism from comprehensive whole-genome sequence annotation.

Fully understanding autism spectrum disorder (ASD) genetics requires whole-genome sequencing (WGS). We present the latest release of the Autism Speaks MSSNG resource, which includes WGS data from 5,100 individuals with ASD and 6,212 non-ASD parents and siblings (total n = 11,312). Examining a wide variety of genetic variants in MSSNG and the Simons Simplex Collection (SSC; n = 9,205), we identified ASD-associated rare variants in 718/5,100 individuals with ASD from MSSNG (14.1%) and 350/2,419 from SSC (14.5%). Considering genomic architecture, 52% were nuclear sequence-level variants, 46% were nuclear structural variants (including copy-number variants, inversions, large insertions, uniparental isodisomies, and tandem repeat expansions), and 2% were mitochondrial variants. Our study provides a guidebook for exploring genotype-phenotype correlations in families who carry ASD-associated rare variants and serves as an entry point to the expanded studies required to dissect the etiology in the ∼85% of the ASD population that remain idiopathic.

Genome-wide rare variant score associates with morphological subtypes of autism spectrum disorder.

Defining different genetic subtypes of autism spectrum disorder (ASD) can enable the prediction of developmental outcomes. Based on minor physical and major congenital anomalies, we categorize 325 Canadian children with ASD into dysmorphic and nondysmorphic subgroups. We develop a method for calculating a patient-level, genome-wide rare variant score (GRVS) from whole-genome sequencing (WGS) data. GRVS is a sum of the number of variants in morphology-associated coding and non-coding regions, weighted by their effect sizes. Probands with dysmorphic ASD have a significantly higher GRVS compared to those with nondysmorphic ASD (P = 0.03). Using the polygenic transmission disequilibrium test, we observe an over-transmission of ASD-associated common variants in nondysmorphic ASD probands (P = 2.9 × 10-3). These findings replicate using WGS data from 442 ASD probands with accompanying morphology data from the Simons Simplex Collection. Our results provide support for an alternative genomic classification of ASD subgroups using morphology data, which may inform intervention protocols.

Mutations in trpγ, the homologue of TRPC6 autism candidate gene, causes autism-like behavioral deficits in Drosophila.

Autism Spectrum Disorder (ASD) is characterized by impaired social communication, restricted interests, and repetitive and stereotyped behaviors. The TRPC6 (transient receptor potential channel 6) represents an ASD candidate gene under an oligogenic/multifactorial model based on the initial description and cellular characterization of an individual with ASD bearing a de novo heterozygous mutation disrupting TRPC6, together with the enrichment of disruptive TRPC6 variants in ASD cases as compared to controls. Here, we perform a clinical re-evaluation of the initial non-verbal patient, and also present eight newly reported individuals ascertained for ASD and bearing predicted loss-of-function mutations in TRPC6. In order to understand the consequences of mutations in TRPC6 on nervous system function, we used the fruit fly, Drosophila melanogaster, to show that null mutations in transient receptor gamma (trpγ; the fly gene most similar to TRPC6), cause a number of behavioral defects that mirror features seen in ASD patients, including deficits in social interactions (based on courtship behavior), impaired sleep homeostasis (without affecting the circadian control of sleep), hyperactivity in both young and old flies, and defects in learning and memory. Some defects, most notably in sleep, differed in severity between males and females and became normal with age. Interestingly, hyperforin, a TRPC6 agonist and the primary active component of the St. John's wort antidepressant, attenuated many of the deficits expressed by trpγ mutant flies. In summary, our results provide further evidence that the TRPC6 gene is a risk factor for ASD. In addition, they show that the behavioral defects caused by mutations in TRPC6 can be modeled in Drosophila, thereby establishing a paradigm to examine the impact of mutations in other candidate genes.

Analysis of recent shared ancestry in a familial cohort identifies coding and noncoding autism spectrum disorder variants.

Autism spectrum disorder (ASD) is a collection of neurodevelopmental disorders characterized by deficits in social communication and restricted, repetitive patterns of behavior or interests. ASD is highly heritable, but genetically and phenotypically heterogeneous, reducing the power to identify causative genes. We performed whole genome sequencing (WGS) in an ASD cohort of 68 individuals from 22 families enriched for recent shared ancestry. We identified an average of 3.07 million variants per genome, of which an average of 112,512 were rare. We mapped runs of homozygosity (ROHs) in affected individuals and found an average genomic homozygosity of 9.65%, consistent with expectations for multiple generations of consanguineous unions. We identified potentially pathogenic rare exonic or splice site variants in 12 known (including KMT2C, SCN1A, SPTBN1, SYNE1, ZNF292) and 12 candidate (including CHD5, GRB10, PPP1R13B) ASD genes. Furthermore, we annotated noncoding variants in ROHs with brain-specific regulatory elements and identified putative disease-causing variants within brain-specific promoters and enhancers for 5 known ASD and neurodevelopmental disease genes (ACTG1, AUTS2, CTNND2, CNTNAP4, SPTBN4). We also identified copy number variants in two known ASD and neurodevelopmental disease loci in two affected individuals. In total we identified potentially etiological variants in known ASD or neurodevelopmental disease genes for ~61% (14/23) of affected individuals. We combined WGS with homozygosity mapping and regulatory element annotations to identify candidate ASD variants. Our analyses add to the growing number of ASD genes and variants and emphasize the importance of leveraging recent shared ancestry to map disease variants in complex neurodevelopmental disorders.

Delineation of a novel neurodevelopmental syndrome associated with PAX5 haploinsufficiency.

PAX5 is a transcription factor associated with abnormal posterior midbrain and cerebellum development in mice. PAX5 is highly loss-of-function intolerant and missense constrained, and has been identified as a candidate gene for autism spectrum disorder (ASD). We describe 16 individuals from 12 families who carry deletions involving PAX5 and surrounding genes, de novo frameshift variants that are likely to trigger nonsense-mediated mRNA decay, a rare stop-gain variant, or missense variants that affect conserved amino acid residues. Four of these individuals were published previously but without detailed clinical descriptions. All these individuals have been diagnosed with one or more neurodevelopmental phenotypes including delayed developmental milestones (DD), intellectual disability (ID), and/or ASD. Seizures were documented in four individuals. No recurrent patterns of brain magnetic resonance imaging (MRI) findings, structural birth defects, or dysmorphic features were observed. Our findings suggest that PAX5 haploinsufficiency causes a neurodevelopmental disorder whose cardinal features include DD, variable ID, and/or ASD.

A recurrent SHANK3 frameshift variant in Autism Spectrum Disorder.

Autism Spectrum Disorder (ASD) is genetically complex with ~100 copy number variants and genes involved. To try to establish more definitive genotype and phenotype correlations in ASD, we searched genome sequence data, and the literature, for recurrent predicted damaging sequence-level variants affecting single genes. We identified 18 individuals from 16 unrelated families carrying a heterozygous guanine duplication (c.3679dup; p.Ala1227Glyfs*69) occurring within a string of 8 guanines (genomic location [hg38]g.50,721,512dup) affecting SHANK3, a prototypical ASD gene (0.08% of ASD-affected individuals carried the predicted p.Ala1227Glyfs*69 frameshift variant). Most probands carried de novo mutations, but five individuals in three families inherited it through somatic mosaicism. We scrutinized the phenotype of p.Ala1227Glyfs*69 carriers, and while everyone (17/17) formally tested for ASD carried a diagnosis, there was the variable expression of core ASD features both within and between families. Defining such recurrent mutational mechanisms underlying an ASD outcome is important for genetic counseling and early intervention.

Single-cell transcriptome identifies molecular subtype of autism spectrum disorder impacted by de novo loss-of-function variants regulating glial cells.

BACKGROUND: In recent years, several hundred autism spectrum disorder (ASD) implicated genes have been discovered impacting a wide range of molecular pathways. However, the molecular underpinning of ASD, particularly from the point of view of 'brain to behaviour' pathogenic mechanisms, remains largely unknown. METHODS: We undertook a study to investigate patterns of spatiotemporal and cell type expression of ASD-implicated genes by integrating large-scale brain single-cell transcriptomes (> million cells) and de novo loss-of-function (LOF) ASD variants (impacting 852 genes from 40,122 cases). RESULTS: We identified multiple single-cell clusters from three distinct developmental human brain regions (anterior cingulate cortex, middle temporal gyrus and primary visual cortex) that evidenced high evolutionary constraint through enrichment for brain critical exons and high pLI genes. These clusters also showed significant enrichment with ASD loss-of-function variant genes (p < 5.23 × 10-11) that are transcriptionally highly active in prenatal brain regions (visual cortex and dorsolateral prefrontal cortex). Mapping ASD de novo LOF variant genes into large-scale human and mouse brain single-cell transcriptome analysis demonstrate enrichment of such genes into neuronal subtypes and are also enriched for subtype of non-neuronal glial cell types (astrocyte, p < 6.40 × 10-11, oligodendrocyte, p < 1.31 × 10-09). CONCLUSION: Among the ASD genes enriched with pathogenic de novo LOF variants (i.e. KANK1, PLXNB1), a subgroup has restricted transcriptional regulation in non-neuronal cell types that are evolutionarily conserved. This association strongly suggests the involvement of subtype of non-neuronal glial cells in the pathogenesis of ASD and the need to explore other biological pathways for this disorder.

Genome-wide detection of tandem DNA repeats that are expanded in autism.

Tandem DNA repeats vary in the size and sequence of each unit (motif). When expanded, these tandem DNA repeats have been associated with more than 40 monogenic disorders1. Their involvement in disorders with complex genetics is largely unknown, as is the extent of their heterogeneity. Here we investigated the genome-wide characteristics of tandem repeats that had motifs with a length of 2-20 base pairs in 17,231 genomes of families containing individuals with autism spectrum disorder (ASD)2,3 and population control individuals4. We found extensive polymorphism in the size and sequence of motifs. Many of the tandem repeat loci that we detected correlated with cytogenetic fragile sites. At 2,588 loci, gene-associated expansions of tandem repeats that were rare among population control individuals were significantly more prevalent among individuals with ASD than their siblings without ASD, particularly in exons and near splice junctions, and in genes related to the development of the nervous system and cardiovascular system or muscle. Rare tandem repeat expansions had a prevalence of 23.3% in children with ASD compared with 20.7% in children without ASD, which suggests that tandem repeat expansions make a collective contribution to the risk of ASD of 2.6%. These rare tandem repeat expansions included previously undescribed ASD-linked expansions in DMPK and FXN, which are associated with neuromuscular conditions, and in previously unknown loci such as FGF14 and CACNB1. Rare tandem repeat expansions were associated with lower IQ and adaptive ability. Our results show that tandem DNA repeat expansions contribute strongly to the genetic aetiology and phenotypic complexity of ASD.

Segregating patterns of copy number variations in extended autism spectrum disorder (ASD) pedigrees.

Autism spectrum disorder (ASD) is a relatively common childhood onset neurodevelopmental disorder with a complex genetic etiology. While progress has been made in identifying the de novo mutational landscape of ASD, the genetic factors that underpin the ASD's tendency to run in families are not well understood. In this study, nine extended pedigrees each with three or more individuals with ASD, and others with a lesser autism phenotype, were phenotyped and genotyped in an attempt to identify heritable copy number variants (CNVs). Although these families have previously generated linkage signals, no rare CNV segregated with these signals in any family. A small number of clinically relevant CNVs were identified. Only one CNV was identified that segregated with ASD phenotype; namely, a duplication overlapping DLGAP2 in three male offspring each with an ASD diagnosis. This gene encodes a synaptic scaffolding protein, part of a group of proteins known to be pathologically implicated in ASD. On the whole, however, the heritable nature of ASD in the families studied remains poorly understood.

A large data resource of genomic copy number variation across neurodevelopmental disorders.

Copy number variations (CNVs) are implicated across many neurodevelopmental disorders (NDDs) and contribute to their shared genetic etiology. Multiple studies have attempted to identify shared etiology among NDDs, but this is the first genome-wide CNV analysis across autism spectrum disorder (ASD), attention deficit hyperactivity disorder (ADHD), schizophrenia (SCZ), and obsessive-compulsive disorder (OCD) at once. Using microarray (Affymetrix CytoScan HD), we genotyped 2,691 subjects diagnosed with an NDD (204 SCZ, 1,838 ASD, 427 ADHD and 222 OCD) and 1,769 family members, mainly parents. We identified rare CNVs, defined as those found in <0.1% of 10,851 population control samples. We found clinically relevant CNVs (broadly defined) in 284 (10.5%) of total subjects, including 22 (10.8%) among subjects with SCZ, 209 (11.4%) with ASD, 40 (9.4%) with ADHD, and 13 (5.6%) with OCD. Among all NDD subjects, we identified 17 (0.63%) with aneuploidies and 115 (4.3%) with known genomic disorder variants. We searched further for genes impacted by different CNVs in multiple disorders. Examples of NDD-associated genes linked across more than one disorder (listed in order of occurrence frequency) are NRXN1, SEH1L, LDLRAD4, GNAL, GNG13, MKRN1, DCTN2, KNDC1, PCMTD2, KIF5A, SYNM, and long non-coding RNAs: AK127244 and PTCHD1-AS. We demonstrated that CNVs impacting the same genes could potentially contribute to the etiology of multiple NDDs. The CNVs identified will serve as a useful resource for both research and diagnostic laboratories for prioritization of variants.

CNTN5-/+or EHMT2-/+human iPSC-derived neurons from individuals with autism develop hyperactive neuronal networks.

Induced pluripotent stem cell (iPSC)-derived neurons are increasingly used to model Autism Spectrum Disorder (ASD), which is clinically and genetically heterogeneous. To study the complex relationship of penetrant and weaker polygenic risk variants to ASD, 'isogenic' iPSC-derived neurons are critical. We developed a set of procedures to control for heterogeneity in reprogramming and differentiation, and generated 53 different iPSC-derived glutamatergic neuronal lines from 25 participants from 12 unrelated families with ASD. Heterozygous de novo and rare-inherited presumed-damaging variants were characterized in ASD risk genes/loci. Combinations of putative etiologic variants (GLI3/KIF21A or EHMT2/UBE2I) in separate families were modeled. We used a multi-electrode array, with patch-clamp recordings, to determine a reproducible synaptic phenotype in 25% of the individuals with ASD (other relevant data on the remaining lines was collected). Our most compelling new results revealed a consistent spontaneous network hyperactivity in neurons deficient for CNTN5 or EHMT2. The biobank of iPSC-derived neurons and accompanying genomic data are available to accelerate ASD research. Editorial note: This article has been through an editorial process in which authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).

Complete Disruption of Autism-Susceptibility Genes by Gene Editing Predominantly Reduces Functional Connectivity of Isogenic Human Neurons.

Autism spectrum disorder (ASD) is phenotypically and genetically heterogeneous. We present a CRISPR gene editing strategy to insert a protein tag and premature termination sites creating an induced pluripotent stem cell (iPSC) knockout resource for functional studies of ten ASD-relevant genes (AFF2/FMR2, ANOS1, ASTN2, ATRX, CACNA1C, CHD8, DLGAP2, KCNQ2, SCN2A, TENM1). Neurogenin 2 (NGN2)-directed induction of iPSCs allowed production of excitatory neurons, and mutant proteins were not detectable. RNA sequencing revealed convergence of several neuronal networks. Using both patch-clamp and multi-electrode array approaches, the electrophysiological deficits measured were distinct for different mutations. However, they culminated in a consistent reduction in synaptic activity, including reduced spontaneous excitatory postsynaptic current frequencies in AFF2/FMR2-, ASTN2-, ATRX-, KCNQ2-, and SCN2A-null neurons. Despite ASD susceptibility genes belonging to different gene ontologies, isogenic stem cell resources can reveal common functional phenotypes, such as reduced functional connectivity.

A genome-wide linkage study of autism spectrum disorder and the broad autism phenotype in extended pedigrees.

BACKGROUND: Although several genetic variants for autism spectrum disorder (ASD) have now been identified, these largely occur sporadically or are de novo. Much less progress has been made in identifying inherited variants, even though the disorder itself is familial in the majority of cases. The objective of this study was to identify chromosomal regions that harbor inherited variants increasing the risk for ASD using an approach that examined both ASD and the broad autism phenotype (BAP) among a unique sample of extended pedigrees. METHODS: ASD and BAP were assessed using standardized tools in 28 pedigrees from Canada and the USA, each with at least three ASD-diagnosed individuals from two nuclear families. Genome-wide linkage analysis was performed using the posterior probability of linkage (PPL) statistic, a quasi-Bayesian method that provides strength of evidence for or against linkage in an essentially model-free manner, with outcomes on the probability scale. RESULTS: The results confirm appreciable interfamilial heterogeneity as well as a high level of intrafamilial heterogeneity. Both ASD and combined ASD/BAP specific loci are apparent. CONCLUSIONS: Inclusion of subclinical phenotypes such as BAP should be more widely employed in genetic studies of ASD as a way of identifying inherited genetic variants for the disorder. Moreover, the results underscore the need for approaches to identifying genetic risk factors in extended pedigrees that are robust to high levels of inter/intrafamilial locus and allelic heterogeneity.

Mutations in RAB39B in individuals with intellectual disability, autism spectrum disorder, and macrocephaly.

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.