Context Neuronal dysfunction in cerebral cortex and other brain regions could contribute to the cognitive and behavioral defects in autism.
Objective To characterize epigenetic signatures of autism in prefrontal cortex neurons.
Design We performed fluorescence-activated sorting and separation of neuronal and nonneuronal nuclei from postmortem prefrontal cortex, digested the chromatin with micrococcal nuclease, and deeply sequenced the DNA from the mononucleosomes with trimethylated H3K4 (H3K4me3), a histone mark associated with transcriptional regulation. Approximately 15 billion base pairs of H3K4me3-enriched sequences were collected from 32 brains.
Setting Academic medical center.
Participants A total of 16 subjects diagnosed as having autism and 16 control subjects ranging in age from 0.5 to 70 years.
Main Outcome Measures Identification of genomic loci showing autism-associated H3K4me3 changes in prefrontal cortex neurons.
Results Subjects with autism showed no evidence for generalized disruption of the developmentally regulated remodeling of the H3K4me3 landscape that defines normal prefrontal cortex neurons in early infancy. However, excess spreading of H3K4me3 from the transcription start sites into downstream gene bodies and upstream promoters was observed specifically in neuronal chromatin from 4 of 16 autism cases but not in controls. Variable subsets of autism cases exhibit altered H3K4me3 peaks at numerous genes regulating neuronal connectivity, social behaviors, and cognition, often in conjunction with altered expression of the corresponding transcripts. Autism-associated H3K4me3 peaks were significantly enriched in genes and loci implicated in neurodevelopmental diseases.
Conclusions Prefrontal cortex neurons from subjects with autism show changes in chromatin structures at hundreds of loci genome-wide, revealing considerable overlap between genetic and epigenetic risk maps of developmental brain disorders.
Autism spectrum disorders comprise a group of complex and etiologically heterogeneous illnesses. The genetic risk architecture remains unknown for most patients, and fewer than 10% of subjects on the autism spectrum harbor rare structural DNA variations and mutations with strong penetrance. Neurons residing in the prefrontal cortex (PFC) and other cortical association areas in autistic subjects are affected by subtle defects in connectivity patterns, cytoarchitecture, and other structural alterations.1-3 While the role of such disordered neural circuitry in the context of neurodevelopmental disease is well understood, including impairment of higher cognition and social communication, very little is known about the molecular pathology of PFC neurons in the autistic brain.
Epigenetic dysregulation of DNA methylation and histone modifications could play a prominent role in the pathophysiology of autism and related disease.4-11 This hypothesis is in part based on the link between autism and deleterious mutations in genes considered of pivotal importance for regulation of chromatin structure and function, such as MECP27 and other methyl-CpG-binding proteins,12,13 the histone deacetylase HDAC4,14 the histone H3 lysine 9 (H3K9)–specific methyltransferase EHMT1/KMT1D/GLP,15 and the H3K4 demethylase JARID1C/KDM5C/SMCX.16 In particular, the trimethylated form of H3K4 (H3K4me3) is primarily located at transcription start sites (TSSs) and linked to the serine 5-phosphorylated, initiation form of RNA polymerase II, thereby providing a docking site at the 5′ end of genes for chromatin remodeling complexes that mostly facilitate (but in some cases repress) transcription.17,18 Epigenetic fine-tuning of H3K4me3 appears to be particularly important for neuronal health. For example, neuronal differentiation is dependent on H3K4 trimethylation mediated by the mixed-lineage leukemia (MLL) methyltransferase and transcriptional activation of RE-1 silencing transcription factor (REST)–sensitive genes.19,20 Furthermore, hippocampal learning and memory require MLL1 -mediated H3K4 trimethylation at growth- and plasticity-regulating genes.21,22
Recently, we presented the first neuronal-specific epigenomes from human brain and provided evidence that nuclei of PFC neurons from normal infants (aged<1 year) showed an excess of H3K4me3 at hundreds of loci compared with older brains.23 This large-scale remodeling of the H3K4me3 landscape in chromatin of young neurons speaks for epigenetic vulnerability of the immature PFC, an intriguing hypothesis when viewed from the neurodevelopmental perspective of autism. Therefore, the aim of this study was to compare, for the first time to our knowledge, the H3K4me3 epigenomes of the PFC neurons of autistic individuals with a panel of controls across a wide age range from infancy to 70 years. Our findings indicate that a subset of autistic individuals are affected by loss or excess of H3K4me3 at hundreds of loci, in conjunction with dysregulated expression of transcripts implicated in neuronal communication and social and other higher-order behaviors. We propose a new disease model for PFC neurons that involves changes in the H3K4me3 landscape in hundreds of loci. These changes are highly variable between different patients and point to a complex interaction between the genetic and epigenetic risk architectures that affect autism spectrum disorders.
Chromatin immunoprecipitation with deep sequencing
The eAppendix has a detailed description of the methods and techniques used. Sorted neuronal (NeuN+) and nonneuronal (NeuN−) nuclei were processed for anti-H3K4me3 chromatin immunoprecipitation (ChIP), and resulting DNA libraries were deep sequenced with an Illumina GAII platform (Illumina, Inc, San Diego, California). Quantification of uniquely mappable sequence tags was based on MACS (Model-Based Analysis of ChIP-Seq) version 1.3.5 software,24 and signals were expressed as parts (tags) per million. Additional bioinformatics approaches to analyze disease-associated changes are provided in the eAppendix.
GENOME-WIDE MAPPING OF H3K4me3 IN HUMAN PFC NEURONS
We previously developed a technique of fluorescence-activated cell sorting of nuclei according to immunoreactivity for the NeuN antigen.25 This allows us to compare changes in the histone methylation landscapes of terminally differentiated neurons across the lifespan, including developmental periods and disease conditions such as autism,26,27 without being affected by shifts in the neuron-glia ratio, which could greatly confound interpretation of conventional chromatin studies from tissue homogenates. We then digested chromatin with micrococcal nuclease, extracted the mononucleosome fraction, and performed ChIP with an antibody against H3K4me3 followed by deep sequencing. We obtained genome-wide maps of the H3K4me3 mark for NeuN+ nuclei from the PFC of 16 subjects diagnosed as having autism spectrum disorders (aged 2-60 years; mean age, 17.6 years; 4 females; supplemental Table S1, http://zlab.umassmed.edu/zlab/publications/ShulhaAGP2011.html) and 10 age-matched control subjects (aged 2.8-69 years; mean age, 19.7 years; 3 females). In addition, we included 4 control subjects younger than 2 years (1 female) to investigate whether H3K4me3 profiles of autistic subjects were more similar to those of normal infants than to those of older control subjects. Nine of the 14 control samples were generated in our earlier study.23 We also sequenced an input library from the neuronal nuclei of 1 control subject, following all the steps of ChIP with deep sequencing but omitting the anti-H3K4me3 antibody. Altogether, 333 million 36-nucleotide reads were obtained for 31 neuronal samples, among which 283 million reads mapped to unique locations in the reference human genome and an additional 32 million reads mapped to multiple locations in the genome. Furthermore, we performed ChIP followed by deep sequencing on NeuN− chromatin of 4 of the autism cases and 2 additional control subjects (supplemental Table S2).
GENOME-WIDE CORRELATION OF H3K4me3 OCCUPANCY IN ANNOTATED PROMOTERS
Consistent with the knowledge that H3K4me3 is a histone mark sharply enriched around the 5′ end of genes, a mean (SD) of 59% (9%) of reads in the 14 control samples mapped within RefSeq promoters (http://www.ncbi.nlm.nih.gov/RefSeq/). Similarly, a mean (SD) of 57% (11%) of reads in the autism samples mapped in proximal promoters (supplemental Table S1). There is no statistical difference in the percentage of tags mapping to promoters between autism and control samples, indicating that autistic patients do not show global displacement of H3K4me3 away from promoters. In sharp contrast, fewer than 4% of reads from the input library mapped to promoters. Postmortem confounds—tissue pH and postmortem interval—do not show significant correlations with the percentage of TSS-proximal tags (correlation coefficients,−0.12 and−0.19, respectively; t tests, 2-tailed P = .55 and .31, respectively).
To investigate the global H3K4me3 occupancy in promoters, we computed the correlation coefficient for the numbers of tags in the entire set of 21 084 annotated promoters between any 2 samples. Figure 1 shows that all 30 samples are highly similar, with correlations ranging from 0.93 to 0.99. The autism samples do not show appreciable difference from controls: the mean (SD) correlation between 16 autism samples and 10 controls of the same age range is 0.969 (0.011), while the mean (SD) correlations are 0.971 (0.012) within autism samples and 0.974 (0.008) within controls. The most distinct samples are from the 3 youngest infants C1-C3 (aged<1 year), indicated by the dark shades of the first 3 columns in Figure 1. They are equally distinct from autism samples and older control subjects; the mean (SD) correlation coefficients are 0.951 (0.009) and 0.954 (0.009), respectively. Our results suggest that autistic patients do not exhibit global difference in H3K4me3 occupancy at annotated promoters. Furthermore, all autism samples included in this study appear to have an age-appropriate promoter H3K4me3 landscape, more similar to control subjects older than 2 years than to control infants.
H3K4me3 PROFILES AROUND THE TSS
As reported previously, the nucleosomes around TSSs are highly enriched in H3K4me3, with the nucleosome immediately downstream from the TSS (the +1 nucleosome) showing the strongest signal.28 Indeed, we observed this phenomenon and further noticed that some TSSs have strong H3K4me3 signals in multiple downstream nucleosomes. To identify these TSSs systematically, we performed K-means clustering (K = 7) on all TSSs of a control sample (C7) according to the H3K4me3 profiles in a±2–kilobase (kb) window centered on the TSS (clustering with other control samples led to highly similar results). Figure 2A shows the heatmap of the resulting clusters sorted inversely by average H3K4me3 signal.
We analyzed whether genes in each cluster were enriched in Gene Ontology terms. Significantly enriched Gene Ontology terms are shown in supplemental Table S3. Cluster 1 contains genes with the strongest average H3K4me3 signal. In addition, genes in cluster 1 exhibit the largest number of H3K4me3-enriched nucleosomes compared with genes in clusters 2 through 6. Six hundred two genes in cluster 1 were among the 1721 genes we previously defined by contrasting the levels of H3K4me3 in PFC neurons and in lymphocytes,23 and the overlap of the 2 sets of genes is highly significant (χ2 test, P < 2.2 × 10−16). Accordingly, genes in cluster 1 are highly enriched in neuronal functions. The most enriched Gene Ontology terms of these genes include nervous system development, generation of neurons, synapse, etc (false discovery rate < 3.8 × 10−16). In contrast, clusters 2 through 6 are enriched in various housekeeping Gene Ontology terms, while cluster 7, enriched in immune response, etc, showed only very low H3K4me3 levels in PFC neurons.
We plotted the H3K4me3 heatmaps of the other 29 samples with TSSs in the same order as in Figure 2A. To compare H3K4me3 profiles across individuals, we generated line graphs, 1 line per individual, representing average H3K4me3 profiles of all genes in each cluster (Figure 2B). Four autism samples (A5, A8, A9, and A13) stood out as showing H3K4me3 profiles spreading away from the TSSs and encroaching into gene bodies and upstream sequences. The heatmap of sample A5 is shown in Figure 2C, and the remaining 28 samples are shown in supplemental Figure S1. No clinical variables appeared to be common to these 4 cases (supplemental Table S1). Figure 2B (red lines) clearly shows that the average H3K4me3 profiles from these 4 samples spread and plateau beneath maximal peak heights of all controls and the remaining autism samples. These spreading profiles are seen for genes in clusters 1 through 6 but not for those in cluster 7 (Figure 2C) and include multiple susceptibility genes previously implicated in autism (Figure 3 and supplemental Table S4).
To evaluate the significance of spreading, we used cluster 1 as the representative and tested the average H3K4me3 levels at 4 positions sampled around the H3K4me3 peak:−2 kb, 0 kb, +1 kb, and +3 kb from the TSS. There was a significant difference between the 4 spreading autism cases and the remaining 12 cases (t test, 2-tailed P = 4.0 × 10−3, 9.5 × 10−4, 9.5 × 10−6, and .01) and all control cases (P = 4.3 × 10−3, 5.0 × 10−3, 8.4 × 10−5, and .02) but not between the 12 nonspreading cases and controls (P > .12). Spreading of the H3K4me3 signal in cases A5, A8, A9, and A13 was not due to an unbalanced number of tags mapping to the 2 genomic strands, because a mean (SD) of 49.97% (0.02%) of genome mapping reads map to the Watson (+) strand in all cases and controls (supplemental Table S5).
Spreading can lead to additional H3K4me3 peaks (genomic regions enriched in H3K4me3 are called H3K4me3 peaks). For example, case A5 shows an additional peak near the NRCAM gene, which is absent from all other samples (Figure 3). We performed ChIP–quantitative polymerase chain reaction on this peak for case A5 and 3 controls (C8, C10, and C12), using the housekeeping gene B2M as control. Indeed, we detected an elevated H3K4me3 signal in case A5: log10(NRCAM /B2M) = −3.67, compared with controls C8 (−4.54), C10 (−3.91), and C12 (−4.29). Excess spreading of H3K4me3 was further verified in cases A5 and A13 by probing AHCYL1 TSSs on chromosome 1 with ChIP–quantitative polymerase chain reaction (supplemental Figure S2).
To test whether the spreading was a general phenomenon affecting all PFC cells in the 4 autism cases, we conducted additional experiments of H3K4me3 ChIP with deep sequencing on NeuN− nuclei of cases A5, A8, A9, and A13 and 2 additional controls (supplemental Table S2) chosen to match the ages of the 4 spreading autism cases. The NeuN− nuclei are primarily from nonneuronal cells, including glia, microglia, and endothelium. The H3K4me3 heatmaps of these 6 samples with TSSs in the same order as in Figure 2A, as well as the line graphs of the 7 clusters as defined in Figure 2, are plotted in supplemental Figure S3. The H3K4me3 profiles of the NeuN− chromatin from cases A5, A8, A9, and A13 are clustered around those of the 2 control samples and do not show any spreading behavior (supplemental Figure S3). We conclude that neuronal chromatin from a subset of autism cases exhibits excessive spreading of the H3K4me3 mark into nucleosomes positioned farther away from the TSS.
FINGERPRINTING SAMPLES USING GENOME-WIDE H3K4me3 PEAKS
Having analyzed H3K4me3 occupancies (Figure 1) and profiles (Figure 2) at annotated promoters, we next asked whether the autism samples differ from controls in nonpromoter regions. We performed unsupervised learning using H3K4me3-based fingerprints computed as follows. For each sample, we determined the total number of genomic regions (not restricted to promoters) enriched in H3K4me3 compared with another sample as the background. Such regions are called H3K4me3 peaks (details in the eAppendix). We used all 30 samples (16 autism samples and 14 control samples); thus, the fingerprint of each sample is a vector of 30 numbers, each being the number of H3K4me3 peaks in the sample using another sample as the background.
We performed principal components analysis on the fingerprints of the 30 samples. Principal components analysis is a mathematical method that reduces the dimensionality of data by identifying directions, called principal components, that maximize the variation in the data. Samples can then be plotted along these principal components to visually assess whether they form clusters. Principal components analysis is an unsupervised learning method, unaware of sample identities (autism or control). Figure 4 shows our 30 samples plotted against the first 3 principal components. Strikingly, the 4 youngest control subjects, C1 through C4, ranging in age from 0.5 to 1.3 years, are located at one end of the graph, far from the space densely occupied by the older control subjects and most autism cases. The next 3 youngest control subjects (C5-C7) and the youngest autistic subject (A1; aged 2 years) bridge the space that separates infants from the remaining subjects. We can draw 2 conclusions from these results. First, the H3K4me3 landscape in PFC neurons undergoes significant remodeling during at least the first 12 months after birth, as previously reported in a smaller cohort.23 Second, no autistic subjects of our study showed evidence for developmental arrest in this preprogrammed remodeling of PFC neuron chromatin.
H3K4me3 PEAKS ENRICHED OR DEPLETED IN SUBSETS OF AUTISM SAMPLES
Autism is a brain disorder with considerable disease heterogeneity in terms of clinical presentation, molecular and cellular pathology, and genetics.29 Therefore, we asked whether H3K4me3 alterations, if present in our clinical cohort, show interindividual variability and affect multiple cases. To explore this, we screened the H3K4me3 data sets from the 16 autistic subjects for peaks that consistently showed greater than a 2-fold upregulation or downregulation in comparison with each of the 10 controls. We thus identified 503 loci as increased and 208 loci as decreased in a subset of autistic subjects (supplemental Table S4, Figure S4, and Figure S5). We call these autism-up and autism-down H3K4me3 loci, respectively, and 330 autism-up loci and 139 autism-down loci overlap with promoters. We then applied Poisson statistics; the significant peaks (P < .001 after Bonferroni multiple-testing correction) are included in supplemental Table S4. Most (but not all) of the autism-up loci were derived from the 4 cases with spreading H3K4me3 profiles (cases A5, A8, A9, and A13 described earlier); furthermore, cases A5 and A13 had the largest shares of the autism-down loci (supplemental Figure S4 and Figure S5). We further validated the autism-associated H3K4me3 changes of 2 more genes by quantitative polymerase chain reaction (supplemental Figure S2).
The H3K4me3 functions as a docking site for a diverse set of transcriptional regulators that could activate or repress gene expression. It has been shown that developmental or disease-associated changes in brain histone methylation are linked to alterations of the corresponding gene transcripts.30-32 Therefore, we asked whether autism cases with altered H3K4me3 levels also showed altered levels of the corresponding messenger RNA (mRNA). We chose 5 genes with decreased H3K4me3 over the TSSs in at least 2 autism cases from a larger set of genes proximal to autism-up H3K4me3 loci, which are implicated in neurodevelopment and cognition based on preclinical and/or clinical research (supplemental Table S4). The following 5 genes were included: (1) the synaptic protein ARC, which interacts with the ubiquitin ligase and autism/Angelman syndrome risk gene UBE3A to regulate synaptic AMPA glutamate receptor trafficking33; (2) complexin 1 (CPLX1), a protein involved in exocytosis and synaptic vesicle functions, which is essential for social behavior and motivational behaviors in mice and is expressed at altered levels in various neuropsychiatric diseases34; (3) the orphan nuclear receptor NR4A1 genetically linked to involuntary movement disorders35; (4) USPx, a dendritic protein implicated in neuronal and synaptic connectivity36; and (5) the neuropeptide VGF, which regulates neuronal plasticity and is expressed at decreased levels in postmortem brain of subjects diagnosed as having manic-depressive illness.37 Strikingly, all autistic subjects with greater than a 2-fold decrease in TSS-associated H4K4me3 peak density also exhibited some of the lowest transcript levels in comparison with controls and the remaining autism cases (Figure 5A and supplemental Figure S6). These changes were highly specific because the nearest neighboring peaks were not altered in the affected autism cases. Furthermore, RNA quantifications were similar independent of the choice of the housekeeping or normalization gene(s) (supplemental Figure S6). In comparing demographic, clinical, and postmortem indices between cases with a decline in specific H3K4me3 peaks and levels of corresponding mRNAs, we did not observe any factor (including seizures, medication histories, or postmortem confounds) as being common to these cases other than the diagnosis of autism (supplemental Table S1).
Next, we probed RNA levels for 3 autism-up loci. For 2 of the 3 transcripts tested, IFI6, which encodes inducible interferonα6, and ALDH1A2, which is a schizophrenia susceptibility gene linked to retinoic acid signaling,38 cases with high H3K4me3 signal tended to show a robust increase in mRNA levels (Figure 5B and supplemental Figure S6). We conclude that in PFC neurons, aberrant H3K4 methylation at a specific TSS indeed is a strong predictor for transcriptional dysregulation, which affects variable subsets of subjects on the autism spectrum.
Overlap with previously annotated autism risk loci
We then asked whether some of the 711 loci with altered H3K4me3 signal (the aforementioned autism-up and autism-down peaks) were in the vicinity of any of the established autism susceptibility loci or genes with strong penetrance of disease risk when mutated. We examined 3 databases: the Simons Foundation Autism Research Initiative database (SFARI) for autism risk genes,39 the Database of Chromosomal Imbalance and Phenotype in Humans Using Ensembl Resources (DECIPHER) for copy number variations (CNVs) associated with diseased neurodevelopment,40 and the Human Unidentified Gene-Encoded protein database (HuGE) for autism-associated polymorphisms.41 Indeed, 89 of 729 loci matched loci in these databases. The list included at least 10 autism risk genes (ASTN2, CACNA1C, CACNA1H, JMJD1C, MEF2C, NRCAM, RAI1, RIMS3, RYR2, and SEM A5A) with abnormal H3K4me3 levels in their promoters and additional disease genes showing H3K4me3 changes in other portions of the gene, such as AUTS2, PARK2, and VGF (supplemental Table S4). For the various subsets of autism cases affected by the epigenetic dysregulation at these genetic risk loci, the altered H3K4me3 profile was specific for neuronal chromatin from PFC, while NeuN− chromatin from PFC appeared unaffected (Figure 3). For the entire set of 16 autism cases, enrichment in the overlap between the autism-up and autism-down peaks and the loci in DECIPHER (58 loci in total) is 1.87-fold and 2.44-fold, respectively.40 Examples include portions of 1p36 associated with one of the most frequent genetic causes of non–X-linked mental retardation42 and a portion of 17p associated with deletion (Smith-Magenis) or duplication (Potocki-Lupski) syndromes conferring developmental disability and autism43 (supplemental Table S4). There is also a significant overlap between the autism-up H3K4me3 peaks and autism-implicated genes in SFARI (χ2 test, P = .009) as well as between autism-down H3K4me3 peaks and autism-implicated genes in HuGE (P = .04). We conclude that epigenetic alterations in PFC neurons of autistic subjects include a significant subset of previously identified disease-associated genes.
Because our mapping algorithm allowed for 1 mismatch per 36-nucleotide read, it is unlikely that the changes in H3K4me3 profiles at these autism risk genes are due to single-nucleotide polymorphisms of the H3K4me3-enriched DNA. To explore whether the observed H3K4me3 changes in our autism samples are the result of CNVs, we screened published data sets from 8 of our autism cases for microdeletions, microduplications, and other CNVs44 (cases A1, A9, A10, and A12-16). Remarkably, from 178 H3K4me3 peaks with a decrease or loss in various subsets of these 8 cases, only 4 peaks (2.2%) were located to a microdeletion. The deletion in case A16 involved 22q11.23, a region frequently affected by segmental duplication upstream of the 22q11 DiGeorge syndrome/velocardiofacial syndrome risk locus45 (supplemental Figure S7). Two cases (A14 and A16) were affected by a partial loss of 1p36.33, representing one form of a subtelomeric deletion syndrome (supplemental Table S4). In most of the affected cases, the syndrome (with a prevalence of approximately 1 in 10 000 individuals) involves several hundred kb or megabases (Mb) of gene-rich portions within the subtelomeric region of chromosome 1 and frequently results in developmental delay and neurological disease.42,46 In our cohort of 16 autistic subjects, we detected 14 autism-up and 8 autism-down peaks in a 28-Mb portion of 1p36.32-33, including multiple cases with no detectable CNV at these loci (Figure 6). Furthermore, none of the 342 peaks with excess H3K4me3 (supplemental Table S4) were associated with CNVs. We conclude that structural variations including CNVs are unlikely to play a role for most of the observed H3K4me3 changes in our cohort of autism cases.
To confirm that some of the H3K4me3 changes in these autism risk genes are associated with dysregulated gene expression, we measured the RNA for RIMS3, which encodes a synaptic vesicle protein associated with very rare, autism-associated structural DNA variations, including the 1p34.2 microdeletion and, in some families, changes in coding exon sequences.47RIMS3 also shows altered expression in lymphoblastoid cell lines derived from autistic patients.48 Indeed, like the 5 previously described genes with a robust decrease in H3K4me3, autistic subjects with the lowest H3K4me3 signal at the RIMS3 TSS also expressed low levels of RIMS3 transcript (supplemental Figure S6).
Chromatin changes in autism
To our knowledge, this study provides the first insights into chromatin structures of PFC neurons from subjects in the autism spectrum. Remarkably, none of our 16 disease cases showed evidence for a generalized disruption of the genome-wide redistribution process of H3K4me3 that defines the normal transition from early infancy (age<1 year after birth) to older ages. Instead, we identified 503 increased loci and 208 decreased loci in the genome, including but not limited to annotated TSSs, that were affected by at least a 2-fold increase or decrease in H3K4me3 peak densities in subsets of autism cases compared with controls. While each of these altered loci affected variable subsets of cases, the majority of them are from 4 diseased brains with an abnormal H3K4me3 profile around the TSS, defined by lower H3K4me3 intensity at the TSS and spreading away from the TSS. Because diagnosis of autism was the only identifiable factor common to these cases, it is likely that this abnormal H3K4me3 profile is related to autism. One spreading case and 1 nonspreading case were exposed to valproic acid, a drug that upregulates histone acetylation in the brain (supplemental Table S1).31,49,50 Nonetheless, additional work with larger cohorts will be necessary to clarify the biological significance of the H3K4me3 spreading and the implications for overall nucleosomal organization around the TSS.
Our observations add further evidence that epigenetic fine-tuning of H3K4 methylation is particularly important for orderly brain development and function19-22,31 and suggest maladaptive chromatin remodeling with imbalances in H3K4me3 distribution as a molecular phenotype of cortical neurons in autism. Indeed, several lines of evidence from this study suggest that these alterations in H3K4me3 profiles of PFC neurons are directly involved in the molecular pathology of the disease in the affected individuals. For example, as exemplified by ARC, CPLX1, USPx, RIMS3, and VGF— genes with a key role in the neurobiology of cognition and higher-order behaviors—most of the subjects affected by a robust decrease in H3K4me3 at the specific TSSs also showed a deficit in levels of the corresponding transcripts. Likewise, increased TSS-bound H3K4me3 was in some cases accompanied by a dramatic increase in transcript levels (compare IFI6 RNA and H3K4me3 in subject A6, Figure 5 and supplemental Figure S6). Given that the focus of this study was on neuronal chromatin, it might have been expected that immune-related genes did not play a major role. Alterations in the immune system of autistic subjects51 could reflect activation of microglia and other nonneuronal constituents residing in diseased tissue.27,52
The finding that genes and loci previously shown to confer genetic risk for autism and related disease significantly overlap with the 711 differential H3K4me3 peaks in our disease cohort implies that there is significant overlap between the genetic and epigenetic risk architectures in autism. It is noteworthy that these epigenetic changes are not limited to the 4 brains with the aforementioned spreading H3K4me3 profile, because the remaining 12 autism samples exhibit altered H3K4me3 at many of the same loci. There are 3 possible explanations for the H3K4me3 changes observed in subsets of autism samples: (1) variation in the underlying genomic DNA, which as discussed earlier is unlikely to account for a significant portion of the disease-associated H3K4me3 peaks given that the genome mapping algorithm allowed only a single base mismatch per tag; (2) heritable transmission of histone modifications53; and (3) a molecular adaptation in a pathophysiological cascade ultimately leading to cortical dysfunction.
Genetic models for an increasing number of neuropsychiatric disorders, including autism, invoke a complex mixture of disease-associated common variants of a select number of genes, each contributing to a small fraction of the disease risk, and rare variations such as microdeletions and duplications of key genomic loci and candidate genes that, when affected, carry a high penetrance of disease.29,54,55 Heterogeneity of disease was also apparent in this study because the 711 epigenetic risk loci were affected in variable subsets of autistic individuals. Of note, the list of H3K4me3 peaks in our small cohort of 16 autistic subjects included the synaptic vesicle gene RIMS3, the retinoic acid signaling–regulated gene RAI1,56 the histone demethylase JMJD1C,57 the astrotactin ASTN2,58 the adhesion molecules NRCAM59-61 and SEMA5,62,63 the ubiquitin ligase PARKIN2 (PARK2),58,64 and other genes for which rare structural DNA variations could carry high disease penetrance. Intriguingly, for most of these high-risk genes, we observed that H3K4me3 changes occurred selectively in PFC neurons but not their surrounding nonneuronal cells (Figure 3). Our findings imply a new disease model suggesting that a subset of susceptibility genes carrying strong penetrance for disease could be epigenetically dysregulated in a cell (type)–specific manner, independent of the occurrence of any rare structural DNA variations so far implicated in neurological disease. Identification of the precise combinations of these aberrantly regulated gene sets in a particular disease case, including a potential association with the severity of specific disease symptoms, can bear great promise toward better understanding of the underlying neurobiology and, perhaps more importantly, could set the foundation for novel therapies specifically tailored toward individual patients.
None of the 16 autism cases in our cohort showed evidence for a generalized disruption of developmentally regulated chromatin remodeling as observed in (normal) infant PFC neurons. Hundreds of loci across the genome show disease-associated alterations in H3K4me3 levels. These defects occur in variable subsets of autistic individuals, are often highly specific for PFC neurons, and are frequently accompanied by dysregulated expression of the corresponding gene transcripts in the same subjects or brain tissues. Autism-associated H3K4me3 alterations are enriched for susceptibility loci conferring genetic risk for neurodevelopmental disease. Fewer than 5% of autism-associated H3K4me3 changes are directly related to CNVs. A subset of subjects with autism showed evidence for a more generalized alteration of their neuronal H3K4me3 landscapes as evidenced by excess spreading around many TSSs. The biological significance of this phenomenon requires further investigation.
Correspondence: Zhiping Weng, PhD, Program in Bioinformatics and Integrative Biology, University of Massachusetts Medical School, 364 Plantation St, LRB 1010, Worcester, MA 01605 (zhiping.weng@umassmed.edu).
Submitted for Publication: June 10, 2011; final revision received August 19, 2011 accepted September 11, 2011.
Published Online: November 7, 2011. doi:10.1001/archgenpsychiatry.2011.151
Author Contributions: Drs Shulha and Cheung contributed equally. Drs Akbarian and Weng are co–senior authors, had full access to all of the data in the study, and take responsibility for the integrity of the data and the accuracy of the data analysis.
Financial Disclosure: None reported.
Funding/Support: This work was supported by Autism Speaks, the International Mental Health Research Organization, the National Alliance for Research on Schizophrenia and Depression, and grants 5 R01MH071476 and RC1MH088047 from the National Institute of Mental Health (Dr Akbarian). Dr Weng was supported in part by grant DBI-0850008 from the National Science Foundation.
Role of the Sponsors: The sponsors had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; or preparation, review, or approval of the manuscript.
Additional Information: This article is dedicated to the memory of Edward G. Jones, MD, PhD.
Additional Contributions: Ron Zielke, MD, and staff from the Brain and Tissue Bank of the University of Maryland, Francine M. Benes, MD, PhD, and staff from the Harvard Brain Tissue Resource Center and the Autism Tissue Program (under director Jane Pickett, MD), William E. Bunney Jr, MD, from the University of California, Irvine, and Edward G. Jones, MD, PhD, from the University of California, Davis, supplied some of the postmortem brain tissue used in this study, and Ellen Kittler, PhD, and Maria Zapp, PhD, from the University of Massachusetts Medical School Deep Sequencing Core, Richard Konz, PhD, and staff from the University of Massachusetts Medical School Flow Cytometry Core, and Eric Morrow, PhD, MD, provided helpful discussions.
1.Zikopoulos B, Barbas H. Changes in prefrontal axons may disrupt the network in autism.
J Neurosci. 2010;30(44):14595-1460921048117
PubMedGoogle ScholarCrossref 2.Simms ML, Kemper TL, Timbie CM, Bauman ML, Blatt GJ. The anterior cingulate cortex in autism: heterogeneity of qualitative and quantitative cytoarchitectonic features suggests possible subgroups.
Acta Neuropathol. 2009;118(5):673-68419590881
PubMedGoogle ScholarCrossref 3.Casanova MF, Buxhoeveden DP, Brown C. Clinical and macroscopic correlates of minicolumnar pathology in autism.
J Child Neurol. 2002;17(9):692-69512503647
PubMedGoogle ScholarCrossref 4.Thatcher KN, LaSalle JM. Dynamic changes in histone H3 lysine 9 acetylation localization patterns during neuronal maturation require MeCP2.
Epigenetics. 2006;1(1):24-3117464364
PubMedGoogle ScholarCrossref 5.Fradin D, Cheslack-Postava K, Ladd-Acosta C, Newschaffer C, Chakravarti A, Arking DE, Feinberg A, Fallin MD. Parent-of-origin effects in autism identified through genome-wide linkage analysis of 16 000 SNPs.
PLoS One. 2010;5(9):e1251320824079
PubMedGoogle ScholarCrossref 6.Akbarian S, Huang HS. Epigenetic regulation in human brain-focus on histone lysine methylation.
Biol Psychiatry. 2009;65(3):198-20318814864
PubMedGoogle ScholarCrossref 8.Bale TL, Baram TZ, Brown AS, Goldstein JM, Insel TR, McCarthy MM, Nemeroff CB, Reyes TM, Simerly RB, Susser ES, Nestler EJ. Early life programming and neurodevelopmental disorders.
Biol Psychiatry. 2010;68(4):314-31920674602
PubMedGoogle ScholarCrossref 9.Iwamoto K, Kato T. Epigenetic profiling in schizophrenia and major mental disorders.
Neuropsychobiology. 2009;60(1):5-1119684418
PubMedGoogle ScholarCrossref 10.Grayson DR, Chen Y, Dong E, Kundakovic M, Guidotti A. From trans-methylation to cytosine methylation: evolution of the methylation hypothesis of schizophrenia.
Epigenetics. 2009;4(3):144-14919395859
PubMedGoogle ScholarCrossref 11.Petronis A. Epigenetics as a unifying principle in the aetiology of complex traits and diseases.
Nature. 2010;465(7299):721-72720535201
PubMedGoogle ScholarCrossref 12.Cukier HN, Rabionet R, Konidari I, Rayner-Evans MY, Baltos ML, Wright HH, Abramson RK, Martin ER, Cuccaro ML, Pericak-Vance MA, Gilbert JR. Novel variants identified in methyl-CpG-binding domain genes in autistic individuals.
Neurogenetics. 2010;11(3):291-30319921286
PubMedGoogle ScholarCrossref 13.Li H, Yamagata T, Mori M, Yasuhara A, Momoi MY. Mutation analysis of methyl-CpG binding protein family genes in autistic patients.
Brain Dev. 2005;27(5):321-32515967618
PubMedGoogle ScholarCrossref 14.Williams SR, Aldred MA, Der Kaloustian VM, Halal F, Gowans G, McLeod DR, Zondag S, Toriello HV, Magenis RE, Elsea SH. Haploinsufficiency of
HDAC4 causes brachydactyly mental retardation syndrome, with brachydactyly type E, developmental delays, and behavioral problems.
Am J Hum Genet. 2010;87(2):219-22820691407
PubMedGoogle ScholarCrossref 15.Balemans MC, Huibers MM, Eikelenboom NW, Kuipers AJ, van Summeren RC, Pijpers MM, Tachibana M, Shinkai Y, van Bokhoven H, Van der Zee CE. Reduced exploration, increased anxiety, and altered social behavior: autistic-like features of euchromatin histone methyltransferase 1 heterozygous knockout mice.
Behav Brain Res. 2010;208(1):47-5519896504
PubMedGoogle ScholarCrossref 16.Adegbola A, Gao H, Sommer S, Browning M. A novel mutation in
JARID1C/SMCX in a patient with autism spectrum disorder (ASD).
Am J Med Genet A. 2008;146A(4):505-51118203167
PubMedGoogle ScholarCrossref 17.Shilatifard A. Molecular implementation and physiological roles for histone H3 lysine 4 (H3K4) methylation.
Curr Opin Cell Biol. 2008;20(3):341-34818508253
PubMedGoogle ScholarCrossref 18.Guenther MG, Levine SS, Boyer LA, Jaenisch R, Young RA. A chromatin landmark and transcription initiation at most promoters in human cells.
Cell. 2007;130(1):77-8817632057
PubMedGoogle ScholarCrossref 19.Wynder C, Hakimi MA, Epstein JA, Shilatifard A, Shiekhattar R. Recruitment of MLL by HMG-domain protein iBRAF promotes neural differentiation.
Nat Cell Biol. 2005;7(11):1113-111716227968
PubMedGoogle ScholarCrossref 20.Lim DA, Huang YC, Swigut T, Mirick AL, Garcia-Verdugo JM, Wysocka J, Ernst P, Alvarez-Buylla A. Chromatin remodelling factor
Mll1 is essential for neurogenesis from postnatal neural stem cells.
Nature. 2009;458(7237):529-53319212323
PubMedGoogle ScholarCrossref 21.Gupta S, Kim SY, Artis S, Molfese DL, Schumacher A, Sweatt JD, Paylor RE, Lubin FD. Histone methylation regulates memory formation.
J Neurosci. 2010;30(10):3589-359920219993
PubMedGoogle ScholarCrossref 22.Kim SY, Levenson JM, Korsmeyer S, Sweatt JD, Schumacher A. Developmental regulation of Eed complex composition governs a switch in global histone modification in brain.
J Biol Chem. 2007;282(13):9962-997217259173
PubMedGoogle ScholarCrossref 23.Cheung I, Shulha HP, Jiang Y, Matevossian A, Wang J, Weng Z, Akbarian S. Developmental regulation and individual differences of neuronal H3K4me3 epigenomes in the prefrontal cortex.
Proc Natl Acad Sci U S A. 2010;107(19):8824-882920421462
PubMedGoogle ScholarCrossref 24.Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, Nusbaum C, Myers RM, Brown M, Li W, Liu XS. Model-based analysis of ChIP-Seq (MACS).
Genome Biol. 2008;9(9):R13718798982
PubMedGoogle ScholarCrossref 25.Connor C, Cheung I, Simon A, Jakovcevski M, Weng Z, Akbarian S. A simple method for improving the specificity of anti-methyl histone antibodies.
Epigenetics. 2010;5(5):392-39520458167
PubMedGoogle ScholarCrossref 26.Vargas DL, Nascimbene C, Krishnan C, Zimmerman AW, Pardo CA. Neuroglial activation and neuroinflammation in the brain of patients with autism.
Ann Neurol. 2005;57(1):67-8115546155
PubMedGoogle ScholarCrossref 27.Morgan JT, Chana G, Pardo CA, Achim C, Semendeferi K, Buckwalter J, Courchesne E, Everall IP. Microglial activation and increased microglial density observed in the dorsolateral prefrontal cortex in autism.
Biol Psychiatry. 2010;68(4):368-37620674603
PubMedGoogle ScholarCrossref 28.Barski A, Cuddapah S, Cui K, Roh TY, Schones DE, Wang Z, Wei G, Chepelev I, Zhao K. High-resolution profiling of histone methylations in the human genome.
Cell. 2007;129(4):823-83717512414
PubMedGoogle ScholarCrossref 30.Stadler F, Kolb G, Rubusch L, Baker SP, Jones EG, Akbarian S. Histone methylation at gene promoters is associated with developmental regulation and region-specific expression of ionotropic and metabotropic glutamate receptors in human brain.
J Neurochem. 2005;94(2):324-33615998284
PubMedGoogle ScholarCrossref 31.Huang HS, Matevossian A, Whittle C, Kim SY, Schumacher A, Baker SP, Akbarian S. Prefrontal dysfunction in schizophrenia involves mixed-lineage leukemia 1-regulated histone methylation at GABAergic gene promoters.
J Neurosci. 2007;27(42):11254-1126217942719
PubMedGoogle ScholarCrossref 32.Ernst C, Chen ES, Turecki G. Histone methylation and decreased expression of TrkB.T1 in orbital frontal cortex of suicide completers.
Mol Psychiatry. 2009;14(9):830-83219696771
PubMedGoogle ScholarCrossref 33.Greer PL, Hanayama R, Bloodgood BL, Mardinly AR, Lipton DM, Flavell SW, Kim TK, Griffith EC, Waldon Z, Maehr R, Ploegh HL, Chowdhury S, Worley PF, Steen J, Greenberg ME. The Angelman syndrome protein Ube3A regulates synapse development by ubiquitinating arc.
Cell. 2010;140(5):704-71620211139
PubMedGoogle ScholarCrossref 34.Drew CJ, Kyd RJ, Morton AJ. Complexin 1 knockout mice exhibit marked deficits in social behaviours but appear to be cognitively normal.
Hum Mol Genet. 2007;16(19):2288-230517652102
PubMedGoogle ScholarCrossref 35.Novak G, Gallo A, Zai CC, Meltzer HY, Lieberman JA, Potkin SG, Voineskos AN, Remington G, Kennedy JL, Levesque D, Le Foll B. Association of the orphan nuclear receptor NR4A1 with tardive dyskinesia.
Psychiatr Genet. 2010;20(1):39-4320010315
PubMedGoogle ScholarCrossref 36.Lee T, Marticke S, Sung C, Robinow S, Luo L. Cell-autonomous requirement of the USP/EcR-B ecdysone receptor for mushroom body neuronal remodeling in
Drosophila.
Neuron. 2000;28(3):807-81811163268
PubMedGoogle ScholarCrossref 37.Thakker-Varia S, Jean YY, Parikh P, Sizer CF, Jernstedt Ayer J, Parikh A, Hyde TM, Buyske S, Alder J. The neuropeptide VGF is reduced in human bipolar postmortem brain and contributes to some of the behavioral and molecular effects of lithium.
J Neurosci. 2010;30(28):9368-938020631166
PubMedGoogle Scholar 38.Wan C, Shi Y, Zhao X, Tang W, Zhang M, Ji B, Zhu H, Xu Y, Li H, Feng G, He L. Positive association between ALDH1A2 and schizophrenia in the Chinese population.
Prog Neuropsychopharmacol Biol Psychiatry. 2009;33(8):1491-149519703508
PubMedGoogle ScholarCrossref 39.Fischbach GD, Lord C. The Simons Simplex Collection: a resource for identification of autism genetic risk factors.
Neuron. 2010;68(2):192-19520955926
PubMedGoogle ScholarCrossref 40.Firth HV, Richards SM, Bevan AP, Clayton S, Corpas M, Rajan D, Van Vooren S, Moreau Y, Pettett RM, Carter NP. DECIPHER: Database of Chromosomal Imbalance and Phenotype in Humans Using Ensembl Resources.
Am J Hum Genet. 2009;84(4):524-53319344873
PubMedGoogle ScholarCrossref 42.Battaglia A, Hoyme HE, Dallapiccola B, Zackai E, Hudgins L, McDonald-McGinn D, Bahi-Buisson N, Romano C, Williams CA, Brailey LL, Zuberi SM, Carey JC. Further delineation of deletion 1p36 syndrome in 60 patients: a recognizable phenotype and common cause of developmental delay and mental retardation.
Pediatrics. 2008;121(2):404-41018245432
PubMedGoogle ScholarCrossref 43.Potocki L, Bi W, Treadwell-Deering D, Carvalho CMB, Eifert A, Friedman EM, Glaze D, Krull K, Lee JA, Lewis RA, Mendoza-Londono R, Robbins-Furman P, Shaw C, Shi X, Weissenberger G, Withers M, Yatsenko SA, Zackai EH, Stankiewicz P, Lupski JR. Characterization of Potocki-Lupski syndrome (dup(17)(p11.2p11.2)) and delineation of a dosage-sensitive critical interval that can convey an autism phenotype.
Am J Hum Genet. 2007;80(4):633-64917357070
PubMedGoogle ScholarCrossref 44.Wintle RF, Lionel AC, Hu P, Ginsberg SD, Pinto D, Thiruvahindrapduram B, Wei J, Marshall CR, Pickett J, Cook EH, Scherer SW. A genotype resource for postmortem brain samples from the Autism Tissue Program.
Autism Res. 2011;4(2):89-9721254448
PubMedGoogle ScholarCrossref 45.Zhao Y, Marotta M, Eichler EE, Eng C, Tanaka H. Linkage disequilibrium between two high-frequency deletion polymorphisms: implications for association studies involving the glutathione-S transferase (
GST) genes.
PLoS Genet. 2009;5(5):e100047219424424
PubMedGoogle ScholarCrossref 46.Campeau PM, Ah Mew N, Cartier L, Mackay KL, Shaffer LG, Der Kaloustian VM, Thomas MA. Prenatal diagnosis of monosomy 1p36: a focus on brain abnormalities and a review of the literature.
Am J Med Genet A. 2008;146A(23):3062-306919006213
PubMedGoogle ScholarCrossref 47.Kumar RA, Sudi J, Babatz TD, Brune CW, Oswald D, Yen M, Nowak NJ, Cook EH, Christian SL, Dobyns WB. A de novo 1p34.2 microdeletion identifies the synaptic vesicle gene
RIMS3 as a novel candidate for autism.
J Med Genet. 2010;47(2):81-9019546099
PubMedGoogle ScholarCrossref 48.Nishimura Y, Martin CL, Vazquez-Lopez A, Spence SJ, Alvarez-Retuerto AI, Sigman M, Steindler C, Pellegrini S, Schanen NC, Warren ST, Geschwind DH. Genome-wide expression profiling of lymphoblastoid cell lines distinguishes different forms of autism and reveals shared pathways.
Hum Mol Genet. 2007;16(14):1682-169817519220
PubMedGoogle ScholarCrossref 49.Marinova Z, Ren M, Wendland JR, Leng Y, Liang MH, Yasuda S, Leeds P, Chuang DM. Valproic acid induces functional heat-shock protein 70 via class I histone deacetylase inhibition in cortical neurons: a potential role of Sp1 acetylation.
J Neurochem. 2009;111(4):976-98719765194
PubMedGoogle ScholarCrossref 50.Dong E, Nelson M, Grayson DR, Costa E, Guidotti A. Clozapine and sulpiride but not haloperidol or olanzapine activate brain DNA demethylation.
Proc Natl Acad Sci U S A. 2008;105(36):13614-1361918757738
PubMedGoogle ScholarCrossref 51.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(7351):380-38421614001
PubMedGoogle ScholarCrossref 53.Hammoud SS, Nix DA, Zhang H, Purwar J, Carrell DT, Cairns BR. Distinctive chromatin in human sperm packages genes for embryo development.
Nature. 2009;460(7254):473-47819525931
PubMedGoogle Scholar 55.Uher R. The role of genetic variation in the causation of mental illness: an evolution-informed framework.
Mol Psychiatry. 2009;14(12):1072-108219704409
PubMedGoogle ScholarCrossref 56.Nakamine A, Ouchanov L, Jiménez P, Manghi ER, Esquivel M, Monge S, Fallas M, Burton BK, Szomju B, Elsea SH, Marshall CR, Scherer SW, McInnes LA. Duplication of 17(p11.2p11.2) in a male child with autism and severe language delay.
Am J Med Genet A. 2008;146A(5):636-64317334992
PubMedGoogle ScholarCrossref 57.Castermans D, Vermeesch JR, Fryns JP, Steyaert JG, Van de Ven WJ, Creemers JW, Devriendt K. Identification and characterization of the
TRIP8 and
REEP3 genes on chromosome 10q21.3 as novel candidate genes for autism.
Eur J Hum Genet. 2007;15(4):422-43117290275
PubMedGoogle ScholarCrossref 58.Glessner JT, Wang K, Cai G, Korvatska O, Kim CE, Wood S, Zhang H, Estes A, Brune CW, Bradfield JP, Imielinski M, Frackelton EC, Reichert J, Crawford EL, Munson J, Sleiman PM, Chiavacci R, Annaiah K, Thomas K, Hou C, Glaberson W, Flory J, Otieno F, Garris M, Soorya L, Klei L, Piven J, Meyer KJ, Anagnostou E, Sakurai T, Game RM, Rudd DS, Zurawiecki D, McDougle CJ, Davis LK, Miller J, Posey DJ, Michaels S, Kolevzon A, Silverman JM, Bernier R, Levy SE, Schultz RT, Dawson G, Owley T, McMahon WM, Wassink TH, Sweeney JA, Nurnberger JI, Coon H, Sutcliffe JS, Minshew NJ, Grant SF, Bucan M, Cook EH, Buxbaum JD, Devlin B, Schellenberg GD, Hakonarson H. Autism genome-wide copy number variation reveals ubiquitin and neuronal genes.
Nature. 2009;459(7246):569-57319404257
PubMedGoogle ScholarCrossref 59.Sakurai T, Ramoz N, Reichert JG, Corwin TE, Kryzak L, Smith CJ, Silverman JM, Hollander E, Buxbaum JD. Association analysis of the NrCAM gene in autism and in subsets of families with severe obsessive-compulsive or self-stimulatory behaviors.
Psychiatr Genet. 2006;16(6):251-25717106428
PubMedGoogle ScholarCrossref 60.Bonora E, Lamb JA, Barnby G, Sykes N, Moberly T, Beyer KS, Klauck SM, Poustka F, Bacchelli E, Blasi F, Maestrini E, Battaglia A, Haracopos D, Pedersen L, Isager T, Eriksen G, Viskum B, Sorensen EU, Brondum-Nielsen K, Cotterill R, Engeland H, Jonge M, Kemner C, Steggehuis K, Scherpenisse M, Rutter M, Bolton PF, Parr JR, Poustka A, Bailey AJ, Monaco AP.International Molecular Genetic Study of Autism Consortium. Mutation screening and association analysis of six candidate genes for autism on chromosome 7q.
Eur J Hum Genet. 2005;13(2):198-20715523497
PubMedGoogle ScholarCrossref 61.Petek E, Windpassinger C, Vincent JB, Cheung J, Boright AP, Scherer SW, Kroisel PM, Wagner K. Disruption of a novel gene (
IMMP2L) by a breakpoint in 7q31 associated with Tourette syndrome.
Am J Hum Genet. 2001;68(4):848-85811254443
PubMedGoogle ScholarCrossref 62.Weiss LA, Arking DE, Daly MJ, Chakravarti A.Gene Discovery Project of Johns Hopkins and the Autism Consortium. A genome-wide linkage and association scan reveals novel loci for autism.
Nature. 2009;461(7265):802-80819812673
PubMedGoogle ScholarCrossref 63.Melin M, Carlsson B, Anckarsater H, Rastam M, Betancur C, Isaksson A, Gillberg C, Dahl N. Constitutional downregulation of SEMA5A expression in autism.
Neuropsychobiology. 2006;54(1):64-6917028446
PubMedGoogle ScholarCrossref 64.Scheuerle A, Wilson K.
PARK2 copy number aberrations in two children presenting with autism spectrum disorder: further support of an association and possible evidence for a new microdeletion/microduplication syndrome.
Am J Med Genet B Neuropsychiatr Genet. 2011;156B(4):413-42021360662
PubMedGoogle Scholar