Abstract

Introduction. Rising studies indicate that the apolipoprotein E (APOE) gene is related to the susceptibility of ischemic stroke (IS). However, certain consensus is limited by the lack of a large sample size of researches. This meta-analysis was performed to explore the potential association between the APOE gene and IS. Methods. To identify relevant case control studies in English publications by October 2020, we searched PubMed, Embase, Web of Science, and the Cochrane Library. Pooled odds ratios (ORs) with fixed- or random-effect models and corresponding 95% confidence intervals (CIs) were calculated to analyze potential associations. Results. A total of 55 researches from 32 countries containing 12207 IS cases and 27742 controls were included. The association between APOE gene ε4 mutation and IS was confirmed (ε4 vs. ε3 allele: pooled , 95% CI, 1.214-1.556; ε2/ε4 vs. ε3/ε3: pooled , 95% CI, 1.056-1.440; ε3/ε4 vs. ε3/ε3: pooled , 95% CI, 1.165-1.542; ε4/ε4 vs. ε3/ε3: pooled , 95% CI, 1.542-2.179; and APOE ε4 carriers vs. non-ε4 carriers: pooled ; 95% CI, 1.203-1.576). Interestingly, APOE ε4 mutation showed a dose-response correlation with IS risk (ε4/ε4 vs. ε2/ε4: pooled ; 95% CI, 1.281-2.060; ε4/ε4 vs. ε3/ε4: pooled ; 95% CI, 1.077-1.571). Similar conclusions were drawn in the small artery disease (SAD) subtype, but not in large artery atherosclerosis (LAA) or in cardioaortic embolism (CE), by subgroup analysis. Conclusions. These observations reveal that specific APOE ε4 mutation was significantly associated with the risk of IS in a dose-dependent manner, while APOE ε4 mutation was related to SAD subtype onset without a cumulative effect.

1. Introduction

Ischemic stroke (IS) is a disturbing problem worldwide, which is attributable to its leading role in disability and mortality worldwide, regardless of age, ethnicity, or gender [1]. Uncovering the etiology of IS is crucial for recognition and prevention of this disorder. Genetic elements and environmental components positively contribute to this multifactorial disease [2, 3]. Genetic inheritance provides a guide to the identification of high-risk individual. It deserves to investigate candidate gene polymorphisms in IS pathophysiological pathways. The apolipoprotein E (APOE) gene locates on chromosome 19q13.2. Two single polymorphisms (rs7412 and rs729358), three common alleles (ε2, ε3, and ε4), and six genotypes (ε2/ε2, ε2/ε3, ε2/ε4, ε3/ε3, ε3/ε4, and ε4/ε4) generate in populations [4]. The product of the APOE gene is a polymorphic protein named apolipoprotein E, which modulates the translocation of the cholesterol and other lipids among highly diverse cells [5], involved with neuroinflammation [6] and myelin integrity maintenance [7]. A study indicated that the activated CypA–MMP9 pathway in APOE4 carriers facilitated pericyte injury, which caused blood vessel dysfunction [8]. APOE polymorphisms and its risk associations with coronary artery disease [9], hypertension [10], diabetes [11], and carotid arterial atherosclerosis [12] are widely debated. The abovementioned diseases place individuals at a potential serious risk of IS. Individual studies of the association between IS and APOE polymorphisms have been explored extensively. Clinical differences, ethnic diversities, and small sample sizes restricted the present finding to an inconsistent and controversial one. Previous meta-analyses concerning to this issue have been published several years ago [13] or limited to specific ethnicity [14, 15]. Accordingly, researches from 32 countries are qualified to form our meta-analysis to clarify how APOE genotypes are associated with IS. Moreover, we firstly revealed the correlation of the APOE gene and three IS subtypes (large artery atherosclerosis (LAA), small artery disease (SAD), and cardioaortic embolism (CE)).

2. Materials and Methods

We followed the rules of the preferred reporting items for systematic reviews and meta-analyses (PRISMA) statement to make this meta-analysis [16].

2.1. Data Availability

The data that contribute to the findings in our study are available and the corresponding authors can be contacted for data access.

2.2. Literature Search

Online databases (PubMed, Embase, Web of Science, and the Cochrane Library) were comprehensively searched for studies potentially involved and published in English publications and prior to October 30, 2020. We used a combination of some search terms relevant for IS (stroke, cerebral infarct, brain infarct, ischemic stroke, cerebral ischemia, transient ischemic attack, and cerebrovascular accident) and for the APOE gene (apolipoprotein E, APOE polymorphisms, apolipoprotein E polymorphisms, apolipoprotein E gene, rs429358, rs7412, apolipoprotein E epsilon 4, APOE e4, apolipoprotein E epsilon 2, and APOE e2). The detailed search strategies were showed next.

2.3. Selection Criteria

The selection of the studies was independently completed by two investigators, and any difference was resolved by discussion until an agreement was reached. We carefully selected case control studies that evaluated the relationship of the APOE gene and IS with definite IS diagnoses (using computed tomography, magnetic resonance, or autopsy) regardless of the ethnic background. The detailed inclusion criteria were (1) high-quality studies which explore the relationship between the APOE gene and IS, (2) explicit IS diagnostic criteria, (3) nonstroke individuals as the control group, and (4) original data including independent and sufficient APOE genotype data, to compute ORs and 95% CIs. The newest and largest studies were chosen to avoid duplicate or overlapped data information.

2.4. Data Extraction

Two investigators separately finished full-text reading to extract the needed information from each selected study and resolved the controversial items through serious discussion. The extracted information was (1) research characteristics, including the first author’s name, year of publication, and geographical location of the study; (2) participant details, such as the sex ratio, mean age, and the sample size of case and control groups; (3) diagnostic criteria for IS; (4) determination methods of the APOE gene; (5) each genotype frequency; (6) the sample sizes of IS subtypes according to TOAST norms and respective genotype frequency; and (7) HWE in controls.

2.5. Quality Assessment

We performed the quality assessment through the Newcastle-Ottawa Scale (NOS) score considering selection, comparability, and exposure. It ranged from 0 (worst) to 9 (best) and high-quality studies were known as with a NOS .

2.6. Statistical Analysis

We performed Stata 14.0 to complete all data analyses. The chi-square test was used to examine the Hardy-Weinberg equilibrium (HWE) in control groups. An overt deviation from HWE was regarded as . The compositive ORs and 95% CI were calculated. We explored five genetic models to generate the respective pooled ORs: (1) allele comparisons (ε2 allele vs. ε3 allele; ε4 allele vs. ε3 allele); (2) genotype comparisons (ε2/ε2 vs. ε3/ε3; ε2/ε3 vs. ε3/ε3; ε2/ε4 vs. ε3/ε3; ε3/ε4 vs. ε3/ε3; ε4/ε4 vs. ε3/ε3); (3) APOE ε4 carrier comparisons: we defined three ε4-containing genotypes () as APOE ε4 carriers and the other genotypes () as non-APOE ε4 carriers; (4) APOE ε2 carrier comparisons: similar comparisons of ε2-containing genotypes () vs. non-ε2-containing genotypes (); and (5) comparisons between APOE ε4 homozygosis and ε4 heterozygote (ε4/ε4 vs. ε2/ε4; ε4/ε4 vs. ε3/ε4). The statistic and Cochran’s test were applied to measure the heterogeneity between studies [17]. We selected the random effect model (DerSimonian-Laird method) when heterogeneity was found between studies ( and fixed-effect model (Mantel-Haenszel method) when no heterogeneity existed (). Subgroup analysis was conducted to confirm the relationship between the APOE polymorphisms and the risk of different IS subgroups. Sensitivity analysis was performed by successively removing a single study one by one to verify the stability and reliability of our conclusions. Meta-regression analysis was operated to recognize sources of heterogeneity. Funnel plots and quantified Egger’s tests were accomplished to test publication bias. Significant publication bias was considered as the value of Egger’s test less than 0.10 or obvious asymmetric funnel plot.

2.7. The Result of Trial Sequential Analysis (TSA)

Insufficient sample size, continuous updating, and repeating “ significance testing” could increase the risk of type I errors. Therefore, traditional meta-analysis that focuses on the specific topic may suffer an increased risk of random error. Trial sequential analysis (TSA) was used to reduce the risk of type I error and obtain important information regarding the required sample size for such trials. Set the time sequence of a single study as the research node, and then, perform an interim analysis between the new study that will be included in meta-analysis and existing data accumulation. The required information size (RIS), trial sequential monitoring boundary, and futility boundary are estimated using the TSA. As the sample size of meta-analysis reaching the RIS or the -curve crossing the trial sequential monitoring boundary, we can conclude that the results of meta-analysis are quite stable and further studies were not needed. We accomplished TSA following the guidelines of the user manual and previous article [18] by setting a significance of 5% for type I error, a relative risk reduction of 20%, and a statistical test power of 80% with TSA software (TSA, version 0.9 beta; Copenhagen Trial Unit, Copenhagen, Denmark).

3. Results

3.1. Characteristics of Eligible Studies

We collect a total of 55 studies from 32 countries containing 12207 IS cases and 27742 controls to make the meta-analysis [1973]. Figure 1 showed the detailed selection process. The selected studies and their main characteristics were exhibited in Table 1. Fifteen of the studies provided data about different subtypes (grouped by classification of cerebrovascular diseases III or TOAST classification) of IS: large artery atherosclerosis (LAA), small artery disease (SAD), and cardioaortic embolism (CE). We extracted them independently and specific information was showed in supplementary material table 1. There were seven studies (Koopal et al. 2016, Lai et al. 2007, Chowdhury et al. 2001, Kokubo et al. 2000, Ji et al. 1998, Couderc et al. 1993, Saidi et al. 2009) which deviated HWE obviously, and one study (Schneider et al. 2005) did not contain enough data to obtain HWE. Forty-eight studies used PCR-based method and seven researches (Slowik et al. 2003, Karttunen et al. 2002, Hachinski et al. 1996, Couderc et al. 1993, Brewin et al. 2020, Aalto-Setala et al. 1998, Schneider et al. 2005) used other methods to identify APOE genotypes. These studies used computed tomography or magnetic resonance to diagnose IS except that one research which used autopsy (Schneider et al. 2005). The NOS score mean value was 7.509, which suggested that the quality of included studies was reliable (supplementary material Table 2). PRISMA2020 checklist was provided to present our meta-analysis items (supplementary material Table 3).

3.2. Main Results of the Comparisons in the Abovementioned Five Genetic Models
3.2.1. Allele Comparisons

In comparison with the ε3 allele, the ε2 allele did not show association of the risk of IS (pooled , 95% CI, 0.867-1.115, ) (as showed in Table 2), while the ε4 allele contributed to an obviously increased risk of IS with the pooled (95% CI, 1.214-1.556, ) (Figure 2(d)).

3.2.2. Genotype Comparisons

When compared with the ε3/ε3 genotype, the pooled effects of the APOE genotype in the meta-analysis were as follows: for the ε2/ε2 genotype, pooled , 95% CI, 0.653-1.486, , and for the ε2/ε3 genotype, pooled , 95% CI, 0.900-1.066, ; those two genotypes presented no association with the risk of IS (as showed in Table 2). Genotypes ε2/ε4, ε3/ε4, and ε4/ε4 were related to a higher risk of IS than ε3/ε3. The respective IS risk ORs were 1.233 (95% CI, 1.056-1.440, ) (Figure 2(a)), 1.340 (95% CI, 1.165-1.542, ) (Figure 2(b)), and 1.833 (95% CI, 1.542-2.179, ) (Figure 2(c)). The above results could be found in Table 2. A conclusion was drawn: every genotype which contained APOE ε4 mutation increased the risk of IS.

3.2.3. APOE ε4 Carrier Comparisons

Compared with the non-ε4 carriers, we confirmed that the ε4 carriers were associated with the increased risk of IS; the pooled outcome was pooled (95% CI, 1.203-1.576, ) (Figure 2(e)).

3.2.4. APOE ε2 Carrier Comparisons

In the genetic model of ε2 carriers vs. non-ε2 carriers, there was no association with the IS risk (pooled , 95% CI 0.841-1.086, ) (Table 2).

3.2.5. APOE ε4 Homozygosis versus APOE ε4 Heterozygote Comparisons

Given the above, the APOE ε4 mutation was linked to IS risk. To identify whether there is a dose-response relationship between the ε4 allele and IS or not, we implemented the comparisons between the ε4/ε4 genotype and ε4 heterozygotes (ε2/ε4 or ε3/ε4 genotype). Compared with the ε2/ε4 and ε3/ε4 genotypes, the IS risk ORs for ε4/ε4 genotypes were 1.625 (95% CI, 1.281-2.060, ) and 1.301 (95% CI, 1.077-1.571, ), respectively (Figures 2(f) and 2(g)); this part provided evidence that ε4 homozygosis might generate a higher risk of IS than ε4 heterozygotes.

3.3. Main Results of the Relationship between APOE Gene and Three IS Subtypes

We further investigated on the correlation of APOE gene polymorphisms and risks of IS subtypes by making comparisons in five genetic models, with a particular focus on the APOE ε4 mutation. Subgroup analyses showed that APOE ε4 mutation significantly increased SAD risk (ε4 allele vs. ε3 allele: pooled , 95% CI, 1.073-1.618, (Figure 3(d)); ε3/ε4 vs.ε3/ε3: pooled , 95% CI, 1.097-1.767, (Figure 3(b)); ε4/ε4 vs. ε3/ε3: pooled , 95%, CI 1.030-3.175, (Figure 3(c)); and APOE ε4 carriers vs. non-APOE ε4 carriers: pooled , 95% CI, 1.064-1.661, (Figure 3(e))). But genotype ε2/ε4 did not increase the risk of SAD onset (Figure 3(a)). The result of APOE ε4 homozygosis versus ε4 heterozygote comparisons (ε4/ε4 vs. ε2/ε4 and ε4/ε4 vs. ε3/ε4) was a matter of concern: APOE ε4 mutation could not cause a cumulative effect in generating higher risk of SAD onset, as showed in Figures 3(f) and 3(g).

3.4. Sensitivity Analysis

Sensitivity analysis was performed by removing studies one by one to check the effect of the individual study on overall ORs. No single study influenced on the pooled ORs and 95% CIs in all genetic model comparisons as our data showed (supplementary material table 4).

3.5. Publication Bias

We carried out publication bias analysis by using funnel plots as qualitative description and Egger’s regression tests as quantitative outcome. Funnel plots of all genetic model comparisons did not exhibit apparent asymmetry (several funnel plots were showed in supplementary material figure 1 and 2). In addition to subtype analysis of ε2/ε2 vs. ε3/3, all the Egger’s regression test outcomes indicated that there existed no evident publication bias with all values exceeding 0.1 (supplementary material table 5). The above results showed that publication bias of our meta-analysis was not significant.

3.6. Regression Analysis

Meta-regression analysis was then performed to explore sources of heterogeneity as shown in supplementary material table 5, considering the year of publication, region, sample size, genotyping method, HWE, NOS score, and source of control. However, the value of each factor affecting overall heterogeneity was not statistically significant in comparisons of ε3/ε4 vs. ε3/3, ε4 vs. non-ε4, ε2 vs. non-ε2, ε4allele vs. ε3allele, and ε2allele vs. ε3allele (supplementary material figure 3). Heterogeneity sources were unascertainable.

3.7. The Result of Trial Sequential Analysis (TSA)

The RIS was 8901 samples and the sample size of our meta-analysis reached it. Moreover, the cumulative -curve crossed the trial sequential monitoring boundary before reaching the RIS as showed in Figure 4. The result of TSA guaranteed the stability of our meta-analysis results. Our sample size was proved to be enough for evaluating the relationship between APOE polymorphisms and IS risk.

4. Discussion

Recently, scholars explored more how gene polymorphisms were contributing to the occurrence and prognosis of diseases. And several previous publications had well explored how gene polymorphisms related to diseases onset and potential mechanisms [74, 75]. As a heterogeneous multifactorial disorder, ischemic stroke could be regulated by certain gene synthesis and specific gene products. The genes involved in the pathological process of stroke are also worth of attention. Apolipoprotein E has been proven to affect atherosclerosis, neurodegeneration, and the process of nerve damage repair. That is why we explored the relationship between APOE gene polymorphisms and ischemic stroke risk.

APOE is a 299-amino acid protein encoded by the APOE gene of three common polymorphisms, ε2, ε3, and ε4. The correlation of APOE gene polymorphisms and the risk of cerebral vascular and degenerative diseases have been investigated a lot, especially in Alzheimer’s disease (AD) and cerebral amyloid angiopathy (CAA) [76]. APOE ε4 is associated with increased risk for AD whereas APOE ε2 is associated with decreased risk [77]. Mirza et al. performed a meta-analysis to find that greater WMH volume was associated with worse performance on all cognitive domains in APOE ε4 carriers only in AD [78]. Charidimou et al. proved that the APOE ε2 allele might be associated with the pathophysiology and severity of cortical superficial siderosis in CAA [79]. As to IS, there existed quite many researches with inconsistent conclusions. Besides method differences, ethnic difference and unclarified pathophysiological mechanisms are probable reasons of the inconsistency.

In a meta-analysis in 1999, McCarron et al. found that the ε4 allele and carriers were more frequent among patients with ischemic cerebrovascular disease, compared with control subjects (27% versus 18%; odds ratio, 1.73; 95% CI, 1.34-2.23; ) [13]. In another meta-analysis based on Chinese population, the ε4 allele is associated with an increased risk of developing cerebral infarction, in which the adjusted risk estimate for the ε4 allele versus ε3 allele was significant (, 95% CI 1.59-2.53, ) [14]. Our estimates seemed to be coinciding with the above ones. Compared with the ε3 allele, the ε4 allele showed a higher risk of IS. Compared with ε3/ε3, both ε4 heterozygote (ε2/ε4, ε3/ε4) and ε4 homozygosis (ε4/ε4) exhibited a significant correlation with an increased risk of IS. Notably, OR in ε4 homozygosis (ε4/ε4 vs. ε3/3: 1.833 (95% CI 1.542-2.179)) was higher than those in ε4 heterozygotes (ε2/ε4 vs. ε3/3: 1.233 (95% CI 1.056-1.440) and ε3/ε4 vs. ε3/3: 1.340 (95% CI 1.165-1.542)), which implied that the ε4 allele might possess a cumulative effect. Then, we performed comparisons between ε4/ε4 and ε2/ε4 or ε3/ε4; there existed significant differences between ε4 homozygosis and ε4 heterozygote. The OR between ε4/ε4 and ε2/ε4 was 1.625 (95% CI 1.281-2.060, ); the OR between ε4/ε4 and ε3/ε4 was 1.301 (95% CI 1.077-1.571, ), giving a hint that ε4 homozygosis might bring a higher risk of IS than ε4 heterozygotes.

There are tremendous researches and discussions focusing on the pathogenicity of ε4. An Indian research reported that VLDL and triglycerides levels were found to be significantly associated with ε2/ε4 and ε3/ε4 genotypes; the ε4 allele exerted a higher influence than the ε3 allele in plasma cholesterol levels [22]. As a lipid transport protein, APOE3 and APOE2 preferentially bind to the smaller, more phospholipid-enriched high-density lipoproteins (HDL), while APOE4 preferentially binds to the larger, triglyceride-rich very low-density lipoproteins (VLDL). Miyata and Smith demonstrated an antioxidant activity in the order , and other researchers also reported similar results that APOE4 was associated with increased oxidative stress [25, 80], which might play a role in atherosclerosis and lead to increased risk of ischemic vascular diseases. Besides the above reasons, APOE4 was proved to be neurotoxic by assuming an abnormal conformation (the unique domain interaction between Arg-61 and Glu-255) which was highly susceptible to neuron specific proteolysis and generating neurotoxic fragments that escaped the secretory pathway and entered the cytosol [81]. Totally, from pathophysiological mechanisms to clinical research results, it seems that APOE4 is indeed related to a higher risk of IS, compared with other isoforms, both in ε4 heterozygote and homozygous. ε2 allele appears to be unclear and controversial in stroke [13]. In a meta-analysis of Martínez-González et al., compared with ε3/ε3, APOE ε2 was associated with intracerebral hemorrhage (; 95% CI, 1.01-1.74); meanwhile, APOE ε2 was more related to lobar hemorrhage than deep hemorrhage [82]. As to the association of IS with APOE based on previous investigation, it is uncertain. Our estimates showed that both ε2/ε2 and ε2/ε3 genotypes exhibited no significant effects on IS risk, compared with ε3/ε3. Also, no differences were found in comparisons of ε2 allele vs. ε3 allele and ε2 vs. non-ε2 carriers. This result remained consistent with another meta-analysis in 2013 [14]. Interestingly, in subtype analysis, ε2/ε2 displayed significances in the CE group (; 95% CI, 1.917-9.600; ) and SAD group (; 95% CI, 1.037-3.134; ). The largest meta-analysis of the APOE genotype with IS showed a positive linear association of increasing risk when ordered from ε2/ε2, ε2/ε3, ε2/ε4, ε3/ε3, ε3/ε4, and ε4/ε4 in European ancestry population [83]. The conclusion might explain why APOE4 brings a higher risk of IS but could not clarify that the CE and SAD subgroups in comparison of ε2/ε2 with ε3/ε3 show significances. It is well known that all patients with type III hyperlipidemia (dysbetalipoproteinemia) were APOE ε2 homozygous, whereas most ε2/ε2 subjects (>90%) were normolipidemic or even hypolipidemic, owing to reductions in LDL or HDL or both. Therefore, the APOE ε2 allele has both increased and decreased risks for atherosclerosis, which induced a comprehensive and undetermined result [84].

As to our subtype analyses, all LAA groups showed no significant difference among comparisons, which raised a question why isoforms of APOE, a lipid transport protein, seemed not to be related with IS caused by large artery atherosclerosis. Besides lipid metabolism and atherosclerosis, there might exist some other pathways underlying the relationships between APOE and risk of IS. Our estimates displayed that APOE isoforms were associated to risk of IS especially in the SAD subgroup. Hypertension was known to be an independent risk factor of SAD. Atherosclerosis, dyslipidemia, and hypertension have a complex interaction, and the causations with APOE need further investigation.

Our meta-analysis has several limitations. First, just as the abovementioned, heterogeneity between studies remains undeterminable. Second, results of our meta-analysis based on case control studies cannot provide a causal relationship, but only an association. Third, age variable and ethnicity can influence APOE frequencies in a population; we cannot obtain sufficient related information to perform further subdivided subgroup analyses. Fourth, other pathogenic factors about IS, a multifactorial disease, such as plasma lipid levels, hypertension, life-style, BMI, and gene-environment interactions, were unachievable. Fifth, the controls in accessible studies were not strictly defined; some were selected from healthy populations and others were from nonstroke people. The expected genotype distribution in controls was not in accordance with HWE in seven studies. Population selection in control groups failed to avoid certain diseases which might have a relation with the APOE gene, such as dyslipidemia, hypertension, other vascular diseases, and diabetes. Sixth, the case groups were not selected by a prospective process and the design of case control studies often caused abnormal gene frequency.

5. Conclusions

In conclusion, our meta-analysis provides rational evidence that APOE ε4 mutation is a genetic risk factor for IS. Prospective studies of a large sample size, which concerns gene-gene and gene-environment interactions, should be carried out in the future to reach a more comprehensive outcome about the association of APOE gene polymorphisms and IS. What is more, future researches should be designed to elucidate the mechanism by which APOE ε4 mutation adds the risk of IS.

Data Availability

Data presented within the paper and the supplementary materials contributed to the findings in our study. They are all are available from our corresponding author for reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

Authors’ Contributions

The conceptualization was done by S.-Y. Q., K. S., Y.-H.C., and X.C.; the methodology was done by D.-S. T., D.-J. P., and C. Q.; K. S., Y-H.C., H.-H. Y, and X.C. took care of the software; meta-analysis was done by D.-S. T., D.-J. P., C. Q., S.-Y. Q., K. S., and X.C.; writing—original draft preparation—was done by S.-Y. Q., K. S., Y.-H.C., and H.-H. Y.; writing—review and editing—was done by D.-S. T., D.-J. P., and C. Q. All authors have read and agreed to the published version of the manuscript. Su-Ya Qiao and Ke Shang contributed equally to this work.

Acknowledgments

We are deeply grateful for all our colleagues of the Department of Neurology in Tongji Hospital.

Supplementary Materials

Supplementary material Table 1: fifteen of the included studies provide data about different subtypes of IS: LAA, SAD, and CE. Supplementary material Table 2: Newcastle-Ottawa Scale (NOS) score of included studies. Supplementary material Table 3: PRISMA list of our meta-analysis. Supplementary material Table 4: sensitivity analysis of the association between ApoE gene polymorphisms and IS. Supplementary material Table 5: publication bias and heterogeneity of our meta-analysis. Supplementary material Figure 1: funnel plots for studies included in Figures 2A–G. Supplementary material Figure 2: funnel plots for studies included in Figures 3A–G. Supplementary material Figure 3: results of meta-regression. (Supplementary Materials)