A multitrait genetic study of hemostatic factors and hemorrhagic transformation after stroke treatment

BACKGROUND

Thrombolytic recombinant tissue plasminogen activator (r-tPA) treatment is the only pharmacologic intervention available in the ischemic stroke acute phase. This treatment is associated with an increased risk of intracerebral hemorrhages, known as hemorrhagic transformations (HTs), which worsen the patient's prognosis.

OBJECTIVES

To investigate the association between genetically determined natural hemostatic factors' levels and increased risk of HT after r-tPA treatment.

METHODS

Using data from genome-wide association studies on the risk of HT after r-tPA treatment and data on 7 hemostatic factors (factor [F]VII, FVIII, von Willebrand factor [VWF], FXI, fibrinogen, plasminogen activator inhibitor-1, and tissue plasminogen activator), we performed local and global genetic correlation estimation multitrait analyses and colocalization and 2-sample Mendelian randomization analyses between hemostatic factors and HT.

RESULTS

Local correlations identified a genomic region on chromosome 16 with shared covariance: fibrinogen-HT, P = 2.45 × 10-11. Multitrait analysis between fibrinogen-HT revealed 3 loci that simultaneously regulate circulating levels of fibrinogen and risk of HT: rs56026866 (PLXND1), P = 8.80 × 10-10; rs1421067 (CHD9), P = 1.81 × 10-14; and rs34780449, near ROBO1 gene, P = 1.64 × 10-8. Multitrait analysis between VWF-HT showed a novel common association regulating VWF and risk of HT after r-tPA at rs10942300 (ZNF366), P = 1.81 × 10-14. Mendelian randomization analysis did not find significant causal associations, although a nominal association was observed for FXI-HT (inverse-variance weighted estimate [SE], 0.07 [-0.29 to 0.00]; odds ratio, 0.87; 95% CI, 0.75-1.00; raw P = .05).

CONCLUSION

We identified 4 shared loci between hemostatic factors and HT after r-tPA treatment, suggesting common regulatory mechanisms between fibrinogen and VWF levels and HT. Further research to determine a possible mediating effect of fibrinogen on HT risk is needed.

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.

Causal role of lipid metabolome on the risk of ischemic stroke, its etiological subtypes, and long-term outcome: A Mendelian randomization study

BACKGROUND AND AIMS

The lipid profile is consistently associated with coronary artery disease (CAD) and ischemic stroke (IS). However, the lipoprotein subfractions have not been deeply explored in stroke subtypes, especially in IS outcome.

METHODS

We performed two-sample Mendelian randomization (MR) analysis using 92 lipid traits measured by nuclear magnetic resonance in 115,000 subjects from the UK Biobank. Data for genetic associations with IS, its subtypes, and long-term outcome (LTO) were obtained from three cohorts of European ancestry: GIGASTROKE (73,652 cases, 1,234,808 controls), GODS (n = 1791) and GISCOME (n = 6165). Results obtained using CARDIoGRAMPlusC4D were used to identify differences with CAD.

RESULTS

Genetically determined low concentration of medium high-density lipoprotein (HDL) particles (odds ratio (OR) = 0.92, 95% CI 0.88-0.96; p = 3.6 × 10-4) and its cholesterol content (OR = 0.92, 95% CI 0.88-0.96; p = 1.9 × 10-4) showed causal associations with an increased risk of stroke. Genetic predisposition to high apolipoprotein (apo)B to apoA-I ratio was causally associated with an increased risk of IS (OR = 1.12, 95% CI 1.06-1.18, p = 1.1 × 10-4), and a highly suggestive association was found between non-esterified cholesterol in low-density lipoprotein (LDL) and increased risk of atherothrombotic stroke (LAS) (OR = 1.35, 95% CI 1.10-1.66; p = 4.0 × 10-3). Low cholesterol in small and medium LDL was suggestively associated with poor LTO.

CONCLUSIONS

Our results support that low medium HDL concentration was causally associated with an increased stroke risk, while high levels of non-esterified cholesterol in LDL were suggestively associated with an increased risk of LAS and with a better LTO.

Ischaemic stroke patients present sex differences in gut microbiota

BACKGROUND

Gut microbiota plays a role in the pathophysiology of ischaemic stroke (IS) through the bidirectional gut-brain axis. Nevertheless, little is known about sex-specific microbiota signatures in IS occurrence.

METHODS

A total of 89 IS patients and 12 healthy controls were enrolled. We studied the taxonomic differences of the gut microbiota between men and women with IS by shotgun metagenomic sequencing. To evaluate the causal effect of several bacteria on IS risk, we performed a two-sample Mendelian randomisation (MR) with inverse-variance weighting (IVW) using genome-wide association analysis (GWAS) summary statistics from two cohorts of 5959 subjects with genetic and microbiota data and 1,296,908 subjects with genetic and IS data, respectively.

RESULTS

α-Diversity analysis measured using Observed Species (p = 0.017), Chao1 (p = 0.009) and Abundance-based Coverage Estimator (p = 0.012) indexes revealed that IS men have a higher species richness compared with IS women. Moreover, we found sex-differences in IS patients in relation to the phylum Fusobacteria, class Fusobacteriia, order Fusobacteriales and family Fusobacteriaceae (all Bonferroni-corrected p < 0.001). MR confirmed that increased Fusobacteriaceae levels in the gut are causally associated with an increased risk of IS (IVW p = 0.02, β = 0.32).

CONCLUSIONS

Our study is the first to indicate that there are gut microbiome differences between men and women with IS, identifying high levels of Fusobacteriaceae in women as a specific risk factor for IS. Incorporating sex stratification analysis is important in the design, analysis and interpretation of studies on stroke and the gut microbiota.

Genetic Architecture of Ischaemic Strokes after COVID-19 Shows Similarities with Large Vessel Strokes

We aimed to analyse whether patients with ischaemic stroke (IS) occurring within eight days after the onset of COVID-19 (IS-COV) are associated with a specific aetiology of IS. We used SUPERGNOVA to identify genome regions that correlate between the IS-COV cohort (73 IS-COV cases vs. 701 population controls) and different aetiological subtypes. Polygenic risk scores (PRSs) for each subtype were generated and tested in the IS-COV cohort using PRSice-2 and PLINK to find genetic associations. Both analyses used the IS-COV cohort and GWAS from MEGASTROKE (67,162 stroke patients vs. 454,450 population controls), GIGASTROKE (110,182 vs. 1,503,898), and the NINDS Stroke Genetics Network (16,851 vs. 32,473). Three genomic regions were associated (p-value < 0.05) with large artery atherosclerosis (LAA) and cardioembolic stroke (CES). We found four loci targeting the genes PITX2 (rs10033464, IS-COV beta = 0.04, p-value = 2.3 × 10-2, se = 0.02), previously associated with CES, HS6ST1 (rs4662630, IS-COV beta = -0.04, p-value = 1.3 × 10-3, se = 0.01), TMEM132E (rs12941838 IS-COV beta = 0.05, p-value = 3.6 × 10-4, se = 0.01), and RFFL (rs797989 IS-COV beta = 0.03, p-value = 1.0 × 10-2, se = 0.01). A statistically significant PRS was observed for LAA. Our results suggest that IS-COV cases are genetically similar to LAA and CES subtypes. Larger cohorts are needed to assess if the genetic factors in IS-COV cases are shared with the general population or specific to viral infection.

A Polygenic Risk Score Based on a Cardioembolic Stroke Multitrait Analysis Improves a Clinical Prediction Model for This Stroke Subtype

Background

Occult atrial fibrillation (AF) is one of the major causes of embolic stroke of undetermined source (ESUS). Knowing the underlying etiology of an ESUS will reduce stroke recurrence and/or unnecessary use of anticoagulants. Understanding cardioembolic strokes (CES), whose main cause is AF, will provide tools to select patients who would benefit from anticoagulants among those with ESUS or AF. We aimed to discover novel loci associated with CES and create a polygenetic risk score (PRS) for a more efficient CES risk stratification.

Methods

Multitrait analysis of GWAS (MTAG) was performed with MEGASTROKE-CES cohort (n = 362,661) and AF cohort (n = 1,030,836). We considered significant variants and replicated those variants with MTAG p-value < 5 × 10⁻⁸ influencing both traits (GWAS-pairwise) with a p-value < 0.05 in the original GWAS and in an independent cohort (n = 9,105). The PRS was created with PRSice-2 and evaluated in the independent cohort.

Results

We found and replicated eleven loci associated with CES. Eight were novel loci. Seven of them had been previously associated with AF, namely, CAV1, ESR2, GORAB, IGF1R, NEURL1, WIPF1, and ZEB2. KIAA1755 locus had never been associated with CES/AF, leading its index variant to a missense change (R1045W). The PRS generated has been significantly associated with CES improving discrimination and patient reclassification of a model with age, sex, and hypertension.

Conclusion

The loci found significantly associated with CES in the MTAG, together with the creation of a PRS that improves the predictive clinical models of CES, might help guide future clinical trials of anticoagulant therapy in patients with ESUS or AF.

Biological Age Acceleration Is Lower in Women With Ischemic Stroke Compared to Men

BACKGROUND

Stroke onset in women occurs later in life compared with men. The underlying mechanisms of these differences have not been established. Epigenetic clocks, based on DNA methylation (DNAm) profiles, are the most accurate biological age estimate. Epigenetic age acceleration (EAA) measures indicate whether an individual is biologically younger or older than expected. Our aim was to analyze whether sexual dichotomy at age of stroke onset is conditioned by EAA.

METHODS

We used 2 DNAm datasets from whole blood samples of case-control genetic studies of ischemic stroke (IS), a discovery cohort of 374 IS patients (N women=163, N men=211), from GRECOS (Genotyping Recurrence Risk of Stroke) and SEDMAN (Dabigatran Study in the Early Phase of Stroke, New Neuroimaging Markers and Biomarkers) studies and a replication cohort of 981 IS patients (N women=411, N men=570) from BASICMAR register. We compared chronological age, 2 DNAm-based biomarkers of aging and intrinsic and extrinsic epigenetic age acceleration EAA (IEAA and extrinsic EAA, respectively), in IS as well as in individual IS etiologic subtypes. Horvath and Hannum epigenetic clocks were used to assess the aging rate. A proteomic study using the SOMAScan multiplex assay was performed on 26 samples analyzing 1305 proteins.

RESULTS

Women present lower Hannum-extrinsic EAA values, whereas men have higher Hannum-extrinsic EAA values (women=-0.64, men=1.24, P=1.34×10^-2); the same tendency was observed in the second cohort (women=-0.57, men=0.79, P=0.02). These differences seemed to be specific to cardioembolic and undetermined stroke subtypes. Additionally, 42 blood protein levels were associated with Hannum-extrinsic EAA (P<0.05), belonging to the immune effector process (P=1.54×10^-6) and platelet degranulation (P<8.74×10^-6) pathways.

CONCLUSIONS

This study shows that sex-specific underlying biological mechanisms associated with stroke onset could be due to differences in biological age acceleration between men and women.

Pharmacogenetics studies in stroke patients treated with rtPA: a review of the most interesting findings

Recombinant tissue-plasminogen activator (rtPA) is the only drug used during the acute phase of stroke. Despite its important benefits, a percentage of patients suffer symptomatic hemorrhagic transformations or a lack of early recanalization rates. These undesirable effects are associated with acute neurological and long-term functional deterioration. For the past 20 years, pharmacogenetic studies have tried to find the genetic risk factors associated with rtPA response. Most of these studies have used a gene-candidate strategy; however, recent genome-wide association studies have emerged indicating that genetic predisposition could modulate rtPA response. This review summarizes the most interesting findings in this field, including which genes and genetic variations are associated with hemorrhagic transformations and recanalization rates after thrombolytic therapy.

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.

RP11-362K2.2:RP11-767I20.1 Genetic Variation Is Associated with Post-Reperfusion Therapy Parenchymal Hematoma. A GWAS Meta-Analysis

Stroke is one of the most common causes of death and disability. Reperfusion therapies are the only treatment available during the acute phase of stroke. Due to recent clinical trials, these therapies may increase their frequency of use by extending the time-window administration, which may lead to an increase in complications such as hemorrhagic transformation, with parenchymal hematoma (PH) being the more severe subtype, associated with higher mortality and disability rates. Our aim was to find genetic risk factors associated with PH, as that could provide molecular targets/pathways for their prevention/treatment and study its genetic correlations to find traits sharing genetic background. We performed a GWAS and meta-analysis, following standard quality controls and association analysis (fastGWAS), adjusting age, NIHSS, and principal components. FUMA was used to annotate, prioritize, visualize, and interpret the meta-analysis results. The total number of patients in the meta-analysis was 2034 (216 cases and 1818 controls). We found rs79770152 having a genome-wide significant association (beta 0.09, p-value 3.90 × 10⁻⁸) located in the RP11-362K2.2:RP11-767I20.1 gene and a suggestive variant (rs13297983: beta 0.07, p-value 6.10 × 10⁻⁸) located in PCSK5 associated with PH occurrence. The genetic correlation showed a shared genetic background of PH with Alzheimer’s disease and white matter hyperintensities. In addition, genes containing the ten most significant associations have been related to aggregated amyloid-β, tau protein, white matter microstructure, inflammation, and matrix metalloproteinases.

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

[Figure: see text].

The copy number variation and stroke (CaNVAS) risk and outcome study

BACKGROUND AND PURPOSE

The role of copy number variation (CNV) variation in stroke susceptibility and outcome has yet to be explored. The Copy Number Variation and Stroke (CaNVAS) Risk and Outcome study addresses this knowledge gap.

METHODS

Over 24,500 well-phenotyped IS cases, including IS subtypes, and over 43,500 controls have been identified, all with readily available genotyping on GWAS and exome arrays, with case measures of stroke outcome. To evaluate CNV-associated stroke risk and stroke outcome it is planned to: 1) perform Risk Discovery using several analytic approaches to identify CNVs that are associated with the risk of IS and its subtypes, across the age-, sex- and ethnicity-spectrums; 2) perform Risk Replication and Extension to determine whether the identified stroke-associated CNVs replicate in other ethnically diverse datasets and use biomarker data (e.g. methylation, proteomic, RNA, miRNA, etc.) to evaluate how the identified CNVs exert their effects on stroke risk, and lastly; 3) perform outcome-based Replication and Extension analyses of recent findings demonstrating an inverse relationship between CNV burden and stroke outcome at 3 months (mRS), and then determine the key CNV drivers responsible for these associations using existing biomarker data.

RESULTS

The results of an initial CNV evaluation of 50 samples from each participating dataset are presented demonstrating that the existing GWAS and exome chip data are excellent for the planned CNV analyses. Further, some samples will require additional considerations for analysis, however such samples can readily be identified, as demonstrated by a sample demonstrating clonal mosaicism.

CONCLUSION

The CaNVAS study will cost-effectively leverage the numerous advantages of using existing case-control data sets, exploring the relationships between CNV and IS and its subtypes, and outcome at 3 months, in both men and women, in those of African and European-Caucasian descent, this, across the entire adult-age spectrum.

Genome-wide transcriptome study in skin biopsies reveals an association of E2F4 with cadasil and cognitive impairment

CADASIL is a small vessel disease caused by mutations in NOTCH3 that lead to an odd number of cysteines in the EGF-like repeat domain, causing protein misfolding and aggregation. The main symptoms are migraine, psychiatric disturbances, recurrent strokes and dementia, being executive function characteristically impaired. The molecular pathways altered by this receptor aggregation need to be studied further. A genome-wide transcriptome study (four cases paired with three healthy siblings) was carried out, in addition to a qRT-PCR for validation purposes (ten new cases and eight new controls). To study the expression profile by cell type of the significant mRNAs found, we performed an in situ hybridization (ISH) (nine cases and eight controls) and a research in the Single-nuclei Brain RNA-seq expression browser (SNBREB). Pathway analysis enrichment was carried out with Gene Ontology and Reactome. Neuropsychological tests were performed in five of the qRT-PCR cases. The two most significant differentially expressed mRNAs (BANP, p-value = 7.23 × 10-4 and PDCD6IP, p-value = 8.36 × 10-4) were selected for the validation study by qRT-PCR. Additionally, we selected two more mRNAs (CAMK2G, p-value = 4.52 × 10-3 and E2F4, p-value = 4.77 × 10-3) due to their association with ischemic neuronal death. E2F4 showed differential expression in the genome-wide transcriptome study and in the qRT-PCR (p = 1.23 × 10-3), and it was upregulated in CADASIL cases. Furthermore, higher E2F4 expression was associated with worse executive function (p = 2.04 × 10-2) and attention and information processing speed (IPS) (p = 8.73 × 10-2). In situ hibridization showed E2F4 expression in endothelial and vascular smooth vessel cells. In silico studies indicated that E2F4 is also expressed in brain endothelial cells. Among the most significant pathways analyzed, there was an enrichment of vascular development, cell adhesion and vesicular machinery terms and autophagy process. E2F4 is more highly expressed in the skin biopsy of CADASIL patients compared to controls, and its expression is present in endothelial cells and VSMCs. Further studies are needed to understand whether E2F4 could be useful as a biomarker, to monitor the disease or be used as a therapeutic target.

A genome-wide methylation study reveals X chromosome and childhood trauma methylation alterations associated with borderline personality disorder

Borderline personality disorder (BPD) is a severe and highly prevalent psychiatric disorder, more common in females than in males and with notable differences in presentation between genders. Recent studies have shown that epigenetic modifications such as DNA methylation may modulate gene × environment interactions and impact on neurodevelopment. We conducted an epigenome wide study (Illumina Infinium HumanMethylation450k beadchip) in a group of BPD patients with (N = 49) and without (N = 47) childhood traumas and in a control group (N = 44). Results were confirmed in a replication cohort (N = 293 BPD patients and N = 114 controls) using EpiTYPER assays. Differentially methylated CpG sites were observed in several genes and intragenic regions in the X chromosome (PQBP1, ZNF41, RPL10, cg07810091 and cg24395855) and in chromosome 6 (TAP2). BPD patients showed significantly lower methylation levels in these CpG sites than healthy controls. These differences seemed to be increased by the existence of childhood trauma. Comparisons between BPD patients with childhood trauma and patients and controls without revealed significant differences in four genes (POU5F1, GGT6, TNFRSF13C and FAM113B), none of them in the X chromosome. Gene set enrichment analyses revealed that epigenetic alterations were more frequently found in genes controlling oestrogen regulation, neurogenesis and cell differentiation. These results suggest that epigenetic alterations in the X chromosome and oestrogen-regulation genes may contribute to the development of BPD and explain the differences in presentation between genders. Furthermore, childhood trauma events may modulate the magnitude of the epigenetic alterations contributing to BPD.

A parsimonious score with a free web tool for predicting disability after an ischemic stroke: the Parsifal Score

BACKGROUND

Most of the models to predict prognosis after an ischemic stroke include complex mathematical equations or too many variables, making them difficult to use in the daily clinic. We want to predict disability 3 months after an ischemic stroke in an independent patient not receiving recanalization treatment within the first 24 h, using a minimum set of variables and an easy tool to facilitate its implementation. As a secondary aim, we calculated the capacity of the score to predict an excellent/devastating outcome and mortality.

METHODS

Eight hundred and forty-four patients were evaluated. A multivariable ordinal logistic regression was used to obtain the score. The Modified Rankin Scale (mRS) was used to estimate disability at the third month. The results were replicated in another independent cohort (378 patients). The "polr" function of R was used to perform the regression, stratifying the sample into seven groups with different cutoffs (from mRS 0 to 6).

RESULTS

The Parsifal score was generated with: age, previous mRS, initial NIHSS, glycemia on admission, and dyslipidemia. This score predicts disability with an accuracy of 80-76% (discovery-replication cohorts). It has an AUC of 0.86 in the discovery and replication cohort. The specificity was 90-80% (discovery-replication cohorts); while, the sensitivity was 64-74% (discovery-replication cohorts). The prediction of an excellent or devastating outcome, as well as mortality, obtained good discrimination with AUC > 0.80.

CONCLUSIONS

The Parsifal Score is a model that predicts disability at the third month, with only five variables, with good discrimination and calibration, and being replicated in an independent cohort.

DNA methylation of MMPs and TIMPs in atherothrombosis process in carotid plaques and blood tissues

Background and Purpose: Polymorphisms and serum levels of Matrix Metalloproteinases (MMP) and Tissue Inhibitor of Metalloproteinases (TIMP) have been studied with regard to atheromatous plaques and ischemic stroke, while no studies of DNA methylation (DNAm) patterns of MMP or TIMP have been performed to that end. Here, we evaluate DNAm levels of the MMP and TIMP gene families in human carotid plaques and blood samples of atherothrombotic stroke patients.

Methods: We profiled the DNAm status of stable and ulcerated atherosclerotic plaques obtained as pair sets from three patients who underwent carotid endarterectomy surgery. We selected 415 CpG sites, mapping into MMPs and TIMPs genes for further study. Secondly, the statistically associated CpG sites were analyzed in blood samples from two separate atherothrombotic stroke cohorts (total sample size = 307), ischemic stroke-cohort 1 (ISC-1): 37 atherothrombotic patients and 6 controls, ischemic stroke-cohort 2 (ISC-2): 80 atherothrombotic patients and 184 controls. DNAm levels from plaque tissue and blood samples were evaluated using a high-density microarray Infinium, HumanMethylation450 BeadChip and Infinium MethylationEPIC BeadChip.

Results: Three CpG sites were statistically significantly associated with unstable plaque portions; cg02969624, q-value = 0.035 (TIMP2), and cg04316754, q-value = 0.037 (MMP24) were hypermethylated, while cg24211657 q-value = 0.035 (TIMP2) was hypomethylated. Association of cg04316754 (MMP24) methylation levels with atherothrombotic risk was also observed in blood tissue: ISC-1 p-values = 0.03, ISC-2 p-value = 1.9 × 10^-04.

Conclusions: The results suggest different DNAm status of MMP24 between stable and unstable atherothrombotic carotid plaques, and between atherothrombotic stroke and controls in blood samples.

Association study of polymorphisms within inflammatory genes and methylation status in treatment response in major depression

Background:

Although pharmacogenetics for major depressive disorder (MDD) is gaining momentum, the role of genetics in differences in response to antidepressant treatment is controversial, as they depend on multifactorial and polygenic phenotypes. Previous studies focused on the genes of the serotonergic system, leaving apart other pathological factors such as the inflammatory pathway. The main objective of the study was to assess whether treatment response might be associated with specific inflammation-related genetic variants or their methylation status.

Methods:

41 SNPs in 8 inflammatory genes: interleukin (IL) 1-β, IL2, IL6, IL6R, IL10, IL18, tumor necrosis factor (TNF)-α and interferon (IFN)-γ were genotyped in 153 patients with MDD, who were evaluated with the Mausdley Staging Method to determine treatment response profiles. Pyrosequencing reactions and methylation quantification were performed in a PyroMark Q24 in 5 selected CpG islands of IL1- β, IL6 and IL6R. Linear and logistic regression analyses were conducted, including age and gender as covariates using PLINK 1.07.

Results:

Allelic distribution of IL1- β rs1143643 was significantly associated with MSM scores (FDR corrected p = 0.04). Allelic distribution of IL6R rs57569414 showed a trend towards significance with MSM scores (p = 0.002; FDR corrected p = 0.07). Haplotype analyses showed associations between allelic combinations of IL1-β and IL10 with treatment response (FDR corrected p < 0.01). Methylation percentage of treatment responders was only higher in an IL6R CpG island (p < 0.05).

Conclusions:

These exploratory findings suggest that IL1-β and, marginally, IL6R polymorphisms may affect treatment response in major depression. If confirmed, these results may account for the heterogeneous phenotypes of major depression that underlie differences in treatment response.

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.

Genome-Wide Association Study of White Blood Cell Counts in Patients With Ischemic Stroke

Background and Purpose- Immune cells play a key role in the first 24h poststroke (acute phase), being associated with stroke outcome. We aimed to find genetic risk factors associated with leukocyte counts during the acute phase of stroke. Methods- Ischemic stroke patients with leukocyte counts data during the first 24h were included. Genome-wide association study and gene expression studies were performed. Results- Our genome-wide association study, which included 2064 (Discovery) and 407 (Replication) patients, revealed a new locus (14q24.3) associated with leukocyte counts. After Joint analysis (n=2471) 5 more polymorphisms reached genome-wide significance (P<5×10^-8). The 14q24.3 locus was associated with acute stroke outcome (rs112809786, P=0.036) and with ACOT1 and PTGR2 gene expression. Previous polymorphisms associated with leukocyte counts in general-population did not show any significance in our study. Conclusions- We have found the first locus associated with leukocyte counts in ischemic stroke, also associated with acute outcome. Genetic analysis of acute endophenotypes could be useful to find the genetic factors associated with stroke outcome. Our findings suggested a different modulation of immune cells in stroke compared with healthy conditions.

Clinical Variables and Genetic Risk Factors Associated with the Acute Outcome of Ischemic Stroke: A Systematic Review

Stroke is a complex disease and one of the main causes of morbidity and mortality among the adult population. A huge variety of factors is known to influence patient outcome, including demographic variables, comorbidities or genetics. In this review, we expound what is known about the influence of clinical variables and related genetic risk factors on ischemic stroke outcome, focusing on acute and subacute outcome (within 24 to 48 hours after stroke and until day 10, respectively), as they are the first indicators of stroke damage. We searched the PubMed data base for articles that investigated the interaction between clinical variables or genetic factors and acute or subacute stroke outcome. A total of 61 studies were finally included in this review. Regarding the data collected, the variables consistently associated with acute stroke outcome are: glucose levels, blood pressure, presence of atrial fibrillation, prior statin treatment, stroke severity, type of acute treatment performed, severe neurological complications, leukocyte levels, and genetic risk factors. Further research and international efforts are required in this field, which should include genome-wide association studies.

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

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.

Systematic Review of Cysteine-Sparing NOTCH3 Missense Mutations in Patients with Clinical Suspicion of CADASIL

CADASIL (cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy) is caused by mutations in the NOTCH3 gene, affecting the number of cysteines in the extracellular domain of the receptor, causing protein misfolding and receptor aggregation. The pathogenic role of cysteine-sparing NOTCH3 missense mutations in patients with typical clinical CADASIL syndrome is unknown. The aim of this article is to describe these mutations to clarify if any could be potentially pathogenic. Articles on cysteine-sparing NOTCH3 missense mutations in patients with clinical suspicion of CADASIL were reviewed. Mutations were considered potentially pathogenic if patients had: (a) typical clinical CADASIL syndrome; (b) diffuse white matter hyperintensities; (c) the 33 NOTCH3 exons analyzed; (d) mutations that were not polymorphisms; and (e) Granular osmiophilic material (GOM) deposits in the skin biopsy. Twenty-five different mutations were listed. Four fulfill the above criteria: p.R61W; p.R75P; p.D80G; and p.R213K. Patients carrying these mutations had typical clinical CADASIL syndrome and diffuse white matter hyperintensities, mostly without anterior temporal pole involvement. Cysteine-sparing NOTCH3 missense mutations are associated with typical clinical CADASIL syndrome and typical magnetic resonance imaging (MRI) findings, although with less involvement of the anterior temporal lobe. Hence, these mutations should be further studied to confirm their pathological role in CADASIL.

GRECOS Project (Genotyping Recurrence Risk of Stroke): The Use of Genetics to Predict the Vascular Recurrence After Stroke

BACKGROUND AND PURPOSE

Vascular recurrence occurs in 11% of patients during the first year after ischemic stroke (IS) or transient ischemic attack. Clinical scores do not predict the whole vascular recurrence risk; therefore, we aimed to find genetic variants associated with recurrence that might improve the clinical predictive models in IS.

METHODS

We analyzed 256 polymorphisms from 115 candidate genes in 3 patient cohorts comprising 4482 IS or transient ischemic attack patients. The discovery cohort was prospectively recruited and included 1494 patients, 6.2% of them developed a new IS during the first year of follow-up. Replication analysis was performed in 2988 patients using SNPlex or HumanOmni1-Quad technology. We generated a predictive model using Cox regression (GRECOS score [Genotyping Reurrence Risk of Stroke]) and generated risk groups using a classification tree method.

RESULTS

The analyses revealed that rs1800801 in the MGP gene (hazard ratio, 1.33; P=9×10^-03), a gene related to artery calcification, was associated with new IS during the first year of follow-up. This polymorphism was replicated in a Spanish cohort (n=1.305); however, it was not significantly associated in a North American cohort (n=1.683). The GRECOS score predicted new IS (P=3.2×10^-09) and could classify patients, from low risk of stroke recurrence (1.9%) to high risk (12.6%). Moreover, the addition of genetic risk factors to the GRECOS score improves the prediction compared with previous Stroke Prognosis Instrument-II score (P=0.03).

CONCLUSIONS

The use of genetics could be useful to estimate vascular recurrence risk after IS. Genetic variability in the MGP gene was associated with vascular recurrence in the Spanish population.

PPM1A Methylation Is Associated With Vascular Recurrence in Aspirin-Treated Patients

BACKGROUND AND PURPOSE

Despite great efforts by pharmacogenetic studies, the causes of aspirin failure to prevent the recurrence of ischemic events remain unclear. Our aim was to study whether epigenetics could be associated with the risk of vascular recurrence in aspirin-treated stroke patients.

METHODS

We performed an epigenetic joint analysis study in 327 patients treated with aspirin. In the discovery stage, we performed a nested case-control study in 38 matched ischemic stroke patients in whom 450 000 methylation sites were analyzed. Nineteen patients presented vascular recurrence after stroke, and 19 matched patients did not present vascular recurrence during the first year of follow-up. In a second stage, 289 new patients were analyzed by EpiTYPER.

RESULTS

The following 3 differentially methylated candidate CpG sites, were identified in the discovery stage and analyzed in the second stage: cg26039762 (P=9.69×10(-06), RAF1), cg04985020 (P=3.47×10(-03), PPM1A), and cg08419850 (P=3.47×10(-03), KCNQ1). Joint analysis identified an epigenome-wide association for cg04985020 (PPM1A; P=1.78×10(-07)), with vascular recurrence in patients treated with aspirin.

CONCLUSIONS

The pattern of differential methylation in PPM1A is associated with vascular recurrence in aspirin-treated stroke patients.

TRAF3 Epigenetic Regulation Is Associated With Vascular Recurrence in Patients With Ischemic Stroke

Background and Purpose—

Clopidogrel is one of the most used antiplatelet drugs in patients with cardiovascular disease. However, 16% to 50% of patients have a high on-clopidogrel platelet reactivity and an increased risk of ischemic events. The pathogenesis of high on-treatment platelet reactivity in patients with stroke is only partially explained by genetic variations. This study aims to find differentially methylated sites across the genome associated with vascular recurrence in ischemic stroke patients treated with clopidogrel.

Methods—

From a cohort of 1900 patients with ischemic stroke, we selected 42 patients treated with clopidogrel, including 21 with a recurrent vascular event and 21 without vascular recurrence during the first year of follow-up. Over 480 000 DNA methylation sites were analyzed across the genome. Differentially methylated CpG sites were identified by nonparametric testing using R. Replication analysis was performed in a new cohort of 191 subjects and results were correlated with platelet reactivity in a subset of 90 subjects using light transmission aggregometry.

Results—

A total of 73 differentially methylated CpG sites ( P <1×10 −05 ) were identified; 3 of them were selected for further replication: cg03548645 ( P =1.42×10 −05 , TRAF3 ), cg09533145 ( P =7.81×10 −06 , ADAMTS2 ), and cg15107336 ( P =1.89×10 −05 , XRCC1 ). The cg03548645 CpG remained significant in the replication study ( P =0.034), a deep analysis of this region revealed another methylation site associated with vascular recurrence, P =0.037. Lower cg03548645 ( TRAF3 ) DNA methylation levels were correlated with an increased platelet aggregation (ρ=−0.29, P =0.0075).

Conclusions—

This study suggests for the first time that epigenetics may significantly contribute to the variability of clopidogrel response and recurrence of ischemic events in patients with stroke.