Lack of association between ALOX5AP genetic polymorphisms and risk of ischemic stroke: evidence from meta-analyses

Abstract

Background

In recent years, there has been substantial research evaluating the relationship between arachidonate 5-lipoxygenase-activating protein (ALOX5AP) polymorphisms and ischemic stroke (IS). The objective of this study was to systematically review and analyze the existing evidence.

Methods

A comprehensive search of major electronic databases for studies published between 1990 and 2018 was carried out. Data were synthesized as OR and 95% CI using fixed-effects and random-effects models.

Results

A total of 30 studies were available for analysis. The aggregate sample size across all studies was 32,782 (16,294 cases and 16,488 controls). We found no association of the ALOX5AP rs10507391 (OR=1.03 for A allele vs T allele; 95% CI: 0.93–1.14; P=0.557), rs4769874 (OR=1.13 for A allele vs G allele; 95% CI: 1.00–1.28; P=0.050), rs9551963 (OR=1.03 for A allele vs C allele; 95% CI: 0.96–1.11; P=0.372), rs17222814 (OR=1.09 for A allele vs G allele; 95% CI: 0.96–1.24; P=0.195), rs17222919 (OR=0.89 for G allele vs T allele; 95% CI: 0.75–1.06; P=0.175), and rs4073259 (OR=1.20 for A allele vs G allele; 95% CI: 1.00–1.45; P=0.056) polymorphisms with IS risk. Haplotype analysis also did not yield significant findings for the HapA (rs17222814G–rs10507391T–rs4769874G–rs9551963A; OR=1.20; 95% CI: 0.91–1.56; P=0.192) and HapB (rs17216473A–rs10507391A–rs9315050A–rs17222842G; OR=1.11; 95% CI: 0.90–1.38; P=0.339) haplotypes.

Conclusion

Current evidence does not support an association of rs10507391, rs4769874, rs9551963, rs17222814, rs17222919, rs4073259, and HapA and HapB with IS risk.

Keywords:

• ischemic stroke
• ALOX5AP
• genetic polymorphism
• haplotype

Introduction

Stroke is the leading cause of mortality and serious long-term disability in China and worldwide, with ~15 million people throughout the world suffering a stroke every year.¹ Ischemic stroke (IS) accounts for ~85%–90% of all cases and is characterized by the sudden loss of blood circulation to an area of the brain. Strong evidence from twin, family, and animal model studies has consistently suggested a genetic influence on IS risk and prognosis, and in recent years, there has been substantial research evaluating specific genetic risk factors.²

Leukotrienes are short-lived lipid mediators derived from the nuclear membrane of cells that are produced and excreted in response to various immune stimuli. The leukotriene pathway starts with oxidation of arachidonic acid to leukotriene A4 (LTA4) by lipoxygenase-activating protein (ALOX5AP).³ The resultant LTA4 can then be converted to leukotriene B4 (LTB4) by leukotriene A4 hydrolase or it can be conjugated with reduced glutathione by the leukotriene C4.³ As a crucial regulator of biosynthesis of proinflammatory leukotriene lipid mediators, ALOX5AP plays an important role in the development of atherosclerosis. Pharmacological targeting of ALOX5AP significantly reduced atherosclerosis burden and vascular as well as adipose tissue T-cell inflammation.⁴ In addition, knockout of the ALOX5AP gene was associated with ceased leukotriene production and amelioration of stroke damage in a mouse model of middle cerebral artery occlusion.⁵ The findings from clinical trials were consistent with those from animal experiments. Hakonarson et al found that the ALOX5AP inhibitor significantly decreased the concentration of the C-reactive protein, a contributor to the risk of myocardial infarction and IS.⁶ Moreover, ALOX5AP mRNA and protein increased significantly in polymorphonuclear cells from OSA patients vs controls and were associated with carotid luminal diameter and intima–media thickness.⁷

The human ALOX5AP is located on chromosome 13q12-13 and consists of five exons and four introns. Genetic variants in ALOX5AP gene have been studied in stroke cases, but case–control studies performed in different ethnic populations obtain contradictory results, indicating that the existing evidence regarding the association between ALOX5AP polymorphisms and IS risk needs to be systematically reviewed and analyzed. To reconcile inconsistencies across individual studies, we performed a quantitative meta-analysis of the effects of ALOX5AP polymorphisms on IS.

Methods

Search strategy

We searched Embase, Medline, PubMed, and China Knowledge Resource Integrated in July 2018 for studies investigating the relationship between ALOX5AP polymorphisms and IS risk. Searches were not limited by date restrictions and language. Search terms included: (“ALOX5AP polymorphism”, or “ALOX5AP variant”, or “ALOX5AP genotype”), AND (“cerebral infarction” OR “ischemic stroke”). The “AND” operator was used to combine these terms in varying combinations. Once articles had been collected, bibliographies were then hand-searched for additional references. The review followed the Preferred Reporting Items for Systematic Reviews and Meta-Analyses guidelines.⁸

Inclusion and exclusion criteria

Two investigators independently assessed titles and abstracts for relevance and full reports for inclusion. Discrepancies were resolved by discussion. Studies were included in our meta-analysis if they met the following criteria: 1) reported on case–control studies in adult humans; 2) published in peer-reviewed journals; and 3) reported genotypic and/or allelic frequencies. The exclusion criteria were as follows: 1) family-based studies; 2) case-only studies; 3) no information on genotypic and/or allelic frequencies; and 4) editorials, narrative reviews or other manuscripts not reporting primary data. If the paper, or author correspondence, suggested overlapping studies, we included only the most comprehensive study for meta-analysis.

Data extraction and quality assessment

Data were tabled in a standardized Excel sheet, and each group comparison was checked by two investigators to confirm accuracy of inclusion. The following information was abstracted from each study: journal, first author, year of publication, participant characteristics, geographical location, dominant ancestry of sample, diagnostic method for IS, numbers of patients and controls, DNA extraction and genotyping methods, allele frequency, and Hardy–Weinberg equilibrium (HWE) in controls. Authors were not contacted to request missing/additional data. For assessing the methodological quality of the primary studies, a quantitative Newcastle–Ottawa Scale (NOS) score was calculated for each study, which was based on the selection of the study groups, the comparability of the groups, and the ascertainment of the exposure.⁹ Each study was graded as either low (scores 0–5) or high quality (scores 6–9).

Statistical analyses

All statistical analyses were performed using Stata 12.0 software (Stata, College Station, TX, USA). Since the majority of studies reported allele frequency instead of genotype data, the relationship between ALOX5AP polymorphisms and IS risk was assessed under allele contrast. All associations were presented as OR with the corresponding 95% CI. The significance of the pooled OR was determined using a Z-test. Heterogeneity was assessed by Cochrane’s Q test of heterogeneity and the I² statistic. The I² values of 25%, 50%, and 75% corresponded to low, medium, and high between-study heterogeneity. Summary ORs were calculated using random-effects models (DerSimonian and Laird) if there was significant between-study heterogeneity.¹⁰ Random-effects modeling assumes a genuine diversity in the results of various studies and takes into account both within-study and between-study variances. Stratified analyses by sample ancestry were performed in order to assess the potential moderating effect of ancestry. HWE deviation was assessed by using the chi-squared goodness-of-fit test. Sensitivity analysis was conducted by omitting studies with HWE deviation. Funnel plots, Begg’s and Egger’s tests were used to assess publication bias. The standard error of log (OR) of each study was plotted against its log (OR) in funnel plots. The funnel plot symmetry was evaluated by Begg’s and Egger’s tests, which used the correlation between the ranks of effect sizes and the ranks of their variances. A P-value of less than 0.05 was considered statistically significant, except for Cochrane’s Q test, in which a P-value less than 0.10 was applied. No corrections for multiple comparisons were performed.

Results

Study characteristics

A total of 302 studies were identified by the search strategy. After duplicates were removed, 176 studies remained. On examination of the full text, 30 papers were eligible for inclusion in the review.^11–40 Study selection is shown in Figure 1. These 30 studies were reported between 2005 and 2018 and included 32,782 participants. Sample sizes ranged from 200 to 3,727 participants per study and represented Caucasian and Asian populations. Imaging (magnetic resonance imaging/computerized tomography) examination of the brain to support the clinical diagnosis was conducted in all studies. The distribution of genotypes in the controls of all studies was in agreement with HWE except for the studies by Cheng et al, Zhao et al, Xu et al, and Diakite et al.^21,26,29,37 Table 1 displays the study characteristics for the identified studies.

Figure 1 Flow diagram of study selection.

Abbreviation: CNKI, China Knowledge Resource Integrated.

Table 1 Summary of included studies

Author	Year	Sample size		Diagnostic method	Polymorphisms evaluated	Female (%)		Age		Ethnicity	NOS score
		Cases	Controls			Cases	Controls	Cases	Controls
Helgadottir et al^Citation11	2004	702	624	CT/MRI confirmed	HapA and HapB	46.9	NM	NM	20–90	Caucasians	8
Helgadottir et al^Citation12	2005	450	710	CT/MRI confirmed	HapA and HapB	42.0	49.0	66.8±0.6	67.2±0.4	Caucasians	7
Meschia et al^Citation13	2005	377	263	CT/MRI confirmed	rs17222814, rs10507391, rs4769874 and rs9551963	46.4	62.0	64.8±15.0	60.0±14.7	Caucasians	8
Lõhmussaar et al^Citation14	2005	639	736	CT/MRI confirmed	rs17222814, rs10507391, rs4769874, rs9551963 and HapA	36.9	39.3	65.0±18.2	62.0±11.7	Caucasians	7
Zee et al^Citation15	2006	259	259	CT/MRI confirmed	rs17222814, rs10507391, rs4769874, rs9551963, HapA and HapB	0.0	0.0	60.8±0.3	61.0±0.3	Caucasians	7
Zhang et al^Citation16	2006	1,285	1,713	CT/MRI confirmed	rs10507391 and rs4769874	36.7	42.7	60.4±9.2	59.6±8.5	Asians	6
Gao and Qc^Citation17	2008	100	100	CT/MRI confirmed	rs10507391 and rs4769874	35.0	40.0	57.3±6.9	60.0±7.5	Asians	7
Lövkvist et al^Citation18	2008	932	396	CT/MRI confirmed	rs17222814, rs10507391 and rs9551963	44.0	43.2	Age above 65 years: 76.2%	Age above 65 years: 78.8%	Caucasians	7
Quarta et al^Citation19	2009	294	235	CT/MRI confirmed	rs9551963 and rs4769874	40.1	43.8	68–83	60–79	Caucasians	7
He et al^Citation20	2009	412	368	CT/MRI confirmed	rs10507391	42.7	38.0	71.5±7.8	70.3±6.9	Asians	8
Cheng et al^Citation21	2010	280	258	CT/MRI confirmed	rs4073259	37.1	39.5	59.6±12.4	58.5±13.6	Asians
Domingues-Montanari et al^Citation22	2010	1,092	781	CT/MRI confirmed	rs10507391	59.3	54.0	Spain cohort: 70.6±11.9 Portugal cohort: 52.4±9.3	Spain cohort: 72.1±6.9 Portugal cohort: 63.0±6.8	Caucasians	8
Gao et al^Citation23	2010	380	425	CT/MRI confirmed	rs10507391, rs17222814, rs4769874and rs9551963	37.1	48.2	63.4±23.9	56.8±21.8	Asians	7
Kim et al^Citation24	2011	117	398	CT/MRI confirmed	rs17222919	45.3	51.3	65.5±12.1	53.8±15.2	Asians	6
Lee et al^Citation25	2011	291	278	CT/MRI confirmed	rs10507391 and rs4293222	19.9	41.0	65.4±10.6	64.1±11.1	Asians	7
Zhao et al^Citation26	2012	682	598	CT/MRI confirmed	rs17222814, rs10507391, rs4769874 and rs9551963,	40.9	43.8	62.1±9.4	61.8±10.1	Asians	7
Wang et al^Citation27	2012	658	704	CT/MRI confirmed	rs17222919	40.0	44.9	69.4±10.6	69.2±9.0	Asians	6
Sharma et al^Citation28	2013	610	610	CT/MRI confirmed	rs10507391	32.0	32.3	49.3±17.3	49.0±16.8	Asians	8
Xu et al^Citation29	2013	547	794	CT/MRI confirmed	rs9551963 and rs10507391	NM	NM	NM	NM	Asians	7
Zhang et al^Citation30	2014	501	497	CT/MRI confirmed	rs4073259	33.9	39.6	58.7±10.2	53.1±8.6	Asians	8
Yao et al^Citation31	2014	420	488	CT/MRI confirmed	rs10507391, rs4769874 and rs9551963	33.3	35.2	64.7±9.1	64.7±9.3	Asians	7
Yang et al^Citation32	2014	492	490	CT/MRI confirmed	rs10507391, rs4769874, rs9551963 and HapB	41.1	45.3	56.7±8.3	56.2±8.9	Asians	7
Fan et al^Citation33	2015	1,913	1,814	CT/MRI confirmed	rs17222919	46.6	48.4	Population 1: 56.1±10.6 Population 2: 60.5±7.8	Population 1: 55.3±10.3 Population 2: 60.5±7.8	Asians	8
Shao et al^Citation34	2015	459	462	CT/MRI confirmed	rs9551963, rs4769874, rs4769060 and rs10507391	39.9	40.3	48.6±11.0	63.8±9.2	Asians	7
Wang et al^Citation35	2015	396	300	CT/MRI confirmed	rs4769874 and rs10507391	40.9	35.3	56–71	54–75	Asians	7
Qu et al^Citation36	2015	456	452	CT/MRI confirmed	rs10507391, rs4769874, and rs9551963	34.9	33.4	59.7±11.5	53.7±8.2	Asians	7
Diakite et al^Citation37	2016	175	201	CT/MRI confirmed	rs10507391	43.4	NM	57.1±2.0	54.6±2.0	Caucasians	6
Hu et al^Citation38	2017	479	580	CT/MRI confirmed	rs10507391	11.1	12.4	65.2±5.3	65.1±5.7	Asians	7
Liu et al^Citation39	2017	346	399	CT/MRI confirmed	rs4293222	38.2	43.9	63.1±10.7	61.9±10.7	Asians	6
Shi et al^Citation40	2018	550	555	CT/MRI confirmed	rs4073259	49.8	49.9	Population 1: 60.7±12.8 Population 2: 60.5±12.5	Population 1: 60.9±13.3 Population 2: 62.5±13.2	Asians	7

Notes: HapA = SG13S25G–SG13S114T–SG13S89G–SG13S32A. HapB = SG13S377A–SG13S114A–SG13S41A–SG13S35G.

Abbreviations: CT, computed tomography; MRI, magnetic resonance imaging; NM, not mentioned; NOS, Newcastle–Ottawa Scale.

Data synthesis

Twenty studies provided summaries of data on the ALOX5AP rs10507391 polymorphism.^{13–18,20,22,23,25,26,28,29,31,32,34–38} Random-effects meta-analysis indicated no association of the polymorphism with IS in the overall population (OR=1.03 for A allele vs T allele; 95% CI: 0.93–1.14; P=0.557) (Figure 2A; Table 2). There was significant between-study heterogeneity (I²=81.6%, P<0.001). In the subgroup analysis, according to ethnicity, no association between the ALOX5AP rs10507391 polymorphism and IS risk was found in Asian (OR=1.10 for A allele vs T allele; 95% CI: 0.98–1.23; P=0.109) and Caucasian (OR=0.89 for A allele vs T allele; 95% CI: 0.76–1.02; P=0.101) populations (Figure 2A; Table 2).

Figure 2 Forest plots for the association of ALOX5AP polymorphisms with IS risk.

Notes: (A) Meta-analysis of the ALOX5AP rs10507391 polymorphism and risk of IS using the random-effects model (A allele vs T allele). (B) Meta-analysis of the ALOX5AP rs4769874 polymorphism and risk of IS using the fixed-effects model (A allele vs G allele). (C) Meta-analysis of the ALOX5AP rs9551963 polymorphism and risk of IS using the random-effects model (A allele vs C allele). (D) Meta-analysis of the ALOX5AP rs17222814 polymorphism and risk of IS using the fixed-effects model (A allele vs G allele). Weights are from random-effects analysis.
Abbreviation: IS, ischemic stroke.

Table 2 Meta-analysis of association between ALOX5AP polymorphisms and IS risk

Polymorphism	Subgroup	Number of studies	Test of heterogeneity		Test of association		Test of publication
			I^2 (%)	P-value	OR (95% CI)	P-value	Begg’s test	Egger’s test
rs10507391	Total	20	81.6	,0.001	1.03 (0.93–1.14)	0.557	0.974	0.857
	Asians	14	81.2	,0.001	1.10 (0.98–1.23)	0.109	NA	NA
	Caucasians	6	68.7	0.008	0.89 (0.76–1.02)	0.101	NA	NA
	HWE	19	80.8	,0.001	1.05 (0.95–1.16)	0.318	NA	NA
rs4769874	Total	13	28.2	0.160	1.13 (1.00–1.28)	0.050	1.000	0.909
	Asians	9	10.9	0.344	1.10 (0.95–1.28)	0.210	NA	NA
	Caucasians	4	59.2	0.061	1.19 (0.96–1.48)	0.106	NA	NA
rs9551963	Total	12	44.7	0.047	1.03 (0.96–1.11)	0.372	0.451	0.365
	Asians	7	58.6	0.025	1.07 (0.96–1.19)	0.230	NA	NA
	Caucasians	5	0	0.484	0.98 (0.90–1.07)	0.656	NA	NA
	HWE	10	53.3	0.023	1.04 (0.95–1.14)	0.381	NA	NA
rs17222814	Total	6	25.3	0.244	1.09 (0.96–1.24)	0.195	0.256	0.137
	Asians	2	80.4	0.025	1.12 (0.89–1.41)	0.354	NA	NA
	Caucasians	4	0	0.704	1.08 (0.92–1.26)	0.348	NA	NA
rs17222919	Total	4	67.9	0.025	0.89 (0.75–1.06)	0.175	0.308	0.541
	Asians	4	67.9	0.025	0.89 (0.75–1.06)	0.175	0.308	0.541
rs4073259	Total	4	63.7	0.041	1.20 (1.00–1.45)	0.056	0.238	0.156
	Asians	4	63.7	0.041	1.20 (1.00–1.45)	0.056	0.238	0.156
HapA	Total	4	65.6	0.033	1.20 (0.91–1.56)	0.192	0.362	0.621
	Caucasians	4	65.6	0.033	1.20 (0.91–1.56)	0.192	0.362	0.621
HapB	Total	4	0	0.835	1.11 (0.90–1.38)	0.339	0.734	0.251
	Asians	1	NA	NA	1.16 (0.84–1.62)	0.370	NA	NA
	Caucasians	3	0	0.692	1.07 (0.80–1.44)	0.627	NA	NA

Abbreviations: HWE, Hardy–Weinberg equilibrium; IS, ischemic stroke; NA, not applicable.

Thirteen studies contributed to the pooled analysis of the ALOX5AP rs4769874 polymorphism (A allele vs G allele).^{13–17,19,23,26,31,32,34–36} Fixed-effect meta-analysis gave an estimated OR of 1.13 (95% CI: 1.00–1.28; P=0.050), indicating no significant association with IS (Figure 2B; Table 2). There was no significant between-study heterogeneity (I²=10.9%, P=0.160). Neither Asian (OR=1.10, 95% CI 0.95–1.28, P=0.210) nor Caucasian (OR=1.19, 95% CI 0.96–1.48, P=0.106) studies found an association with IS (Figure 2B; Table 2).

Twelve studies examined the ALOX5AP rs9551963 polymorphism.^{13–15,18,19,23,26,29,31,32,34,36} In the primary pooled analysis (Figure 2C; Table 2), we identified no significant association between the polymorphism and IS risk (OR=1.03 for A allele vs C allele; 95% CI: 0.96–1.11; P=0.372). There was significant between-study heterogeneity (I²=44.7%, P=0.047). In the subgroup analysis based on ethnicity, studies with patients of Asian descent or Caucasian descent yielded similar findings (Asian population: OR=1.07, 95% CI 0.96–1.19, P=0.230; Caucasian population: OR=0.98, 95% CI 0.90–1.07, P=0.656), indicating no association with IS (Figure 2C; Table 2).

Six studies provided data for the ALOX5AP rs17222814 polymorphism.^{13–15,18,23,26} Fixed-effect meta-analysis gave an estimated OR of 1.09 (95% CI: 0.96–1.24; P=0.195), indicating no association with IS in the overall population (Figure 2D; Table 2). There was no significant between-study heterogeneity (I²=25.3%, P=0.244). In the subgroup analysis by ethnicity, neither Asian nor Caucasian studies yielded significant findings (Table 2).

Four Asian studies in three publications investigated the role of the ALOX5AP rs17222919 polymorphism.^24,27,33 Random-effects meta-analysis indicated no association of the SNP with IS (OR=0.89 for G allele vs T allele; 95% CI: 0.75–1.06; P=0.175) (Table 2). Similarly, four Asian studies providing summaries of data on the ALOX5AP rs4073259 polymorphism did not find any association with IS risk (OR=1.20 for A allele vs G allele; 95% CI: 1.00–1.45; P=0.056) (Table 2).^21,30,40

In addition to individual polymorphisms, we performed combined analysis for haplotypes. HapA was defined by rs17222814, rs10507391, rs4769874, and rs9551963 polymorphisms, with alleles G, T, G, and A. No association between HapA and IS risk was found (OR=1.20, 95% CI: 0.91–1.56, P=0.192) (Table 2). The pooled effect estimate for HapB (rs17216473A–rs10507391A–rs9315050A–rs17222842G) was similar, indicating no significant association with IS (OR=1.11, 95% CI: 0.90–1.38, P=0.339) (Table 2).

Publication bias and sensitivity analysis

Funnel plots did not suggest publication bias (Figure 3). The results of Begg’s and Egger’s tests for small-study effect biases were not significant, indicating no evidence of publication bias (Table 2). For the ALOX5AP rs10507391 and rs9551963 polymorphisms, sensitivity analyses were performed to assess the stability of the results by omission of studies, which were not in HWE. The significance of pooled ORs was not influenced.

Figure 3 Funnel plots to detect publication bias.

Notes: (A) Funnel plot analysis to detect publication bias for the ALOX5AP rs10507391 polymorphism and risk of IS. (B) Funnel plot analysis to detect publication bias for the ALOX5AP rs4769874 polymorphism and risk of IS. (C) Funnel plot analysis to detect publication bias for the ALOX5AP rs9551963 polymorphism and risk of IS. (D) Funnel plot analysis to detect publication bias for the ALOX5AP rs17222814 polymorphism and risk of IS.
Abbreviations: IS, ischemic stroke; SE, standard error.

Discussion

IS is the leading cause of morbidity and mortality in industrialized countries. Given the limited therapeutic options for IS patients, there remains great interest in identifying novel therapeutic target that may prevent stroke or improve stroke recovery.

The ALOX5AP gene encodes a major regulator for 5-lipoxygenase that catalyzes the oxidation of arachidonic acid to LTA4, which is released by inflammatory cells at injured sites and thus plays an important role in atherosclerosis and other vascular damage.⁴¹ In the literature, there are apparently conflicting associations reported between ALOX5AP polymorphisms and IS risk. The conflicting results were possibly because of small effects of ALOX5AP polymorphisms on IS risk or the relatively low statistical power of published studies. To summarize the available evidence, we combined studies involving a total of 32,782 participants. Six single polymorphisms and two haplotypes were analyzed. Our analysis did not support an association of the ALOX5AP rs10507391, rs4769874, rs9551963, rs17222814, rs17222919, and rs4073259 polymorphisms with IS risk in Asian and Caucasian populations. Haplotype analysis also did not yield significant findings for the HapA (rs17222814G– rs10507391T–rs4769874G–rs9551963A) and HapB (rs17216473A–rs10507391A–rs9315050A–rs17222842G) haplotypes. To the best of our knowledge, this is the largest and most comprehensive meta-analysis on the association between ALOX5AP polymorphisms and IS risk.

Our results are in agreement with a previous meta-analysis published in 2009. Zintzaras et al analyzed the relationship between ALOX5AP polymorphisms and risk of stroke with a total of 5,194 cases and 4,566 controls, finding no association of rs17222814, rs10507391, rs4769874, rs9551963, HapA, and HapB with stroke.⁴² Although the two meta-analyses obtained similar findings, there were still some differences. First, both hemorrhagic and IS patients were included in the Zintzaras et al’s study, whereas we only took into account IS patients. Second, Zintzaras et al’s meta-analysis overtook some Asian studies^16,17,20 because they excluded the non-English papers. Third, Zintzaras et al’s analysis arose from 12 studies, considerably less in number (5,194 cases and 4,566 controls) than the overall 30 in our study (16,294 cases and 16,488 controls). Besides Zintzaras et al’s analysis, the Wang et al’s⁴³ meta-analysis involving 10 studies and the Chen et al’s⁴⁴ meta-analysis combining seven studies investigated rs10507391.^43,44 The positive association between rs10507391 and IS risk reported by these two small size meta-analyses was not confirmed by our study applying a large sample.

It should be mentioned that much effort has been performed to conduct appropriately the meta-analysis and avoid any possible source of bias. To avoid the local literature bias, we included both English and non-English articles. We properly assessed deviations from HWE in the control population and performed sensitivity analysis excluding studies where the sample violated HWE to ensure stability of the pooled ORs. In addition, we evaluated small-study effects and publication bias through Begg’s and Egger’s tests; there was no evidence of systematic missingness of scientific results from the literature. Finally, subgroup analysis by ethnicity confirmed the lack of association between the ALOX5AP polymorphisms and IS risk not only in the overall populations but also in each ethnic group.

Previous studies demonstrated the presence of strong linkage disequilibrium in the ALOX5AP gene.^14,15,32 It is believed that haplotype analysis can provide more information than single-locus analysis because it allows us to explore potential interactions among alleles. In this meta-analysis, we evaluated the relationship of the HapA (rs17222814G– rs10507391T–rs4769874G–rs9551963A) and HapB (rs17216473A–rs10507391A–rs9315050A–rs17222842G) haplotypes with IS risk. Although no associations of these haplotypes with IS were identified, we could not exclude that other haplotypes may play a role in the development of IS. The study by Zhao et al reported a statistically significant association between the haplotype GCGA constructed by rs17222814, rs10507391, rs4769874, and rs9551963 and IS risk.²⁶ In addition, Yao et al found that the haplotype TGC constructed using rs10507391, rs4769874, and rs9551963 was associated with the increased risk of IS (OR=1.60, 95% CI: 1.28–1.98, P<0.001), while the haplotype AGA was protective (OR=0.66, 95% CI: 0.53–0.81) in a case–control study, involving a total of 420 cases and 488 controls.³¹ Considering the large sample sizes applied by these studies, further studies are warranted to confirm their findings.

Our study has some limitations that should be considered. First, the relationship between the ALOX5AP polymorphisms and the levels of leukotriene lipid mediators in IS patients was not evaluated due to lack of published data. Increased levels of products of the 5-lipoxygenase (5-LO)/5-LO-activating protein (FLAP) pathway including LTB4 had been reported in IS patients.⁴⁵ The evaluation of potential interactions between ALOX5AP polymorphisms and leukotriene lipid mediators in IS patients could provide valuable insights into the development of IS. Second, because there was a paucity of data, we did not assess gender differences in the effects of ALOX5AP polymorphisms on IS risk. It is widely known that gender has a complex and interactive effect on IS risk and there are gender differences in the prevalence of IS.⁴⁶ More research is needed to evaluate gender differences in the association of ALOX5AP polymorphisms with IS development.

In conclusion, our meta-analysis of 30 studies involving 16,294 IS patients and 16,488 controls indicates no significant association of the ALOX5AP polymorphisms rs10507391, rs4769874, rs9551963, rs17222814, rs17222919, and rs4073259 with IS risk. In addition, the HapA and HapB haplotypes also did not show any association with IS.

Acknowledgments

This work was supported, in part, by the National Natural Science Foundation of China (81660776 and 81360535), the Natural Science Foundation of Guangxi Zhuang Autonomous Region (2016JJA140115), and the Science Research and Technology Development Project of Guangxi Zhuang Autonomous Region (1598012-55).

Disclosure

The authors report no conflicts of interest in this work.