Genomic analysis of 116 autism families strengthens known risk genes and highlights promising candidates.

Autism spectrum disorder (ASD) is a complex neurodevelopmental condition with a strong genetic component in which rare variants contribute significantly to risk. We performed whole genome and/or exome sequencing (WGS and WES) and SNP-array analysis to identify both rare sequence and copy number variants (SNVs and CNVs) in 435 individuals from 116 ASD families. We identified 37 rare potentially damaging de novo SNVs (pdSNVs) in the cases (n = 144). Interestingly, two of them (one stop-gain and one missense variant) occurred in the same gene, BRSK2. Moreover, the identification of 8 severe de novo pdSNVs in genes not previously implicated in ASD (AGPAT3, IRX5, MGAT5B, RAB8B, RAP1A, RASAL2, SLC9A1, YME1L1) highlighted promising candidates. Potentially damaging CNVs (pdCNVs) provided support to the involvement of inherited variants in PHF3, NEGR1, TIAM1 and HOMER1 in neurodevelopmental disorders (NDD), although mostly acting as susceptibility factors with incomplete penetrance. Interpretation of identified pdSNVs/pdCNVs according to the ACMG guidelines led to a molecular diagnosis in 19/144 cases, although this figure represents a lower limit and is expected to increase thanks to further clarification of the role of likely pathogenic variants in ASD/NDD candidate genes not yet established. In conclusion, our study highlights promising ASD candidate genes and contributes to characterize the allelic diversity, mode of inheritance and phenotypic impact of de novo and inherited risk variants in ASD/NDD genes.

Whole genome analysis of rare deleterious variants adds further evidence to BRSK2 and other risk genes in Autism Spectrum Disorder.

Autism spectrum disorder (ASD) is a complex neurodevelopmental condition with a strong genetic component in which rare variants contribute significantly to risk. We have performed whole genome and/or exome sequencing (WGS and WES) and SNP-array analysis to identify both rare sequence and copy number variants (SNVs and CNVs) in 435 individuals from 116 ASD families. We identified 37 rare potentially damaging de novo SNVs (pdSNVs) in cases (n = 144). Interestingly, two of them (one stop-gain and one missense variant) occurred in the same gene, BRSK2. Moreover, the identification of 9 severe de novo pdSNVs in genes not previously implicated in ASD (RASAL2, RAP1A, IRX5, SLC9A1, AGPAT3, MGAT3, RAB8B, MGAT5B, YME1L1), highlighted novel candidates. Potentially damaging CNVs (pdCNVs) provided support to the involvement of inherited variants in PHF3, NEGR1, TIAM1 and HOMER1 in neurodevelopmental disorders (NDD), although mostly acting as susceptibility factors with incomplete penetrance. Interpretation of identified pdSNVs/pdCNVs according to the ACMG guidelines led to a molecular diagnosis in 19/144 cases, but this figure represents a lower limit and is expected to increase thanks to further clarification of the role of likely pathogenic variants in new ASD/NDD candidates. In conclusion, our study strengthens the role of BRSK2 and other neurodevelopmental genes in ASD risk, highlights novel candidates and contributes to characterize the allelic diversity, mode of inheritance and phenotypic impact of de novo and inherited risk variants in ASD/NDD genes.

Dissecting the multifaceted contribution of the mitochondrial genome to autism spectrum disorder.

Autism spectrum disorder (ASD) is a clinically heterogeneous class of neurodevelopmental conditions with a strong, albeit complex, genetic basis. The genetic architecture of ASD includes different genetic models, from monogenic transmission at one end, to polygenic risk given by thousands of common variants with small effects at the other end. The mitochondrial DNA (mtDNA) was also proposed as a genetic modifier for ASD, mostly focusing on maternal mtDNA, since the paternal mitogenome is not transmitted to offspring. We extensively studied the potential contribution of mtDNA in ASD pathogenesis and risk through deep next generation sequencing and quantitative PCR in a cohort of 98 families. While the maternally-inherited mtDNA did not seem to predispose to ASD, neither for haplogroups nor for the presence of pathogenic mutations, an unexpected influence of paternal mtDNA, apparently centered on haplogroup U, came from the Italian families extrapolated from the test cohort (n = 74) when compared to the control population. However, this result was not replicated in an independent Italian cohort of 127 families and it is likely due to the elevated paternal age at time of conception. In addition, ASD probands showed a reduced mtDNA content when compared to their unaffected siblings. Multivariable regression analyses indicated that variants with 15%-5% heteroplasmy in probands are associated to a greater severity of ASD based on ADOS-2 criteria, whereas paternal super-haplogroups H and JT were associated with milder phenotypes. In conclusion, our results suggest that the mtDNA impacts on ASD, significantly modifying the phenotypic expression in the Italian population. The unexpected finding of protection induced by paternal mitogenome in term of severity may derive from a role of mtDNA in influencing the accumulation of nuclear de novo mutations or epigenetic alterations in fathers' germinal cells, affecting the neurodevelopment in the offspring. This result remains preliminary and needs further confirmation in independent cohorts of larger size. If confirmed, it potentially opens a different perspective on how paternal non-inherited mtDNA may predispose or modulate other complex diseases.

Contribution of CACNA1H Variants in Autism Spectrum Disorder Susceptibility.

Autism Spectrum Disorder (ASD) is a highly heterogeneous neuropsychiatric disorder with a strong genetic component. The genetic architecture is complex, consisting of a combination of common low-risk and more penetrant rare variants. Voltage-gated calcium channels (VGCCs or Cav) genes have been implicated as high-confidence susceptibility genes for ASD, in accordance with the relevant role of calcium signaling in neuronal function. In order to further investigate the involvement of VGCCs rare variants in ASD susceptibility, we performed whole genome sequencing analysis in a cohort of 105 families, composed of 124 ASD individuals, 210 parents and 58 unaffected siblings. We identified 53 rare inherited damaging variants in Cav genes, including genes coding for the principal subunit and genes coding for the auxiliary subunits, in 40 ASD families. Interestingly, biallelic rare damaging missense variants were detected in the CACNA1H gene, coding for the T-type Cav3.2 channel, in ASD probands from two different families. Thus, to clarify the role of these CACNA1H variants on calcium channel activity we performed electrophysiological analysis using whole-cell patch clamp technology. Three out of four tested variants were shown to mildly affect Cav3.2 channel current density and activation properties, possibly leading to a dysregulation of intracellular Ca2+ ions homeostasis, thus altering calcium-dependent neuronal processes and contributing to ASD etiology in these families. Our results provide further support for the role of CACNA1H in neurodevelopmental disorders and suggest that rare CACNA1H variants may be involved in ASD development, providing a high-risk genetic background.

An increased burden of rare exonic variants in NRXN1 microdeletion carriers is likely to enhance the penetrance for autism spectrum disorder.

Autism spectrum disorder (ASD) is characterized by a complex polygenic background, but with the unique feature of a subset of cases (~15%-30%) presenting a rare large-effect variant. However, clinical interpretation in these cases is often complicated by incomplete penetrance, variable expressivity and different neurodevelopmental trajectories. NRXN1 intragenic deletions represent the prototype of such ASD-associated susceptibility variants. From chromosomal microarrays analysis of 104 ASD individuals, we identified an inherited NRXN1 deletion in a trio family. We carried out whole-exome sequencing and deep sequencing of mitochondrial DNA (mtDNA) in this family, to evaluate the burden of rare variants which may contribute to the phenotypic outcome in NRXN1 deletion carriers. We identified an increased burden of exonic rare variants in the ASD child compared to the unaffected NRXN1 deletion-transmitting mother, which remains significant if we restrict the analysis to potentially deleterious rare variants only (P = 6.07 × 10-5 ). We also detected significant interaction enrichment among genes with damaging variants in the proband, suggesting that additional rare variants in interacting genes collectively contribute to cross the liability threshold for ASD. Finally, the proband's mtDNA presented five low-level heteroplasmic mtDNA variants that were absent in the mother, and two maternally inherited variants with increased heteroplasmic load. This study underlines the importance of a comprehensive assessment of the genomic background in carriers of large-effect variants, as penetrance modulation by additional interacting rare variants to might represent a widespread mechanism in neurodevelopmental disorders.

Brain Magnetic Resonance Findings in 117 Children with Autism Spectrum Disorder under 5 Years Old.

We examined the potential benefits of neuroimaging measurements across the first 5 years of life in detecting early comorbid or etiological signs of autism spectrum disorder (ASD). In particular, we analyzed the prevalence of neuroradiologic findings in routine magnetic resonance imaging (MRI) scans of a group of 117 ASD children younger than 5 years old. These data were compared to those reported in typically developing (TD) children. MRI findings in children with ASD were analyzed in relation to their cognitive level, severity of autistic symptoms, and the presence of electroencephalogram (EEG) abnormalities. The MRI was rated abnormal in 55% of children with ASD with a significant prevalence in the high-functioning subgroup compared to TD children. We report significant incidental findings of mega cisterna magna, ventricular anomalies and abnormal white matter signal intensity in ASD without significant associations between these MRI findings and EEG features. Based on these results we discuss the role that brain MRI may play in the diagnostic procedure of ASD.

An integrated analysis of rare CNV and exome variation in Autism Spectrum Disorder using the Infinium PsychArray.

Autism spectrum disorder (ASD) is a neurodevelopmental condition with a complex and heterogeneous genetic etiology. While a proportion of ASD risk is attributable to common variants, rare copy-number variants (CNVs) and protein-disrupting single-nucleotide variants (SNVs) have been shown to significantly contribute to ASD etiology. We analyzed a homogeneous cohort of 127 ASD Italian families genotyped with the Illumina PsychArray, to perform an integrated analysis of CNVs and SNVs and to assess their contribution to ASD risk. We observed a higher burden of rare CNVs, especially deletions, in ASD individuals versus unaffected controls. Furthermore, we identified a significant enrichment of rare CNVs intersecting ASD candidate genes reported in the SFARI database. Family-based analysis of rare SNVs genotyped by the PsychArray also indicated an increased transmission of rare SNV variants from heterozygous parents to probands, supporting a multigenic model of ASD risk with significant contributions of both variant types. Moreover, our study reinforced the evidence for a significant role of VPS13B, WWOX, CNTNAP2, RBFOX1, MACROD2, APBA2, PARK2, GPHN, and RNF113A genes in ASD susceptibility. Finally, we showed that the PsychArray, besides providing useful genotyping data in psychiatric disorders, is a valuable and cost-efficient tool for genic CNV detection, down to 10 kb.

ELMOD3-SH2D6 gene fusion as a possible co-star actor in autism spectrum disorder scenario.

Autism spectrum disorder (ASD) is a group of neurodevelopmental disorders characterized by high heritability. It is known that genetic factors contribute to ASD pathogenesis. In particular, copy number variants (CNVs) are involved in ASD susceptibility and can affect gene expression regulation. 2p11.2 microdeletions encompassing ELMOD3, CAPG and SH2D6 genes have been described in four unrelated ASD families. The present study revealed that this microdeletion is responsible for the production of a chimeric transcript generated from the fusion between ELMOD3 and SH2D6. The identified transcript showed significantly higher expression levels in subjects carrying the deletion compared to control subjects, suggesting that it is not subjected to nonsense-mediated decay and might encode for a chimeric protein. In conclusion, this study suggests the possible involvement of this gene fusion, together with the other previously identified variants, in ASD.

Contribution of ultrarare variants in mTOR pathway genes to sporadic focal epilepsies.

Objective: We investigated the contribution to sporadic focal epilepsies (FE) of ultrarare variants in genes coding for the components of complexes regulating mechanistic Target Of Rapamycin (mTOR)complex 1 (mTORC1). Methods: We collected genetic data of 121 Italian isolated FE cases and 512 controls by Whole Exome Sequencing (WES) and single-molecule Molecular Inversion Probes (smMIPs) targeting 10 genes of the GATOR1, GATOR2, and TSC complexes. We collapsed "qualifying" variants (ultrarare and predicted to be deleterious or loss of function) across the examined genes and sought to identify their enrichment in cases compared to controls. Results: We found eight qualifying variants in cases and nine in controls, demonstrating enrichment in FE patients (P = 0.006; exact unconditional test, one-tailed). Pathogenic variants were identified in DEPDC5 and TSC2, both major genes for Mendelian FE syndromes. Interpretation: Our findings support the contribution of ultrarare variants in genes in the mTOR pathway complexes GATOR and TSC to the risk of sporadic FE and a shared genetic basis between rare and common epilepsies. The identification of a monogenic etiology in isolated cases, most typically encountered in clinical practice, may offer to a broader community of patients the perspective of precision therapies directed by the underlying genetic cause.

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.

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.

Genetic variation in CHRNA7 and CHRFAM7A is associated with nicotine dependence and response to varenicline treatment.

The role of nicotinic acetylcholine receptors (nAChR) in nicotine dependence (ND) is well established; CHRNA7, encoding the α7 subunit, has a still uncertain role in ND, although it is implicated in a wide range of neuropsychiatric conditions. CHRFAM7A, a hybrid gene containing a partial duplication of CHRNA7, is possibly involved in modulating α7 nAChR function. The aim of this study was to investigate the role of CHRNA7 and CHRFAM7A genetic variants in ND and to test the hypothesis that α7 nAChR variation may modulate the efficacy of varenicline treatment in smoking cessation. We assessed CHRNA7 and CHRFAM7A copy number, CHRFAM7A exon 6 ∆2 bp polymorphism, and sequence variants in the CHRNA7 proximal promoter in an Italian sample of 408 treatment-seeking smokers. We conducted case-control and quantitative association analyses using two smoking measures (cigarettes per day, CPD, and Fagerström Test for Nicotine Dependence, FTND). Next, driven by the hypothesis that varenicline may exert some of its therapeutic effects through activation of α7 nAChRs, we restricted the analysis to a subgroup of 142 smokers who received varenicline treatment. The CHRNA7 promoter variant rs28531779 showed association with both smoking quantitative measures (FNTD p = 0.026, ß = 0.89, 95% CI 0.11-1.67; CPD p = 0.006, ß = 4.82 95% CI 1.42-8.22). Moreover, in the varenicline-treated subgroup we observed association of CHRFAM7A copy number with 6 months smoking abstinence (p = 0.035, OR = 3.18, 95% CI = 1.09-9.30). Thus, our study points to a possible role of genetic variation in CHRNA7 and CHRFAM7A in tobacco addiction mechanisms and response to varenicline treatment.

Integrated DNA methylation analysis identifies topographical and tumoral biomarkers in pilocytic astrocytomas.

Pilocytic astrocytoma (PA) is the most common glioma in pediatric patients and occurs in different locations. Chromosomal alterations are mostly located at chromosome 7q34 comprising the BRAF oncogene with consequent activation of the mitogen-activated protein kinase pathway. Although genetic and epigenetic alterations characterizing PA from different localizations have been reported, the role of epigenetic alterations in PA development is still not clear. The aim of this study was to investigate whether distinctive methylation patterns may define biologically relevant groups of PAs. Integrated DNA methylation analysis was performed on 20 PAs and 4 normal brain samples by Illumina Infinium HumanMethylation27 BeadChips. We identified distinct methylation profiles characterizing PAs from different locations (infratentorial vs supratentorial) and tumors with onset before and after 3 years of age. These results suggest that PA may be related to the specific brain site where the tumor arises from region-specific cells of origin. We identified and validated in silico the methylation alterations of some CpG islands. Furthermore, we evaluated the expression levels of selected differentially methylated genes and identified two biomarkers, one, IRX2, related to the tumor localization and the other, TOX2, as tumoral biomarker.

A genome-wide analysis in cluster headache points to neprilysin and PACAP receptor gene variants.

BACKGROUND: Cluster Headache (CH) is a severe primary headache, with a poorly understood pathophysiology. Complex genetic factors are likely to play a role in CH etiology; however, no confirmed gene associations have been identified. The aim of this study is to identify genetic variants influencing risk to CH and to explore the potential pathogenic mechanisms. METHODS: We have performed a genome-wide association study (GWAS) in a clinically well-defined cohort of 99 Italian patients with CH and in a control sample of 360 age-matched sigarette smoking healthy individuals, using the Infinium PsychArray (Illumina), which combines common highly-informative genome-wide tag SNPs and exonic SNPs. Genotype data were used to carry out a genome-wide single marker case-control association analysis using common SNPs, and a gene-based association analysis focussing on rare protein altering variants in 745 candidate genes with a putative role in CH. RESULTS: Although no single variant showed statistically significant association at the genome-wide threshold, we identified an interesting suggestive association (P = 9.1 × 10-6) with a common variant of the PACAP receptor gene (ADCYAP1R1). Furthermore, gene-based analysis provided significant evidence of association (P = 2.5 × 10-5) for a rare potentially damaging missense variant in the MME gene, encoding for the membrane metallo-endopeptidase neprilysin. CONCLUSIONS: Our study represents the first genome-wide association study of common SNPs and rare exonic variants influencing risk for CH. The most interesting results implicate ADCYAP1R1 and MME gene variants in CH susceptibility and point to a role for genes involved in pain processing. These findings provide new insights into the pathogenesis of CH that need further investigation and replication in larger CH samples.

Analysis of CHRNA7 rare variants in autism spectrum disorder susceptibility.

Chromosome 15q13.3 recurrent microdeletions are causally associated with a wide range of phenotypes, including autism spectrum disorder (ASD), seizures, intellectual disability, and other psychiatric conditions. Whether the reciprocal microduplication is pathogenic is less certain. CHRNA7, encoding for the alpha7 subunit of the neuronal nicotinic acetylcholine receptor, is considered the likely culprit gene in mediating neurological phenotypes in 15q13.3 deletion cases. To assess if CHRNA7 rare variants confer risk to ASD, we performed copy number variant analysis and Sanger sequencing of the CHRNA7 coding sequence in a sample of 135 ASD cases. Sequence variation in this gene remains largely unexplored, given the existence of a fusion gene, CHRFAM7A, which includes a nearly identical partial duplication of CHRNA7. Hence, attempts to sequence coding exons must distinguish between CHRNA7 and CHRFAM7A, making next-generation sequencing approaches unreliable for this purpose. A CHRNA7 microduplication was detected in a patient with autism and moderate cognitive impairment; while no rare damaging variants were identified in the coding region, we detected rare variants in the promoter region, previously described to functionally reduce transcription. This study represents the first sequence variant analysis of CHRNA7 in a sample of idiopathic autism.