KCTD13 is a major driver of mirrored neuroanatomical phenotypes of the 16p11.2 copy number variant

Gene expansions contributing to human brain evolution

Genomic drivers of human-specific neurological traits remain largely undiscovered. Duplicated genes expanded uniquely in the human lineage likely contributed to brain evolution, including the increased complexity of synaptic connections between neurons and the dramatic expansion of the neocortex. Discovering duplicate genes is challenging because the similarity of paralogs makes them prone to sequence-assembly errors. To mitigate this issue, we analyzed a complete telomere-to-telomere human genome sequence (T2T-CHM13) and identified 213 duplicated gene families likely containing human-specific paralogs (>98% identity). Positing that genes important in universal human brain features should exist with at least one copy in all modern humans and exhibit expression in the brain, we narrowed in on 362 paralogs with at least one copy across thousands of ancestrally diverse genomes and present in human brain transcriptomes. Of these, 38 paralogs co-express in gene modules enriched for autism-associated genes and potentially contribute to human language and cognition. We narrowed in on 13 duplicate gene families with human-specific paralogs that are fixed among modern humans and show convincing brain expression patterns. Using long-read DNA sequencing revealed hidden variation across 200 modern humans of diverse ancestries, uncovering signatures of selection not previously identified, including possible balancing selection of CD8B . To understand the roles of duplicated genes in brain development, we generated zebrafish CRISPR "knockout" models of nine orthologs and transiently introduced mRNA-encoding paralogs, effectively "humanizing" the larvae. Morphometric, behavioral, and single-cell RNA-seq screening highlighted, for the first time, a possible role for GPR89B in dosage-mediated brain expansion and FRMPD2B function in altered synaptic signaling, both hallmark features of the human brain. Our holistic approach provides important insights into human brain evolution as well as a resource to the community for studying additional gene expansion drivers of human brain evolution.

Abstract short

Duplicated genes expanded in the human lineage likely contributed to brain evolution, yet challenges exist in their discovery due to sequence-assembly errors. We used a complete telomere-to-telomere genome sequence to identify 213 human-specific gene families. From these, 362 paralogs were found in all modern human genomes tested and brain transcriptomes, making them top candidates contributing to human-universal brain features. Choosing a subset of paralogs, we used long-read DNA sequencing of hundreds of modern humans to reveal previously hidden signatures of selection. To understand their roles in brain development, we generated zebrafish CRISPR "knockout" models of nine orthologs and introduced mRNA-encoding paralogs, effectively "humanizing" larvae. Our findings implicate two new genes in possibly contributing to hallmark features of the human brain: GPR89B in dosage-mediated brain expansion and FRMPD2B in altered synapse signaling. Our holistic approach provides new insights and a comprehensive resource for studying gene expansion drivers of human brain evolution.

The pleiotropic spectrum of proximal 16p11.2 CNVs

Recurrent genomic rearrangements at 16p11.2 BP4-5 represent one of the most common causes of genomic disorders. Originally associated with increased risk for autism spectrum disorder, schizophrenia, and intellectual disability, as well as adiposity and head circumference, these CNVs have since been associated with a plethora of phenotypic alterations, albeit with high variability in expressivity and incomplete penetrance. Here, we comprehensively review the pleiotropy associated with 16p11.2 BP4-5 rearrangements to shine light on its full phenotypic spectrum. Illustrating this phenotypic heterogeneity, we expose many parallels between findings gathered from clinical versus population-based cohorts, which often point to the same physiological systems, and emphasize the role of the CNV beyond neuropsychiatric and anthropometric traits. Revealing the complex and variable clinical manifestations of this CNV is crucial for accurate diagnosis and personalized treatment strategies for carrier individuals. Furthermore, we discuss areas of research that will be key to identifying factors contributing to phenotypic heterogeneity and gaining mechanistic insights into the molecular pathways underlying observed associations, while demonstrating how diversity in affected individuals, cohorts, experimental models, and analytical approaches can catalyze discoveries.

Optimized CRISPR inhibition and activation opens key avenues for systematic biological exploration in zebrafish

The application of CRISPR interference (CRISPRi) and CRISPR activation (CRISPRa) technologies in zebrafish has the potential to expand its capacity for the study of gene function significantly. We have developed a codon optimized CRISPRi/a for zebrafish; here we provide proof-of-principle data across established pigmentary and growth phenotypes. We established a zebrafish codon-optimized cas9 gene, harboring mutations D10A and D839A to render the protein catalytically inactive ( dCas9 ). Similarly, codon-optimized Krüppel associated box (KRAB) and methylated CP2 (MeCP2) inactivating domains or VP64 activator domain were cloned downstream from dCas9 for CRISPRi and CRISPRa, respectively. To validate CRISPRi, we targeted key genes in melanocyte differentiation ( sox10, mitfa, and mitfb) ; and melanin production (tyrosinase; tyr ). Microinjection of CRISPRi mRNA and single guide RNAs (sgRNAs) targeting the tyr promoter or 5'-UTR resulted in larvae with hypopigmented epidermal melanocytes. Transcription factors mitfa and mitfb similarly direct differentiation and pigmentation of epidermal melanocytes ( mitfa ) and retinal pigment epithelium (RPE, mitfb ). CRISPRi-mediated targeting of their promoters or 5'-UTR also results in pronounced hypopigmentation of epidermal melanocytes ( mitfa ), and RPE ( mitfb ). So too, targeting CRISPRi to the sox10 promoter results in hypopigmentation of both epidermal melanocytes and RPE consistent with its role upstream of mitfa and mitfb , and tyr . Finally, we asked whether CRISPRi and CRISPRa could be used to modulate a single gene, to yield hypomorphic and hypermorphic effects, selecting mrap2a , as our target. This gene regulates energy homeostasis and somatic growth via inhibition of the melanocortin 4 receptor gene ( mc4r ). We demonstrate that targeting the mrap2a 5'-UTR with CRISPRa or CRISPRi significantly increases or decreases larval body length, respectively. We demonstrate the utility of CRISPRi/a for modulating control of zebrafish gene expression, facilitating efficient assay of candidate gene function and disease relevance.

Chromosome 16p11.2 microdeletion syndrome with microcephaly and Dandy-Walker malformation spectrum: expanding the known phenotype

BACKGROUND

Chromosome 16p11.2 deletions and duplications were found to be the second most common copy number variation (CNV) reported in cases with clinical presentation suggestive of chromosomal syndromes. Chromosome 16p11.2 deletion syndrome shows remarkable phenotypic heterogeneity with a wide variability of presentation extending from normal development and cognition to severe phenotypes. The clinical spectrum ranges from neurocognitive and global developmental delay (GDD), intellectual disability, and language defects (dysarthria /apraxia) to neuropsychiatric and autism spectrum disorders. Other presentations include dysmorphic features, congenital malformations, insulin resistance, and a tendency for obesity. Our study aims to narrow the gap of knowledge in Saudi Arabia and the Middle Eastern and Northern African (MENA) region about genetic disorders, particularly CNV-associated disorders. Despite their rarity, genetic studies in the MENA region revealed high potential with remarkable genetic and phenotypic novelty.

RESULTS

We identified a heterozygous de novo recurrent proximal chromosome 16p11.2 microdeletion by microarray (arr[GRCh38]16p11.2(29555974_30166595)x1) [(arr[GRCh37]16p11.2(29567295_30177916)x1)] and confirmed by whole exome sequencing (arr[GRCh37]16p11.2(29635211_30199850)x1). We report a Saudi girl with severe motor and cognitive disability, myoclonic epilepsy, deafness, and visual impairment carrying the above-described deletion. Our study broadens the known phenotypic spectrum associated with recurrent proximal 16p11.2 microdeletion syndrome to include developmental dysplasia of the hip, optic atrophy, and a flat retina. Notably, the patient exhibited a rare combination of microcephaly, features consistent with the Dandy-Walker spectrum, and a thin corpus callosum (TCC), which are extremely infrequent presentations in patients with the 16p11.2 microdeletion. Additionally, the patient displayed areas of skin and hair hypopigmentation, attributed to a homozygous hypomorphic allele in the TYR gene.

CONCLUSION

This report expands on the clinical phenotype associated with proximal 16p11.2 microdeletion syndrome, highlighting the potential of genetic research in Saudi Arabia and the MENA region. It underscores the importance of similar future studies.

Characterization of alternative splicing during mammalian brain development reveals the extent of isoform diversity and potential effects on protein structural changes

Regulation of gene expression is critical for fate commitment of stem and progenitor cells during tissue formation. In the context of mammalian brain development, a plethora of studies have described how changes in the expression of individual genes characterize cell types across ontogeny and phylogeny. However, little attention has been paid to the fact that different transcripts can arise from any given gene through alternative splicing (AS). Considered a key mechanism expanding transcriptome diversity during evolution, assessing the full potential of AS on isoform diversity and protein function has been notoriously difficult. Here, we capitalize on the use of a validated reporter mouse line to isolate neural stem cells, neurogenic progenitors and neurons during corticogenesis and combine the use of short- and long-read sequencing to reconstruct the full transcriptome diversity characterizing neurogenic commitment. Extending available transcriptional profiles of the mammalian brain by nearly 50,000 new isoforms, we found that neurogenic commitment is characterized by a progressive increase in exon inclusion resulting in the profound remodeling of the transcriptional profile of specific cortical cell types. Most importantly, we computationally infer the biological significance of AS on protein structure by using AlphaFold2, revealing how radical protein conformational changes can arise from subtle changes in isoforms sequence. Together, our study reveals that AS has a greater potential to impact protein diversity and function than previously thought, independently from changes in gene expression.

A critical review of the impact of candidate copy number variants on autism spectrum disorder

Autism spectrum disorder (ASD) is a complex neurodevelopmental disorder (NDD) influenced by genetic, epigenetic, and environmental factors. Recent advancements in genomic analysis have shed light on numerous genes associated with ASD, highlighting the significant role of both common and rare genetic mutations, as well as copy number variations (CNVs), single nucleotide polymorphisms (SNPs) and unique de novo variants. These genetic variations disrupt neurodevelopmental pathways, contributing to the disorder's complexity. Notably, CNVs are present in 10 %-20 % of individuals with autism, with 3 %-7 % detectable through cytogenetic methods. While the role of submicroscopic CNVs in ASD has been recently studied, their association with genomic loci and genes has not been thoroughly explored. In this review, we focus on 47 CNV regions linked to ASD, encompassing 1632 genes, including protein-coding genes and long non-coding RNAs (lncRNAs), of which 659 show significant brain expression. Using a list of ASD-associated genes from SFARI, we detect 17 regions harboring at least one known ASD-related protein-coding gene. Of the remaining 30 regions, we identify 24 regions containing at least one protein-coding gene with brain-enriched expression and a nervous system phenotype in mouse mutants, and one lncRNA with both brain-enriched expression and upregulation in iPSC to neuron differentiation. This review not only expands our understanding of the genetic diversity associated with ASD but also underscores the potential of lncRNAs in contributing to its etiology. Additionally, the discovered CNVs will be a valuable resource for future diagnostic, therapeutic, and research endeavors aimed at prioritizing genetic variations in ASD.

Understanding copy number variations through their genes: a molecular view on 16p11.2 deletion and duplication syndromes

Neurodevelopmental disorders (NDDs) include a broad spectrum of pathological conditions that affect >4% of children worldwide, share common features and present a variegated genetic origin. They include clinically defined diseases, such as autism spectrum disorders (ASD), attention-deficit/hyperactivity disorder (ADHD), motor disorders such as Tics and Tourette's syndromes, but also much more heterogeneous conditions like intellectual disability (ID) and epilepsy. Schizophrenia (SCZ) has also recently been proposed to belong to NDDs. Relatively common causes of NDDs are copy number variations (CNVs), characterised by the gain or the loss of a portion of a chromosome. In this review, we focus on deletions and duplications at the 16p11.2 chromosomal region, associated with NDDs, ID, ASD but also epilepsy and SCZ. Some of the core phenotypes presented by human carriers could be recapitulated in animal and cellular models, which also highlighted prominent neurophysiological and signalling alterations underpinning 16p11.2 CNVs-associated phenotypes. In this review, we also provide an overview of the genes within the 16p11.2 locus, including those with partially known or unknown function as well as non-coding RNAs. A particularly interesting interplay was observed between MVP and MAPK3 in modulating some of the pathological phenotypes associated with the 16p11.2 deletion. Elucidating their role in intracellular signalling and their functional links will be a key step to devise novel therapeutic strategies for 16p11.2 CNVs-related syndromes.

Copy number variant risk loci for schizophrenia converge on the BDNF pathway

OBJECTIVES

Schizophrenia genetics is intricate, with common and rare variants' contributions not fully understood. Certain copy number variations (CNVs) elevate risk, pivotal for understanding mental disorder models. Despite CNVs' genome-wide distribution and variable gene and protein effects, we must explore beyond affected genes to interaction partners and molecular pathways.

METHODS

In this study, we developed machine-readable interactive pathways to enable analysis of functional effects of genes within CNV loci and identify ten common pathways across CNVs with high schizophrenia risk using the WikiPathways database, schizophrenia risk gene collections from GWAS studies, and a gene-disease association database.

RESULTS

For CNVs that are pathogenic for schizophrenia, we found overlapping pathways, including BDNF signalling, cytoskeleton, and inflammation. Common schizophrenia risk genes identified by different studies are found in all CNV pathways, but not enriched.

CONCLUSIONS

Our findings suggest that specific pathways - BDNF signalling - are critical contributors to schizophrenia risk conferred by rare CNVs. Our approach highlights the importance of not only investigating deleted or duplicated genes within pathogenic CNV loci, but also study their direct interaction partners, which may explain pleiotropic effects of CNVs on schizophrenia risk and offer a broader field for interventions.

Using rare genetic mutations to revisit structural brain asymmetry

Asymmetry between the left and right hemisphere is a key feature of brain organization. Hemispheric functional specialization underlies some of the most advanced human-defining cognitive operations, such as articulated language, perspective taking, or rapid detection of facial cues. Yet, genetic investigations into brain asymmetry have mostly relied on common variants, which typically exert small effects on brain-related phenotypes. Here, we leverage rare genomic deletions and duplications to study how genetic alterations reverberate in human brain and behavior. We designed a pattern-learning approach to dissect the impact of eight high-effect-size copy number variations (CNVs) on brain asymmetry in a multi-site cohort of 552 CNV carriers and 290 non-carriers. Isolated multivariate brain asymmetry patterns spotlighted regions typically thought to subserve lateralized functions, including language, hearing, as well as visual, face and word recognition. Planum temporale asymmetry emerged as especially susceptible to deletions and duplications of specific gene sets. Targeted analysis of common variants through genome-wide association study (GWAS) consolidated partly diverging genetic influences on the right versus left planum temporale structure. In conclusion, our gene-brain-behavior data fusion highlights the consequences of genetically controlled brain lateralization on uniquely human cognitive capacities.

KCTD10 regulates brain development by destabilizing brain disorder-associated protein KCTD13

KCTD10 belongs to the KCTD (potassiumchannel tetramerization domain) family, many members of which are associated with neuropsychiatric disorders. However, the biological function underlying the association with brain disorders remains to be explored. Here, we reveal that Kctd10 is highly expressed in neuronal progenitors and layer V neurons throughout brain development. Kctd10 deficiency triggers abnormal proliferation and differentiation of neuronal progenitors, reduced deep-layer (especially layer V) neurons, increased upper-layer neurons, and lowered brain size. Mechanistically, we screened and identified a unique KCTD10-interacting protein, KCTD13, associated with neurodevelopmental disorders. KCTD10 mediated the ubiquitination-dependent degradation of KCTD13 and KCTD10 ablation resulted in a considerable increase of KCTD13 expression in the developing cortex. KCTD13 overexpression in neuronal progenitors led to reduced proliferation and abnormal cell distribution, mirroring KCTD10 deficiency. Notably, mice with brain-specific Kctd10 knockout exhibited obvious motor deficits. This study uncovers the physiological function of KCTD10 and provides unique insights into the pathogenesis of neurodevelopmental disorders.

Pervasive alterations of intra-axonal volume and network organization in young children with a 16p11.2 deletion

Reciprocal Copy Number Variants (CNVs) at the 16p11.2 locus confer high risk for autism spectrum disorder (ASD) and other neurodevelopmental disorders (NDDs). Morphometric MRI studies have revealed large and pervasive volumetric alterations in carriers of a 16p11.2 deletion. However, the specific neuroanatomical mechanisms underlying such alterations, as well as their developmental trajectory, are still poorly understood. Here we explored differences in microstructural brain connectivity between 24 children carrying a 16p11.2 deletion and 66 typically developing (TD) children between 2 and 8 years of age. We found a large pervasive increase of intra-axonal volume widespread over a high number of white matter tracts. Such microstructural alterations in 16p11.2 deletion children were already present at an early age, and led to significant changes in the global efficiency and integration of brain networks mainly associated to language, motricity and socio-emotional behavior, although the widespread pattern made it unlikely to represent direct functional correlates. Our results shed light on the neuroanatomical basis of the previously reported increase of white matter volume, and align well with analogous evidence of altered axonal diameter and synaptic function in 16p11.2 mice models. We provide evidence of a prevalent mechanistic deviation from typical maturation of brain structural connectivity associated with a specific biological risk to develop ASD. Future work is warranted to determine how this deviation contributes to the emergence of symptoms observed in young children diagnosed with ASD and other NDDs.

A Comprehensive Analysis of the Structural Recognition between KCTD Proteins and Cullin 3

KCTD ((K)potassium Channel Tetramerization Domain-containing) proteins constitute an emerging class of proteins involved in fundamental physio-pathological processes. In these proteins, the BTB domain, which represents the defining element of the family, may have the dual role of promoting oligomerization and favoring functionally important partnerships with different interactors. Here, by exploiting the potential of recently developed methodologies for protein structure prediction, we report a comprehensive analysis of the interactions of all KCTD proteins with their most common partner Cullin 3 (Cul3). The data here presented demonstrate the impressive ability of this approach to discriminate between KCTDs that interact with Cul3 and those that do not. Indeed, reliable and stable models of the complexes were only obtained for the 15 members of the family that are known to interact with Cul3. The generation of three-dimensional models for all KCTD-Cul3 complexes provides interesting clues on the determinants of the structural basis of this partnership as clear structural differences emerged between KCTDs that bind or do not bind Cul3. Finally, the availability of accurate three-dimensional models for KCTD-Cul3 interactions may be valuable for the ad hoc design and development of compounds targeting specific KCTDs that are involved in several common diseases.

Parallelized engineering of mutational models using piggyBac transposon delivery of CRISPR libraries

New technologies and large-cohort studies have enabled novel variant discovery and association at unprecedented scale, yet functional characterization of these variants remains paramount to deciphering disease mechanisms. Approaches that facilitate parallelized genome editing of cells of interest or induced pluripotent stem cells (iPSCs) have become critical tools toward this goal. Here, we developed an approach that incorporates libraries of CRISPR-Cas9 guide RNAs (gRNAs) together with inducible Cas9 into a piggyBac (PB) transposon system to engineer dozens to hundreds of genomic variants in parallel against isogenic cellular backgrounds. This method empowers loss-of-function (LoF) studies through the introduction of insertions or deletions (indels) and copy-number variants (CNVs), though generating specific nucleotide changes is possible with prime editing. The ability to rapidly establish high-quality mutational models at scale will facilitate the development of isogenic cellular collections and catalyze comparative functional genomic studies investigating the roles of hundreds of genes and mutations in development and disease.

RNA methyltransferase SPOUT1/CENP-32 links mitotic spindle organization with the neurodevelopmental disorder SpADMiSS

SPOUT1/CENP-32 encodes a putative SPOUT RNA methyltransferase previously identified as a mitotic chromosome associated protein. SPOUT1/CENP-32 depletion leads to centrosome detachment from the spindle poles and chromosome misalignment. Aided by gene matching platforms, we identified 24 individuals with neurodevelopmental delays from 18 families with bi-allelic variants in SPOUT1/CENP-32 detected by exome/genome sequencing. Zebrafish spout1/cenp-32 mutants showed reduction in larval head size with concomitant apoptosis likely associated with altered cell cycle progression. In vivo complementation assays in zebrafish indicated that SPOUT1/CENP-32 missense variants identified in humans are pathogenic. Crystal structure analysis of SPOUT1/CENP-32 revealed that most disease-associated missense variants mapped to the catalytic domain. Additionally, SPOUT1/CENP-32 recurrent missense variants had reduced methyltransferase activity in vitro and compromised centrosome tethering to the spindle poles in human cells. Thus, SPOUT1/CENP-32 pathogenic variants cause an autosomal recessive neurodevelopmental disorder: SpADMiSS ( SPOUT1 Associated Development delay Microcephaly Seizures Short stature) underpinned by mitotic spindle organization defects and consequent chromosome segregation errors.

Integrating non-mammalian model organisms in the diagnosis of rare genetic diseases in humans

Next-generation sequencing technology has rapidly accelerated the discovery of genetic variants of interest in individuals with rare diseases. However, showing that these variants are causative of the disease in question is complex and may require functional studies. Use of non-mammalian model organisms - mainly fruitflies (Drosophila melanogaster), nematode worms (Caenorhabditis elegans) and zebrafish (Danio rerio) - enables the rapid and cost-effective assessment of the effects of gene variants, which can then be validated in mammalian model organisms such as mice and in human cells. By probing mechanisms of gene action and identifying interacting genes and proteins in vivo, recent studies in these non-mammalian model organisms have facilitated the diagnosis of numerous genetic diseases and have enabled the screening and identification of therapeutic options for patients. Studies in non-mammalian model organisms have also shown that the biological processes underlying rare diseases can provide insight into more common mechanisms of disease and the biological functions of genes. Here, we discuss the opportunities afforded by non-mammalian model organisms, focusing on flies, worms and fish, and provide examples of their use in the diagnosis of rare genetic diseases.

PLASTAMINATION: Outcomes on the Central Nervous System and Reproduction

BACKGROUND

Environmental exposures to non-biodegradable and biodegradable plastics are unavoidable. Microplastics (MPs) and nanoplastics (NPs) from the manufacturing of plastics (primary sources) and the degradation of plastic waste (secondary sources) can enter the food chain directly or indirectly and, passing biological barriers, could target both the brain and the gonads. Hence, the worldwide diffusion of environmental plastic contamination (PLASTAMINATION) in daily life may represent a possible and potentially serious risk to human health.

OBJECTIVE

This review provides an overview of the effects of non-biodegradable and the more recently introduced biodegradable MPs and NPs on the brain and brain-dependent reproductive functions, summarizing the molecular mechanisms and outcomes on nervous and reproductive organs. Data from in vitro, ex vivo, non-mammalian and mammalian animal models and epidemiological studies have been reviewed and discussed.

RESULTS

MPs and NPs from non-biodegradable plastics affect organs, tissues and cells from sensitive systems such as the brain and reproductive organs. Both MPs and NPs induce oxidative stress, chronic inflammation, energy metabolism disorders, mitochondrial dysfunction and cytotoxicity, which in turn are responsible for neuroinflammation, dysregulation of synaptic functions, metabolic dysbiosis, poor gamete quality, and neuronal and reproductive toxicity. In spite of this mechanistic knowledge gained from studies of non-biodegradable plastics, relatively little is known about the adverse effects or molecular mechanisms of MPs and NPs from biodegradable plastics.

CONCLUSION

The neurological and reproductive health risks of MPs/NPs exposure warrant serious consideration, and further studies on biodegradable plastics are recommended.

Prenatal prevalence and postnatal manifestations of 16p11.2 deletions: A new insights into neurodevelopmental disorders based on clinical investigations combined with multi-omics analysis

BACKGROUND

The 16p11.2 deletion is one of the most common genetic aetiologies of neurodevelopmental disorders (NDDs). The prenatal phenotype of 16p11.2 deletion and the potential mechanism associated with postnatal clinical manifestations were largely unknow. We revealed the developmental trajectories of 16p11.2 deletion from the prenatal to postnatal periods and to identify key signaling pathways and candidate genes contributing to neurodevelopmental abnormalities.

METHODS

In this 5-y retrospective cohort study, women with singleton pregnancies who underwent amniocentesis for chromosomal abnormalities were included. Test of copy-number variations (CNVs) involved single nucleotide polymorphism-array and CNV-seq was performed to detected 16p11.2 deletion. For infants born carrying the 16p11.2 deletion, neurological and intellectual evaluations using the Chinese version of the Gesell Development Scale. For patients observed to have vertebral malformations, Sanger sequencing for T-C-A haplotype of TBX6 was performed. For those infants with clinical manifestations, whole-exome sequencing was consecutively performed in trios to rule out single-gene diseases, and transcriptomics combined with untargeted metabolomics were performed.

RESULTS

The prevalence of 16p11.2 deletion was 0.063% (55/86,035) in the prenatal period. Up to 80% (20/25) of the 16p11.2 deletions were proven de novo by parental confirmation. Approximately half of 16p11.2 deletions (28/55) were detected with prenatal abnormal ultrasound findings. Vertebral malformations were identified as the most distinctive structural malformations and were enriched in fetuses with 16p11.2 deletions compared with controls (90.9‰ [5/55] vs. 8.4‰ [72/85,980]; P < 0.001). All 5 fetuses with vertebral malformations were confirmed to have the TBX6 haplotype of T-C-A. Overall, 47.6% (10/21) infants birthed were diagnosed with NDDs of different degrees. Language impairment was the predominant manifestation (7/10; 70.0%), followed by motor delay (5/10; 50%). Multi-omics analysis indicated that MAPK3 was the central hub of the differentially expressed gene (DEG) network. We firstly reported that histidine-associated metabolism may be the core metabolic pathway related to the 16p11.2 deletion.

CONCLUSION

We demonstrated the prenatal presentation, incomplete penetrance and variable expressivity of the 16p11.2 deletion. We identified vertebral malformations were the most distinctive prenatal phenotypes, and language impairment was the predominant postnatal manifestation. Most of the 16p11.2 deletion was de novo. Meanwhile, we suggested that MAPK3 and histidine-associated metabolism may contribute to neurodevelopmental abnormalities of 16p11.2 deletion.

Multi-omics analysis of zebrafish response to tick saliva reveals biological processes associated with alpha-Gal syndrome

The alpha-Gal syndrome (AGS) is a tick-borne allergy. A multi-omics approach was used to determine the effect of tick saliva and mammalian meat consumption on zebrafish gut transcriptome and proteome. Bioinformatics analysis using R software was focused on significant biological and metabolic pathway changes associated with AGS. Ortholog mapping identified highly concordant human ortholog genes for the detection of disease-enriched pathways. Tick saliva treatment increased zebrafish mortality, incidence of hemorrhagic type allergic reactions and changes in behavior and feeding patterns. Transcriptomics analysis showed downregulation of biological and metabolic pathways correlated with anti-alpha-Gal IgE and allergic reactions to tick saliva affecting blood circulation, cardiac and vascular smooth muscle contraction, behavior and sensory perception. Disease enrichment analysis revealed downregulated orthologous genes associated with human disorders affecting nervous, musculoskeletal, and cardiovascular systems. Proteomics analysis revealed suppression of pathways associated with immune system production of reactive oxygen species and cardiac muscle contraction. Underrepresented proteins were mainly linked to nervous and metabolic human disorders. Multi-omics data revealed inhibition of pathways associated with adrenergic signaling in cardiomyocytes, and heart and muscle contraction. Results identify tick saliva-related biological pathways supporting multisystemic organ involvement and linking α-Gal sensitization with other illnesses for the identification of potential disease biomarkers.

Dissecting the autism-associated 16p11.2 locus identifies multiple drivers in neuroanatomical phenotypes and unveils a male-specific role for the major vault protein

BACKGROUND

Using mouse genetic studies and systematic assessments of brain neuroanatomical phenotypes, we set out to identify which of the 30 genes causes brain defects at the autism-associated 16p11.2 locus.

RESULTS

We show that multiple genes mapping to this region interact to regulate brain anatomy, with female mice exhibiting far fewer brain neuroanatomical phenotypes. In male mice, among the 13 genes associated with neuroanatomical defects (Mvp, Ppp4c, Zg16, Taok2, Slx1b, Maz, Fam57b, Bola2, Tbx6, Qprt, Spn, Hirip3, and Doc2a), Mvp is the top driver implicated in phenotypes pertaining to brain, cortex, hippocampus, ventricles, and corpus callosum sizes. The major vault protein (MVP), the main component of the vault organelle, is a conserved protein found in eukaryotic cells, yet its function is not understood. Here, we find MVP expression highly specific to the limbic system and show that Mvp regulates neuronal morphology, postnatally and specifically in males. We also recapitulate a previously reported genetic interaction and show that Mvp+/-;Mapk3+/- mice exhibit behavioral deficits, notably decreased anxiety-like traits detected in the elevated plus maze and open field paradigms.

CONCLUSIONS

Our study highlights multiple gene drivers in neuroanatomical phenotypes, interacting with each other through complex relationships. It also provides the first evidence for the involvement of the major vault protein in the regulation of brain size and neuroanatomy, specifically in male mice.

Genomic and Reverse Translational Analysis Discloses a Role for Small GTPase RhoA Signaling in the Pathogenesis of Schizophrenia: Rho-Kinase as a Novel Drug Target

Schizophrenia is one of the most serious psychiatric disorders and is characterized by reductions in both brain volume and spine density in the frontal cortex. RhoA belongs to the RAS homolog (Rho) family and plays critical roles in neuronal development and structural plasticity via Rho-kinase. RhoA activity is regulated by GTPase-activating proteins (GAPs) and guanine nucleotide exchange factors (GEFs). Several variants in GAPs and GEFs associated with RhoA have been reported to be significantly associated with schizophrenia. Moreover, several mouse models carrying schizophrenia-associated gene variants involved in RhoA/Rho-kinase signaling have been developed. In this review, we summarize clinical evidence showing that variants in genes regulating RhoA activity are associated with schizophrenia. In the last half of the review, we discuss preclinical evidence indicating that RhoA/Rho-kinase is a potential therapeutic target of schizophrenia. In particular, Rho-kinase inhibitors exhibit anti-psychotic-like effects not only in Arhgap10 S490P/NHEJ mice, but also in pharmacologic models of schizophrenia (methamphetamine- and MK-801-treated mice). Accordingly, we propose that Rho-kinase inhibitors may have antipsychotic effects and reduce cognitive deficits in schizophrenia despite the presence or absence of genetic variants in small GTPase signaling pathways.

Microglial cannabinoid receptor type 1 mediates social memory deficits in mice produced by adolescent THC exposure and 16p11.2 duplication

Adolescent cannabis use increases the risk for cognitive impairments and psychiatric disorders. Cannabinoid receptor type 1 (Cnr1) is expressed not only in neurons and astrocytes, but also in microglia, which shape synaptic connections during adolescence. However, the role of microglia in mediating the adverse cognitive effects of delta-9-tetrahydrocannabinol (THC), the principal psychoactive constituent of cannabis, is not fully understood. Here, we report that in mice, adolescent THC exposure produces microglial apoptosis in the medial prefrontal cortex (mPFC), which was exacerbated in a model of 16p11.2 duplication, a representative copy number variation (CNV) risk factor for psychiatric disorders. These effects are mediated by microglial Cnr1, leading to reduction in the excitability of mPFC pyramidal-tract neurons and deficits in social memory in adulthood. Our findings suggest the microglial Cnr1 may contribute to adverse effect of cannabis exposure in genetically vulnerable individuals.

Screening of Tnfaip1-Interacting Proteins in Zebrafish Embryonic cDNA Libraries Using a Yeast Two-Hybrid System

TNFAIP1 regulates cellular biological functions, including DNA replication, DNA repair, and cell cycle, by binding to target proteins. Identification of Tnfaip1-interacting proteins contributes to the understanding of the molecular regulatory mechanisms of their biological functions. In this study, 48 hpf, 72 hpf, and 96 hpf wild-type zebrafish embryo mRNAs were used to construct yeast cDNA library. The library titer was 1.12 × 107 CFU/mL, the recombination rate was 100%, and the average length of the inserted fragments was greater than 1000 bp. A total of 43 potential interacting proteins of Tnfaip1 were identified using zebrafish Tnfaip1 as a bait protein. Utilizing GO functional annotation and KEGG signaling pathway analysis, we found that these interacting proteins are mainly involved in translation, protein catabolic process, ribosome assembly, cytoskeleton formation, amino acid metabolism, and PPAR signaling pathway. Further yeast spotting analyses identified four interacting proteins of Tnfaip1, namely, Ubxn7, Tubb4b, Rpl10, and Ybx1. The Tnfaip1-interacting proteins, screened from zebrafish embryo cDNA in this study, increased our understanding of the network of Tnfaip1-interacting proteins during the earliest embryo development and provided a molecular foundation for the future exploration of tnfaip1's biological functions.

Complexity in Genetic Epilepsies: A Comprehensive Review

Epilepsy is a highly prevalent neurological disorder, affecting between 5-8 per 1000 individuals and is associated with a lifetime risk of up to 3%. In addition to high incidence, epilepsy is a highly heterogeneous disorder, with variation including, but not limited to the following: severity, age of onset, type of seizure, developmental delay, drug responsiveness, and other comorbidities. Variable phenotypes are reflected in a range of etiologies including genetic, infectious, metabolic, immune, acquired/structural (resulting from, for example, a severe head injury or stroke), or idiopathic. This review will focus specifically on epilepsies with a genetic cause, genetic testing, and biomarkers in epilepsy.

oFlowSeq: a quantitative approach to identify protein coding mutations affecting cell type enrichment using mosaic CRISPR-Cas9 edited cerebral organoids

Cerebral organoids are comprised of diverse cell types found in the developing human brain, and can be leveraged in the identification of critical cell types perturbed by genetic risk variants in common, neuropsychiatric disorders. There is great interest in developing high-throughput technologies to associate genetic variants with cell types. Here, we describe a high-throughput, quantitative approach (oFlowSeq) by utilizing CRISPR-Cas9, FACS sorting, and next-generation sequencing. Using oFlowSeq, we found that deleterious mutations in autism-associated gene KCTD13 resulted in increased proportions of Nestin+ cells and decreased proportions of TRA-1-60+ cells within mosaic cerebral organoids. We further identified that a locus-wide CRISPR-Cas9 survey of another 18 genes in the 16p11.2 locus resulted in most genes with > 2% maximum editing efficiencies for short and long indels, suggesting a high feasibility for an unbiased, locus-wide experiment using oFlowSeq. Our approach presents a novel method to identify genotype-to-cell type imbalances in an unbiased, high-throughput, quantitative manner.

Microglial cannabinoid receptor type 1 mediates social memory deficits produced by adolescent THC exposure and 16p11.2 duplication

Adolescent cannabis use increases the risk for cognitive impairments and psychiatric disorders. Cannabinoid receptor type 1 (Cnr1) is expressed not only in neurons and astrocytes, but also in microglia, which shape synaptic connections during adolescence. Nonetheless, until now, the role of microglia in mediating the adverse cognitive effects of delta-9-tetrahydrocannabinol (THC), the principal psychoactive constituent of cannabis, has been unexplored. Here, we report that adolescent THC exposure produces microglial apoptosis in the medial prefrontal cortex (mPFC), which was exacerbated in the mouse model of 16p11.2 duplication, a representative copy number variation (CNV) risk factor for psychiatric disorders. These effects are mediated by microglial Cnr1, leading to reduction in the excitability of mPFC pyramidal-tract neurons and deficits in social memory in adulthood. Our findings highlight the importance of microglial Cnr1 to produce the adverse effect of cannabis exposure in genetically vulnerable individuals.

Genome-wide identification and phenotypic characterization of seizure-associated copy number variations in 741,075 individuals

Copy number variants (CNV) are established risk factors for neurodevelopmental disorders with seizures or epilepsy. With the hypothesis that seizure disorders share genetic risk factors, we pooled CNV data from 10,590 individuals with seizure disorders, 16,109 individuals with clinically validated epilepsy, and 492,324 population controls and identified 25 genome-wide significant loci, 22 of which are novel for seizure disorders, such as deletions at 1p36.33, 1q44, 2p21-p16.3, 3q29, 8p23.3-p23.2, 9p24.3, 10q26.3, 15q11.2, 15q12-q13.1, 16p12.2, 17q21.31, duplications at 2q13, 9q34.3, 16p13.3, 17q12, 19p13.3, 20q13.33, and reciprocal CNVs at 16p11.2, and 22q11.21. Using genetic data from additional 248,751 individuals with 23 neuropsychiatric phenotypes, we explored the pleiotropy of these 25 loci. Finally, in a subset of individuals with epilepsy and detailed clinical data available, we performed phenome-wide association analyses between individual CNVs and clinical annotations categorized through the Human Phenotype Ontology (HPO). For six CNVs, we identified 19 significant associations with specific HPO terms and generated, for all CNVs, phenotype signatures across 17 clinical categories relevant for epileptologists. This is the most comprehensive investigation of CNVs in epilepsy and related seizure disorders, with potential implications for clinical practice.

Structural basis for the ubiquitination of G protein βγ subunits by KCTD5/Cullin3 E3 ligase

G protein-coupled receptor (GPCR) signaling is precisely controlled to avoid overstimulation that results in detrimental consequences. Gβγ signaling is negatively regulated by a Cullin3 (Cul3)-dependent E3 ligase, KCTD5, which triggers ubiquitination and degradation of free Gβγ. Here, we report the cryo-electron microscopy structures of the KCTD5-Gβγ fusion complex and the KCTD7-Cul3 complex. KCTD5 in pentameric form engages symmetrically with five copies of Gβγ through its C-terminal domain. The unique pentameric assembly of the KCTD5/Cul3 E3 ligase places the ubiquitin-conjugating enzyme (E2) and the modification sites of Gβγ in close proximity and allows simultaneous transfer of ubiquitin from E2 to five Gβγ subunits. Moreover, we show that ubiquitination of Gβγ by KCTD5 is important for fine-tuning cyclic adenosine 3´,5´-monophosphate signaling of GPCRs. Our studies provide unprecedented insights into mechanisms of substrate recognition by unusual pentameric E3 ligases and highlight the KCTD family as emerging regulators of GPCR signaling.

KCTD13-mediated ubiquitination and degradation of GluN1 regulates excitatory synaptic transmission and seizure susceptibility

Temporal lobe epilepsy (TLE) is the most common and severe form of epilepsy in adults; however, its underlying pathomechanisms remain elusive. Dysregulation of ubiquitination is increasingly recognized to contribute to the development and maintenance of epilepsy. Herein, we observed for the first time that potassium channel tetramerization domain containing 13 (KCTD13) protein, a substrate-specific adapter for cullin3-based E3 ubiquitin ligase, was markedly down-regulated in the brain tissue of patients with TLE. In a TLE mouse model, the protein expression of KCTD13 dynamically changed during epileptogenesis. Knockdown of KCTD13 in the mouse hippocampus significantly enhanced seizure susceptibility and severity, whereas overexpression of KCTD13 showed the opposite effect. Mechanistically, GluN1, an obligatory subunit of N-methyl-D-aspartic acid receptors (NMDARs), was identified as a potential substrate protein of KCTD13. Further investigation revealed that KCTD13 facilitates lysine-48-linked polyubiquitination of GluN1 and its degradation through the ubiquitin-proteasome pathway. Besides, the lysine residue 860 of GluN1 is the main ubiquitin site. Importantly, dysregulation of KCTD13 affected membrane expression of glutamate receptors and impaired glutamate synaptic transmission. Systemic administration of the NMDAR inhibitor memantine significantly rescued the epileptic phenotype aggravated by KCTD13 knockdown. In conclusion, our results demonstrated an unrecognized pathway of KCTD13-GluN1 in epilepsy, suggesting KCTD13 as a potential neuroprotective therapeutic target for epilepsy.

Combinations of genes at the 16p11.2 and 22q11.2 CNVs contribute to neurobehavioral traits

The 16p11.2 and 22q11.2 copy number variants (CNVs) are associated with neurobehavioral traits including autism spectrum disorder (ASD), schizophrenia, bipolar disorder, obesity, and intellectual disability. Identifying specific genes contributing to each disorder and dissecting the architecture of CNV-trait association has been difficult, inspiring hypotheses of more complex models, such as multiple genes acting together. Using multi-tissue data from the GTEx consortium, we generated pairwise expression imputation models for CNV genes and then applied these elastic net models to GWAS for: ASD, bipolar disorder, schizophrenia, BMI (obesity), and IQ (intellectual disability). We compared the variance in these five traits explained by gene pairs with the variance explained by single genes and by traditional interaction models. We also modeled polygene region-wide effects using summed predicted expression ranks across many genes to create a regionwide score. We found that in all CNV-trait pairs except for bipolar disorder at 22q11.2, pairwise effects explain more variance than single genes. Pairwise model superiority was specific to the CNV region for all 16p11.2 traits and ASD at 22q11.2. We identified novel individual genes over-represented in top pairs that did not show single-gene signal. We also found that BMI and IQ have significant regionwide association with both CNV regions. Overall, we observe that genetic architecture differs by trait and region, but 9/10 CNV-trait combinations demonstrate evidence for multigene contribution, and for most of these, the importance of combinatorial models appears unique to CNV regions. Our results suggest that mechanistic insights for CNV pathology may require combinational models.

CRISPR/Cas9-Mediated Knockout of tnfaip1 in Zebrafish Plays a Role in Early Development

TNF α-induced protein 1 (TNFAIP1) was first identified in human umbilical vein endothelial cells and can be induced by tumor necrosis factor α (TNFα). Early studies have found that TNFAIP1 is involved in the development of many tumors and is closely associated with the neurological disorder Alzheimer's disease. However, little is known about the expression pattern of TNFAIP1 under physiological conditions and its function during embryonic development. In this study, we used zebrafish as a model to illustrate the early developmental expression pattern of tnfaip1 and its role in early development. First, we examined the expression pattern of tnfaip1 during early zebrafish development using quantitative real-time PCR and whole mount in situ hybridization and found that tnfaip1 was highly expressed in early embryonic development and, subsequently, expression became localized to anterior embryonic structures. To investigate the function of tnfaip1 during early development, we constructed a model of a stably inherited tnfaip1 mutant using the CRISPR/Cas9 system. Tnfaip1 mutant embryos showed significant developmental delays as well as microcephaly and microphthalmia. At the same time, we found decreased expression of the neuronal marker genes tuba1b, neurod1, and ccnd1 in tnfaip1 mutants. Analysis of transcriptome sequencing data revealed altered expression of the embryonic development related genes dhx40, hspa13, tnfrsf19, nppa, lrp2b, hspb9, clul1, zbtb47a, cryba1a, and adgrg4a in the tnfaip1 mutants. These findings suggest an important role for tnfaip1 in the early development of zebrafish.

Translational Neuroscience Approaches to Understanding Autism

While autism spectrum disorder affects nearly 2% of children in the United States, little is known with certainty concerning the etiologies and brain systems involved. This is due, in part, to the substantial heterogeneity in the presentation of the core symptoms of autism as well as the great number of co-occurring conditions that are common in autistic individuals. Understanding the neurobiology of autism is further hampered by the limited availability of postmortem brain tissue to determine the cellular and molecular alterations that take place in the autistic brain. Animal models therefore provide great translational value in helping to define the neural systems that constitute the social brain and mediate repetitive behaviors or interests. If they are based on genetic or environmental factors that contribute to autism, organisms from flies to nonhuman primates may serve as models of the neural structure or function of the autistic brain. Ultimately, successful models can also be employed to test the safety and effectiveness of potential therapeutics. This is an overview of the major animal species that are currently used as models of autism, including an appraisal of the advantages and limitations of each.

Calpain activity is negatively regulated by a KCTD7-Cullin-3 complex via non-degradative ubiquitination

Calpains are a class of non-lysosomal cysteine proteases that exert their regulatory functions via limited proteolysis of their substrates. Similar to the lysosomal and proteasomal systems, calpain dysregulation is implicated in the pathogenesis of neurodegenerative disease and cancer. Despite intensive efforts placed on the identification of mechanisms that regulate calpains, however, calpain protein modifications that regulate calpain activity are incompletely understood. Here we show that calpains are regulated by KCTD7, a cytosolic protein of previously uncharacterized function whose pathogenic mutations result in epilepsy, progressive ataxia, and severe neurocognitive deterioration. We show that KCTD7 works in complex with Cullin-3 and Rbx1 to execute atypical, non-degradative ubiquitination of calpains at specific sites (K398 of calpain 1, and K280 and K674 of calpain 2). Experiments based on single-lysine mutants of ubiquitin determined that KCTD7 mediates ubiquitination of calpain 1 via K6-, K27-, K29-, and K63-linked chains, whereas it uses K6-mediated ubiquitination to modify calpain 2. Loss of KCTD7-mediated ubiquitination of calpains led to calpain hyperactivation, aberrant cleavage of downstream targets, and caspase-3 activation. CRISPR/Cas9-mediated knockout of Kctd7 in mice phenotypically recapitulated human KCTD7 deficiency and resulted in calpain hyperactivation, behavioral impairments, and neurodegeneration. These phenotypes were largely prevented by pharmacological inhibition of calpains, thus demonstrating a major role of calpain dysregulation in KCTD7-associated disease. Finally, we determined that Cullin-3-KCTD7 mediates ubiquitination of all ubiquitous calpains. These results unveil a novel mechanism and potential target to restrain calpain activity in human disease and shed light on the molecular pathogenesis of KCTD7-associated disease.

The Oncosuppressive Properties of KCTD1: Its Role in Cell Growth and Mobility

The KCTD protein family is traditionally regarded as proteins that play key roles in neurological physiopathology. However, new studies are increasingly demonstrating their involvement in many other biological processes, including cancers. This is particularly evident for KCTD proteins not involved in protein ubiquitination and degradation, such as KCTD1. We explored the role of KCTD1 in colorectal cancer by knocking down this protein in the human colon adenocarcinoma cell line, SW480. We re-assessed its ability to downregulate β-catenin, a central actor in the WNT/β-catenin signalling pathway. Interestingly, opposite effects are observed when the protein is upregulated in CACO2 colorectal cancer cells. Moreover, interrogation of the TCGA database indicates that KCTD1 downregulation is associated with β-catenin overexpression in colorectal cancer patients. Indeed, knocking down KCTD1 in SW480 cells led to a significant increase in their motility and stemness, two important tumorigenesis traits, suggesting an oncosuppressor role for KCTD1. It is worth noting that similar effects are induced on colorectal cancer cells by the misregulation of KCTD12, a protein that is distantly related to KCTD1. The presented results further expand the spectrum of KCTD1 involvement in apparently unrelated physiopathological processes. The similar effects produced on colorectal cancer cell lines by KCTD1 and KCTD12 suggest novel, previously unreported analogous activities among members of the KCTD protein family.

Structure-based developmental toxicity and ASD-phenotypes of bisphenol A analogues in embryonic zebrafish

Autism spectrum disorder (ASD) is a neurodevelopmental disorder that has become more prevalent in recent years. Environmental endocrine disruptor bisphenol A (BPA) has been linked to ASD. BPA analogues (BPs) are structure-modified substitutes widely used as safer alternatives in consumer products, yet few studies have explored the developmental neurotoxicity (DNT) of BPA analogues. In the present study, we used the larval zebrafish model to assess the DNT effects of BPA and its analogues. Our results showed that many BPA analogues are more toxic than BPA in the embryonic zebrafish assay regarding teratogenic effect and mortality, which may partially due to differences in lipophilicity and/or different substitutes of structural function groups such as CF3, benzene, or cyclohexane. At sublethal concentrations, zebrafish embryos exposed to BPA or BPs also displayed reduced prosocial behavior in later larval development, evidenced by increased nearest neighbor distance (NND) and the interindividual distance (IID) in shoaling, which appears to be structurally independent. An in-depth analysis of BPA, bisphenol F (BPF), and bisphenol S (BPS) revealed macrocephaly and ASD-like behavioral deficits resulting from exposures to sublethal concentrations of these chemicals. The ASD-like behavioral deficits were characterized by hyperactivity, increased anxiety-like behavior, and decreased social contact. Mechanistically, accelerated neurogenesis that manifested by increased cell proliferation, the proportion of newborn mature neurons, and the number of neural stem cells in proliferation, as well as upregulated genes related to the K⁺ channels, may have contributed to the observed ASD-like morphological and behavioral alterations. Our findings indicate that BPF and BPS may also pose significant risks to ASD development in humans and highlight the importance of a comprehensive assessment of DNT effects for all BPA analogues in the future.

Convergent coexpression of autism-associated genes suggests some novel risk genes may not be detectable in large-scale genetic studies

Autism spectrum disorder (ASD) is a heritable neurodevelopmental disorder characterized by deficits in social interactions and communication. Protein-altering variants in many genes have been shown to contribute to ASD; however, understanding the convergence across many genes remains a challenge. We demonstrate that coexpression patterns from 993 human postmortem brains are significantly correlated with the transcriptional consequences of CRISPR perturbations in human neurons. Across 71 ASD risk genes, there was significant tissue-specific convergence implicating synaptic pathways. Tissue-specific convergence was further demonstrated across schizophrenia and atrial fibrillation risk genes. The degree of ASD convergence was significantly correlated with ASD association from rare variation and differential expression in ASD brains. Positively convergent genes showed intolerance to functional mutations and had shorter coding lengths than known risk genes even after removing association with ASD. These results indicate that convergent coexpression can identify potentially novel genes that are unlikely to be discovered by sequencing studies.

16p11.2 deletion accelerates subpallial maturation and increases variability in human iPSC-derived ventral telencephalic organoids

Inhibitory interneurons regulate cortical circuit activity, and their dysfunction has been implicated in autism spectrum disorder (ASD). 16p11.2 microdeletions are genetically linked to 1% of ASD cases. However, few studies investigate the effects of this microdeletion on interneuron development. Using ventral telencephalic organoids derived from human induced pluripotent stem cells, we have investigated the effect of this microdeletion on organoid size, progenitor proliferation and organisation into neural rosettes, ganglionic eminence marker expression at early developmental timepoints, and expression of the neuronal marker NEUN at later stages. At early stages, deletion organoids exhibited greater variations in size with concomitant increases in relative neural rosette area and the expression of the ventral telencephalic marker COUPTFII, with increased variability in these properties. Cell cycle analysis revealed an increase in total cell cycle length caused primarily by an elongated G1 phase, the duration of which also varied more than normal. At later stages, deletion organoids increased their NEUN expression. We propose that 16p11.2 microdeletions increase developmental variability and may contribute to ASD aetiology by lengthening the cell cycle of ventral progenitors, promoting premature differentiation into interneurons.

Stage-Related Neurotoxicity of BPA in the Development of Zebrafish Embryos

Bisphenol A (BPA) is one of the most widely produced chemicals in the world used in the production of epoxy resins and polycarbonate plastics. BPA is easily migrated from the outer packaging to the contents. Due to the lipophilic property, BPA is easily accumulated in organisms. Perinatal low-dose BPA exposure alters brain neural development in later generations. In this study, after BPA treatment, the spontaneous movement of zebrafish larvae from the cleavage period to the segmentation period (1-24 hpf) was significantly decreased, with speed decreasing by 18.97% and distance decreasing between 18.4 and 29.7% compared to controls. Transcriptomics analysis showed that 131 genes were significantly differentially expressed in the exposed group during the 1-24 hpf period, among which 39 genes were significantly upregulated and 92 genes were significantly downregulated. The GO enrichment analysis, gene function analysis and real-time quantitative PCR of differentially expressed genes showed that the mRNA level of guanine deaminase (cypin) decreased significantly in the 1-24 hpf period. Moreover, during the 1-24 hpf period, BPA exposure reduced guanine deaminase activity. Therefore, we confirmed that cypin is a key sensitive gene for BPA during this period. Finally, the cypin mRNA microinjection verified that the cypin level of zebrafish larvae was restored, leading to the restoration of the locomotor activity. Taken together, the current results show that the sensitive period of BPA to zebrafish embryos is from the cleavage period to the segmentation period (1-24 hpf), and cypin is a potential target for BPA-induced neurodevelopmental toxicity. This study provides a potential sensitive period and a potential target for the deep understanding of neurodevelopmental toxicity mechanisms caused by BPA.

Identifying foetal forebrain interneurons as a target for monogenic autism risk factors and the polygenic 16p11.2 microdeletion

BACKGROUND

Autism spectrum condition or 'autism' is associated with numerous genetic risk factors including the polygenic 16p11.2 microdeletion. The balance between excitatory and inhibitory neurons in the cerebral cortex is hypothesised to be critical for the aetiology of autism making improved understanding of how risk factors impact on the development of these cells an important area of research. In the current study we aim to combine bioinformatics analysis of human foetal cerebral cortex gene expression data with anatomical and electrophysiological analysis of a 16p11.2^+/- rat model to investigate how genetic risk factors impact on inhibitory neuron development.

METHODS

We performed bioinformatics analysis of single cell transcriptomes from gestational week (GW) 8-26 human foetal prefrontal cortex and anatomical and electrophysiological analysis of 16p11.2^+/- rat cerebral cortex and hippocampus at post-natal day (P) 21.

RESULTS

We identified a subset of human interneurons (INs) first appearing at GW23 with enriched expression of a large fraction of risk factor transcripts including those expressed from the 16p11.2 locus. This suggests the hypothesis that these foetal INs are vulnerable to mutations causing autism. We investigated this in a rat model of the 16p11.2 microdeletion. We found no change in the numbers or position of either excitatory or inhibitory neurons in the somatosensory cortex or CA1 of 16p11.2^+/- rats but found that CA1 Sst INs were hyperexcitable with an enlarged axon initial segment, which was not the case for CA1 pyramidal cells.

LIMITATIONS

The human foetal gene expression data was acquired from cerebral cortex between gestational week (GW) 8 to 26. We cannot draw inferences about potential vulnerabilities to genetic autism risk factors for cells not present in the developing cerebral cortex at these stages. The analysis 16p11.2^+/- rat phenotypes reported in the current study was restricted to 3-week old (P21) animals around the time of weaning and to a single interneuron cell-type while in human 16p11.2 microdeletion carriers symptoms likely involve multiple cell types and manifest in the first few years of life and on into adulthood.

CONCLUSIONS

We have identified developing interneurons in human foetal cerebral cortex as potentially vulnerable to monogenic autism risk factors and the 16p11.2 microdeletion and report interneuron phenotypes in post-natal 16p11.2^+/- rats.

Zebrafish: A Model Deciphering the Impact of Flavonoids on Neurodegenerative Disorders

Over the past century, advances in biotechnology, biochemistry, and pharmacognosy have spotlighted flavonoids, polyphenolic secondary metabolites that have the ability to modulate many pathways involved in various biological mechanisms, including those involved in neuronal plasticity, learning, and memory. Moreover, flavonoids are known to impact the biological processes involved in developing neurodegenerative diseases, namely oxidative stress, neuroinflammation, and mitochondrial dysfunction. Thus, several flavonoids could be used as adjuvants to prevent and counteract neurodegenerative disorders such as Alzheimer's and Parkinson's diseases. Zebrafish is an interesting model organism that can offer new opportunities to study the beneficial effects of flavonoids on neurodegenerative diseases. Indeed, the high genome homology of 70% to humans, the brain organization largely similar to the human brain as well as the similar neuroanatomical and neurochemical processes, and the high neurogenic activity maintained in the adult brain makes zebrafish a valuable model for the study of human neurodegenerative diseases and deciphering the impact of flavonoids on those disorders.

Loss of autism-candidate CHD8 perturbs neural crest development and intestinal homeostatic balance

Individuals with mutations in CHD8 present with gastrointestinal complaints, yet the underlying mechanisms are understudied. Here, using a stable constitutive chd8 mutant zebrafish model, we found that the loss of chd8 leads to a reduced number of vagal neural crest cells (NCCs), enteric neural and glial progenitors, emigrating from the neural tube, and that their early migration capability was altered. At later stages, although the intestinal colonization by NCCs was complete, we found the decreased numbers of both serotonin-producing enterochromaffin cells and NCC-derived serotonergic neurons, suggesting an intestinal hyposerotonemia in the absence of chd8 Furthermore, transcriptomic analyses revealed an altered expression of key receptors and enzymes in serotonin and acetylcholine signaling pathways. The tissue examination of chd8 mutants revealed a thinner intestinal epithelium accompanied by an accumulation of neutrophils and the decreased numbers of goblet cells and eosinophils. Last, single-cell sequencing of whole intestines showed a global disruption of the immune balance with a perturbed expression of inflammatory interleukins and changes in immune cell clusters. Our findings propose a causal developmental link between chd8, NCC development, intestinal homeostasis, and autism-associated gastrointestinal complaints.

The polymerase δ-interacting protein family and their emerging roles in diseases

The polymerase δ-interacting protein (POLDIP) family is a new family that can interact with DNA polymerase δ (delta). The members of the POLDIP family include POLDIP1, POLDIP2, and POLDIP3. Screened by the two-hybrid method, POLDIP1, POLDIP2, and POLDIP3 were initially discovered and named for their ability to bind to the p50 subunit of DNA polymerase δ. Recent studies have confirmed that POLDIPs are involved in the regulation of signal transduction pathways in neurodevelopment, neuropsychiatric diseases, cardiovascular diseases, tumors, and other diseases. However, each protein participates in different signaling pathways. In this review, we elucidate upon the family in terms of their genes and protein structures, their biological functions, in addition to the pathways that they are involved in during the development of diverse diseases. Finally, to provide new insights to the scientific community, we used the TCGA database to analyze and summarize the gene expressions of POLDIP family members in various tumors, as well as the correlations between their expressions and the overall survival times of tumor patients. Our data summary will give researchers working on cancer new concepts.

Statistical and functional convergence of common and rare genetic influences on autism at chromosome 16p

The canonical paradigm for converting genetic association to mechanism involves iteratively mapping individual associations to the proximal genes through which they act. In contrast, in the present study we demonstrate the feasibility of extracting biological insights from a very large region of the genome and leverage this strategy to study the genetic influences on autism. Using a new statistical approach, we identified the 33-Mb p-arm of chromosome 16 (16p) as harboring the greatest excess of autism's common polygenic influences. The region also includes the mechanistically cryptic and autism-associated 16p11.2 copy number variant. Analysis of RNA-sequencing data revealed that both the common polygenic influences within 16p and the 16p11.2 deletion were associated with decreased average gene expression across 16p. The transcriptional effects of the rare deletion and diffuse common variation were correlated at the level of individual genes and analysis of Hi-C data revealed patterns of chromatin contact that may explain this transcriptional convergence. These results reflect a new approach for extracting biological insight from genetic association data and suggest convergence of common and rare genetic influences on autism at 16p.

Gene dosage changes in KCTD13 result in penile and testicular anomalies via diminished androgen receptor function

Despite the high prevalence of hypospadias and cryptorchidism, the genetic basis for these conditions is only beginning to be understood. Using array-comparative-genomic-hybridization (aCGH), potassium-channel-tetramerization-domain-containing-13 (KCTD13) encoded at 16p11.2 was identified as a candidate gene involved in hypospadias, cryptorchidism and other genitourinary (GU) tract anomalies. Copy number variants (CNVs) at 16p11.2 are among the most common syndromic genomic variants identified to date. Many patients with CNVs at this locus exhibit GU and/or neurodevelopmental phenotypes. KCTD13 encodes a substrate-specific adapter of a BCR (BTB-CUL3-RBX1) E3-ubiquitin-protein-ligase complex (BCR (BTB-CUL3-RBX1) E3-ubiquitin-protein-ligase complex (B-cell receptor (BCR) [BTB (the BTB domain is a conserved motif involved in protein-protein interactions) Cullin3 complex RING protein Rbx1] E3-ubiqutin-protein-ligase complex), which has essential roles in the regulation of cellular cytoskeleton, migration, proliferation, and neurodevelopment; yet its role in GU development is unknown. The prevalence of KCTD13 CNVs in patients with GU anomalies (2.58%) is significantly elevated when compared with patients without GU anomalies or in the general population (0.10%). KCTD13 is robustly expressed in the developing GU tract. Loss of KCTD13 in cell lines results in significantly decreased levels of nuclear androgen receptor (AR), suggesting that loss of KCTD13 affects AR sub-cellular localization. Kctd13 haploinsufficiency and homozygous deletion in mice cause a significant increase in the incidence of cryptorchidism and micropenis. KCTD13-deficient mice exhibit testicular and penile abnormalities together with significantly reduced levels of nuclear AR and SOX9. In conclusion, gene-dosage changes of murine Kctd13 diminish nuclear AR sub-cellular localization, as well as decrease SOX9 expression levels which likely contribute in part to the abnormal GU tract development in Kctd13 mouse models and in patients with CNVs in KCTD13.

Tissue- and cell-type-specific molecular and functional signatures of 16p11.2 reciprocal genomic disorder across mouse brain and human neuronal models

Chromosome 16p11.2 reciprocal genomic disorder, resulting from recurrent copy-number variants (CNVs), involves intellectual disability, autism spectrum disorder (ASD), and schizophrenia, but the responsible mechanisms are not known. To systemically dissect molecular effects, we performed transcriptome profiling of 350 libraries from six tissues (cortex, cerebellum, striatum, liver, brown fat, and white fat) in mouse models harboring CNVs of the syntenic 7qF3 region, as well as cellular, transcriptional, and single-cell analyses in 54 isogenic neural stem cell, induced neuron, and cerebral organoid models of CRISPR-engineered 16p11.2 CNVs. Transcriptome-wide differentially expressed genes were largely tissue-, cell-type-, and dosage-specific, although more effects were shared between deletion and duplication and across tissue than expected by chance. The broadest effects were observed in the cerebellum (2,163 differentially expressed genes), and the greatest enrichments were associated with synaptic pathways in mouse cerebellum and human induced neurons. Pathway and co-expression analyses identified energy and RNA metabolism as shared processes and enrichment for ASD-associated, loss-of-function constraint, and fragile X messenger ribonucleoprotein target gene sets. Intriguingly, reciprocal 16p11.2 dosage changes resulted in consistent decrements in neurite and electrophysiological features, and single-cell profiling of organoids showed reciprocal alterations to the proportions of excitatory and inhibitory GABAergic neurons. Changes both in neuronal ratios and in gene expression in our organoid analyses point most directly to calretinin GABAergic inhibitory neurons and the excitatory/inhibitory balance as targets of disruption that might contribute to changes in neurodevelopmental and cognitive function in 16p11.2 carriers. Collectively, our data indicate the genomic disorder involves disruption of multiple contributing biological processes and that this disruption has relative impacts that are context specific.

A Mini-Review Regarding the Modalities to Study Neurodevelopmental Disorders-Like Impairments in Zebrafish—Focussing on Neurobehavioural and Psychological Responses

Neurodevelopmental disorders (NDDs) are complex disorders which can be associated with many comorbidities and exhibit multifactorial-dependent phenotypes. An important characteristic is represented by the early onset of the symptoms, during childhood or young adulthood, with a great impact on the socio-cognitive functioning of the affected individuals. Thus, the aim of our review is to describe and to argue the necessity of early developmental stages zebrafish models, focusing on NDDs, especially autism spectrum disorders (ASD) and also on schizophrenia. The utility of the animal models in NDDs or schizophrenia research remains quite controversial. Relevant discussions can be opened regarding the specific characteristics of the animal models and the relationship with the etiologies, physiopathology, and development of these disorders. The zebrafish models behaviors displayed as early as during the pre-hatching embryo stage (locomotor activity prone to repetitive behavior), and post-hatching embryo stage, such as memory, perception, affective-like, and social behaviors can be relevant in ASD and schizophrenia research. The neurophysiological processes impaired in both ASD and schizophrenia are generally highly conserved across all vertebrates. However, the relatively late individual development and conscious social behavior exhibited later in the larval stage are some of the most important limitations of these model animal species.

A cross-disorder dosage sensitivity map of the human genome

Rare copy-number variants (rCNVs) include deletions and duplications that occur infrequently in the global human population and can confer substantial risk for disease. In this study, we aimed to quantify the properties of haploinsufficiency (i.e., deletion intolerance) and triplosensitivity (i.e., duplication intolerance) throughout the human genome. We harmonized and meta-analyzed rCNVs from nearly one million individuals to construct a genome-wide catalog of dosage sensitivity across 54 disorders, which defined 163 dosage sensitive segments associated with at least one disorder. These segments were typically gene dense and often harbored dominant dosage sensitive driver genes, which we were able to prioritize using statistical fine-mapping. Finally, we designed an ensemble machine-learning model to predict probabilities of dosage sensitivity (pHaplo & pTriplo) for all autosomal genes, which identified 2,987 haploinsufficient and 1,559 triplosensitive genes, including 648 that were uniquely triplosensitive. This dosage sensitivity resource will provide broad utility for human disease research and clinical genetics.

Autism NPCs from both idiopathic and CNV 16p11.2 deletion patients exhibit dysregulation of proliferation and mitogenic responses

Neural precursor cell (NPC) dysfunction has been consistently implicated in autism. Induced pluripotent stem cell (iPSC)-derived NPCs from two autism groups (three idiopathic [I-ASD] and two 16p11.2 deletion [16pDel]) were used to investigate if proliferation is commonly disrupted. All five individuals display defects, with all three macrocephalic individuals (two 16pDel, one I-ASD) exhibiting hyperproliferation and the other two I-ASD subjects displaying hypoproliferation. NPCs were challenged with bFGF, and all hyperproliferative NPCs displayed blunted responses, while responses were increased in hypoproliferative cells. mRNA expression studies suggest that different pathways can result in similar proliferation phenotypes. Since 16pDel deletes MAPK3, P-ERK was measured. P-ERK is decreased in hyperproliferative but increased in hypoproliferative NPCs. While these P-ERK changes are not responsible for the phenotypes, P-ERK and bFGF response are inversely correlated with the defects. Finally, we analyzed iPSCs and discovered that 16pDel displays hyperproliferation, while idiopathic iPSCs were normal. These data suggest that NPC proliferation defects are common in ASD.

Orgo-Seq integrates single-cell and bulk transcriptomic data to identify cell type specific-driver genes associated with autism spectrum disorder

Cerebral organoids can be used to gain insights into cell type specific processes perturbed by genetic variants associated with neuropsychiatric disorders. However, robust and scalable phenotyping of organoids remains challenging. Here, we perform RNA sequencing on 71 samples comprising 1,420 cerebral organoids from 25 donors, and describe a framework (Orgo-Seq) to integrate bulk RNA and single-cell RNA sequence data. We apply Orgo-Seq to 16p11.2 deletions and 15q11-13 duplications, two loci associated with autism spectrum disorder, to identify immature neurons and intermediate progenitor cells as critical cell types for 16p11.2 deletions. We further applied Orgo-Seq to identify cell type-specific driver genes. Our work presents a quantitative phenotyping framework to integrate multi-transcriptomic datasets for the identification of cell types and cell type-specific co-expressed driver genes associated with neuropsychiatric disorders.

Clinical Implications of Chromosome 16 Copy Number Variation

Chromosome 16 is one of the gene-rich chromosomes; however, approximately 10% of the chromosome 16 sequence is composed of segmental copies, which renders this chromosome instable and predisposes it to rearrangements via frequent nonallelic homologous recombination. Microarray technologies have enabled the analysis of copy number variations (CNV), which may be associated with the risk of developing complex diseases. Through comparative genomic hybridisation in 1,298 patients, we detected 18 cases with chromosome 16 CNV. We identified 2recurrent CNV regions, including 1 at 16p13.11 in 4 patients and another at 16p11.2 in 7 patients. We also detected atypical chromosome 16 rearrangements in 7 patients. Furthermore, we noted an increased frequency of co-occurring genomic changes, supporting the two-hit hypothesis to explain the phenotypic variability in the clinical presentation of CNV syndromes. Our findings can contribute to the creation of a chromosome 16 disease map based on regions that may be associated with disease development.

Imaging neural circuit pathology of autism spectrum disorders: autism-associated genes, animal models and the application of in vivo two-photon imaging

Recent advances in human genetics identified genetic variants involved in causing autism spectrum disorders (ASDs). Mouse models that mimic mutations found in patients with ASD exhibit behavioral phenotypes consistent with ASD symptoms. These mouse models suggest critical biological factors of ASD etiology. Another important implication of ASD genetics is the enrichment of ASD risk genes in molecules involved in developing synapses and regulating neural circuit function. Sophisticated in vivo imaging technologies applied to ASD mouse models identify common synaptic impairments in the neocortex, with genetic-mutation-specific defects in local neural circuits. In this article, we review synapse- and circuit-level phenotypes identified by in vivo two-photon imaging in multiple mouse models of ASD and discuss the contributions of altered synapse properties and neural circuit activity to ASD pathogenesis.