Direct Measure of the De Novo Mutation Rate in Autism and Schizophrenia Cohorts

The role of de novo mutations (DNMs) in common diseases remains largely unknown. Nonetheless, the rate of de novo deleterious mutations and the strength of selection against de novo mutations are critical to understanding the genetic architecture of a disease. Discovery of high-impact DNMs requires substantial high-resolution interrogation of partial or complete genomes of families via resequencing. We hypothesized that deleterious DNMs may play a role in cases of autism spectrum disorders (ASD) and schizophrenia (SCZ), two etiologically heterogeneous disorders with significantly reduced reproductive fitness. We present a direct measure of the de novo mutation rate (μ) and selective constraints from DNMs estimated from a deep resequencing data set generated from a large cohort of ASD and SCZ cases (n = 285) and population control individuals (n = 285) with available parental DNA. A survey of ∼430 Mb of DNA from 401 synapse-expressed genes across all cases and 25 Mb of DNA in controls found 28 candidate DNMs, 13 of which were cell line artifacts. Our calculated direct neutral mutation rate (1.36 × 10⁻⁸) is similar to previous indirect estimates, but we observed a significant excess of potentially deleterious DNMs in ASD and SCZ individuals. Our results emphasize the importance of DNMs as genetic mechanisms in ASD and SCZ and the limitations of using DNA from archived cell lines to identify functional variants.     

A Higher Mutational Burden in Females Supports a “Female Protective Model” in Neurodevelopmental Disorders

Increased male prevalence has been repeatedly reported in several neurodevelopmental disorders (NDs), leading to the concept of a “female protective model.” We investigated the molecular basis of this sex-based difference in liability and demonstrated an excess of deleterious autosomal copy-number variants (CNVs) in females compared to males (odds ratio [OR] = 1.46, p = 8 × 10⁻¹⁰) in a cohort of 15,585 probands ascertained for NDs. In an independent autism spectrum disorder (ASD) cohort of 762 families, we found a 3-fold increase in deleterious autosomal CNVs (p = 7 × 10⁻⁴) and an excess of private deleterious single-nucleotide variants (SNVs) in female compared to male probands (OR = 1.34, p = 0.03). We also showed that the deleteriousness of autosomal SNVs was significantly higher in female probands (p = 0.0006). A similar bias was observed in parents of probands ascertained for NDs. Deleterious CNVs (>400 kb) were maternally inherited more often (up to 64%, p = 10⁻¹⁵) than small CNVs < 400 kb (OR = 1.45, p = 0.0003). In the ASD cohort, increased maternal transmission was also observed for deleterious CNVs and SNVs. Although ASD females showed higher mutational burden and lower cognition, the excess mutational burden remained, even after adjustment for those cognitive differences. These results strongly suggest that females have an increased etiological burden unlinked to rare deleterious variants on the X chromosome. Carefully phenotyped and genotyped cohorts will be required for identifying the symptoms, which show gender-specific liability to mutational burden.     

Next-Generation Sequencing and In Vitro Expression Study of ADAMTS13 Single Nucleotide Variants in Deep Vein Thrombosis

BackgroundDeep vein thrombosis (DVT) genetic predisposition is partially known.ObjectivesThis study aimed at assessing the functional impact of nine ADAMTS13 single nucleotide variants (SNVs) previously reported to be associated as a group with DVT in a burden test and the individual association of selected variants with DVT risk in two replication studies.MethodsWild-type and mutant recombinant ADAMTS13 were transiently expressed in HEK293 cells. Antigen and activity of recombinant ADAMTS13 were measured by ELISA and FRETS-VWF73 assays, respectively. The replication studies were performed in an Italian case-control study (Milan study; 298/298 patients/controls) using a next-generation sequencing approach and in a Dutch case-control study (MEGA study; 4306/4887 patients/controls) by TaqMan assays.ResultsIn vitro results showed reduced ADAMTS13 activity for three SNVs (p.Val154Ile [15%; 95% confidence interval [CI] 14–16], p.Asp187His [19%; 95%[CI] 17–21], p.Arg421Cys [24%; 95%[CI] 22–26]) similar to reduced plasma ADAMTS13 levels of patients carriers for these SNVs. Therefore these three SNVs were interrogated for risk association. The first replication study identified 3 heterozygous carriers (2 cases, 1 control) of p.Arg421Cys (odds ratio [OR] 2, 95%[CI] 0.18–22.25). The second replication study identified 2 heterozygous carriers (1 case, 1 control) of p.Asp187His ([OR] 1.14, 95%[CI] 0.07–18.15) and 10 heterozygous carriers (4 cases, 6 controls) of p.Arg421Cys ([OR] 0.76, 95%[CI] 0.21–2.68).ConclusionsThree SNVs (p.Val154Ile, p.Asp187His and p.Arg421Cys) showed reduced ex vivo and in vitro ADAMTS13 levels. However, the low frequency of these variants makes it difficult to confirm their association with DVT.     

Joint Linkage and Association Analysis with Exome Sequence Data Implicates SLC25A40 in Hypertriglyceridemia

Hypertriglyceridemia (HTG) is a heritable risk factor for cardiovascular disease. Investigating the genetics of HTG may identify new drug targets. There are ∼35 known single-nucleotide variants (SNVs) that explain only ∼10% of variation in triglyceride (TG) level. Because of the genetic heterogeneity of HTG, a family study design is optimal for identification of rare genetic variants with large effect size because the same mutation can be observed in many relatives and cosegregation with TG can be tested. We considered HTG in a five-generation family of European American descent (n = 121), ascertained for familial combined hyperlipidemia. By using Bayesian Markov chain Monte Carlo joint oligogenic linkage and association analysis, we detected linkage to chromosomes 7 and 17. Whole-exome sequence data revealed shared, highly conserved, private missense SNVs in both SLC25A40 on chr7 and PLD2 on chr17. Jointly, these SNVs explained 49% of the genetic variance in TG; however, only the SLC25A40 SNV was significantly associated with TG (p = 0.0001). This SNV, c.374A>G, causes a highly disruptive p.Tyr125Cys substitution just outside the second helical transmembrane region of the SLC25A40 inner mitochondrial membrane transport protein. Whole-gene testing in subjects from the Exome Sequencing Project confirmed the association between TG and SLC25A40 rare, highly conserved, coding variants (p = 0.03). These results suggest a previously undescribed pathway for HTG and illustrate the power of large pedigrees in the search for rare, causal variants.     

Resequencing and Association Analysis of CLN8 with Autism Spectrum Disorder in a Japanese Population

Rare variations contribute substantially to autism spectrum disorder (ASD) liability. We recently performed whole-exome sequencing in two families with affected siblings and then carried out a follow-up study and identified ceroid-lipofuscinosis neuronal 8 (epilepsy, progressive with mental retardation) (CLN8) as a potential genetic risk factor for ASD. To further investigate the role of CLN8 in the genetic etiology of ASD, we performed resequencing and association analysis of CLN8 with ASD in a Japanese population. Resequencing the CLN8 coding region in 256 ASD patients identified five rare missense variations: g.1719291G>A (R24H), rs201670636 (F39L), rs116605307 (R97H), rs143701028 (T108M) and rs138581191 (N152S). These variations were genotyped in 568 patients (including the resequenced 256 patients) and 1017 controls. However, no significant association between these variations and ASD was identified. This study does not support a contribution of rare missense CLN8 variations to ASD susceptibility in the Japanese population.     

Population-Based Resequencing of APOA1 in 10,330 Individuals: Spectrum of Genetic Variation,Phenotype, and Comparison with Extreme Phenotype Approach

Rare genetic variants, identified by in-detail resequencing of loci, may contribute to complex traits. We used the apolipoprotein A-I gene (APOA1), a major high-density lipoprotein (HDL) gene, and population-based resequencing to determine the spectrum of genetic variants, the phenotypic characteristics of these variants, and how these results compared with results based on resequencing only the extremes of the apolipoprotein A-I (apoA-I) distribution. First, we resequenced APOA1 in 10,330 population-based participants in the Copenhagen City Heart Study. The spectrum and distribution of genetic variants was determined as a function of the number of individuals resequenced. Second, apoA-I and HDL cholesterol phenotypes were determined for nonsynonymous (NS) and synonymous (S) variants and were validated in the Copenhagen General Population Study (n = 45,239). Third, observed phenotypes were compared with those predicted using an extreme phenotype approach based on the apoA-I distribution. Our results are as follows: First, population-based resequencing of APOA1 identified 40 variants of which only 7 (18%) had minor allele frequencies >1%, and most were exceedingly rare. Second, 0.27% of individuals in the general population were heterozygous for NS variants which were associated with substantial reductions in apoA-I (up to 39 mg/dL) and/or HDL cholesterol (up to 0.9 mmol/L) and, surprisingly, 0.41% were heterozygous for variants predisposing to amyloidosis. NS variants associated with a hazard ratio of 1.72 (1.09–2.70) for myocardial infarction (MI), largely driven by A164S, a variant not associated with apoA-I or HDL cholesterol levels. Third, using the extreme apoA-I phenotype approach, NS variants correctly predicted the apoA-I phenotype observed in the population-based resequencing. However, using the extreme approach, between 79% (screening 0–1^st percentile) and 21% (screening 0–20^th percentile) of all variants were not identified; among these were variants previously associated with amyloidosis. Population-based resequencing of APOA1 identified a majority of rare NS variants associated with reduced apoA-1 and HDL cholesterol levels and/or predisposing to amyloidosis. In addition, NS variants associated with increased risk of MI.     

No Evidence for Association of Autism with Rare Heterozygous Point Mutations in Contactin-Associated Protein-Like 2 (CNTNAP2), or in Other Contactin-Associated Proteins or Contactins

Contactins and Contactin-Associated Proteins, and Contactin-Associated Protein-Like 2 (CNTNAP2) in particular, have been widely cited as autism risk genes based on findings from homozygosity mapping, molecular cytogenetics, copy number variation analyses, and both common and rare single nucleotide association studies. However, data specifically with regard to the contribution of heterozygous single nucleotide variants (SNVs) have been inconsistent. In an effort to clarify the role of rare point mutations in CNTNAP2 and related gene families, we have conducted targeted next-generation sequencing and evaluated existing sequence data in cohorts totaling 2704 cases and 2747 controls. We find no evidence for statistically significant association of rare heterozygous mutations in any of the CNTN or CNTNAP genes, including CNTNAP2, placing marked limits on the scale of their plausible contribution to risk.     

Genetic and functional analyses of SHANK2 mutations suggest a multiple hit model of autism spectrum disorders

Autism spectrum disorders (ASD) are a heterogeneous group of neurodevelopmental disorders with a complex inheritance pattern. While many rare variants in synaptic proteins have been identified in patients with ASD, little is known about their effects at the synapse and their interactions with other genetic variations. Here, following the discovery of two de novo SHANK2 deletions by the Autism Genome Project, we identified a novel 421 kb de novo SHANK2 deletion in a patient with autism. We then sequenced SHANK2 in 455 patients with ASD and 431 controls and integrated these results with those reported by Berkel et al. 2010 (n = 396 patients and n = 659 controls). We observed a significant enrichment of variants affecting conserved amino acids in 29 of 851 (3.4%) patients and in 16 of 1,090 (1.5%) controls (P = 0.004, OR = 2.37, 95% CI = 1.23–4.70). In neuronal cell cultures, the variants identified in patients were associated with a reduced synaptic density at dendrites compared to the variants only detected in controls (P = 0.0013). Interestingly, the three patients with de novo SHANK2 deletions also carried inherited CNVs at 15q11–q13 previously associated with neuropsychiatric disorders. In two cases, the nicotinic receptor CHRNA7 was duplicated and in one case the synaptic translation repressor CYFIP1 was deleted. These results strengthen the role of synaptic gene dysfunction in ASD but also highlight the presence of putative modifier genes, which is in keeping with the “multiple hit model” for ASD. A better knowledge of these genetic interactions will be necessary to understand the complex inheritance pattern of ASD.     

Truncating mutations in NRXN2 and NRXN1 in autism spectrum disorders and schizophrenia

Growing genetic evidence is converging in favor of common pathogenic mechanisms for autism spectrum disorders (ASD), intellectual disability (ID or mental retardation) and schizophrenia (SCZ), three neurodevelopmental disorders affecting cognition and behavior. Copy number variations and deleterious mutations in synaptic organizing proteins including NRXN1 have been associated with these neurodevelopmental disorders, but no such associations have been reported for NRXN2 or NRXN3. From resequencing the three neurexin genes in individuals affected by ASD (n = 142), SCZ (n = 143) or non-syndromic ID (n = 94), we identified a truncating mutation in NRXN2 in a patient with ASD inherited from a father with severe language delay and family history of SCZ. We also identified a de novo truncating mutation in NRXN1 in a patient with SCZ, and other potential pathogenic ASD mutations. These truncating mutations result in proteins that fail to promote synaptic differentiation in neuron coculture and fail to bind either of the established postsynaptic binding partners LRRTM2 or NLGN2 in cell binding assays. Our findings link NRXN2 disruption to the pathogenesis of ASD for the first time and further strengthen the involvement of NRXN1 in SCZ, supporting the notion of a common genetic mechanism in these disorders.     

A Cytogenetic Abnormality and Rare Coding Variants Identify ABCA13 as a Candidate Gene in Schizophrenia, Bipolar Disorder, and Depression

Schizophrenia and bipolar disorder are leading causes of morbidity across all populations, with heritability estimates of ∼80% indicating a substantial genetic component. Population genetics and genome-wide association studies suggest an overlap of genetic risk factors between these illnesses but it is unclear how this genetic component is divided between common gene polymorphisms, rare genomic copy number variants, and rare gene sequence mutations. We report evidence that the lipid transporter gene ABCA13 is a susceptibility factor for both schizophrenia and bipolar disorder. After the initial discovery of its disruption by a chromosome abnormality in a person with schizophrenia, we resequenced ABCA13 exons in 100 cases with schizophrenia and 100 controls. Multiple rare coding variants were identified including one nonsense and nine missense mutations and compound heterozygosity/homozygosity in six cases. Variants were genotyped in additional schizophrenia, bipolar, depression (n > 1600), and control (n > 950) cohorts and the frequency of all rare variants combined was greater than controls in schizophrenia (OR = 1.93, p = 0.0057) and bipolar disorder (OR = 2.71, p = 0.00007). The population attributable risk of these mutations was 2.2% for schizophrenia and 4.0% for bipolar disorder. In a study of 21 families of mutation carriers, we genotyped affected and unaffected relatives and found significant linkage (LOD = 4.3) of rare variants with a phenotype including schizophrenia, bipolar disorder, and major depression. These data identify a candidate gene, highlight the genetic overlap between schizophrenia, bipolar disorder, and depression, and suggest that rare coding variants may contribute significantly to risk of these disorders.     

Identification of rare DNA variants in mitochondrial disorders with improved array-based sequencing

A common goal in the discovery of rare functional DNA variants via medical resequencing is to incur a relatively lower proportion of false positive base-calls. We developed a novel statistical method for resequencing arrays (SRMA, sequence robust multi-array analysis) to increase the accuracy of detecting rare variants and reduce the costs in subsequent sequence verifications required in medical applications. SRMA includes single and multi-array analysis and accounts for technical variables as well as the possibility of both low- and high-frequency genomic variation. The confidence of each base-call was ranked using two quality measures. In comparison to Sanger capillary sequencing, we achieved a false discovery rate of 2% (false positive rate 1.2 × 10⁻⁵, false negative rate 5%), which is similar to automated second-generation sequencing technologies. Applied to the analysis of 39 nuclear candidate genes in disorders of mitochondrial DNA (mtDNA) maintenance, we confirmed mutations in the DNA polymerase gamma POLG in positive control cases, and identified novel rare variants in previously undiagnosed cases in the mitochondrial topoisomerase TOP1MT, the mismatch repair enzyme MUTYH, and the apurinic-apyrimidinic endonuclease APEX2. Some patients carried rare heterozygous variants in several functionally interacting genes, which could indicate synergistic genetic effects in these clinically similar disorders.     

Whole Exome Sequencing of Distant Relatives in Multiplex Families Implicates Rare Variants in Candidate Genes for Oral Clefts

A dozen genes/regions have been confirmed as genetic risk factors for oral clefts in human association and linkage studies, and animal models argue even more genes may be involved. Genomic sequencing studies should identify specific causal variants and may reveal additional genes as influencing risk to oral clefts, which have a complex and heterogeneous etiology. We conducted a whole exome sequencing (WES) study to search for potentially causal variants using affected relatives drawn from multiplex cleft families. Two or three affected second, third, and higher degree relatives from 55 multiplex families were sequenced. We examined rare single nucleotide variants (SNVs) shared by affected relatives in 348 recognized candidate genes. Exact probabilities that affected relatives would share these rare variants were calculated, given pedigree structures, and corrected for the number of variants tested. Five novel and potentially damaging SNVs shared by affected distant relatives were found and confirmed by Sanger sequencing. One damaging SNV in CDH1, shared by three affected second cousins from a single family, attained statistical significance (P = 0.02 after correcting for multiple tests). Family-based designs such as the one used in this WES study offer important advantages for identifying genes likely to be causing complex and heterogeneous disorders.     

Whole-Exome Sequencing in a South American Cohort Links ALDH1A3, FOXN1 and Retinoic Acid Regulation Pathways to Autism Spectrum Disorders

Autism spectrum disorders (ASDs) are a range of complex neurodevelopmental conditions principally characterized by dysfunctions linked to mental development. Previous studies have shown that there are more than 1000 genes likely involved in ASD, expressed mainly in brain and highly interconnected among them. We applied whole exome sequencing in Colombian—South American trios. Two missense novel SNVs were found in the same child: ALDH1A3 (RefSeq NM_000693: c.1514T>C (p.I505T)) and FOXN1 (RefSeq NM_003593: c.146C>T (p.S49L)). Gene expression studies reveal that Aldh1a3 and Foxn1 are expressed in ~E13.5 mouse embryonic brain, as well as in adult piriform cortex (PC; ~P30). Conserved Retinoic Acid Response Elements (RAREs) upstream of human ALDH1A3 and FOXN1 and in mouse Aldh1a3 and Foxn1 genes were revealed using bioinformatic approximation. Chromatin immunoprecipitation (ChIP) assay using Retinoid Acid Receptor B (Rarb) as the immunoprecipitation target suggests RA regulation of Aldh1a3 and Foxn1 in mice. Our results frame a possible link of RA regulation in brain to ASD etiology, and a feasible non-additive effect of two apparently unrelated variants in ALDH1A3 and FOXN1 recognizing that every result given by next generation sequencing should be cautiously analyzed, as it might be an incidental finding.     

Lipoprotein lipase gene sequencing and plasma lipid profile

Lipoprotein lipase (LPL) plays a crucial role in lipid metabolism by hydrolyzing triglyceride (TG)-rich particles and affecting HDL cholesterol (HDL-C) levels. In this study, the entire LPL gene plus flanking regions were resequenced in individuals with extreme HDL-C/TG levels (n = 95), selected from a population-based sample of 623 US non-Hispanic White (NHW) individuals. A total of 176 sequencing variants were identified, including 28 novel variants. A subset of 64 variants [common tag single nucleotide polymorphisms (tagSNP) and selected rare variants] were genotyped in the total sample, followed by association analyses with major lipid traits. A gene-based association test including all genotyped variants revealed significant association with HDL-C (P = 0.024) and TG (P = 0.006). Our single-site analysis revealed seven independent signals (P < 0.05; r² < 0.40) with either HDL-C or TG. The most significant association was for the SNP rs295 exerting opposite effects on TG and HDL-C levels with P values of 7.5.10⁻⁴ and 0.002, respectively. Our work highlights some common variants and haplotypes in LPL with significant associations with lipid traits; however, the analysis of rare variants using burden tests and SKAT-O method revealed negligible effects on lipid traits. Comprehensive resequencing of LPL in larger samples is warranted to further test the role of rare variants in affecting plasma lipid levels.     

Biallelic SUN5 Mutations Cause Autosomal-Recessive Acephalic Spermatozoa Syndrome

Acephalic spermatozoa syndrome is a rare and severe form of teratozoospermia characterized by a predominance of headless spermatozoa in the ejaculate. Family clustering and consanguinity suggest a genetic origin; however, causative mutations have yet to be identified. We performed whole-exome sequencing in two unrelated infertile men and subsequent variant filtering identified one homozygous (c.824C>T [p.Thr275Met]) and one compound heterozygous (c.1006C>T [p.Arg356Cys] and c.485T>A [p.Met162Lys]) SUN5 (also named TSARG4) variants. Sanger sequencing of SUN5 in 15 additional unrelated infertile men revealed four compound heterozygous (c.381delA [p.Val128Serfs^∗7] and c.824C>T [p.Thr275Met]; c.381delA [p.Val128Serfs^∗7] and c.781G>A [p.Val261Met]; c.216G>A [p.Trp72^∗] and c.1043A>T [p.Asn348Ile]; c.425+1G>A/c.1043A>T [p.Asn348Ile]) and two homozygous (c.851C>G [p.Ser284^∗]; c.350G>A [p.Gly114Arg]) variants in six individuals. These 10 SUN5 variants were found in 8 of 17 unrelated men, explaining the genetic defect in 47.06% of the affected individuals in our cohort. These variants were absent in 100 fertile population-matched control individuals. SUN5 variants lead to absent, significantly reduced, or truncated SUN5, and certain variants altered SUN5 distribution in the head-tail junction of the sperm. In summary, these results demonstrate that biallelic SUN5 mutations cause male infertility due to autosomal-recessive acephalic spermatozoa syndrome.     

Pooled Sequencing of 531 Genes in Inflammatory Bowel Disease Identifies an Associated Rare Variant in BTNL2 and Implicates Other Immune Related Genes

The contribution of rare coding sequence variants to genetic susceptibility in complex disorders is an important but unresolved question. Most studies thus far have investigated a limited number of genes from regions which contain common disease associated variants. Here we investigate this in inflammatory bowel disease by sequencing the exons and proximal promoters of 531 genes selected from both genome-wide association studies and pathway analysis in pooled DNA panels from 474 cases of Crohn’s disease and 480 controls. 80 variants with evidence of association in the sequencing experiment or with potential functional significance were selected for follow up genotyping in 6,507 IBD cases and 3,064 population controls. The top 5 disease associated variants were genotyped in an extension panel of 3,662 IBD cases and 3,639 controls, and tested for association in a combined analysis of 10,147 IBD cases and 7,008 controls. A rare coding variant p.G454C in the BTNL2 gene within the major histocompatibility complex was significantly associated with increased risk for IBD (p = 9.65x10⁻¹⁰, OR = 2.3[95% CI = 1.75–3.04]), but was independent of the known common associated CD and UC variants at this locus. Rare (<1%) and low frequency (1–5%) variants in 3 additional genes showed suggestive association (p<0.005) with either an increased risk (ARIH2 c.338-6C>T) or decreased risk (IL12B p.V298F, and NICN p.H191R) of IBD. These results provide additional insights into the involvement of the inhibition of T cell activation in the development of both sub-phenotypes of inflammatory bowel disease. We suggest that although rare coding variants may make a modest overall contribution to complex disease susceptibility, they can inform our understanding of the molecular pathways that contribute to pathogenesis.     

Deletion 17q12 is a recurrent copy number variant that confers high risk of autism and schizophrenia

Autism spectrum disorders (ASD) and schizophrenia are neurodevelopmental disorders for which recent evidence indicates an important etiologic role for rare copy number variants (CNVs) and suggests common genetic mechanisms. We performed cytogenomic array analysis in a discovery sample of patients with neurodevelopmental disorders referred for clinical testing. We detected a recurrent 1.4 Mb deletion at 17q12, which harbors HNF1B, the gene responsible for renal cysts and diabetes syndrome (RCAD), in 18/15,749 patients, including several with ASD, but 0/4,519 controls. We identified additional shared phenotypic features among nine patients available for clinical assessment, including macrocephaly, characteristic facial features, renal anomalies, and neurocognitive impairments. In a large follow-up sample, the same deletion was identified in 2/1,182 ASD/neurocognitive impairment and in 4/6,340 schizophrenia patients, but in 0/47,929 controls (corrected p = 7.37 × 10⁻⁵). These data demonstrate that deletion 17q12 is a recurrent, pathogenic CNV that confers a very high risk for ASD and schizophrenia and show that one or more of the 15 genes in the deleted interval is dosage sensitive and essential for normal brain development and function. In addition, the phenotypic features of patients with this CNV are consistent with a contiguous gene syndrome that extends beyond RCAD, which is caused by HNF1B mutations only.     

Sequencing ASMT Identifies Rare Mutations in Chinese Han Patients with Autism

Melatonin is involved in the regulation of circadian and seasonal rhythms and immune function. Prior research reported low melatonin levels in autism spectrum disorders (ASD). ASMT located in pseudo-autosomal region 1 encodes the last enzyme of the melatonin biosynthesis pathway. A previous study reported an association between ASD and single nucleotide polymorphisms (SNPs) rs4446909 and rs5989681 located in the promoter of ASMT. Furthermore, rare deleterious mutations were identified in a subset of patients. To investigate the association between ASMT and autism, we sequenced all ASMT exons and its neighboring region in 398 Chinese Han individuals with autism and 437 healthy controls. Although our study did not detect significant differences of genotypic distribution and allele frequencies of the common SNPs in ASMT between patients with autism and healthy controls, we identified new rare coding mutations of ASMT. Among these rare variants, 4 were exclusively detected in patients with autism including a stop mutation (p.R115W, p.V166I, p.V179G, and p.W257X). These four coding variants were observed in 6 of 398 (1.51%) patients with autism and none in 437 controls (Chi-Square test, Continuity Correction p = 0.032, two-sided). Functional prediction of impact of amino acid showed that p.R115W might affect protein function. These results indicate that ASMT might be a susceptibility gene for autism. Further studies in larger samples are needed to better understand the degree of variation in this gene as well as to understand the biochemical and clinical impacts of ASMT/melatonin deficiency.     

Rare inherited missense variants of POGZ associate with autism risk and disrupt neuronal development

Excess de novo likely gene-disruptive and missense variants within dozens of genes have been identified in autism spectrum disorder(ASD)and other neurodevelopmental disorders.However,many rare inherited missense variants of these high-risk genes have not been thoroughly evaluated.In this study,we analyzed the rare missense variant burden of POGZ in a large cohort of ASD patients from the Autism Clinical and Genetic Resources in China(ACGC)and further dissected the functional effect of diseaseassociated missense variants on neuronal development.Our results showed a significant burden of rare missense variants in ASD patients compared to the control population(P=4.6×10-5,OR=3.96),and missense variants in ASD patients showed more severe predicted functional outcomes than those in controls.Furthermore,by leveraging published large-scale sequencing data of neurodevelopmental disorders(NDDs)and sporadic case reports,we identified 8 de novo missense variants of POGZ in NDD patients.Functional analysis revealed that two inherited,but not de novo,missense variants influenced the cellular localization of POGZ and failed to rescue the defects in neurite and dendritic spine development caused by Pogz knockdown in cultured mouse primary cortical neurons.Significantly,L1CAM,an autism candidate risk gene,is differentially expressed in POGZ deficient cell lines.Reduced expression of L1cam was able to partially rescue the neurite length defects caused by Pogz knockdown.Our study showed the important roles of rare inherited missense variants of POGZ in ASD risk and neuronal development and identified the potential downstream targets of POGZ,which are important for further molecular mechanism studies.