Kasra Honarmand Tamizkar1
Elham Badrlou1
Termeh Aslani2
Serge Brand3,4,5,6,7
Shahram Arsang-Jang8
Soudeh Ghafouri-Fard9*
Mohammad Taheri10*- 1Department of Medical Genetics, School of Medicine, Shahid Beheshti University of Medical Sciences, Tehran, Iran
- 2Department of Genetics, Tehran Medical Sciences Branch, Islamic Azad University, Tehran, Iran
- 3Center for Affective, Stress and Sleep Research, Psychiatric Clinics, University of Basel, Basel, Switzerland
- 4Division of Sport Science and Psychosocial Health, Department of Sport, Exercise and Health, University of Basel, Basel, Switzerland
- 5Substance Abuse Prevention Research Center, Kermanshah University of Medical Sciences, Kermanshah, Iran
- 6Sleep Disorder Research Center, Kermanshah University of Medical Sciences, Kermanshah, Iran
- 7School of Medicine, Tehran University of Medical Sciences, Tehran, Iran
- 8Cancer Gene therapy Research Center, Zanjan University of Medical Science, Zanjan, Iran
- 9Men's Health and Reproductive Health Research Center, Shahid Beheshti University of Medical Sciences, Tehran, Iran
- 10Skull Base Research Center, Loghman Hakim Hospital, Shahid Beheshti University of Medical Sciences, Tehran, Iran
Autism spectrum disorder (ASD) is a long-standing neurodevelopmental condition with prominent effects on social behavior of affected children. This disorder has been linked with neuroinflammatory responses. NF-κB has been shown to affect these responses in the orbitofrontal cortex of patients with ASD, thus being implicated in the pathogenesis of ASD. We measured expression of some NF-κB-associated lncRNAs and mRNAs (DILC, ANRIL, PACER, CHAST, ADINR, DICER1-AS1, HNF1A-AS1, NKILA, ATG5 and CEBPA) in the peripheral blood of ASD kids vs. healthy children. Expression quantities of ADINR, ANRIL, DILC, NKILA and CHAST were meaningfully higher in ASD cases compared with healthy kids (Posterior Beta = 1.402, P value < 0.0001; Posterior Beta = 2.959, P value < 0.0001; Posterior Beta = 0.882, P value = 0.012; Posterior Beta = 1.461, P value < 0.0001; Posterior Beta = 0.541, P value = 0.043, respectively). The Bonferroni corrected P values for these lncRNAs remained significant except for CHAST and DILC. Expression levels of other genes were not considerably different between cases and controls. Expressions of ATG5, DICER-AS1 and DILC were correlated with age of ASD patients (P < 0.0001). Among ASD cases, the most robust correlation has been detected between ADINR and NKILA (r = 0.87, P < 0.0001). Expression of none of genes has been correlated with age of healthy children. Among this group of children, expression levels of ADINR and CHAST were robustly correlated (r = 0.83, P < 0.0001). ANRIL had the greatest AUC value (AUC = 0.857), thus the best diagnostic power among the assessed genes. NKILA ranked the second position in this regard (AUC = 0.757). Thus, NF-κB-associated lncRNAs might partake in the pathogenesis of ASD.
Introduction
Autism spectrum disorder (ASD) is a long-standing neurodevelopmental condition pigeonholed by defects in social abilities and speaking communication and the occurrence of stereotypic behaviors and interests (American Psychiatric Association, 2013). A wide range of genetic, environmental, neurological, and immune-related parameters partake in the etiology of ASD (Neuhaus et al., 2010). Growing evidence suggests the influence of inflammatory responses in the anterior areas of the neocortex in the pathophysiology of ASD (Pardo et al., 2005; Vargas et al., 2005). Moreover, stimulation of microglia and astrocytes in the brain regions correlated with cognitive activities has led to neuroinflammatory responses in these patients (Anderson et al., 2008). Nuclear factor κ-light-chain-enhancer of activated B cells (NF-κB) has been shown to regulate response to extracellular stress and expression amounts of pro-inflammatory cytokines (Pahl, 1999; Perkins, 2004). A previous study has reported aberrant levels of NF-κB in orbitofrontal cortex of patients with ASD especially in extremely activated microglia and its role in the molecular cascade resulting in the neuroinflammation, particularly in inhabitant immune cells in brain areas linked with the behavioral changes in ASD (Young et al., 2011). Moreover, expression of IKKα kinase which phosphorylates IκBα thus enhances the inhibitory impact of this factor on NF-κB, has been remarkably elevated in the cerebellum of patients with ASD. Besides, the expression of NF-κB and its phosphorylation at Ser536 have been considerably altered in the cerebellum and cortex of ASD patients and animal model of this disorder (Malik et al., 2011). We have newly identified dysregulation of some of the NF-κB-associated genes and lncRNAs in the peripheral blood of patients with schizophrenia vs. healthy controls (Safa et al., 2020a). In the current project, we measured expression of these lncRNAs and mRNAs (DILC, ANRIL, PACER, CHAST, ADINR, DICER1-AS1, HNF1A-AS1, NKILA, ATG5 and CEBPA) in the peripheral blood of patients with ASD vs. healthy children. LncRNAs are a group of regulatory transcripts with sizes more than 200 nt to several kbs. These transcripts have influential effects on expression of genes through modulating chromatin structure, serving as enhancers for transcription, serving as molecular decoys to limit availability of other regulatory molecules and making scaf-folds for recruitment of other biomolecules (Fang and Fullwood, 2016). A bulk of evidence has shown particular involvement of these transcripts in the pathoetiology of neurodevelopmental diseases (Roberts et al., 2014).
Materials and Methods
ASD Patients and Normally Developed Children
The current research project was accomplished using the blood samples acquired from 30 ASD children (male/female ratio: 19/11) with mean age of 6 years (Standard deviation = 1.39). Forty-one age- and sex-matched healthy children (male/female ratio: 30/11) from the same ethnic group were selected as controls. ASD status was diagnosed by a psychiatric according to the DSM-V (American Psychiatric Association, 2013). Children having other neuropsychiatric, metabolic or immune-related conditions were omitted from the project. None of the ASD patients had a comorbid disorder such as attention deficit hyperactive disorder. Children registered as controls had no personal or familial history of neuropsychiatric disorders or developmental problem. Written informed consent forms were signed by all parents. The study protocol was approved by the ethical committee of Shahid Beheshti University of Medical Science and all methods were performed in accordance with the relevant guidelines and regulations (IR.SBMU.RETECH.REC.1399.740).
Evaluation of Transcript Levels of Genes
RNA was retrieved from blood specimens using the extraction kit (Hybrid-RTM Blood RNA) purchased from GeneAll Company (Seoul, South Korea). Next, 10 μg of extracted RNA specimens was converted to cDNA using the BioFact RT kit (Seoul, South Korea). Expression of genes was quantified in all specimens using the RealQ Plus 2x Master Mix (Amplicon, Denmark) and the primers which are shown in Table 1. All experiments were performed in duplicate. Each PCR run included a negative control (no template control) for each set of primers. Based on our previous studies validating the constant levels of B2M gene in the blood (Fallah et al., 2019), expression quantities of B2M gene were measured to normalize the expression data of target genes.
Table 1. Information about primers and the corresponding amplified region.
Statistical Strategies
Statistical methods were executed in Stan, “ggplot2,” “brms” and pROC packages in the R v.4 environment. The Spearman correlation coefficients were calculated to judge the correlation between expressions of the selected genes. Expression of these genes were compared between ASD children and normal children using the Bayesian regression model, the asymmetric Laplace family prior with 2,000 burnouts and 5,000 iterations. Rhat, posterior predictive plots and Loo were used. P values were corrected for multiple comparisons. To obtain the Bonferroni corrected P value, the original α-values were divided by the quantities of analyses on the dependent variable. Receiver operating characteristic (ROC) curves were illustrated to evaluate the diagnostic power of genes.
Results
Expression quantities of ADINR, ANRIL, DILC, NKILA and CHAST were meaningfully higher in ASD cases compared with healthy kids (Posterior Beta = 1.402, P value < 0.0001; Posterior Beta = 2.959, P value < 0.0001; Posterior Beta = 0.882, P value = 0.012; Posterior Beta = 1.461, P value < 0.0001; Posterior Beta = 0.541, P value = 0.043, respectively). The Bonferroni corrected P values for these lncRNAs remained significant except for CHAST and DILC.
Expression levels of other genes were not considerably different between cases and controls. Table 2 reports the detailed parameters obtained from Bayesian regression model for evaluation of gene expression between ASD children and control subjects.
Table 2. Detailed parameters obtained from Bayesian regression model for comparison of gene expression between ASD children and control subjects (Expression Ratio: case/control; Gender: girl/boy; P-Value was calculated from Median regression model using bootstrap).
Figure 1 demonstrates the levels of assessed genes among ASD children and normal children.
Figure 1. Relative levels of assessed genes among ASD children and healthy children.
Correlation Between Expressions of Genes
Expressions of ATG5, DICER-AS1 and DILC were correlated with age of ASD patients (P < 0.0001). Among ASD cases, the most considerable correlation has been detected between ADINR and NKILA (r = 0.87, P < 0.0001) (Figure 2).
Figure 2. Correlation between expressions of genes among patients with ASD.
Expression of none of genes has been correlated with age of healthy children. Among this group of children, expression levels of ADINR and CHAST were robustly correlated (r = 0.83, P < 0.0001) (Figure 3).
Figure 3. Correlation between expressions of genes among healthy controls.
ROC Curves
ANRIL had the greatest AUC value (AUC = 0.857), thus the best diagnostic power among the assessed genes. NKILA ranked the second position in this regard (AUC = 0.757). Figure 4 displays the plotted ROC curves for distinguishing ASD from healthy status.
Figure 4. Receiver Operating Characteristics (ROC) curves for appraisal of diagnostic value of ADINR, ANRIL, DILC, NKILA and CHAST in ASD.
Discussion
NF-κB has an essential role in the regulation of immune responses which partake in the pathogenesis of ASD. Based on our recent experience in patients with schizophrenia regarding the possible role of NF-κB-related genes, we hypothesized that these genes might also be involved in ASD. This hypothesis has been supported by the existence of inflammatory responses at the cross point of shared underlying mechanisms for schizophrenia and autism (Meyer et al., 2011). Notably, we detected over-expression of ADINR, ANRIL, DILC, NKILA and CHAST in the peripheral blood of ASD cases compared with controls. However, the Bonferroni corrected P values for CHAST and DILC did not remain significant. ADINR is an lncRNA which is transcribed from a region upstream of the CEBPA gene, regulating its expression. This lncRNA binds with PA1 and recruits histone methyl-transferase complexes to alter histone epigenetic marks in the CEBPA locus (Xiao et al., 2015). Although the role of CEBPA in the pathogenesis of ASD has not been elucidated yet, CEBPA has functional interactions with the ASD-related protein CHD8 (Kita et al., 2018). Moreover, CEBPA has been shown to modulate nfkb1 expression and relocate histone deacetylases from NF-κB p50 homodimers to activate expression of NF-κB target genes (Paz-Priel et al., 2011). Therefore, ADINR might participate in the pathogenesis of ASD through modulating CEBPA or CHD8 activity or expression levels. We have recently demonstrated a trend toward association between some ANRIL haplotypes and risk of ASD among Iranians (Safa et al., 2020b). Moreover, ANRIL has an acknowledged role in the regulation of immune responses as an element of NF-κB pathway (Zhou et al., 2016). DILC has an appreciated role in the pathogenesis a the immune-mediated disorder rheumatoid arthritis through modulating IL-6 expression (Wang et al., 2019). IL-6 is possibly involved in the pathogenesis of ASD since its expression has been elevated IL-6 is increased in the cerebellum of patients with ASD and it changes adhesion and migration of neurons and modulates development of synapses (Wei et al., 2011). NKILA has been shown to decline TNF-α-associated inflammatory responses (Han et al., 2020). Notably, cytokine profiling has recognized TNF-α as an important dysregulated cytokine in patients with ASD (Xie et al., 2017). Finally, CHAST over-expression has also been detected in patients with schizophrenia (Safa et al., 2020a). Expression levels of other genes were not considerably different between cases and controls.
It is worth mentioning that levels of lncRNAs circulating in peripheral blood do not necessarily correlate with their levels in the brain. Since brain tissues are not available from these patients, we could not evaluate correlations between these two sets of samples.
Notably, expressions of ATG5, DICER-AS1 and DILC were correlated with age of ASD patients, but not healthy subjects, indicating different impact of age-related factors among patients and healthy controls. Previous studies have reported age-dependent changes in genes expressions in the brain tissues of ASD cases. For instance, Chow et al. have shown altered expressions of pathways controlling cell number, cortical patterning, and differentiation in prefrontal cortex samples of young ASD cases. On the other hand, signaling and repair pathways have been found to be dysregulated in prefrontal cortex samples of adult ASD cases. They have concluded that age-dependent alterations in gene signature in ASD might represent distinctive abnormal processes in different phases of neurodevelopment (Chow et al., 2012). Moreover, experiments in lymphoblastoid cell lines have shown an obvious inconsistency between neuroanatomical and cellular aberrations detected in ASD cases at younger ages and molecular abnormalities at higher age (Ansel et al., 2017). Thus, identification of age-related changes in gene expression might facilitate recognition of the most important biomedical pathways in pathoetiology of ASD in each age.
The most robust correlations have been detected between ADINR and NKILA among cases, and between ADINR and CHAST among healthy persons, indicating the impact of ASD on modulation of the correlation between NF-κB-associated genes.
Finally, we appraised the suitability of genes in distinguishing between ASD and healthy status. ANRIL had the greatest AUC value, thus the best diagnostic power among the assessed genes. NKILA ranked the second position in this regard. Thus, NF-κB-associated lncRNAs might partake in the pathogenesis of ASD and can be used as diagnostic markers in ASD. However, the AUC data need to be replicated in a larger sample in which ASD diagnosis is established using a standardized diagnostic tool such as the ADOS or ADI-R before it can be clinically useful. Moreover, functional experiments are desired to appraise the molecular cascade of involvement of NF-κB-associated lncRNAs in ASD.
Data Availability Statement
The original contributions presented in the study are included in the article/supplementary files, further inquiries can be directed to the corresponding authors.
Ethics Statement
The studies involving human participants were reviewed and approved by the Ethical Committee of Shahid Beheshti University of Medical Science and all methods were performed in accordance with the relevant guidelines and regulations (IR.SBMU.RETECH.REC.1399.740). Written informed consent to participate in this study was provided by the participants' legal guardian/next of kin.
Author Contributions
EB, TA, and MT performed the experiment. SA-J analyzed the data. SG-F and KH designed and supervised the study. SG-F wrote and revised the draft. All the authors are contributed equally and approved the submitted version.
Funding
The current study was supported by a grant from Shahid Beheshti University of Medical Sciences.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
American Psychiatric Association (2013). Diagnostic and statistical Manual of Mental Disorders (DSM-5®). Washington, DC: American Psychiatric Publishing. doi: 10.1176/appi.journals.9780890425596
Anderson, M. P., Hooker, B. S., and Herbert, M. R., (2008). Bridging from cells to cognition in autism pathophysiology: biological pathways to defective brain function and plasticity. Am. J. Biochem. Biotechnol. 4, 167–176. doi: 10.3844/ajbbsp.2008.167.176
Ansel, A., Rosenzweig, J. P., Zisman, P. D., Melamed, M., and Gesundheit, B. (2017). Variation in gene expression in autism spectrum disorders: an extensive review of transcriptomic studies. Front. Neurosci. 10:601. doi: 10.3389/fnins.2016.00601
Chow, M. L., Pramparo, T., Winn, M. E., Barnes, C. C., Li, H. R., Weiss, L., et al. (2012). Age-dependent brain gene expression and copy number anomalies in autism suggest distinct pathological processes at young versus mature ages. PLoS Genet. 8:e1002592. doi: 10.1371/journal.pgen.1002592
Fallah, H., Azari, I., Neishabouri, S. M., Oskooei, V. K., Taheri, M., and Ghafouri-Fard, S. (2019). Sex-specific up-regulation of lncRNAs in peripheral blood of patients with schizophrenia. Sci. Rep. 9:12737. doi: 10.1038/s41598-019-49265-z
Fang, Y., and Fullwood, M. J. (2016). Roles, functions, and mechanisms of long non-coding RNAs in cancer. Genomics Proteomics Bioinformatics 14, 42–54. doi: 10.1016/j.gpb.2015.09.006
Han, Y., Cai, Y., Lai, X., Wang, Z., Wei, S., Tan, K., et al. (2020). lncRNA RMRP prevents mitochondrial dysfunction and cardiomyocyte apoptosis via the miR-1-5p/hsp70 axis in LPS-induced sepsis mice. Inflammation 43:605–18. doi: 10.1007/s10753-019-01141-8
Kita, Y., Katayama, Y., Shiraishi, T., Oka, T., Sato, T., Suyama, M., et al. (2018). The autism-related protein CHD8 cooperates with C/EBPβ to regulate adipogenesis. Cell Rep. 23, 1988–2000. doi: 10.1016/j.celrep.2018.04.050
Malik, M., Tauqeer, Z., Sheikh, A. M., Wen, G., Nagori, A., Yang, K., et al. (2011). NF-κB signaling in the brain of autistic subjects. Mediators Inflamm. 2011:785265. doi: 10.1155/2011/785265
Meyer, U., Feldon, J., and Dammann, O. (2011). Schizophrenia and autism: both shared and disorder-specific pathogenesis via perinatal inflammation? Pediatr. Res. 69(5 Pt 2), 26R−33R. doi: 10.1203/PDR.0b013e318212c196
Neuhaus, E., Beauchaine, T. P., and Bernier, R. (2010). Neurobiological correlates of social functioning in autism. Clin. Psychol. Rev. 30, 733–748. doi: 10.1016/j.cpr.2010.05.007
Pahl, H. L. (1999). Activators and target genes of Rel/NF-kappaB transcription factors. Oncogene 18, 6853–6866. doi: 10.1038/sj.onc.1203239
Pardo, C. A., Vargas, D. L., and Zimmerman, A. W. (2005). Immunity, neuroglia and neuroinflammation in autism. Int. Rev. Psychiatry (Abingdon, England) 17, 485–495. doi: 10.1080/02646830500381930
Paz-Priel, I., Houng, S., Dooher, J., and Friedman, A. D. (2011). C/EBPα and C/EBPα oncoproteins regulate nfkb1 and displace histone deacetylases from NF-κB p50 homodimers to induce NF-κB target genes. Blood 117, 4085–4094. doi: 10.1182/blood-2010-07-294470
Perkins, N., (2004). NF-kB: tumour promoter or tumour suppressor? Br. J. Cancer. 91, pp. S7–S7. doi: 10.1016/j.tcb.2003.12.004
Roberts, T. C., Morris, K. V., and Wood, M. J. A. (2014). The role of long non-coding RNAs in neurodevelopment, brain function and neurological disease. Philos. Trans. R. Soc. Lond. B. Biol. Sci. 369:20130507. doi: 10.1098/rstb.2013.0507
Safa, A., Badrlou, E., Arsang-Jang, S., Sayad, A., Taheri, M., and Ghafouri-Fard, S. (2020a). Expression of NF-κB associated lncRNAs in schizophrenia. Sci. Rep. 10:18105. doi: 10.1038/s41598-020-75333-w
Safa, A., Noroozi, R., Taheri, M., and Ghafouri-Fard, S. (2020b). Association analysis of ANRIL polymorphisms and haplotypes with autism spectrum disorders. J. Mol. Neurosci. 71, 1–6. doi: 10.1007/s12031-020-01657-x
Vargas, D. L., Nascimbene, C., Krishnan, C., Zimmerman, A. W., and Pardo, C. A. (2005). Neuroglial activation and neuroinflammation in the brain of patients with autism. Ann. Neurol. 57, 67–81. doi: 10.1002/ana.20315
Wang, G., Tang, L., Zhang, X., and Li, Y. (2019). LncRNA DILC participates in rheumatoid arthritis by inducing apoptosis of fibroblast-like synoviocytes and down-regulating IL-6. Biosci. Rep. 39:BSR20182374. doi: 10.1042/BSR20182374
Wei, H., Zou, H., Sheikh, A. M., Malik, M., Dobkin, C., Brown, W. T., et al. (2011). IL-6 is increased in the cerebellum of autistic brain and alters neural cell adhesion, migration and synaptic formation. J. Neuroinflamm. 8:52. doi: 10.1186/1742-2094-8-52
Xiao, T., Liu, L., Li, H., Sun, Y., Luo, H., Li, T., et al. (2015). Long noncoding RNA ADINR regulates adipogenesis by transcriptionally activating C/EBPα. Stem Cell Rep. 5, 856–865. doi: 10.1016/j.stemcr.2015.09.007
Xie, J., Huang, L., Li, X., Li, H., Zhou, Y., Zhu, H., et al. (2017). Immunological cytokine profiling identifies TNF-α as a key molecule dysregulated in autistic children. Oncotarget 8, 82390–82398. doi: 10.18632/oncotarget.19326
Young, A. M. H., Campbell, E., Lynch, S., Suckling, J., and Powis, S. J. (2011). Aberrant NF-kappaB expression in autism spectrum condition: a mechanism for neuroinflammation. Front. Psychiatry 2:27. doi: 10.3389/fpsyt.2011.00027
Keywords: autism spectrum disorder, ADINR, ANRIL, DILC, NKILA, CHAST
Citation: Honarmand Tamizkar K, Badrlou E, Aslani T, Brand S, Arsang-Jang S, Ghafouri-Fard S and Taheri M (2021) Dysregulation of NF-κB-Associated LncRNAs in Autism Spectrum Disorder. Front. Mol. Neurosci. 14:747785. doi: 10.3389/fnmol.2021.747785
Received: 26 July 2021; Accepted: 01 September 2021;
Published: 28 September 2021.
Edited by:
Mukesh K. Pandey, Mayo Clinic, United StatesReviewed by:
Supriya Behl, Mayo Clinic, United StatesAdeline Lau, Flinders University, Australia
Sunil Mehta, Mayo Clinic, United States
Copyright © 2021 Honarmand Tamizkar, Badrlou, Aslani, Brand, Arsang-Jang, Ghafouri-Fard and Taheri. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Soudeh Ghafouri-Fard, s.ghafourifard@sbmu.ac.ir; Mohammad Taheri, mohammad_823@yahoo.com