Identification of Epigenetically Altered Genes in Sporadic Amyotrophic Lateral Sclerosis

1Amyotrophic lateral sclerosis (ALS) is a progressive and terminal neurodegenerative disease characterized by the selective degeneration of motor neurons within the motor cortex, brainstem and spinal cord. In the United States, approximately 14 cases of ALS are diagnosed each day and 30,000 people are living with the disease. The average time from disease onset to death is 3 years and no treatment that substantially improves the clinical course of the disease is currently available.

2Proposed pathogenic mechanisms of ALS include oxidative stress, glutamate excitotoxicity, impaired axonal transport, neurotrophic deprivation, neuroinflammation, apoptosis, altered protein turnover, and mitochondrial dysfunction,. Moreover, influences from astrocytes and microglia in the motor neuron microenvironment contribute to pathogenesis. In the last 20 years, a search for genetic factors has identified several genes associated with familial ALS (fALS) and a few with sporadic ALS (sALS) –. Because fALS only accounts for 5–10% of all cases of ALS, the causes leading to the vast majority of ALS (sALS) are poorly understood.

3Environmental exposure to toxins, excessive physical activity, dietary factors, and changes in immunity increase the risk of developing sALS. These factors may drive epigenetic changes, which are well suited to explain disease onset and progression in sALS, as they may be acquired throughout life. Epigenetic modifications, including covalent modifications of DNA and histones as well as RNA editing, dynamically regulate gene expression without altering the genetic code,. These modifications are important in chromosome integrity, cellular differentiation, development, and aging,. Two such modifications, 5-methylcytosine (5 mC) and 5-hydroxymethylcytosine (5 HmC) are associated with repression or activation of gene expression, respectively, in response to environmental and developmental factors linked to age-related diseases.

45mC at CpG (cytosine nucleotide separated by a phosphate from a guanine nucleotide) sites is a reversible mechanism facilitated by DNA (cytosine-5)-methyltransferases (DNMTs). Conversely, the Fe(II) and α-ketoglutarate (α-KG)-dependent ten-eleven translocation (TET) family of proteins catalyze oxidation and decarboxylation reactions of 5mC leading to 5-hydroxymethylcytosine (5 HmC), 5-formylcytosine (5 fC) and 5-carboxylcytosine (5 caC),. 5 HmC may be an intermediate for passive (during DNA replication) and active demethylation and/or serve as a docking site for proteins with high affinity for 5 HmC, thereby dissociating interactions between the transcriptional repression machinery and 5 mC.

5In addition to the identification of alterations in global 5 HmC associated with sALS, this study represents one of the first methylation assessments of sALS by integrating methylome and transcriptome profiles of postmortem frozen human spinal cord samples. We identified differentially methylated sALS spinal cord genes exhibiting concordant mRNA expression overrepresented in functional categories implicated in sALS. These data support a role for epigenetic regulation in sALS and it may provide a better understanding of disease pathogenesis and facilitate the discovery of new therapeutic targets.

6A workflow of our data analysis is provided in Fig. 1.

7Chestnut et al. recently reported an increase in DNMTs and 5mC immunoreactivity in ALS brain and spinal cord, suggesting that a global increase in 5mC is associated with the pathogenesis of ALS. We assessed global 5mC of genomic DNA extracted from postmortem human spinal cord samples (sALS, n = 11; matching controls, n = 11; Tables 1, S1) using a colorimetric ELISA approach. We observed a modest but significant 1.4-fold increase in global 5mC in sALS (3.58±0.18) compared to controls (2.56±0.18) (p = 0.0006, Fig. 2), confirming previous observations.

8To identify loci-specific differentially methylated genes (DMGs) in sALS, postmortem spinal cord tissue (Tables 1, S1) was subjected to an unbiased genome-wide methylation screen using the Infinium Human Methylation27 DNA BeadChip (HM27K) array. We identified 4,261 significant differentially methylated autosomal CpG sites, representing 3,574 genes. Functional enrichment analysis identified biological categories of extracellular region (p = 3.5E-30), defense response (p = 2.3E-18), cytokine activity (p = 2.7E-16), immune response (p = 6.9E-15), JAK-STAT signaling pathway (p = 2.2E-7), steroid hormone biosynthesis (p = 3.6E-7), and drug metabolism (p = 3.1E-3) were over-represented in spinal cord DMGs (Fig. 3). This suggests a role for epigenetic regulation in molecular mechanisms known to drive sALS pathogenesis –. Furthermore, the identified DMGs included previously reported epigenetically regulated genes such as runt-related transcription factor 3 (RUNX3), TNF-related apoptosis-inducing ligand (TRAIL/TNFSF10), H19, neuritin 1 (NRN1) and signal transducer and activator of transcription 5A (STAT5A) –, validating our approach. The ALS-dependent differential methylation of three genes encoding the chemokines CKLF-like MARVEL transmembrane domain-containing proteins 2 and 3 (CMTM2 and CMTM3) and the chemokine (C-X-C motif) ligand 12 (CXCL12), as well as the two transcription factors, STAT5A and the CCAAT/enhancer binding protein beta (C/EBPB) was assessed by pyrosequencing (Fig. S1). These results indicate parallel direction of differential methylation between the two assays.

9An increase in promoter 5mC (hyper-methylation) is associated with gene silencing, while a decrease (hypo-methylation) reflects up-regulation of gene expression. To identify functionally relevant epigenes, we performed genome-wide gene expression profiling of total RNA from sALS (n = 12) and control (n = 10) spinal cord (Tables 1, S1). We found 1,182 DEGs in sALS enriched with immune response-related biological categories such as defense response (p = 2.3E-20), inflammatory response (p = 1.3E-18), T-cell activation (p = 3.0E-9), leukocyte activation (p = 1.6E-8), cytokine binding (p = 1.7E-7) as well as carbohydrate binding (p = 7.2E-6), transmembrane receptor protein tyrosine kinase signaling (p = 2.6E-5), cell death (p = 1.5E-3), and JAK-STAT signaling (p = 3.6E-2) (Fig. 3). Several of these biological categories overlap with the 5mC profile, suggesting methylation may play a role in the regulation of these DEGs.

10DMGs and DEGs were analyzed for their direction of change to select genes demonstrating concordant and potentially true regulation of gene expression by methylation. Between 3,574 DMGs and 1,182 DEGs, 251 genes were common (Table S2), of which approximately 70% were hyper-methylated as observed with global 5mC. Among these common genes, 112 epigenes demonstrated concordant direction of methylation and expression changes (hypo-methylation/up-regulation (51 genes; Table 2) or hyper-methylation/down-regulation (61 genes; Table 3) (Fig. 4). These concordant epigenes were enriched in biological categories including immune response (p = 7.1E-5), plasma membrane part (p = 9.6E-5), defense response (p = 0.03), and neuron adhesion (p = 0.05) (Fig. 3).

11Gene co-citation network analysis revealed potential associations between some of the concordant genes identified in our screens. Among the 112 concordant epigenes, 53 genes were co-cited at least once in PubMed abstracts at the sentence or abstract-level. Within the co-citation network, the epigenes were classified into representative biological functions such as immune response, antigen presentation, tumor/suppressor related, and extracellular matrix repair (Fig. 5). In this network, the transcription factors C/EBPB and STAT5A were highly connected, suggesting a role in sALS pathogenesis. Moreover, some of these connections were directly or indirectly associated with concordant genes previously related to neuronal development such as muscle segment homeobox (msh) 2 (MSX2), megalencephalic leukoencephalopathy with subcortical cysts 1 (MLC1), RUNX3, catenin delta-2 (CTNND2), receptor tyrosine kinase (AXL), neuronatin (NNAT), and NRN1 – (Fig. 5).

12SciMiner identified 4,128 genes from ALS-related publications (as of 7/23/2012), which were compared to our 112 concordant epigenes. Fourteen genes were identified in two or more ALS-related publications with frequencies that were significantly different from those in over 20 million abstracts in PubMed (p<0.05). Fifty-one genes demonstrated ≥2-fold altered expression, including the chitinase 3-like protein 2 (CHI3L2), the triggering receptor expressed on myeloid cells-2 (TREM2), cathepsin Z (CTSZ), the lumican precursor protein (LUM), H19, and TRAIL/TNFSF10 (Tables 2, 3). Thus, bioinformatics evaluation of the concordant epigenes identified by integrating methylomic and transcriptomic analyses detected both novel and previously known ALS-related genes.

13Expression of 14 concordant epigenes selected either from the ALS-related literature or from the expression array data was confirmed by RT-PCR (Figs. 6 and S2, Table 4). Of the genes previously related to the ALS literature, NRN1, FMO1, and the lumican precursor protein (LUM) were under-expressed, while the lysosomal protease CTSZ was over-expressed in sALS. No significant difference between sALS and control subjects was observed for FES-upstream region (FURIN). Novel sALS-associated epigenes such as STAT5A, TREM2, the high-affinity IgE receptor (FCER1G), CHI3L2, and the proton-couple divalent metal ion transporter solute carrier family 11 member 1A (SLC11A1) were over-expressed in sALS. SLC11A1 presented the highest increase by 12.7-fold. Down-regulation was validated for gap junction ß-2 (GJB2)/Connexin-26 as well as imprinted genes such as H19, NNAT, and the paternally expressed 10 (PEG10). In summary, the RT-PCR expression data indicate high concordance with the microarray expression data, validating our results.

145 HmC, an alternate epigenetic modification of DNA, is increased in brain compared to other human tissues, and alterations in global 5 HmC are associated with age-related neurodegenerative disorders, suggesting an important role of 5 HmC in neuronal the function –. We measured global 5 HmC for sALS and control spinal cord samples previously analyzed for global 5mC. We observed an approximately 3.0-fold increase in global 5 HmC in sALS (0.31±0.02) compared to controls (0.11±0.03) (p<0.0001) (Fig. 7). This is the first report of aberrant levels of global 5 HmC associated with the pathogenesis of sALS.

15High correlation of epigenetic marks in spinal cord and blood may be useful for diagnostic and therapeutic application in ALS. We investigated whether global 5mC and 5 HmC would be altered in sALS whole blood similarly to spinal cord. Whole blood genomic DNA from a different cohort (Tables 5, S1) was subjected to global 5mC and global 5 HmC by ELISA. The levels of percent global 5mC and 5 HmC in whole blood were 10 fold lower compared to spinal cord, in agreement with recent reports,. Contrary to spinal cord, no differential percent global 5mC (controls, 0.408±0.300; ALS, 0.405±0.027 p = 0.941) and global 5 HmC (controls, 0.033±0.005; ALS, 0.032±0.004; p = 0.401) (Fig. 8) were observed in whole blood.

16Although several genes have been implicated in the pathogenesis of ALS, the causes leading to most cases remain unknown. Environmental factors may be associated with the onset and development of sALS by altering epigenetic regulation,. The aim of this study was to identify sALS-associated epigenetic marks resulting in aberrant gene expression. Abnormal 5mC patterns of repetitive elements such as Alu and LINE1, as well as altered function of methylation regulators such as the DNMTs, lead to changes in global 5mC or 5 HmC associated with neurodegeneration,. We demonstrate increased global methylation in sALS spinal cord, perhaps due to an increase in DNMT activity. Furthermore, we report for the first time an increase in global 5 HmC in sALS spinal cord. Increased 5mC and 5 HmC may be due to 5mC providing more substrate for the TET proteins, TETs are not differentially expressed in spinal cord sALS according to our microarray data (data not shown). TET should decrease the amount of 5mC only if 5mC is not increasing at a faster rate than the oxidation reaction. Although normal aging leads to increased global 5 HmC in mouse hippocampal DNA independently of increased levels of oxidative stress, in ALS, increased oxidative DNA damage and free radicals may contribute to global 5 HmC dysregulation. The base excision repair (BER) pathway responsible for oxidative DNA damage restoration and one of the active demethylation pathways, is deficient in ALS,,.

17Methylomics and transcriptomics analyses identified potential biologically relevant epigenes in postmortem sALS spinal cord. These epigenes were enriched with biological functions related to inflammation and the immune responses, previously linked to ALS –. Our data suggest that alterations in gene expression of immune-related genes in sALS may be regulated by methylation. Immune-related concordant epigenes including TREM2, chemokine (C-C motif) receptor 1/RANTES receptor (CCR1), SLC11A1, the transmembrane receptor C-type lectin domain family 4 member A isoform 1 (CLEC4A), and the IgE receptor (FCER1G) were found to be over-expressed in sALS. Our findings suggest an infiltration of myeloid cells, mast cells, or natural killer cells to the damaged area and/or activation of resident microglia –. Supporting our observations, neuro-inflammation was recently associated with systemic macrophage activation independent of T-cell activation and the recruitment of activated inflammatory monocytes to the spinal cord in ALS,. Although immunosuppressive and anti-inflammatory therapies have shown to delay disease onset in ALS animal models, clinical trials have not revealed a major effect on disease progression or survival –. This suggests that continuous activation of microglia leading to neuronal damage surpasses the capacity of the nervous system to respond to immunosuppressive and anti-inflammatory therapies at later stages of ALS, implicating a need for biomarkers identifying early immune-related changes in sALS.

18Co-citation network and literature mining approaches identified connections among novel and previously implicated ALS-related epigenes and pathways,. The transcription factors STAT5A and C/EBPB are highly connected in our co-citation network and their interplay promote activation of various genes including interleukin-6 (IL-6). Moreover, recent reports implicate C/EBPB and STAT5A in ALS pathogenesis and neurodegeneration. For instance, expression of C/EBPB in ALS microglia from spinal cord suggests an important role of C/EBPB in the regulation of neurotoxic genes in the ALS neuronal microenvironment,. Furthermore, changes in STAT5A expression may reflect an altered inflammatory response contributing to the pathogenesis of ALS. Over-expression of STAT5A reduces neuronal degeneration associated with spinal muscular atrophy, a neurodegenerative disease with similar pathogenesis as ALS, and it provides oligodendrocyte protection, which in turn favors neuronal environment preservation –. Whether positive regulation of STAT5A in sALS is due to an anti-apoptotic response to compensate for the degeneration of the nervous system, or its over-expression is responsible, in part, for the pathogenesis of the disease remains to be determined. Interestingly, we observed potential transcription factor binding sites (TFBSs) for STAT5A and C/EBPB in 40% and 48% of the promoters of our identified DEGs, respectively; the binding sites for STAT5A and C/EBPB are 1.2 (p = 4.1E-12) and 1.3 (p = 3.8E-13) times more frequent in the DEGs than in the vertebrate promoters, respectively. Our observations suggest epigenetic mechanisms, in part, drive the expression of central regulators of downstream targets in sALS.

19Our study identified ALS-dependent methylation dysregulation of several genes previously implicated in neuronal development, differentiation, and proliferation, including Slit-Robo Rho GTPase activating protein 1 (SRGAP1), crumbs homolog 1 (CRB1), MSX2, MLC1, CTNND2, AXL, RUNX3, NNAT, and NRN1 –,,. Altered expression and/or mutations in CRB1, MLC1, and CTNND2 are associated with mental retardation and neurodegeneration,,,. Moreover, we found CRB1 to be hypo-methylated and up-regulated in sALS, in agreement with overexpression observed by Kudo et al. in spinal cord from a mouse model of familial ALS (SOD1 G93A) as well as human sALS. Interestingly, most of these genes were identified by our literature-based association network of concordant epigenes and were connected to C/EBPB and STAT5A. Analysis of the promoter region of these genes indicates a high incidence of potential TFBSs for these two transcription factors, suggesting a potential role of STAT5A and C/EBPB in the regulation of neuronal genes in the pathogenesis of sALS. Our observations suggest sALS-related alterations in methylation may lead to aberrant expression of genes required for neuronal homeostasis. Nevertheless, more studies need to be done to address the role of methylation, STAT5A, and C/EBPB in the regulation of neuronal genes.

20Another sALS-related epigene, CTSZ warrants further investigation since its aberrant expression is associated with neurodegeneration by promoting neurotoxin elimination in the damaged cellular environment. Expression of CTSZ as well as two other members of the cathepsin family, cathepsins B and D, increases in human and rodent ALS spinal cord and mutant SOD1 (G86R, G93A) mouse skeletal muscle suggesting they play an important role in ALS,.

21Except for optineurin (OPTN), which was identified as a hypo-methylated DMG without demonstrating changes in gene expression, loci known to be mutated in fALS were not present in our concordant epigenes. This agrees with recent studies indicating that promoter regions of SOD1, VEGF, and metallothioneins I and II are not differentially methylated in sALS,. When compared with the ALS Online genetics Database (ALSoD)-reported genes and other ALS-dependent methylation/gene expression of profiling studies,, –, we observed a modest overlap of four concordant epigenes; Purkinje cell protein 4 (PCP4), catenin (CTNNAL1), fibroblast growth factor 18 (FGF18), and flavin containing monooxygenase 1 (FMO1). Furthermore, five of our concordant genes presented opposite direction of expression when compared to known ALS-dependent differentially expressed gene. Our data indicate that epigenetic mechanisms are potential regulators of these key genes in ALS.

22Based on the large number of genes identified in the methylation (3,574 genes) and expression (1,182 genes) arrays, relatively few sALS-associated genes presented concordant direction between methylation and gene expression. The low occurrence of a small subset of genes potentially regulated by CpG modification in such a way that hyper-methylation promotes gene silencing and hypo-methylation promotes gene expression has been previously documented. 5mC within promoter regions is associated with repression of gene expression by interfering with transcription factor binding or by providing a binding site for transcriptional repressors. Interestingly, over half (55%) of the 251 common DMGs/DEGs presented same direction of 5mC and expression. In some cases, 5mC positively regulates gene transcription by promoting transcription factor binding at promoter regions or, more commonly, by modifying intragenic CpG sites facilitating transcription efficiency, histone conformation, and regulating levels of sense and antisense mRNA. Furthermore, 5 HmC, a highly enriched modification in brain, correlates with increased gene expression. HM27K does not differentiate between 5mC and 5 HmC; therefore, some of the common epigenes presenting same direction of methylation and expression may be regulated by 5 HmC.

23Although the high incidence of same direction sALS concordant epigenes parallels the high levels of global 5 HmC in spinal cord, loci specific 5 HmC modifications associated with sALS remain to be identified. Gene expression of non-common DEG (non-DMGs) could be determined, in part, by 5mC-dependent regulation of transcription factors. In addition to STAT5 and C/EBPB, we identified several transcription factors as concordant genes such as the transcription factor 7 (TCF7), RUNX3, IKAROS family zinc finger 1 (IKZF1), MSX2, and hypoxia inducible factor 3, alpha subunit (HIF3A). Furthermore, regulation of gene expression is a dynamic and complex mechanism and the interplay of several epigenetic pathways has been reported to modulate adult neurogenesis. Therefore, alterations to epigenetic networks in conjunction with genetic predisposition may result in the development of sALS.

24The prospect of identifying sALS epigenetic biomarkers in blood is exciting as it provides a minimally invasive alternative for sALS diagnostic and prognostic assessments. Although we did not detect significant global 5mC and 5 HmC differences in blood and inflammation-related epigene biomarkers may reflect systemic inflammatory changes rather than neuronal changes, further investigation of individual loci may provide potential epigenetic biomarkers for sALS.

25There were several limitations to our study. First, a relatively small number of samples were analyzed and loci-specific 5 HmC analysis is still needed. Nevertheless, this is an initial step towards identifying epigenetic mechanisms altering key pathways leading to sALS, which will be validated in larger cohorts. Second, sALS postmortem tissue reflects the terminal disease stage rather than the pathogenic mechanisms leading to disease onset and progression. As sALS-affected motor neurons deteriorate at the terminal stage and heterogeneous tissue consisting of both gray and white matter was analyzed, our results may represent epigenetic regulation of the neuronal microenvironment, including microglia activation and the scarce neurons surviving the degenerative process,. This may explain, in part the discrepancy in the direction of expression of common and concordant genes reported here with other sALS genome-wide expression profiles, as well as the heavily represented inflammation-related genes, in our concordant epigenes, which are not differentially expressed specifically in sALS motor neurons or ventral horns. Finally, more studies are needed to concretely identify whether or not the genes identified in this study are involved in ALS pathogenesis.

26Advances in identifying epigenetic regulators in disease states have led to new therapeutic approaches. Interestingly, demethylating agents have been extensively studied to reverse aberrant epigenetic changes associated with cancer and more recently, histone deacetylase inhibitors have shown to have neuroprotective properties in animal models of neurodegenerative diseases. These observations suggest reversible epigenetic modifications carry the potential for therapeutic treatment in sALS. We contend that environmental life exposures result in failure to maintain epigenetic homeostasis in the nervous system microenvironment leading to global and loci specific aberrant regulation of gene expression in sALS-affected tissue. Ascertaining the role of epigenetic regulation may provide a better understanding of the pathogenesis of sALS and new therapeutic targets.

27Frozen human spinal cord samples from 12 Caucasian sALS subjects and 11 age and gender-matched neurologically-normal controls were obtained from the National Center for Child Health and Human Development (NICHD) Brain and Tissue Bank for Developmental Disorders at the University of Maryland, Baltimore, MD (Table 1). Whole blood was collected in EDTA tubes from a different cohort of 11 Caucasian sALS and 11 age- (5±years) and gender-matched neurologically-normal control subjects at the University of Michigan ALS Consortium (Table 4). Table S1 summarizes the samples used for each assay.

28The participants donating blood reviewed and signed a written informed consent under a protocol reviewed and approved by the University of Michigan Medical School Institutional Review Board (HUM00028826).

29Genomic DNA was extracted from 50 mg of frozen postmortem spinal cord tissue (mostly grey matter including the ventral horn, but some white matter included) using the Promega Maxwell 16 Tissue DNA Purification kit and a Maxwell instrument (Promega Co, Madison, WI). Genomic DNA (1 µg) was bisulfite converted with an EZ DNA Methylation Kit (Zymo Research, Irvine, CA) according to the manufacturer’s instructions. Total RNA was extracted from the same tissue for methylation profiling (Table S1) using the RNeasy kit and treated with RNAse-free DNase1 according to the manufacturer’s instructions (Qiagen, Valencia, CA). Automated genomic DNA extraction from whole blood was performed at the Michigan Institute for Clinical & Health Research (MICHR) at the University of Michigan using Autogen FlexStar (Autogen, Holliston, MA) and Qiagen Flexigene reagents. Nucleotide concentration was assessed using a Nanodrop 2000 (Thermo Scientific) and RNA integrity was determined by microfluid electrophoresis with a 21000 Bioanalyzer (Agilent Technologies, Palo Alto, CA).

30Differences in genomic DNA global methylation (Global 5mC) and hydroxymethylation (Global 5 HmC) from sALS and control spinal cord or whole blood were determined in duplicate using the colorimetric enzyme-linked immuno-sorbent assay (ELISA) MethylFlash (Methylated or Hydroxymethylated) DNA Quantification Kits according to the manufacturer’s directions (Epigentek Group Inc., Farmingdale, NY). The absorbance at 450 nm was captured in a Fluoroskan Ascent microplate reader (Labsystems, Vienna, Va). The percentage of Global 5mC and 5 HmC is expressed as mean ± standard error mean (SEM). The two-tailed t-test was used for statistical comparison. Graphs and statistical analysis were obtained with GraphPad Prism 5.

31For high-throughput methylation profiling, 200 ng of bisulfite-converted DNA was whole-genome amplified (WGA), enzymatically fragmented, purified, and hybridized to the Infinium Human Methylation27 DNA BeadChip array (HM27K; Illumina, Inc., San Diego, CA) following the manufacturer’s instructions at the University of Michigan Sequencing Core. The HM27K quantitatively determines DNA methylation for 27,578 CpG sites spanning 14,495 genes. DMGs were identified using Illumina’s GenomeStudio software. Single-base resolution corresponding to DNA methylation levels for each locus was reported and the methylation level is given by a beta (ß) value describing the percentage of the degree of methylation ranging from 0 (no methylation) to 1 (complete methylation). Any methylation value with a detection P-value >0.05 was excluded. Differential methylation of the selected CpG target regions of autosomal chromosomes between sALS and control groups were tested using Illumina Custom algorithm with multiple testing corrections applied. DiffScore, GenomeStudio’s statistical significance score for differential methylation, of >13 for hyper-methylation or <-13 for hypo-methylation, equivalent to False Discovery Rate (FDR) <5%, were used.

32Microarray gene expression analysis was performed as previously described in our published protocols. Briefly, RNA samples with an RNA integrity number (RIN) >6.4 were used for further microarray and real-time PCR analysis. Total RNA (75 ng) was amplified and biotin-labeled using the Ovation Biotin-RNA Amplification and Labeling System (NuGEN Technologies, Inc., San Carlos, CA) according to the manufacturer’s instructions. Amplification and hybridization was performed at the University of Michigan DNA Sequencing Core Affymetrix and Microarray Core Group (Ann Arbor, MI) using the Affymetrix GeneChip Human Genome U133 Plus 2.0 Array measuring over 47,000 transcripts representing over 20,000 human genes.

33Affymetrix CEL files were analyzed using a local version of the GenePattern genomic analysis platform from the Broad Institute. Samples were Robust Multi-array Average (RMA) normalized using the BrainArray Custom CDF HGU133Plus2_Hs_ENTREZG version 14. Microarray quality was assessed as previously published. Briefly, probe-level modeling (PLM) and quality metrics provided by the BioConductor affy package were used to identify low-quality arrays –. Outlier arrays, skewed away from other arrays, identified by Principal Component Analysis (PCA) were excluded from further analyses. Intensity-Based Moderated T-statistic (IBMT) was employed to identify DMGs with a 10% FDR cut-off between sALS and control samples.

34Concordant epigenes are those exhibiting significant differential methylation (hyper- or hypo-methylation) and a parallel change of gene expression (under- or over-expression, respectively) between sALS and control. Differentially methylated (DMGs)/expressed (DEGs) genes were subjected to bioinformatics analyses.

35Database for Annotation, Visualization and Integrated Discovery (DAVID; http://david.abcc.ncifcrf.gov/), was used to identify enriched molecular biological functions and ALS-relevant pathways of concordant epigenes. A Benjamini-Hochberg corrected P-value of 0.05 was used as the cut-off for statistically significant over-representation.

36A literature mining approach was used to obtain a comprehensive list of potential ALS-associated targets (genes/proteins). SciMiner, a web-based literature mining tool,, retrieves, processes documents, and identifies potential ALS-associated targets from the ALS-related literature, defined by a PubMed-style query of “Amyotrophic Lateral Sclerosis”. The concordant epigenes were compared against the literature-derived ALS-associated targets that were observed in at least 2 or more papers and whose frequency (in terms of the number of papers) was significantly different from the background. Fisher’s exact test (p-value <0.05) was used to determine whether each gene’s frequency was significantly different from the complete collection of abstracts of over 20 million papers in PubMed. The concordant genes identified by the high-throughput arrays were compared with these literature-derived ALS-related genes to identify which known disease-relevant genes are most highly methylated/expressed and, consequently, likely involved in disease pathogenesis. The resulting genes were designated as literature-derived ALS-associated epigenes.

37To elucidate the functional relationships among the concordant epigenes, we generated transcriptional networks using Genomatix Pathway Systems (GePS; Genomatix Software GmbH, Munich, Germany) with a sentence-level co-citation filter. Two genes co-cited at the sentence level in the literature are linked, resulting in a co-citation network. Additionally, transcriptional regulatory information of predicted transcription factor binding sites (TFBS) in promoter regions of genes could be further incorporated. The network allows the visualization of concordant epigenes, their potential associations, and transcriptional regulation with each other. Therefore, it helps in the identification of key genes that are highly connected to genes, and which play potentially important roles in the pathogenesis of sALS. Potential TFBSs of two highly connected genes in the network, STAT5A and C/EBPB, were searched among the promoters of the concordant epigenes using MatInspector (Genomatix).

38To validate the HM27K arrays, we assessed gene-specific methylation of three selected cytokine genes based on the fact that immune response is associated with the pathogenesis of sALS and two transcription factors. Amplicons of the promoter regions of the genes coding for the CKLF-like MARVEL transmembrane domain-containing proteins 2 and 3 (CMTM2 and CMTM3), the chemokine (C-X-C motif) ligand 12 (CXCL12), signal transducer and activator of transcription 5A (STAT5A), and CCAAT/enhancer binding protein beta (C/EBPB) were generated in 30 ml reactions using the PyroMark kit (Qiagen, Valencia, CA) with 4.8 pmol of the forward non-biotinylated primer, 2.4 pmol of the reverse biotinylated primer (Table S3), and 25 ng of bisulfite converted genomic DNA as previously described. PCR conditions: 95°C for 15 min, 50 cycles [95°C for 30 s, 40–50°C for 30 s, 72°C for 20 s], 72°C for 10 min. Ten ml of the amplicon was Streptavidin Sepharose (Amersham Bioscience, Uppsala, Sweden) were purified, denatured with 0.2 M NaOH, and pyrosequenced using 0.5 mM of sequencing primer in a PSQ96 HS System (Qiagen) following the manufacturer’s protocol. Percent methylation of the region analyzed containing the identified Illumina methylation site or individual sites are presented as mean ± SEM with a two-sample equal variance t-test using GraphPad Prism 5.

39cDNA was generated by reverse transcription from total RNA isolated for microarray analysis using an iScript cDNA Synthesis Kit (Bio-Rad, Hercules, CA). RT-PCR was performed in triplicate using sequence-specific primers (Table S4) with SYBR Green PCR reagents (Bio-Rad, Hercules, CA). The PCR amplification profile was as follows: 95°C for 5 min, [denaturation at 95°C for 30 s, annealing at 55–60°C for 60 s, and extension at 72°C for 30 s] x40 cycles, and a final phase of 72°C for 5 min. The fluorescence threshold CT value, representing mRNA expressed in sALS samples, was calculated by the iCycler iQ system software. mRNA levels were normalized to an endogenous reference (ΔCT) and then relative to the control group (ΔΔCT). Levels of PCR products are demonstrated as mean ± SEM and a two-sample equal variance t-test was performed using GraphPad Prism 5 to confirm that mRNA levels were significantly different between sALS and control.