Effect of gene-gene and gene-environment interaction on the risk of first-ever stroke and poststroke death

Integration of Mendelian randomisation and systems biology models to identify novel blood-based biomarkers for stroke

Stroke is the second largest cause of mortality in the world. Genome-wide association studies (GWAS) have identified some genetic variants associated with stroke risk, but their putative functional causal genes are unknown. Hence, we aimed to identify putative functional causal gene biomarkers of stroke risk. We used a summary-based Mendelian randomisation (SMR) approach to identify the pleiotropic associations of genetically regulated traits (i.e., gene expression and DNA methylation) with stroke risk. Using SMR approach, we integrated cis-expression quantitative loci (cis-eQTLs) and cis-methylation quantitative loci (cis-mQTLs) data with GWAS summary statistics of stroke. We also utilised heterogeneity in dependent instruments (HEIDI) test to distinguish pleiotropy from linkage from the observed associations identified through SMR analysis. Our integrative SMR analyses and HEIDI test revealed 45 candidate biomarker genes (FDR < 0.05; P_HEIDI>0.01) that were pleiotropically or potentially causally associated with stroke risk. Of those candidate biomarker genes, 10 genes (HTRA1, PMF1, FBN2, C9orf84, COL4A1, BAG4, NEK6, SH2B3, SH3PXD2A, ACAD10) were differentially expressed in genome-wide blood transcriptomics data from stroke and healthy individuals (FDR<0.05). Functional enrichment analysis of the identified candidate biomarker genes revealed gene ontologies and pathways involved in stroke, including "cell aging", "metal ion binding" and "oxidative damage". Based on the evidence of genetically regulated expression of genes through SMR and directly measured expression of genes in blood, our integrative analysis suggests ten genes as blood biomarkers of stroke risk. Furthermore, our study provides a better understanding of the influence of DNA methylation on the expression of genes linked to stroke risk.

Large-Scale Targeted Sequencing Study of Ischemic Stroke in the Han Chinese Population

Background

Ischemic stroke is likely caused by interactions of multiple genes and environmental determinants. However, large‐scale sequencing studies to discern functional genetic variants and their interactions with clinical and lifestyle risk factors on ischemic stroke are limited.

Methods and Results

We sequenced functional regions of 740 previously identified genes associated with atherosclerotic disease among 999 ischemic stroke cases and 1001 controls of Chinese ancestry. Multiple logistic regression models were used to examine the associations between variants and ischemic stroke and test interactions between variants and clinical and lifestyle risk factors. Functional variants achieving suggestive significance were replicated in an independent sample of 4724 ischemic stroke cases and 5029 controls. Driven by variant main effects, each minor allele of the correlated rs174535, rs174545, and rs3834458 variants at MYRF‐FADS1‐FADS2 conferred an average 0.83‐fold (95% CI, 0.78–0.88) decreased odds of stroke. Significant main effects of MTHFR rs1801133 missense variant were also observed, with each copy of the A allele associated with a 1.20‐fold (95% CI, 1.13–1.27) higher odds of ischemic stroke. The functional ALDH2 rs671 variant was identified in interaction analyses with alcohol drinking (Meta‐P=3.39×10⁻¹⁷). Each minor allele conferred a 0.54‐fold (95% CI, 0.45–0.64) decreased odds of stroke among drinkers and a 0.89‐fold (95% CI, 0.83–0.97) decreased odds among nondrinkers.

Conclusions

Significant associations at MYRF‐FADS1‐FADS2 indicate that genetically elevated polyunsaturated fatty acids may decrease ischemic stroke risk in East Asians. Significant associations at MTHFR and ALDH2 robustly confirm deleterious effects of genetically elevated homocysteine and alcohol intake, respectively, on ischemic stroke.

New Drug Targets to Prevent Death Due to Stroke: A Review Based on Results of Protein-Protein Interaction Network, Enrichment, and Annotation Analyses

This study used established biomarkers of death from ischemic stroke (IS) versus stroke survival to perform network, enrichment, and annotation analyses. Protein-protein interaction (PPI) network analysis revealed that the backbone of the highly connective network of IS death consisted of IL6, ALB, TNF, SERPINE1, VWF, VCAM1, TGFB1, and SELE. Cluster analysis revealed immune and hemostasis subnetworks, which were strongly interconnected through the major switches ALB and VWF. Enrichment analysis revealed that the PPI immune subnetwork of death due to IS was highly associated with TLR2/4, TNF, JAK-STAT, NOD, IL10, IL13, IL4, and TGF-β1/SMAD pathways. The top biological and molecular functions and pathways enriched in the hemostasis network of death due to IS were platelet degranulation and activation, the intrinsic pathway of fibrin clot formation, the urokinase-type plasminogen activator pathway, post-translational protein phosphorylation, integrin cell-surface interactions, and the proteoglycan-integrin extracellular matrix complex (ECM). Regulation Explorer analysis of transcriptional factors shows: (a) that NFKB1, RELA and SP1 were the major regulating actors of the PPI network; and (b) hsa-mir-26-5p and hsa-16-5p were the major regulating microRNA actors. In conclusion, prevention of death due to IS should consider that current IS treatments may be improved by targeting VWF, the proteoglycan-integrin-ECM complex, TGF-β1/SMAD, NF-κB/RELA and SP1.

Biobank Scale Pharmacogenomics Informs the Genetic Underpinnings of Simvastatin Use

Studying drug-metabolizing enzymes, encoded by pharmacogenes, may inform biological mechanisms underlying the diseases for which a medication is prescribed. Until recently, pharmacogenes could not be studied at biobank scale. In 7,649 unrelated African-ancestry (AFR) and 326,214 unrelated European-ancestry (EUR) participants from the UK Biobank, we associated pharmacogene haplotypes from 50 genes with 265 (EUR) and 17 (AFR) medication use phenotypes using generalized linear models. In EUR, N-acetyltransferase 2 (NAT2) metabolizer phenotype and activity score were associated with simvastatin use. The dose of NAT2*1 was associated with simvastatin use when compared with NAT2*5 (the most common haplotype). This association was robust to effects of low-density lipoprotein cholesterol (LDL-C) concentration (NAT2*1 odds ratio (OR) = 1.07, 95% CI: 1.05-1.09, P = 1.14 × 10^-8 ) and polygenic risk for LDL-C concentration (NAT2*1 OR = 1.09, 95% CI: 1.04-1.14, P = 2.26 × 10^-4 ). Interactive effects between NAT2*1 and simvastatin use on LDL-C concentration (OR = 0.957, 95% CI: 0.916-0.998, P = 0.045) were replicated in the electronic Medical Records and Genomics Pharmacogenetic Sequencing Pilot (eMERGE-PGx) cohort (OR = 0.987, 95% CI: 0.976-0.998, P = 0.029). We used biobank-scale data to uncover and replicate an association between NAT2 locus variation and better response to statin therapy. Testing NAT2 alleles may be useful for making clinical decisions regarding the potential benefit (e.g., absolute risk reduction) in LDL-C concentration prior to statin treatment.

Genetic polymorphisms in pri-let-7a-2 are associated with ischemic stroke risk in a Chinese Han population from Liaoning, China: a case-control study

Ischemic stroke is a complicated disease, and its pathogenesis has been attributed to the occurrence of genetic polymorphisms. Evidence has suggested that the microRNA let-7a is involved in the pathogenesis of ischemic stroke. Pri-miRNA is the primary transcript, which undergoes several processing steps to generate pre-miRNA and, later, mature miRNAs. In this case-control study, we analyzed the distribution of pri-let-7a-2 variants in patients at a high risk for ischemic stroke and the interactions of pri-let-7a-2 variants and environmental factors. Blood samples and clinical information were collected from 1086 patients with ischemic stroke and 836 healthy controls between December 2013 and December 2015 at the First Affiliated Hospital of China Medical University. We found that the rs1143770 CC genotype and the C allele were associated with a decreased risk of ischemic stroke, whereas the rs629367 CC genotype was associated with an increased risk for ischemic stroke. Moreover, these two single-nucleotide polymorphisms were in linkage disequilibrium in this study sample. We analyzed gene-environment interactions and found that rs1143770 exerted a combined effect on the pathogenesis of ischemic stroke, together with alcohol use, smoking, and a history of hypertension. Therefore, the detection of pri-let-7a-2 polymorphisms may increase the awareness of ischemic stroke risk. This study was approved by the Institutional Ethics Committee of the First Affiliated Hospital of China Medical University, China (approval No. 2012-38-1) on February 20, 2012, and was registered with the Chinese Clinical Trial Registry (registration number: ChiCTR-COC-17013559) on December 27, 2017.

Interactive Effects of a Combination of the HDAC3 and HDAC9 Genes with Diabetes Mellitus on the Risk of Ischemic Stroke

BACKGROUND AND AIM

 Previous studies indicated that the HDAC3 and HDAC9 genes play critical roles in atherosclerosis and ischemic stroke (IS). The purpose of this study was to investigate the association of combined single-nucleotide polymorphisms in the HDAC3 and HDAC9 genes with the susceptibility to IS.

METHODS

 A case-control study was conducted including 863 IS patients and 863 age- and gender-matched healthy participants. A polygenic score was developed to estimate the contribution of a combination of the HDAC3 and HDAC9 genes to the risk of IS. The interactive effects of traditional risk factors of stroke and the polygenic score on the risk of IS were explored. Additionally, the association between the polygenic score and the progression of atherosclerosis, a potential risk factor of IS, was examined in our healthy controls.

RESULTS

 Subjects with a higher polygenic score had an increased risk of IS (odds ratio: 1.83; 95% confidence interval: 1.38-2.43) after adjusting for covariates compared with individuals with a lower polygenic score. An interactive effect of diabetes mellitus and the polygenic score on the risk of IS was observed. A significant positive correlation between the polygenic score and a change in the plaque score (standardized β = 0.42, p = 0.0235) in healthy controls with diabetes mellitus was found.

CONCLUSION

 Our results suggested that the combination of the HDAC3 and HDAC9 genes with a history of diabetes mellitus could exacerbate the deterioration of atherosclerosis, thereby increasing the risk of IS. Further studies are warranted to explore our results in other populations.

Temporal Trends and Familial Clustering of Ideal Cardiovascular Health in Parents and Offspring Over the Life Course: An Investigation Using The Framingham Heart Study

Background

Evidence suggests familial aggregation and intergenerational associations for individual cardiovascular health (CVH) metrics. Over a 53‐year life course, we examined trends and association of CVH between parents and their offspring at similar mean ages.

Methods and Results

We conducted a series of cross‐sectional analyses of the FHS (Framingham Heart Study). Parent‐offspring pairs were assessed at exams where their mean age distributions were similar. Ideal CVH was defined using 5 CVH metrics: blood pressure (<120/<80 mm Hg), fasting blood glucose (<100 mg/dL), blood cholesterol (<200 mg/dL), body mass index (<25 kg/m²), and non‐smoking. Joinpoint regression and Chi‐squared test were used to assess linear trend; proportional‐odds regression was used to examine the association between parents and offspring CVH. A total of 2637 parents were paired with 3119 biological offspring throughout 6 exam cycles. Similar patterns of declining ideal CVH with advancing age were observed in parents and offspring. Small proportions of parents (4%) and offspring (17%) achieved 5 CVH metrics at ideal levels (P‐trend <0.001). Offspring of parents with poor CVH had more than twice the odds of having poor CVH (pooled odds ratio, 2.59; 95% CI, 1.98–3.40). Over time, elevated glucose levels and obesity doubled among the offspring and were the main drivers for declining ideal CVH trends.

Conclusions

Parental CVH was positively associated with offspring CVH. However, intergenerational CVH gains from declining smoking rates, cholesterol, and blood pressure were offset by rising offspring obesity and elevated glucose levels. This suggests an intergenerational phenotypic shift of risk factors and the need for a family‐centered approach to cardiovascular care.

Effect of gene-gene and gene-environment interaction on the risk of first-ever stroke and poststroke death

BACKGROUND

Multiple genetic and environmental factors contribute to the individual-level heterogeneity in stroke. This study aimed to assess how the genetic interactions confer risk of stroke.

METHODS

In a Chinese case-control study including 1,405 strokes and 1,263 controls who were followed up (range, 0.1-6.0 years), eight genes, including apolipoprotein(a) (APOA1), methylenetetrahydrofolate reductase (MTHFR), vitamin K epoxide reductase complex subunit 1 (VKORC1), arachidonate 5-lipoxygenase-activating protein (ALOX5AP), NOTCH3, chromosome 9p21.3(Chr.9p21.3), vascular endothelial growth factor (VEGFA), and kinase insert domain-containing receptor (KDR), were analyzed for interactions by the generalized multifactor dimensionality reduction method and validated by the multivariate logistic regression models. The genetic associations with carotid artery intima-media thickness (IMT) were examined.

RESULTS

The interaction of VKORC1 and Chr.9p21.3 was identified for stroke and its worse prognosis, and subjects having the VKORC1 rs2359612C and Chr.9p21.3 rs10757274G alleles had higher risks for stroke (OR = 1.83, 95% CI = 1.32-2.52) as well as for stroke recurrence (HR = 1.84, 95% CI = 1.24-2.73), cardiovascular events (HR = 1.65, 95% CI = 1.15-2.38), and cardiovascular mortality (HR = 2.16, 95% CI = 1.24-3.79). Supporting, they were associated with higher IMT. Hypertension or physical inactivity increased the risk effect. The interaction of VEGFA rs833061C and KDR rs2305948T was identified for hemorrhagic stroke.

CONCLUSIONS

Our findings identified two novel genetic interactions of VKORC1 and Chr.9p21.3 and of VEGFA and KDR for risk of stroke and subtypes as well as future stroke prognosis.