Alternative splicing is prevalent in the mammalian brain. To interrogate the functional role of alternative splicing in neural development, we analyzed purified neural progenitor cells (NPCs) and neurons from developing cerebral cortices, revealing hundreds of differentially spliced exons that preferentially alter key protein domains-especially in cytoskeletal proteins-and can harbor disease-causing mutations. We show that Ptbp1 and Rbfox proteins antagonistically govern the NPC-to-neuron transition by regulating neuron-specific exons. Whereas Ptbp1 maintains apical progenitors partly through suppressing a poison exon of Flna in NPCs, Rbfox proteins promote neuronal differentiation by switching Ninein from a centrosomal splice form in NPCs to a non-centrosomal isoform in neurons. We further uncover an intronic human mutation within a PTBP1-binding site that disrupts normal skipping of the FLNA poison exon in NPCs and causes a brain-specific malformation. Our study indicates that dynamic control of alternative splicing governs cell fate in cerebral cortical development.

The genetic architecture of autism spectrum disorder involves the interplay of common and rare variants and their impact on hundreds of genes. Using exome sequencing, here we show that analysis of rare coding variation in 3,871 autism cases and 9,937 ancestry-matched or parental controls implicates 22 autosomal genes at a false discovery rate (FDR) < 0.05, plus a set of 107 autosomal genes strongly enriched for those likely to affect risk (FDR < 0.30). These 107 genes, which show unusual evolutionary constraint against mutations, incur de novo loss-of-function mutations in over 5% of autistic subjects. Many of the genes implicated encode proteins for synaptic formation, transcriptional regulation and chromatin-remodelling pathways. These include voltage-gated ion channels regulating the propagation of action potentials, pacemaking and excitability-transcription coupling, as well as histone-modifying enzymes and chromatin remodellers-most prominently those that mediate post-translational lysine methylation/demethylation modifications of histones.

FSHD2 is a rare form of facioscapulohumeral muscular dystrophy (FSHD) characterized by the absence of a contraction in the D4Z4 macrosatellite repeat region on chromosome 4q35 that is the hallmark of FSHD1. However, hypomethylation of this region is common to both subtypes. Recently, mutations in SMCHD1 combined with a permissive 4q35 allele were reported to cause FSHD2. We identified a novel p.Lys275del SMCHD1 mutation in a family affected with FSHD2 using whole-exome sequencing and linkage analysis. This mutation alters a highly conserved amino acid in the ATPase domain of SMCHD1. Subject III-11 is a male who developed asymmetrical muscle weakness characteristic of FSHD at 13 years. Physical examination revealed marked bilateral atrophy at biceps brachii, bilateral scapular winging, some asymmetrical weakness at tibialis anterior and peroneal muscles, and mild lower facial weakness. Biopsy of biceps brachii in subject II-5, the father of III-11, demonstrated lobulated fibers and dystrophic changes. Endomysial and perivascular inflammation was found, which has been reported in FSHD1 but not FSHD2. Given the previous report of SMCHD1 mutations in FSHD2 and the clinical presentations consistent with the FSHD phenotype, we conclude that the SMCHD1 mutation is the likely cause of the disease in this family.

The current study characterizes a cohort of limb-girdle muscular dystrophy (LGMD) in the United States using whole-exome sequencing. Fifty-five families affected by LGMD were recruited using an institutionally approved protocol. Exome sequencing was performed on probands and selected parental samples. Pathogenic mutations and cosegregation patterns were confirmed by Sanger sequencing. Twenty-two families (40%) had novel and previously reported pathogenic mutations, primarily in LGMD genes, and also in genes for Duchenne muscular dystrophy, facioscapulohumeral muscular dystrophy, congenital myopathy, myofibrillar myopathy, inclusion body myopathy and Pompe disease. One family was diagnosed via clinical testing. Dominant mutations were identified in COL6A1, COL6A3, FLNC, LMNA, RYR1, SMCHD1 and VCP, recessive mutations in ANO5, CAPN3, GAA, LAMA2, SGCA and SGCG, and X-linked mutations in DMD. A previously reported variant in DMD was confirmed to be benign. Exome sequencing is a powerful diagnostic tool for LGMD. Despite careful phenotypic screening, pathogenic mutations were found in other muscle disease genes, largely accounting for the increased sensitivity of exome sequencing. Our experience suggests that broad sequencing panels are useful for these analyses because of the phenotypic overlap of many neuromuscular conditions. The confirmation of a benign DMD variant illustrates the potential of exome sequencing to help determine pathogenicity.

Although autism has a clear genetic component, the high genetic heterogeneity of the disorder has been a challenge for the identification of causative genes. We used homozygosity analysis to identify probands from nonconsanguineous families that showed evidence of distant shared ancestry, suggesting potentially recessive mutations. Whole-exome sequencing of 16 probands revealed validated homozygous, potentially pathogenic recessive mutations that segregated perfectly with disease in 4/16 families. The candidate genes (UBE3B, CLTCL1, NCKAP5L, ZNF18) encode proteins involved in proteolysis, GTPase-mediated signaling, cytoskeletal organization, and other pathways. Furthermore, neuronal depolarization regulated the transcription of these genes, suggesting potential activity-dependent roles in neurons. We present a multidimensional strategy for filtering whole-exome sequence data to find candidate recessive mutations in autism, which may have broader applicability to other complex, heterogeneous disorders.

Adult mesenchymal stem cells, including preadipocytes, possess a cellular sensory organelle called the primary cilium. Ciliated preadipocytes abundantly populate perivascular compartments in fat and are activated by a high-fat diet. Here, we sought to understand whether preadipocytes use their cilia to sense and respond to external cues to remodel white adipose tissue. Abolishing preadipocyte cilia in mice severely impairs white adipose tissue expansion. We discover that TULP3-dependent ciliary localization of the omega-3 fatty acid receptor FFAR4/GPR120 promotes adipogenesis. FFAR4 agonists and ω-3 fatty acids, but not saturated fatty acids, trigger mitosis and adipogenesis by rapidly activating cAMP production inside cilia. Ciliary cAMP activates EPAC signaling, CTCF-dependent chromatin remodeling, and transcriptional activation of PPARγ and CEBPα to initiate adipogenesis. We propose that dietary ω-3 fatty acids selectively drive expansion of adipocyte numbers to produce new fat cells and store saturated fatty acids, enabling homeostasis of healthy fat tissue.

Proteins related to the ubiquitin-proteasome system play an important role during the differentiation and ciliogenesis of photoreceptor cells. Mutations in several genes involved in ubiquitination and proteostasis have been identified as causative of inherited retinal dystrophies (IRDs) and ciliopathies. USP48 is a deubiquitinating enzyme whose role in the retina is still unexplored although previous studies indicate its relevance for neurosensory organs. In this work, we describe that a pool of endogenous USP48 localises to the basal body in retinal cells and provide data that supports the function of USP48 in the photoreceptor cilium. We also demonstrate that USP48 interacts with the IRD-associated proteins ARL3 and UNC119a, and stabilise their protein levels using different mechanisms. Our results suggest that USP48 may act in the regulation/stabilisation of key ciliary proteins for photoreceptor function, in the modulation of intracellular protein transport, and in ciliary trafficking to the photoreceptor outer segment.

Cells inherit two centrioles, the older of which is uniquely capable of generating a cilium. Using proteomics and superresolved imaging, we identify a module that we term DISCO (distal centriole complex). The DISCO components CEP90, MNR, and OFD1 underlie human ciliopathies. This complex localizes to both distal centrioles and centriolar satellites, proteinaceous granules surrounding centrioles. Cells and mice lacking CEP90 or MNR do not generate cilia, fail to assemble distal appendages, and do not transduce Hedgehog signals. Disrupting the satellite pools does not affect distal appendage assembly, indicating that it is the centriolar populations of MNR and CEP90 that are critical for ciliogenesis. CEP90 recruits the most proximal known distal appendage component, CEP83, to root distal appendage formation, an early step in ciliogenesis. In addition, MNR, but not CEP90, restricts centriolar length by recruiting OFD1. We conclude that DISCO acts at the distal centriole to support ciliogenesis by restraining centriole length and assembling distal appendages, defects in which cause human ciliopathies.

During the development of the cerebral cortex, neurons are generated directly from radial glial cells or indirectly via basal progenitors. The balance between these division modes determines the number and types of neurons formed in the cortex thereby affecting cortical functioning. Here, we investigate the role of primary cilia in controlling the decision between forming neurons directly or indirectly. We show that a mutation in the ciliary gene Inpp5e leads to a transient increase in direct neurogenesis and subsequently to an overproduction of layer V neurons in newborn mice. Loss of Inpp5e also affects ciliary structure coinciding with reduced Gli3 repressor levels. Genetically restoring Gli3 repressor rescues the decreased indirect neurogenesis in Inpp5e mutants. Overall, our analyses reveal how primary cilia determine neuronal subtype composition of the cortex by controlling direct versus indirect neurogenesis. These findings have implications for understanding cortical malformations in ciliopathies with INPP5E mutations.

Craniofacial defects are among the most common phenotypes caused by ciliopathies, yet the developmental and molecular etiology of these defects is poorly understood. We investigated multiple mouse models of human ciliopathies (including Tctn2, Cc2d2a, and Tmem231 mutants) and discovered that each displays hypotelorism, a narrowing of the midface. As early in development as the end of gastrulation, Tctn2 mutants displayed reduced activation of the Hedgehog (HH) pathway in the prechordal plate, the head organizer. This prechordal plate defect preceded a reduction of HH pathway activation and Shh expression in the adjacent neurectoderm. Concomitant with the reduction of HH pathway activity, Tctn2 mutants exhibited increased cell death in the neurectoderm and facial ectoderm, culminating in a collapse of the facial midline. Enhancing HH signaling by decreasing the gene dosage of a negative regulator of the pathway, Ptch1, decreased cell death and rescued the midface defect in both Tctn2 and Cc2d2a mutants. These results reveal that ciliary HH signaling mediates communication between the prechordal plate and the neurectoderm to provide cellular survival cues essential for development of the facial midline.

There is tremendous potential for genome sequencing to improve clinical diagnosis and care once it becomes routinely accessible, but this will require formalizing research methods into clinical best practices in the areas of sequence data generation, analysis, interpretation and reporting. The CLARITY Challenge was designed to spur convergence in methods for diagnosing genetic disease starting from clinical case history and genome sequencing data. DNA samples were obtained from three families with heritable genetic disorders and genomic sequence data were donated by sequencing platform vendors. The challenge was to analyze and interpret these data with the goals of identifying disease-causing variants and reporting the findings in a clinically useful format. Participating contestant groups were solicited broadly, and an independent panel of judges evaluated their performance.

Primary cilia are required for Smoothened to transduce vertebrate Hedgehog signals, but how Smoothened accumulates in cilia and is activated is incompletely understood. Here, we identify cilia-associated oxysterols that promote Smoothened accumulation in cilia and activate the Hedgehog pathway. Our data reveal that cilia-associated oxysterols bind to two distinct Smoothened domains to modulate Smoothened accumulation in cilia and tune the intensity of Hedgehog pathway activation. We find that the oxysterol synthase HSD11β2 participates in the production of Smoothened-activating oxysterols and promotes Hedgehog pathway activity. Inhibiting oxysterol biosynthesis impedes oncogenic Hedgehog pathway activation and attenuates the growth of Hedgehog pathway-associated medulloblastoma, suggesting that targeted inhibition of Smoothened-activating oxysterol production may be therapeutically useful for patients with Hedgehog-associated cancers.

Cilia have a unique diffusion barrier ("gate") within their proximal region, termed transition zone (TZ), that compartmentalises signalling proteins within the organelle. The TZ is known to harbour two functional modules/complexes (Meckel syndrome [MKS] and Nephronophthisis [NPHP]) defined by genetic interaction, interdependent protein localisation (hierarchy), and proteomic studies. However, the composition and molecular organisation of these modules and their links to human ciliary disease are not completely understood. Here, we reveal Caenorhabditis elegans CEP-290 (mammalian Cep290/Mks4/Nphp6 orthologue) as a central assembly factor that is specific for established MKS module components and depends on the coiled coil region of MKS-5 (Rpgrip1L/Rpgrip1) for TZ localisation. Consistent with a critical role in ciliary gate function, CEP-290 prevents inappropriate entry of membrane-associated proteins into cilia and keeps ARL-13 (Arl13b) from leaking out of cilia via the TZ. We identify a novel MKS module component, TMEM-218 (Tmem218), that requires CEP-290 and other MKS module components for TZ localisation and functions together with the NPHP module to facilitate ciliogenesis. We show that TZ localisation of TMEM-138 (Tmem138) and CDKL-1 (Cdkl1/Cdkl2/Cdkl3/Cdlk4 related), not previously linked to a specific TZ module, similarly depends on CEP-290; surprisingly, neither TMEM-138 or CDKL-1 exhibit interdependent localisation or genetic interactions with core MKS or NPHP module components, suggesting they are part of a distinct, CEP-290-associated module. Lastly, we show that families presenting with Oral-Facial-Digital syndrome type 6 (OFD6) have likely pathogenic mutations in CEP-290-dependent TZ proteins, namely Tmem17, Tmem138, and Tmem231. Notably, patient fibroblasts harbouring mutated Tmem17, a protein not yet ciliopathy-associated, display ciliogenesis defects. Together, our findings expand the repertoire of MKS module-associated proteins--including the previously uncharacterised mammalian Tmem80--and suggest an MKS-5 and CEP-290-dependent assembly pathway for building a functional TZ.

The Meckel syndrome (MKS) complex functions at the transition zone, located between the basal body and axoneme, to regulate the localization of ciliary membrane proteins. We investigated the role of Tmem231, a two-pass transmembrane protein, in MKS complex formation and function. Consistent with a role in transition zone function, mutation of mouse Tmem231 disrupts the localization of proteins including Arl13b and Inpp5e to cilia, resulting in phenotypes characteristic of MKS such as polydactyly and kidney cysts. Tmem231 and B9d1 are essential for each other and other complex components such as Mks1 to localize to the transition zone. As in mouse, the Caenorhabditis elegans orthologue of Tmem231 localizes to and controls transition zone formation and function, suggesting an evolutionarily conserved role for Tmem231. We identified TMEM231 mutations in orofaciodigital syndrome type 3 (OFD3) and MKS patients that compromise transition zone function. Thus, Tmem231 is critical for organizing the MKS complex and controlling ciliary composition, defects in which cause OFD3 and MKS.

While next-generation sequencing has accelerated the discovery of human disease genes, progress has been largely limited to the "low hanging fruit" of mutations with obvious exonic coding or canonical splice site impact. In contrast, the lack of high-throughput, unbiased approaches for functional assessment of most noncoding variants has bottlenecked gene discovery. We report the integration of transcriptome sequencing (RNA-seq), which surveys all mRNAs to reveal functional impacts of variants at the transcription level, into the gene discovery framework for a unique human disease, microcephaly-micromelia syndrome (MMS). MMS is an autosomal recessive condition described thus far in only a single First Nations population and causes intrauterine growth restriction, severe microcephaly, craniofacial anomalies, skeletal dysplasia, and neonatal lethality. Linkage analysis of affected families, including a very large pedigree, identified a single locus on Chromosome 21 linked to the disease (LOD > 9). Comprehensive genome sequencing did not reveal any pathogenic coding or canonical splicing mutations within the linkage region but identified several nonconserved noncoding variants. RNA-seq analysis detected aberrant splicing in DONSON due to one of these noncoding variants, showing a causative role for DONSON disruption in MMS. We show that DONSON is expressed in progenitor cells of embryonic human brain and other proliferating tissues, is co-expressed with components of the DNA replication machinery, and that Donson is essential for early embryonic development in mice as well, suggesting an essential conserved role for DONSON in the cell cycle. Our results demonstrate the utility of integrating transcriptomics into the study of human genetic disease when DNA sequencing alone is not sufficient to reveal the underlying pathogenic mutation.

Genes associated with human microcephaly, a condition characterized by a small brain, include critical regulators of proliferation, cell fate and DNA repair. We describe a syndrome of congenital microcephaly and diverse defects in cerebral cortical architecture. Genome-wide linkage analysis in two families identified a 7.5-Mb locus on chromosome 19q13.12 containing 148 genes. Targeted high throughput sequence analysis of linked genes in each family yielded > 4,000 DNA variants and implicated a single gene, WDR62, as harboring potentially deleterious changes. We subsequently identified additional WDR62 mutations in four other families. Magnetic resonance imaging and postmortem brain analysis supports important roles for WDR62 in the proliferation and migration of neuronal precursors. WDR62 is a WD40 repeat-containing protein expressed in neuronal precursors as well as in postmitotic neurons in the developing brain and localizes to the spindle poles of dividing cells. The diverse phenotypes of WDR62 suggest it has central roles in many aspects of cerebral cortical development.

Despite significant heritability of autism spectrum disorders (ASDs), their extreme genetic heterogeneity has proven challenging for gene discovery. Studies of primarily simplex families have implicated de novo copy number changes and point mutations, but are not optimally designed to identify inherited risk alleles. We apply whole-exome sequencing (WES) to ASD families enriched for inherited causes due to consanguinity and find familial ASD associated with biallelic mutations in disease genes (AMT, PEX7, SYNE1, VPS13B, PAH, and POMGNT1). At least some of these genes show biallelic mutations in nonconsanguineous families as well. These mutations are often only partially disabling or present atypically, with patients lacking diagnostic features of the Mendelian disorders with which these genes are classically associated. Our study shows the utility of WES for identifying specific genetic conditions not clinically suspected and the importance of partial loss of gene function in ASDs.

We have undertaken a genome-wide analysis of rare copy-number variation (CNV) in 1124 autism spectrum disorder (ASD) families, each comprised of a single proband, unaffected parents, and, in most kindreds, an unaffected sibling. We find significant association of ASD with de novo duplications of 7q11.23, where the reciprocal deletion causes Williams-Beuren syndrome, characterized by a highly social personality. We identify rare recurrent de novo CNVs at five additional regions, including 16p13.2 (encompassing genes USP7 and C16orf72) and Cadherin 13, and implement a rigorous approach to evaluating the statistical significance of these observations. Overall, large de novo CNVs, particularly those encompassing multiple genes, confer substantial risks (OR = 5.6; CI = 2.6-12.0, p = 2.4 × 10(-7)). We estimate there are 130-234 ASD-related CNV regions in the human genome and present compelling evidence, based on cumulative data, for association of rare de novo events at 7q11.23, 15q11.2-13.1, 16p11.2, and Neurexin 1.

The transcription factor NRF2 is a master regulator of cellular antioxidant and detoxification responses, but it also regulates other processes such as autophagy and pluripotency. In human embryonic stem cells (hESCs), NRF2 antagonizes neuroectoderm differentiation, which only occurs after NRF2 is repressed via a Primary Cilia-Autophagy-NRF2 (PAN) axis. However, the functional connections between NRF2 and primary cilia, microtubule-based plasma membrane protrusions that function as cellular antennae, remain poorly understood. For instance, nothing is known about whether NRF2 affects cilia, or whether cilia regulation of NRF2 extends beyond hESCs. Here, we show that NRF2 and primary cilia reciprocally regulate each other. First, we demonstrate that fibroblasts lacking primary cilia have higher NRF2 activity, which is rescued by autophagy-activating mTOR inhibitors, indicating that the PAN axis also operates in differentiated cells. Furthermore, NRF2 controls cilia formation and function. NRF2-null cells grow fewer and shorter cilia and display impaired Hedgehog signaling, a cilia-dependent pathway. These defects are not due to increased oxidative stress or ciliophagy, but rather to NRF2 promoting expression of multiple ciliogenic and Hedgehog pathway genes. Among these, we focused on GLI2 and GLI3, the transcription factors controlling Hh pathway output. Both their mRNA and protein levels are reduced in NRF2-null cells, consistent with their gene promoters containing consensus ARE sequences predicted to bind NRF2. Moreover, GLI2 and GLI3 fail to accumulate at the ciliary tip of NRF2-null cells upon Hh pathway activation. Given the importance of NRF2 and ciliary signaling in human disease, our data may have important biomedical implications.

Motor, sensory, and integrative activities of the brain are coordinated by a series of midline-bridging neuronal commissures whose development is tightly regulated. Here we report a new human syndrome in which these commissures are widely disrupted, thus causing clinical manifestations of horizontal gaze palsy, scoliosis, and intellectual disability. Affected individuals were found to possess biallelic loss-of-function mutations in the gene encoding the axon-guidance receptor 'deleted in colorectal carcinoma' (DCC), which has been implicated in congenital mirror movements when it is mutated in the heterozygous state but whose biallelic loss-of-function human phenotype has not been reported. Structural MRI and diffusion tractography demonstrated broad disorganization of white-matter tracts throughout the human central nervous system (CNS), including loss of all commissural tracts at multiple levels of the neuraxis. Combined with data from animal models, these findings show that DCC is a master regulator of midline crossing and development of white-matter projections throughout the human CNS.