Rare complete knockouts in humans: population distribution and significant role in autism spectrum disorders

Single nucleotide variations encoding missense mutations in G protein-coupled receptors may contribute to autism

Autism is a neurodevelopmental condition with a range of symptoms that vary in intensity and severity from person to person. Genetic sequencing has identified thousands of genes containing mutations in autistic individuals, which may contribute to the development of autistic symptoms. Several of these genes encode G protein-coupled receptors (GPCRs), which are cell surface expressed proteins that transduce extracellular messages to the intracellular space. Mutations in GPCRs can impact their function, resulting in aberrant signalling within cells and across neurotransmitter systems in the brain. This review summarises the current knowledge on autism-associated single nucleotide variations encoding missense mutations in GPCRs and the impact of these genetic mutations on GPCR function. For some autism-associated mutations, changes in GPCR expression levels, ligand affinity, potency and efficacy have been observed. However, for many the functional consequences remain unknown. Thus, further work to characterise the functional impacts of the genetically identified mutations is required. LINKED ARTICLES: This article is part of a themed issue Therapeutic Targeting of G Protein-Coupled Receptors: hot topics from the Australasian Society of Clinical and Experimental Pharmacologists and Toxicologists 2021 Virtual Annual Scientific Meeting. To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v181.14/issuetoc.

Expanding the Clinical Spectrum of DRP2-Associated Charcot-Marie-Tooth Disease

BACKGROUND AND OBJECTIVES

Germline truncating variants in the DRP2 gene (encoding dystrophin-related protein 2) cause the disruption of the periaxin-DRP2-dystroglycan complex and have been linked to Charcot-Marie-Tooth disease. However, the causality and the underlying phenotype of the genetic alterations are not clearly defined.

METHODS

This cross-sectional retrospective observational study includes 9 patients with Charcot-Marie-Tooth disease (CMT) with DRP2 germline variants evaluated at 6 centers throughout Spain.

RESULTS

We identified 7 Spanish families with 4 different DRP2 likely pathogenic germline variants. In agreement with an X-linked inheritance, men harboring hemizygous DRP2 variants presented with an intermediate form of CMT, whereas heterozygous women were asymptomatic. Symptom onset was variable (36.6 ± 16 years), with lower limb weakness and multimodal sensory loss producing a mild-to-moderate functional impairment. Nerve echography revealed an increase in the cross-sectional area of nerve roots and proximal nerves. Lower limb muscle magnetic resonance imaging confirmed the presence of a length-dependent fatty infiltration. Immunostaining in intradermal nerve fibers demonstrated the absence of DRP2 and electron microscopy revealed abnormal myelin thickness that was also detectable in the sural nerve sections.

DISCUSSION

Our findings support the causality of DRP2 pathogenic germline variants in CMT and further define the phenotype as a late-onset sensory and motor length-dependent neuropathy, with intermediate velocities and thickening of proximal nerve segments.

Loss of interleukin 1 signaling causes impairment of microglia- mediated synapse elimination and autistic-like behaviour in mice

In the last years, the hypothesis that elevated levels of proinflammatory cytokines contribute to the pathogenesis of neurodevelopmental diseases has gained popularity. IL-1 is one of the main cytokines found to be elevated in Autism spectrum disorder (ASD), a complex neurodevelopmental condition characterized by defects in social communication and cognitive impairments. In this study, we demonstrate that mice lacking IL-1 signaling display autistic-like defects associated with an excessive number of synapses. We also show that microglia lacking IL-1 signaling at early neurodevelopmental stages are unable to properly perform the process of synapse engulfment and display excessive activation of mammalian target of rapamycin (mTOR) signaling. Notably, even the acute inhibition of IL-1R1 by IL-1Ra is sufficient to enhance mTOR signaling and reduce synaptosome phagocytosis in WT microglia. Finally, we demonstrate that rapamycin treatment rescues the defects in IL-1R deficient mice. These data unveil an exclusive role of microglial IL-1 in synapse refinement via mTOR signaling and indicate a novel mechanism possibly involved in neurodevelopmental disorders associated with defects in the IL-1 pathway.

Massively parallel disruption of enhancers active in human neural stem cells

Changes in gene regulation have been linked to the expansion of the human cerebral cortex and to neurodevelopmental disorders, potentially by altering neural progenitor proliferation. However, the effects of genetic variation within regulatory elements on neural progenitors remain obscure. We use sgRNA-Cas9 screens in human neural stem cells (hNSCs) to disrupt 10,674 genes and 26,385 conserved regions in 2,227 enhancers active in the developing human cortex and determine effects on proliferation. Genes with proliferation phenotypes are associated with neurodevelopmental disorders and show biased expression in specific fetal human brain neural progenitor populations. Although enhancer disruptions overall have weaker effects than gene disruptions, we identify enhancer disruptions that severely alter hNSC self-renewal. Disruptions in human accelerated regions, implicated in human brain evolution, also alter proliferation. Integrating proliferation phenotypes with chromatin interactions reveals regulatory relationships between enhancers and their target genes contributing to neurogenesis and potentially to human cortical evolution.

Evolution of the nervous system by acquisition of retrovirus-derived genes in mammals

In the course of evolution, the most highly developed organ is likely the brain, which has become more complex over time and acquired diverse forms and functions in different species. In particular, mammals have developed complex and high-functioning brains, and it has been reported that several genes derived from retroviruses were involved in mammalian brain evolution, that is, generating the complexity of the nervous system. Especially, the sushi-ichi-related retrotransposon homolog (SIRH)/retrotransposon gag-like (RTL) genes have been suggested to play a role in the evolutionary processes shaping brain morphology and function in mammals. Genetic mutation and altered expression of genes are linked to neurological disorders, highlighting how the acquisition of virus-derived genes in mammals has both driven brain evolution and imposed a susceptibility to diseases. This review provides an overview of the functions, diversity, evolution and diseases associated with SIRH/RTL genes in the nervous system. The contribution of retroviruses to brain evolution is an important research topic in evolutionary biology and neuroscience, and further insights are expected to be gained through future studies.

Intellectual Disability and Behavioral Deficits Linked to CYFIP1 Missense Variants Disrupting Actin Polymerization

BACKGROUND

15q11.2 deletions and duplications have been linked to autism spectrum disorder, schizophrenia, and intellectual disability. Recent evidence suggests that dysfunctional CYFIP1 (cytoplasmic FMR1 interacting protein 1) contributes to the clinical phenotypes observed in individuals with 15q11.2 deletion/duplication syndrome. CYFIP1 plays crucial roles in neuronal development and brain connectivity, promoting actin polymerization and regulating local protein synthesis. However, information about the impact of single nucleotide variants in CYFIP1 on neurodevelopmental disorders is limited.

METHODS

Here, we report a family with 2 probands exhibiting intellectual disability, autism spectrum disorder, spastic tetraparesis, and brain morphology defects and who carry biallelic missense point mutations in the CYFIP1 gene. We used skin fibroblasts from one of the probands, the parents, and typically developing individuals to investigate the effect of the variants on the functionality of CYFIP1. In addition, we generated Drosophila knockin mutants to address the effect of the variants in vivo and gain insight into the molecular mechanism that underlies the clinical phenotype.

RESULTS

Our study revealed that the 2 missense variants are in protein domains responsible for maintaining the interaction within the wave regulatory complex. Molecular and cellular analyses in skin fibroblasts from one proband showed deficits in actin polymerization. The fly model for these mutations exhibited abnormal brain morphology and F-actin loss and recapitulated the core behavioral symptoms, such as deficits in social interaction and motor coordination.

CONCLUSIONS

Our findings suggest that the 2 CYFIP1 variants contribute to the clinical phenotype in the probands that reflects deficits in actin-mediated brain development processes.

Rare X-linked variants carry predominantly male risk in autism, Tourette syndrome, and ADHD

Autism spectrum disorder (ASD), Tourette syndrome (TS), and attention-deficit/hyperactivity disorder (ADHD) display strong male sex bias, due to a combination of genetic and biological factors, as well as selective ascertainment. While the hemizygous nature of chromosome X (Chr X) in males has long been postulated as a key point of "male vulnerability", rare genetic variation on this chromosome has not been systematically characterized in large-scale whole exome sequencing studies of "idiopathic" ASD, TS, and ADHD. Here, we take advantage of informative recombinations in simplex ASD families to pinpoint risk-enriched regions on Chr X, within which rare maternally-inherited damaging variants carry substantial risk in males with ASD. We then apply a modified transmission disequilibrium test to 13,052 ASD probands and identify a novel high confidence ASD risk gene at exome-wide significance (MAGEC3). Finally, we observe that rare damaging variants within these risk regions carry similar effect sizes in males with TS or ADHD, further clarifying genetic mechanisms underlying male vulnerability in multiple neurodevelopmental disorders that can be exploited for systematic gene discovery.

Diversity and Classification of Genetic Variations in Autism Spectrum Disorder

Autism spectrum disorder (ASD) is a neurodevelopmental condition with symptoms that affect the whole personality and all aspects of life. Although there is a high degree of heterogeneity in both its etiology and its characteristic behavioral patterns, the disorder is well-captured along the autistic triad. Currently, ASD status can be confirmed following an assessment of behavioral features, but there is a growing emphasis on conceptualizing autism as a spectrum, which allows for establishing a diagnosis based on the level of support need, free of discrete categories. Since ASD has a high genetic predominance, the number of genetic variations identified in the background of the condition is increasing exponentially as genetic testing methods are rapidly evolving. However, due to the huge amount of data to be analyzed, grouping the different DNA variations is still challenging. Therefore, in the present review, a multidimensional classification scheme was developed to accommodate most of the currently known genetic variants associated with autism. Genetic variations have been grouped according to six criteria (extent, time of onset, information content, frequency, number of genes involved, inheritance pattern), which are themselves not discrete categories, but form a coherent continuum in line with the autism spectrum approach.

The autism spectrum disorder risk gene NEXMIF over-synchronizes hippocampal CA1 network and alters neuronal coding

Mutations in autism spectrum disorder (ASD) risk genes disrupt neural network dynamics that ultimately lead to abnormal behavior. To understand how ASD-risk genes influence neural circuit computation during behavior, we analyzed the hippocampal network by performing large-scale cellular calcium imaging from hundreds of individual CA1 neurons simultaneously in transgenic mice with total knockout of the X-linked ASD-risk gene NEXMIF (neurite extension and migration factor). As NEXMIF knockout in mice led to profound learning and memory deficits, we examined the CA1 network during voluntary locomotion, a fundamental component of spatial memory. We found that NEXMIF knockout does not alter the overall excitability of individual neurons but exaggerates movement-related neuronal responses. To quantify network functional connectivity changes, we applied closeness centrality analysis from graph theory to our large-scale calcium imaging datasets, in addition to using the conventional pairwise correlation analysis. Closeness centrality analysis considers both the number of connections and the connection strength between neurons within a network. We found that in wild-type mice the CA1 network desynchronizes during locomotion, consistent with increased network information coding during active behavior. Upon NEXMIF knockout, CA1 network is over-synchronized regardless of behavioral state and fails to desynchronize during locomotion, highlighting how perturbations in ASD-implicated genes create abnormal network synchronization that could contribute to ASD-related behaviors.

Ankyrin B promotes developmental spine regulation in the mouse prefrontal cortex

Postnatal regulation of dendritic spine formation and refinement in cortical pyramidal neurons is critical for excitatory/inhibitory balance in neocortical networks. Recent studies have identified a selective spine pruning mechanism in the mouse prefrontal cortex mediated by class 3 Semaphorins and the L1 cell adhesion molecules, neuron-glia related cell adhesion molecule, Close Homolog of L1, and L1. L1 cell adhesion molecules bind Ankyrin B, an actin-spectrin adaptor encoded by Ankyrin2, a high-confidence gene for autism spectrum disorder. In a new inducible mouse model (Nex1Cre-ERT2: Ank2flox: RCE), Ankyrin2 deletion in early postnatal pyramidal neurons increased spine density on apical dendrites in prefrontal cortex layer 2/3 of homozygous and heterozygous Ankyrin2-deficient mice. In contrast, Ankyrin2 deletion in adulthood had no effect on spine density. Sema3F-induced spine pruning was impaired in cortical neuron cultures from Ankyrin B-null mice and was rescued by re-expression of the 220 kDa Ankyrin B isoform but not 440 kDa Ankyrin B. Ankyrin B bound to neuron-glia related CAM at a cytoplasmic domain motif (FIGQY1231), and mutation to FIGQH inhibited binding, impairing Sema3F-induced spine pruning in neuronal cultures. Identification of a novel function for Ankyrin B in dendritic spine regulation provides insight into cortical circuit development, as well as potential molecular deficiencies in autism spectrum disorder.

Genomic architecture of autism spectrum disorder in Qatar: The BARAKA-Qatar Study

BACKGROUND

Autism spectrum disorder (ASD) is a neurodevelopmental condition characterized by impaired social and communication skills, restricted interests, and repetitive behaviors. The prevalence of ASD among children in Qatar was recently estimated to be 1.1%, though the genetic architecture underlying ASD both in Qatar and the greater Middle East has been largely unexplored. Here, we describe the first genomic data release from the BARAKA-Qatar Study-a nationwide program building a broadly consented biorepository of individuals with ASD and their families available for sample and data sharing and multi-omics research.

METHODS

In this first release, we present a comprehensive analysis of whole-genome sequencing (WGS) data of the first 100 families (372 individuals), investigating the genetic architecture, including single-nucleotide variants (SNVs), copy number variants (CNVs), tandem repeat expansions (TREs), as well as mitochondrial DNA variants (mtDNA) segregating with ASD in local families.

RESULTS

Overall, we identify potentially pathogenic variants in known genes or regions in 27 out of 100 families (27%), of which 11 variants (40.7%) were classified as pathogenic or likely-pathogenic based on American College of Medical Genetics (ACMG) guidelines. Dominant variants, including de novo and inherited, contributed to 15 (55.6%) of these families, consisting of SNVs/indels (66.7%), CNVs (13.3%), TREs (13.3%), and mtDNA variants (6.7%). Moreover, homozygous variants were found in 7 families (25.9%), with a sixfold increase in homozygous burden in consanguineous versus non-consanguineous families (13.6% and 1.8%, respectively). Furthermore, 28 novel ASD candidate genes were identified in 20 families, 23 of which had recurrent hits in MSSNG and SSC cohorts.

CONCLUSIONS

This study illustrates the value of ASD studies in under-represented populations and the importance of WGS as a comprehensive tool for establishing a molecular diagnosis for families with ASD. Moreover, it uncovers a significant role for recessive variation in ASD architecture in consanguineous settings and provides a unique resource of Middle Eastern genomes for future research to the global ASD community.

Retrovirus-Derived RTL9 Plays an Important Role in Innate Antifungal Immunity in the Eutherian Brain

Retrotransposon Gag-like (RTL) genes play a variety of essential and important roles in the eutherian placenta and brain. It has recently been demonstrated that RTL5 and RTL6 (also known as sushi-ichi retrotransposon homolog 8 (SIRH8) and SIRH3) are microglial genes that play important roles in the brain's innate immunity against viruses and bacteria through their removal of double-stranded RNA and lipopolysaccharide, respectively. In this work, we addressed the function of RTL9 (also known as SIRH10). Using knock-in mice that produce RTL9-mCherry fusion protein, we examined RTL9 expression in the brain and its reaction to fungal zymosan. Here, we demonstrate that RTL9 plays an important role, degrading zymosan in the brain. The RTL9 protein is localized in the microglial lysosomes where incorporated zymosan is digested. Furthermore, in Rtl9 knockout mice expressing RTL9ΔC protein lacking the C-terminus retroviral GAG-like region, the zymosan degrading activity was lost. Thus, RTL9 is essentially engaged in this reaction, presumably via its GAG-like region. Together with our previous study, this result highlights the importance of three retrovirus-derived microglial RTL genes as eutherian-specific constituents of the current brain innate immune system: RTL9, RTL5 and RTL6, responding to fungi, viruses and bacteria, respectively.

Roles of retrovirus-derived PEG10 and PEG11/RTL1 in mammalian development and evolution and their involvement in human disease

PEG10 and PEG11/RTL1 are paternally expressed, imprinted genes that play essential roles in the current eutherian developmental system and are therefore associated with developmental abnormalities caused by aberrant genomic imprinting. They are also presumed to be retrovirus-derived genes with homology to the sushi-ichi retrotransposon GAG and POL, further expanding our comprehension of mammalian evolution via the domestication (exaptation) of retrovirus-derived acquired genes. In this manuscript, we review the importance of PEG10 and PEG11/RTL1 in genomic imprinting research via their functional roles in development and human disease, including neurodevelopmental disorders of genomic imprinting, Angelman, Kagami-Ogata and Temple syndromes, and the impact of newly inserted DNA on the emergence of newly imprinted regions. We also discuss their possible roles as ancestors of other retrovirus-derived RTL/SIRH genes that likewise play important roles in the current mammalian developmental system, such as in the placenta, brain and innate immune system.

Retrovirus-Derived RTL/SIRH: Their Diverse Roles in the Current Eutherian Developmental System and Contribution to Eutherian Evolution

Eutherians have 11 retrotransposon Gag-like (RTL)/sushi-ichi retrotransposon homolog (SIRH) genes presumably derived from a certain retrovirus. Accumulating evidence indicates that the RTL/SIRH genes play a variety of roles in the current mammalian developmental system, such as in the placenta, brain, and innate immune system, in a eutherian-specific manner. It has been shown that the functional role of Paternally Expressed 10 (PEG10) in placental formation is unique to the therian mammals, as are the eutherian-specific roles of PEG10 and PEG11/RTL1 in maintaining the fetal capillary network and the endocrine regulation of RTL7/SIRH7 (aka Leucine Zipper Down-Regulated in Cancer 1 (LDOCK1)) in the placenta. In the brain, PEG11/RTL1 is expressed in the corticospinal tract and hippocampal commissure, mammalian-specific structures, and in the corpus callosum, a eutherian-specific structure. Unexpectedly, at least three RTL/SIRH genes, RTL5/SIRH8, RTL6/SIRH3, and RTL9/SIRH10, play important roles in combating a variety of pathogens, namely viruses, bacteria, and fungi, respectively, suggesting that the innate immunity system of the brain in eutherians has been enhanced by the emergence of these new components. In this review, we will summarize the function of 10 out of the 11 RTL/SIRH genes and discuss their roles in eutherian development and evolution.

TMCO3, a Putative K+ :Proton Antiporter at the Golgi Apparatus, Is Important for Longitudinal Growth in Mice and Humans

Isolated short stature, defined as short stature without any other abnormalities, is a common heterogeneous condition in children. Exome sequencing identified the homozygous nonsense variant c.1832G>A/p.(Trp611*) in TMCO3 in two sisters with isolated short stature. Radiological studies, biochemical measurements, assessment of the skeletal status, and three-dimensional bone microarchitecture revealed no relevant skeletal and bone abnormalities in both sisters. The homozygous TMCO3 variant segregated with short stature in the family. TMCO3 transcript levels were reduced by ~50% in leukocyte-derived RNA of both sisters compared with controls, likely due to nonsense-mediated mRNA decay. In primary urinary cells of heterozygous family members, we detected significantly reduced TMCO3 protein levels. TMCO3 is functionally uncharacterized. We ectopically expressed wild-type TMCO3 in HeLa and ATDC5 chondrogenic cells and detected TMCO3 predominantly at the Golgi apparatus, whereas the TMCO3W611* mutant did not reach the Golgi. Coordinated co-expression of TMCO3W611* -HA and EGFP in HeLa cells confirmed intrinsic instability and/or degradation of the mutant. Tmco3 is expressed in all relevant mouse skeletal cell types. Highest abundance of Tmco3 was found in chondrocytes of the prehypertrophic zone in mouse and minipig growth plates where it co-localizes with a Golgi marker. Knockdown of Tmco3 in differentiated ATDC5 cells caused reduced and increased expression of Pthlh and Ihh, respectively. Measurement of long bones in Tmco3tm1b(KOMP)Wtsi knockout mice revealed significant shortening of forelimbs and hindlimbs. TMCO3 is a potential member of the monovalent cation:proton antiporter 2 (CPA2) family. By in silico tools and homology modeling, TMCO3 is predicted to have an N-terminal secretory signal peptide, forms a dimer localized to the membrane, and is organized in a dimerization and a core domain. The core domain contains the CPA2 motif essential for K+ binding and selectivity. Collectively, our data demonstrate that loss of TMCO3 causes growth defects in both humans and mice. © 2023 American Society for Bone and Mineral Research (ASBMR).

Ankyrin B Promotes Developmental Spine Regulation in the Mouse Prefrontal Cortex

Postnatal regulation of dendritic spine formation and refinement in cortical pyramidal neurons is critical for excitatory/inhibitory balance in neocortical networks. Recent studies have identified a selective spine pruning mechanism in the mouse prefrontal cortex (PFC) mediated by class 3 Semaphorins and the L1-CAM cell adhesion molecules Neuron-glia related CAM (NrCAM), Close Homolog of L1 (CHL1), and L1. L1-CAMs bind Ankyrin B (AnkB), an actin-spectrin adaptor encoded by Ankyrin2 ( ANK2 ), a high confidence gene for autism spectrum disorder (ASD). In a new inducible mouse model (Nex1Cre-ERT2: Ank2 flox : RCE), Ank2 deletion in early postnatal pyramidal neurons increased spine density on apical dendrites in PFC layer 2/3 of homozygous and heterozygous Ank2 -deficient mice. In contrast, Ank2 deletion in adulthood had no effect on spine density. Sema3F-induced spine pruning was impaired in cortical neuron cultures from AnkB-null mice and was rescued by re-expression of the 220 kDa AnkB isoform but not 440 kDa AnkB. AnkB bound to NrCAM at a cytoplasmic domain motif (FIGQY 1231 ), and mutation to FIGQH inhibited binding, impairing Sema3F-induced spine pruning in neuronal cultures. Identification of a novel function for AnkB in dendritic spine regulation provides insight into cortical circuit development, as well as potential molecular deficiencies in ASD.

Molecular diagnosis of 405 individuals with autism spectrum disorder

Autism spectrum disorder (ASD) is caused by combined genetic and environmental factors. Genetic heritability in ASD is estimated as 60-90%, and genetic investigations have revealed many monogenic factors. We analyzed 405 patients with ASD using family-based exome sequencing to detect disease-causing single-nucleotide variants (SNVs), small insertions and deletions (indels), and copy number variations (CNVs) for molecular diagnoses. All candidate variants were validated by Sanger sequencing or quantitative polymerase chain reaction and were evaluated using the American College of Medical Genetics and Genomics/Association for Molecular Pathology guidelines for molecular diagnosis. We identified 55 disease-causing SNVs/indels in 53 affected individuals and 13 disease-causing CNVs in 13 affected individuals, achieving a molecular diagnosis in 66 of 405 affected individuals (16.3%). Among the 55 disease-causing SNVs/indels, 51 occurred de novo, 2 were compound heterozygous (in one patient), and 2 were X-linked hemizygous variants inherited from unaffected mothers. The molecular diagnosis rate in females was significantly higher than that in males. We analyzed affected sibling cases of 24 quads and 2 quintets, but only one pair of siblings shared an identical pathogenic variant. Notably, there was a higher molecular diagnostic rate in simplex cases than in multiplex families. Our simulation indicated that the diagnostic yield is increasing by 0.63% (range 0-2.5%) per year. Based on our simple simulation, diagnostic yield is improving over time. Thus, periodical reevaluation of ES data should be strongly encouraged in undiagnosed ASD patients.

Ribonuclease inhibitor 1 (RNH1) deficiency cause congenital cataracts and global developmental delay with infection-induced psychomotor regression and anemia

Ribonuclease inhibitor 1, also known as angiogenin inhibitor 1, encoded by RNH1, is a ubiquitously expressed leucine-rich repeat protein, which is highly conserved in mammalian species. Inactivation of rnh1 in mice causes an embryonically lethal anemia, but the exact biological function of RNH1 in humans remains unknown and no human genetic disease has so far been associated with RNH1. Here, we describe a family with two out of seven siblings affected by a disease characterized by congenital cataract, global developmental delay, myopathy and psychomotor deterioration, seizures and periodic anemia associated with upper respiratory tract infections. A homozygous splice-site variant (c.615-2A > C) in RNH1 segregated with the disease. Sequencing of RNA derived from patient fibroblasts and cDNA analysis of skeletal muscle mRNA showed aberrant splicing with skipping of exon 7. Western blot analysis revealed a total lack of the RNH1 protein. Functional analysis revealed that patient fibroblasts were more sensitive to RNase A exposure, and this phenotype was reversed by transduction with a lentivirus expressing RNH1 to complement patient cells. Our results demonstrate that loss-of-function of RNH1 in humans is associated with a multiorgan developmental disease with recessive inheritance. It may be speculated that the infection-induced deterioration resulted from an increased susceptibility toward extracellular RNases and/or other inflammatory responses normally kept in place by the RNase inhibitor RNH1.

LoFTK: a framework for fully automated calculation of predicted Loss-of-Function variants and genes

BACKGROUND

Loss-of-Function (LoF) variants in human genes are important due to their impact on clinical phenotypes and frequent occurrence in the genomes of healthy individuals. The association of LoF variants with complex diseases and traits may lead to the discovery and validation of novel therapeutic targets. Current approaches predict high-confidence LoF variants without identifying the specific genes or the number of copies they affect. Moreover, there is a lack of methods for detecting knockout genes caused by compound heterozygous (CH) LoF variants.

RESULTS

We have developed the Loss-of-Function ToolKit (LoFTK), which allows efficient and automated prediction of LoF variants from genotyped, imputed and sequenced genomes. LoFTK enables the identification of genes that are inactive in one or two copies and provides summary statistics for downstream analyses. LoFTK can identify CH LoF variants, which result in LoF genes with two copies lost. Using data from parents and offspring we show that 96% of CH LoF genes predicted by LoFTK in the offspring have the respective alleles donated by each parent.

CONCLUSIONS

LoFTK is a command-line based tool that provides a reliable computational workflow for predicting LoF variants from genotyped and sequenced genomes, identifying genes that are inactive in 1 or 2 copies. LoFTK is an open software and is freely available to non-commercial users at https://github.com/CirculatoryHealth/LoFTK .

Fatty Acid Amide Hydrolase: An Integrative Clinical Perspective

Introduction: Fatty acid amide hydrolase (FAAH) is one of the main terminating enzymes of the endocannabinoid system (ECS). Since being discovered in 1996, the modulation of FAAH has been viewed as a compelling alternative strategy to obtain the beneficial effect of the ECS. With a considerable amount of FAAH-related publication over time, the next step would be to comprehend the proximity of this evidence for clinical application. Objective: This review intends to highlight the rationale of FAAH modulation and provide the latest evidence from clinical studies. Methods: Publication searches were conducted to gather information focused on FAAH-related clinical evidence with an extension to the experimental research to understand the biological plausibility. The subtopics were selected to be multidisciplinary to offer more perspective on the current state of the arts. Discussion: Experimental and clinical studies have demonstrated that FAAH was highly expressed not only in the central nervous system but also in the peripheral tissues. As the key regulator of endocannabinoid signaling, it would appear that FAAH plays a role in the modulation of mood and emotional response, reward system, pain perception, energy metabolism and appetite regulation, inflammation, and other biological processes. Genetic variants may be associated with some conditions such as substance/alcohol use disorders, obesity, and eating disorder. The advancement of functional neuroimaging has enabled the evaluation of the neurochemistry of FAAH in brain tissues and this can be incorporated into clinical trials. Intriguingly, the application of FAAH inhibitors in clinical trials seems to provide less striking results in comparison with the animal models, although some potential still can be seen. Conclusion: Modulation of FAAH has an immense potential to be a new therapeutic candidate for several disorders. Further exploration, however, is still needed to ensure who is the best candidate for the treatment strategy.

Genomic strategies to untangle the etiology of autism: A primer

Autism spectrum disorder (ASD) is a neurodevelopmental disorder characterized by deficits in communication, diminished social skills, and restrictive and repetitive behaviors and interests. ASD affects approximately 2.3% of the population and is highly heterogeneous, both phenotypically and genetically. As genomic technologies advance, our understanding of the genetic architecture of ASD is becoming clearer, encompassing spontaneous and inherited alterations throughout the genome, and delineating alterations that are either rare or common in the population. This commentary provides an overview of the genomic strategies and resulting major findings of genetic alterations associated with ASD.

Diverse monogenic subforms of human spermatogenic failure

Non-obstructive azoospermia (NOA) is the most severe form of male infertility and typically incurable. Defining the genetic basis of NOA has proven challenging, and the most advanced classification of NOA subforms is not based on genetics, but simple description of testis histology. In this study, we exome-sequenced over 1000 clinically diagnosed NOA cases and identified a plausible recessive Mendelian cause in 20%. We find further support for 21 genes in a 2-stage burden test with 2072 cases and 11,587 fertile controls. The disrupted genes are primarily on the autosomes, enriched for undescribed human "knockouts", and, for the most part, have yet to be linked to a Mendelian trait. Integration with single-cell RNA sequencing data shows that azoospermia genes can be grouped into molecular subforms with synchronized expression patterns, and analogs of these subforms exist in mice. This analysis framework identifies groups of genes with known roles in spermatogenesis but also reveals unrecognized subforms, such as a set of genes expressed across mitotic divisions of differentiating spermatogonia. Our findings highlight NOA as an understudied Mendelian disorder and provide a conceptual structure for organizing the complex genetics of male infertility, which may provide a rational basis for disease classification.

Can the "female protective effect" liability threshold model explain sex differences in autism spectrum disorder?

Male sex is a strong risk factor for autism spectrum disorder (ASD). The leading theory for a "female protective effect" (FPE) envisions males and females have "differing thresholds" under a "liability threshold model" (DT-LTM). Specifically, this model posits that females require either a greater number or larger magnitude of risk factors (i.e., greater liability) to manifest ASD, which is supported by the finding that a greater proportion of females with ASD have highly penetrant genetic mutations. Herein, we derive testable hypotheses from the DT-LTM for ASD, investigating heritability, familial recurrence, correlation between ASD penetrance and sex ratio, population traits, clinical features, the stability of the sex ratio across diagnostic changes, and highlight other key prerequisites. Our findings reveal that several key predictions of the DT-LTM are not supported by current data, requiring us to establish a different conceptual framework for evaluating alternate models that explain sex differences in ASD.

The people behind the papers - Masahito Irie, Fumitoshi Ishino and Tomoko Kaneko-Ishino

Viral-derived genes have had a huge impact during mammalian evolution, with many of the exapted genes being expressed in the placenta. Now, new research published in Development describes the importance of two genes with retroviral origins in microglia, the innate immune cells of the brain, which are derived from another extra-embryonic tissue, the yolk sac. We caught up with the first author, Masahito Irie, and the corresponding authors, Fumitoshi Ishino, Professor at Tokyo Medical and Dental University, and Tomoko Kaneko-Ishino, Professor at Tokai University, to hear about more about their research.

A phenotypic spectrum of autism is attributable to the combined effects of rare variants, polygenic risk and sex

The genetic etiology of autism spectrum disorder (ASD) is multifactorial, but how combinations of genetic factors determine risk is unclear. In a large family sample, we show that genetic loads of rare and polygenic risk are inversely correlated in cases and greater in females than in males, consistent with a liability threshold that differs by sex. De novo mutations (DNMs), rare inherited variants and polygenic scores were associated with various dimensions of symptom severity in children and parents. Parental age effects on risk for ASD in offspring were attributable to a combination of genetic mechanisms, including DNMs that accumulate in the paternal germline and inherited risk that influences behavior in parents. Genes implicated by rare variants were enriched in excitatory and inhibitory neurons compared with genes implicated by common variants. Our results suggest that a phenotypic spectrum of ASD is attributable to a spectrum of genetic factors that impact different neurodevelopmental processes.

Genomics, convergent neuroscience and progress in understanding autism spectrum disorder

More than a hundred genes have been identified that, when disrupted, impart large risk for autism spectrum disorder (ASD). Current knowledge about the encoded proteins - although incomplete - points to a very wide range of developmentally dynamic and diverse biological processes. Moreover, the core symptoms of ASD involve distinctly human characteristics, presenting challenges to interpreting evolutionarily distant model systems. Indeed, despite a decade of striking progress in gene discovery, an actionable understanding of pathobiology remains elusive. Increasingly, convergent neuroscience approaches have been recognized as an important complement to traditional uses of genetics to illuminate the biology of human disorders. These methods seek to identify intersection among molecular-level, cellular-level and circuit-level functions across multiple risk genes and have highlighted developing excitatory neurons in the human mid-gestational prefrontal cortex as an important pathobiological nexus in ASD. In addition, neurogenesis, chromatin modification and synaptic function have emerged as key potential mediators of genetic vulnerability. The continued expansion of foundational 'omics' data sets, the application of higher-throughput model systems and incorporating developmental trajectories and sex differences into future analyses will refine and extend these results. Ultimately, a systems-level understanding of ASD genetic risk holds promise for clarifying pathobiology and advancing therapeutics.

Novel genes bearing mutations in rare cases of early-onset ataxia with cerebellar hypoplasia

We propose an approach for the identification of mutant genes for rare diseases in single cases of unknown etiology. All genes with rare biologically significant variants sorted from individual exome data are tested further for profiling of their spatial-temporal and cell/tissue specific expression compared to that of their paralogs. We developed a simple bioinformatics tool ("Essential Paralogue by Expression" (EPbE)) for such analysis. Here, we present rare clinical forms of early ataxia with cerebellar hypoplasia. Using whole-exome sequencing and the EPbE tool, we identified two novel mutant genes previously not associated with congenital human diseases. In Family I, the unique missense mutation (p.Lys258Glu) was found in the LRCH2 gene inherited in an X-linked manner. p.Lys258Glu occurs in the evolutionarily invariant site of the leucine-rich repeat domain of LRCH2. In Family II and Family III, the identical genetic variant was found in the CSMD1 gene inherited as an autosomal-recessive trait. The variant leads to amino acid substitution p.Gly2979Ser in a highly conserved region of the complement-interacting domain of CSMD1. The LRCH2 gene for Family I patients (in which congenital cerebellar hypoplasia was associated with demyelinating polyneuropathy) is expressed in Schwann and precursor Schwann cells and predominantly over its paralogous genes in the developing cerebellar cortex. The CSMD1 gene is predominantly expressed over its paralogous genes in the cerebellum, specifically in the period of late childhood. Thus, the comparative spatial-temporal expression of the selected genes corresponds to the neurological manifestations of the disease.

A family study implicates GBE1 in the etiology of autism spectrum disorder

Autism spectrum disorders (ASD) are neurodevelopmental disorders with an estimated heritability of >60%. Family-based genetic studies of ASD have generally focused on multiple small kindreds, searching for de novo variants of major effect. We hypothesized that molecular genetic analysis of large multiplex families would enable the identification of variants of milder effects. We studied a large multigenerational family of European ancestry with multiple family members affected with ASD or the broader autism phenotype (BAP). We identified a rare heterozygous variant in the gene encoding 1,4-ɑ-glucan branching enzyme 1 (GBE1) that was present in seven of seven individuals with ASD, nine of ten individuals with the BAP, and none of four tested unaffected individuals. We genotyped a community-acquired cohort of 389 individuals with ASD and identified three additional probands. Cascade analysis demonstrated that the variant was present in 11 of 13 individuals with familial ASD/BAP and neither of the two tested unaffected individuals in these three families, also of European ancestry. The variant was not enriched in the combined UK10K ASD cohorts of European ancestry but heterozygous GBE1 deletion was overrepresented in large ASD cohorts, collectively suggesting an association between GBE1 and ASD.

Somatic Mosaicism and Autism Spectrum Disorder

Autism spectrum disorder (ASD) is a genetically heterogenous neurodevelopmental disorder. In the early years of next-generation sequencing, de novo germline variants were shown to contribute to ASD risk. These germline mutations are present in all of the cells of an affected individual and can be detected in any tissue, including clinically accessible DNA sources such as blood or saliva. In recent years, studies have also implicated de novo somatic variants in ASD risk. These somatic mutations arise postzygotically and are present in only a subset of the cells of an affected individual. Depending on the developmental time and progenitor cell in which a somatic mutation occurs, it may be detectable in some tissues and not in others. Somatic mutations detectable at relatively low sequencing coverage in clinically accessible tissues are suggested to contribute to 3-5% of simplex ASD diagnoses, and "brain limited" somatic mutations have been identified in postmortem ASD brain tissue. Somatic mutations likely represent the genetic diagnosis in a proportion of otherwise unexplained individuals with ASD, and brain limited somatic mutations can be used as markers to discover risk genes, cell types, brain regions, and cellular pathways important for ASD pathogenesis and to potentially target for therapeutics.

Lethal variants in humans: lessons learned from a large molecular autopsy cohort

BACKGROUND

Molecular autopsy refers to DNA-based identification of the cause of death. Despite recent attempts to broaden its scope, the term remains typically reserved to sudden unexplained death in young adults. In this study, we aim to showcase the utility of molecular autopsy in defining lethal variants in humans.

METHODS

We describe our experience with a cohort of 481 cases in whom the cause of premature death was investigated using DNA from the index or relatives (molecular autopsy by proxy). Molecular autopsy tool was typically exome sequencing although some were investigated using targeted approaches in the earlier stages of the study; these include positional mapping, targeted gene sequencing, chromosomal microarray, and gene panels.

RESULTS

The study includes 449 cases from consanguineous families and 141 lacked family history (simplex). The age range was embryos to 18 years. A likely causal variant (pathogenic/likely pathogenic) was identified in 63.8% (307/481), a much higher yield compared to the general diagnostic yield (43%) from the same population. The predominance of recessive lethal alleles allowed us to implement molecular autopsy by proxy in 55 couples, and the yield was similarly high (63.6%). We also note the occurrence of biallelic lethal forms of typically non-lethal dominant disorders, sometimes representing a novel bona fide biallelic recessive disease trait. Forty-six disease genes with no OMIM phenotype were identified in the course of this study. The presented data support the candidacy of two other previously reported novel disease genes (FAAH2 and MSN). The focus on lethal phenotypes revealed many examples of interesting phenotypic expansion as well as remarkable variability in clinical presentation. Furthermore, important insights into population genetics and variant interpretation are highlighted based on the results.

CONCLUSIONS

Molecular autopsy, broadly defined, proved to be a helpful clinical approach that provides unique insights into lethal variants and the clinical annotation of the human genome.

Autism Spectrum Disorder and Duchenne Muscular Dystrophy: A Clinical Case on the Potential Role of the Dystrophin in Autism Neurobiology

A diagnosis of autism spectrum disorder is reported in up to 19% of dystrophinopathies. However, over the last ten years, only a few papers have been published on this topic. Therefore, further studies are required to analyze this association in depth and ultimately to understand the role of the brain dystrophin isoform in the pathogenesis of ASD and other neurodevelopmental disorders. In this paper, we report a clinical case of a patient affected by ASD and Duchenne muscular dystrophy, who carries a large deletion of the dystrophin gene. Then we present a brief overview of the literature about similar cases and about the potential role of the dystrophin protein in the neurobiology of autism spectrum disorder.

Molecular Dysregulation in Autism Spectrum Disorder

Autism Spectrum Disorder (ASD) comprises a heterogeneous group of neurodevelopmental disorders with a strong heritable genetic component. At present, ASD is diagnosed solely by behavioral criteria. Advances in genomic analysis have contributed to numerous candidate genes for the risk of ASD, where rare mutations and s common variants contribute to its susceptibility. Moreover, studies show rare de novo variants, copy number variation and single nucleotide polymorphisms (SNPs) also impact neurodevelopment signaling. Exploration of rare and common variants involved in common dysregulated pathways can provide new diagnostic and therapeutic strategies for ASD. Contributions of current innovative molecular strategies to understand etiology of ASD will be explored which are focused on whole exome sequencing (WES), whole genome sequencing (WGS), microRNA, long non-coding RNAs and CRISPR/Cas9 models. Some promising areas of pharmacogenomic and endophenotype directed therapies as novel personalized treatment and prevention will be discussed.

Genetic architecture of autism spectrum disorder: Lessons from large-scale genomic studies

Autism spectrum disorder (ASD) is a neurodevelopmental disorder with a strong genetic component. Recently developed genomic technologies, including microarray and next-generation sequencing (NGS), have enabled researchers to genetic analyses aimed at identifying genetic variations associated with ASD and to elucidate the genetic architecture of the disorder. Large-scale microarray, exome sequencing analyses, and robust statistical methods have resulted in successful gene discovery and identification of high-confidence ASD genes from among de novo and inherited variants. Efforts have been made to understand the genetic architecture of ASD using whole-genome sequencing and genome-wide association studies aimed at identifying noncoding mutations and common variants associated with ASD. In addition, the development of systems biology approaches has resulted in the integration of genetic findings with functional genomic datasets, thereby providing a unique insight into the functional convergence of ASD risk genes and their neurobiology. In this review, we summarize the latest findings of ASD genetic studies involving large cohorts and discuss their implications in ASD neurobiology and in clinical practice.

Targeted sequencing and integrative analysis to prioritize candidate genes in neurodevelopmental disorders

Neurodevelopmental disorders (NDDs) are a group of diseases characterized by high heterogeneity and frequently co-occurring symptoms. The mutational spectrum in patients with NDDs is largely incomplete. Here, we sequenced 547 genes from 1102 patients with NDDs and validated 1271 potential functional variants, including 108 de novo variants (DNVs) in 78 autosomal genes and seven inherited hemizygous variants in six X chromosomal genes. Notably, 36 of these 78 genes are the first to be reported in Chinese patients with NDDs. By integrating our genetic data with public data, we prioritized 212 NDD candidate genes with FDR < 0.1, including 17 novel genes. The novel candidate genes interacted or were co-expressed with known candidate genes, forming a functional network involved in known pathways. We highlighted MSL2, which carried two de novo protein-truncating variants (p.L192Vfs*3 and p.S486Ifs*11) and was frequently connected with known candidate genes. This study provides the mutational spectrum of NDDs in China and prioritizes 212 NDD candidate genes for further functional validation and genetic counseling.

Targeted sequencing and integrative analysis of 3,195 Chinese patients with neurodevelopmental disorders prioritized 26 novel candidate genes

Neurodevelopmental disorders (NDDs) are a set of complex disorders characterized by diverse and co-occurring clinical symptoms. The genetic contribution in patients with NDDs remains largely unknown. Here, we sequence 519 NDD-related genes in 3,195 Chinese probands with neurodevelopmental phenotypes and identify 2,522 putative functional mutations consisting of 137 de novo mutations (DNMs) in 86 genes and 2,385 rare inherited mutations (RIMs) with 22 X-linked hemizygotes in 13 genes, 2 homozygous mutations in 2 genes and 23 compound heterozygous mutations in 10 genes. Furthermore, the DNMs of 16,807 probands with NDDs are retrieved from public datasets and combine in an integrated analysis with the mutation data of our Chinese NDD probands by taking 3,582 in-house controls of Chinese origin as background. We prioritize 26 novel candidate genes. Notably, six of these genes - ITSN1, UBR3, CADM1, RYR3, FLNA, and PLXNA3 - preferably contribute to autism spectrum disorders (ASDs), as demonstrated by high co-expression and/or interaction with ASD genes confirmed via rescue experiments in a mouse model. Importantly, these genes are differentially expressed in the ASD cortex in a significant manner and involved in ASD-associated networks. Together, our study expands the genetic spectrum of Chinese NDDs, further facilitating both basic and translational research.

Identification of primary copy number variations reveal enrichment of Calcium, and MAPK pathways sensitizing secondary sites for autism

Background

Autism is a neurodevelopmental condition with genetic heterogeneity. It is characterized by difficulties in reciprocal social interactions with strong repetitive behaviors and stereotyped interests. Copy number variations (CNVs) are genomic structural variations altering the genomic structure either by duplication or deletion. De novo or inherited CNVs are found in 5–10% of autistic subjects with a size range of few kilobases to several megabases. CNVs predispose humans to various diseases by altering gene regulation, generation of chimeric genes, and disruption of the coding region or through position effect. Although, CNVs are not the initiating event in pathogenesis; additional preceding mutations might be essential for disease manifestation. The present study is aimed to identify the primary CNVs responsible for autism susceptibility in healthy cohorts to sensitize secondary-hits. In the current investigation, primary-hit autism gene CNVs are characterized in 1715 healthy cohorts of varying ethnicities across 12 populations using Affymetrix high-resolution array study. Thirty-eight individuals from twelve families residing in Karnataka, India, with the age group of 13–73 years are included for the comparative CNV analysis. The findings are validated against global 179 autism whole-exome sequence datasets derived from Simons Simplex Collection. These datasets are deposited at the Simons Foundation Autism Research Initiative (SFARI) database.

Results

The study revealed that 34.8% of the subjects carried 2% primary-hit CNV burden with 73 singleton-autism genes in different clusters. Of these, three conserved CNV breakpoints were identified with ARHGAP11B, DUSP22, and CHRNA7 as the target genes across 12 populations. Enrichment analysis of the population-specific autism genes revealed two signaling pathways—calcium and mitogen-activated protein kinases (MAPK) in the CNV identified regions. These impaired pathways affected the downstream cascades of neuronal function and physiology, leading to autism behavior. The pathway analysis of enriched genes unravelled complex protein interaction networks, which sensitized secondary sites for autism. Further, the identification of miRNA targets associated with autism gene CNVs added severity to the condition.

Conclusion

These findings contribute to an atlas of primary-hit genes to detect autism susceptibility in healthy cohorts, indicating their impact on secondary sites for manifestation.

Linking Autism Risk Genes to Disruption of Cortical Development

Autism spectrum disorder (ASD) is a prevalent neurodevelopmental disorder characterized by impairments in social communication and social interaction, and the presence of repetitive behaviors and/or restricted interests. In the past few years, large-scale whole-exome sequencing and genome-wide association studies have made enormous progress in our understanding of the genetic risk architecture of ASD. While showing a complex and heterogeneous landscape, these studies have led to the identification of genetic loci associated with ASD risk. The intersection of genetic and transcriptomic analyses have also begun to shed light on functional convergences between risk genes, with the mid-fetal development of the cerebral cortex emerging as a critical nexus for ASD. In this review, we provide a concise summary of the latest genetic discoveries on ASD. We then discuss the studies in postmortem tissues, stem cell models, and rodent models that implicate recently identified ASD risk genes in cortical development.

Less Is More, Natural Loss-of-Function Mutation Is a Strategy for Adaptation

Gene gain and loss are crucial factors that shape the evolutionary success of diverse organisms. In the past two decades, more attention has been paid to the significance of gene gain through gene duplication or de novo genes. However, gene loss through natural loss-of-function (LoF) mutations, which is prevalent in the genomes of diverse organisms, has been largely ignored. With the development of sequencing techniques, many genomes have been sequenced across diverse species and can be used to study the evolutionary patterns of gene loss. In this review, we summarize recent advances in research on various aspects of LoF mutations, including their identification, evolutionary dynamics in natural populations, and functional effects. In particular, we discuss how LoF mutations can provide insights into the minimum gene set (or the essential gene set) of an organism. Furthermore, we emphasize their potential impact on adaptation. At the genome level, although most LoF mutations are neutral or deleterious, at least some of them are under positive selection and may contribute to biodiversity and adaptation. Overall, we highlight the importance of natural LoF mutations as a robust framework for understanding biological questions in general.

The Evolution of the Mammalian ABCA6-like Genes: Analysis of Phylogenetic, Expression, and Population Genetic Data Reveals Complex Evolutionary Histories

ABC membrane transporters are a large and complex superfamily of ATP-binding cassette transporters that are present in all domains of life. Both their essential function and complexity are reflected by their retention across large expanses of organismal diversity and by the extensive expansion of individual members and subfamilies during evolutionary history. This expansion has resulted in the diverse ABCA transporter family that has in turn evolved into multiple subfamilies. Here, we focus on the ABCA6-like subfamily of ABCA transporters with the goal of understanding their evolutionary history including potential functional changes in, or loss of, individual members. Our analysis finds that ABCA6-like genes, consisting of ABCA6, 8, 9, and 10, are absent from representatives of both monotremes and marsupials and thus the duplications that generated these families most likely occurred at the base of the Eutherian or placental mammals. We have found evidence of both positive and relaxed selection among the ABCA6-like genes, suggesting dynamic changes in function and the potential of gene redundancy. Analysis of the ABCA10 genes further suggests that this gene has undergone relaxed selection only within the human lineage. These findings are complemented by human population data, where we observe an excess of deactivating homozygous mutations. We describe the complex evolutionary history of this ABCA transporter subfamily and demonstrate through the combination of evolutionary and population genetic analysis that ABCA10 is undergoing pseudogenization within humans.

Functional relationships between recessive inherited genes and genes with de novo variants in autism spectrum disorder

Background

Both de novo variants and recessive inherited variants were associated with autism spectrum disorder (ASD). This study aimed to use exome data to prioritize recessive inherited genes (RIGs) with biallelically inherited variants in autosomes or X-linked inherited variants in males and investigate the functional relationships between RIGs and genes with de novo variants (DNGs).

Methods

We used a bioinformatics pipeline to analyze whole-exome sequencing data from 1799 ASD quads (containing one proband, one unaffected sibling, and their parents) from the Simons Simplex Collection and prioritize candidate RIGs with rare biallelically inherited variants in autosomes or X-linked inherited variants in males. The relationships between RIGs and DNGs were characterized based on different genetic perspectives, including genetic variants, functional networks, and brain expression patterns.

Results

Among the biallelically or hemizygous constrained genes that were expressed in the brain, ASD probands carried significantly more biallelically inherited protein-truncating variants (PTVs) in autosomes (p = 0.038) and X-linked inherited PTVs in males (p = 0.026) than those in unaffected siblings. We prioritized eight autosomal, and 13 X-linked candidate RIGs, including 11 genes already associated with neurodevelopmental disorders. In total, we detected biallelically inherited variants or X-linked inherited variants of these 21 candidate RIGs in 26 (1.4%) of 1799 probands. We then integrated previously reported known or candidate genes in ASD, ultimately obtaining 70 RIGs and 87 DNGs for analysis. We found that RIGs were less likely to carry multiple recessive inherited variants than DNGs were to carry multiple de novo variants. Additionally, RIGs and DNGs were significantly co-expressed and interacted with each other, forming a network enriched in known functional ASD clusters, although RIGs were less likely to be enriched in these functional clusters compared with DNGs. Furthermore, although RIGs and DNGs presented comparable expression patterns in the human brain, RIGs were less likely to be associated with prenatal brain regions, the middle cortical layers, and excitatory neurons than DNGs.

Limitations

The RIGs analyzed in this study require functional validation, and the results should be replicated in more patients with ASD.

Conclusions

ASD RIGs were functionally associated with DNGs; however, they exhibited higher heterogeneity than DNGs.

Homozygous deletions implicate non-coding epigenetic marks in Autism spectrum disorder

More than 98% of the human genome is made up of non-coding DNA, but techniques to ascertain its contribution to human disease have lagged far behind our understanding of protein coding variations. Autism spectrum disorder (ASD) has been mostly associated with coding variations via de novo single nucleotide variants (SNVs), recessive/homozygous SNVs, or de novo copy number variants (CNVs); however, most ASD cases continue to lack a genetic diagnosis. We analyzed 187 consanguineous ASD families for biallelic CNVs. Recessive deletions were significantly enriched in affected individuals relative to their unaffected siblings (17% versus 4%, p < 0.001). Only a small subset of biallelic deletions were predicted to result in coding exon disruption. In contrast, biallelic deletions in individuals with ASD were enriched for overlap with regulatory regions, with 23/28 CNVs disrupting histone peaks in ENCODE (p < 0.009). Overlap with regulatory regions was further demonstrated by comparisons to the 127-epigenome dataset released by the Roadmap Epigenomics project, with enrichment for enhancers found in primary brain tissue and neuronal progenitor cells. Our results suggest a novel noncoding mechanism of ASD, describe a powerful method to identify important noncoding regions in the human genome, and emphasize the potential significance of gene activation and regulation in cognitive and social function.

Biallelic loss-of-function variants in TBC1D2B cause a neurodevelopmental disorder with seizures and gingival overgrowth

The family of Tre2-Bub2-Cdc16 (TBC)-domain containing GTPase activating proteins (RABGAPs) is not only known as key regulatorof RAB GTPase activity but also has GAP-independent functions. Rab GTPases are implicated in membrane trafficking pathways, such as vesicular trafficking. We report biallelic loss-of-function variants in TBC1D2B, encoding a member of the TBC/RABGAP family with yet unknown function, as the underlying cause of cognitive impairment, seizures, and/or gingival overgrowth in three individuals from unrelated families. TBC1D2B messenger RNA amount was drastically reduced, and the protein was absent in fibroblasts of two patients. In immunofluorescence analysis, ectopically expressed TBC1D2B colocalized with vesicles positive for RAB5, a small GTPase orchestrating early endocytic vesicle trafficking. In two independent TBC1D2B CRISPR/Cas9 knockout HeLa cell lines that serve as cellular model of TBC1D2B deficiency, epidermal growth factor internalization was significantly reduced compared with the parental HeLa cell line suggesting a role of TBC1D2B in early endocytosis. Serum deprivation of TBC1D2B-deficient HeLa cell lines caused a decrease in cell viability and an increase in apoptosis. Our data reveal that loss of TBC1D2B causes a neurodevelopmental disorder with gingival overgrowth, possibly by deficits in vesicle trafficking and/or cell survival.

The role of rare compound heterozygous events in autism spectrum disorder

The identification of genetic variants underlying autism spectrum disorders (ASDs) may contribute to a better understanding of their underlying biology. To examine the possible role of a specific type of compound heterozygosity in ASD, namely, the occurrence of a deletion together with a functional nucleotide variant on the remaining allele, we sequenced 550 genes in 149 individuals with ASD and their deletion-transmitting parents. This approach allowed us to identify additional sequence variants occurring in the remaining allele of the deletion. Our main goal was to compare the rate of sequence variants in remaining alleles of deleted regions between probands and the deletion-transmitting parents. We also examined the predicted functional effect of the identified variants using Combined Annotation-Dependent Depletion (CADD) scores. The single nucleotide variant-deletion co-occurrence was observed in 13.4% of probands, compared with 8.1% of parents. The cumulative burden of sequence variants (n = 68) in pooled proband sequences was higher than the burden in pooled sequences from the deletion-transmitting parents (n = 41, X2 = 6.69, p = 0.0097). After filtering for those variants predicted to be most deleterious, we observed 21 of such variants in probands versus 8 in their deletion-transmitting parents (X2 = 5.82, p = 0.016). Finally, cumulative CADD scores conferred by these variants were significantly higher in probands than in deletion-transmitting parents (burden test, β = 0.13; p = 1.0 × 10-5). Our findings suggest that the compound heterozygosity described in the current study may be one of several mechanisms explaining variable penetrance of CNVs with known pathogenicity for ASD.

Recent Advances in Understanding the Genetic Architecture of Autism

Recent advances in understanding the genetic architecture of autism spectrum disorder have allowed for unprecedented insight into its biological underpinnings. New studies have elucidated the contributions of a variety of forms of genetic variation to autism susceptibility. While the roles of de novo copy number variants and single-nucleotide variants-causing loss-of-function or missense changes-have been increasingly recognized and refined, mosaic single-nucleotide variants have been implicated more recently in some cases. Moreover, inherited variants (including common variants) and, more recently, rare recessive inherited variants have come into greater focus. Finally, noncoding variants-both inherited and de novo-have been implicated in the last few years. This work has revealed a convergence of diverse genetic drivers on common biological pathways and has highlighted the ongoing importance of increasing sample size and experimental innovation. Continuing to synthesize these genetic findings with functional and phenotypic evidence and translating these discoveries to clinical care remain considerable challenges for the field.

5' splice site GC>GT and GT>GC variants differ markedly in terms of their functionality and pathogenicity

In the human genome, most 5' splice sites (~99%) employ the canonical GT dinucleotide whereas a small minority (~1%) use the noncanonical GC dinucleotide. The functionality and pathogenicity of 5' splice site GT>GC (+2T>C) variants have been extensively studied but we know very little about 5' splice site GC>GT (+2C>T) variants. Herein, we have addressed this deficiency by performing a meta-analysis of reported +2C>T "pathogenic" variants together with a functional analysis of engineered +2C>T substitutions using a cell culture-based full-length gene splicing assay. Our results establish proof of concept that +2C>T variants are qualitatively different from +2T>C variants in terms of their functionality and suggest that, in sharp contrast to +2T>C variants, most if not all +2C>T variants have no pathological relevance. Our findings have important implications for interpreting the clinical relevance of +2C>T variants and understanding the evolutionary switching between GT and GC 5' splice sites in mammalian genomes.

Novel genetic susceptibility loci identified by family based whole exome sequencing in Han Chinese schizophrenia patients

Schizophrenia (SCZ) is a highly heritable psychiatric disorder that affects approximately 1% of population around the world. However, early relevant studies did not reach clear conclusions of the genetic mechanisms of SCZ, suggesting that additional susceptibility loci that exert significant influence on SCZ are yet to be revealed. So, in order to identify novel susceptibility genes that account for the genetic risk of SCZ, we performed a systematic family-based study using whole exome sequencing (WES) in 65 Han Chinese families. The analysis of 51 SCZ trios with both unaffected parents identified 22 exonic and 1 splice-site de novo mutations (DNMs) on a total of 23 genes, and showed that 12 genes carried rare protein-altering compound heterozygous mutations in more than one trio. In addition, we identified 26 exonic or splice-site single nucleotide polymorphisms (SNPs) on 18 genes with nominal significance (P < 5 × 10-4) using a transmission disequilibrium test (TDT) in all the families. Moreover, TDT result confirmed a SCZ susceptibility locus on 3p21.1, encompassing the multigenetic region NEK4-ITIH1-ITIH3-ITIH4. Through several different strategies to predict the potential pathogenic genes in silico, we revealed 4 previous discovered susceptibility genes (TSNARE1, PBRM1, STAB1 and OLIG2) and 4 novel susceptibility loci (PSEN1, TLR5, MGAT5B and SSPO) in Han Chinese SCZ patients. In summary, we identified a list of putative candidate genes for SCZ using a family-based WES approach, thus improving our understanding of the pathology of SCZ and providing critical clues to future functional validation.

Genetic evidence of gender difference in autism spectrum disorder supports the female-protective effect

Autism spectrum disorder (ASD) is a complex neurodevelopmental disorder with a male-to-female prevalence of 4:1. However, the genetic mechanisms underlying this gender difference remain unclear. Mutation burden analysis, a TADA model, and co-expression and functional network analyses were performed on de novo mutations (DNMs) and corresponding candidate genes. We found that the prevalence of putative functional DNMs (loss-of-function and predicted deleterious missense mutations) in females was significantly higher than that in males, suggesting that a higher genetic load was required in females to reach the threshold for a diagnosis. We then prioritized 174 candidate genes, including 60 shared genes, 91 male-specific genes, and 23 female-specific genes. All of the three subclasses of candidate genes were significantly more frequently co-expressed in female brains than male brains, suggesting that compensation effects of the deficiency of ASD candidate genes may be more likely in females. Nevertheless, the three subclasses of candidate genes were co-expressed with each other, suggesting a convergent functional network of male and female-specific genes. Our analysis of different aspects of genetic components provides suggestive evidence supporting the female-protective effect in ASD. Moreover, further study is needed to integrate neuronal and hormonal data to elucidate the underlying gender difference in ASD.

Double hits in schizophrenia

The co-occurrence of a copy number variant (CNV) and a functional variant on the other allele may be a relevant genetic mechanism in schizophrenia. We hypothesized that the cumulative burden of such double hits-in particular those composed of a deletion and a coding single-nucleotide variation (SNV)-is increased in patients with schizophrenia. We combined CNV data with coding variants data in 795 patients with schizophrenia and 474 controls. To limit false CNV-detection, only CNVs called by two algorithms were included. CNV-affected genes were subsequently examined for coding SNVs, which we termed "CNV-SNVs." Correcting for total queried sequence, we assessed the CNV-SNV-burden and the combined predicted deleterious effect. We estimated P-values by permutation of the phenotype. We detected 105 CNV-SNVs; 67 in duplicated and 38 in deleted genic sequence. Although the difference in CNV-SNVs rates was not significant, the combined deleteriousness inferred by CNV-SNVs in deleted sequence was almost 4-fold higher in cases compared with controls (nominal P = 0.009). This effect may be driven by a higher number of CNV-SNVs and/or by a higher degree of predicted deleteriousness of CNV-SNVs. No such effect was observed for duplications. We provide early evidence that deletions co-occurring with a functional variant may be relevant, albeit of modest impact, for the genetic etiology of schizophrenia. Large-scale consortium studies are required to validate our findings. Sequence-based analyses would provide the best resolution for detection of CNVs as well as coding variants genome-wide.

Affected Sib-Pair Analyses Identify Signaling Networks Associated With Social Behavioral Deficits in Autism

Autism spectrum disorders (ASDs) are characterized by deficits in three core behavioral domains: reciprocal social interactions, communication, and restricted interests and/or repetitive behaviors. Several hundreds of risk genes for autism have been identified, however, it remains a challenge to associate these genes with specific core behavioral deficits. In multiplex autism families, affected sibs often show significant differences in severity of individual core phenotypes. We hypothesize that a higher mutation burden contributes to a larger difference in the severity of specific core phenotypes between affected sibs. We tested this hypothesis on social behavioral deficits in autism. We sequenced synaptome genes (n = 1,886) in affected male sib-pairs (n = 274) in families from the Autism Genetics Research Exchange (AGRE) and identified rare (MAF ≤ 1%) and predicted functional variants. We selected affected sib-pairs with a large (≥10; n = 92 pairs) or a small (≤4; n = 108 pairs) difference in total cumulative Autism Diagnostic Interview-Revised (ADI-R) social scores (SOCT_CS). We compared burdens of unshared variants present only in sibs with severe social deficits and found a higher burden in SOCT_CS≥10 compared to SOCT_CS ≤ 4 (SOCT_CS≥10: 705.1 ± 16.2; SOCT_CS ≤ 4, 668.3 ± 9.0; p = 0.025). Unshared SOCT_CS≥10 genes only in sibs with severe social deficits are significantly enriched in the SFARI gene set. Network analyses of these genes using InWeb_IM, molecular signatures database (MSigDB), and GeNetMeta identified enrichment for phosphoinositide 3-kinase (PI3K)-AKT-mammalian target of rapamycin (mTOR) (Enrichment Score [eScore] p value = 3.36E−07; n = 8 genes) and Nerve growth factor (NGF) (eScore p value = 8.94E−07; n = 9 genes) networks. These studies support a key role for these signaling networks in social behavioral deficits and present a novel approach to associate risk genes and signaling networks with core behavioral domains in autism.