Immune function genes CD99L2, JARID2 and TPO show association with autism spectrum disorder
© Ramos et al.; licensee BioMed Central Ltd 2012
Received: 16 December 2011
Accepted: 9 June 2012
Published: 9 June 2012
A growing number of clinical and basic research studies have implicated immunological abnormalities as being associated with and potentially responsible for the cognitive and behavioral deficits seen in autism spectrum disorder (ASD) children. Here we test the hypothesis that immune-related gene loci are associated with ASD.
We identified 2,012 genes of known immune-function via Ingenuity Pathway Analysis. Family-based tests of association were computed on the 22,904 single nucleotide polymorphisms (SNPs) from the 2,012 immune-related genes on 1,510 trios available at the Autism Genetic Resource Exchange (AGRE) repository. Several SNPs in immune-related genes remained statistically significantly associated with ASD after adjusting for multiple comparisons. Specifically, we observed significant associations in the CD99 molecule-like 2 region (CD99L2, rs11796490, P = 4.01 × 10-06, OR = 0.68 (0.58-0.80)), in the jumonji AT rich interactive domain 2 (JARID2) gene (rs13193457, P = 2.71 × 10-06, OR = 0.61 (0.49-0.75)), and in the thyroid peroxidase gene (TPO) (rs1514687, P = 5.72 × 10-06, OR = 1.46 (1.24-1.72)).
This study suggests that despite the lack of a general enrichment of SNPs in immune function genes in ASD children, several novel genes with known immune functions are associated with ASD.
KeywordsAutism Spectrum Disorders immune loci family-based association analysis single-nucleotide polymorphism
Autism spectrum disorders (ASDs) encompass a heterogeneous group of clinical descriptors whose core features include deficits in cognition, communication and social acuity, coupled with stereotypical behaviors. Currently, 1 in 88 children in the United States have an ASD diagnosis, with boys being affected at a ratio of 5:1 compared to girls . The etiology of ASDs, with very few exceptions, is unknown although a growing number of clinical and basic research studies have implicated immunological abnormalities as being associated with and potentially responsible for the cognitive and behavioral deficits seen in ASD children [2–4].
A judicious and efficient way to identify genetic variation predisposing to complex diseases is the application of a hypothesis-driven framework that incorporates prior biological knowledge. Unlike agnostic genome wide association studies (GWAS), this approach narrows the hypothesis space to provide a more focused and powerful examination of the data. Indeed, several studies have shown that a hypothesis-driven approach of selecting candidate genes serves to increase the reliability and likelihood of finding genes that are truly associated with disease [5, 6]. In accord with these considerations, we tested for associations between immune-related genes and ASD children.
We used Ingenuity Pathway Analysis (IPA) , a bioinformatics database, to search for molecules with known immune function. IPA offers the largest and most comprehensive set of functional annotations, integrating manually curated data from other databases, and a broad range of scientific literature. Briefly, we searched for “immune” under Functions and Diseases, exported the molecule annotations, and used MatchMiner  to get the positions of the corresponding 2,012 immune-function genes (Additional file 1: Table S1).
We then used the genotypic data available at the Autism Genetic Resource Exchange (AGRE)  repository to perform a family-based association analysis between variants in those 2,012 immune loci and ASD. These data have previously been used in several GWAS [10–12], and the families used in our analysis are, with the exception of minor differences due to slightly different quality control criteria, the same ones as reported in Ma et al., Wang et al. and Weiss et al.. These data and families also overlap with those reported by Voineagu et al., who, included in their expression profiling study, used the AGRE data to perform an analysis searching for evidence of a genetic enrichment of immune loci. We searched for SNPs in a region including 5 kb up- and downstream of each immune gene. After applying quality control filters, we identified 22,904 SNPs genotyped in the AGRE collection of 1,510 trios on the Illumina Hap550 platform. We tested these SNPs for association with ASD using the standard Transmission Disequilibrium Test (TDT), as implemented in PLINK v1.07 . In addition, given that the 1,510 trios are part of 1,057 independent nuclear families, we also used the Pedigree Disequilibrium Test (PDT), which uses data from related nuclear families and discordant sibships from extended pedigrees . The following exclusion filters were applied: minor allele frequency (MAF) < 0.01, SNP genotyping missing rate > 10%, individuals > 10% missing genotypes, Hardy-Weinberg Equilibrium (HWE) P-value < 0.001 (founders), Mendelian error rate > 5% per family and > 4% per SNP.
Since we hypothesized that ASD children may be enriched for risk loci in immune genes, we tested for a potential enrichment of significant P-values for SNPs in immune-related genes. We tested for an enrichment of association signals in these genes using INRICH , which corrects the empirical gene-set P-value using Bootstrapping-based re-sampling. In addition, we assessed the distribution of the linkage-disequilibrium (LD)-pruned immune-gene-set with random LD-pruned SNPs, using an r2 > 0.2 for LD-based SNP pruning. We did not observe an enrichment of immune-gene associations with either test (P = 1.0 and P = 0.45, respectively).
Most significant immune loci from the association analyses
OR (95% CI)
Finally, in order to fine map the novel loci, we used BEAGLE v3.3.2  and the 1000 Genomes Project Phase 1 (version 3; http://www.1000genomes.org/) reference panel to impute all SNPs in the regions of the top genes.
Novel immune genes
The lack of an enrichment of significant associations in immune genes corroborates the findings by Voineagu et al., who used a complimentary approach to test for evidence of a genetic component for the up-regulation of immune response genes in the autistic brain using the same AGRE genotype data. Nevertheless, we uncovered a few associations in immune-related genes that met an FDR-adjustment.
A variant in the first intron of the jumonji AT rich interactive domain 2 (JARID2) gene showed association with ASD (rs13193457, P = 2.71 × 10-06, OR = 0.61 (0.49-0.75) (Table 1). Despite having a borderline Mendelian error rate, this SNP is of interest. This gene has been previously implicated in schizophrenia [19, 20] and, more recently, it has been reported in a GWAS for ASD that also included the same AGRE dataset herein used . The variant reported by Weiss et al. (rs7766973) does not show LD with rs13193457 (r2 = 0.13 in CEU from 1000 Genomes Pilot 1). The jumonji-domain functions by removing methyl marks on histones that are associated with gene regulation . JARID2 is required for neural tube formation and is essential for normal heart development and function, as well as acting as a transcriptional repressor of ANF by binding to both GATA4 and NKX2-5 and repressing their transcriptional activator activities.
An intronic variant in the thyroid peroxidase gene (TPO) has also shown association with ASD (rs1514687, P = 5.72 × 10-06, OR = 1.46 (1.24-1.72)). This association is corroborated by those of neighboring SNPs (Figure 1). This gene has been associated with hypothyroidism . This gene encodes a membrane-bound glycoprotein which acts as an enzyme and plays a central role in thyroid gland function. The protein functions in the iodination of tyrosine residues in thyroglobulin and phenoxy-ester formation between pairs of iodinated tyrosines to generate the thyroid hormones, thyroxine and triiodothyronine. Mutations in this gene are associated with several disorders of thyroid hormonogenesis, including congenital hypothyroidism, congenital goiter and thyroid hormone organification defect IIA. Multiple transcript variants encoding distinct isoforms have been identified for this gene.
In spite of its reported high heritability, the identification of genetic risk factors predisposing to ASDs has been difficult. To date, only a small number of loci are considered established (that is, have met genome-wide significance and have been replicated in independent cohorts). GWAS are by design agnostic. Increased statistical power can be achieved with more focused hypotheses. Given the increasing speculation of a significant role for immune genes in the etiology and pathogenesis of ASD, we undertook a search for genetic variation associated with ASD in genes with immune functions. This hypothesis-driven candidate gene approach based on prior knowledge was selected because it serves to increase the power, reliability and likelihood of finding genes that are associated with disease [5, 6].
A lack of enrichment of genetic associations in immune-related genes has recently been reported  in the same dataset we used, and this conclusion is supported here using another statistical approach. The associations we uncovered in immune-related genes and the observed dysregulation of immune genes in ASD  suggest that variation in a limited number of immune-function genes may be responsible for observed up-regulation of their immune downstream targets.
The transmission/disequilibrium test (TDT) is largely robust to population substructure as it is a within-family comparison. The individuals in these data are predominately of European descent (76%), with Asian (3%) and African (2%) ancestry represented; the Autism Genetic Resource Exchange (AGRE) repository lists 14% and 5% of the individuals of unknown or admixed ancestry. One limitation of using the TDT in our study is the lack of total independence between the trios analyzed: the 1,510 trios used in the TDT are part of 1,101 independent nuclear families. The correlation among children in the same family can inflate the type I error probability. We, thus, also applied the PDT, also a within family test that is robust to population substructure, and the PDT analysis corroborates the TDT analysis.
SNPs within the CD99 molecule-like 2 gene, the jumonji AT rich interactive domain 2 (JARID2) gene, and the thyroid peroxidase gene were associated with ASD after multiple comparisons adjustment. Understanding how these genetic factors might contribute to pathogenesis should ultimately lead to important opportunities for discovering therapeutic targets that can be used to treat ASD.
Autism Genetic Resource Exchange
Autism Spectrum Disorders
False Discovery Rate
genome wide association studies
Ingenuity Pathway Analysis
jumonji AT rich interactive domain 2
minor allele frequency
Pedigree Disequilibrium Test
False Discovery Rate-adjusted P-value
single nucleotide polymorphism
Transmission Disequilibrium Test
thyroid peroxidase gene.
We gratefully acknowledge the resources provided by the Autism Genetic Resource Exchange (AGRE) Consortium* and the participating AGRE families. The Autism Genetic Resource Exchange is a program of Autism Speaks and is supported, in part, by grant 1U24MH081810 from the National Institute of Mental Health to Clara M. Lajonchere (PI). We also acknowledge funding from the Jane Botsford Johnson Foundation, The Thoughtful House Center for Children, and the Wake Forest School of Medicine Center for Public Health Genomics.
- CDC: Prevalence of autism spectrum disorders - Autism and Developmental Disabilities Monitoring Network, United States. Morb Mortal Wkly Rep Surveill Summ. 2012, 61: 1-19.Google Scholar
- Ashwood P, Krakowiak P, Hertz-Picciotto I, Hansen R, Pessah I, Van de Water J: Elevated plasma cytokines in autism spectrum disorders provide evidence of immune dysfunction and are associated with impaired behavioral outcome. Brain Behav Immun. 2011, 25: 40-45. 10.1016/j.bbi.2010.08.003.PubMed CentralView ArticlePubMedGoogle Scholar
- Jyonouchi H, Geng L, Streck DL, Toruner GA: Children with autism spectrum disorders (ASD) who exhibit chronic gastrointestinal (GI) symptoms and marked fluctuation of behavioral symptoms exhibit distinct innate immune abnormalities and transcriptional profiles of peripheral blood (PB) monocytes. J Neuroimmunol. 2011, 238: 73-80. 10.1016/j.jneuroim.2011.07.001.View ArticlePubMedGoogle Scholar
- Onore C, Careaga M, Ashwood P: The role of immune dysfunction in the pathophysiology of autism. Brain Behav Immun. 2012, 26: 383-392. 10.1016/j.bbi.2011.08.007.PubMed CentralView ArticlePubMedGoogle Scholar
- Armstrong DL, Reiff A, Myones BL, Quismorio FP, Klein-Gitelman M, McCurdy D, Wagner-Weiner L, Silverman E, Ojwang JO, Kaufman KM, Kelly JA, Merrill JT, Harley JB, Bae SC, Vyse TJ, Gilkeson GS, Gaffney PM, Moser KL, Putterman C, Edberg JC, Brown EE, Ziegler J, Langefeld CD, Zidovetzki R, Jacob CO: Identification of new SLE-associated genes with a two-step Bayesian study design. Genes Immun. 2009, 10: 446-456. 10.1038/gene.2009.38.PubMed CentralView ArticlePubMedGoogle Scholar
- Gaulton KJ, Willer CJ, Li Y, Scott LJ, Conneely KN, Jackson AU, Duren WL, Chines PS, Narisu N, Bonnycastle LL, Luo J, Tong M, Sprau AG, Pugh EW, Doheny KF, Valle TT, Abecasis GR, Tuomilehto J, Bergman RN, Collins FS, Boehnke M, Mohlke KL: Comprehensive association study of type 2 diabetes and related quantitative traits with 222 candidate genes. Diabetes. 2008, 57: 3136-3144. 10.2337/db07-1731.PubMed CentralView ArticlePubMedGoogle Scholar
- Ingenuity Pathway Analysis.http://www.ingenuity.com,
- Autism Genetic Resource Exchange.http://research.agre.org,
- Ma D, Salyakina D, Jaworski JM, Konidari I, Whitehead PL, Andersen AN, Hoffman JD, Slifer SH, Hedges DJ, Cukier HN, Griswold AJ, McCauley JL, Beecham GW, Wright HH, Abramson RK, Martin ER, Hussman JP, Gilbert JR, Cuccaro ML, Haines JL, Pericak-Vance MA: A genome-wide association study of autism reveals a common novel risk locus at 5p14.1. Ann Hum Genet. 2009, 73: 263-273. 10.1111/j.1469-1809.2009.00523.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang K, Zhang H, Ma D, Bucan M, Glessner JT, Abrahams BS, Salyakina D, Imielinski M, Bradfield JP, Sleiman PM, Kim CE, Hou C, Frackelton E, Chiavacci R, Takahashi N, Sakurai T, Rappaport E, Lajonchere CM, Munson J, Estes A, Korvatska O, Piven J, Sonnenblick LI, Alvarez Retuerto AI, Herman EI, Dong H, Hutman T, Sigman M, Ozonoff S, Klin A: Common genetic variants on 5p14.1 associate with autism spectrum disorders. Nature. 2009, 459: 528-533. 10.1038/nature07999.PubMed CentralView ArticlePubMedGoogle Scholar
- Weiss LA, Arking DE, Daly MJ, Chakravarti A: A genome-wide linkage and association scan reveals novel loci for autism. Nature. 2009, 461: 802-808. 10.1038/nature08490.PubMed CentralView ArticlePubMedGoogle Scholar
- Voineagu I, Wang X, Johnston P, Lowe JK, Tian Y, Horvath S, Mill J, Cantor RM, Blencowe BJ, Geschwind DH: Transcriptomic analysis of autistic brain reveals convergent molecular pathology. Nature. 2011, 474: 380-384. 10.1038/nature10110.PubMed CentralView ArticlePubMedGoogle Scholar
- Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Maller J, Sklar P, de Bakker PI, Daly MJ, Sham PC: PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007, 81: 559-575. 10.1086/519795.PubMed CentralView ArticlePubMedGoogle Scholar
- Martin ER, Monks SA, Warren LL, Kaplan NL: A test for linkage and association in general pedigrees: the pedigree disequilibrium test. Am J Hum Genet. 2000, 67: 146-154. 10.1086/302957.PubMed CentralView ArticlePubMedGoogle Scholar
- Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc B. 1995, 57: 289-300.Google Scholar
- Browning BL, Browning SR: A unified approach to genotype imputation and haplotype-phase inference for large data sets of trios and unrelated individuals. Am J Hum Genet. 2009, 84: 210-223. 10.1016/j.ajhg.2009.01.005.PubMed CentralView ArticlePubMedGoogle Scholar
- Pedrosa E, Ye K, Nolan KA, Morrell L, Okun JM, Persky AD, Saito T, Lachman HM: Positive association of schizophrenia to JARID2 gene. Am J Med Genet B Neuropsychiatr Genet. 2007, 144B: 45-51. 10.1002/ajmg.b.30386.View ArticlePubMedGoogle Scholar
- Liu Y, Chen G, Norton N, Liu W, Zhu H, Zhou P, Luan M, Yang S, Chen X, Carroll L, Williams NM, O'Donovan MC, Kirov G, Owen MJ: Whole genome association study in a homogenous population in Shandong peninsula of China reveals JARID2 as a susceptibility gene for schizophrenia. J Biomed Biotechnol. 2009, 2009: 536918-PubMed CentralPubMedGoogle Scholar
- Cloos PA, Christensen J, Agger K, Helin K: Erasing the methyl mark: histone demethylases at the center of cellular differentiation and disease. Genes Dev. 2008, 22: 1115-1140. 10.1101/gad.1652908.PubMed CentralView ArticlePubMedGoogle Scholar
- Avbelj M, Tahirovic H, Debeljak M, Kusekova M, Toromanovic A, Krzisnik C, Battelino T: High prevalence of thyroid peroxidase gene mutations in patients with thyroid dyshormonogenesis. Eur J Endocrinol. 2007, 156: 511-519. 10.1530/EJE-07-0037.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.