Next Article in Journal
Collagen Fiber Array of Peritumoral Stroma Influences Epithelial-to-Mesenchymal Transition and Invasive Potential of Mammary Cancer Cells
Previous Article in Journal
Discrepancies in Persistent Dry Eye Signs and Symptoms in Bilateral Pseudophakic Patients
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
JCM
Volume 8
Issue 2
10.3390/jcm8020212
jcm-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of a Sardinian Multiplex Family with Autism Spectrum Disorder Points to Post-Synaptic Density Gene Variants and Identifies CAPG as a Functionally Relevant Candidate Gene

by
Elena Bacchelli
1,*,†,
Eleonora Loi
2,†,
Cinzia Cameli
1,†,
Loredana Moi
2,
Ana Florencia Vega Benedetti
2,
Sylvain Blois
2,
Antonio Fadda
2,
Elena Bonora
3,
Sandra Mattu
4,
Roberta Fadda
5,
Rita Chessa
6,
Elena Maestrini
1,
Giuseppe Doneddu
6,† and
Patrizia Zavattari
2,*,†
1
Department of Pharmacy and Biotechnology, University of Bologna, 40126 Bologna, Italy
2
Department of Biomedical Sciences, Unit of Biology and Genetics, University of Cagliari, 09042 Cagliari, Italy
3
Department of Medical and Surgical Sciences, DIMEC, St. Orsola-Malpighi Hospital, University of Bologna, 40138 Bologna, Italy
4
Department of Biomedical Sciences, Unit of Oncology and Molecular Pathology, University of Cagliari, 09124 Cagliari, Italy
5
Department of Pedagogy, Psychology, Philosophy, University of Cagliari, 09123 Cagliari, Italy
6
Center for Pervasive Developmental Disorders, AO Brotzu, 09134 Cagliari, Italy
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
J. Clin. Med. 2019, 8(2), 212; https://doi.org/10.3390/jcm8020212
Submission received: 15 January 2019 / Revised: 4 February 2019 / Accepted: 5 February 2019 / Published: 7 February 2019
(This article belongs to the Section Clinical Laboratory Medicine)
Browse Figures
Versions Notes

Abstract

:
Autism spectrum disorders (ASDs) are a group of neurodevelopmental disorders with high heritability, although their underlying genetic factors are still largely unknown. Here we present a comprehensive genetic characterization of two ASD siblings from Sardinia by genome-wide copy number variation analysis and whole exome sequencing (WES), to identify novel genetic alterations associated with this disorder. Single nucleotide polymorphism (SNP) array data revealed a rare microdeletion involving CAPG, ELMOD3, and SH2D6 genes, in both siblings. CAPG encodes for a postsynaptic density (PSD) protein known to regulate spine morphogenesis and synaptic formation. The reduced CAPG mRNA and protein expression levels in ASD patients, in the presence of hemizygosity or a particular genetic and/or epigenetic background, highlighted the functional relevance of CAPG as a candidate gene for ASD. WES analysis led to the identification in both affected siblings of a rare frameshift mutation in VDAC3, a gene intolerant to loss of function mutation, encoding for a voltage-dependent anion channel localized on PSD. Moreover, four missense damaging variants were identified in genes intolerant to loss of function variation encoding for PSD proteins: PLXNA2, KCTD16, ARHGAP21, and SLC4A1. This study identifies CAPG and VDAC3 as candidate genes and provides additional support for genes encoding PSD proteins in ASD susceptibility.

Graphical Abstract

1. Introduction

Autism spectrum disorder (ASD) (Mendelian Inheritance in Man (MIM) 209850) is a neurodevelopmental disorder characterized by early onset in childhood, atypical and repetitive behaviors and impairments in communication and reciprocal social interaction. The phenotype of ASDs is extremely heterogeneous, with individual differences in a wide range of symptoms and severity [1,2,3,4].
Moreover, comorbidity with other neurological and behavioral disorders, such as developmental delay/intellectual disability, obsessive–compulsive disorder, attention deficit hyperactivity disorder, depression, anxiety disorder, social phobia, bipolar disorder, seizures, and others, is frequently observed.
The estimated prevalence of ASD in the worldwide population is 1 in 160 children and appears to be increasing [5], making it a major burden to society. Autism is more prevalent in male individuals, with an overall male-to-female ratio close to 3:1 [6].
Multiple lines of evidence support the strong role of genetics and heritability in the etiology of ASDs [2,3,4]. However, in line with the highly heterogeneous phenotypic manifestations, the underlying genetic architecture of ASD is also complex, with a combined contribution of both common and rare forms of genetic variations [7].
The major advances in the field of molecular diagnostic technology, such as the development of chromosome microarray technology (CMA) and whole-exome sequencing (WES), have allowed to identify more than 800 genes and loci, harboring rare de novo and inherited copy number variants (CNVs) and single nucleotide variants (SNVs), predisposing to ASDs. Even if the contribution of each gene/locus to ASD risk is very low, most of them converge on a smaller set of specific pathways such as synaptic function, transcriptional regulation and chromatin remodeling. In particular, alteration of synaptic homeostasis may be a common biological process associated with ASDs [8]. Indeed, many genes involved in synaptic function, and notably postsynaptic density (PSD) genes, a huge protein complex crucial for synaptic transmission and plasticity, were shown to play a role in ASD and other neurodevelopmental disorders [9].
Despite the major advances that have been made in recent years in our understanding of the genetic architecture of ASD, causative mutations remain unknown in the majority of ASD patients.
The Sardinian population has often proved to be fundamental for genetic studies of association in monogenic but also multifactorial diseases. The complex multi-locus and multi-allelic nature of certain diseases has indeed been unveiled studying Sardinian cohorts [10,11,12,13,14], other times showing founder effects only detectable in genetic isolates [15,16].
In the last years, CNV analysis by microarray technology has become the first-tiered genetic investigation for ASD, in particular when the affected individual shows also intellectual disability. However, CNV analysis can explain only a portion of ASD genetic architecture, given that the ASD phenotype may be the result of multiple variants of different size and effect, possibly affecting different genes involved in the same biological pathway. Therefore, even in presence of relevant CNVs, combining CNV analysis with sequencing data is warranted in order to better understand the role of different variants in the ASD phenotype [17].
Here we present a comprehensive genetic characterization of two ASD sibling from Sardinia. We identified a rare 2p11.2 deletion involving the whole coding region of CAPG, a member of postsynaptic density genes, and part of ELMOD3 and SH2D6 genes. We evaluated the impact of the deletion on CAPG expression. Finally, we identified deleterious variants in genes intolerant to loss of function mutations that could have a role in ASD susceptibility in combination with the 2p11.2 deletion.

2. Material and Methods

2.1. Participants

The Sardinian family participating in this study consists of father (AUT003.1), mother (AUT003.2), and two affected brothers (AUT003.3 and AUT003.4). Genomic DNA obtained from blood was available for all members of the pedigree.
Both siblings were diagnosed using the Autism Diagnostic Observation Schedule—Generic (ADOS) protocols [18] plus clinical evaluation, according to the VI Diagnostic and Statistical Manual of Mental Disorders (DSM-IV) criteria.
The older brother, named AUT003.3, received a diagnosis of Pervasive Developmental Disorder—Not otherwise Specified (PDD—NOS) when he was 7 years old, according to the diagnostic criteria of the DSM-IV. The results of the module II of the ADOS showed a deficit in the reciprocal social interaction with restricted and repetitive interests, difficulty in the changes of the routine of daily life, significant global motor clumsiness and deficit of eye contact. At that time, his non-verbal intellectual functioning was in the norm (Leiter-r Intelligence Quotient, IQ = 106).
The younger brother, named AUT003.4, received a diagnosis of Autism Disorder when he was 4 years old (DSM-IV). The results of the module I of the ADOS showed impairment in communication and language refers to a framework of autism, and in reciprocal social interaction, corresponding to a spectrum disorder. He also presented restricted and repetitive interests and stereotypies. At that time, his non-verbal intellectual functioning was in the norm (Leiter-r Intelligence Quotient, IQ = 87). He also showed a significant delay in receptive and expressive language, with aspects of immediate echolalia and poor spontaneous verbalization. At the age of 10, he showed a significant cognitive delay (IQ = 46). The ADOS at the age of 10 indicates difficulties in the areas of communication and reciprocal social interaction, as well as restricted and repetitive interests and stereotypies.
The Social Responsiveness Scale: Adult Research Version (SRS: ARV) and the Broad Autism Phenotype Questionnaire (BAPQ) were used in order to assess the parents’ Broader Autism Phenotype (BAP) (Table 1). Both parents did not exceed clinical cutoffs in BAP total measures (BAPQ Total and SRS: ARV Total). Only the father exceeded BAPQ score cutoffs for the “Social Behavior” (Aloof) subscale. In general, the father shows higher scores than the mother in all the measures of the BAPQ and the SRS: ARV. In addition, to exclude a possible diagnosis of ASD, parents were also assessed with the Autism Diagnostic Observation Schedule (ADOS—module IV), that resulted to be negative for both parents.
Simultaneously with the AUT003 family, blood samples were taken from two unaffected children of the same age of probands AUT003.3 and AUT003.4, to be used as experimental controls in the expression studies and in the sequencing of the CAPG promoter.
Subsequently, blood samples of 13 ASD probands and 22 healthy control subjects were collected to test the expression levels of CAPG mRNA. The 13 ASD probands were recruited at the Stella Maris Clinical Research Institute for Child and Adolescent Neuropsychiatry (Calambrone, Pisa, Italy). ASD diagnosis was based on the Autism Diagnostic Interview-Revised (ADI-R) and ADOS. A clinical evaluation was undertaken in order to exclude known syndromes associated with autism. Standard karyotyping, fragile-X testing, electroencephalogram (EEG), and array-based comparative genomic hybridization (aCGH) were obtained for all probands. The control sample consists of 22 unrelated Italian individuals with no psychiatric disorders.
This study was approved by the local Ethical Committee and took place in observation of the declaration of Helsinki.

2.2. DNA and RNA Samples Extraction

Genomic DNA was extracted from peripheral whole blood lymphocytes using salting-out procedure. DNA was quantified with NanoDrop (NanoDrop Products, Thermo Scientific, Wilmington, DE, USA) and by fluorometric reading (Quant-iT™ PicoGreen® dsDNA Assay Kit, Thermo Fisher Scientific, Waltham, MA, USA).
Total RNA from all members of family AUT003, 13 ASD probands, and 24 healthy control subjects was extracted from PBMCs (Peripheral Blood Mononuclear Cells, isolated over a Ficoll-Hypaque density gradient) using Qiagen RNAeasy mini kit (Qiagen, Hilden, Germany) and quantified with NanoDrop. cDNA was synthesized using the High Capacity Kit (Applied Biosystems, Carlsbad, CA, USA). Then, 40 ng of cDNA were used to test CAPG and ELMOD3 gene expression by quantitative reverse transcription PCR (qRT-PCR) using SYBR Green in all individuals.

2.3. CNV Analysis

Genotyping was performed using Human1M-Duo DNA Analysis BeadChip® (Illumina, San Diego, CA, USA), which interrogate 1.2 million loci, according to the manufacturer’s instructions.
The clustering algorithm implemented in GenomeStudio (V2010.3) was used to cluster the data. GenomeStudio was also used to evaluate all genotypes using a quantitative genotype quality score called GenCall (GC) score, ranging from 0 to 1 with 1 being the best. The quality score cutoff was set at 0.15. Single nucleotide polymorphism (SNP) quality control (QC) was performed according to recommended guidelines [19].
Genome-wide CNVs were detected using a multiple-algorithm approach to maximize sensitivity and specificity of CNV calling. CNVs required calling by at least two different CNV detection algorithms out of these three: PennCNV [20], QuantiSNP [21], and CNVPartition (Illumina, San Diego, CA, USA). The required data for CNV analysis, (i.e., between-sample normalized fluorescence (Log R ratios (LRR) and B allele frequency (BAF) values) and genotypes for each sample), were exported directly from Illumina’s GenomeStudio software. To reduce false positive calls, we excluded CNVs if failed stringent quality-control criteria: spanning <10 probes, low confidence score (confidence score <15), more than 50% of their length overlapping segmental duplications and if they were within centromere proximal cytobands. As a final step, we joined CNVs that appeared to be artificially split by either of the calling algorithms.
The physical coordinates of each CNV, base pair (bp) position along the chromosome, were based on Build 37 (hg19) of the human genome assembly. CNVs passing QC were finally annotated to RefSeq genes, within 10 kb beyond the 5′- and 3′-untranslated regions (UTRs), to include CNVs extending to neighboring genes or overlapping potential regulatory regions. Similarly, we annotated CNVs overlapping exons.
We further filtered these stringent genic CNVs to identify a subset of “rare” CNVs, defined as CNV calls overlapping (>75%) copy number stable regions of the genome [22].
We manually inspected all rare CNVs by plotting their Log R ratios (LRR) and B allele frequencies (BAF).

2.4. CNV Characterisation and Segregation Analysis

In order to validate rare interesting CNVs in family AUT003, quantitative PCR experiments were performed using Fast SYBR Green (Thermo Fisher Scientific, Waltham, MA, USA). Each assay was conducted in triplicate, with at least three sets of primers corresponding to the region of interest and another mapping to a control region on FOXP2 gene at 7q31.1. The number of copies of each amplified fragment was calculated using the ΔΔCt method. All family members (the parents and both siblings) were tested for confirming inheritance and segregation of CNVs.
To help resolve the microdeletion breakpoints, long range PCR across the 2p11.2 deletion was then carried out using the KAPA HiFi HotStart DNA Polymerase following the manufacturer’s suggested protocol (primers GGTGAAAATCTCCATTCCTTGA and CGACCCAGGCTTTAGTTATCAG).
The 241-bp deletion-spanning amplicon was purified and sequenced using Sanger method.

2.5. Gene Expression Analysis

CAPG and ELMOD3 transcript levels were tested in family AUT003, in a commercially available cDNA panel (Human MTC Panel I, Clontech, Mountain View, CA, USA) including five different tissues (brain, testis, kidney, colon, prostate, thyroid) and in cDNA derived from blood of all member of family AUT003 and two unaffected controls by qRT-PCR with SYBR Green using 40 ng of cDNA. CAPG transcript levels were also assessed in 13 ASD probands and 22 additional healthy control subjects by qRT-PCR. To amplify CAPG fragment, we designed the forward primer in exon 4 and reverse primer in exon 5 (primers: CAGGGCAATGAGTCTGACCT and GGCTCCTGTGGAGGTCTTGT; 102 bp amplicon length; NM_001256140.1); for ELMOD3 fragment, we designed forward primer in exon 5 and reverse primer in exon 6 (primers: GGTGAAAGATTGTCTGCTGGA and AGAATGCCATGGTTCTTCAA; 104 bp amplicon length; NM_001135022.1). All the data were normalized using the housekeeping gene TFRC as reference gene.

2.6. Western Blot

Proteins from lymphocytes were isolated using the Membrane Protein Extraction Kit (Mem-PERTM Plus, Thermo Fisher Scientific, Waltham, MA, USA) following the manufacturer’s instructions. Protein extracts were fractioned by 10% Sodium dodecyl sulfate (SDS) polyacrylamide gel and separated at 100 V during 90 min. Then, the proteins were transferred onto 0.45 µm nitrocellulose membrane, blocked in 5% milk overnight at 4 °C and then incubated 2 h with primary antibody anti-CAPG (1:1000, 10194-1-AP, Proteintech, Rosemont, IL, USA). After hybridization with primary antibody, the membrane was washed with Tris-buffered saline containing Tween-20 three times, then incubated with horseradish peroxidase-labeled secondary antibody (Jackson ImmunoResearch, Cambridge House, UK) for 1 h at room temperature and washed with Tris-buffered saline containing Tween-20 three times. Final detection was performed with Luminol-based Enhanced Chemiluminescent (ECL) Western blotting reagents (BIO-RAD, Hercules, CA, USA). An antibody against GAPDH (1:500, MAB374, Millipore, Burlington, MA, USA) was used for gel-loading control. Western blot signals were quantified using ImageJ program. The intensity of each band was normalized respect to that of GAPDH. The assay was performed in duplicate.

2.7. Whole Exome Analysis

DNA from the two affected siblings was captured using the Nextera exome enrichment kit (Illumina Inc., San Diego, CA, USA) according to the manufacturer’s protocol and sequenced as 100-bp paired-end reads on the Illumina HiSeq2500 platform (Illumina Inc., San Diego, CA, USA).
Exomes had a read depth (DP) of 10× or more for 90% of the total exome coverage and 20× or more for 80%. The raw data generated by WES were analyzed using CoVaCS, a recently developed pipeline for Next Generation Sequencing data analysis [23]. CoVaCS generates a final set of variants using a consensus call-set approach from three different algorithm (GATK, Varscan and Freebayes).
VCF (Variant Call Format) files of the CoVaCS consensus calls were annotated using Annovar [24]. The gene based annotation (position, nomenclature, gene name, gene function) was performed using RefSeq.
In order to remove low-quality variants called genotypes were required to have DP ≥10, and Genome Quality (GQ) ≥20. We selected only rare variants with MAF ≤1% in gnomeAD exome, gnomeAD genome [25] and the 1000 Genomes Project [26]. We considered in our analysis only variants located at exons or splicing sites; furthermore, we selected only Likely Gene Disrupting (LGD) and damaging missense mutations. LGD mutations included stop-gain, stop-loss, frameshift and splicing mutations.
Damaging missense mutations were defined as deleterious by at least two of the following criteria: SIFT (Sorting Intolerant From Tolerant) score ≤0.05, Polyphen2 (HDIV) score ≥0.95, Mutation Assessor ≥2, Phred transformed CADD (Combined Annotation-Dependent Depletion) score ≥15, placental mammal PhyloP ≥2.4, vertebrate PhyloP ≥4 [27].
Synaptic genes and PSD genes were defined as previously described [28,29].

3. Results

3.1. Identification and Characterisation of CAPG Deletion in Family AUT003

During a genome-wide CNV scan using a multi-algorithm approach on one multiplex family with ASD from Sardinia, we identified a rare genic heterozygous deletion of ~37 kb, involving the whole CAPG gene, the last coding exons of ELMOD3 and the first non-coding exon of SH2D6.
We validated the microdeletion, confirming that it was inherited from the unaffected mother to both the affected siblings (Figure 1A). Analysis of chromosome 2 SNPs (mapping in a region of 138 kb including the deletion) indicates that the non-deleted paternal copy of CAPG was also shared identical-by-descent in all two affected children.
To determine whether this rare CNV may have been mediated by repetitive elements, we mapped the breakpoints at base pair resolution. At the deletion junction we identified the presence of 21 bp mapping both in chr2:85,610,292-85,610,313 and in chr2:85,647,268-85,647,289 inside two Alu repeats with high homology (AluSx1, chr2:85610243-85610455 and AluSx, chr2:85647216-85647540) (NCBI, National Center for Biotechnology Information, build 37; Figure 1B), suggesting that this deletion is likely to be generated through non-allelic homologous recombination. We resolved the size of the deletion of 36,977 bp.
CNV analysis of this family allowed us also to identify two additional rare exonic CNVs: a chr3q26.32 deletion including exon 1 of KCNMB2, transmitted from the unaffected mother to affected son AUT003.3, and a chr6q21 duplication involving the first four exons of BEND3 and the last four exons of PDSS2, inherited from the unaffected father by affected son AUT003.3 (Supplementary Figure S1).

3.2. ELMOD3-CAPG-SH2D6 Deletion Frequency in ASD Cases and Controls

In order to assess the relevance and test the frequency of ELMOD3-CAPG-SH2D6 in ASD susceptibility in autism and control populations, we exploited large-cohort whole-genome CNV data already available. We identified two rare deletions encompassing ELMOD3-CAPG-SH2D6 in two simplex ASD families from the AGP genome-wide study including CNV data of 2147 European ASD families [30], and in both cases the deletion was inherited from the father (Sample_ID: 5424_3 and 5536_3) [30]. Recently, a homozygous deletion affecting the last three exons of ELMOD3, the entire CAPG gene, and the first non-coding exon of SH2D6 has been identified in a Tunisian ASD proband with ID and hearing impairment. The deletion is inherited from consanguineous parents, both heterozygous for the deletion [31]. No CNVs in this genomic region are present in 6639 European controls [32,33,34,35,36]. Indeed, this genomic region is considered a copy number stable region, according to the CNV map of the human genome [22].
To exclude the hypothesis that the identified deletion might be relatively common in the Sardinian population, we analyzed 80 Sardinian control samples finding that no control sample showed the deletion.

3.3. CAPG and ELMOD3 Gene Expression Analysis

In order to obtain a global overview of the CAPG and ELMOD3 expression, we evaluated their expression levels in different human tissues and two unaffected controls. We observed the presence of CAPG and ELMOD3 transcripts in all tested tissues, including brain (Figure 2).
We did not test SH2D6 expression since the deletion involves only SH2D6 non-coding exon 1. Moreover, from the GTExPortal gene expression data, SH2D6 is reported to be expressed at very low levels in brain and its expression is barely detectable in whole blood. Then we tested CAPG and ELMOD3 expression levels in all members of the AUT003 family and in two healthy controls. CAPG expression was significantly reduced (p <0.001) in the mother and in the two ASD siblings carrying the heterozygous deletion, compared to two healthy controls who do not carry the deletion (Figure 3A). Conversely, we did not observe a reduction of ELMOD3 levels in the two affected sibling and their mother (Figure 3B) but the expression levels were highly variable between the tested individuals.
We also performed a quantitative characterization of CAPG expression in an independent cohort of 13 Italian ASD probands and 22 unaffected individuals. We observed a significant reduction of CAPG expression in ASD probands compared with healthy control subjects (Mann–Whitney test p = 0.02647) (Figure 4). Since all ASD probands were between 4 and 18 years of age while in our control samples there were also 12 older subjects (age range: 21–47 years), we repeated the analysis only considering the 10 controls aged 8 months–12 years. The difference between ASD and control children was still significant (Mann–Whitney test p = 0.03493), excluding a possible age effect influencing CAPG expression.

3.4. CAPG Protein Expression Analysis

The expression protein level of CAPG was determined in the two siblings affected by ASD, their unaffected mother and in two healthy control children. The Western blot analysis showed significant lower CAPG expression levels in the AUT003 family members with the deletion compared to healthy control subjects (Figure 5).

3.5. Whole Exome Sequencing

We performed a WES analysis for AUT003.3 and AUT003.4 by Illumina Nextera technology.
The 81.3–87.5% and 94.6–98.7% of target regions showed at least 20× and 10× coverage for AUT003.3 and AUT003.4, respectively.
In the first step of WES analysis, we focused on LGD variants given their known deleterious effect on protein function (Table 2). In particular, we restricted the analysis on variants in genes intolerant to loss-of-function (LoF) variants based on the pLi score, i.e. probability of being loss-of-function (LoF) intolerant (pLi) [37]. Among the 33 identified LGD variants, three are in genes strongly intolerant to loss of function variants (pLi score ≥0.9). Two of these variants, located in SIK1 and VDAC3, are shared by both affected siblings. The variant in SIK1 disrupts a splice site and it is inherited from the unaffected mother, while the variant in VDAC3 is a frameshift deletion transmitted by the unaffected father. The effect of the VDAC3 2 bp-frameshift deletion is the introduction of five novel amino acids followed by premature termination codon probably inducing nonsense-mediated decay. According to GTExPortal gene expression data, SIK1 is expressed at low levels in brain, while VDAC3 is highly expressed in brain and, interestingly, it encodes for a postsynaptic density member, CAPG.
To comprehensively inspect the load of potentially damaging variant shared by both affected siblings, we also investigated the presence of damaging missense variants in genes with a pLi score ≥0.9. Among the 20 shared variants identified in genes intolerant to mutation, four affect PSD genes: PLXNA2, KCTD16, ARHGAP21, and SLC4A1 (Table 3).
No damaging variants were identified in the hemizygous CAPG genomic region, thereby excluding a complete CAPG knock-out in both affected siblings.

4. Discussion

Genome-wide CNV analysis in a Sardinian multiplex family led to the identification of a rare genic heterozygous deletion on chromosome 2, involving the whole CAPG gene, the last coding exons of ELMOD3 and the first non-coding exon of SH2D6, inherited by both affected siblings from the unaffected mother. The deletion encompasses a CNV stable region of genome [22] and has not been reported in controls, while similar deletions spanning CAPG-ELMOD3-SH2D6 have been detected in three independent ASD families, in the heterozygous [30] and the homozygous [31] state, suggesting that this deletion might be involved in ASD susceptibility.
Interestingly, CAPG expression analysis in AUT003 family showed a significant reduction of mRNA (p <0.001) and protein expression (p <0.05), due to hemizygosity, in all carriers of the identified deletion (mother and both affected siblings). Moreover, we observed a significant reduction of CAPG mRNA expression levels in 13 ASD Italian probands compared to 22 unaffected controls (p = 0.026), suggesting that decreased CAPG levels might represent a risk factor with a low effect size for ASD. We could speculate that differences in CAPG gene expression could likely be due to the presence of different genetic and/or epigenetic background in ASD subjects compared to controls.
CAPG, capping actin protein gelsolin like, is an excellent candidate for neurologic development abnormalities, given that it is essential for dendritic spine development and synapse formation [38]. Abnormality in spine number and shape has been observed in a number of neurological disorders, including autism, and contributes to brain dysfunction [39]. Recently, it has been showed that CAPG is present in dendritic spines of cultured hippocampal neurons and the branched actin filament network containing CAPG appears to be a prominent feature of the spine head [40]. Capg knockdown in rat hippocampal neurons cultures resulted in a marked decline in spine density accompanied by increased filopodia-like protrusions and consequently a reduced number of functional synapses [38]. Moreover, through a computational approach to improve the reconstruction of the PSD network it has been suggested that CAPG interacts with SRC at the postsynaptic density [41]. Overall, these findings indicate that capping of actin filaments by CAPG represents an essential step for the remodeling of the actin architecture underlying spine morphogenesis and synaptic formation during development.
Interestingly, CAPG has also been implicated in Rett syndrome, a genetic neurodevelopmental disorder caused by variants in MECP2. In particular, transcriptomic data analysis obtained from human cells with impaired MECP2 function suggests that MECP2 regulates CAPG expression and that CAPG could be involved in altered cytoskeleton organization observed in Rett syndrome [42]. Moreover, Mecp2-null mice showed Capg upregulation that might be involved in the molecular mechanisms leading to altered neuronal architecture observed in Rett syndrome [43].
Given the presence of CAPG deletion also in the unaffected mother, we can hypothesize that only a recessive loss of function is sufficient to cross the risk threshold to develop ASD, while the presence of hemizygous CAPG, or even its reduced levels of expression, could confer a smaller increase in susceptibility and lead to the ASD only in association with other genetic risk variants.
Therefore, a WES analysis has been performed for both affected siblings, in order to search for additional genetic alterations that could act in combination with CAPG deletion to determine ASD phenotype. Among all rare variants identified by WES analysis, we focused first on LGD variants given their harmful effect on protein function. Moreover, since it has been reported that genes have a different tolerance for deleterious variants, we considered only LGD variants in genes reported to be intolerant to loss-of-function (LoF) variants based on pLi score [37]. The most interesting variant is a frameshift deletion in VDAC3 gene, shared by the two siblings and inherited from the unaffected father. The VDAC3 2 bp-frameshift deletion determines the introduction of five novel amino acids followed by premature termination codon but also should lead to mRNA degradation owing to nonsense-mediated mRNA decay. VDAC3 encodes the voltage-dependent anion channel 3 and it is a member of human postsynaptic density (hPSD) genes highly expressed in several brain regions including the hypothalamus and caudate nucleus. VDACs are a family of proteins present on the mitochondrial outer membrane, on which they supply diffusion pores for small hydrophilic molecules and are involved in the release of cytochrome c at the onset of apoptotic cell death. Interestingly, recent evidence suggests that VDAC could be involved in the pathogenesis of several syndromes, including autism [44,45,46]. In particular, Gonzalez-Gronow and colleagues identified the presence of autoantibodies against VDAC proteins in autistic individuals, suggesting their possible causal role in the neurologic pathogenesis of autism [44]. Moreover, Crane and colleagues have hypothesized that the beneficial effect recorded in patients with autism following treatment with coenzyme Q, or other agents influencing the transport of electrons, may be due to the control of such molecules on the porin channels [45]. Therefore, the VDAC3 frameshift mutation identified in this study might contribute to the ASD phenotype in the two siblings.
Recent studies suggested that the autism phenotype may be influenced by the combined effects of variants in multiple genes involved in common pathways. Thus, it is worthy of note that the two affected siblings share disrupting mutations (a microdeletion and a frameshift mutation) in two different genes, both encoding for proteins associated to the PSD, a huge synaptic protein complex linked to various neuropsychiatric diseases including neurodevelopmental disorders such as ASD [29].
Further evidence for a possible involvement of PSD genes in modulating the ASD phenotype in this family is provided by the identification of damaging missense mutations shared by the two affected siblings in four additional LoF intolerant genes (pLi score ≥0.9) encoding for PSD proteins: PLXNA2, KCTD16, ARHGAP21, and SLC4A1.
Of particular interest is the damaging mutation involving PLXNA2, a gene specifically expressed in fetal and adult brain [47]. This gene encodes the plexin A2, a member of the plexin-A family of semaphorin co-receptors, a large family of secreted or membrane-bound proteins that mediate repulsive effects on axon pathfinding during nervous system development. Interestingly, reduced PLXNA4 expression, another plexin-A family member, have been observed in individual with ASD [48], while a genetic variant in PLXNA2 has been associated with schizophrenia [49]. Moreover, it has been suggested that PlxnA2 is an inhibitor of the GTPase Rap1B-regulating granule cell (GC) progenitors in the hippocampal neuroepithelium in a Rap1-dependent manner [50]. Interestingly, similar to PlxnA2 GAP activity, Shank3, one of the stronger ASD candidate genes, inhibits Rap1 GTPases. In particular, two different SHANK3 mutations associated to ASD led to Rap1 activity reduction [51]. In addition, Plxna2−/− mice showed deficits in associative learning, sociability, and sensorimotor gating in behavioral studies [50]. Altogether, this evidence suggests that variants in PLXNA2 contribute to neuropsychiatric disorders.
The other three shared variants in synaptic genes were identified in KCTD16, ARHGAP21, and SLC4A1. KCTD16 encodes for an auxiliary subunit of GABA-B receptors that modulate both pre and postsynaptic GABA-B receptor response [52,53,54,55]. Interestingly, a KCTD16 polymorphism has been associated with epilepsy by a genetic linkage study [56].
SLC4A1 encodes for a member of the ion transporter SLC4 family, that mediate the pH regulation and transepithelial ion transport, critical for several biological processes, including modulation of neuronal excitability [57,58,59,60]. Dysfunctions of SLC4 transporters are associated with a broad range of diseases, including central nervous system disorders, such as intellectual disability, migraine, epilepsy, autism and drug addiction [61,62,63,64,65,66].
Finally, ARHGAP21 is a widely expressed gene with high levels in brain and muscle [67]. It encodes for a Rho GTPase activating-protein, localized to the Golgi complex through binding to GTP-bound ARF1. Arhgap21 protein acts as a Rho-GAP for Cdc42 and regulates the Arp2/3 complex and actin cytoskeleton remodeling, playing a role in Golgi complex organization [68]. Interestingly, CAPG plays a central role in the assembly of the actin-related protein 2/3 (Arp2/3)-generated branched actin networks, by capping the actin filament barbed end, which promotes the Arp2/3-dependent filament nucleation and optimal branching [69,70,71].
Therefore, both ARHGAP21 and CAPG have a crucial role in actin assembly control, through the Arp2/3 complex regulation.
In conclusion, this study identified CAPG as an excellent candidate gene in ASD susceptibility and allowed the identification of new promising candidate genes encoding for PSD proteins, thus contributing to the comprehensive framework for understanding this complex and heterogeneous disorder. Indeed, we can hypothesize that the injurious effect of inherited CAPG deletions can be modified by the specific genetic background of the carrier individual. Therefore, in the Tunisian ASD proband [31], the presence of CAPG homozygous deletion could be sufficient to determine the autistic phenotype, while, in the two ASD Sardinian siblings here described, the CAPG deletion might act in combination with the VDAC3 2 bp-frameshift deletion and/or with the other identified deleterious variants in the PSD genes (PLXNA2, KCTD16, ARHGAP21, and SLC4A1). However, further studies are warranted to better clarify the role of CAPG, alone and in combination with the other identified PSD genes, in the impairment of synaptic homeostasis that is considered a common biological process associated with ASDs.

Supplementary Materials

The following are available online at https://www.mdpi.com/2077-0383/8/2/212/s1, Figure S1: Validation of chr3q26.32 (KCNMB2) deletion and chr6q21 (BEND3-PDSS2) duplication in family AUT003. (A) The genomic position of KCNMB2 deletion, chr3:178246465-178511664 (NCBI build 37 coordinates), is highlighted in blue. The microdeletion was confirmed in the mother and the affected son AUT003_3 by qPCR using four different probes. (B) The genomic position of BEND3-PDSS2 duplication, chr6:107408585-107546262 (NCBI build 37 coordinates), is highlighted in blue. The microduplication was confirmed in the father and the affected son AUT003_3 by qPCR using four different probes.

Author Contributions

Conceptualization, G.D. and P.Z.; Data curation, E.B., C.C., L.M., S.B., A.F., R.C., G.D., and P.Z.; Formal analysis, E.B., E.L., C.C., A.F.V.B., E.B., S.M., R.F., and P.Z.; Funding acquisition, E.M., G.D., and P.Z.; Investigation, E.B., E.L., C.C., L.M., A.F.V.B., R.F., E.M., G.D., and P.Z.; Methodology, E.B., R.F., G.D., and P.Z.; Project administration, P.Z.; Resources, G.D.; Software, E.B., C.C., and A.F.; Supervision, E.B., E.M., G.D., and P.Z.; Validation, E.B.; Writing—original draft, E.B., E.L., C.C., A.F.V.B., S.B., R.F., and P.Z.; Writing—review and editing, E.B., E.L., C.C., E.M., G.D., and P.Z.

Funding

This work was supported by grants from Fondazione Banco di Sardegna (FBS) (1175/2009.0470) to GD, from Associazione Italiana Studio Malformazioni (ASM) (2010) to PZ, from the Kyulan Family Foundation to EBa and from CINECA (computational resources to EBa). Partly supported also by grants from Fondo per la Ricerca Locale (ex 60%), Università di Cagliari, to PZ and from “Fondo per la Ricerca Fondamentale Orientata (ex quota 60%)”, University of Bologna, to EM.

Acknowledgments

We gratefully acknowledge the AUT003 family and all the subjects who have participated in the study. We thank the Centro Interdipartimentale di Ricerche sul Cancro “Giorgio Prodi” (CIRC), University of Bologna for Illumina Sequencing Service. We wish to thank Manuela and Francesco Zavattari for their skillful contribution to the drawing and graphic design of the graphical abstract.

Conflicts of Interest

The authors declare no potential conflicts of interest.

References

  1. Schaefer, G.B.; Mendelsohn, N.J. Professional practice and guidelines committee clinical genetics evaluation in identifying the etiology of autism spectrum disorders: 2013 guideline revisions. Genet. Med. 2013, 15, 399–407. [Google Scholar] [CrossRef] [PubMed]
  2. Caglayan, A.O. Genetic causes of syndromic and non-syndromic autism. Dev. Med. Child Neurol. 2010, 52, 130–138. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. El-Fishawy, P.; State, M.W. The genetics of autism: Key issues, recent findings and clinical implications. Psychiatry Clin North Am. 2011, 33, 83–105. [Google Scholar] [CrossRef] [PubMed]
  4. Geschwind, D.H. Genetics of autism spectrum disorders. Trends Cogn. Sci. 2011, 15, 409–416. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. World Health Organization: WHO. Available online: https://www.who.int/news-room/fact-sheets/detail/autism-spectrum-disorders (accessed on 7 December 2018).
  6. Loomes, R.; Hull, L.; Mandy, W.P.L. What is the male-to-female ratio in autism spectrum disorder? A systematic review and meta-analysis. J. Am. Acad. Child Adolesc. Psychiatry 2017, 56, 466–474. [Google Scholar] [CrossRef] [PubMed]
  7. Yin, J.; Schaaf, C.P. Autism genetics—An overview. Prenat. Diagn. 2016, 5, 14–30. [Google Scholar] [CrossRef] [PubMed]
  8. Toro, R.; Konyukh, M.; Delorme, R.; Leblond, C.; Chaste, P.; Fauchereau, F.; Coleman, M.; Leboyer, M.; Gillberg, C.; Bourgeron, T. Key role for gene dosage and synaptic homeostasis in autism spectrum disorders. Trends Genet. 2010, 26, 363–372. [Google Scholar] [CrossRef]
  9. Gilman, S.R.; Iossifov, I.; Levy, D.; Ronemus, M.; Wigler, M.; Vitkup, D. Rare de novo variants associated with autism implicate a large functional network of genes involved in formation and function of synapses. Neuron 2011, 70, 898–907. [Google Scholar] [CrossRef]
  10. Zavattari, P.; Lampis, R.; Mulargia, A.; Loddo, M.; Angius, E.; Todd, J.A.; Cucca, F. Confirmation of the DRB1-DQB1 loci as the major component of IDDM1 in the isolated founder population of Sardinia. Hum. Mol. Genet. 2000, 9, 2967–2972. [Google Scholar] [CrossRef] [Green Version]
  11. Contu, D.; Morelli, L.; Zavattari, P.; Lampis, R.; Angius, E.; Frongia, P.; Murru, D.; Maioli, M.; Francalacci, P.; Todd, J.A.; et al. Sex-related bias and exclusion mapping of the nonrecombinant portion of chromosome Y in human type 1 diabetes in the isolated founder population of Sardinia. Diabetes 2002, 51, 3573–3576. [Google Scholar] [CrossRef]
  12. Pitzalis, M.; Zavattari, P.; Murru, R.; Deidda, E.; Zoledziewska, M.; Murru, D.; Moi, L.; Motzo, C.; Orrù, V.; Costa, G.; et al. Genetic loci linked to Type 1 Diabetes and Multiple Sclerosis families in Sardinia. BMC Med. Genet. 2008, 9, 3. [Google Scholar] [CrossRef] [PubMed]
  13. Zavattari, P.; Loche, A.; Civolani, P.; Pilia, S.; Moi, L.; Casini, M.R.; Minerba, L.; Loche, S. An INSIG2 polymorphism affects glucose homeostasis in Sardinian obese children and adolescents. Ann. Hum. Genet. 2010, 74, 381–386. [Google Scholar] [CrossRef] [PubMed]
  14. Zavattari, P.; Loche, A.; Pilia, S.; Ibba, A.; Moi, L.; Guzzetti, C.; Casini, M.R.; Loche, S. rs9939609 in the FTO gene is associated with obesity but not with several biochemical parameters in Sardinian obese children. Ann. Hum. Genet. 2011, 75, 648–654. [Google Scholar] [CrossRef] [PubMed]
  15. Loudianos, G.; Dessi, V.; Lovicu, M.; Angius, A.; Cao, A.; Pirastu, M. The -75A-->C substitution in the 5′ UTR of the Wilson disease gene is a sequence polymorphism in the Mediterranean population. Am. J. Hum. Genet. 1998, 62, 484–485. [Google Scholar] [CrossRef]
  16. Crisponi, L.; Crisponi, G.; Meloni, A.; Toliat, M.R.; Nürnberg, G.; Usala, G.; Uda, M.; Masala, M.; Höhne, W.; Becker, C.; et al. Crisponi syndrome is caused by mutations in the CRLF1 gene and is allelic to cold-induced sweating syndrome Type 1. Am. J. Hum. Genet. 2007, 80, 971–981. [Google Scholar] [CrossRef] [PubMed]
  17. Tammimies, K.; Marshall, C.R.; Walker, S.; Kaur, G.; Thiruvahindrapuram, B.; Lionel, A.C.; Yuen, R.K.C.; Uddin, M.; Roberts, W.; Weksberg, R.; et al. Molecular diagnostic yield of chromosomal microarray analysis and whole-exome sequencing in children with autism spectrum disorder. JAMA 2015, 314, 895. [Google Scholar] [CrossRef] [PubMed]
  18. Lord, C.; Rutter, M.; DiLavore, P.C.; Risi, S.; Gotham, K.B.S. (ADOS-2) Manual (Part I): Modules 1–4. In Autism Diagnostic Observation Schedule, 2nd ed.; Western Psychological Services: Torrance, CA, USA, 2012. [Google Scholar]
  19. Guo, Y.; He, J.; Zhao, S.; Wu, H.; Zhong, X.; Sheng, Q.; Samuels, D.C.; Shyr, Y.; Long, J. Illumina human exome genotyping array clustering and quality control. Nat. Protoc. 2014, 9, 2643–2662. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Wang, K.; Li, M.; Hadley, D.; Liu, R.; Glessner, J.; Grant, S.F.A.; Hakonarson, H.; Bucan, M. PennCNV: An integrated hidden Markov model designed for high-resolution copy number variation detection in whole-genome SNP genotyping data. Genome Res. 2007, 17, 1665–1674. [Google Scholar] [CrossRef] [Green Version]
  21. Colella, S.; Yau, C.; Taylor, J.M.; Mirza, G.; Butler, H.; Clouston, P.; Bassett, A.S.; Seller, A.; Holmes, C.C.; Ragoussis, J. QuantiSNP: An objective bayes hidden-markov model to detect and accurately map copy number variation using SNP genotyping data. Nucleic Acids Res. 2007, 35, 2013–2025. [Google Scholar] [CrossRef]
  22. Zarrei, M.; MacDonald, J.R.; Merico, D.; Scherer, S.W. A copy number variation map of the human genome. Nat. Rev. Genet. 2015, 16, 172–183. [Google Scholar] [CrossRef]
  23. Chiara, M.; Gioiosa, S.; Chillemi, G.; D’Antonio, M.; Flati, T.; Picardi, E.; Zambelli, F.; Horner, D.S.; Pesole, G.; Castrignanò, T. CoVaCS: A consensus variant calling system. BMC Genomics 2018, 19, 120. [Google Scholar] [CrossRef] [PubMed]
  24. Wang, K.; Li, M.; Hakonarson, H. ANNOVAR: Functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010, 38, e164. [Google Scholar] [CrossRef] [PubMed]
  25. Genome aggregation database (gnomAD). Available online: http://gnomad.broadinstitute.org/ (accessed on 8 February 2018).
  26. 1000 Genomes Project. Available online: http://www.1000genomes.org/ (accessed on 8 February 2018).
  27. Yuen, R.K.; Merico, D.; Cao, H.; Pellecchia, G.; Alipanahi, B.; Thiruvahindrapuram, B.; Tong, X.; Sun, Y.; Cao, D.; Zhang, T.; et al. Genome-wide characteristics of de novo mutations in autism. NPJ Genomic Med. 2016, 1, 16027. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Jansen, A.; Dieleman, G.C.; Smit, A.B.; Verhage, M.; Verhulst, F.C.; Polderman, T.J.C.; Posthuma, D. Gene-set analysis shows association between FMRP targets and autism spectrum disorder. Eur. J. Hum. Genet. 2017, 25, 863–868. [Google Scholar] [CrossRef] [PubMed]
  29. Bayés, A.; van de Lagemaat, L.N.; Collins, M.O.; Croning, M.D.R.; Whittle, I.R.; Choudhary, J.S.; Grant, S.G.N. Characterization of the proteome, diseases and evolution of the human postsynaptic density. Nat. Neurosci. 2011, 14, 19–21. [Google Scholar] [CrossRef] [PubMed]
  30. Pinto, D.; Delaby, E.; Merico, D.; Barbosa, M.; Merikangas, A.; Klei, L.; Thiruvahindrapuram, B.; Xu, X.; Ziman, R.; Wang, Z.; et al. Convergence of genes and cellular pathways dysregulated in autism spectrum disorders. Am. J. Hum. Genet. 2014, 94, 677–694. [Google Scholar] [CrossRef] [PubMed]
  31. Lahbib, S.; Leblond, C.S.; Hamza, M.; Regnault, B.; Lemée, L.; Mathieu, A.; Jaouadi, H.; Mkaouar, R.; Youssef-Turki, I.; Belhadj, A.B.; et al. Homozygous 2p11.2 deletion supports the implication of ELMOD3 in hearing loss and reveals the potential association of CAPG with ASD/ID etiology. J. Appl. Genet. 2018, 60, 46–56. [Google Scholar] [CrossRef] [PubMed]
  32. Bierut, L.J.; Agrawal, A.; Bucholz, K.K.; Doheny, K.F.; Laurie, C.; Pugh, E.; Fisher, S.; Fox, L.; Howells, W.; Bertelsen, S.; et al. A genome-wide association study of alcohol dependence. Proc. Natl. Acad. Sci. USA 2010, 107, 5082–5087. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Krawczak, M.; Nikolaus, S.; von Eberstein, H.; Croucher, P.J.P.; El Mokhtari, N.E.; Schreiber, S. PopGen: Population-based recruitment of patients and controls for the analysis of complex genotype-phenotype relationships. Public Health Genomics 2006, 9, 55–61. [Google Scholar] [CrossRef] [PubMed]
  34. Shaikh, T.H.; Gai, X.; Perin, J.C.; Glessner, J.T.; Xie, H.; Murphy, K.; O’Hara, R.; Casalunovo, T.; Conlin, L.K.; D’Arcy, M.; et al. High-resolution mapping and analysis of copy number variations in the human genome: A data resource for clinical and research applications. Genome Res. 2009, 19, 1682–1690. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Stewart, A.F.R.; Dandona, S.; Chen, L.; Assogba, O.; Belanger, M.; Ewart, G.; LaRose, R.; Doelle, H.; Williams, K.; Wells, G.A.; et al. Kinesin family member 6 variant Trp719Arg does not associate with angiographically defined coronary artery disease in the Ottawa heart genomics study. J. Am. Coll. Cardiol. 2009, 53, 1471–1472. [Google Scholar] [CrossRef] [PubMed]
  36. Figueiredo, J.C.; Lewinger, J.P.; Song, C.; Campbell, P.T.; Conti, D.V.; Edlund, C.K.; Duggan, D.J.; Rangrej, J.; Lemire, M.; Hudson, T.; et al. Genotype-environment interactions in microsatellite stable/microsatellite instability-low colorectal cancer: Results from a genome-wide association study. Cancer Epidemiol. Biomark. Prev. 2011, 20, 758–766. [Google Scholar] [CrossRef] [PubMed]
  37. Lek, M.; Karczewski, K.J.; Minikel, E.V.; Samocha, K.E.; Banks, E.; Fennell, T.; O’Donnell-Luria, A.H.; Ware, J.S.; Hill, A.J.; Cummings, B.B.; et al. Analysis of protein-coding genetic variation in 60,706 humans. Nature 2016, 536, 285–291. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Fan, Y.; Tang, X.; Vitriol, E.; Chen, G.; Zheng, J.Q. Actin capping protein is required for dendritic spine development and synapse formation. J. Neurosci. 2011, 31, 10228–10233. [Google Scholar] [CrossRef] [PubMed]
  39. Van Spronsen, M.; Hoogenraad, C.C. Synapse pathology in psychiatric and neurologic disease. Curr. Neurol. Neurosci. Rep. 2010, 10, 207–214. [Google Scholar] [CrossRef] [PubMed]
  40. Korobova, F.; Svitkina, T. Molecular architecture of synaptic actin cytoskeleton in hippocampal neurons reveals a mechanism of dendritic spine morphogenesis. Mol. Biol. Cell 2010, 21, 165–176. [Google Scholar] [CrossRef] [PubMed]
  41. Bar-shira, O.; Chechik, G. Predicting protein–protein interactions in the post synaptic density. Mol. Cell. Neurosci. 2013, 56, 128–139. [Google Scholar] [CrossRef]
  42. Ehrhart, F.; Coort, S.L.; Eijssen, L.; Cirillo, E.; Smeets, E.; Sangani, N.B.; Evelo, C.; Curfs, L.M.G. Integrated analysis of human transcriptome data for Rett syndrome finds a network of involved genes. BioRxiv 2018. [Google Scholar]
  43. Matarazzo, V.; Cohen, D.; Palmer, A.M.; Simpson, P.J.; Khokhar, B.; Pan, S.-J.; Ronnett, G.V. The transcriptional repressor Mecp2 regulates terminal neuronal differentiation. Mol. Cell. Neurosci. 2004, 27, 44–58. [Google Scholar] [CrossRef]
  44. Gonzalez-Gronow, M.; Cuchacovich, M.; Francos, R.; Cuchacovich, S.; del Fernandez, M.P.; Blanco, A.; Bowers, E.V.; Kaczowka, S.; Pizzo, S.V. Antibodies against the voltage-dependent anion channel (VDAC) and its protective ligand hexokinase-I in children with autism. J. Neuroimmunol. 2010, 227, 153–161. [Google Scholar] [CrossRef] [Green Version]
  45. Gvozdjáková, A.; Low, H.; Sun, I.; Navas, P.; Crane, F. Plasma membrane coenzyme Q: Evidence for a role in autism. Biol. Targets Ther. 2014, 8, 199. [Google Scholar] [CrossRef] [PubMed]
  46. Thinnes, F.P. Plasmalemmal VDAC-1 corroborated as amyloid Aß-receptor. Front. Aging Neurosci. 2015, 7, 188. [Google Scholar] [CrossRef] [PubMed]
  47. Su, A.I.; Wiltshire, T.; Batalov, S.; Lapp, H.; Ching, K.A.; Block, D.; Zhang, J.; Soden, R.; Hayakawa, M.; Kreiman, G.; et al. A gene atlas of the mouse and human protein-encoding transcriptomes. Proc. Natl. Acad. Sci. USA 2004, 101, 6062–6067. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Suda, S.; Iwata, K.; Shimmura, C.; Kameno, Y.; Anitha, A.; Thanseem, I.; Nakamura, K.; Matsuzaki, H.; Tsuchiya, K.J.; Sugihara, G.; et al. Decreased expression of axon-guidance receptors in the anterior cingulate cortex in autism. Mol. Autism 2011, 2, 14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Mah, S.; Nelson, M.R.; DeLisi, L.E.; Reneland, R.H.; Markward, N.; James, M.R.; Nyholt, D.R.; Hayward, N.; Handoko, H.; Mowry, B.; et al. Identification of the semaphorin receptor PLXNA2 as a candidate for susceptibility to schizophrenia. Mol. Psychiatry 2006, 11, 471–478. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Zhao, X.-F.; Kohen, R.; Parent, R.; Duan, Y.; Fisher, G.L.; Korn, M.J.; Ji, L.; Wan, G.; Jin, J.; Püschel, A.W.; et al. PlexinA2 forward signaling through Rap1 GTPases regulates dentate gyrus development and schizophrenia-like behaviors. Cell Rep. 2018, 22, 456–470. [Google Scholar] [CrossRef] [PubMed]
  51. Lilja, J.; Zacharchenko, T.; Georgiadou, M.; Jacquemet, G.; Franceschi, N.D.; Peuhu, E.; Hamidi, H.; Pouwels, J.; Martens, V.; Nia, F.H.; et al. SHANK proteins limit integrin activation by directly interacting with Rap1 and R-Ras. Nat. Cell Biol. 2017, 19, 292–305. [Google Scholar] [CrossRef]
  52. Booker, S.A.; Althof, D.; Gross, A.; Loreth, D.; Müller, J.; Unger, A.; Fakler, B.; Varro, A.; Watanabe, M.; Gassmann, M.; et al. KCTD12 auxiliary proteins modulate kinetics of GABA B receptor-mediated inhibition in cholecystokinin-containing interneurons. Cereb. Cortex 2016, 27, bhw090. [Google Scholar] [CrossRef]
  53. Fritzius, T.; Turecek, R.; Seddik, R.; Kobayashi, H.; Tiao, J.; Rem, P.D.; Metz, M.; Kralikova, M.; Bouvier, M.; Gassmann, M.; et al. KCTD hetero-oligomers confer unique kinetic properties on hippocampal GABAB receptor-induced K+ Currents. J. Neurosci. 2017, 37, 1162–1175. [Google Scholar] [CrossRef]
  54. Schwenk, J.; Pérez-Garci, E.; Schneider, A.; Kollewe, A.; Gauthier-Kemper, A.; Fritzius, T.; Raveh, A.; Dinamarca, M.C.; Hanuschkin, A.; Bildl, W.; et al. Modular composition and dynamics of native GABAB receptors identified by high-resolution proteomics. Nat. Neurosci. 2016, 19, 233–242. [Google Scholar] [CrossRef]
  55. Schwenk, J.; Metz, M.; Zolles, G.; Turecek, R.; Fritzius, T.; Bildl, W.; Tarusawa, E.; Kulik, A.; Unger, A.; Ivankova, K.; et al. Native GABAB receptors are heteromultimers with a family of auxiliary subunits. Nature 2010, 465, 231–235. [Google Scholar] [CrossRef] [PubMed]
  56. Angelicheva, D.; Tournev, I.; Guergueltcheva, V.; Mihaylova, V.; Azmanov, D.N.; Morar, B.; Radionova, M.; Smith, S.J.; Zlatareva, D.; Stevens, J.M.; et al. Partial epilepsy syndrome in a Gypsy family linked to 5q31.3-q32. Epilepsia 2009, 50, 1679–1688. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Hentschke, M.; Wiemann, M.; Hentschke, S.; Kurth, I.; Hermans-Borgmeyer, I.; Seidenbecher, T.; Jentsch, T.J.; Gal, A.; Hübner, C.A. Mice with a targeted disruption of the Cl-/HCO3- exchanger AE3 display a reduced seizure threshold. Mol. Cell. Biol. 2006, 26, 182–191. [Google Scholar] [CrossRef] [PubMed]
  58. Jacobs, S.; Ruusuvuori, E.; Sipila, S.T.; Haapanen, A.; Damkier, H.H.; Kurth, I.; Hentschke, M.; Schweizer, M.; Rudhard, Y.; Laatikainen, L.M.; et al. Mice with targeted Slc4a10 gene disruption have small brain ventricles and show reduced neuronal excitability. Proc. Natl. Acad. Sci. USA 2008, 105, 311–316. [Google Scholar] [CrossRef] [PubMed]
  59. Sinning, A.; Liebmann, L.; Hübner, C.A. Disruption of Slc4a10 augments neuronal excitability and modulates synaptic short-term plasticity. Front. Cell. Neurosci. 2015, 9, 223. [Google Scholar] [CrossRef]
  60. Sinning, A.; Liebmann, L.; Kougioumtzes, A.; Westermann, M.; Bruehl, C.; Hübner, C.A. Synaptic glutamate release is modulated by the Na+ -driven Cl-/HCO₃− exchanger Slc4a8. J. Neurosci. 2011, 31, 7300–7311. [Google Scholar] [CrossRef]
  61. Sander, T.; Toliat, M.R.; Heils, A.; Leschik, G.; Becker, C.; Rüschendorf, F.; Rohde, K.; Mundlos, S.; Nürnberg, P. Association of the 867Asp variant of the human anion exchanger 3 gene with common subtypes of idiopathic generalized epilepsy. Epilepsy Res. 2002, 51, 249–255. [Google Scholar] [CrossRef]
  62. Demirci, F.Y.K.; Chang, M.-H.; Mah, T.S.; Romero, M.F.; Gorin, M.B. Proximal renal tubular acidosis and ocular pathology: A novel missense mutation in the gene (SLC4A4) for sodium bicarbonate cotransporter protein (NBCe1). Mol. Vis. 2006, 12, 324–330. [Google Scholar]
  63. Sebat, J.; Lakshmi, B.; Malhotra, D.; Troge, J.; Lese-Martin, C.; Walsh, T.; Yamrom, B.; Yoon, S.; Krasnitz, A.; Kendall, J.; et al. Strong association of de novo copy number mutations with autism. Science 2007, 316, 445–449. [Google Scholar] [CrossRef]
  64. Gurnett, C.A.; Veile, R.; Zempel, J.; Blackburn, L.; Lovett, M.; Bowcock, A. Disruption of sodium bicarbonate transporter SLC4A10 in a patient with complex partial epilepsy and mental retardation. Arch. Neurol. 2008, 65, 550. [Google Scholar] [CrossRef]
  65. Krepischi, A.C.V.; Knijnenburg, J.; Bertola, D.R.; Kim, C.A.; Pearson, P.L.; Bijlsma, E.; Szuhai, K.; Kok, F.; Vianna-Morgante, A.M.; Rosenberg, C. Two distinct regions in 2q24.2-q24.3 associated with idiopathic epilepsy. Epilepsia 2010, 51, 2457–2460. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Suzuki, M.; Van Paesschen, W.; Stalmans, I.; Horita, S.; Yamada, H.; Bergmans, B.A.; Legius, E.; Riant, F.; De Jonghe, P.; Li, Y.; et al. Defective membrane expression of the Na+ -HCO3- cotransporter NBCe1 is associated with familial migraine. Proc. Natl. Acad. Sci. USA 2010, 107, 15963–15968. [Google Scholar] [CrossRef] [PubMed]
  67. Bassères, D.S.; Tizzei, E.V.; Duarte, A.A.; Costa, F.F.; Saad, S.T.O. ARHGAP10, a novel human gene coding for a potentially cytoskeletal Rho-GTPase activating protein. Biochem. Biophys. Res. Commun. 2002, 294, 579–585. [Google Scholar] [CrossRef]
  68. Dubois, T.; Paléotti, O.; Mironov, A.A.; Fraisier, V.; Stradal, T.E.B.; De Matteis, M.A.; Franco, M.; Chavrier, P. Golgi-localized GAP for Cdc42 functions downstream of ARF1 to control Arp2/3 complex and F-actin dynamics. Nat. Cell Biol. 2005, 7, 353–364. [Google Scholar] [CrossRef] [PubMed]
  69. Akin, O.; Mullins, R.D. Capping protein increases the rate of actin-based motility by promoting filament nucleation by the Arp2/3 complex. Cell 2008, 133, 841–851. [Google Scholar] [CrossRef] [PubMed]
  70. Iwasa, J.H.; Mullins, R.D. Spatial and temporal relationships between actin-filament nucleation, capping, and disassembly. Curr. Biol. 2007, 17, 395–406. [Google Scholar] [CrossRef] [PubMed]
  71. Mejillano, M.R.; Kojima, S.; Applewhite, D.A.; Gertler, F.B.; Svitkina, T.M.; Borisy, G.G. Lamellipodial versus filopodial mode of the actin nanomachinery. Cell 2004, 118, 363–373. [Google Scholar] [CrossRef]
Jcm 08 00212 g001 550
Figure 1. Validation and breakpoint definition of chr2p11.2 (ELMOD3-CAPG-SH2D6) deletion in family AUT003. (A) The microdeletion was confirmed in the unaffected mother and the two affected siblings by quantitative PCR using four different probes. The genomic position of the deletion is indicated with a red rectangle; (B) Sequence electropherogram: the 21-bp breakpoint-spanning sequence, common to both ends, is pinpointed by a blue rectangle (NCBI build 37 coordinates). DEL: assay located in the deletion; NONdel: assay located outside the deletion.
Figure 1. Validation and breakpoint definition of chr2p11.2 (ELMOD3-CAPG-SH2D6) deletion in family AUT003. (A) The microdeletion was confirmed in the unaffected mother and the two affected siblings by quantitative PCR using four different probes. The genomic position of the deletion is indicated with a red rectangle; (B) Sequence electropherogram: the 21-bp breakpoint-spanning sequence, common to both ends, is pinpointed by a blue rectangle (NCBI build 37 coordinates). DEL: assay located in the deletion; NONdel: assay located outside the deletion.
Jcm 08 00212 g001
Jcm 08 00212 g002 550
Figure 2. CAPG and ELMOD3 gene expression levels in different tissues. (A) CAPG expression levels; (B) ELMOD3 gene expression levels Bar plots represent fold change with upper/lower limits relative to TFRC reference gene.
Figure 2. CAPG and ELMOD3 gene expression levels in different tissues. (A) CAPG expression levels; (B) ELMOD3 gene expression levels Bar plots represent fold change with upper/lower limits relative to TFRC reference gene.
Jcm 08 00212 g002
Jcm 08 00212 g003 550
Figure 3. CAPG and ELMOD3 gene expression levels in AUT003 family members and in controls. (A) CAPG expression levels; (B) ELMOD3 gene expression levels. Bar plots represent fold change with upper/lower limits relative to controls. The statistical comparisons were conducted between AUT003 family members and the average expression level between the two controls. Asterisks indicate statistically significant differences (** p <0.01, *** p <0.001).
Figure 3. CAPG and ELMOD3 gene expression levels in AUT003 family members and in controls. (A) CAPG expression levels; (B) ELMOD3 gene expression levels. Bar plots represent fold change with upper/lower limits relative to controls. The statistical comparisons were conducted between AUT003 family members and the average expression level between the two controls. Asterisks indicate statistically significant differences (** p <0.01, *** p <0.001).
Jcm 08 00212 g003
Jcm 08 00212 g004 550
Figure 4. CAPG gene expression levels in autism spectrum disorder (ASD) cases and in controls. Boxplot of the relative CAPG expression of 13 ASD cases and of 22 healthy controls. The expression values were normalized using TFRC reference gene. Mean values are represented by crosses. The asterisk indicates a statistically significant difference (p = 0.026).
Figure 4. CAPG gene expression levels in autism spectrum disorder (ASD) cases and in controls. Boxplot of the relative CAPG expression of 13 ASD cases and of 22 healthy controls. The expression values were normalized using TFRC reference gene. Mean values are represented by crosses. The asterisk indicates a statistically significant difference (p = 0.026).
Jcm 08 00212 g004
Jcm 08 00212 g005 550
Figure 5. CAPG protein expression levels in AUT003 family members and controls. (A) Representative blots of CAPG protein expression in autism spectrum disorder (ASD) cases and controls. GAPDH was used as loading control; (B) CAPG expression. Bars represent mean ± SD (standard deviation). The statistical comparisons were made between AUT003 family members and the average between the two controls. Asterisks indicate statistically significant differences (* p <0.05, ** p <0.01).
Figure 5. CAPG protein expression levels in AUT003 family members and controls. (A) Representative blots of CAPG protein expression in autism spectrum disorder (ASD) cases and controls. GAPDH was used as loading control; (B) CAPG expression. Bars represent mean ± SD (standard deviation). The statistical comparisons were made between AUT003 family members and the average between the two controls. Asterisks indicate statistically significant differences (* p <0.05, ** p <0.01).
Jcm 08 00212 g005
Table
Table 1. Assessment of parents’ autism phenotype using the Broad Autism Phenotype Questionnaire, the Social Responsiveness Scale, and the Autism Diagnostic Observation Schedule.
Table 1. Assessment of parents’ autism phenotype using the Broad Autism Phenotype Questionnaire, the Social Responsiveness Scale, and the Autism Diagnostic Observation Schedule.
MeasureMotherFather
BAPQ Scores
Aloof2.43.5
Pragmatic Language2.22.7
Rigidity2.22.7
BAPQ Total2.43
SRS: ARV scores
Awareness108
Mannerism810
Communication921
Cognition410
Motivation310
SRS:ARV Total3459
ADOS—IV Module Scores
Communication00
Reciprocal Social Interaction00
BAPQ: Broad Autism Phenotype Questionnaire; SRS: Social Responsiveness Scale; ARV: Adult Research Version; ADOS: Autism Diagnostic Observation Schedule–Generic. BAPQ, SRS: ARV, and ADOS–IV module scores are shown.
Table
Table 2. LGD variants.
Table 2. LGD variants.
Identified VariantpLi ScoreSynaptic Genes [28,29]Genotype
Gene Base ChangeAminoAcid ChangeGeneAUT003.3AUT003.4
chr1:g.19559158delA NM_001271428:c.1739delTp.F580Sfs*10EMC1 0/10/0
chr1:g.24779922delT NM_001322854:c.565delTp.L189Cfs*15NIPAL3 0/10/0
chr1:g.108771659G>A NM_001143989:c.C1543Tp.R515XNBPF4 0/10/1
chr1:g.158299688A>T NM_001764:c.T561Ap.Y187XCD1B 0/10/1
chr2:g.30961265C>T NM_144575:c.1594+1G>A CAPN13 0/10/0
chr2:g.58386928-58386929insTAAT NM_001114636:c.1114_1115insATTAp.T372Nfs*12FANCL 0/10/1
chr2:g.120439282C>T NM_001105198:c.C853Tp.R285XTMEM177 0/10/1
chr3:g.38770224G>A NM_001293306:c.C2449Tp.R817XSCN10A 0/00/1
chr3:g.188327410T>A NM_001167671:c.T891Ap.Y297XLPP 0/00/1
chr4:g.8235241G>A NM_001318480:c.3054+1G>A SH3TC1 0/10/1
chr6:43306235-43306238delTCTT NM_014345:c.5498_5501delp.K1833Rfs*6ZNF318 0/00/1
chr7:g.150028118C>T NM_138434:c.C625Tp.Q209XZBED6CL 0/10/0
chr7:g.150864372G>T NM_001098834:c.C264Ap.C88XGBX1 0/00/1
chr8:g.42259530-42259531delAT NM_001135694:c.551_552delp.H184Rfs*5VDAC30.97PSD; Excitability0/10/1
chr8:g.145140986delC NM_003801:c.1824delCp.C610Afs*43GPAA1 0/00/1
chr10:g.102240866-102240872del NM_003393:c.353_359delp.G120Wfs*37WNT8B 0/10/0
chr10:g.123845143-123845144delCC NM_001291876:c.3128_3129delp.P1044Tfs*4TACC2 0/00/1
chr11:g.59577372dupC NM_017840:c.76dupGp.A26Gfs*14MRPL16 0/10/0
chr11:g.116701354G>A NM_000040:c.55+1G>A APOC3 0/00/1
chr12:g.25705804C>T NM_001145728:c.89+1G>A LMNTD1 0/00/1
chr12:g.54757526delC NM_020370:c.110delGp.G37Afs*22GPR84 0/10/0
chr12:g.64436690C>T NM_001346201:c.C610TR204XSRGAP10.99 0/10/0
chr12:g.133198246dupA NM_012226:c.889dupAp.N298Kfs*6P2RX2 0/00/1
chr15:g.75641682T>C NM_024608:c.434+2T>C NEIL1 0/00/1
chr16:g.107113-107114insTC NM_024571:c.369_370insTCp.F125Pfs*5SNRNP25 0/10/1
chr18:g.71822640delT NM_014177:c.462+2T>- TIMM21 0/00/1
chr18:g.77917844dupG NM_032510:c.940dupCp.Q314Pfs*72PARD6G 0/00/1
chr19:g.10204081dupA NM_001321411:c.1165dupTp.S389Ffs*4ANGPTL6 0/00/1
chr19:g.43859812G>C NM_020406:c.380-1G>C CD177 0/10/1
chr20:g.61591931delT NM_022082:c.473delTp.I158Tfs*11SLC17A9 0/00/1
chr21:g.44839738C>A NM_001320643:c.1119+1G>T SIK10.99 0/10/1
chr22:g.25124143C>T NM_001255975:c.1905+1G>A PIWIL3 0/00/1
chrX:g.23019507dupT NM_182699:c.1334dupTp.I446Nfs*35DDX53 1/11/1
LGD: Likely Gene Disrupting; PSD: postsynaptic density; pLi: probability of being loss-of-function (LoF) intolerant (pLi).
Table
Table 3. Rare damaging missense variants, shared by both affected siblings, in genes with pLi score ≥0.9.
Table 3. Rare damaging missense variants, shared by both affected siblings, in genes with pLi score ≥0.9.
Identified VariantpLi ScoreSynaptic Genes [28,29]Genotype
Gene Base ChangeAminoAcid ChangeGeneAUT003.3AUT003.4
chr1:g.109801543A>G NM_001408:c.A3800Gp.H1267RCELSR20.9999992 0/10/1
chr1:g.117146504G>A NM_001007237:c.C1366Tp.R456CIGSF30.987084919 0/10/1
chr1:g.162725022G>A NM_006182:c.G494Ap.R165QDDR20.990992372 0/10/1
chr1:g.208390903G>T NM_025179:c.C365Ap.S122YPLXNA20.993981182PSD0/10/1
chr2:g.96956078C>T NM_014014:c.G2728Ap.V910ISNRNP2001 0/10/1
chr2:g.239975199C>T NM_006037:c.G3172Ap.A1058THDAC40.999989727 0/10/1
chr4:g.126238072T>C NM_001291285:c.T506Cp.V169AFAT40.999999998 0/10/1
chr5:g.141052986C>T NM_022481:c.G954Ap.M318IARAP30.967505029 0/10/1
chr5:g.143587022G>A NM_020768:c.G745Ap.A249TKCTD160.959972608PSD; Excitability0/10/1
chr7:g.100491451G>C NM_000665:c.C403Gp.P135AACHE0.996674942 0/10/1
chr9:g.88938095T>A NM_001185074:c.A2201Tp.E734VZCCHC60.999999563 0/10/1
chr10:g.70726833C>A NM_001256910:c.C760Ap.H254NDDX210.998039216 0/10/1
chr10:g.79581043G>A NM_004747:c.C3199Tp.R1067CDLG50.999998542 0/10/1
chr10:g.79590510C>G NM_004747:c.G1870Cp.E624QDLG50.999998542 0/10/1
chr10:g.24874455G>T NM_020824:c.C4763Ap.S1588YARHGAP210.999946673PSD0/10/1
chr11:g.94918634T>G NM_144665:c.A548Cp.N183TSESN30.994296083 0/10/1
chr11:g.114027109G>A NM_001018011:c.G1319Ap.R440QZBTB160.980087306 0/10/1
chr12:g.112699118C>T NM_001109662:c.G2432Ap.R811HHECTD41 0/10/1
chr17:g.42335485C>T NM_000342:c.G1151Ap.R384HSLC4A10.909484295PSD; Excitability0/10/1
chr19:g.50100192G>A NM_020719:c.G2600Ap.R867HPRR120.999969577 0/10/1
PSD: postsynaptic density; pLi: probability of being loss-of-function (LoF) intolerant (pLi).

Share and Cite

MDPI and ACS Style

Bacchelli, E.; Loi, E.; Cameli, C.; Moi, L.; Vega Benedetti, A.F.; Blois, S.; Fadda, A.; Bonora, E.; Mattu, S.; Fadda, R.; et al. Analysis of a Sardinian Multiplex Family with Autism Spectrum Disorder Points to Post-Synaptic Density Gene Variants and Identifies CAPG as a Functionally Relevant Candidate Gene. J. Clin. Med. 2019, 8, 212. https://doi.org/10.3390/jcm8020212

AMA Style

Bacchelli E, Loi E, Cameli C, Moi L, Vega Benedetti AF, Blois S, Fadda A, Bonora E, Mattu S, Fadda R, et al. Analysis of a Sardinian Multiplex Family with Autism Spectrum Disorder Points to Post-Synaptic Density Gene Variants and Identifies CAPG as a Functionally Relevant Candidate Gene. Journal of Clinical Medicine. 2019; 8(2):212. https://doi.org/10.3390/jcm8020212

Chicago/Turabian Style

Bacchelli, Elena, Eleonora Loi, Cinzia Cameli, Loredana Moi, Ana Florencia Vega Benedetti, Sylvain Blois, Antonio Fadda, Elena Bonora, Sandra Mattu, Roberta Fadda, and et al. 2019. "Analysis of a Sardinian Multiplex Family with Autism Spectrum Disorder Points to Post-Synaptic Density Gene Variants and Identifies CAPG as a Functionally Relevant Candidate Gene" Journal of Clinical Medicine 8, no. 2: 212. https://doi.org/10.3390/jcm8020212

APA Style

Bacchelli, E., Loi, E., Cameli, C., Moi, L., Vega Benedetti, A. F., Blois, S., Fadda, A., Bonora, E., Mattu, S., Fadda, R., Chessa, R., Maestrini, E., Doneddu, G., & Zavattari, P. (2019). Analysis of a Sardinian Multiplex Family with Autism Spectrum Disorder Points to Post-Synaptic Density Gene Variants and Identifies CAPG as a Functionally Relevant Candidate Gene. Journal of Clinical Medicine, 8(2), 212. https://doi.org/10.3390/jcm8020212

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop