Novel methods for integration and visualization of genomics and genetics data in Alzheimer's disease.

INTRODUCTION: Numerous omics studies have been conducted to understand the molecular networks involved in Alzheimer's disease (AD), but the pathophysiology is still not completely understood; new approaches that enable neuroscientists to better interpret the results of omics analysis are required. METHODS: We have developed advanced methods to analyze and visualize publicly-available genomics and genetics data. The tools include a composite clinical-neuropathological score for defining AD, gene expression maps in the brain, and networks integrating omics data to understand the impact of polymorphisms on AD pathways. RESULTS: We have analyzed over 50 public human gene expression data sets, spanning 19 different brain regions and encompassing three separate cohorts. We integrated genome-wide association studies with expression data to identify important genes in the pathophysiology of AD, which provides further insight into the calcium signaling and calcineurin pathways. DISCUSSION: Biologists can use these freely-available tools to obtain a comprehensive, information-rich view of the pathways in AD.

24-hour rhythms of DNA methylation and their relation with rhythms of RNA expression in the human dorsolateral prefrontal cortex.

Circadian rhythms modulate the biology of many human tissues, including brain tissues, and are driven by a near 24-hour transcriptional feedback loop. These rhythms are paralleled by 24-hour rhythms of large portions of the transcriptome. The role of dynamic DNA methylation in influencing these rhythms is uncertain. While recent work in Neurospora suggests that dynamic site-specific circadian rhythms of DNA methylation may play a role in modulating the fungal molecular clock, such rhythms and their relationship to RNA expression have not, to our knowledge, been elucidated in mammalian tissues, including human brain tissues. We hypothesized that 24-hour rhythms of DNA methylation exist in the human brain, and play a role in driving 24-hour rhythms of RNA expression. We analyzed DNA methylation levels in post-mortem human dorsolateral prefrontal cortex samples from 738 subjects. We assessed for 24-hour rhythmicity of 420,132 DNA methylation sites throughout the genome by considering methylation levels as a function of clock time of death and parameterizing these data using cosine functions. We determined global statistical significance by permutation. We then related rhythms of DNA methylation with rhythms of RNA expression determined by RNA sequencing. We found evidence of significant 24-hour rhythmicity of DNA methylation. Regions near transcription start sites were enriched for high-amplitude rhythmic DNA methylation sites, which were in turn time locked to 24-hour rhythms of RNA expression of nearby genes, with the nadir of methylation preceding peak transcript expression by 1-3 hours. Weak ante-mortem rest-activity rhythms were associated with lower amplitude DNA methylation rhythms as were older age and the presence of Alzheimer's disease. These findings support the hypothesis that 24-hour rhythms of DNA methylation, particularly near transcription start sites, may play a role in driving 24-hour rhythms of gene expression in the human dorsolateral prefrontal cortex, and may be affected by age and Alzheimer's disease.

A 17q12 allele is associated with altered NK cell subsets and function.

NK cells play an important role in innate immunity. A previous genome-wide association study demonstrated an association between a 17q12 allele (rs9916629(C)) and lower frequency of CD3(-)CD56(+) NK cells in peripheral blood. We performed an analysis that not only replicates the original result of the genome-wide association study (p = 0.036) but also defines the specific cell subpopulations and functions that are modulated by the rs9916629 polymorphism in a cohort of 96 healthy adult subjects using targeted multiparameter flow cytometric profiling of NK cell phenotypes and functions. We found that rs9916629(C) is associated with alterations in specific NK cell subsets, including lower frequency of predominantly cytotoxic CD56(dim) NK cells (p = 0.011), higher frequency of predominantly regulatory CD56(bright) NK cells (p = 0.019), and a higher proportion of NK cells expressing the inhibitory NKG2A receptor (p = 0.0002). Functionally, rs9916629(C) is associated with decreased secretion of macrophage inflammatory protein-1ß by NK cells in the context of Ab-dependent cell-mediated cytotoxicity (p = 0.039) and increased degranulation in response to MHC class I-deficient B cells (p = 0.017). Transcriptional profiling of NK cells suggests that rs9916629 influences the expression of transcription factors such as TBX21, which has a role in NK cell differentiation, offering a possible mechanism for the phenotypic and functional differences between the different alleles. The rs9916629(C) allele therefore has a validated effect on the proportion of NK cells in peripheral blood and skews NK cells toward a specific phenotypic and functional profile, potentially influencing the impact that these innate immune cells have on infection and autoimmunity.

Broad proteomics analysis of seeding-induced aggregation of α-synuclein in M83 neurons reveals remodeling of proteostasis mechanisms that might contribute to Parkinson's disease pathogenesis.

Aggregation of misfolded α-synuclein (α-syn) is a key characteristic feature of Parkinson's disease (PD) and related synucleinopathies. The nature of these aggregates and their contribution to cellular dysfunction is still not clearly elucidated. We employed mass spectrometry-based total and phospho-proteomics to characterize the underlying molecular and biological changes due to α-syn aggregation using the M83 mouse primary neuronal model of PD. We identified gross changes in the proteome that coincided with the formation of large Lewy body-like α-syn aggregates in these neurons. We used protein-protein interaction (PPI)-based network analysis to identify key protein clusters modulating specific biological pathways that may be dysregulated and identified several mechanisms that regulate protein homeostasis (proteostasis). The observed changes in the proteome may include both homeostatic compensation and dysregulation due to α-syn aggregation and a greater understanding of both processes and their role in α-syn-related proteostasis may lead to improved therapeutic options for patients with PD and related disorders.

A link between SIN1 (MAPKAP1) and poly(rC) binding protein 2 (PCBP2) in counteracting environmental stress.

When SIN1 (MAPKAP1) was used as the bait in a two-hybrid screen of a human bone marrow cDNA library, its most frequent partner was poly(rC) binding protein 2 (PCBP2/hnRNP-E2), which associates with the N-terminal domain of SIN1 and can be coimmunoprecipitated with SIN1 and the cytoplasmic domain of the IFN receptor IFNAR2 from HeLa cells. SIN1, but not PCBP2, also associates with the receptors that bind TNFalpha. PCBP2 is known to bind pyrimidine-rich repeats within the 3' UTR of mRNAs and has been implicated in control of RNA stability and translation and selective cap-independent transcription. RNAi silencing of either SIN1 or PCBP2 renders cells sensitive to basal and stress-induced apoptosis. Stress in the form of TNFalpha and H(2)O(2) treatments rapidly raises the cell content of SIN1 and PCBP2, an effect reversible by inhibiting MAPK14. A meta analysis of human microarray information with an algorithm that discerns similarities in gene-regulatory profiles shows that SIN1 and PCBP2 are generally coregulated with large numbers of genes implicated in both cell survival and death and in cellular stress responses, including RNA translation and processing. We predict that SIN1 is a scaffold protein that organizes antiapoptotic responses in stressed cells, whereas PCBP2, its binding partner, provides for the selective expression of cell survival factors through posttranslational events.

PRIMEGENS-v2: genome-wide primer design for analyzing DNA methylation patterns of CpG islands.

MOTIVATION: DNA methylation plays important roles in biological processes and human diseases, especially cancers. High-throughput bisulfite genomic sequencing based on new generation of sequencers, such as the 454-sequencing system provides an efficient method for analyzing DNA methylation patterns. The successful implementation of this approach depends on the use of primer design software capable of performing genome-wide scan for optimal primers from in silico bisulfite-treated genome sequences. We have developed a method, which fulfills this requirement and conduct primer design for sequences including regions of given promoter CpG islands. RESULTS: The developed method has been implemented using the C and JAVA programming languages. The primer design results were tested in the PCR experiments of 96 selected human DNA sequences containing CpG islands in the promoter regions. The results indicate that this method is efficient and reliable for designing sequence-specific primers. AVAILABILITY: The sequence-specific primer design for DNA meth-ylated sequences including CpG islands has been integrated into the second version of PRIMEGENS as one of the primer design features. The software is freely available for academic use at http://digbio.missouri.edu/primegens/.

The Role of MAPT Haplotype H2 and Isoform 1N/4R in Parkinsonism of Older Adults.

BACKGROUND AND OBJECTIVE: Recently, we have shown that the Parkinson's disease (PD) susceptibility locus MAPT (microtubule associated protein tau) is associated with parkinsonism in older adults without a clinical diagnosis of PD. In this study, we investigated the relationship between parkinsonian signs and MAPT transcripts by assessing the effect of MAPT haplotypes on alternative splicing and expression levels of the most common isoforms in two prospective clinicopathologic studies of aging. MATERIALS AND METHODS: using regression analysis, controlling for age, sex, study and neuropathology, we evaluated 976 subjects with clinical, genotyping and brain pathology data for haplotype analysis. For transcript analysis, we obtained MAPT gene and isoform-level expression from the dorsolateral prefrontal cortex for 505 of these subjects. RESULTS: The MAPT H2 haplotype was associated with lower total MAPT expression (p = 1.2x10-14) and global parkinsonism at both study entry (p = 0.001) and proximate to death (p = 0.050). Specifically, haplotype H2 was primarily associated with bradykinesia in both assessments (p<0.001 and p = 0.008). MAPT total expression was associated with age and decreases linearly with advancing age (p<0.001). Analysing MAPT alternative splicing, the expression of 1N/4R isoform was inversely associated with global parkinsonism (p = 0.008) and bradykinesia (p = 0.008). Diminished 1N/4R isoform expression was also associated with H2 (p = 0.001). CONCLUSIONS: Overall, our results suggest that age and H2 are associated with higher parkinsonism score and decreased total MAPT RNA expression. Additionally, we found that H2 and parkinsonism are associated with altered expression levels of specific isoforms. These findings may contribute to the understanding of the association between MAPT locus and parkinsonism in elderly subjects and in some extent to age-related neurodegenerative diseases.

Epigenomics of Alzheimer's disease.

Alzheimer's disease (AD) is a large and growing public health problem. It is characterized by the accumulation of amyloid ß peptides and abnormally phosphorylated tau proteins that are associated with cognitive decline and dementia. Much has been learned about the genomics of AD from linkage analyses and, more recently, genome-wide association studies. Several but not all aspects of the genomic landscape are involved in amyloid ß metabolism. The moderate concordance of disease among twins suggests other factors, potentially epigenomic factors, are related to AD. We are at the earliest stages of examining the relation of the epigenome to the clinical and pathologic phenotypes that characterize AD. Our literature review suggests that there is some evidence of age-related changes in human brain methylation. Unfortunately, studies of AD have been relatively small with limited coverage of methylation sites and microRNA, let alone other epigenomic marks. We are in the midst of 2 large studies of human brains including coverage of more than 420,000 autosomal cytosine-guanine dinucleotides with the Illumina Infinium HumanMethylation450 BeadArray, and histone acetylation with chromatin immunoprecipitation sequencing. We present descriptive data to help inform other researchers what to expect from these approaches to better design and power their studies. We then discuss future directions to inform on the epigenomic architecture of AD.

Association of Brain DNA methylation in SORL1, ABCA7, HLA-DRB5, SLC24A4, and BIN1 with pathological diagnosis of Alzheimer disease.

IMPORTANCE: Recent large-scale genome-wide association studies have discovered several genetic variants associated with Alzheimer disease (AD); however, the extent to which DNA methylation in these AD loci contributes to the disease susceptibility remains unknown. OBJECTIVE: To examine the association of brain DNA methylation in 28 reported AD loci with AD pathologies. DESIGN, SETTING, AND PARTICIPANTS: Ongoing community-based clinical pathological cohort studies of aging and dementia (the Religious Orders Study and the Rush Memory and Aging Project) among 740 autopsied participants 66.0 to 108.3 years old. EXPOSURES: DNA methylation levels at individual CpG sites generated from dorsolateral prefrontal cortex tissue using a bead assay. MAIN OUTCOMES AND MEASURES: Pathological diagnosis of AD by National Institute on Aging-Reagan criteria following a standard postmortem examination. RESULTS: Overall, 447 participants (60.4%) met the criteria for pathological diagnosis of AD. Brain DNA methylation in SORL1, ABCA7, HLA-DRB5, SLC24A4, and BIN1 was associated with pathological AD. The association was robustly retained after replacing the binary trait of pathological AD with 2 quantitative and molecular specific hallmarks of AD, namely, Aß load and paired helical filament tau tangle density. Furthermore, RNA expression of transcripts of SORL1 and ABCA7 was associated with paired helical filament tau tangle density, and the expression of BIN1 was associated with Aß load. CONCLUSIONS AND RELEVANCE: Brain DNA methylation in multiple AD loci is associated with AD pathologies. The results provide further evidence that disruption of DNA methylation is involved in the pathological process of AD.

Association of DNA methylation in the brain with age in older persons is confounded by common neuropathologies.

DNA methylation plays a crucial role in the regulation of gene expression, cell differentiation and development. Previous studies have reported age-related alterations of methylation levels in the human brain across the lifespan, but little is known about whether the observed association with age is confounded by common neuropathologies among older persons. Using genome-wide DNA methylation data from 740 postmortem brains, we interrogated 420,132 CpG sites across the genome in a cohort of individuals with ages from 66 to 108 years old, a range of ages at which many neuropathologic indices become quite common. We compared the association of DNA methylation prior to and following adjustment for common neuropathologies using a series of linear regression models. In the simplest model adjusting for technical factors including batch effect and bisulfite conversion rate, we found 8156 CpGs associated with age. The number of CpGs associated with age dropped by more than 10% following adjustment for sex. Notably, after adjusting for common neuropathologies, the total number of CpGs associated with age was reduced by approximately 40%, compared to the sex-adjusted model. These data illustrate that the association of methylation changes in the brain with age is inflated if one does not account for age-related brain pathologies. This article is part of a Directed Issue entitled: Epigenetics dynamics in development and disease.

Polarization of the effects of autoimmune and neurodegenerative risk alleles in leukocytes.

To extend our understanding of the genetic basis of human immune function and dysfunction, we performed an expression quantitative trait locus (eQTL) study of purified CD4(+) T cells and monocytes, representing adaptive and innate immunity, in a multi-ethnic cohort of 461 healthy individuals. Context-specific cis- and trans-eQTLs were identified, and cross-population mapping allowed, in some cases, putative functional assignment of candidate causal regulatory variants for disease-associated loci. We note an over-representation of T cell-specific eQTLs among susceptibility alleles for autoimmune diseases and of monocyte-specific eQTLs among Alzheimer's and Parkinson's disease variants. This polarization implicates specific immune cell types in these diseases and points to the need to identify the cell-autonomous effects of disease susceptibility variants.