The Endosome Localized Arf-GAP AGAP1 Modulates Dendritic Spine Morphology Downstream of the Neurodevelopmental Disorder Factor Dysbindin

AGAP duplicons associate with structural diversity at Chromosome 10q11.22

The 10q11.22 chromosomal region is a duplication-rich interval of the human genome and one of the last to be fully assembled. It carries copy number-variable genes associated with intellectual disability, bipolar disorder, and obesity. In this study, we characterized the structural diversity at this locus by analyzing 64 haploid assemblies produced by the Human Pangenome Reference Consortium. We identified 11 alternative haplotypes that differ in the copy number and/or orientation of large genomic segments, ranging from hundreds of kilobase pairs (kbp) to over one megabase pair (Mbp). We uncovered a 2.4 Mbp size difference between the shortest and longest haplotypes. Breakpoint analysis revealed that genomic instability results from nonallelic homologous recombination between segmental duplication (SD) pairs with varying similarity (94.4%-99.6%). Nonetheless, these pairs generally recombine at positions where their identity is higher (>99.6%). Recurrent inversions occur with different breakpoints within the same inverted SD pair. Inversion polymorphisms shuffle the entire SD arrangement, creating new predispositions to copy-number variations. The SD architecture is associated with a catarrhine-specific subgroup of the AGAP gene family, which likely triggered the accumulation of SDs at this locus over the past 25 million years of human evolution. Our results reveal extensive structural diversity and genomic instability at the 10q11.22 locus, and expand the general understanding of the mutational mechanisms behind SD-mediated rearrangements.

An engineered trafficking biosensor reveals a role for DNAJC13 in DOR downregulation

Trafficking of G protein-coupled receptors (GPCRs) through the endosomal-lysosomal pathway is critical to homeostatic regulation of GPCRs following activation with agonist. Identifying the genes involved in GPCR trafficking is challenging due to the complexity of sorting operations and the large number of cellular proteins involved in the process. Here, we developed a high-sensitivity biosensor for GPCR expression and agonist-induced trafficking to the lysosome by leveraging the ability of the engineered peroxidase APEX2 to activate the fluorogenic substrate Amplex UltraRed (AUR). We used the GPCR-APEX2/AUR assay to perform a genome-wide CRISPR interference screen focused on identifying genes regulating expression and trafficking of the δ-opioid receptor (DOR). We identified 492 genes consisting of both known and new regulators of DOR function. We demonstrate that one new regulator, DNAJC13, controls trafficking of multiple GPCRs, including DOR, through the endosomal-lysosomal pathway by regulating the composition of the endosomal proteome and endosomal homeostasis.

Social isolation-induced transcriptomic changes in mouse hippocampus impact the synapse and show convergence with human genetic risk for neurodevelopmental phenotypes

Early life stress (ELS) can impact brain development and is a risk factor for neurodevelopmental disorders such as schizophrenia. Post-weaning social isolation (SI) is used to model ELS in animals, using isolation stress to disrupt a normal developmental trajectory. We aimed to investigate how SI affects the expression of genes in mouse hippocampus and to investigate how these changes related to the genetic basis of neurodevelopmental phenotypes. BL/6J mice were exposed to post-weaning SI (PD21-25) or treated as group-housed controls (n = 7-8 per group). RNA sequencing was performed on tissue samples from the hippocampus of adult male and female mice. Four hundred and 1,215 differentially-expressed genes (DEGs) at a false discovery rate of < 0.05 were detected between SI and control samples for males and females respectively. DEGS for both males and females were significantly overrepresented in gene ontologies related to synaptic structure and function, especially the post-synapse. DEGs were enriched for common variant (SNP) heritability in humans that contributes to risk of neuropsychiatric disorders (schizophrenia, bipolar disorder) and to cognitive function. DEGs were also enriched for genes harbouring rare de novo variants that contribute to autism spectrum disorder and other developmental disorders. Finally, cell type analysis revealed populations of hippocampal astrocytes that were enriched for DEGs, indicating effects in these cell types as well as neurons. Overall, these data suggest a convergence between genes dysregulated by the SI stressor in the mouse and genes associated with neurodevelopmental disorders and cognitive phenotypes in humans.

AGAP1-associated endolysosomal trafficking abnormalities link gene-environment interactions in neurodevelopmental disorders

AGAP1 is an Arf1 GTPase-activating protein that regulates endolysosomal trafficking. Damaging variants have been linked to cerebral palsy and autism. We report three new cases in which individuals had microdeletion variants in AGAP1. The affected individuals had intellectual disability (3/3), autism (3/3), dystonia with axial hypotonia (1/3), abnormalities of brain maturation (1/3), growth impairment (2/3) and facial dysmorphism (2/3). We investigated mechanisms potentially underlying AGAP1 variant-mediated neurodevelopmental impairments using the Drosophila ortholog CenG1a. We discovered reduced axon terminal size, increased neuronal endosome abundance and elevated autophagy compared to those in controls. Given potential incomplete penetrance, we assessed gene-environment interactions. We found basal elevation in the phosphorylation of the integrated stress-response protein eIF2α (or eIF2A) and inability to further increase eIF2α phosphorylation with subsequent cytotoxic stressors. CenG1a-mutant flies had increased lethality from exposure to environmental insults. We propose a model wherein disruption of AGAP1 function impairs endolysosomal trafficking, chronically activating the integrated stress response and leaving AGAP1-deficient cells susceptible to a variety of second-hit cytotoxic stressors. This model may have broader applicability beyond AGAP1 in instances where both genetic and environmental insults co-occur in individuals with neurodevelopmental disorders.

A neurodevelopmental disorder associated with an activating de novo missense variant in ARF1

ADP-ribosylation factor 1 (ARF1) is a small GTPase that regulates membrane traffic at the Golgi apparatus and endosomes through recruitment of several coat proteins and lipid-modifying enzymes. Here, we report a pediatric patient with an ARF1-related disorder because of a monoallelic de novo missense variant (c.296 G > A; p.R99H) in the ARF1 gene, associated with developmental delay, hypotonia, intellectual disability and motor stereotypies. Neuroimaging revealed a hypoplastic corpus callosum and subcortical white matter abnormalities. Notably, this patient did not exhibit periventricular heterotopias previously observed in other patients with ARF1 variants (including p.R99H). Functional analysis of the R99H-ARF1 variant protein revealed that it was expressed at normal levels and properly localized to the Golgi apparatus; however, the expression of this variant caused swelling of the Golgi apparatus, increased the recruitment of coat proteins such as coat protein complex I, adaptor protein complex 1 and GGA3 and altered the morphology of recycling endosomes. In addition, we observed that the expression of R99H-ARF1 prevented dispersal of the Golgi apparatus by the ARF1-inhibitor brefeldin A. Finally, protein interaction analyses showed that R99H-ARF1 bound more tightly to the ARF1-effector GGA3 relative to wild-type ARF1. These properties were similar to those of the well-characterized constitutively active Q71L-ARF1 mutant, indicating that the pathogenetic mechanism of the R99H-ARF1 variant involves constitutive activation with resultant Golgi and endosomal alterations. The absence of periventricular nodular heterotopias in this R99H-ARF1 subject also indicates that this finding may not be a consistent phenotypic expression of all ARF1-related disorders.

AGAP1-associated endolysosomal trafficking abnormalities link gene-environment interactions in a neurodevelopmental disorder

AGAP1 is an Arf1 GAP that regulates endolysosomal trafficking. Damaging variants have been linked to cerebral palsy and autism. We report 3 new individuals with microdeletion variants in AGAP1 . Affected individuals have intellectual disability (3/3), autism (3/3), dystonia with axial hypotonia (1/3), abnormalities of brain maturation (1/3), growth impairment (2/3) and facial dysmorphism (2/3). We investigated mechanisms potentially underlying AGAP1 neurodevelopmental impairments using the Drosophila ortholog, CenG1a . We discovered reduced axon terminal size, increased neuronal endosome abundance, and elevated autophagy at baseline. Given potential incomplete penetrance, we assessed gene-environment interactions. We found basal elevation in phosphorylation of the integrated stress-response protein eIF2α and inability to further increase eIF2α-P with subsequent cytotoxic stressors. CenG1a -mutant flies have increased lethality from exposure to environmental insults. We propose a model wherein disruption of AGAP1 function impairs endolysosomal trafficking, chronically activating the integrated stress response, and leaving AGAP1-deficient cells susceptible to a variety of second hit cytotoxic stressors. This model may have broader applicability beyond AGAP1 in instances where both genetic and environmental insults co-occur in individuals with neurodevelopmental disorders.

Summary statement

We describe 3 additional patients with heterozygous AGAP1 deletion variants and use a loss of function Drosophila model to identify defects in synaptic morphology with increased endosomal sequestration, chronic autophagy induction, basal activation of eIF2α-P, and sensitivity to environmental stressors.

Differential Methylation Profile in Fragile X Syndrome-Prone Offspring Mice after in Utero Exposure to Lactobacillus Reuteri

Environmental factors such as diet, gut microbiota, and infections have proven to have a significant role in epigenetic modifications. It is known that epigenetic modifications may cause behavioral and neuronal changes observed in neurodevelopmental disabilities, including fragile X syndrome (FXS) and autism (ASD). Probiotics are live microorganisms that provide health benefits when consumed, and in some cases are shown to decrease the chance of developing neurological disorders. Here, we examined the epigenetic outcomes in offspring mice after feeding of a probiotic organism, Lactobacillus reuteri (L. reuteri), to pregnant mother animals. In this study, we tested a cohort of Western diet-fed descendant mice exhibiting a high frequency of behavioral features and lower FMRP protein expression similar to what is observed in FXS in humans (described in a companion manuscript in this same GENES special topic issue). By investigating 17,735 CpG sites spanning the whole mouse genome, we characterized the epigenetic profile in two cohorts of mice descended from mothers treated and non-treated with L. reuteri to determine the effect of prenatal probiotic exposure on the prevention of FXS-like symptoms. We found several genes involved in different neurological pathways being differentially methylated (p ≤ 0.05) between the cohorts. Among the key functions, synaptogenesis, neurogenesis, synaptic modulation, synaptic transmission, reelin signaling pathway, promotion of specification and maturation of neurons, and long-term potentiation were observed. The results of this study are relevant as they could lead to a better understanding of the pathways involved in these disorders, to novel therapeutics approaches, and to the identification of potential biomarkers for early detection of these conditions.

Crystal structure of the GTP-binding protein-like domain of AGAP1

AGAP1 is often considered to regulate membrane trafficking, protein transport and actin cytoskeleton dynamics. Recent studies have shown that aberrant expression of AGAP1 is associated with many diseases, including neurodevelopmental disorders and acute lymphoblastic leukemia. It has been proposed that the GTP-binding protein-like domain (GLD) is involved in the binding of cofactors and thus regulates the catalytic activity of AGAP1. To obtain a better understanding of the pathogenic mechanism underpinning AGAP1-related diseases, it is essential to obtain structural information. Here, the GLD (residues 70-235) of AGAP1 was overexpressed in Escherichia coli BL21 (DE3) cells. Affinity and gel-filtration chromatography were used to obtain AGAP1GLD with high purity for crystallization. Using the hanging-drop vapor-diffusion method with the protein at a final concentration of 20 mg ml-1, AGAP1GLD protein crystals of suitable size were obtained. The crystals were found to diffract to 3.0 Å resolution and belonged to space group I4, with unit-cell parameters a = 100.39, b = 100.39, c = 48.08 Å. The structure of AGAP1GLD exhibits the highly conserved functional G1-G5 loops and is generally similar to other characterized ADP-ribosylation factor (Arf) GTPase-activating proteins (GAPs), implying an analogous function to Arf GAPs. Additionally, this study indicates that AGAP1 could be classified as a type of NTPase, the activity of which might be regulated by protein partners or by its other domains. Taken together, these results provide insight into the regulatory mechanisms of AGAP1 in cell signaling.

The small GTPase ARF6 regulates GABAergic synapse development

ADP ribosylation factors (ARFs) are a family of small GTPases composed of six members (ARF1-6) that control various cellular functions, including membrane trafficking and actin cytoskeletal rearrangement, in eukaryotic cells. Among them, ARF1 and ARF6 are the most studied in neurons, particularly at glutamatergic synapses, but their roles at GABAergic synapses have not been investigated. Here, we show that a subset of ARF6 protein is localized at GABAergic synapses in cultured hippocampal neurons. In addition, we found that knockdown (KD) of ARF6, but not ARF1, triggered a reduction in the number of GABAergic synaptic puncta in mature cultured neurons in an ARF activity-dependent manner. ARF6 KD also reduced GABAergic synaptic density in the mouse hippocampal dentate gyrus (DG) region. Furthermore, ARF6 KD in the DG increased seizure susceptibility in an induced epilepsy model. Viewed together, our results suggest that modulating ARF6 and its regulators could be a therapeutic strategy against brain pathologies involving hippocampal network dysfunction, such as epilepsy.

Targeted resequencing identifies genes with recurrent variation in cerebral palsy

A growing body of evidence points to a considerable and heterogeneous genetic aetiology of cerebral palsy (CP). To identify recurrently variant CP genes, we designed a custom gene panel of 112 candidate genes. We tested 366 clinically unselected singleton cases with CP, including 271 cases not previously examined using next-generation sequencing technologies. Overall, 5.2% of the naïve cases (14/271) harboured a genetic variant of clinical significance in a known disease gene, with a further 4.8% of individuals (13/271) having a variant in a candidate gene classified as intolerant to variation. In the aggregate cohort of individuals from this study and our previous genomic investigations, six recurrently hit genes contributed at least 4% of disease burden to CP: COL4A1, TUBA1A, AGAP1, L1CAM, MAOB and KIF1A. Significance of Rare VAriants (SORVA) burden analysis identified four genes with a genome-wide significant burden of variants, AGAP1, ERLIN1, ZDHHC9 and PROC, of which we functionally assessed AGAP1 using a zebrafish model. Our investigations reinforce that CP is a heterogeneous neurodevelopmental disorder with known as well as novel genetic determinants.

Epigenome-wide Association Study of Attention-Deficit/Hyperactivity Disorder Symptoms in Adults

BACKGROUND

Previous studies have reported associations between attention-deficit/hyperactivity disorder symptoms and DNA methylation in children. We report the first epigenome-wide association study meta-analysis of adult attention-deficit/hyperactivity disorder symptoms, based on peripheral blood DNA methylation (Infinium HumanMethylation450K array) in three population-based adult cohorts.

METHODS

An epigenome-wide association study was performed in the Netherlands Twin Register (N = 2258, mean age 37 years), Dunedin Multidisciplinary Health and Development Study (N = 800, age 38 years), and Environmental Risk Longitudinal Twin Study (N = 1631, age 18 years), and results were combined through meta-analysis (total sample size N = 4689). Region-based analyses accounting for the correlation between nearby methylation sites were also performed.

RESULTS

One epigenome-wide significant differentially methylated position was detected in the Dunedin study, but meta-analysis did not detect differentially methylated positions that were robustly associated across cohorts. In region-based analyses, six significant differentially methylation regions (DMRs) were identified in the Netherlands Twin Register, 19 in the Dunedin study, and none in the Environmental Risk Longitudinal Twin Study. Of these DMRs, 92% were associated with methylation quantitative trait loci, and 68% showed moderate to large blood-brain correlations for DNA methylation levels. DMRs included six nonoverlapping DMRs (three in the Netherlands Twin Register, three in the Dunedin study) in the major histocompatibility complex, which were associated with expression of genes in the major histocompatibility complex, including C4A and C4B, previously implicated in schizophrenia.

CONCLUSIONS

Our findings point at new candidate loci involved in immune and neuronal functions that await further replication. Our work also illustrates the need for further research to examine to what extent epigenetic associations with psychiatric traits depend on characteristics such as age, comorbidities, exposures, and genetic background.

Genome-wide methylation analysis of a large population sample shows neurological pathways involvement in chronic widespread musculoskeletal pain

Chronic widespread musculoskeletal pain (CWP), has a considerable heritable component, which remains to be explained. Epigenetic factors may contribute to and account for some of the heritability estimate. We analysed epigenome-wide methylation using MeDIPseq in whole blood DNA from 1708 monozygotic and dizygotic Caucasian twins having CWP prevalence of 19.9%. Longitudinally stable methylation bins (lsBINs), were established by testing repeated measurements conducted ≥3 years apart, n = 292. DNA methylation variation at lsBINs was tested for association with CWP in a discovery set of 50 monozygotic twin pairs discordant for CWP, and in an independent dataset (n = 1608 twins), and the results from the 2 samples were combined using Fisher method. Functional interpretation of the most associated signals was based on functional genomic annotations, gene ontology, and pathway analyses. Of 723,029 signals identified as lsBINs, 26,399 lsBINs demonstrated the same direction of association in both discovery and replication datasets at nominal significance (P ≤ 0.05). In the combined analysis across 1708 individuals, whereas no lsBINs showed genome-wide significance (P < 10-8), 24 signals reached p≤9E-5, and these included association signals mapping in or near to IL17A, ADIPOR2, and TNFRSF13B. Bioinformatics analyses of the associated methylation bins showed enrichment for neurological pathways in CWP. We estimate that the variance explained by epigenetic factors in CWP is 6%. This, the largest study to date of DNA methylation in CWP, points towards epigenetic modification of neurological pathways in CWP and provides proof of principle of this method in teasing apart the complex risk factors for CWP.

The importance of data structure in statistical analysis of dendritic spine morphology

BACKGROUND

Dendritic spine morphology is heterogeneous and highly dynamic. To study the changing or aberrant morphology in test setups, often spines from several neurons from a few experimental units e.g. mice or primary neuronal cultures are measured. This strategy results in a multilevel data structure, which, when not properly addressed, has a high risk of producing false positive and false negative findings.

METHODS

We used mixed-effects models to deal with data with a multilevel data structure and compared this method to analyses at each level. We apply these statistical tests to a dataset of dendritic spine morphology parameters to illustrate advantages of multilevel mixed-effects model, and disadvantages of other models.

RESULTS

We present an application of mixed-effects models for analyzing dendritic spine morphology datasets while correcting for the data structure.

COMPARISON WITH EXISTING METHODS

We further show that analyses at spine level and aggregated levels do not adequately account for the data structure, and that they may lead to erroneous results.

CONCLUSION

We highlight the importance of data structure in dendritic spine morphology analyses and highly recommend the use of mixed-effects models or other appropriate statistical methods to deal with multilevel datasets. Mixed-effects models are easy to use and superior to commonly used methods by including the data structure and the addition of other explanatory variables, for example sex, and age, etc., as well as interactions between variables or between variables and level identifiers.

Next Generation Proteomic Pipeline for Chromosome-Based Proteomic Research Using NeXtProt and GENCODE Databases

Human Proteome Project aims to map all human proteins including missing proteins as well as proteoforms with post translational modifications, alternative splicing variants (ASVs), and single amino acid variants (SAAVs). neXtProt and Ensemble databases are usually used to provide curated information on human coding genes. However, to find these proteoforms, we (Chr #11 team) first introduce a streamlined pipeline using customized and concatenated neXtProt and GENCODE originated from Ensemble, with controlled false discovery rate (FDR). Because of large sized databases used in this pipeline, we found more stringent FDR filtering (0.1% at the peptide level and 1% at the protein level) to claim novel findings, such as GENCODE ASVs and missing proteins, from human hippocampus data set (MSV000081385) and ProteomeXchange (PXD007166). Using our next generation proteomic pipeline (nextPP) with neXtProt and GENCODE databases, two missing proteins such as activity-regulated cytoskeleton-associated protein (ARC, Chr 8) and glutamate receptor ionotropic, kainite 5 (GRIK5, Chr 19) were additionally identified with two or more unique peptides from human brain tissues. Additionally, by applying the pipeline to human brain related data sets such as cortex (PXD000067 and PXD000561), spinal cord, and fetal brain (PXD000561), seven GENCODE ASVs such as ACTN4-012 (Chr.19), DPYSL2-005 (Chr.8), MPRIP-003 (Chr.17), NCAM1-013 (Chr.11), EPB41L1-017 (Chr.20), AGAP1-004 (Chr.2), and CPNE5-005 (Chr.6) were identified from two or more data sets. The identified peptides of GENCODE ASVs were mapped onto novel exon insertions, alternative translations at 5'-untranslated region, or novel protein coding sequence. Applying the pipeline to male reproductive organ related data sets, 52 GENCODE ASVs were identified from two testis (PXD000561 and PXD002179) and a spermatozoa (PXD003947) data sets. Four out of 52 GENCODE ASVs such as RAB11FIP5-008 (Chr. 2), RP13-347D8.7-001 (Chr. X), PRDX4-002 (Chr. X), and RP11-666A8.13-001 (Chr. 17) were identified in all of the three samples.