DNA methylation and stroke prognosis: an epigenome-wide association study.

BACKGROUND AND AIMS: Stroke is the leading cause of adult-onset disability. Although clinical factors influence stroke outcome, there is a significant variability among individuals that may be attributed to genetics and epigenetics, including DNA methylation (DNAm). We aimed to study the association between DNAm and stroke prognosis. METHODS AND RESULTS: To that aim, we conducted a two-phase study (discovery-replication and meta-analysis) in Caucasian patients with ischemic stroke from two independent centers (BasicMar [discovery, N = 316] and St. Pau [replication, N = 92]). Functional outcome was assessed using the modified Rankin Scale (mRS) at three months after stroke, being poor outcome defined as mRS > 2. DNAm was determined using the 450K and EPIC BeadChips in whole-blood samples collected within the first 24 h. We searched for differentially methylated positions (DMPs) in 370,344 CpGs, and candidates below p-value < 10-5 were subsequently tested in the replication cohort. We then meta-analyzed DMP results from both cohorts and used them to identify differentially methylated regions (DMRs). After doing the epigenome-wide association study, we found 29 DMPs at p-value < 10-5 and one of them was replicated: cg24391982, annotated to thrombospondin-2 (THBS2) gene (p-valuediscovery = 1.54·10-6; p-valuereplication = 9.17·10-4; p-valuemeta-analysis = 6.39·10-9). Besides, four DMRs were identified in patients with poor outcome annotated to zinc finger protein 57 homolog (ZFP57), Arachidonate 12-Lipoxygenase 12S Type (ALOX12), ABI Family Member 3 (ABI3) and Allantoicase (ALLC) genes (p-value < 1·10-9 in all cases). DISCUSSION: Patients with poor outcome showed a DMP at THBS2 and four DMRs annotated to ZFP57, ALOX12, ABI3 and ALLC genes. This suggests an association between stroke outcome and DNAm, which may help identify new stroke recovery mechanisms.

Role of PATJ in stroke prognosis by modulating endothelial to mesenchymal transition through the Hippo/Notch/PI3K axis.

Through GWAS studies we identified PATJ associated with functional outcome after ischemic stroke (IS). The aim of this study was to determine PATJ role in brain endothelial cells (ECs) in the context of stroke outcome. PATJ expression analyses in patient's blood revealed that: (i) the risk allele of rs76221407 induces higher expression of PATJ, (ii) PATJ is downregulated 24 h after IS, and (iii) its expression is significantly lower in those patients with functional independence, measured at 3 months with the modified Rankin scale ((mRS) ≤2), compared to those patients with marked disability (mRS = 4-5). In mice brains, PATJ was also downregulated in the injured hemisphere at 48 h after ischemia. Oxygen-glucose deprivation and hypoxia-dependent of Hypoxia Inducible Factor-1α also caused PATJ depletion in ECs. To study the effects of PATJ downregulation, we generated PATJ-knockdown human microvascular ECs. Their transcriptomic profile evidenced a complex cell reprogramming involving Notch, TGF-ß, PI3K/Akt, and Hippo signaling that translates in morphological and functional changes compatible with endothelial to mesenchymal transition (EndMT). PATJ depletion caused loss of cell-cell adhesion, upregulation of metalloproteases, actin cytoskeleton remodeling, cytoplasmic accumulation of the signal transducer C-terminal transmembrane Mucin 1 (MUC1-C) and downregulation of Notch and Hippo signaling. The EndMT phenotype of PATJ-depleted cells was associated with the nuclear recruitment of MUC1-C, YAP/TAZ, ß-catenin, and ZEB1. Our results suggest that PATJ downregulation 24 h after IS promotes EndMT, an initial step prior to secondary activation of a pro-angiogenic program. This effect is associated with functional independence suggesting that activation of EndMT shortly after stroke onset is beneficial for stroke recovery.

Genome-wide association and Mendelian randomization study of fibroblast growth factor 21 reveals causal associations with hyperlipidemia and possibly NASH.

BACKGROUND: Fibroblast growth factor 21 (FGF21) is a hepatokine that produces metabolic benefits, such as improvements of lipid profile. We performed a genome-wide association study (GWAS) to identify genetic variants associated with circulating FGF21 and investigated the causal effects of FGF21 on pertinent outcomes using Mendelian randomization (MR). METHODS: We conducted a GWAS testing ∼7.8 million DNA sequence variants with circulating FGF21 in a discovery cohort of 6259 Swedish adults with replication in 4483 Swedish women. We then performed two-sample MR analyses of genetically predicted circulating FGF21 in relation to alcohol and nutrient intake, cardiovascular and metabolic biomarkers and diseases, and liver function biomarkers using publicly available GWAS summary statistics data. RESULTS: Our GWAS identified multiple single-nucleotide polymorphisms with genome-wide significant associations (P < 5 × 10-8) with circulating FGF21 on chromosomes 2 and 19 in or near the GCKR and FGF21 genes, respectively. The strongest signal at the FGF21 locus (rs2548957, ß = 0.181, P < 2.18 × 10-42) displayed in two-sample MR analyses robust associations with lower alcohol intake, lower circulating low-density lipoprotein cholesterol, apolipoprotein B, C-reactive protein, gamma-glutamyl transferase, and galectin-3 concentrations, and higher circulating insulin-like growth factor-I and alkaline phosphatase concentrations after correcting for multiple testing (P < 0.0018) whereas associations with fat mass, type 2 diabetes, and cardiovascular disease were largely null. CONCLUSIONS: We identified robust associations of certain genetic variants in or near the GCKR and FGF21 genes with circulating FGF21 concentrations. Furthermore, our results support a strong causal effect of FGF21 on improved lipid profile, reduced alcohol consumption and C-reactive protein concentrations, and liver function biomarkers including fibrosis. We found largely null or weak positive associations with fat mass, diabetes, and cardiovascular disease as well as higher insulin-like growth factor-I concentrations, which could indicate a compensatory increase to regulate the above FGF21 resistant states in humans.

Multi-ancestry GWAS reveals excitotoxicity associated with outcome after ischaemic stroke.

During the first hours after stroke onset, neurological deficits can be highly unstable: some patients rapidly improve, while others deteriorate. This early neurological instability has a major impact on long-term outcome. Here, we aimed to determine the genetic architecture of early neurological instability measured by the difference between the National Institutes of Health Stroke Scale (NIHSS) within 6 h of stroke onset and NIHSS at 24 h. A total of 5876 individuals from seven countries (Spain, Finland, Poland, USA, Costa Rica, Mexico and Korea) were studied using a multi-ancestry meta-analyses. We found that 8.7% of NIHSS at 24 h of variance was explained by common genetic variations, and also that early neurological instability has a different genetic architecture from that of stroke risk. Eight loci (1p21.1, 1q42.2, 2p25.1, 2q31.2, 2q33.3, 5q33.2, 7p21.2 and 13q31.1) were genome-wide significant and explained 1.8% of the variability suggesting that additional variants influence early change in neurological deficits. We used functional genomics and bioinformatic annotation to identify the genes driving the association from each locus. Expression quantitative trait loci mapping and summary data-based Mendelian randomization indicate that ADAM23 (log Bayes factor = 5.41) was driving the association for 2q33.3. Gene-based analyses suggested that GRIA1 (log Bayes factor = 5.19), which is predominantly expressed in the brain, is the gene driving the association for the 5q33.2 locus. These analyses also nominated GNPAT (log Bayes factor = 7.64) ABCB5 (log Bayes factor = 5.97) for the 1p21.1 and 7p21.1 loci. Human brain single-nuclei RNA-sequencing indicates that the gene expression of ADAM23 and GRIA1 is enriched in neurons. ADAM23, a presynaptic protein and GRIA1, a protein subunit of the AMPA receptor, are part of a synaptic protein complex that modulates neuronal excitability. These data provide the first genetic evidence in humans that excitotoxicity may contribute to early neurological instability after acute ischaemic stroke.

Causal Effect of MMP-1 (Matrix Metalloproteinase-1), MMP-8, and MMP-12 Levels on Ischemic Stroke: A Mendelian Randomization Study.

Background and Purpose: MMP (matrix metalloproteinase) levels have been widely associated with ischemic stroke risk and poststroke outcome. However, their role as a risk factor or as a subeffect because of ischemia is uncertain. Methods: We performed a literature search of genome-wide studies that evaluate serum/plasma levels of MMPs. We used a 2-sample Mendelian randomization approach to evaluate the causality of MMP levels on ischemic stroke risk or poststroke outcome, using 2 cohorts: MEGASTROKE (n=440 328) and GODs (n=1791). Results: Genome-wide association studies of MMP-1, MMP-8, and MMP-12 plasma/serum levels were evaluated. A significant association, which was also robust in the sensitivity analysis, was found with all ischemic strokes: MMP-12 (odds ratio=0.90 [95% CI, 0.86­0.94]; q value=7.43×10−5), and with subtypes of stroke, large-artery atherosclerosis: MMP-1 (odds ratio=0.95 [95% CI, 0.92­0.98]; q value=0.01) and MMP-12 (odds ratio=0.71 [95% CI, 0.65­0.77]; q value=5.11×10−14); small-vessel occlusion: MMP-8 (odds ratio=1.24 [95% CI, 1.06­1.45]; q value=0.03). No associations were found in relation to stroke outcome. Conclusions: Our study suggests a causal link between lower serum levels of MMP-12 and the risk of ischemic stroke, lower serum levels of MMP-1 and MMP-12 and the risk of large-artery stroke and higher serum levels of MMP-8 and the risk of lacunar stroke.

Single nucleotide variations in ZBTB46 are associated with post-thrombolytic parenchymal haematoma.

Haemorrhagic transformation is a complication of recombinant tissue-plasminogen activator treatment. The most severe form, parenchymal haematoma, can result in neurological deterioration, disability, and death. Our objective was to identify single nucleotide variations associated with a risk of parenchymal haematoma following thrombolytic therapy in patients with acute ischaemic stroke. A fixed-effect genome-wide meta-analysis was performed combining two-stage genome-wide association studies (n = 1904). The discovery stage (three cohorts) comprised 1324 ischaemic stroke individuals, 5.4% of whom had a parenchymal haematoma. Genetic variants yielding a P-value < 0.05 1 × 10-5 were analysed in the validation stage (six cohorts), formed by 580 ischaemic stroke patients with 12.1% haemorrhagic events. All participants received recombinant tissue-plasminogen activator; cases were parenchymal haematoma type 1 or 2 as defined by the European Cooperative Acute Stroke Study (ECASS) criteria. Genome-wide significant findings (P < 5 × 10-8) were characterized by in silico functional annotation, gene expression, and DNA regulatory elements. We analysed 7 989 272 single nucleotide polymorphisms and identified a genome-wide association locus on chromosome 20 in the discovery cohort; functional annotation indicated that the ZBTB46 gene was driving the association for chromosome 20. The top single nucleotide polymorphism was rs76484331 in the ZBTB46 gene [P = 2.49 × 10-8; odds ratio (OR): 11.21; 95% confidence interval (CI): 4.82-26.55]. In the replication cohort (n = 580), the rs76484331 polymorphism was associated with parenchymal haematoma (P = 0.01), and the overall association after meta-analysis increased (P = 1.61 × 10-8; OR: 5.84; 95% CI: 3.16-10.76). ZBTB46 codes the zinc finger and BTB domain-containing protein 46 that acts as a transcription factor. In silico studies indicated that ZBTB46 is expressed in brain tissue by neurons and endothelial cells. Moreover, rs76484331 interacts with the promoter sites located at 20q13. In conclusion, we identified single nucleotide variants in the ZBTB46 gene associated with a higher risk of parenchymal haematoma following recombinant tissue-plasminogen activator treatment.

Multi-ancestry genetic study in 5,876 patients identifies an association between excitotoxic genes and early outcomes after acute ischemic stroke.

During the first hours after stroke onset neurological deficits can be highly unstable: some patients rapidly improve, while others deteriorate. This early neurological instability has a major impact on long-term outcome. Here, we aimed to determine the genetic architecture of early neurological instability measured by the difference between NIH stroke scale (NIHSS) within six hours of stroke onset and NIHSS at 24h (ΔNIHSS). A total of 5,876 individuals from seven countries (Spain, Finland, Poland, United States, Costa Rica, Mexico and Korea) were studied using a multi-ancestry meta-analyses. We found that 8.7% of ΔNIHSS variance was explained by common genetic variations, and also that early neurological instability has a different genetic architecture than that of stroke risk. Seven loci (2p25.1, 2q31.2, 2q33.3, 4q34.3, 5q33.2, 6q26 and 7p21.1) were genome-wide significant and explained 2.1% of the variability suggesting that additional variants influence early change in neurological deficits. We used functional genomics and bioinformatic annotation to identify the genes driving the association from each loci. eQTL mapping and SMR indicate that ADAM23 (log Bayes Factor (LBF)=6.34) was driving the association for 2q33.3. Gene based analyses suggested that GRIA1 (LBF=5.26), which is predominantly expressed in brain, is the gene driving the association for the 5q33.2 locus. These analyses also nominated PARK2 (LBF=5.30) and ABCB5 (LBF=5.70) for the 6q26 and 7p21.1 loci. Human brain single nuclei RNA-seq indicates that the gene expression of ADAM23 and GRIA1 is enriched in neurons. ADAM23 , a pre-synaptic protein, and GRIA1 , a protein subunit of the AMPA receptor, are part of a synaptic protein complex that modulates neuronal excitability. These data provides the first evidence in humans that excitotoxicity may contribute to early neurological instability after acute ischemic stroke. RESEARCH INTO CONTEXT: Evidence before this study: No previous genome-wide association studies have investigated the genetic architecture of early outcomes after ischemic stroke.Added Value of this study: This is the first study that investigated genetic influences on early outcomes after ischemic stroke using a genome-wide approach, revealing seven genome-wide significant loci. A unique aspect of this genetic study is the inclusion of all of the major ethnicities by recruiting from participants throughout the world. Most genetic studies to date have been limited to populations of European ancestry.Implications of all available evidence: The findings provide the first evidence that genes implicating excitotoxicity contribute to human acute ischemic stroke, and demonstrates proof of principle that GWAS of acute ischemic stroke patients can reveal mechanisms involved in ischemic brain injury.

Validation of a clinical-genetics score to predict hemorrhagic transformations after rtPA.

OBJECTIVE: To validate the Genot-PA score, a clinical-genetic logistic regression score that stratifies the thrombolytic therapy safety, in a new cohort of patients with stroke. METHODS: We enrolled 1,482 recombinant tissue plasminogen activator (rtPA)-treated patients with stroke in Spain and Finland from 2003 to 2016. Cohorts were analyzed on the basis of ethnicity and therapy: Spanish patients treated with IV rtPA within 4.5 hours of onset (cohort A and B) or rtPA in combination with mechanical thrombectomy within 6 hours of onset (cohort C) and Finnish participants treated with IV rtPA within 4.5 hours of onset (cohort D). The Genot-PA score was calculated, and hemorrhagic transformation (HT) and parenchymal hematoma (PH) risks were determined for each score stratum. RESULTS: Genot-PA score was tested in 1,324 (cohort A, n = 726; B, n = 334; C, n = 54; and D, n = 210) patients who had enough information to complete the score. Of these, 213 (16.1%) participants developed HT and 85 (6.4%) developed PH. In cohorts A, B, and D, HT occurrence was predicted by the score (p = 2.02 × 10-6, p = 0.023, p = 0.033); PH prediction was associated in cohorts A through C (p = 0.012, p = 0.034, p = 5.32 × 10-4). Increased frequency of PH events from the lowest to the highest risk group was found (cohort A 4%-15.7%, cohort B 1.5%-18.2%, cohort C 0%-100%). The best odds ratio for PH prediction in the highest-risk group was obtained in cohort A (odds ratio 5.16, 95% confidence interval 1.46-18.08, p = 0.009). CONCLUSION: The Genot-PA score predicts HT in patients with stroke treated with IV rtPA. Moreover, in an exploratory study, the score was associated with PH risk in mechanical thrombectomy-treated patients.

PATJ Low Frequency Variants Are Associated With Worse Ischemic Stroke Functional Outcome.

RATIONALE: Ischemic stroke is among the leading causes of adult disability. Part of the variability in functional outcome after stroke has been attributed to genetic factors but no locus has been consistently associated with stroke outcome. OBJECTIVE: Our aim was to identify genetic loci influencing the recovery process using accurate phenotyping to produce the largest GWAS (genome-wide association study) in ischemic stroke recovery to date. METHODS AND RESULTS: A 12-cohort, 2-phase (discovery-replication and joint) meta-analysis of GWAS included anterior-territory and previously independent ischemic stroke cases. Functional outcome was recorded using 3-month modified Rankin Scale. Analyses were adjusted for confounders such as discharge National Institutes of Health Stroke Scale. A gene-based burden test was performed. The discovery phase (n=1225) was followed by open (n=2482) and stringent joint-analyses (n=1791). Those cohorts with modified Rankin Scale recorded at time points other than 3-month or incomplete data on previous functional status were excluded in the stringent analyses. Novel variants in PATJ (Pals1-associated tight junction) gene were associated with worse functional outcome at 3-month after stroke. The top variant was rs76221407 (G allele, ß=0.40, P=1.70×10-9). CONCLUSIONS: Our results identify a set of common variants in PATJ gene associated with 3-month functional outcome at genome-wide significance level. Future studies should examine the role of PATJ in stroke recovery and consider stringent phenotyping to enrich the information captured to unveil additional stroke outcome loci.

Biological Age is a predictor of mortality in Ischemic Stroke.

Age and stroke severity are the main mortality predictors after ischemic stroke. However, chronological age and biological age are not exactly concordant. Age-related changes in DNA methylation in multiple CpG sites across the genome can be used to estimate biological age, which is influenced by lifestyle, environmental factors, and genetic variation. We analyzed the impact of biological age on 3-month mortality in ischemic stroke. We assessed 594 patients with acute ischemic stroke in a cohort from Hospital del Mar (Barcelona) and validated the results in an independent cohort. Demographic and clinical data, including chronological age, vascular risk factors, initial stroke severity (NIHSS score), recanalization treatment, and previous modified Rankin scale were registered. Biological age was estimated with an algorithm based on DNA methylation in 71 CpGs. Biological age was predictive of 3-month mortality (p = 0.041; OR = 1.05, 95% CI 1.00-1.10), independently of NIHSS score, chronological age, TOAST, vascular risk factors, and blood cell composition. Stratified by TOAST classification, biological age was associated with mortality only in large-artery atherosclerosis etiology (p = 0.004; OR = 1.14, 95% CI 1.04-1.25). As estimated by DNA methylation, biological age is an independent predictor of 3-month mortality in ischemic stroke regardless of chronological age, NIHSS, previous modified Rankin scale, and vascular risk factors.

Biological age is better than chronological as predictor of 3-month outcome in ischemic stroke.

OBJECTIVE: To analyze the effect of age-related DNA methylation changes in multiple cytosine-phosphate-guanine (CpG) sites (biological age [b-age]) on patient outcomes at 3 months after an ischemic stroke. METHODS: We included 511 patients with first-ever acute ischemic stroke assessed at Hospital del Mar (Barcelona, Spain) as the discovery cohort. Demographic and clinical data, including chronological age (c-age), vascular risk factors, initial stroke severity, recanalization treatment, and previous and 3-month modified Rankin Scale (p-mRS and 3-mRS, respectively) were registered. B-age was estimated with an algorithm, based on DNA methylation in 71 CpGs. Bivariate analysis determined variables associated with 3-mRS for inclusion in ordinal multivariate analysis. RESULTS: After ordinal regressions for 3-month ischemic stroke outcome (3-mRS), b-age was associated with outcome (odds ratio 1.04 [95% confidence interval 1.01-1.07]), nullifying c-age. Stepwise regression kept b-age, basal NIH Stroke Scale, sex, p-mRS, and recanalization treatment as better explanatory variables, instead of c-age. These results were successfully replicated in an independent cohort. CONCLUSIONS: B-age, estimated by DNA methylation, is an independent predictor of ischemic stroke outcome regardless of chronological years.

Ischemic stroke patients are biologically older than their chronological age.

Ischemic stroke is associated with aging. It is possible to predict chronological age by measuring age-related changes in DNA methylation from multiple CpG sites across the genome, known as biological age. The difference between biological age and actual chronological age would indicate an individual's level of aging. Our aim was to determine the biological age of ischemic stroke patients and compare their aging with controls of the same chronological age. A total of 123 individuals, 41 controls and 82 patients with ischemic stroke were paired by chronological age, ranging from 39 to 82 years. Illumina HumanMethylation450 BeadChip array was used to measure DNA methylation in CpG sites in both groups, and biological age was estimated using methylation values of specific CpGs. Ischemic stroke patients were biologically an average 2.5 years older than healthy controls (p-value=0.010). Stratified by age tertiles, younger stroke patients (≤57 years old) were biologically older than controls (OR=1.19; 95%CI 1.00-1.41, p-value=0.046). The older groups showed no biological age differences between cases and controls, but were close to reaching the significance level. Ischemic stroke patients are biologically older than controls. Biological age should be considered as a potential new biomarker of stroke risk.

Whole exome sequencing analysis reveals TRPV3 as a risk factor for cardioembolic stroke.

Genetic studies suggest that hundreds of genes associated with stroke remain unidentified. Exome sequencing proves useful for finding new genes associated with stroke. We aimed to find new genetic risk factors for cardioembolic stroke by analysing exome sequence data using new strategies. For discovery, we analysed 42 cardioembolic stroke cases and controls with extreme phenotypes (cohort 1), and for replication, 32 cardioembolic stroke cases and controls (cohort 2) using the SeqCapExome capture kit. We then analysed the replicated genes in two new cohorts that comprised 834 cardioembolic strokes and controls (cohort 3) and 64,373 cardioembolic strokes and controls (cohort 4). Transcriptomic in-silico functional analyses were also performed. We found 26 coding regions with a higher frequency of mutations in cardioembolic strokes after correcting for the number of mutations found in the whole exome of every patient. The TRPV3 gene was associated with cardioembolic stroke after replication of exome sequencing analysis (p-value-discovery: 0.018, p-value-replication: 0.014). The analysis of the TRPV3 gene using polymorphisms in cohort 3 and 4 revealed two polymorphisms associated with cardioembolic stroke in both cohorts, the most significant polymorphism being rs151091899 (p-value: 3.1 × 10-05; odds ratio: 5.4) in cohort 3. The genotype of one polymorphism of TRPV3 was associated with a differential expression of genes linked to cardiac malformations. In conclusion, new strategies using exome sequence data have revealed TRPV3 as a new gene associated with cardioembolic stroke. This strategy among others might be useful in finding new genes associated with complex genetic diseases.

Epigenome-wide association study identifies TXNIP gene associated with type 2 diabetes mellitus and sustained hyperglycemia.

Type 2 diabetes mellitus (DM) is an established risk factor for a wide range of vascular diseases, including ischemic stroke (IS). Glycated hemoglobin A1c (HbA1c), a marker for average blood glucose levels over the previous 12 weeks, is used as a measure of glycemic control and also as a diagnostic criterion for diabetes (HbA1c levels ≥ 6.5%). Epigenetic mechanisms, such as DNA methylation, may be associated with aging processes and with modulation of the risk of various pathologies, such as DM. Specifically, DNA methylation could be one of the mechanisms mediating the relation between DM and environmental exposures. Our goal was to identify new CpG methylation sites associated with DM. We performed a genome-wide methylation study in whole-blood DNA from an IS patient cohorts. Illumina HumanMethylation450 BeadChip array was used to measure DNA methylation in CpG sites. All statistical analyses were adjusted for sex, age, hyperlipidemia, body mass index (BMI), smoking habit and cell count. Findings were replicated in two independent cohorts, an IS cohort and a population-based cohort, using the same array. In the discovery phase (N = 355), we identified a CpG site, cg19693031 (located in the TXNIP gene) that was associated with DM (P = 1.17 × 10(-12)); this CpG was replicated in two independent cohorts (N = 167 and N = 645). Methylation of TXNIP was inversely and intensely associated with HbA1c levels (P = 7.3 × 10(-16)), specifically related to diabetic patients with poor control of glucose levels. We identified an association between the TXNIP gene and DM through epigenetic mechanisms, related to sustained hyperglycemia levels (HbA1c ≥ 7%).