Study selection
Of the 5645 articles screened for eligibility, 23 studies with 15,382,537 individuals were included in the final meta-analyses [4,5,6,7,8, 10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27] (Fig. 1). From these 23 studies, 21 studies were for singular cardiovascular events [4,5,6,7,8, 10,11,12,13,14,15,16,17,18,19,20,21, 23, 24, 26, 27], 11 for cardiovascular death [4, 6, 7, 10,11,12, 14, 20, 23, 24, 26], 8 for MI events [4, 5, 13,14,15,16,17, 21], 7 for any stroke event [5, 6, 14,15,16,17, 21], and 10 for composite CVD outcomes [4,5,6, 8, 11, 12, 14, 22, 25, 27]. Of these 23 studies, 17 investigated body weight variability [4,5,6, 11, 13,14,15,16,17,18,19,20,21,22, 25,26,27], whilst 7 investigated BMI variability [7, 8, 10, 12, 20, 23, 24]. The study performed by Nam et al. [20] contributed data towards both body weight and BMI variability analyses [20]. Despite meeting inclusion criteria, estimated RRs from the 2019 study by Oh et al. were not included in the final analyses due to only recording estimated RRs per + 1 increase in average successive variability (ASV), and as such not compatible with our analysis [31].
Study characteristics
The recorded average (mean or median) age of the participants within studies ranged from 35 to 72 years. The average weight of participants ranged from 63.5 kg to 92.5 kg. The average BMI of participants ranged from 22.1 to 33.2 kg/m2. Almost every study included in this analysis had > 50% male participation, except for 3 [8, 11, 22]. Of the included studies, 12 had > 50% White participants [4, 5, 8, 11, 13,14,15, 22,23,24, 26, 27], and 11 had > 50% East Asian participants [6, 7, 10, 12, 16,17,18,19,20,21, 25]. The participants of 6 of the studies had been previously diagnosed with T2D [13,14,15, 18, 21, 27], whilst 1 more provided separate statistics for the participants within the study who had been previously diagnosed with T2D [12]. The average follow up time for the included studies ranged from 3.7 to 32 years (See Additional file 1: Table S2). Of the 23 studies, 16 were judged as having a high quality (≥ 7) based on the Newcastle–Ottawa scale (See Additional file 1: Table S4). Sponholtz et al. separated their participant population based upon whether they were obese, and as such generated 2 reports per cardiovascular outcome. Similarly, Wannamethee et al. stratified their population based upon whether their participants initially lost or gained weight [24], Lissner et al. stratified their population by gender [11], and Youk et al. stratified by diabetes status [12], and as such each of these studies produced 2 reports per cardiovascular outcome. In total, 58 reports regarding cardiovascular outcomes were collected from the 23 studies. Of these 58 reports, 47 reported on a singular (i.e. non-compound) recorded CV event, and 11 reported composite CVD outcomes. Of the 47 reports of singular CV events, 17 reports were for CV death, 8 for MI, and 7 for stroke.
Body weight and BMI variability were associated with increased risk of any cardiovascular event
A total of 21 studies (15 for body weight and 6 for BMI) consisting of 15,141,102 individuals investigated the association between weight/BMI variability and any CV outcomes. Compared to the least variable group, the summary RR for any CV outcomes for people in the most variable group of body weight was 1.27 (95% CI 1.17–1.38; P < 0.0001; I2 = 97.28%; P < 0.0001 for heterogeneity; Fig. 2). A similar summary RR of 1.39 (95% CI 1.17–1.64; P < 0.0001; I2 = 76.39%; P < 0.0001 for heterogeneity; See Additional file 1: Figure S2a) was found when the variability was defined using BMI. The summary RR estimate for the association between per + 1 SD increase in unit of body weight variability and any cardiovascular event was 1.16 (95% CI 1.06–1.26; P < 0.0001; I2 = 94.70%; P = 0.0013 for heterogeneity; See Additional file 1: Figure S1a).
Body weight and BMI variability were associated with increased risk of cardiovascular death
A total of 11 studies (5 for body weight variability and 6 for BMI variability) consisting of 633,592 participants investigated the association between weight/BMI variability and risk of CV death. Compared to the least variable group, the summary RR for CV deaths for the most variable group of body weight and BMI was 1.29 (95% CI 1.03–1.60; P < 0.0001; I2 = 55.16%; P = 0.062 for heterogeneity; Fig. 3) and 1.27 (95% CI 1.09–1.49; P = 0.0027; I2 = 68.51%; P = 0.002 for heterogeneity; See Additional file 1: Figure S2b), respectively. The summary RR for CV Deaths per + 1 SD increase in body weight variability was 1.11 (95% CI 1.02–1.21; P = 0.0132; I2 = 49.66%; P for heterogeneity = 0.1359; See Additional file 1: Figure S1b).
Body weight variability was associated with increased risk of myocardial infarction
There were eight studies with a total population of 5,742,933 that investigated the association between body weight variability and MI. The summary RR for MI associated with being in the most variable strata of body weight compared to the least variable strata was 1.32 (95% CI 1.09–1.59; P = 0.0037; I2 = 97.14%; P for heterogeneity < 0.0001; Fig. 3), and the summary RR per + 1 SD increase in body weight variability was 1.14 (95% CI 0.92–1.42; P = 0.2234; I2 = 82.32%; P for heterogeneity = 0.0174; Additional file 1: Figure S1b). No study that investigated BMI variability reported RR for MI.
Body weight variability was associated with increased risk of stroke
There were 7 studies consisting of 5,779,027 subjects that investigated the association between body weight variability and risk of stroke. The summary RR for stroke associated with being in the most variable strata of body weight compared to the least variable strata was 1.21 (95% CI 1.19–1.24; P < 0.0001; Fig. 3). Significant heterogeneity was detected in this analysis (I2 = 0.06%; P for heterogeneity = 0.0073). No study that investigated BMI variability recorded RR estimates for stroke outcomes. Similarly, it was not possible to perform a meta-analysis on the risk of stroke per + 1 SD increase in body weight variability, as there is only 1 study that reported RR for stroke [5].
Body weight variability and risk of composite cardiovascular outcomes
Eight studies consisting of 339,566 participants reported association between body weight variability and composite CV outcomes. The RR for composite CVD outcomes associated with being in the most variable body weight group compared to the least variable was 1.36 (95% CI 1.08–1.73; P = 0.01; I2 = 92.41%; P for heterogeneity < 0.0001; Fig. 3), and the RR of composite CVD outcomes associated per + 1 SD increase in body weight variability was 1.14 (95% CI 1.04–1.25; P = 0.0047; I2 = 91.77%; P for heterogeneity < 0.0001; See Additional file 1: Figure S1b). There was only one study that investigated the association between BMI variability and composite CV outcome [8] and thus a meta-analysis is not performed.
The association between weight variability and CV outcomes was not modified by ethnicity or diabetes status
Given ethnicity and diabetes status are known risk factors for CV outcomes [32, 33], we performed subgroup analyses stratified by ethnicity and diabetes status. Due to a lack of data on other ethnicities, this analysis was performed in White Europeans and East Asians. We observed that both ethnicity and diabetic status generally had no significant effect on the observed association between body weight variability and CV events (See Additional file 1: Figures S3 and S4). Compared to the group with the least variability, a significantly higher risk of any CV event was observed in the group with the highest degree of body weight variability in both the White (RR = 1.42; 95% CI 1.25–1.62; P < 0.0001; See Additional file 1: Figure S3a) and East Asian populations (RR = 1.16; 95% CI 1.12–1.19; P < 0.0001;
See Additional file 1: Figure S3b). The RR for CV death was not statistically significant in either the Whites (RR = 1.33; 95% CI 0.97–1.83; P = 0.0741; See Additional file 1: Figure S3c) or the East Asians (RR = 1.22; 95% CI 0.90–1.66; P = 0.2022; See Additional file 1: Figure S3d). The RR of composite CVD outcome in East Asians was also not significant (RR = 1.11; 95% CI 0.81–1.52; P = 0.5154; See Additional file 1: Figure S3d). This is most likely due to a lack of power after stratification.
Similar results were obtained after stratifying by diabetes status. The RR for any CV event in the most variable group compared to the least variable group with diabetes was 1.25 (95% CI 1.13–1.38; P < 0.0001; I2 = 98.03%; P for heterogeneity < 0.0001) and in non-diabetics it was 1.29 (95% CI 1.14–1.46; P < 0.0001; I2 = 98.03%; P for heterogeneity < 0.0001). No differential association by diabetes status was also observed between weight variability and MI, however the association between weight variability and stroke was observed to be insignificant in the non-diabetic population (RR = 1.31; 95% CI 0.99–1.72; P = 0.0566; I2 = 98.35%; P for heterogeneity = 0.0029) (See Additional file 1: Figure S4c and d). Only one of the papers included in this analysis investigated cardiovascular death or composite cardiovascular outcomes in individuals with type II diabetes [14], and as such these outcomes were not suitable for subgroup meta-analysis. Therefore, the effect of diabetes status on the risk of these outcomes cannot be analysed.
Sensitivity analyses
Given different studies capture weight variability using different metrics (i.e. ASV, SD, CoV, VIM, RMSE), we performed sensitivity analysis after stratification by exposure definition. Overall, exposure definition had little impact on the results, with a few notable exceptions (See Additional file 1: Figure S5a–i). When ASV was used, the summary RRs associated with the most variable strata of body weight (compared to the least variable) for both MI and stroke were significantly higher, with the RR for MI 1.97 (95% CI 1.60–2.44; P < 0.0001; See Additional file 1: Figure S5b) and stroke 2.17 (95% CI 1.57–3.00; P < 0.0001; See Additional file 1: Figure S5b). However, when ASV was used as a measure of variability in the per + 1 SD increase in body weight variability analysis, the summary RR for MI became insignificant (See Additional file 1: Figure S5i). Similarly, when ASV was used as a measure of variability, the summary RR for CV death became insignificant (See Additional file 1: Figure S5b and i). These changes are most likely explained by a lack of power in these sub-analyses. Of note however is the general trend of decreased heterogeneity of results observed after stratifying by metric of effect, which suggests that this difference in study methodology is a large contributor to the difference in result we observe between studies.
An important question to address when interpreting these results is how much the specific quantile used in the most variable group versus the least variable group affects the observed association. To explore this, we stratified the most variable versus least variable body weight analysis based on the quantile used by the included studies. Of the 15 included studies that compared quantiles, 6 compared quintiles, 7 compared quartiles, 1 compared tertiles, and 1 compared medians (See Additional file 1: Table S3). We therefore stratified the analysis by studies that investigated quintiles or quartiles. Due to a lack of included studies investigating BMI variability, this sensitivity analysis was only performed on studies that investigated body weight variability. We also removed studies that did not compare the highest quantile to the lowest quantile, and instead compared the top quantile against the combined lower quantiles [18, 20]. This results in 3 studies being included in the quintile analysis [5, 14, 25], and 5 included in the quartile analysis [15,16,17, 19, 21]. Compared to the quintile with the least variability, a significantly higher risk of any CV event was observed in the quintile with the highest degree of body weight variability (RR = 1.86; 95% CI 1.41–2.44; P < 0.001; I2 = 54.40%; P for heterogeneity = 0.0725). A similar result was found in studies that compared quartiles (RR = 1.18; 95% CI 1.11–1.25; P < 0.0001; I2 = 95.30%; P for heterogeneity < 0.0001) (Additional file 1: Figure S6a and b). When comparing the risks of secondary outcomes, all outcomes remained significant in both the quintile and quartile strata, however a general trend was found where the risk observed in the quintile group was higher than that of the same outcome in quartile group. While this difference in RR between studies that compared quintiles and quartiles could reflect the true impact of weight variability on CV outcomes in the top 20% of the population, it could also be due to differences in sample sizes or unseen bias as studies that compared quartiles are largely from similar studies. Therefore, further investigation in a well powered prospective study is warranted (Additional file 1: Figure S6c and d).
One of the key questions when investigating weight variability is to what degree is the increased risk of CVD observed due to a general increase or decrease in body weight rather than a variability in weight. To explore this question, we stratified our analysis into studies that adjusted for average BMI or change in BMI versus studies that did not control for these covariates. When compared to the group with the least variability, a significantly higher risk of any CV event was observed in the group with the highest degree of body weight variability in studies that adjusted for BMI (RR = 1.60; 95% CI 1.37–1.87; P < 0.0001; I2 = 52.50%; P for heterogeneity = 0.0104; See Additional file 1: Figure S8a) and studies that did not (RR = 1.16; 95% CI 1.09–1.23; P < 0.0001; I2 = 95.63%; P for heterogeneity < 0.0001;
See Additional file 1: Figure S8c). When comparing secondary outcomes between these two strata, we found that the observed relative risk were similar, with a few exception (Additional file 1: Figure S8b and d). After stratification, the relative risk of CV death associated with weight variability was found to be insignificant in both groups, however this is most likely due to reduced power post-stratification. Similarly, the risk of composite cardiovascular outcomes was found to be insignificant in the unadjusted strata, again most likely due to reduced power post-stratification. An interesting observation is that the relative risks reported in the adjusted group were typically higher than those in the unadjusted group, however, this could be mostly explained by other differences between these two groups. For example, all of the studies in the adjusted analysis included participants of White Europeans and 80% of the participants in the unadjusted analysis had East Asian ancestry. In addition, the adjusted studies were generally smaller than the unadjusted studies, with a total adjusted population of 50,095 and an unadjusted population of 14,650,699.
Of the studies included in this analysis, several were performed on populations with pre-existing CVD. It is therefore important to investigate whether weight variability has a differential effect on risk of CVD among populations with high risk of CVD at baseline versus populations with no known CVD. Compared to the group with the least variability, a significantly higher risk of any CV event was observed in the group with the highest degree of body weight variability in both populations with pre-existing CVD (RR = 1.59; 95% CI 1.29–1.96; P < 0.0001; I2 = 65.75%; P for heterogeneity = 0.0035; See Additional file 1: Figure S7a) and no known CVD (1.19; 95% CI 1.11–1.26; P < 0.0001; I2 = 95.49%; P for heterogeneity < 0.0001; See Additional file 1: Figure S7c). The observed relative risk of secondary outcomes between these two groups were largely similar, with a few exceptions (See Additional file 1: Figure S7b and d). The relative risk of CV death associated with weight variability was found to be insignificant in the group with pre-existing CVD risk. The risk of composite cardiovascular outcomes associated with weight variability in both groups was found to be insignificant after stratification. However, these differences could be explained by a reduced power after stratification.
To further explore this question, we performed a univariate meta-regression to investigate whether the observed relative risk of CVD associated with weight variability was correlated on the average age of the study populations (See Additional file 1: Figure S9a–e). The relative risk of any CV event associated with weight variability was found to be significantly positively correlated with average age (β = 0.0081; P = 0.0345; R2 = 0.028). No significant correlation between average age and relative risk associated with weight variability was observed for any of the secondary outcomes.
Heterogeneity and bias analysis
As heterogeneity was significant in the analysis of the RRs of the primary and secondary outcomes associated with being in the top strata of body weight variability, it was important to investigate whether this heterogeneity was due to publication bias. As such, Egger’s regression test and funnel plots were created for these analyses. Egger’s regression found no funnel plot asymmetry for the CV event (z = 1.7567; P = 0.079; See Additional file 1: Figure S10a), CV death (z = −1.0027; P = 0.316; See Additional file 1: Figure S10b), MI (z = 1.1849; P = 0.236; See Additional file 1: Figure S10c), or the most composite CVD outcomes (z = −1.7294; P = 0.0837; See Additional file 1: Figure S10e) analyses, however, significant asymmetry was found for the analysis of stroke (z = 2.9287; P = 0.0034; See Additional file 1: Figure S10d). As such, the Duval and Tweedie trim-and-fill method was employed in order to estimate the effect that hypothetical missing publications would have on the summary RR estimate. This analysis found no change in the estimated summary RR of any cardiovascular event associated with degree of body weight variability (RR = 1.21; 95% CI 1.19–1.24; P < 0.0001; See Additional file 1: Figure S10f).
To assess the effect that low quality papers had on the results, the 7 studies that scored < 7 using NOS were removed and then the primary and secondary outcomes were analysed a second time. The effect sizes remain the same after the removal of low-quality studies (See Additional file 1: Figure S11).