- Primary research
- Open Access
Analysis of the SNP rs3747333 and rs3747334 in NLGN4X gene in autism spectrum disorder: a meta-analysis
Annals of General Psychiatry volume 18, Article number: 6 (2019)
The SNP rs3747333 and rs3747334 in Neuroligin 4X (NLGN4X) gene have been demonstrated to be associated with the susceptibility to Autism spectrum disorder (ASDs; MIM 209850), but the results are inconsistent. Therefore, a meta-analysis of eligible studies reporting the association between rs3747333 and rs3747334 and ASD was carried out to enhance the reliability of published results.
A systematic literature search was performed using PubMed, Web of Science, Cochrane Library to search English articles concerning the relation between rs3747333, rs3747334 and ASD up to Sep. 21th, 2017. Summary odds ratios (OR) and 95% confidence interval (CI) were used to evaluate the risk of rs3747333, rs3747334 in the ASD. The heterogeneity and publication bias of the eligible studies were also evaluated.
Six eligible studies involving 1284 subjects (735 patients and 549 healthy controls) were included in this meta-analysis. Overall, the results indicated that there was no significant risk elevation between rs3747333, rs3747334 variants and ASD (OR = 0.39, 95% CI 0.10–1.60). Furthermore, sensitivity analysis and publication bias analysis confirmed this result.
In conclusion, our meta-analysis suggests that the rs3747333, rs3747334 in NLGN4X gene are not frequent causes of ASD.
Introduction of autism spectrum disorder
Autism spectrum disorder (ASDs; MIM 209850) is a heterogeneous group of neurodevelopmental diseases defined by significant language delays, social, and communication deficits, and highly repetitive and stereotypic interest and/or behaviors . The global prevalence of autism has risen to 6.2% . Autism patients in China are estimated in the range of 0.6–1.8 million according to the WHO’s report. According to the latest report released by Autistica, people with ASD die on average 18 years before the general population . ASD is a result of complex gene–environment interactions, with strong and clear genetic influences . Although the etiologic and clinical heterogeneity are universally recognized, in practice, many studies still fail to take this into account. As of December, 2015, genetic studies have pointed to 791 potential candidate genes in ASD (https://gene.sfari.org/autdb/HG_Home.do), making it difficult to find common underlying pathogenic mechanisms.
The risk of the rs3747333 and rs3747334 in the autism spectrum disorder
Neuroligin 4X (NLGN4X) is a member of neuroligins expressed in the postsynaptic neurons and mediate transsynaptic signaling by interacting their ligand, neurexins . Several studies identified significant variants in the NLGN4X gene that were replicated in two or more independent samples and were associated with specific phenotypes of ASD [6,7,8]. However, other studies found no variants in the NLGN4X gene in autistic patients in their cohorts [9,10,11]. To date, many common single-nucleotide polymorphisms (SNPs) of NLGN4X such as rs72413786, rs2290487, rs7049300, rs3747333, rs3747334, have been studied in extensive research on ASD. Of these, the SNP rs3747333 and rs3747334 have received more favorable observations concerning on the ASD risk [11,12,13,14,15,16]. A case–control study of Chinese descendants found that the common SNPs (rs3747333 and rs3747334) in NLGN4X gene significantly associated with risk for autism , a previous work showed that seven boys were carriers of mutations at three SNPs including the synonymous SNP rs7049300, the non-synonymous SNP rs3747333 (L593) and the synonymous SNP rs3747334 (L593), but the co-occurrence of rs3747333C>T leading to a substitution of leucine to phenylalanine and rs3747334C>G entails a double substitution of CTC to TTG in codon 593 and thus protect the amino acid leucine . However, a work from German showed that the polymorphisms rs 3747333 and rs 3747334 did not cause any amino acid substitutions in the total of the eight detected carriers . Those inconsistencies could be attributed to the phenotypic heterogeneity of ASD and to relatively small sample sizes.
Given the increasing attention directed to the inconsistencies results, to further clarify the inconclusive association between the two SNPs and ASD, we conducted a meta-analysis with available data published up to Sep. 21th, 2017.
Materials and methods
Literature search and selection strategy
A systematic literature search was conducted to collect studies relevant for the current meta-analysis. First, we searched three online databases (PubMed, Web of Science, Cochrane Library) for articles published in English on associations between the rs3747333, rs3747334 in NLGN4X and ASD up to Sep. 21th, 2017. Combinations of terms identifying the gene (rs3747333, rs3747334, NLGN4X, X-linked Neuroligin 4 gene) and terms identifying the phenotype of interest (autism, autism spectrum disorder, Asperger, childhood disintegrative disorder, pervasive developmental disorder) were entered as Medical Subject Heading terms and text words. After identifying a set of publications returned from each search, we collated the results and eliminated duplicates. We then reviewed the titles, abstracts, keywords and texts of all studies to exclude studies that were clearly irrelevant. Hand-searching was carried out to determine other potential eligible studies through scanning the references cited in the retrieved articles. Finally, we applied criteria to the set of studies to determine their inclusion or exclusion in the present meta-analysis. If the two reviewers (Hongli Sun and Ying Yang) disagreed, all the authors critically evaluated the studies to judge of the inclusion or exclusion of a certain study.
All the eligible articles meet the inclusion criteria:
The study was peer-reviewed to ensure high levels of methodological adequacy and to avoid the inevitable bias caused by dependence on investigators agreeing to provide data from unpublished studies, as suggested by the Cochrane Group .
The study was adopted as “Original Research”.
Assessment of the association between the rs3747333 or rs3747334 or both of them in NLGN4X and ASD.
Written in English.
The two authors independently extracted the following data: (1) first author, (2) year of publication, (3) location, (4) predominant ethnicity of participants, (5) study design, (6) number of participants, (7) diagnostic criteria. When the study did not report sufficient data to calculate effect size, we contacted corresponding authors via e-mail to request additional information. If the contact with the author was not successful, we could not include the study in the meta-analysis.
The association of the rs3747333, rs3747334 in NLGN4X and ASD was assessed by odds ratio (OR) with 95% confidence interval (CI). Subgroup (according to study design) analysis was evaluated by study type. Publication bias was considered present when P-value was less than 0.05. Sensitivity analysis was also carried out to evaluate the stability of the meta-analysis. Statistical heterogeneity among studies was evaluated by the Chi-square-based Q-test and the I2 index. I2 index < 50% and/or P-value > 0.10 for Q-test indicated a lack of heterogeneity among the studies. The fixed-effect model (the Mantel–Haenszel method) was utilized to calculate the pooled ORs when the heterogeneity was not significant among the studies. Otherwise, the random-effects model (the DerSimonian and Laird method) was employed. Publication bias was evaluated by the funnel plot. All statistical analyses were carried out with the STATA 12.0 (STATA Corp, College Station, Texas).
Meta-analyses of the rs3747333, rs3747334 in NLGN4X and ASD
After applying inclusion and exclusion criteria (see “Materials and methods”), data from 6 studies with 1284 subjects [9, 11, 14,15,16, 18] were included in the meta-analysis (Fig. 1). We assessed the association between the rs3747333, rs3747334, and ASD. The characteristics of the populations used for meta-analysis are shown in Table 1.
Regarding the rs3747333, rs3747334, no statistically significant correlation with ASD (OR = 0.39, 95% CI 0.10–1.60) was observed under the random-effects model. There was substantial heterogeneity in the estimated effect sizes (I-squared: I2 = 95.1%, P = 0.000). We further performed the subgroup analysis by study design and the result showed that there was no statistical significant association existed between the rs3747333, rs3747334 and the risk of ASD under random-effects model (Screening: OR = 0.13, 95% CI 0.05–0.31; Case–control: OR = 3.68, 95% CI 0.40–33.73). There was substantial heterogeneity in the estimated effect sizes (Screening: I-squared: I2 = 80.5%, P = 0.002; case–control: I-squared: I2 = 92.3%, P = 0.000) (Fig. 2).
Sensitivity analysis and publication bias
The aim of sensitivity analysis was to assess the effects of each study on the pooled OR and ensure that none of solo study was completely responsible for the combined results. The results of sensitivity analysis showed that the pooled OR was not materially affected before and after removal of individual study (see Additional file 1: Figures S1–S6). The studies distribution of the funnel plot was substantially symmetrical about the combined effect size in Fig. 3. Visual inspection of the funnel plot revealed obvious asymmetry.
In the current meta-analysis, pooling the data from 6 published studies indicated the rs3747333 and rs3747334 in NLGN4X were not a risk for ASD. The rs3747333, rs3747334 are linked in exon 6 of NLGN4X, which localized on Xq13.1 and Xp22.33, encode a member of a family of neuronal cell surface proteins serving as splice site specific ligands for neurexins and may be involved in the formation, organization, and remodeling of central nervous system synapses. In our meta-analysis, the most works showed that no evidence for involvement of the rs3747333, rs3747334 in NLGN4X in probands with ASD [9, 11, 14, 15, 18], but Xiaojuan’ work form Chinese sample, conversely, identified that the rs3747333 and rs3747334 in NLGN4X gene significantly associated with the risk for ASD . The inconsistent results of these works are, to a very great extent, determined by ethnicity. However, it is needed to recollect more Chinese sample to identify the risk of the rs3747333, rs3747334 in NLGN4X for ASD.
The main limitation of this meta-analysis is that there was a significant heterogeneity in the estimated effect sizes across studies for the rs3747333, rs3747334 with ASD, suggesting that potential moderator variables, such as ethnicity, may cause the observed differences in effect sizes across studies. Another limitation is that we collected too few studies to reliably test moderators of that heterogeneity, and there may be some potential bias. We performed sensitivity analysis, but the results showed that the pooled OR was not altered materially when ruling out each individual study. A recent meta-analysis using genome-wide association (GWAS) analysis showed that over 16,000 individuals with highlights a novel locus at 10q24.32 and a significant overlap with schizophrenia , which made a great contribution to understand the complicated pathogenesis of autism spectrum disorder. Thus, genome-wide association analysis as a better approach to the study of the genetics of ASD should be applied in the future work.
In summary, our meta-analysis indicates that there was no significant risk elevation between the rs3747333, rs3747334 variants and ASD based on the currently published data. However, the meta-analysis adds new information to the ongoing debate relating to the risk of NLGN4X for ASD. Further investigations of these findings are merited by well-designed and large-sample multi-centric studies.
Availability of data and materials
DSM-5 American Psychiatric Association. Diagnostic and statistical manual of mental disorders. Arlington: American Psychiatric Publishing; 2013.
Gauthier J, Bonnel A, Stonge J, Karemera L, Laurent S, Mottron L, Fombonne E, Joober R, Rouleau GA. NLGN3/NLGN4 gene mutations are not responsible for autism in the Quebec population. Am J Med Genet Part B Neuropsychiatr Genet. 2005;132B:74.
Xu X, Xiong Z, Zhang L, Liu Y, Lu L, Peng Y, Guo H, Zhao J, Xia K, Hu Z. Variations analysis of NLGN3 and NLGN4X gene in Chinese autism patients. Mol Biol Rep. 2014;41:4133–40.
Südhof TC. Neuroligins and neurexins link synaptic function to cognitive disease. Nature. 2008;455:903–11.
Volaki K, Pampanos A, Kitsioutzeli S, Vrettou C, Oikonomakis V, Sofocleous C, Kanavakis E. Mutation screening in the Greek population and evaluation of NLGN3 and NLGN4X genes causal factors for autism. Psychiatr Genet. 2013;23:198.
Elsabbagh M, Divan G, Koh YJ, Kim YS, Kauchali S, Marcín C, Montiel-Nava C, Patel V, Paula CS, Wang C. Global prevalence of autism and other pervasive developmental disorders. Autism Res. 2012;5:160–79.
Liu Y, Du Y, Liu W, Yang C, Liu Y, Wang H, Gong X. Lack of association between NLGN3, NLGN4, SHANK2 and SHANK3 gene variants and autism spectrum disorder in a Chinese population. PLoS ONE. 2013;8:e56639.
Wermter AK, Kamp-Becker I, Strauch K, Schulte-Körne G, Remschmidt H. No evidence for involvement of genetic variants in the X-linked neuroligin genes NLGN3 and NLGN4X in probands with autism spectrum disorder on high functioning level. Am J Med Genet B Neuropsychiatr Genet. 2008;147B:535–7.
Underwood E. People on autism spectrum die 18 years younger than average. Science. 2016.
Francesca B, Bacchelli E, Pesaresi G, Carone S, Bailey AJ, Maestrini E. Absence of coding mutations in the X-linked genes neuroligin 3 and neuroligin 4 in individuals with autism from the IMGSAC collection. Am J Med Genet Part B Neuropsychiatr Genet. 2006;141B:220.
Gauthier J, Bonnel A, St-Onge J, Karemera L, Laurent S, Mottron L, Fombonne E, Joober R, Rouleau GA. NLGN3/NLGN4 gene mutations are not responsible for autism in the Quebec population. Am J Med Genet B Neuropsychiatr Genet. 2005;132B:74–5.
Chaste P, Leboyer M. Autism risk factors: genes, environment, and gene-environment interactions. Dialogues Clin Neurosci. 2012;14:281.
Avdjieva-Tzavella DM, Todorov TP, Todorova AP, Kirov AV, Hadjidekova SP, Rukova BB, Litvinenko IO, Hristova-Naydenova DN, Tincheva RS, Toncheva DI. Analysis of the genes encoding neuroligins NLGN3 and NLGN4 in Bulgarian patients with autism. Genet Couns. 2012;23:505–11.
Mikhailov A, Fennell A, Plong-On O, Sripo T, Hansakunachai T, Roongpraiwan R, Sombuntham T, Ruangdaraganon N, Vincent JB, Limprasert P. Screening of NLGN3 and NLGN4X genes in Thai children with autism spectrum disorder. Psychiatr Genet. 2014;24:42–3.
Higgins JPT, Green S. Cochrane handbook for systematic reviews of interventions, vol. 5. Wiley Online Library; 2008.
Holmans P. Meta-analysis of GWAS of over 16,000 individuals with autism spectrum disorder highlights a novel locus at 10q24.32 and a significant overlap with schizophrenia. Mol Autism. 2017;8:21.
Lawson-Yuen A, Saldivar JS, Sommer S, Picker J. Familial deletion within NLGN4 associated with autism and Tourette syndrome. Eur J Hum Genet. 2008;16:614–8.
Laumonnier F, Bonnet-Brilhault F, Gomot M, Blanc R, David A, Moizard MP, Raynaud M, Ronce N, Lemonnier E, Calvas P, Laudier B, Chelly J, Fryns JP, Ropers HH, Hamel BC, Andres C, Barthelemy C, Moraine C, Briault S. X-linked mental retardation and autism are associated with a mutation in the NLGN4 gene, a member of the neuroligin family. Am J Hum Genet. 2004;74:552–7.
Blasi F, Bacchelli E, Pesaresi G, Carone S, Bailey AJ, Maestrini E. Absence of coding mutations in the X-linked genes neuroligin 3 and neuroligin 4 in individuals with autism from the IMGSAC collection. Am J Med Genet B Neuropsychiatr Genet. 2006;141B:220–1.
We thank Wei Fan for the meta-analysis.
This study was supported by Natural Science Foundation of China (No. 81701331) and Xi’an science and technology project (No.SF1510(4)).
Ethics approval and consent to participate
Consent for publication
The authors provide consent for publication.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.