Association between Lys198Asn polymorphism of endothelin‐1 gene and ischemic stroke: A meta‐analysis

Abstract Background Endothelin (ET)‐1 is a potent vasoconstrictor peptide produced by endothelial cells and associated with vascular dysfunction and cardiovascular disease. Lys198Asn is a single‐nucleotide polymorphism (SNP) of gene encoding ET‐1 (EDN1). It is hypothesized that it might have a role in altering ET‐1 and ultimately leading to vascular dysfunction and ischemic stroke. We therefore conducted a meta‐analysis to investigate the association between Lys198Asn polymorphism of EDN1 gene and susceptibility of ischemic stroke. Methods This meta‐analysis was conducted according to the guidance of the Preferred Reporting Items for Systematic Reviews and Meta‐Analyses (PRISMA) statement. We searched PubMed, Google Scholar, Embase, Web of Science, J‐STAGE, and China National Knowledge Infrastructure (CNKI) for relevant studies. The association between Lys198Asn polymorphism and ischemic stroke susceptibility was evaluated by calculating the pooled ORs and 95% CIs. Results Our analysis included 1,291 cases and 2,513 controls. Meta‐analysis established a significant association between Lys198Asn SNP of EDN1 gene and ischemic stroke when assuming either recessive model (OR: 1.30; 95% CI: 1.02–1.65; p = .03; I 2 = 41%) or dominant model (OR: 1.48; 95% CI: 1.24–1.76; p < .001; I 2 = 61%). There was no evidence of publication bias in either of the recessive model (Egger test: p = .23; Begg test: p = .85) or dominant model (Egger test: p = .79; Begg test: p = .85). A subgroup analysis based on subtypes of ischemic stroke showed that Lys198Asn SNP was only associated with large vessel infarction but not with lacunar infarction caused by small vessel disease. A subgroup analysis based on ethnicity revealed that the Lys198Asn polymorphism of the EDN1 gene was associated with ischemic stroke only in Caucasians. Conclusions The present meta‐analysis suggests that Lys198Asn polymorphism of EDN1 gene is associated with an increased risk for ischemic stroke.


| BACKG ROU N D
Stroke is the second leading cause of death globally and leading cause of long-term disability in both developing and developed countries.
Approximately 70% of new strokes are ischemic in origin, 51% stroke death, and 58% of stroke disability-adjusted life years are because of ischemic stroke (Feigin et al., 2016). A complex interaction between modifiable and non-modifiable conventional risk factors and genetic factors could be behind the pathogenesis of ischemic stroke (Boehme, Esenwa, & Elkind, 2017). Genetic influence on the pathogenesis of ischemic stroke has been established by epidemiological and animal studies (Hassan, 2000;Humphries & Morgan, 2004), and influences are polygenic whereby multiple genes exert a small influence or risk on phenotype.
The basis of ischemic stroke is the reduction in cerebral blood flow caused by alteration in structure and function of the blood vessel wall, which begins with changes in the structure and function of the endothelium (Widmer & Lerman, 2014). Endothelin (ET)-1 is a potent vasoconstrictor peptide produced by endothelial cells and has been implicated as an important factor in the development of vascular dysfunction and cardiovascular disease (Böhm & Pernow, 2007). In addition to vasoconstriction, ET-1 is involved in pro-inflammatory actions, stimulation of free radical formation, platelet activation, and proliferative effects (Böhm & Pernow, 2007). Some studies have found high plasma level (Estrada et al., 2017;Sapira, Cojocaru, Lilios, & Grigorian, 2010;Ziv et al., 1992) and high cerebrospinal fluid (CSF) level (Lampl et al., 1997) of ET-1 in ischemic stroke patients. In a study by Dai et al. (2017), injection of ET-1 in rhesus monkey induced transient ischemic attack. Thus ET-1 might have a role in development of ischemic stroke.

| MATERIAL S AND ME THODS
This current meta-analysis was conducted according to the guidance of the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) statement (Liberati et al., 2009).

| Literature search
We searched PubMed, Google Scholar, Embase, Web of Science, J-STAGE, and China National Knowledge Infrastructure (CNKI) for studies published from January 2000 to September 2018. Searches F I G U R E 1 Figure showing posttranslational modification of pre-proendothelin-1 (PPET) protein encoded by endothelin-1 gene (EDN1). Lys198Asn polymorphism does not fall within endothelin (ET)-1 amino acids, but is located near the carboxy-terminal region of PPET, which is removed from PPET by furin-like enzyme during post-translational processing steps were conducted using the keywords "EDN1 polymorphism" or "Endothelin-1 gene polymorphism" or "rs5370 polymorphism" or "Lys198Asn polymorphism" in combination with "stroke" or "ischemic stroke" or "cerebrovascular disease" or "cerebral infarction." Titles and abstracts were screened and the full text of any that were considered to possibly meet the inclusion criteria were then obtained. The cited references of retrieved articles were manually checked to identify any additional eligible trials. Two authors (GN and RO) screened and retrieved reports and excluded irrelevant studies. When uncertainty about eligibility criteria arose, it was solved by discussion with all the authors.

| Eligibility criteria
We included studies that were conducted on human subjects, and the studies were searched without any limitations on language. In this meta-analysis, all included studies met the following criteria: (a) case-control studies focused on the association between rs5370/ Lys198Asn polymorphism and ischemic stroke susceptibility, (b) patients with ischemic stroke were diagnosed with neuroimaging confirmed by neurologist, and (c) there were sufficient data of the genotypes in the case-control groups to evaluate the ORs and 95% CIs. The exclusion criteria were as follows: (a) publications with overlapping cases and controls from the similar study, (b) no genotypic data available, and (c) publications with patients with cardioembolic stroke or hemorrhagic stroke.

| Data abstraction and assessment of methodological quality
The relevant data from each study were independently extracted by two reviewers (GN and RO) using a standardized, structured form. Extracted data included the first author, type of design, site of study, patient ethnicity, year of publication, sample size of cases and controls, control design, Hardy-Weinberg equilibrium, source of DNA, genotyping method, data of the genotypes in the case-control groups, type of ischemic stroke, mean age of participants, mean BMI of participants, and male to female ratio.
The methodological quality of each study was assessed independently by two reviewers (GN and RO) using the parameters previously used by Kumar et al. (2016) In this scale, following sections were included for assessment of methodological quality: representativeness of cases, source of controls, matching of controls, ascertainment of cases, ascertainment of controls, blinding while genotyping, genotyping methods, Hardy-Weinberg equilibrium, and association assessment. Maximum score of 2 can be provided in each section, except source of controls (maximum 3), ascertainment of controls (maximum 1), blinding while genotyping (maximum 1), and association assessment (maximum 1). Methodological quality of scale score varies from 0 to 16, in which higher score represents the better quality and lower score represents the lower quality. Any discrepancies during data extraction and quality assessment were resolved by discussion with all the authors.

| Statistical analysis
In the majority of selected studies, the authors presented odds ratio (OR) and 95% confidence intervals (CIs) for association be-

| Literature search
The results of the systematic literature search and selection are summarized in Figure 2. We identified 364 articles through database searching and three through reference lists searching. After exclusion of duplicates, 258 articles remained, and after screening titles and abstracts, 65 relevant abstracts remained. After application of inclusion and exclusion criteria, six studies were chosen for final analysis and the data were extracted.

| Study and patient characteristics
Our research included studies from two continents and five different countries. Two studies were from China (Li Wei et al., 2009;Zhang & Sui, 2014) and one study was from Japan (Yamaguchi et al., 2006), the United Kingdom (Gormley et al., 2005), Turkey (Aslan et al., 2017), and Ukraine (Dubovyk et al., 2018) each. The publication date ranged from 2005 (Gormley et al., 2005) to 2018 (Dubovyk et al., 2018). All included studies had case-control designs. In all included studies, Hardy-Weinberg equilibrium was achieved in the control population. Venous blood was the source of DNA in all studies. Each study used a unique form of genotyping. Of the six studies, only three studies matched age and sex for controls (Gormley et al., 2005;Li Wei et al., 2009;Zhang & Sui, 2014). Four studies did inpatient recruitment of patient (Aslan et al., 2017;Dubovyk et al., 2018;Li Wei et al., 2009;Zhang & Sui, 2014) while outpatient recruitment (Yamaguchi et al., 2006) and community recruitment (Gormley et al., 2005) was done by one study each. According to quality assessment tools, the methodological quality of most studies was found to be moderately high. Table 1 lists the details of the study characteristics. The patients included in the three studies were Asian (Li Wei et al., 2009;Yamaguchi et al., 2006;Zhang & Sui, 2014) and in the remaining three were Caucasian (Aslan et al., 2017;Dubovyk et al., 2018;Gormley et al., 2005). Two studies recruited small vessel disease patients and four recruited large vessel disease patients. In the included studies, the controls had no history of cerebrovascular, coronary and other atherosclerotic diseases. Among the cases, the mean age ranged between 60 and 71 years. Among the controls, the mean age ranged between 60 and 68 years. All other patient characteristics can be seen in Table 2. Caucasians. Subgroup analysis is tabulated in Table 3.

| Meta-analysis
Meta-regression analysis was used to explore whether the methodological quality scores of the included studies play an important role in the association between the Lys198Asn polymorphism of the EDN1 gene and the risk of ischemic stroke. In recessive model, we found that quality scores were not significant predictors of the effect size for the association between Lys198Asn polymorphism and ischemic stroke risk (p = .14). However, in the dominant model, quality scores were significant predictors of the effect size (p = .02).

Meta-regression plot for recessive model and dominant model is
shown in Figures 7 and 8, respectively.

| D ISCUSS I ON
There is growing evidence that genetic factors are involved in the development of ischemic stroke. This meta-analysis was designed to determine an accurate estimate of the association between the EDN1 Lys198Asn SNP and the risk of ischemic stroke. Several studies have shown an independent association of Lys198Asn SNP of the EDN1 gene and number of diseases such as hypertension (Barath et al., 2007;Barbeau et al., 2003;Fang, Li, & Ma, 2017;Jin et al., 2003 (Colombo, Ciofini, Paradossi, Bevilacqua, & Biagini, 2006), diabetes (Li et al., 2008), pulmonary hypertension (Benza et al., 2015), intraventricular hemorrhage (SzpechtOctober , Gadzinowski, Seremak-Mrozikiewicz, Kurzawińska, & Szymankiewicz, 2016), and subarachnoid hemorrhage (Foreman et al., 2017;Gallek et al., 2013). Therefore, the Lys198Asn SNP is considered to be a risk factor for cardiovascular disease and cerebrovascular disease. Our metaanalysis suggests a significant association of Lys198Asn SNP of EDN1 with ischemic stroke. Our findings are consistent with the results of most of the included studies.
Lys198Asn SNP does not belong to the nucleotide sequence encoding the ET-1 amino acid, but is located near the nucleotide sequence encoding the amino acid adjacent to the carboxy-terminal region of PPET (Figure 1), which is removed from PPET by furin-like enzyme during post-translational processing steps (Tanaka et al., 2005). Considering the location of this SNP, it is unlikely that structural changes will be made in ET-1. However, variation in gene regulation is more important in disease causality than in structural changes in specific proteins. Although Lys198Asn does not structurally alter the ET-1 protein, studies have demonstrated that silent SNP like Lys198Asn can affect human mRNA synthesis, splicing, processing, stability, transport, translation, and protein folding (Robert & Pelletier, 2018;Shen & Basilion, 1999). Furthermore, Lys198Asn SNP can also be coinherited together with another gene highly associated with ischemic stroke due to physical proximity or with a highly functional SNP within the EDN1 gene itself, for the same reasons as above, and can structurally alter ET-1 (Tanaka et al., 2005). In other words, above-mentioned SNP and the causative genetic factor might be in linkage disequilibrium (Tanaka et al., 2005). Such SNP may not have influence on disease susceptibility however they can be an important maker if they coinherit with other highly disease-associated gene or with other highly functional SNP of same gene.  In our subgroup analysis, we found that Lys198Asn SNP was only associated with large vessel infarction but not with lacunar infarction caused by small vessel disease (SVD). This may be due to a number of issues, including sample size, heterogeneous pathophysiology, and silent SVD. Silent SVD is defined as hyperintensities of the white matter that correlate with an asymptomatic cerebral infarction. Many cases of lacunar infarctions due to SVD are silent. As these cases often do not present to the hospital, they are excluded from the analysis. Complex interactions between genetic and nongenetic factors have been reported to play an important role in the progression of ischemic stroke (Boehme et al., 2017). Ethnicity is one of the important nongenetic factors. Therefore, we conducted an ethnicity-based subgroup analysis. In the recessive model, the Lys198Asn SNP of the EDN1 gene was not associated with an ischemic stroke in both Asians and Caucasians. However, in dominant model, the Lys198Asn SNP of the EDN1 gene was associated with ischemic stroke in Caucasians but not in Asians. In Asians, especially in Chinese, lacunar infarctions due to SVD are believed to be responsible for a higher proportion of ischemic stroke than in Caucasians (Tsai, Thomas, & Sudlow, 2013). Since the Lys198Asn polymorphism is not associated with lacunar infarction due to SVD, this could be the reason why the Lys198Asn SNP of the EDN1 gene is not associated with an ischemic stroke in Asians. Yamaguchi et al. (2006) included female patients only for analysis of Lys198Asn SNP and it was significantly associated with ischemic stroke. Dubovyk et al. (2018) and Zhang and Sui (2014)  provided inconsistent results due to low statistical power. Among the included studies, different methodologies were used for genotyping and the selection of controls. The haplotype analysis, which plays a pivotal role in genetic association studies, was not performed. The gene-environment interaction and the gene-gene interaction were not studied in this meta-analysis. Owing to the above limitations, the results of this meta-analysis should be interpreted with caution.

| CON CLUS ION
Our meta-analysis suggests a significant association of Lys198Asn polymorphism of EDN1 gene with ischemic stroke. Lys198Asn polymorphism might influence on disease susceptibility and can also be an important marker of risk for ischemic stroke. Large cohort studies with gene-gene and gene-environment interactions in diverse ethnicity should be conducted to reach more reliable conclusion.

ACK N OWLED G M ENTS
None.

CO N FLI C T O F I NTE R E S T
The authors declare no conflicts of interest.

AUTH O R S ' CO NTR I B UTI O N S
GN and BKY designed the study. GN and RO carried out literature search, review, and selection. GN and HD carried out the statistical analysis and drafted the manuscript. All authors were involved in revising the manuscript critically for important intellectual content. All authors read and approve the final manuscript.

E TH I C A L A PPROVA L
This research does not include human participants.

DATA AVA I L A B I L I T Y S TAT E M E N T
Data sharing is not applicable to this article as no new data were created or analyzed in this study.

I N FO R M E D CO N S E NT
For this type of study formal consent is not required.