Nicole M. Davis Armstrong1
Wei-Min Chen2,3
Michael S. Brewer1
Stephen R. Williams4
Michèle M. Sale2,3
Bradford B. Worrall3,4
Keith L. Keene1,5*- 1Department of Biology, East Carolina University, Greenville, NC, United States
- 2Center for Public Health Genomics, University of Virginia, Charlottesville, VA, United States
- 3Department of Public Health Sciences, University of Virginia, Charlottesville, VA, United States
- 4Department of Neurology, University of Virginia, Charlottesville, VA, United States
- 5Center for Health Disparities, East Carolina University, Greenville, NC, United States
DNA methylation, a well-characterized epigenetic modification that is influenced by both environment and genetic variation, has previously been implicated in a number of complex diseases, including cardiovascular disease and stroke. The goal of this study was to evaluate epigenome-wide associations with recurrent stroke and the folate one-carbon metabolism-related trait, plasma homocysteine (hcy). Differential methylation analyses were performed on 473,864 autosomal CpG loci, using Illumina HumanMethylation 450K array data in 180 ischemic stroke cases from the Vitamin Intervention for Stroke Prevention (VISP) clinical trial. Linear regression was used to assess associations between number of strokes prior to VISP enrollment and measures of hcy with degree of methylation (β-values), while logistic regression was used to evaluate recurrent stroke status and incident recurrent stroke associations. All regression analyses were stratified by race. Two differentially methylated CpG sites exceeded epigenome-wide significance (p ≤ 1.055 × 10−7) for prior number of strokes (PNS) in European Americans. The top locus, cg22812874, was located in the ankyrin repeat and SOCS box containing 10 gene (ASB10; p = 3.4 × 10−9; β = −0.0308; 95% CI = −0.040, −0.002). Methylation locus cg00340919, located in an intron of the tetratricopeptide repeat domain 37 gene, was also statistically significant (TTC37; p = 8.74 × 10−8; β = −0.0517; 95% CI = −0.069, −0.034). An additional 138 CpG sites met our threshold for suggestive significance (p ≤ 5 × 10−5). We evaluated DNA methylation associated with recurrent stroke and hcy phenotypes across the epigenome. Hypermethylation at two CpG sites located in ASB10 and TTC37 was associated with fewer strokes prior to VISP enrollment. Our findings present a foundation for additional epigenome-wide studies, as well as mechanistic studies into epigenetic marks that influence recurrent stroke risk.
Introduction
Epigenetics refers to chemical modifications in DNA that influence and regulate gene expression without altering the sequence. In mammalian cells, DNA methylation occurs with a covalent addition of a methyl group to a cytosine primarily located at a CpG dinucleotide (Li et al., 1992; Robertson, 2011) and is influenced by both environmental and genetic factors (Jirtle and Skinner, 2007; Mathers et al., 2010). Aberrant DNA methylation patterns have been implicated in a number of diseases (Robertson, 2005). The FOCM pathway plays a crucial role in DNA methylation through the conversion of methionine to homocysteine (hcy) with the aid of B vitamins (Fox and Stover, 2008). Elevated levels of plasma hcy have been associated with an increased risk of stroke (Boushey et al., 1995; Wald et al., 2005; Shi et al., 2015).
Annually, stroke contributes to 6.5 million deaths worldwide and is the fifth leading cause of death in the United States (Benjamin et al., 2017). There are roughly 795,000 reported cases of stroke in the United States each year and nearly a quarter are recurrent attacks, which are often more disabling and fatal compared to an initial stroke. Approximately 85% of strokes are ischemic attacks; a complex disorder with both genetic and environmental contributions (Benjamin et al., 2017). After adjusting for the major risk factors contributing to ischemic stroke, an estimated 69% of the population-attributable risk is not accounted for (Jamrozik et al., 1994; Benjamin et al., 2017), suggesting genetic and epigenetic factors may account for a substantial proportion of stroke risk (Jamrozik et al., 1994; Qureshi and Mehler, 2010).
A more mechanistic understanding of the relationship between the FOCM pathway, hcy, DNA methylation, and vascular disease could help identify prevention strategies for stroke. This study performed epigenome-wide associations with recurrent stroke and FOCM-related phenotypes, using methylation data of 180 participants from the VISP clinical trial. Our findings implicate two novel methylation loci associated with the number of strokes suffered prior to trial enrollment.
Materials and Methods
Ethics Approval
The institutional review boards (IRBs) of Wake Forest School of Medicine and the University of North Carolina at Chapel Hill School of Medicine approved the VISP clinical trial study protocol. The local IRB for each of the individual recruiting sites approved the VISP protocol and all participants provided written, informed consent. The genetic and epigenetic components of the VISP study were approved by the IRBs at the University of Virginia School of Medicine and East Carolina University.
Subjects
Vitamin Intervention for Stroke Prevention was a multicenter, double-blind, randomized, and controlled clinical trial designed to determine whether folic acid, pyridoxine (B6), and cyanocobalamin (B12) supplementation, in addition to best medical management, prevented recurrent cerebral infarction and reduced the risk for myocardial infarction or fatal coronary heart disease in participants aged 35 years or older. VISP participants with baseline hcy levels above the 25th percentile were enrolled within 120 days of suffering a NDCI characterized by the sudden onset of a neurologic deficit persisting for at least 24 h and not due to a cardiac embolism (Toole, 2002). Patients (n = 3,860) were randomly assigned to receive daily doses of either the high-dose (n = 1,827; 25 mg B6, 0.4 mg B12, and 2.5 mg folic acid) or the low-dose formulation (n = 1,853; 200 μg B6, 6 μg B12, and 20 μg folic acid). A subset of 2,100 VISP participants provided consent for inclusion in genetic studies, of which methylation data was generated for 204 participants. Although all VISP participants had elevated baseline hcy levels, the 204 participants from the methylation subset were selected at the extremes (highest and lowest) for hcy levels. After filtering, quality control, and normalization, methylation data from 180 individuals were used in subsequent analyses (104 of European descent and 76 of African descent).
FOCM Metabolite Measurements
Baseline fasting and 2 h post methionine-load hcy levels were determined at randomization for all VISP participants. Each participant fasted overnight, gave a whole blood sample for baseline measurements and was subsequently given methionine (0.1 g/kg body weight) as a crystalline powder dissolved in unsweetened juice. Two hours after consumption, another blood sample was taken to determine post-methionine hcy level (Toole et al., 2004; Pettigrew et al., 2008). The absolute difference between the baseline and post-methionine load hcy (DeltaPost) was calculated and provides a measurement of an individual’s ability to convert methionine into homocysteine.
Methylation Data Generation
DNA methylation data was generated using the Illumina Infinium HumanMethylation450K BeadChip microarray (Illumina, Inc., San Diego, CA, United States) per the manufacturer’s protocol. Denaturation and sodium bisulfite conversion of genomic DNA was performed using Zymo EZ DNA Methylation Kits (Zymo Research Corporation, Irvine, CA, United States). The BCD was eluted from the Zymo-Spin column, amplified, enzymatically fragmented, precipitated, and resuspended. This solution was hybridized to the BeadChip for 22 h at 48°C and BeadChips were subsequently washed to remove unhybridized DNA samples. Single base extensions of hybridized primers labeled with biotin underwent staining. The BeadChips were washed, coated, and imaged on an iScan (Illumina, Inc., San Diego, CA, United States) (Maksimovic et al., 2012).
Imaged data was viewed using GenomeStudio (Illumina, Inc., San Diego, CA, United States). The degree of methylation, or β-value, was calculated as the ratio of methylated allele intensity divided by the sum of the intensities of methylated and unmethylated loci for 485,512 CpG sites. Quality control and quantile normalization were performed using the Bioconductor Minfi and watermelon packages, respectively (Touleimat and Tost, 2012; Pidsley et al., 2013; Aryee et al., 2014). Removal of methylation loci located on sex chromosomes resulted in 473,864 autosomal methylation sites that were included in subsequent analyses. Raw DNA methylation data, linked to genome-wide association SNP data, will be available at the NCBI database of Genotypes and Phenotypes (dbGaP Study Accession: phs000343.v3.p1).
Statistical Analyses
Principal components analysis was performed on SNP data acquired from VISP genetic studies using KING software (Manichaikul et al., 2010). The top four and ten principal components were used as covariates in regression analyses to account for population substructure and admixture in EA and AA, respectively.
Cellular heterogeneity due to variation in cell population proportions in whole blood (Reinius et al., 2012) was controlled for using cell proportion estimates based upon the Houseman algorithm (Houseman et al., 2012). Estimated cell counts were generated for B-lymphocytes, CD4+ and CD8+ T-lymphocytes, natural killer cells, granulocytes, and monocytes using the “estimateCellCounts()” function in the R package minfi (Aryee et al., 2014; Jaffe and Irizarry, 2014) and were used as covariates in regression models, as described below.
To identify differentially methylated CpG sites associated with measures of baseline plasma hcy (hcy), post-methionine load hcy (POST), and the absolute difference (DeltaPost), a multiple linear regression stratified by race and adjusting for age, sex, principal components, batch, and estimated cell proportions was performed in R (version 3.3.1).
To identify CpG sites associated with recurrent stroke (experiencing more than one stroke in one’s lifetime, “Recurrent Ever”; incident stroke during the VISP trial, “VISP Recurrent”), a logistic regression stratified by race and adjusting for age, sex, principal components, batch, estimated cell proportions, and trial arm (high- or low-dose supplementation) was performed. Differentially methylated CpG sites associated with the number of strokes prior to enrollment (PNS) were analyzed using a linear regression model adjusted for age, sex, principal components, batch, and estimated cell proportions. Two individuals (one AA and one EA) were excluded in this model due to dependent variable outliers.
Significant associations were corrected for multiple testing using a Bonferroni-adjusted epigenome-wide threshold of 1.055 × 10−7. The beta coefficients and 95% confidence intervals (CIs) for 473,864 methylation loci were calculated in R (Supplementary Tables S1, S2). Manhattan plots were generated using qqman (Turner, 2018) and lambda values were estimated using GenABEL (Aulchenko et al., 2007).
Statistical power was calculated for both AA and EA subpopulations at an effect size of 0.15, significance level of α = 0.05, and degrees of freedom corresponding to race subpopulation size using the pwr package (Cohen, 1988; Champely, 2017). The AA subpopulation resulted in limited power of 0.327, while the EA subpopulation had power of 0.643.
Gene Ontology Term Enrichment
Gene ontology (GO) term enrichment was performed using the GOrilla (Gene Ontology enRIchment anaLysis and visuaLizAtion) tool (Eden et al., 2007, 2009) and the missMethyl package in R (Phipson et al., 2016). The top 86 unique genes associated with differentially methylated CpG sites were ranked and compared to a background list of 15,057 terms for biological processes, and 4,518 terms for molecular function in the GOrilla program (p ≤ 3.321 × 10−6 and p ≤ 1.107 × 10−5, respectively). The same subset of unique genes was used in missMethyl across 20,732 terms (p ≤ 2.411 × 10−6).
Results
DNA methylation data was generated for a subset of 180 individuals from the VISP trial (AA n = 76; EA n = 104; Table 1), using the Illumina 450K array. Linear and logistic regression for recurrent stroke phenotypes (PNS, VISP Recurrent, Recurrent Ever) identified two statistically significant (p ≤ 1.055 × 10−7) CpG sites in EA (Table 2 and Figure 1). The most significant association was observed at cg22812874 (p = 3.40 × 10−9, β = −0.0308; 95% CI = −0.040, −0.002), located in exon 6 of the ankyrin repeat and SOCs box containing 10 gene, ASB10, on chromosome 7. This locus is also found in a TFBS for the CCCTC-binding factor, CTCF. Methylation locus cg00340919, located within intron 4 of the tetratricopeptide repeat domain 37 gene, TTC37, on chromosome 5, also exceeded epigenome-wide significance (p = 8.74 × 10−8, β = −0.0517; 95% CI = −0.069, −0.034). A subset of 82 CpG sites reached an a priori suggestive significance threshold of p ≤ 5 × 10−5 (Supplementary Figures S7–S11 and Supplementary Table S1).
TABLE 1. Demographic summary statistics.
TABLE 2. Statistically significant (p ≤ 1.055 × 10−7) differentially methylated sites associated with the previous number of strokes in European Americans.
FIGURE 1. Manhattan plot depicting results from linear regression of previous number of strokes in European American population. The red line indicates a Bonferroni-adjusted threshold of 1.055 × 10−7; blue line indicates suggestive significance threshold of 5.0 × 10−7. Inflation value, = 1.31.
Regression analyses of hcy measures for both AA and EA (hcy, POST, DeltaPost) identified 56 CpG sites that met the suggestive threshold (Supplementary Figures S1–S6 and Supplementary Table S2). For all hcy related analyses, the two most significant associations were observed at cg08657886 (AA POST; p = 3.26 × 10−7) and cg21649277 (EA POST; p = 4.69 × 10−7), both located on chromosome 19. Methylation locus cg08657886 was located within the promoter region of the Family with sequence similarity 32, member A (FAM32A) gene, while cg21649277 is located in the first intron of the arrestin domain containing 2 gene, ARRDC2 (Supplementary Table S2). Suggestive evidence of association for the POST phenotype at cg14818960, located within the first intron of the folate receptor 1 gene, FOLR1 (p = 9.66 × 10−6) is also notable, given this gene’s function in the FOCM pathway (Supplementary Table S2).
Gene ontology analysis using GOrilla, identified 22 terms with suggestive association (p ≤ 1.0 × 10−3) using the top 86 CpG sites from the PNS regression (Supplementary Table S3). Sulfur compound binding (GO:1901681) was the most significant GO term (p = 5.01 × 10−5), followed by heparin binding (GO:0008201; p = 5.38 × 10−5). The genes ASB10 and TTC37 were both associated at the suggestive threshold with the terms signal transduction (GO:0007165) and (unknown) molecular function (GO:0003674), albeit both of these terms are broadly defined. missMethyl identified 55 suggestive GO terms at p ≤ 5 × 10−3, including proteinaceous ECM (GO:0005578, p = 6.15 × 10−5) and adenylate cyclase-inhibiting G-protein coupled receptor signaling pathway (GO:0007193, p = 1.20 × 10−4) (Supplementary Table S4).
Discussion
Differential methylation analyses were performed for six phenotypes related to hcy and recurrent stroke, resulting in 140 significant or suggestive associations. Two epigenome-wide significant, differentially methylated CpG sites located within the ASB10 and TTC37 genes were associated with PNS in the EA subgroup. These two methylation loci were not statistically significant in the AA subset analyses for PNS (cg22812874 p = 0.554, β = −2.445; 95% CI = −10.682, 5.791; cg00340919 p = 0.397, β = 0.841; 95% CI = −1.135, 2.818).
The ASB10 gene functions by targeting SOCS proteins and the elongin BC regulatory complex for degradation (Kile et al., 2002). There are no studies implicating this gene with stroke risk, however, a synonymous variant in ASB10 has been associated with primary open angle glaucoma (Pasutto et al., 2012). Immunofluorescence and confocal analyses determined that ASB10 is localized intracellularly in human TM cells (Keller et al., 2013) and that TNF-α and IL-1β up-regulate ASB10 in these cells (Keller and Wirtz, 2017). TNF-α and IL-1β are risk factors for ischemic stroke in humans (Luheshi et al., 2011; Bokhari et al., 2014) and result in neurotoxicity in rats after an ischemic event (Clausen et al., 2008). As previously mentioned, the associated methylation locus cg22812874 is located in a TFBS for the CTCF protein as predicted using ORegAnno data from UCSC Genome Browser (Kent et al., 2002; Montgomery et al., 2006; Griffith et al., 2008). This protein is ubiquitously expressed and has been described traditionally as a transcriptional repressor (Lobanenkov et al., 1990; Holwerda and de Laat, 2013). Additionally, CTCF can act as an insulator when positioned between an enhancer and promoter. The insulator activity blocks signaling between the enhancer and promoter and prevents transcriptional activation (Bell et al., 1999). CTCF was shown to be essential in cardiogenesis, as well as in the mediation of cardiomyocyte differentiation. Inactivation of this transcription factor causes severe cardiac defects and death in embryonic development (Gomez-Velazquez et al., 2017). The ASB10 association is likely specific to populations of European decent, as our AA analyses showed no evidence for association. A recent study of ischemic cardiomyopathy in a Spanish population detected significant differential DNA methylation patterns for another member of the ASB family (ASB1 gene), and further highlighted the potential role of DNA methylation of ASB family members for stroke and cardiovascular phenotypes (Ortega et al., 2018).
TTC37 encodes a protein with twenty tetratricopeptide repeats, which is highly expressed in vascular tissues (Fabre et al., 2011) and has been implicated in THES. Symptoms of THES include platelet abnormalities, such as large platelet sizes and thrombocytosis (Hartley et al., 2010), which are risk factors for stroke (Arboix et al., 1995; Posfai et al., 2014; Kato et al., 2015). While the associated methylation locus cg00340919 is not directly located in a regulatory element, it is adjacent to the SNP rs73777218. Combined eQTL analyses from GTEx suggests that rs73777218 is associated with the expression of pseudogene CTD-2538A21.1 in brain-cerebellum tissues (GTEx Consortium, 2013, 2017 Hormozdiari et al., 2014; Ongen et al., 2016). Although limited knowledge is available on the function of this particular pseudogene, it is located within intron 7 of the FAM81B gene. Expression of FAM81B is related to progression free survival of estrogen positive breast cancer when treated with aromatase inhibitors (Ramirez-Ardila et al., 2016). Although debatable, aromatase inhibitors are thought to increase risk of coronary heart disease and stroke in women (Chlebowski et al., 2006).
The most differentially methylated loci in the hcy analyses were cg08657886 and cg21649277. cg08657886 is located within the promoter region of FAM32A, a housekeeping gene whose protein may function in inducing G2 arrest and apoptosis (The UniProt Consortium, 2017). cg21649277 is located in the first intron of ARRDC2, whose protein is a member of the arrestin-related trafficking (ART) family of proteins (Rauch and Martin-Serrano, 2011). ART proteins functions as regulators of G protein-mediated signaling. The most notable locus that reached suggestive significance with functional implications in FOCM or cerebral diseases was cg14818960. This locus is located within FOLR1, the gene encoding the folate receptor-α protein, which regulates folate transport into cells. Folate receptor-α is found within cell membranes, where it binds to 5-MTHF (Steinfeld et al., 2009). Variants in FOLR1 have been independently associated with cerebral folate transport deficiency and brain abnormalities caused by the lack of folate (Steinfeld et al., 2009).
GO analyses, using GOrilla and missMethyl, identified 22 and 55 terms, respectively, that reached a suggestive threshold. The most statistically significant term across both packages was sulfur compound binding including but not limited to organosulfur compounds diallyl sulfide and allyl mercaptan. These compounds have been recognized to decrease total cholesterol, triglycerides, thrombosis, and lower blood pressure (Warshafsky et al., 1993; Kris-Etherton et al., 2002), all of which are risk factors of cardiovascular disease.
Both GOrilla and missMethyl identified heparin binding as a molecular function with a suggestive association. Heparin, a commonly prescribed anticoagulant, is used in ischemic stroke prevention. Studies focusing on the efficacy and safety of heparin alone or in combination therapies have demonstrated that heparin is safe and effective in ischemic attacks (Yi et al., 2014; Dluha et al., 2016); however, novel therapeutics are becoming more promising and could soon be implemented more regularly in favor to traditional heparin treatment (Douketis et al., 2015; Han et al., 2015). The ontology terms that were suggestively associated with ASB10 and TTC37 included signal transduction. Signal transduction mediators have been reported in brain ischemia injury and depending on the mediator, have both neuroprotective and neurotoxic effects (Jin et al., 2013).
missMethyl enrichment identified proteinaceous ECM and G-protein-coupled receptors as suggestive terms. The proteinaceous ECM is a layer that consists predominately of proteins and glycosaminoglycans (Binns et al., 2009; QuickGO, 2017). ECM composition was reported to be altered upon BBB disruption, which directly affects neurological disease progression (Baeten and Akassoglou, 2011). Degradation of these ECM components by MMPs has been reported (Mott and Werb, 2004; Cai et al., 2015). Matrix metalloproteinase-9 (MMP-9) has been associated with ischemic stroke progression in mouse brain models (Zhao et al., 2006). MMP-9 upregulated in the first 3 days after ischemic stroke contributed to BBB permeability and brain inflammation. However, elevated levels of MMP-9 on days 7–14 post-stroke played a beneficial role in angiogenesis and brain recovery (Zhao et al., 2006). Therefore, strategies that modulate MMP function in regulation of ECM components can be utilized for ischemic stroke recovery. In addition, G-protein-coupled receptors have been implicated in neuronal function and protection in human studies, cell-culture, and animal models (Martin et al., 2005). Possible roles for class II G-protein receptors have been revealed in the pathogenesis of neurological and neurodegenerative conditions (Martin et al., 2005).
This study has a number of strengths and limitations. The VISP clinical trial included a mixture of both AA and EA participants and was well-phenotyped from a clinical perspective. Although specific stroke subtypes were not included, based on the inclusion/exclusion criteria, the strokes likely represented a homogenous, small-vessel disease group. Recurrent stroke is a grossly understudied phenotype; therefore the inclusion of recurrent stroke phenotypes is a strength of this study. However, no public datasets are available with both epigenetic (DNA methylation) data and stroke recurrence phenotypes. This limitation prohibits us from validating our findings, yet further emphasizes the need for more studies focused on these traits/datasets. The use of the commercially available Illumina arrays allowed the entire epigenome to be evaluated with confidence and reproducibility.
The VISP clinical trial had a randomized, controlled, double-blind experimental design, and thus does not represent a population sample. Study entry required hcy levels in the top quartile, while the methylation subset included VISP participants with hcy levels at the extremes. All subjects experienced an ischemic stroke within 120 days of enrollment, thus our sample is skewed toward an older age distribution with hypertension and elevated body mass indexes. Despite its relevance for a stroke recurrence phenotype, our hcy results may not be generalizable. We are unable to determine the timing of the epigenetic changes relative to the phenotypes investigated, thus we can only establish association and not causation.
DNA from VISP participants was extracted from whole blood samples. Although not optimal due to cellular heterogeneity, whole blood provides a valuable resource that is potentially available for replication studies. To minimize confounding effects due to cell heterogeneity in whole blood, cell-type proportions were estimated based on methylation signatures and were included as covariates in the regression models. The Houseman method, as implemented in minfi, is generally regarded as the gold standard for cell proportion estimations in DNA methylation studies using whole blood. Alternative programs exist that are capable of generating cell estimations and corresponding cell proportions. Using MeDeCom (Lutsik et al., 2017), we compared latent DNA methylation components (LMCs) with estimated cell proportions from minfi (Supplementary Figure S12 and Supplementary Table S5). Monocytes and NK cells are not fully captured in the LMC estimations. Monocytes are thought to play a major role in ischemic stroke pathology (ElAli and LeBlanc, 2016), therefore the minfi cell proportion estimations which account for monocyte proportions were considered most biologically appropriate, and were utilized in our analysis. Despite having limited power, statistically significant associations between two CpG sites and PNS in EA were detected, but larger, independent sample sizes will be needed to confirm and validate our results. Overall, these analyses present a foundation for further epigenetic investigations of novel loci in stroke-susceptible populations.
Data Availability Statement
The datasets used and analyzed during this study are available from the corresponding author on reasonable request.
Author Contributions
NDA contributed in data analysis, interpretation of the data, and drafting of the manuscript. W-MC contributed in data analysis, interpretation, and revising of the manuscript. MB contributed in data manipulation, quality control and analysis, and revising of the manuscript. SW contributed in the study design and revising of the manuscript. MS and BW contributed in overall conception, design and interpretation of the data, and revising of the manuscript. KK contributed in overall conception, design and interpretation of the data, drafting and revising of the manuscript. All authors read and approved the final manuscript.
Funding
This study was supported by the National Institutes of Health National Human Genome Research Institute (U01HG005160; Supplement U01HG005160S1-MMS/BBW), the American Heart Association Scientist Development Award (KK), and the Walking for My Life 5K Stroke Walk sponsored by Faith Christian Center International (Charlottesville, VA, United States). Study recruitment and datasets for the VISP clinical trial were supported by a research grant (RO1 NS34447; PI James Toole) from the National Institute of Neurological Disorders and Stroke. This study was also supported by American Heart Association Grant: 12SDG9180012.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We would like to thank all individuals who volunteered and participated in the VISP study. We would also like to thank the physicians, study-site investigators, UVA Genome Sciences Laboratory, and technical staff involved in data generation and subsequent analyses.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2018.00358/full#supplementary-material
Abbreviations
5-MTHF, 5-methyl-tetrahydrofolate; AA, African American; ARRDC2, arrestin domain containing 2 gene; ASB10, Ankyrin repeat and SOCs box containing 10; BBB, blood–brain barrier; BCD, bisulfite-converted DNA; DeltaPost, Absolute difference of post-methionine load test and baseline homocysteine levels; EA, European American; ECM, extracellular matrix; eQTL, expression quantitative trait loci; FAM32A, Family with sequence similarity 32, member A gene; FOCM, folate-mediated one-carbon metabolism; FOLR1, folate receptor alpha; GO, gene ontology; GTEx, Genotype-Tissue Expression; GWAS, genome-wide association study; Hcy, Homocysteine; HCY, Measured baseline homocysteine level; IL-1β, interleukin-1 beta; MMP, matrix metalloproteinases; MTHFR, Methylenetetrahydrofolate reductase; NDCI, non-disabling cerebral infarction; PCA, principal component analysis; PNS, previous number of strokes; POST, Post-methionine load test homocysteine level; Recurrent Ever, Recurrent stroke status prior to or during VISP trial; SNP, single nucleotide polymorphism; SOCs, suppressor of cytokine signaling; TFBS, transcription factor binding site; THES, trichohepatoenteric syndrome; TM cells, trabecular meshwork cells; TNF-α, tumor necrosis factor alpha; TTC37, tetratricopeptide repeat domain 37 gene; VISP, Vitamin Intervention for Stroke Prevention; VISP Recurrent, Recurrent stroke during the VISP trial.
References
Arboix, A., Bessess, C., Acin, P., Massons, J. B., Florensa, L., Oliveres, M., et al. (1995). Ischemic stroke as first manifestation of essential thrombocythemia. Report of six cases. Stroke 26, 1463–1466.
Aryee, M. J., Jaffee, A. E., Corrada-Bravo, H., Ladd-Acosta, C., Feinberg, A. P., Hansen, K. D., et al. (2014). Minfi: a flexible and comprehensive bioconductor package for the analysis of infinium DNA methylation microarrays. Bioinformatics 30, 1363–1369. doi: 10.1093/bioinformatics/btu049
Aulchenko, Y. S., Ripke, S., Isaacs, A., and van Duijn, C. M. (2007). GenABEL: an R package for genome-wide association analysis. Bioinformatics 23, 1294–1296. doi: 10.1093/bioinformatics/btm108
Baeten, K. M., and Akassoglou, K. (2011). Extracellular matrix and matrix receptors in blood-brain barrier formation and stroke. Dev. Neurobiol. 71, 1018–1039. doi: 10.1002/dneu.20954
Bell, A. C., West, A. G., and Felsenfeld, G. (1999). The protein CTCF is required for the enhancer blocking activity of vertebrate insulators. Cell 98, 387–396. doi: 10.1016/S0092-8674(00)81967-4
Benjamin, E. J., Blaha, M. J., Chiuve, S. E., Cushman, M., Das, S. R., Deo, R., et al. (2017). Heart disease and stroke statistics-2017 update: a report from the american heart association. Circulation. 135, e146–e603. doi: 10.1161/CIR.0000000000000485
Binns, D., Dimmer, E., Huntley, R., Barrel, D., O’Donovan, C., and Apweiler, R. (2009). QuickGO: a web-based tool for gene ontology searching. Bioinformatics 25, 3045–3046. doi: 10.1093/bioinformatics/btp536
Bokhari, F. A., Shakoori, T. A., Butt, A., and Ghafoor, F. (2014). TNF-alpha: a risk factor for ischemic stroke. J. Ayub. Med. Coll. Abbottabad. 26, 111–114.
Boushey, C. J., Beresford, S. A., Omenn, G. S., and Motulsky, A. G. (1995). A quantitative assessment of plasma homocysteine as a risk for vascular disease. Probable benefits of increasing folic acid intakes. JAMA 274, 1049–1057. doi: 10.1001/jama.1995.03530130055028
Cai, H., Mu, Z., Jian, Z., Wang, Y., Yang, G. Y., and Zhang, Z. (2015). Hypoxia-controlled matrix metalloproteinase-9 hyperexpression promotes behavioral recovery after ischemia. Neurosci. Bull. 31, 550–560. doi: 10.1007/s12264-015-1533-1
Champely, S. (2017). Package ‘pwr’. Available at: http://cran.r-project.org/package=pwr
Chlebowski, R. T., Anderson, G. L., Geller, M., and Col, N. (2006). Coronary heart disease and stroke with aromatase inhibitor, tamoxifen, and menopausal hormone therapy use. Clin. Breast Cancer 6, S58–S64. doi: 10.3816/CBC.2006.s.005
Clausen, B. H., Lambertsen, K. L., Babcock, A. A., Holm, T. H., Dagnaes-Hansen, F., and Finsen, B. (2008). Interleukin-1beta and tumor necrosis factor-alpha are expressed by different subsets of microglia and macrophages after ischemic stroke in mice. J. Neuroinflam. 5:46. doi: 10.1186/1742-2094-5-46
Cohen, J. (1988). Statistical Power Analysis for the Behavioral Sciences, 2nd Edn. Hillsdale, NJ: Lawrence Erlbaum Associates.
Dluha, J., Sivak, S., Kurca, E., Dusenka, R., Kalmarova, K., Turcanova Koprusakova, M., et al. (2016). The safety and efficacy of Heparin and Nadroparin compared to placebo in acute ischemic stroke-pilot study. Biomed. Pap. Med. Fac. Univ. Palacky Olomouc Czech. Repub. 160, 543–548. doi: 10.5507/bp.2016.042
Douketis, J. D., Spyropoulos, A. C., Kaatz, S., Becker, R. C., Caprini, J. A., Dunn, A. S., et al. (2015). Perioperative bridging anticoagulation in patients with atrial fibrillation. N Engl J Med. 373, 823–833. doi: 10.1056/NEJMoa1501035
Eden, E., Lipson, D., Yogev, S., and Yakhini, Z. (2007). Discovering motifs in ranked lists of DNA sequences. PLoS Comp. Biol. 3:e39. doi: 10.1371/journal.pcbi.0030039
Eden, E., Navon, R., Steinfeld, I., Lipson, D., and Yakhini, Z. (2009). GOrilla: a tool for discovery and visualization of enriched GO terms in ranked gene lists. BMC Bioinformatics 10:48. doi: 10.1186/1471-2105-10-48
ElAli, A., and LeBlanc, N. J. (2016). The role of monocytes in ischemic stroke pathobiology: new avenues to explore. Front. Aging Neurosci. 8:29. doi: 10.3389/fnagi.2016.00029
Fabre, A., Martinex-Vinson, C., Roquelaure, B., Missirian, C., Andre, N., Breton, A., et al. (2011). Novel mutations in TTC37 associated with tricho-hepato-enteric syndrome. Hum. Mutat. 32, 277–281. doi: 10.1002/humu.21420
Fox, J. T., and Stover, P. T. (2008). Folate-mediated one-carbon metabolism. Vitam. Horm. 79, 1–44. doi: 10.1016/S0083-6729(08)00401-9
Gomez-Velazquez, M., Badia-Careaga, C., Lechuga-Vieco, A. V., Nieto-Arrellano, R., Tena, J. J., Rollan, I., et al. (2017). CTCF counter-regulates cardiomyocyte development and maturation programs in the embryonic heart. PLoS Genet. 13:e1006985. doi: 10.1371/journal.pgen.1006985
Griffith, O. L., Montgomery, S. B., Bernier, B., Chu, B., Kasaian, K., Aerts, S., et al. (2008). ORegAnno: an open-access community-driven resource for regulatory annotation. Nucleic Acids Res 36, D107–D113.
GTEx Consortium. (2013). The Genotype-Tissue Expression (GTEx) project. Nat. Genet. 45, 580–585. doi: 10.1038/ng.2653
GTEx Consortium. (2017). Genetic effects on gene expression across human tissues. Nature 550, 204–213. doi: 10.1038/nature24277
Han, Y., Guo, J., Zheng, Y., Zang, H., Su, X., Wang, Y., et al. (2015). Bivalirduin vs heparin with or without tirofiban during primary percutaneous coronary intervention in acute myocardial infarction: the BRIGHT randomized clinical trial. JAMA 313, 1336–1346. doi: 10.1001/jama.2015.2323
Hartley, J. L., Zachos, N. C., Dawood, B., Donoqitz, M., Forman, J., Politt, R. J., et al. (2010). Mutations in TTC37 cause trichohepatoenteric syndrome (phenotypic diarrhea of infancy). Gastroenterology 138, 2388–2398. doi: 10.1053/j.gastro.2010.02.010
Holwerda, S. J., and de Laat, W. (2013). CTCF: the protein, the binding partners, the binding sites and their chromatin loops. Philos. Trans. R. Soc. Lond. B Biol. Sci. 368:20120369. doi: 10.1098/rstb.2010.0369
Hormozdiari, F., Kostem, E., Kang, E. Y., Pasaniuc, B., and Eskin, E. (2014). Identifying causal variants at loci with multiple signals of association. Genetics 198, 497–508. doi: 10.1534/genetics.114.167908
Houseman, E. A., Accomando, W. P., Koestler, D. C., Christensen, B. C., Marsit, C. J., Nelson, H. H., et al. (2012). DNA methylation arrays as surrogate measures of cell mixture distribution. BMC Bioinformatics 13:86. doi: 10.1186/1471-2105-13-86
Jaffe, A. E., and Irizarry, R. A. (2014). Accounting for cellular heterogeneity is critical in epigenome-wide association studies. Genome Biol. 15:R31. doi: 10.1186/gb-2014-15-2-r31
Jamrozik, K., Broadhurst, R. J., Anderson, C. S., and Stewart-Wynne, E. G. (1994). The role of lifestyle factors in the etiology of stroke. A population-based case-control study in Perth, Western Australia. Stroke. 25, 51–59. doi: 10.1161/01.STR.25.1.51
Jin, R., Liu, L., Zhang, S., Nanda, A., and Li, G. (2013). Role of inflammation and its mediators in acute ischemic stroke. J. Cardiovasc. Transl. Res. 6, 834–851. doi: 10.1007/s12265-013-9508-6
Jirtle, R. L., and Skinner, M. K. (2007). Environmental epigenomics and disease susceptibility. Nat. Rev. Genet. 8, 253–262. doi: 10.1038/nrg2045
Kato, Y., Hayashi, T., Sehara, Y., Deguchi, I., Fukuonka, T., Maruyama, H., et al. (2015). Ischemic stroke with essential thrombocythemia: a case series. J. Stroke Cerebrovasc. Dis. 24, 890–893. doi: 10.1016/j.jstrokecerebrovasdis.2014.12.012
Keller, K. E., and Wirtz, M. K. (2017). Working your SOCS off: the role of ASB10 and protein degradation pathways in glaucoma. Exp. Eye Res. 158, 154–160. doi: 10.1016/j.exer.2016.06.003
Keller, K. E., Yang, Y. F., Sun, Y. Y., Sykes, R., Acott, T. S., and Wirtz, M. K. (2013). Ankyrin repeat and suppressor of cytokine signaling box containing protein-10 is associated with ubiquitin-mediated degradation pathways in trabecular meshwork cells. Mol. Vis. 19, 1639–1655.
Kent, W. J., Sugnet, C. W., Furey, T. S., Roskin, K. M., Pringle, T. H., Zahler, A. M., et al. (2002). The human genome browser at UCSC. Genome Res. 12, 996–1006. doi: 10.1101/gr.229102
Kile, B. T., Schulman, B. A., Alexander, W. S., Nicola, N. A., Martin, H. M., and Hilton, D. J. (2002). The SOCS box: a tale of destruction and degradation. Trends Biochem. Sci. 27, 235–241. doi: 10.1016/S0968-0004(02)02085-6
Kris-Etherton, P. M., Hecker, K. D., Bonanome, A., Coval, S. M., Binkoski, A. E., Hilpert, K. F., et al. (2002). Bioactive compounds in foods: their role in the prevention of cardiovascular disease and cancer. Am. J. Med. 113, 71–88. doi: 10.1016/S0002-9343(01)00995-0
Li, E., Bestor, T. H., and Jaenisch, R. (1992). Targeted mutation of the DNA methyltransferase gene results in embryonic lethality. Cell 69, 915–926. doi: 10.1016/0092-8674(92)90611-F
Lobanenkov, V. V., Nicolas, R. H., Adler, V. V., Paterson, H., Klenova, E. M., Polotskaja, A. V., et al. (1990). A novel sequence-specific DNA binding protein which interacts with three regularly spaced direct repeats of the CCCTC-motif in the 5’-flanking sequence of the chicken c-myc gene. Oncogene 5, 1743–1753.
Luheshi, N. M., Kovacs, K. J., Lopez-Castejon, G., Brough, D., and Denes, A. (2011). Interleukin-1α expression precedes IL-1β after ischemic brain injury and is localised to areas of focal neuronal loss and penumbral tissues. J. Neuroinflam. 8:186. doi: 10.1186/1742-2094-8-186
Lutsik, P., Slawski, M., Gasparoni, G., Vedeneev, N., Hein, M., and Walter, J. (2017). MeDeCom: discovery and quantification of latent components of heterogeneous methylomes. Genome Biol. 18:55. doi: 10.1186/s13059-017-1182-6
Maksimovic, J., Gordon, L., and Oshlack, A. (2012). SWAN: subset-quantile within array normalization for Illumina Infinium HumanMethylation450 BeadChips. Genome Biol. 13:R44. doi: 10.1186/gb-2012-13-6-r44
Manichaikul, A., Mychaleckyj, J. C., Rich, S. S., Daly, K., Sale, M., and Chen, W. M. (2010). Robust relationship inference in genome-wide association studies. Bioinformatics 26, 2867–2873. doi: 10.1093/bioinformatics/btq559
Martin, B., Lopez de Maturana, R., Brenneman, R., Walent, T., Mattson, M. P., and Maudsley, S. (2005). Class II G protein-coupled receptors and their ligands in neuronal function and protection. Neuromol. Med. 7, 3–36. doi: 10.1385/NMM:7:1-2:003
Mathers, J. C., Strathdee, G., and Relton, C. L. (2010). Induction of epigenetic alterations by dietary and other environmental factors. Adv. Genet. 71, 3–39. doi: 10.1016/B978-0-12-380864-6.00001-8
Montgomery, S. B., Griffith, O. L., Sleumer, M. C., Bergman, C. M., Bilenky, M., Pleasance, E. D., et al. (2006). ORegAnno: an open access database and curation system for literature-derived promoters, transcription factor binding sites and regulatory variation. Bioinformatics 22, 637–640. doi: 10.1093/bioinformatics/btk027
Mott, J. D., and Werb, Z. (2004). Regulation of matrix biology by matrix metalloproteinases. Curr. Opin. Cell Biol. 16, 448–464. doi: 10.1016/j.ceb.2004.07.010
Ongen, H., Buil, A., Brown, A. A., Dermitzakis, E. T., and Delaneau, O. (2016). Fast and efficient QTL mapper for thousands of molecular phenotypes. Bioinformatics. 32, 1479–1485. doi: 10.1093/bioinformatics/btv722
Ortega, A., Tarazon, E., Gil-Cayuela, C., Martinez-Dolz, L., Lago, F., Gonzalez-Juanatey, J., et al. (2018). ASB1 differential methylation in ischaemic cardiomyopathy:relationship with left ventricular performance in end-stage heart failure patients. ESC Heart Fail 5, 732–737. doi: 10.1002/ehf2.12289
Pasutto, F., Keller, K. E., Weisschih, N., Sticht, H., Samples, J. R., Yang, Y. F., et al. (2012). Variants in ASB10 are associated with open-angle glaucoma. Hum. Mol. Genet. 21, 1336–1349. doi: 10.1093/hmg/ddr572
Pettigrew, L. C., Bang, H., Chambless, L. E., Howard, V. J., and Toole, J. F. (2008). VISP Investigators. Assessment of pre- and post-methionine load homocysteine for prediction of recurrent stroke and coronary artery disease in the vitamin intervention for stroke prevention trial. Atherosclerosis 200, 345–349. doi: 10.1016/j.atherosclerosis.2007.11.014
Phipson, B., Maksimovic, J., and Oshlack, A. (2016). missMethyl: an R package for analyzing data from illumina’s human methylation 450 platform. Bioinformatics 32, 286–288. doi: 10.1093/bioinformatics/btv560
Pidsley, R., Wong, Y. C. C., Volta, M., Lunnon, K., Mill, J., and Schalkwyk, L. C. (2013). A data-driven approach to preprocessing Illumina 450K methylation array data. BMC Genomics 14:293. doi: 10.1186/1471-2164-14-293
Posfai, E., Marton, I., Szoke, A., Borbenyi, Z., Vecsei, L., Csomor, A., et al. (2014). Stroke in essential thrombocythemia. J. Neurol. Sci. 336, 260–262. doi: 10.1016/j.jns.2013.10.016
QuickGO. (2017)). Embl-EBI. Available at: http://www.ebi.ac.uk/QuickGO/
Qureshi, I. A., and Mehler, M. F. (2010). Emerging role of epigenetics in stroke: part 1: DNA methylation and chromatin modifications. Arch. Neurol. 67, 1316–1322. doi: 10.1001/archneurol.2010.275
Ramirez-Ardila, D. E., Ruigro-Ritstier, K., Helmijr, J. C., Look, M. P., van Laere, S., Dirix, L., et al. (2016). LRG1 mRNA expression in breast cancer associates with PIK3CA genotype and with aromatase inhibitor therapy outcome. Mol Oncol. 10, 1363–1373. doi: 10.1016/j.molonc.2016.07.004
Rauch, S., and Martin-Serrano, J. (2011). Multiple interactions between the ESCRT machinery and arrestin-related proteins: implications for PPXY-dependent budding. J. Virol. 85, 3546–3556. doi: 10.1128/JVI.02045-10
Reinius, L. E., Acevedo, N., Joerink, M., Pershagen, G., Dahlen, S.-E., Greco, D., et al. (2012). Differential DNA methylation in purified human blood cells: implications for cell lineage and studies on disease susceptibility. PLoS One 7:e41361. doi: 10.1371/journal.pone.0041361
Robertson, K. D. (2005). DNA methylation and human disease. Nat. Rev. Genet. 6, 597–610. doi: 10.1038/ngr1655
Robertson, K. D. (2011). DNA methylation, methyltransferases, and cancer. Oncogene 20, 3139–3155. doi: 10.1038/sj.onc.1204341
Shi, Z., Guan, Y., Hua, Y. R., Liu, S., Zhang, M., Lu, H., et al. (2015). Elevated total homocysteine levels in acute ischemic stroke are associated with long-term mortality. Stroke 46, 2419–2425. doi: 10.1161/STROKEAHA.115.009136
Steinfeld, R., Grapp, M., Kraetzner, R., Dreha-Kulaczewski, S., Helms, G., Dechent, P., et al. (2009). Folate receptor alpha defect causes cerebral folate transport deficiency: a treatable neurodegenerative disorder associated with disturbed myelin metabolism. Am. J. Hum. Genet. 85, 354–363. doi: 10.1016/j.ajhg.2009.08.005
The UniProt Consortium. (2017). UniProt: the universal protein knowledgebase. Nucl. Acids Res. 45, D158–D169. doi: 10.1093/nar/gkw1099
Toole, J. F. (2002). Vitamin intervention for stroke prevention. J. Neurol. Sci. 20, 121–124. doi: 10.1016/S0022-510X(02)00265-4
Toole, J. F., Malinow, M. R., Chambless, L. E., Spence, J. D., Pettigrew, L. C., Howard, V. J., et al. (2004). Lowering homocysteine in patients with ischemic stroke to prevent recurrent stroke, myocardial infarction, and death. The vitamin intervention for stroke prevention (VISP) randomized controlled trial. JAMA 291, 565–575. doi: 10.1001/jama.291.5.565
Touleimat, N., and Tost, J. (2012). Complete pipeline for infiunium® human methylation 450K beadchip data processing using subset quantile normalization for accurate DNA methylation estimation. Epigenomics 4, 325–341. doi: 10.2219/epi.12.21
Turner, S. D. (2018). qqman: an R package for visualizing GWAS results using Q-Q and manhattan plots. J. Open Source Softw. 3:731. doi: 10.21105/joss.00731
Wald, D. S., Law, M., and Morris, J. K. (2005). Homocysteine and cardiovascular disease: evidence on causality from a meta-analysis. BMJ. 325, 1202.
Warshafsky, S., Kamer, R. S., and Sivak, S. L. (1993). Effect of garlic on total serum cholesterol. A meta-analysis. Ann. Intern. Med. 119, 599–605. doi: 10.7326/0003-4819-119-7_Part_1-199310010-00009
Yi, X., Lin, J., Wang, C., Zhang, B., and Chi, W. (2014). Low-molecular-weight heparin is more effective than aspirin in preventing early neurologic deterioration and improving six-month outcome. J. Stroke Cerebrovasc. Dis. 23, 1537–1544. doi: 10.1016/j.jstrokecerebrovasdis.2013.12.036
Keywords: DNA methylation, recurrent stroke, VISP, association, folate one carbon metabolism, homocysteine, epigenome, epigenetics
Citation: Davis Armstrong NM, Chen W-M, Brewer MS, Williams SR, Sale MM, Worrall BB and Keene KL (2018) Epigenome-Wide Analyses Identify Two Novel Associations With Recurrent Stroke in the Vitamin Intervention for Stroke Prevention Clinical Trial. Front. Genet. 9:358. doi: 10.3389/fgene.2018.00358
Received: 03 April 2018; Accepted: 20 August 2018;
Published: 06 September 2018.
Edited by:
Ryan Yuen, Hospital for Sick Children, CanadaReviewed by:
Shicheng Guo, Marshfield Clinic Research Institute, United StatesYanfeng Zhang, HudsonAlpha Institute for Biotechnology, United States
Copyright © 2018 Davis Armstrong, Chen, Brewer, Williams, Sale, Worrall and Keene. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Keith L. Keene, keenek@ecu.edu