Comprehensive multi-omic profiling of somatic mutations in malformations of cortical development.

Malformations of cortical development (MCD) are neurological conditions involving focal disruptions of cortical architecture and cellular organization that arise during embryogenesis, largely from somatic mosaic mutations, and cause intractable epilepsy. Identifying the genetic causes of MCD has been a challenge, as mutations remain at low allelic fractions in brain tissue resected to treat condition-related epilepsy. Here we report a genetic landscape from 283 brain resections, identifying 69 mutated genes through intensive profiling of somatic mutations, combining whole-exome and targeted-amplicon sequencing with functional validation including in utero electroporation of mice and single-nucleus RNA sequencing. Genotype-phenotype correlation analysis elucidated specific MCD gene sets associated with distinct pathophysiological and clinical phenotypes. The unique single-cell level spatiotemporal expression patterns of mutated genes in control and patient brains indicate critical roles in excitatory neurogenic pools during brain development and in promoting neuronal hyperexcitability after birth.

Somatic mosaicism reveals clonal distributions of neocortical development.

The structure of the human neocortex underlies species-specific traits and reflects intricate developmental programs. Here we sought to reconstruct processes that occur during early development by sampling adult human tissues. We analysed neocortical clones in a post-mortem human brain through a comprehensive assessment of brain somatic mosaicism, acting as neutral lineage recorders1,2. We combined the sampling of 25 distinct anatomic locations with deep whole-genome sequencing in a neurotypical deceased individual and confirmed results with 5 samples collected from each of three additional donors. We identified 259 bona fide mosaic variants from the index case, then deconvolved distinct geographical, cell-type and clade organizations across the brain and other organs. We found that clones derived after the accumulation of 90-200 progenitors in the cerebral cortex tended to respect the midline axis, well before the anterior-posterior or ventral-dorsal axes, representing a secondary hierarchy following the overall patterning of forebrain and hindbrain domains. Clones across neocortically derived cells were consistent with a dual origin from both dorsal and ventral cellular populations, similar to rodents, whereas the microglia lineage appeared distinct from other resident brain cells. Our data provide a comprehensive analysis of brain somatic mosaicism across the neocortex and demonstrate cellular origins and progenitor distribution patterns within the human brain.

Natural Selection Shapes Codon Usage in the Human Genome.

Synonymous codon usage has been identified as a determinant of translational efficiency and mRNA stability in model organisms and human cell lines. However, whether natural selection shapes human codon content to optimize translation efficiency is unclear. Furthermore, aside from those that affect splicing, synonymous mutations are typically ignored as potential contributors to disease. Using genetic sequencing data from nearly 200,000 individuals, we uncover clear evidence that natural selection optimizes codon content in the human genome. In deriving intolerance metrics to quantify gene-level constraint on synonymous variation, we discover that dosage-sensitive genes, DNA-damage-response genes, and cell-cycle-regulated genes are particularly intolerant to synonymous variation. Notably, we illustrate that reductions in codon optimality in BRCA1 can attenuate its function. Our results reveal that synonymous mutations most likely play an underappreciated role in human variation.

Agenesis of the putamen and globus pallidus caused by recessive mutations in the homeobox gene GSX2.

Basal ganglia are subcortical grey nuclei that play essential roles in controlling voluntary movements, cognition and emotion. While basal ganglia dysfunction is observed in many neurodegenerative or metabolic disorders, congenital malformations are rare. In particular, dysplastic basal ganglia are part of the malformative spectrum of tubulinopathies and X-linked lissencephaly with abnormal genitalia, but neurodevelopmental syndromes characterized by basal ganglia agenesis are not known to date. We ascertained two unrelated children (both female) presenting with spastic tetraparesis, severe generalized dystonia and intellectual impairment, sharing a unique brain malformation characterized by agenesis of putamina and globi pallidi, dysgenesis of the caudate nuclei, olfactory bulbs hypoplasia, and anomaly of the diencephalic-mesencephalic junction with abnormal corticospinal tract course. Whole-exome sequencing identified two novel homozygous variants, c.26C>A; p.(S9*) and c.752A>G; p.(Q251R) in the GSX2 gene, a member of the family of homeobox transcription factors, which are key regulators of embryonic development. GSX2 is highly expressed in neural progenitors of the lateral and median ganglionic eminences, two protrusions of the ventral telencephalon from which the basal ganglia and olfactory tubercles originate, where it promotes neurogenesis while negatively regulating oligodendrogenesis. The truncating variant resulted in complete loss of protein expression, while the missense variant affected a highly conserved residue of the homeobox domain, was consistently predicted as pathogenic by bioinformatic tools, resulted in reduced protein expression and caused impaired structural stability of the homeobox domain and weaker interaction with DNA according to molecular dynamic simulations. Moreover, the nuclear localization of the mutant protein in transfected cells was significantly reduced compared to the wild-type protein. Expression studies on both patients' fibroblasts demonstrated reduced expression of GSX2 itself, likely due to altered transcriptional self-regulation, as well as significant expression changes of related genes such as ASCL1 and PAX6. Whole transcriptome analysis revealed a global deregulation in genes implicated in apoptosis and immunity, two broad pathways known to be involved in brain development. This is the first report of the clinical phenotype and molecular basis associated to basal ganglia agenesis in humans.

Mutations in PIK3C2A cause syndromic short stature, skeletal abnormalities, and cataracts associated with ciliary dysfunction.

PIK3C2A is a class II member of the phosphoinositide 3-kinase (PI3K) family that catalyzes the phosphorylation of phosphatidylinositol (PI) into PI(3)P and the phosphorylation of PI(4)P into PI(3,4)P2. At the cellular level, PIK3C2A is critical for the formation of cilia and for receptor mediated endocytosis, among other biological functions. We identified homozygous loss-of-function mutations in PIK3C2A in children from three independent consanguineous families with short stature, coarse facial features, cataracts with secondary glaucoma, multiple skeletal abnormalities, neurological manifestations, among other findings. Cellular studies of patient-derived fibroblasts found that they lacked PIK3C2A protein, had impaired cilia formation and function, and demonstrated reduced proliferative capacity. Collectively, the genetic and molecular data implicate mutations in PIK3C2A in a new Mendelian disorder of PI metabolism, thereby shedding light on the critical role of a class II PI3K in growth, vision, skeletal formation and neurological development. In particular, the considerable phenotypic overlap, yet distinct features, between this syndrome and Lowe's syndrome, which is caused by mutations in the PI-5-phosphatase OCRL, highlight the key role of PI metabolizing enzymes in specific developmental processes and demonstrate the unique non-redundant functions of each enzyme. This discovery expands what is known about disorders of PI metabolism and helps unravel the role of PIK3C2A and class II PI3Ks in health and disease.

Improved Pathogenic Variant Localization via a Hierarchical Model of Sub-regional Intolerance.

Different parts of a gene can be of differential importance to development and health. This regional heterogeneity is also apparent in the distribution of disease-associated mutations, which often cluster in particular regions of disease-associated genes. The ability to precisely estimate functionally important sub-regions of genes will be key in correctly deciphering relationships between genetic variation and disease. Previous methods have had some success using standing human variation to characterize this variability in importance by measuring sub-regional intolerance, i.e., the depletion in functional variation from expectation within a given region of a gene. However, the ability to precisely estimate local intolerance was restricted by the fact that only information within a given sub-region is used, leading to instability in local estimates, especially for small regions. We show that borrowing information across regions using a Bayesian hierarchical model stabilizes estimates, leading to lower variability and improved predictive utility. Specifically, our approach more effectively identifies regions enriched for ClinVar pathogenic variants. We also identify significant correlations between sub-region intolerance and the distribution of pathogenic variation in disease-associated genes, with AUCs for classifying de novo missense variants in Online Mendelian Inheritance in Man (OMIM) genes of up to 0.86 using exonic sub-regions and 0.91 using sub-regions defined by protein domains. This result immediately suggests that considering the intolerance of regions in which variants are found may improve diagnostic interpretation. We also illustrate the utility of integrating regional intolerance into gene-level disease association tests with a study of known disease-associated genes for epileptic encephalopathy.

The Burden of Candidate Pathogenic Variants for Kidney and Genitourinary Disorders Emerging From Exome Sequencing.

Background: Exome sequencing is increasingly being used for clinical diagnostics, with an impetus to expand reporting of incidental findings across a wide range of disorders. Analysis of population cohorts can help reduce risk for genetic variant misclassification and resultant unnecessary referrals to subspecialists. Objective: To examine the burden of candidate pathogenic variants for kidney and genitourinary disorders emerging from exome sequencing. Design: Secondary analysis of genetic data. Setting: A tertiary care academic medical center. Patients: A convenience sample of exome sequence data from 7974 self-declared healthy adults. Measurements: Assessment of the prevalence of candidate pathogenic variants in 625 genes associated with Mendelian kidney and genitourinary disorders. Results: Of all participants, 23.3% carried a candidate pathogenic variant, most of which were attributable to previously reported variants that had implausibly high allele frequencies. In particular, 25 genes (discovered before the creation of the Exome Aggregation Consortium, a genetic database comprising data from a large control population) accounted for 67.7% of persons with candidate pathogenic variants. After stringent filtering based on allele frequency, 1.4% of persons still had a candidate pathogenic variant, an excessive rate given the prevalence of monogenic kidney and genitourinary disorders. Manual annotation of a subset of variants showed that the majority would be classified as nonbenign under current guidelines for clinical sequence interpretation and could prompt subspecialty referrals if returned. Limitation: Limited access to health record data prevented comprehensive assessment of the phenotypic concordance with genetic diagnoses. Conclusion: Widespread reporting of incidental genetic findings related to kidney and genitourinary disorders will require stringent curation of clinical variant databases and detailed case-level review to avoid genetic misdiagnosis and unnecessary referrals. These findings motivate similar analyses for genes relevant to other medical subspecialties. Primary Funding Source: National Institute of Diabetes and Digestive and Kidney Diseases and National Human Genome Research Institute.

Diagnostic Utility of Exome Sequencing for Kidney Disease.

BACKGROUND: Exome sequencing is emerging as a first-line diagnostic method in some clinical disciplines, but its usefulness has yet to be examined for most constitutional disorders in adults, including chronic kidney disease, which affects more than 1 in 10 persons globally. METHODS: We conducted exome sequencing and diagnostic analysis in two cohorts totaling 3315 patients with chronic kidney disease. We assessed the diagnostic yield and, among the patients for whom detailed clinical data were available, the clinical implications of diagnostic and other medically relevant findings. RESULTS: In all, 3037 patients (91.6%) were over 21 years of age, and 1179 (35.6%) were of self-identified non-European ancestry. We detected diagnostic variants in 307 of the 3315 patients (9.3%), encompassing 66 different monogenic disorders. Of the disorders detected, 39 (59%) were found in only a single patient. Diagnostic variants were detected across all clinically defined categories, including congenital or cystic renal disease (127 of 531 patients [23.9%]) and nephropathy of unknown origin (48 of 281 patients [17.1%]). Of the 2187 patients assessed, 34 (1.6%) had genetic findings for medically actionable disorders that, although unrelated to their nephropathy, would also lead to subspecialty referral and inform renal management. CONCLUSIONS: Exome sequencing in a combined cohort of more than 3000 patients with chronic kidney disease yielded a genetic diagnosis in just under 10% of cases. (Funded by the National Institutes of Health and others.).

Loss of Protocadherin-12 Leads to Diencephalic-Mesencephalic Junction Dysplasia Syndrome.

OBJECTIVE: To identify causes of the autosomal-recessive malformation, diencephalic-mesencephalic junction dysplasia (DMJD) syndrome. METHODS: Eight families with DMJD were studied by whole-exome or targeted sequencing, with detailed clinical and radiological characterization. Patient-derived induced pluripotent stem cells were derived into neural precursor and endothelial cells to study gene expression. RESULTS: All patients showed biallelic mutations in the nonclustered protocadherin-12 (PCDH12) gene. The characteristic clinical presentation included progressive microcephaly, craniofacial dysmorphism, psychomotor disability, epilepsy, and axial hypotonia with variable appendicular spasticity. Brain imaging showed brainstem malformations and with frequent thinned corpus callosum with punctate brain calcifications, reflecting expression of PCDH12 in neural and endothelial cells. These cells showed lack of PCDH12 expression and impaired neurite outgrowth. INTERPRETATION: DMJD patients have biallelic mutations in PCDH12 and lack of protein expression. These patients present with characteristic microcephaly and abnormalities of white matter tracts. Such pathogenic variants predict a poor outcome as a result of brainstem malformation and evidence of white matter tract defects, and should be added to the phenotypic spectrum associated with PCDH12-related conditions. Ann Neurol 2018;84:646-655.

Biallelic loss of human CTNNA2, encoding αN-catenin, leads to ARP2/3 complex overactivity and disordered cortical neuronal migration.

Neuronal migration defects, including pachygyria, are among the most severe developmental brain defects in humans. Here, we identify biallelic truncating mutations in CTNNA2, encoding αN-catenin, in patients with a distinct recessive form of pachygyria. CTNNA2 was expressed in human cerebral cortex, and its loss in neurons led to defects in neurite stability and migration. The αN-catenin paralog, αE-catenin, acts as a switch regulating the balance between ß-catenin and Arp2/3 actin filament activities1. Loss of αN-catenin did not affect ß-catenin signaling, but recombinant αN-catenin interacted with purified actin and repressed ARP2/3 actin-branching activity. The actin-binding domain of αN-catenin or ARP2/3 inhibitors rescued the neuronal phenotype associated with CTNNA2 loss, suggesting ARP2/3 de-repression as a potential disease mechanism. Our findings identify CTNNA2 as the first catenin family member with biallelic mutations in humans, causing a new pachygyria syndrome linked to actin regulation, and uncover a key factor involved in ARP2/3 repression in neurons.

Correction: Orion: Detecting regions of the human non-coding genome that are intolerant to variation using population genetics.

[This corrects the article DOI: 10.1371/journal.pone.0181604.].

A case-control collapsing analysis identifies epilepsy genes implicated in trio sequencing studies focused on de novo mutations.

Trio exome sequencing has been successful in identifying genes with de novo mutations (DNMs) causing epileptic encephalopathy (EE) and other neurodevelopmental disorders. Here, we evaluate how well a case-control collapsing analysis recovers genes causing dominant forms of EE originally implicated by DNM analysis. We performed a genome-wide search for an enrichment of "qualifying variants" in protein-coding genes in 488 unrelated cases compared to 12,151 unrelated controls. These "qualifying variants" were selected to be extremely rare variants predicted to functionally impact the protein to enrich for likely pathogenic variants. Despite modest sample size, three known EE genes (KCNT1, SCN2A, and STXBP1) achieved genome-wide significance (p<2.68×10-6). In addition, six of the 10 most significantly associated genes are known EE genes, and the majority of the known EE genes (17 out of 25) originally implicated in trio sequencing are nominally significant (p<0.05), a proportion significantly higher than the expected (Fisher's exact p = 2.33×10-17). Our results indicate that a case-control collapsing analysis can identify several of the EE genes originally implicated in trio sequencing studies, and clearly show that additional genes would be implicated with larger sample sizes. The case-control analysis not only makes discovery easier and more economical in early onset disorders, particularly when large cohorts are available, but also supports the use of this approach to identify genes in diseases that present later in life when parents are not readily available.

Orion: Detecting regions of the human non-coding genome that are intolerant to variation using population genetics.

There is broad agreement that genetic mutations occurring outside of the protein-coding regions play a key role in human disease. Despite this consensus, we are not yet capable of discerning which portions of non-coding sequence are important in the context of human disease. Here, we present Orion, an approach that detects regions of the non-coding genome that are depleted of variation, suggesting that the regions are intolerant of mutations and subject to purifying selection in the human lineage. We show that Orion is highly correlated with known intolerant regions as well as regions that harbor putatively pathogenic variation. This approach provides a mechanism to identify pathogenic variation in the human non-coding genome and will have immediate utility in the diagnostic interpretation of patient genomes and in large case control studies using whole-genome sequences.

An AKT3-FOXG1-reelin network underlies defective migration in human focal malformations of cortical development.

Focal malformations of cortical development (FMCDs) account for the majority of drug-resistant pediatric epilepsy. Postzygotic somatic mutations activating the phosphatidylinositol-4,5-bisphosphate-3-kinase (PI3K)-protein kinase B (AKT)-mammalian target of rapamycin (mTOR) pathway are found in a wide range of brain diseases, including FMCDs. It remains unclear how a mutation in a small fraction of cells disrupts the architecture of the entire hemisphere. Within human FMCD-affected brain, we found that cells showing activation of the PI3K-AKT-mTOR pathway were enriched for the AKT3(E17K) mutation. Introducing the FMCD-causing mutation into mouse brain resulted in electrographic seizures and impaired hemispheric architecture. Mutation-expressing neural progenitors showed misexpression of reelin, which led to a non-cell autonomous migration defect in neighboring cells, due at least in part to derepression of reelin transcription in a manner dependent on the forkhead box (FOX) transcription factor FOXG1. Treatments aimed at either blocking downstream AKT signaling or inactivating reelin restored migration. These findings suggest a central AKT-FOXG1-reelin signaling pathway in FMCD and support pathway inhibitors as potential treatments or therapies for some forms of focal epilepsy.

Inactivating mutations in MFSD2A, required for omega-3 fatty acid transport in brain, cause a lethal microcephaly syndrome.

Docosahexanoic acid (DHA) is the most abundant omega-3 fatty acid in brain, and, although it is considered essential, deficiency has not been linked to disease. Despite the large mass of DHA in phospholipids, the brain does not synthesize it. DHA is imported across the blood-brain barrier (BBB) through the major facilitator superfamily domain-containing 2a (MFSD2A) protein. MFSD2A transports DHA as well as other fatty acids in the form of lysophosphatidylcholine (LPC). We identify two families displaying MFSD2A mutations in conserved residues. Affected individuals exhibited a lethal microcephaly syndrome linked to inadequate uptake of LPC lipids. The MFSD2A mutations impaired transport activity in a cell-based assay. Moreover, when expressed in mfsd2aa-morphant zebrafish, mutants failed to rescue microcephaly, BBB breakdown and lethality. Our results establish a link between transport of DHA and LPCs by MFSD2A and human brain growth and function, presenting the first evidence of monogenic disease related to transport of DHA in humans.

Functional genome-wide siRNA screen identifies KIAA0586 as mutated in Joubert syndrome.

Defective primary ciliogenesis or cilium stability forms the basis of human ciliopathies, including Joubert syndrome (JS), with defective cerebellar vermis development. We performed a high-content genome-wide small interfering RNA (siRNA) screen to identify genes regulating ciliogenesis as candidates for JS. We analyzed results with a supervised-learning approach, using SYSCILIA gold standard, Cildb3.0, a centriole siRNA screen and the GTex project, identifying 591 likely candidates. Intersection of this data with whole exome results from 145 individuals with unexplained JS identified six families with predominantly compound heterozygous mutations in KIAA0586. A c.428del base deletion in 0.1% of the general population was found in trans with a second mutation in an additional set of 9 of 163 unexplained JS patients. KIAA0586 is an orthologue of chick Talpid3, required for ciliogenesis and Sonic hedgehog signaling. Our results uncover a relatively high frequency cause for JS and contribute a list of candidates for future gene discoveries in ciliopathies.

Biallelic mutations in SNX14 cause a syndromic form of cerebellar atrophy and lysosome-autophagosome dysfunction.

Pediatric-onset ataxias often present clinically as developmental delay and intellectual disability, with prominent cerebellar atrophy as a key neuroradiographic finding. Here we describe a new clinically distinguishable recessive syndrome in 12 families with cerebellar atrophy together with ataxia, coarsened facial features and intellectual disability, due to truncating mutations in the sorting nexin gene SNX14, encoding a ubiquitously expressed modular PX domain-containing sorting factor. We found SNX14 localized to lysosomes and associated with phosphatidylinositol (3,5)-bisphosphate, a key component of late endosomes/lysosomes. Patient-derived cells showed engorged lysosomes and a slower autophagosome clearance rate upon autophagy induction by starvation. Zebrafish morphants for snx14 showed dramatic loss of cerebellar parenchyma, accumulation of autophagosomes and activation of apoptosis. Our results characterize a unique ataxia syndrome due to biallelic SNX14 mutations leading to lysosome-autophagosome dysfunction.

Biallelic truncating mutations in FMN2, encoding the actin-regulatory protein Formin 2, cause nonsyndromic autosomal-recessive intellectual disability.

Dendritic spines represent the major site of neuronal activity in the brain; they serve as the receiving point for neurotransmitters and undergo rapid activity-dependent morphological changes that correlate with learning and memory. Using a combination of homozygosity mapping and next-generation sequencing in two consanguineous families affected by nonsyndromic autosomal-recessive intellectual disability, we identified truncating mutations in formin 2 (FMN2), encoding a protein that belongs to the formin family of actin cytoskeleton nucleation factors and is highly expressed in the maturing brain. We found that FMN2 localizes to punctae along dendrites and that germline inactivation of mouse Fmn2 resulted in animals with decreased spine density; such mice were previously demonstrated to have a conditioned fear-learning defect. Furthermore, patient neural cells derived from induced pluripotent stem cells showed correlated decreased synaptic density. Thus, FMN2 mutations link intellectual disability either directly or indirectly to the regulation of actin-mediated synaptic spine density.