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.

Rheumatoid arthritis-associated RBPJ polymorphism alters memory CD4+ T cells.

Notch signaling has recently emerged as an important regulator of immune responses in autoimmune diseases. The recombination signal-binding protein for immunoglobulin kappa J region (RBPJ) is a transcriptional repressor, but converts into a transcriptional activator upon activation of the canonical Notch pathway. Genome-wide association studies of rheumatoid arthritis (RA) identified a susceptibility locus, rs874040(CC), which implicated the RBPJ gene. Here, chromatin state mapping generated using the chromHMM algorithm reveals strong enhancer regions containing DNase I hypersensitive sites overlapping the rs874040 linkage disequilibrium block in human memory, but not in naïve CD4(+) T cells. The rs874040 overlapping this chromatin state was associated with increased RBPJ expression in stimulated memory CD4(+) T cells from healthy subjects homozygous for the risk allele (CC) compared with memory CD4(+) T cells bearing the protective allele (GG). Transcriptomic analysis of rs874040(CC) memory T cells showed a repression of canonical Notch target genes IL (interleukin)-9, IL-17 and interferon (IFN)γ in the basal state. Interestingly, activation of the Notch pathway using soluble Notch ligand, Jagged2-Fc, induced IL-9 and IL-17A while delta-like 4Fc, another Notch ligand, induced higher IFNγ expression in the rs874040(CC) memory CD4(+) T cells compared with their rs874040(GG) counterparts. In RA, RBPJ expression is elevated in memory T cells from RA patients compared with control subjects, and this was associated with induced inflammatory cytokines IL-9, IL-17A and IFNγ in response to Notch ligation in vitro. These findings demonstrate that the rs874040(CC) allele skews memory T cells toward a pro-inflammatory phenotype involving Notch signaling, thus increasing the susceptibility to develop RA.

Water-deficit-induced changes in transcription factor expression in maize seedlings.

Plants tolerate water deficits by regulating gene networks controlling cellular and physiological traits to modify growth and development. Transcription factor (TF)-directed regulation of transcription within these gene networks is key to eliciting appropriate responses. In this study, reverse transcription quantitative PCR (RT-qPCR) was used to examine the abundance of 618 transcripts from 536 TF genes in individual root and shoot tissues of maize seedlings grown in vermiculite under well-watered (water potential of -0.02 MPa) and water-deficit conditions (water potentials of -0.3 and -1.6 MPa). A linear mixed model identified 433 TF transcripts representing 392 genes that differed significantly in abundance in at least one treatment, including TFs that intersect growth and development and environmental stress responses. TFs were extensively differentially regulated across stressed maize seedling tissues. Hierarchical clustering revealed TFs with stress-induced increased abundance in primary root tips that likely regulate root growth responses to water deficits, possibly as part of abscisic acid and/or auxin-dependent signaling pathways. Ten of these TFs were selected for validation in nodal root tips of drought-stressed field-grown plants (late V1 to early V2 stage). Changes in abundance of these TF transcripts under a field drought were similar to those observed in the seedling system.

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.

Proteome-wide analysis reveals G protein-coupled receptor-like proteins in rice (Oryza sativa).

G protein-coupled receptors (GPCRs) constitute the largest family of transmembrane proteins in metazoans that mediate the regulation of various physiological responses to discrete ligands through heterotrimeric G protein subunits. The existence of GPCRs in plant is contentious, but their comparable crucial role in various signaling pathways necessitates the identification of novel remote GPCR-like proteins that essentially interact with the plant G protein α subunit and facilitate the transduction of various stimuli. In this study, we identified three putative GPCR-like proteins (OsGPCRLPs) (LOC_Os06g09930.1, LOC_Os04g36630.1, and LOC_Os01g54784.1) in the rice proteome using a stringent bioinformatics workflow. The identified OsGPCRLPs exhibited a canonical GPCR 'type I' 7TM topology, patterns, and biologically significant sites for membrane anchorage and desensitization. Cluster-based interactome mapping revealed that the identified proteins interact with the G protein α subunit which is a characteristic feature of GPCRs. Computational results showing the interaction of identified GPCR-like proteins with G protein α subunit and its further validation by the membrane yeast-two-hybrid assay strongly suggest the presence of GPCR-like 7TM proteins in the rice proteome. The absence of a regulator of G protein signaling (RGS) box in the C- terminal domain, and the presence of signature motifs of canonical GPCR in the identified OsGPCRLPs strongly suggest that the rice proteome contains GPCR-like proteins that might be involved in signal transduction.

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.

Homolog-specific PCR primer design for profiling splice variants.

To study functional diversity of proteins encoded from a single gene, it is important to distinguish the expression levels among the alternatively spliced variants. A variant-specific primer pair is required to amplify each alternatively spliced variant individually. For this purpose, we developed a new feature, homolog-specific primer design (HSPD), in our high-throughput primer and probe design software tool, PRIMEGENS-v2. The algorithm uses a de novo approach to design primers without any prior information of splice variants or close homologs for an input query sequence. It not only designs primer pairs but also finds potential isoforms and homologs of the input sequence. Efficiency of this algorithm was tested for several gene families in soybean. A total of 187 primer pairs were tested under five different abiotic stress conditions with three replications at three time points. Results indicate a high success rate of primer design. Some primer pairs designed were able to amplify all splice variants of a gene. Furthermore, by utilizing combinations within the same multiplex pool, we were able to uniquely amplify a specific variant or duplicate gene. Our method can also be used to design PCR primers to specifically amplify homologs in the same gene family. PRIMEGENS-v2 is available at: http://primegens.org.

Targeted bisulfite sequencing by solution hybrid selection and massively parallel sequencing.

We applied a solution hybrid selection approach to the enrichment of CpG islands (CGIs) and promoter sequences from the human genome for targeted high-throughput bisulfite sequencing. A single lane of Illumina sequences allowed accurate and quantitative analysis of ~1 million CpGs in more than 21,408 CGIs and more than 15,946 transcriptional regulatory regions. Of the CpGs analyzed, 77-84% fell on or near capture probe sequences; 69-75% fell within CGIs. More than 85% of capture probes successfully yielded quantitative DNA methylation information of targeted regions. Differentially methylated regions (DMRs) were identified in the 5'-end regulatory regions, as well as the intra- and intergenic regions, particularly in the X-chromosome among the three breast cancer cell lines analyzed. We chose 46 candidate loci (762 CpGs) for confirmation with PCR-based bisulfite sequencing and demonstrated excellent correlation between two data sets. Targeted bisulfite sequencing of three DNA methyltransferase (DNMT) knockout cell lines and the wild-type HCT116 colon cancer cell line revealed a significant decrease in CpG methylation for the DNMT1 knockout and DNMT1, 3B double knockout cell lines, but not in DNMT3B knockout cell line. We demonstrated the targeted bisulfite sequencing approach to be a powerful method to uncover novel aberrant methylation in the cancer epigenome. Since all targets were captured and sequenced as a pool through a series of single-tube reactions, this method can be easily scaled up to deal with a large number of samples.

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.

Harnessing peripheral DNA methylation differences in the Alzheimer's Disease Neuroimaging Initiative (ADNI) to reveal novel biomarkers of disease.

BACKGROUND: Alzheimer's disease (AD) is a chronic progressive neurodegenerative disease impacting an estimated 44 million adults worldwide. The causal pathology of AD (accumulation of amyloid-beta and tau), precedes hallmark symptoms of dementia by more than a decade, necessitating development of early diagnostic markers of disease onset, particularly for new drugs that aim to modify disease processes. To evaluate differentially methylated positions (DMPs) as novel blood-based biomarkers of AD, we used a subset of 653 individuals with peripheral blood (PB) samples in the Alzheimer's disease Neuroimaging Initiative (ADNI) consortium. The selected cohort of AD, mild cognitive impairment (MCI), and age-matched healthy controls (CN) all had imaging, genetics, transcriptomics, cerebrospinal protein markers, and comprehensive clinical records, providing a rich resource of concurrent multi-omics and phenotypic information on a well-phenotyped subset of ADNI participants. RESULTS: In this manuscript, we report cross-diagnosis differential peripheral DNA methylation in a cohort of AD, MCI, and age-matched CN individuals with longitudinal DNA methylation measurements. Epigenome-wide association studies (EWAS) were performed using a mixed model with repeated measures over time with a P value cutoff of 1 × 10-5 to test contrasts of pairwise differential peripheral methylation in AD vs CN, AD vs MCI, and MCI vs CN. The most highly significant differentially methylated loci also tracked with Mini Mental State Examination (MMSE) scores. Differentially methylated loci were enriched near brain and neurodegeneration-related genes (e.g., BDNF, BIN1, APOC1) validated using the genotype tissue expression project portal (GTex). CONCLUSIONS: Our work shows that peripheral differential methylation between age-matched subjects with AD relative to healthy controls will provide opportunities to further investigate and validate differential methylation as a surrogate of disease. Given the inaccessibility of brain tissue, the PB-associated methylation marks may help identify the stage of disease and progression phenotype, information that would be central to bringing forward successful drugs for AD.

Meta-Analysis of the Alzheimer's Disease Human Brain Transcriptome and Functional Dissection in Mouse Models.

We present a consensus atlas of the human brain transcriptome in Alzheimer's disease (AD), based on meta-analysis of differential gene expression in 2,114 postmortem samples. We discover 30 brain coexpression modules from seven regions as the major source of AD transcriptional perturbations. We next examine overlap with 251 brain differentially expressed gene sets from mouse models of AD and other neurodegenerative disorders. Human-mouse overlaps highlight responses to amyloid versus tau pathology and reveal age- and sex-dependent expression signatures for disease progression. Human coexpression modules enriched for neuronal and/or microglial genes broadly overlap with mouse models of AD, Huntington's disease, amyotrophic lateral sclerosis, and aging. Other human coexpression modules, including those implicated in proteostasis, are not activated in AD models but rather following other, unexpected genetic manipulations. Our results comprise a cross-species resource, highlighting transcriptional networks altered by human brain pathophysiology and identifying correspondences with mouse models for AD preclinical studies.

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 Major Risk Factors for Alzheimer's Disease: Age, Sex, and Genes Modulate the Microglia Response to Aß Plaques.

Gene expression profiles of more than 10,000 individual microglial cells isolated from cortex and hippocampus of male and female AppNL-G-F mice over time demonstrate that progressive amyloid-ß accumulation accelerates two main activated microglia states that are also present during normal aging. Activated response microglia (ARMs) are composed of specialized subgroups overexpressing MHC type II and putative tissue repair genes (Dkk2, Gpnmb, and Spp1) and are strongly enriched with Alzheimer's disease (AD) risk genes. Microglia from female mice progress faster in this activation trajectory. Similar activated states are also found in a second AD model and in human brain. Apoe, the major genetic risk factor for AD, regulates the ARMs but not the interferon response microglia (IRMs). Thus, the ARMs response is the converging point for aging, sex, and genetic AD risk factors.

Reproductive period and epigenetic modifications of the oxidative phosphorylation pathway in the human prefrontal cortex.

PURPOSE: Human females have a unique duration of post-reproductive longevity, during which sex-specific mechanisms ma influence later-life mechanisms of neuronal resilience and vulnerability. The maintenance of energy metabolism, through the oxidative phosphorylation (OXPHOS) apparatus, is essential for brain health. Given the known association between reproductive period (years from menarche to menopause) and cognitive aging, we examined the hypothesis that cumulative estrogen exposure across the lifetime may be associated with differential methylation of genes in the OXPHOS pathway. METHODS: Using DNA methylation patterns in the post-mortem dorsolateral prefrontal cortex (DLPFC) of 426 women prospectively followed until death in the Religious Orders Study and Rush Memory and Aging Project, we examined the relationship between reproductive period (subtracting age at menarche from age at menopause) and DNA methylation of a published set of autosomal OXPHOS genes previously implicated in stroke susceptibility. We then performed an unsupervised analysis of methylation levels across the Hallmark pathways from the Molecular Signatures Database. RESULTS: We observed a strong association between reproductive period and DNA methylation status across OXPHOS CpGs. We replicated this association between reproductive period and DNA methylation in a much larger set of OXPHOS genes in our unsupervised analysis. Here, reproductive period also showed associations with methylation in genes related to E2F, MYC and MTORC1 signaling, fatty acid metabolism and DNA repair. CONCLUSION: This study provides evidence from both a supervised and unsupervised analyses, that lifetime cumulative endogenous steroid exposures may play a role in maintenance of post-menopausal cellular balance, including in brain tissue.

Genome-scale probe and primer design with PRIMEGENS.

This chapter introduces the software package PRIMEGENS for designing gene-specific probes and associated PCR primers on a large scale. Such design is especially useful for constructing cDNA or oligo microarray to minimize cross-hybridization. PRIMEGENS can also be used for designing primers to amplify a segment of a unique target gene using reverse-transcriptase (RT)-PCR. The input to PRIMEGENS is a set of sequences, whose primers need to be designed, and a sequence pool containing all the genes in a genome. It provides options to choose various parameters. PRIMEGENS uses a systematic algorithm for designing gene-specific probes and its primer pair. For a given sequence, PRIMEGENS first searches for the longest gene-specific fragment and then designs best PCR product for this fragment. The 2.0 version of PRIMEGENS provides a graphical user interface (GUI) with additional features. The software is freely available for any users and can be downloaded from 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.

PRIMEGENSw3: a web-based tool for high-throughput primer and probe design.

Highly specific and efficient primer and probe design has been a major hurdle in many high-throughput techniques. Successful implementation of any PCR or probe hybridization technique depends on the quality of primers and probes used in terms of their specificity and cross-hybridization. Here we describe PRIMEGENSw3, a set of web-based utilities for high-throughput primer and probe design. These utilities allow users to select genomic regions and to design primer/probe for selected regions in an interactive, user-friendly, and automatic fashion. The system runs the PRIMEGENS algorithm in the back-end on the high-performance server with the stored genomic database or user-provided custom database for cross-hybridization check. Cross-hybridization is checked not only using BLAST but also by checking mismatch positions and energy calculation of potential hybridization hits. The results can be visualized online and also can be downloaded. The average success rate of primer design using PRIMEGENSw3 is ~90 %. The web server also supports primer design for methylated sequences, which is used in epigenetic studies. Stand-alone version of the software is also available for download at the website.

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.