Persistent DNA Methylation Changes across the First Year of Life and Prenatal NO2 Exposure in a Canadian Prospective Birth Study.

BACKGROUND: Evidence suggests that prenatal air pollution exposure alters DNA methylation (DNAm), which could go on to affect long-term health. It remains unclear whether DNAm alterations present at birth persist through early life. Identifying persistent DNAm changes would provide greater insight into the molecular mechanisms contributing to the association of prenatal air pollution exposure with atopic diseases. OBJECTIVES: This study investigated DNAm differences associated with prenatal nitrogen dioxide (NO2) exposure (a surrogate measure of traffic-related air pollution) at birth and 1 y of age and examined their role in atopic disease. We focused on regions showing persistent DNAm differences from birth to 1 y of age and regions uniquely associated with postnatal NO2 exposure. METHODS: Microarrays measured DNAm at birth and at 1 y of age for an atopy-enriched subset of Canadian Health Infant Longitudinal Development (CHILD) study participants. Individual and regional DNAm differences associated with prenatal NO2 (n=128) were identified, and their persistence at age 1 y were investigated using linear mixed effects models (n=124). Postnatal-specific DNAm differences (n=125) were isolated, and their association with NO2 in the first year of life was examined. Causal mediation investigated whether DNAm differences mediated associations between NO2 and age 1 y atopy or wheeze. Analyses were repeated using biological sex-stratified data. RESULTS: At birth (n=128), 18 regions of DNAm were associated with NO2, with several annotated to HOX genes. Some of these regions were specifically identified in males (n=73), but not females (n=55). The effect of prenatal NO2 across CpGs within altered regions persisted at 1 y of age. No significant mediation effects were identified. Sex-stratified analyses identified postnatal-specific DNAm alterations. DISCUSSION: Regional cord blood DNAm differences associated with prenatal NO2 persisted through at least the first year of life in CHILD participants. Some differences may represent sex-specific alterations, but replication in larger cohorts is needed. The early postnatal period remained a sensitive window to DNAm perturbations. https://doi.org/10.1289/EHP13034.

Wheeze trajectories: Determinants and outcomes in the CHILD Cohort Study.

BACKGROUND: Wheezing in early life is associated with asthma in adulthood; however, the determinants of wheezing trajectories and their associations with asthma and lung function in childhood remain poorly understood. OBJECTIVE: In the CHILD Cohort Study, we aimed to identify wheezing trajectories and examine the associations between these trajectories, risk factors, and clinical outcomes at age 5 years. METHODS: Wheeze data were collected at 8 time points from 3 months to 5 years of age. We used group-based trajectory models to derive wheeze trajectories among 3154 children. Associations with risk factors and clinical outcomes were analyzed by weighted regression models. RESULTS: We identified 4 trajectories: a never/infrequent trajectory, transient wheeze, intermediate-onset (preschool) wheeze, and persistent wheeze. Higher body mass index was a common risk factor for all wheeze trajectories compared with that in the never/infrequent group. The unique predictors for specific wheeze trajectories included male sex, lower respiratory tract infections, and day care attendance for transient wheeze; paternal history of asthma, atopic sensitization, and child genetic risk score of asthma for intermediate wheeze; and maternal asthma for persistent wheeze. Blood eosinophil counts were higher in children with the intermediate wheeze trajectory than in those children with the other trajectories at the ages of 1 and 5 years. All wheeze trajectories were associated with decreased lung function and increased risk of asthma at age 5 years. CONCLUSIONS: We identified 4 distinct trajectories in children from 3 months to 5 years of age, reflecting different phenotypes of early childhood wheeze. These trajectories were characterized by different biologic and physiologic traits and risk factors.

Missense variants in ATP1A3 and FXYD gene family are associated with childhood-onset schizophrenia.

Childhood-onset schizophrenia (COS) is a rare and severe form of schizophrenia defined as onset before age of 13. Here we report on two unrelated cases diagnosed with both COS and alternating hemiplegia of childhood (AHC), and for whom two distinct pathogenic de novo variants were identified in the ATP1A3 gene. ATP1A3 encodes the α-subunit of a neuron-specific ATP-dependent transmembrane sodium-potassium pump. Using whole exome sequencing data derived from a cohort of 17 unrelated COS cases, we also examined ATP1A3 and all of its interactors known to be expressed in the brain to establish if variants could be identified. This led to the identification of a third case with a possibly damaging missense mutation in ATP1A3 and three others cases with predicted pathogenic missense variants in the FXYD gene family (FXYD1, FXYD6, and FXYD6-FXYD2 readthrough). FXYD genes encode proteins that modulate the ATP-dependant pump function. This report is the first to identify variants in the same pathway for COS. Our COS study illustrates the interest of stratifying a complex condition according to the age of onset for the identification of deleterious variants. Whereas ATP1A3 is a replicated gene in rare neuropediatric diseases, this gene has previously been linked with COS in only one case report. The association with rare variants in FXYD gene family is novel and highlights the interest of exploring these genes in COS as well as in pediatric neurodevelopmental disorders.

Genetic architecture and adaptations of Nunavik Inuit.

The Canadian Inuit have a distinct population background that may entail particular implications for the health of its individuals. However, the number of genetic studies examining this Inuit population is limited, and much remains to be discovered in regard to its genetic characteristics. In this study, we generated whole-exome sequences and genomewide genotypes for 170 Nunavik Inuit, a small and isolated founder population of Canadian Arctic indigenous people. Our study revealed the genetic background of Nunavik Inuit to be distinct from any known present-day population. The majority of Nunavik Inuit show little evidence of gene flow from European or present-day Native American peoples, and Inuit living around Hudson Bay are genetically distinct from those around Ungava Bay. We also inferred that Nunavik Inuit have a small effective population size of 3,000 and likely split from Greenlandic Inuit ∼10.5 kya. Nunavik Inuit went through a bottleneck at approximately the same time and might have admixed with a population related to the Paleo-Eskimos. Our study highlights population-specific genomic signatures in coding regions that show adaptations unique to Nunavik Inuit, particularly in pathways involving fatty acid metabolism and cellular adhesion (CPNE7, ICAM5, STAT2, and RAF1). Subsequent analyses in selection footprints and the risk of intracranial aneurysms (IAs) in Nunavik Inuit revealed an exonic variant under weak negative selection to be significantly associated with IA (rs77470587; P = 4.6 × 10-8).

SMPD1 mutations, activity, and α-synuclein accumulation in Parkinson's disease.

BACKGROUND: SMPD1 (acid-sphingomyelinase) variants have been associated with Parkinson's disease in recent studies. The objective of this study was to further investigate the role of SMPD1 mutations in PD. METHODS: SMPD1 was sequenced in 3 cohorts (Israel Ashkenazi Jewish cohort, Montreal/Montpellier, and New York), including 1592 PD patients and 975 controls. Additional data were available for 10,709 Ashkenazi Jewish controls. Acid-sphingomyelinase activity was measured by a mass spectrometry-based assay in the New York cohort. α-Synuclein levels were measured in vitro following CRISPR/Cas9-mediated knockout and siRNA knockdown of SMPD1 in HeLa and BE(2)-M17 cells. Lysosomal localization of acid-sphingomyelinase with different mutations was studied, and in silico analysis of their effect on acid-sphingomyelinase structure was performed. RESULTS: SMPD1 mutations were associated with PD in the Ashkenazi Jewish cohort, as 1.4% of PD patients carried the p.L302P or p.fsP330 mutation, compared with 0.37% in 10,709 Ashkenazi Jewish controls (OR, 3.7; 95%CI, 1.6-8.2; P = 0.0025). In the Montreal/Montpellier cohort, the p.A487V variant was nominally associated with PD (1.5% versus 0.14%; P = 0.0065, not significant after correction for multiple comparisons). Among PD patients, reduced acid-sphingomyelinase activity was associated with a 3.5- to 5.8-year earlier onset of PD in the lowest quartile versus the highest quartile of acid-sphingomyelinase activity (P = 0.01-0.001). We further demonstrated that SMPD1 knockout and knockdown resulted in increased α-synuclein levels in HeLa and BE(2)-M17 dopaminergic cells and that the p.L302P and p.fsP330 mutations impair the traffic of acid-sphingomyelinase to the lysosome. CONCLUSIONS: Our results support an association between SMPD1 variants, acid-sphingomyelinase activity, and PD. Furthermore, they suggest that reduced acid-sphingomyelinase activity may lead to α-synuclein accumulation. © 2019 International Parkinson and Movement Disorder Society.

Exome sequencing of sporadic childhood-onset schizophrenia suggests the contribution of X-linked genes in males.

Childhood-onset schizophrenia (COS) is a rare and severe form of schizophrenia, defined as having an onset before the age of 13. The male COS cases have a slightly younger age of onset than female cases. They also present with a higher rate of comorbid developmental disorders. These sex differences are not explained by the frequency of chromosomal abnormalities, and the contribution of other forms of genetic variations remains unestablished. Using a whole-exome sequencing approach, we examined 12 COS trios where the unaffected parents had an affected male child. The sequencing data enabled us to test if the hemizygous variants, transmitted from the unaffected carrying mother, could mediate the phenotype (X-linked recessive inheritance model). Our results revealed that affected children have a significantly greater number of X-linked rare variants than their unaffected fathers. The variants identified in the male probands were mostly found in genes previously linked to other neuropsychiatric diseases like autism, intellectual disability, and epilepsy, including LUZP4, PCDH19, RPS6KA3, and OPHN1. The level of expression of the genes was assessed at different developmental periods in normal brain using the BrainSpan database. This approach revealed that some of them were expressed earlier in males than in females, consistent with the younger age of onset in male COS. In conclusion, this article suggests that X-linked genes might play a role in the pathophysiology of COS. Candidate genes detailed here could explain the higher level of comorbidities and the earlier age of onset observed in a subset of the male COS cases.

Absence of Mutation Enrichment for Genes Phylogenetically Conserved in the Olivocerebellar Motor Circuitry in a Cohort of Canadian Essential Tremor Cases.

Essential Tremor is a prevalent neurological disorder of unknown etiology. Studies suggest that genetic factors contribute to this pathology. To date, no causative mutations in a gene have been reproducibly reported. All three structures of the olivocerebellar motor circuitry have been linked to Essential Tremor. We postulated that genes enriched for their expression in the olivocerebellar circuitry would be more susceptible to harbor mutations in Essential Tremor patients. A list of 11 candidate genes, enriched for their expression in the olivocerebellar circuitry, was assessed for their variation spectrum and frequency in a cohort of Canadian Essential Tremor cases. Our results from this list of 11 candidate genes do not support an association for Essential Tremor in our cohort of Canadian cases. The heterogenic nature of ET and modest size of the cohort used in this study are two confounding factors that could explain these results.

Full sequencing and haplotype analysis of MAPT in Parkinson's disease and rapid eye movement sleep behavior disorder.

BACKGROUND: MAPT haplotypes are associated with PD, but their association with rapid eye movement sleep behavior disorder is unclear. OBJECTIVE: To study the role of MAPT variants in rapid eye movement sleep behavior disorder. METHODS: Two cohorts were included: (A) PD (n = 600), rapid eye movement sleep behavior disorder (n = 613) patients, and controls (n = 981); (B) dementia with Lewy bodies patients with rapid eye movement sleep behavior disorder (n = 271) and controls (n = 950). MAPT-associated variants and the entire coding sequence of MAPT were analyzed. Age-, sex-, and ethnicity-adjusted analyses were performed to examine the association between MAPT, PD, and rapid eye movement sleep behavior disorder. RESULTS: MAPT-H2 variants were associated with PD (odds ratios: 0.62-0.65; P = 0.010-0.019), but not with rapid eye movement sleep behavior disorder. In PD, the H1 haplotype odds ratio was 1.60 (95% confidence interval: 1.12-2.28; P = 0.009), and the H2 odds ratio was 0.68 (95% confidence interval: 0.48-0.96; P = 0.03). The H2/H1 haplotypes were not associated with rapid eye movement sleep behavior disorder. CONCLUSIONS: Our results confirm the protective effect of the MAPT-H2 haplotype in PD, and define its components. Furthermore, our results suggest that MAPT does not play a major role in rapid eye movement sleep behavior disorder, emphasizing different genetic background than in PD in this locus. © 2018 International Parkinson and Movement Disorder Society.

A rare variant in MLKL confers susceptibility to ApoE ɛ4-negative Alzheimer's disease in Hong Kong Chinese population.

Alzheimer's disease (AD) is the most common neurodegenerative disorders in the elderly. To identify rare genetic factors other than apolipoprotein E ɛ4 allele (ApoE ɛ4) contributing to the pathogenesis of late-onset AD (LOAD), we conducted a whole-exome analysis of 246 ApoE ɛ4-negative LOAD cases and 172 matched controls in Hong Kong Chinese population. LOAD patients showed a significantly higher burden of rare loss-of-function variants in genes related to immune function than healthy controls. Among the genes involved in immune function, we identified a rare stop-gain variant (p.Q48X) in mixed lineage kinase domain like pseudokinase (MLKL) gene present exclusively in 6 LOAD cases. MLKL is expressed in neurons, and the its expression levels in the p.Q48X carriers were significantly lower than that in age-matched wild-type controls. The ratio of Aß42 to Aß40 significantly increased in MLKL knockdown cells compared to scramble controls. MLKL loss-of-function mutation might contribute to late-onset ApoE ɛ4-negative AD in the Hong Kong Chinese population.

Genome-wide association analysis identifies new candidate risk loci for familial intracranial aneurysm in the French-Canadian population.

Intracranial Aneurysm (IA) is a common disease with a worldwide prevalence of 1-3%. In the French-Canadian (FC) population, where there is an important founder effect, the incidence of IA is higher and is frequently seen in families. In this study, we genotyped a cohort of 257 mostly familial FC IA patients and 1,992 FC controls using the Illumina NeuroX SNP-chip. The most strongly associated loci were tested in 34 Inuit IA families and in 32 FC IA patients and 106 FC controls that had been exome sequenced (WES). After imputation, one locus at 3p14.2 (FHIT, rs1554600, p = 4.66 × 10-9) reached a genome-wide significant level of association and a subsequent validation in Nunavik Inuit cohort further confirmed the significance of the FHIT variant association (rs780365, FBAT-O, p = 0.002839). Additionally, among the other promising loci (p < 5 × 10-6), the one at 3q13.2 (rs78125721, p = 4.77 × 10-7), which encompasses CCDC80, also showed an increased mutation burden in the WES data (CCDC80, SKAT-O, p = 0.0005). In this study, we identified two new potential IA loci in the FC population: FHIT, which is significantly associated with hypertensive IA, and CCDC80, which has potential genetic and functional relevance to IA pathogenesis, providing evidence on the additional risk loci for familial IA. We also replicated the previous IA GWAS risk locus 18q11.2, and suggested a potential locus at 8p23.1 that warrants further study.

Teneurin transmembrane protein 4 is not a cause for essential tremor in a Canadian population.

INTRODUCTION: Mutations in teneurin transmembrane protein 4 were reported to be a risk factor for essential tremor, but the relevance of this across different population remains to be examined. The aim of this study was to determine the frequency and spectrum of variations in teneurin transmembrane protein 4 in a cohort of Canadian essential tremor cases. METHODS: The coding portion of teneurin transmembrane protein 4 was sequenced in 269 unrelated essential tremor cases and 288 matched control individuals using a targeted and high-throughput sequencing approach. RESULTS: A total of 157 single nucleotide variations were identified, and from these 99 were a missense or nonsense mutation. A total of 68 cases were carriers of ≥1 rare missense or nonsense mutations, and 39 control individuals were carriers of the same types of variations. Gene-based association tests were used to jointly analyze the single nucleotide variations. CONCLUSIONS: Our results do not support a positive association between teneurin transmembrane protein 4 and the Canadian population. © 2017 International Parkinson and Movement Disorder Society.

No rare deleterious variants from STK32B, PPARGC1A, and CTNNA3 are associated with essential tremor.

OBJECTIVE: To assess the contribution of variants in STK32B, PPARGC1A, and CTNNA3 as essential tremor (ET) predisposing factors following their association in a 2-stage genome-wide association study (GWAS). METHODS: The coding regions of these genes was examined for the presence of rare variants using two approaches: (1) Looking at whole-exome and whole-genome sequencing data of 14 autosomal dominant multiplex ET families. (2) Conducting a targeted massive parallel sequencing to examine the three genes in cohorts of 269 ET cases and 287 control individuals. The cumulative impact of rare variants was assessed using SKAT-O analyses using (1) all variants, (2) only rare variants, and (3) only the rare variants altering the mRNA. RESULTS: Thirty-four variants were identified. No difference emerged regarding the distributions of individual variants (or gene) between cases and controls. CONCLUSION: No rare exonic variants further validated one of these genes as a risk factor for ET. The recent GWAS offers promising avenues, but the genetic heterogeneity of ET is nonetheless challenging for the validation of risk factors, and ultimately larger cohorts of cases should help to overcome this task.

Paternal Age Explains a Major Portion of De Novo Germline Mutation Rate Variability in Healthy Individuals.

De novo mutations (DNM) are an important source of rare variants and are increasingly being linked to the development of many diseases. Recently, the paternal age effect has been the focus of a number of studies that attempt to explain the observation that increasing paternal age increases the risk for a number of diseases. Using disease-free familial quartets we show that there is a strong positive correlation between paternal age and germline DNM in healthy subjects. We also observed that germline CNVs do not follow the same trend, suggesting a different mechanism. Finally, we observed that DNM were not evenly distributed across the genome, which adds support to the existence of DNM hotspots.

RNF213 Is Associated with Intracranial Aneurysms in the French-Canadian Population.

Intracranial aneurysms (IAs) are the result of focal weakness in the artery wall and have a complex genetic makeup. To date, genome-wide association and sequencing studies have had limited success in identifying IA risk factors. Distinct populations, such as the French-Canadian (FC) population, have increased IA prevalence. In our study, we used exome sequencing to prioritize risk variants in a discovery cohort of six FC families affected by IA, and the analysis revealed an increased variation burden for ring finger protein 213 (RNF213). We resequenced RNF213 in a larger FC validation cohort, and association tests on further identified variants supported our findings (SKAT-O, p = 0.006). RNF213 belongs to the AAA+ protein family, and two variants (p.Arg2438Cys and p.Ala2826Thr) unique to affected FC individuals were found to have increased ATPase activity, which could lead to increased risk of IA by elevating angiogenic activities. Common SNPs in RNF213 were also extracted from the NeuroX SNP-chip genotype data, comprising 257 FC IA-affected and 1,988 control individuals. We discovered that the non-ancestral allele of rs6565666 was significantly associated with the affected individuals (p = 0.03), and it appeared as though the frequency of the risk allele had changed through genetic drift. Although RNF213 is a risk factor for moyamoya disease in East Asians, we demonstrated that it might also be a risk factor for IA in the FC population. It therefore appears that the function of RNF213 can be differently altered to predispose distinct populations to dissimilar neurovascular conditions, highlighting the importance of a population's background in genetic studies of heterogeneous disease.

GBA p.T369M substitution in Parkinson disease: Polymorphism or association? A meta-analysis.

The lysosomal enzyme glucocerebrosidase (GCase), encoded by GBA, has an important role in Parkinson disease (PD). GBA mutation carriers have an increased risk for PD, earlier age at onset, faster progression, and various nonmotor symptoms including cognitive decline, REM sleep behavior disorder, hyposmia, and autonomic dysfunction.(1) Furthermore, GCase enzymatic activity is reduced in the peripheral blood(2) and brain(3) of noncarrier, sporadic PD patients. Biallelic GBA mutations, which have been classified as "severe" or "mild," may cause Gaucher disease (GD), a lysosomal storage disorder. Mild mutations may lead to GD type 1, and 2 severe mutations result in neuronopathic GD (type 2 and type 3).(4) There are 2 GBA variants, p.E326K and p.T369M, which do not cause GD in homozygous carriers, but may modify GCase activity and GD phenotype. It is now clear that p.E326K is a risk factor for PD,(5) but whether p.T369M is associated with PD is still controversial. In some studies, the p.T369M substitution was associated with PD,(6) while in others it had similar or increased frequency in controls. Of interest, it was recently demonstrated that the GBA p.T369M substitution was associated with reduced enzymatic activity in patients with PD and controls compared with that in noncarriers (7.64 vs 11.93 µmol/L/h, p < 0.001).(2) Of interest, it was even lower than the average enzymatic activity of the p.E326K substitution, which was 9.81 µmol/L/h. Because clinical trials on GBA-associated PD are ongoing, and because treatment specifically targeting GBA is likely to be available in the future, it is important to determine whether the GBA p.T369M substitution is associated with PD.

Identification of genetic variants of LGI1 and RTN4R (NgR1) linked to schizophrenia that are defective in NgR1-LGI1 signaling.

BACKGROUND: The protein NgR1 is encoded by RTN4R, a gene linked to schizophrenia. We previously reported NgR1 as receptor for the epilepsy-linked protein LGI1. NgR1 regulates synapse number and synaptic plasticity, whereas LGI1 antagonizes NgR1 signaling and promotes synapse formation. Impairments in synapse formation are common in neurological disease and we hypothesized that an LGI1-NgR1 signaling pathway may contribute to the development of schizophrenia. METHODS: We screened two unrelated schizophrenic populations for variants in RTN4R and LGI1 using whole exome sequencing and Sanger sequencing. We tested the ability of LGI1 to bind rare coding variants of NgR1 using a cell surface binding assays and the signaling ability of NgR1 using COS7 cell-spreading assays. RESULTS: We observed a previously reported rare coding variant in RTN4R (c.1195C>T, pR399W). We report the first LGI1 mutations to be identified in individuals with schizophrenia. Three different LGI1 mutations were found, two missense mutations (c.205G>A, p.V69I) and (c.313G>A, V105M), and an intronic variant (g.897T>C) that likely leads to a protein truncation. We found NgR1(R119W) and NgR1(277C) have a reduced ability to bind LGI1 in a cell surface binding assay. COS7 cell-spreading assays reveal that NgR1 mutants are impaired in their ability to mediate RhoA activation. CONCLUSION: Variants in NgR1 and LGI1 may be associated with schizophrenia and variants in NgR1 found in schizophrenic patients have impaired LGI1-NgR1 signaling. Impaired LGI1-NgR1 signaling may contribute to disease progression.

Analysis of DNAJC13 mutations in French-Canadian/French cohort of Parkinson's disease.

DNAJC13 mutations have been suggested to cause Parkinson's disease (PD), yet subsequent studies reported conflicting results on this association. In the present study, we sequenced the coding region of DNAJC13 in a French-Canadian/French cohort of 528 PD patients and 692 controls. A total of 62 (11.7%) carriers of rare DNAJC13 variants were identified among the PD patients compared with 82 (11.8%) among controls (p = 1.0). Two variants that were previously suggested to be associated with PD, p.R1516H and p.L2170W, were identified with similar directions of association as previously reported. The p.R1516H was found in 2 (0.4%) patients versus 6 (0.9%, nonsignificant) controls and the p.L2170W variant was found in 9 (1.7%) patients and 5 (0.7%, nonsignificant) controls. Meta-analysis with previous reports resulted in odds ratios of 0.32 (95% confidence interval = 0.15-0.68, p = 0.0037) and 2.68 (95% confidence interval = 1.32-5.42, p = 0.007), respectively. Our results provide some support for the possibility that specific DNAJC13 variants may play a minor role in PD susceptibility, although studies in additional populations are necessary.

The role of the melanoma gene MC1R in Parkinson disease and REM sleep behavior disorder.

The MC1R gene, suggested to be involved in Parkinson disease (PD) and melanoma, was sequenced in PD patients (n = 539) and controls (n = 265) from New York, and PD patients (n = 551), rapid eye movement sleep behavior disorder (RBD) patients (n = 351), and controls (n = 956) of European ancestry. Sixty-eight MC1R variants were identified, including 7 common variants with frequency > 0.01. None of the common variants was associated with PD or RBD in the different regression models. In a meta-analysis with fixed-effect model, the p.R160W variant was associated with an increased risk for PD (odds ratio = 1.22, 95% confidence interval = 1.02-1.47, p = 0.03) but with significant heterogeneity (p = 0.048). Removing one study that introduced the heterogeneity resulted in nonsignificant association (odds ratio = 1.11, 95% confidence interval, 0.92-1.35, p = 0.27, heterogeneity p = 0.57). Rare variants had similar frequencies in patients and controls (10.54% and 10.15%, respectively, p = 0.75), and no cumulative effect of carrying more than one MC1R variant was found. The present study does not support a role for the MC1R p.R160W and other variants in susceptibility for PD or RBD.

De novo variants in sporadic cases of childhood onset schizophrenia.

Childhood-onset schizophrenia (COS), defined by the onset of illness before age 13 years, is a rare severe neurodevelopmental disorder of unknown etiology. Recently, sequencing studies have identified rare, potentially causative de novo variants in sporadic cases of adult-onset schizophrenia and autism. In this study, we performed exome sequencing of 17 COS trios in order to test whether de novo variants could contribute to this disease. We identified 20 de novo variants in 17 COS probands, which is consistent with the de novo mutation rate reported in the adult form of the disease. Interestingly, the missense de novo variants in COS have a high likelihood for pathogenicity and were enriched for genes that are less tolerant to variants. Among the genes found disrupted in our study, SEZ6, RYR2, GPR153, GTF2IRD1, TTBK1 and ITGA6 have been previously linked to neuronal function or to psychiatric disorders, and thus may be considered as COS candidate genes.