Television Viewing and Incident Cardiovascular Disease: Prospective Associations and Mediation Analysis in the EPIC Norfolk Study

Abstract

Background

Although television viewing time is detrimentally associated with intermediate cardiovascular risk factors, the relationship with incident total (i.e. combined fatal and non-fatal) cardiovascular disease (CVD), non-fatal CVD and coronary heart disease is largely unknown. This study examined whether television viewing time is associated with these three outcomes, independently of physical activity energy expenditure and other confounding variables.

Methodology/Principal Findings

A population-based cohort of 12,608 men and women (aged 61.4±9.0), free from stroke, myocardial infarction and cancer at baseline in 1998–2000 were followed up until 2007 (6.9±1.9 years). Participants self-reported education, smoking, alcohol use, antihypertensive, lipid lowering and antidepressant medication, disease history, total energy intake, sleep duration, physical activity and television viewing. BMI, waist circumference, blood pressure, triglycerides, HDL cholesterol and glycated haemoglobin (HbA1c) were measured by standardized procedures; a clustered metabolic risk score was constructed. Every one hour/day increase in television viewing was associated with an increased hazard for total (HR = 1.06, 95%CI = 1.03–1.08; 2,620 cases), non-fatal CVD (HR = 1.06, 95%CI = 1.03–1.09; 2,134 cases), and coronary heart disease (HR = 1.08, 95%CI = 1.03–1.13; 940 cases), independent of gender, age, education, smoking, alcohol, medication, diabetes status, CVD family history, sleep duration and physical activity energy expenditure. Energy intake, BMI, waist circumference, blood pressure, triglycerides, HDL cholesterol, HbA1c and the clustered metabolic risk score only partially mediated these associations.

Conclusions

These results indicate that the most prevalent leisure time (sedentary) behaviour, television viewing, independently contributes to increased CVD risk. Recommendations on reducing television viewing time should be considered.

About the Authors:

Katrien Wijndaele

Affiliations Medical Research Council Epidemiology Unit, Cambridge, United Kingdom, Research Foundation Flanders, Brussels, Belgium, Department of Movement and Sports Sciences, Ghent University, Ghent, Belgium

Søren Brage

Affiliation: Medical Research Council Epidemiology Unit, Cambridge, United Kingdom

Hervé Besson

Affiliation: Medical Research Council Epidemiology Unit, Cambridge, United Kingdom

Kay-Tee Khaw

Affiliation: Department of Public Health and Primary Care, Institute of Public Health, University of Cambridge, Cambridge, United Kingdom

Stephen J. Sharp

Affiliation: Medical Research Council Epidemiology Unit, Cambridge, United Kingdom

Robert Luben

Affiliation: Department of Public Health and Primary Care, Institute of Public Health, University of Cambridge, Cambridge, United Kingdom

Amit Bhaniani

Affiliation: Department of Public Health and Primary Care, Institute of Public Health, University of Cambridge, Cambridge, United Kingdom

Nicholas J. Wareham

Affiliation: Medical Research Council Epidemiology Unit, Cambridge, United Kingdom

Ulf Ekelund

* E-mail: [email protected]

Affiliation: Medical Research Council Epidemiology Unit, Cambridge, United Kingdom

Introduction

Sedentary behaviour (i.e. excessive sitting, as distinct from insufficient physical activity), occupies a substantial part of the waking day in modern society (e.g. 55% in the US) [1]. Increasing evidence from prospective studies indicates detrimental associations between excessive sitting, and television viewing time in particular, and intermediate cardiovascular risk factors, which are independent of physical activity and other relevant covariates. These include weight gain and incident obesity [2]–[4], dyslipidaemia [5], hypertension [6] and insulin resistance and diabetes [3], [7], [8]. Potential mechanistic explanations for these associations are beginning to be unravelled [9].

Few prospective studies however have so far examined independent sedentary behaviour effects on more downstream incident cardiovascular disease (CVD). Excessive total or non-occupational sitting time increased cardiovascular mortality hazard in Canadian [10] and middle-aged US adults [11] and total CVD hazard in US post-menopausal women [12] following adjustment for leisure time physical activity. Specifically for television viewing, an increased hazard of prolonged viewing time on cardiovascular mortality was found in Australian [13] and UK adults [14], independent of leisure time exercise and total physical activity respectively, but not in US men from the Aerobics Center Longitudinal Study [15]. One recent study in Scottish adults examined the prospective associations between leisure screen time and total (i.e. combined fatal and non-fatal) incident CVD [16]. However, so far it has not been studied yet whether television viewing is independently associated with incident coronary heart disease or other specific types of CVD, such as cardiac failure and stroke. Compared to previous studies focusing on cardiovascular mortality [13]–[15], examining the relationship between television viewing and non-fatal CVD is of additional value, as morbidity and incapacity caused by non-fatal cardiovascular events are major contributors to the economic burden for society and the health care system [17], [18]. In addition, the first CVD event is most often more proximal in the cardiovascular disease process to the exposure of interest, and so may provide a more complete picture of the potential public health hazard associated with television viewing in terms of incident CVD.

Similar to the US and Australia, television viewing is the predominant leisure time sedentary behaviour in British adults and one of the three main activities besides sleeping and working [19]–[21]. Therefore it is likely to be representative of total leisure time spent sitting [22]. Further, it might be more susceptible to voluntary change compared to other sources of excessive sitting, such as occupational sitting. Moreover, compared to sitting in general, television viewing is associated with additional factors potentially mediating associations with CVD, such as increased energy intake by exposure to unhealthy food advertisements [23]. This may result in detrimental associations which differ in effect size from those caused by sitting per se.

In the EPIC Norfolk prospective population study, we examined the independent association between television viewing time and incident total cardiovascular, non-fatal cardiovascular and coronary heart disease events in a population of healthy middle-aged white adults. We determined whether any such associations were independent of total (at home, to work, at work and during leisure-time) physical activity energy expenditure (PAEE) and other relevant confounders. Furthermore, we also examined to what extent these associations were mediated by a range of intermediate cardiovascular risk factors.

Results

Descriptive characteristics

During 87,572 person-years of follow-up (men: 36,846, women: 50,726), 2,620 participants developed incident CVD (1,351 men, 1,269 women). Out of those 2,620 incident CVD cases, 2,134 (81.5%) were non-fatal CVD events (men: 1,064, women: 1,070) and 940 (35.9%) were (fatal or non-fatal) coronary heart disease events (men: 625, women: 315). As shown in Table 1, participants who developed CVD watched television at baseline for on average half an hour/day more than those who did not (any and non-fatal CVD: 0.4, coronary heart disease: 0.5 hours/day). Comparing baseline characteristics by tertiles of the exposure of interest (see Table 2), higher levels of television were associated with gradually less favourable profiles of most covariates in men and women. Television viewing (hours/day) and PAEE (MET*hours/week) were weakly inversely correlated (Spearman correlation: −0.17, P<0.01).

[Figure omitted. See PDF.]

Table 1. Descriptive characteristics (mean (SD) or No. (%)) at baseline between incident cardiovascular event cases and non-cases in 12,608 men and women in EPIC Norfolk, 1998–2007.

https://doi.org/10.1371/journal.pone.0020058.t001

[Figure omitted. See PDF.]

Table 2. Descriptive characteristics (mean (SD) or No. (%)) at baseline by television viewing tertiles (lowest: <2.5, middle: 2.5–3.6, highest: >3.6 hours/day) in 12,608 men and women in EPIC Norfolk.

https://doi.org/10.1371/journal.pone.0020058.t002

Cox proportional hazards regression

Figure 1 shows the increases in age-adjusted CVD event rates per 10,000 person-years by 1-hour increments in television viewing, suggesting a linear association between television viewing time and all 3 outcomes. Adding the quadratic television viewing term to Model C did not result in a significant improvement in model fit for any of the 3 outcomes.

[Figure omitted. See PDF.]

Figure 1. Age-adjusted cardiovascular disease event rates per 10,000 person-years by 1-hour increments of television viewing (hours/day).

N television categories: 0–0.99: 800; 1–1.99: 2,126; 2–2.99: 3,180; 3–3.99: 3,057; 4–4.99: 2,040; 5–5.99: 926; ≥6: 479.

https://doi.org/10.1371/journal.pone.0020058.g001

As shown in Table 3, television viewing (hours/day) was positively associated with incident total CVD, non-fatal CVD and coronary heart disease, independent of age and gender (Model A), education, smoking, alcohol, hypertension, dyslipidaemia and antidepressant medication, baseline diabetes status, family history of CVD, sleep duration (Model B), and total PAEE (Model C). The increase in hazard associated with every one hour increase in television viewing time was 6% for non-fatal and total CVD, and 8% for coronary heart disease. These results were virtually unchanged after excluding participants (141 men and 96 women) who encountered any CVD event within the first year of follow-up (data not shown). Adjustment for social class instead of education level did not change results (data not shown).

[Figure omitted. See PDF.]

Table 3. Hazard ratios (95%CI) for any cardiovascular, non-fatal cardiovascular and coronary heart disease events per hour/day increase in television viewing in 12,608 men and women in EPIC Norfolk, 1998–2007.

https://doi.org/10.1371/journal.pone.0020058.t003

After introducing multiplicative interaction terms to Model C, there were no significant interactions on any of the three outcomes between television viewing time and physical activity, BMI, gender, or education. A significant interaction by age and clustered metabolic risk was found for total CVD (age: P = 0.02; clustered metabolic risk: P = 0.03) and non-fatal CVD (age: P = 0.03; clustered metabolic risk: P<0.01), but not for coronary heart disease. Table 4 illustrates associations between television viewing and total CVD derived from Model C for subgroups. Consistent with the results of the interaction tests, very similar hazard ratios were found between subgroups for most characteristics. Older participants and participants with higher clustered metabolic risk respectively showed a somewhat lower hazard ratio for the association between television viewing time and incident total CVD compared to their younger and lower risk counterparts, although associations were still in the same direction.

[Figure omitted. See PDF.]

Table 4. Hazard ratios (95%CI) for any cardiovascular disease event per hour/day increase in television viewing in subgroups according to physical activity level, body mass index, gender, education level, age and clustered metabolic risk, in 12,608 men and women in EPIC Norfolk, 1998–2007.

https://doi.org/10.1371/journal.pone.0020058.t004

Adding hormone replacement therapy use to model C in women resulted in the same hazard ratio for television viewing time (1.05 (1.01–1.09)), indicating that results in women were also independent of this factor.

As a second aim, potential mediation effects were examined by adding continuous energy intake, BMI, the clustered metabolic risk score and its individual metabolic risk variables to Model C. Preliminary analyses showed significant correlations between continuous television viewing time and these variables (all P<0.001; negative correlation with HDL-cholesterol, positive correlation with all other variables), and significant differences for these variables between incident cases and non-cases for the three outcomes (see Table 1), suggesting a mediation effect by these variables. As shown in Table 5, although an attenuation of the hazard ratios (95%CI) for television viewing was found in most models, they remained significant. This suggests that these characteristics only partially mediated the association between television viewing time and incident CVD. Additional inclusion of intake of fruit and vegetables (grams/day) and saturated fatty acids (% total energy intake) as indicators of an (un)healthy diet in the model examining a potential mediation effect of dietary intake (Model C+total energy intake) resulted in very similar hazard ratios (95%CI) for television viewing time: 1) total CVD: 1.07 (1.03–1.10), P<0.001; 2) non-fatal CVD: 1.07 (1.03–1.11), P<0.001; and 3) coronary heart disease: 1.09 (1.04–1.15), P<0.001.

[Figure omitted. See PDF.]

Table 5. Hazard ratios (95%CI) for any cardiovascular, non-fatal cardiovascular and coronary heart disease events per hour/day increase in television viewing in 12,608 men and women in EPIC Norfolk, 1998–2007, after additional adjustment for potential mediating variables.

https://doi.org/10.1371/journal.pone.0020058.t005

Additionally, we examined the association between television viewing time and incident cardiac failure (229 fatal or non-fatal events; men: 147, women: 82) and stroke (261 fatal or non-fatal events; men: 127, women: 134). Television viewing (hours/day) was positively associated with incident cardiac failure, independent of age and gender (Model A: HR (95%CI): 1.15 (1.05–1.25)), education, smoking, alcohol, hypertension, dyslipidaemia and antidepressant medication, baseline diabetes status, family history of CVD, sleep duration (Model B: HR (95%CI): 1.11 (1.01–1.21)), and total PAEE (Model C: HR (95%CI): 1.10 (1.01–1.20)). Positive but non-significant associations were found for incident stroke (HR (95%CI): Model A: 1.06 (0.98–1.16), Model B: 1.04 (0.95–1.13), Model C: 1.03 (0.95–1.12)).

Discussion

In this population-based sample of healthy middle-aged white adults, higher television viewing time was independently associated with an increased hazard of incident total CVD, non-fatal CVD and coronary heart disease. These associations were independent of a continuous measure of overall (during leisure-time, to work, at work and at home) PAEE. Although PAEE was self-reported, this may suggest that these associations are not just the result of a compensatory reduction in physical activity. Furthermore, associations were equally strong in high and low physically active individuals and may support the need for recommendations on limiting television viewing time, in addition to the established physical activity recommendations [9]. So far, two studies reported independent associations between television viewing time and cardiovascular mortality, the final outcome in the CVD process [13], [14]. The recent findings from the Scottish Health Survey [16] indicated that leisure screen time (viewing television or another type of screen such as computer or video game) was associated with total (fatal and non-fatal combined) incident CVD (215 cases) in 4,512 adults (aged ≥35) followed up for 4.3 (±0.5) years. Despite minor differences in the exposure variable and adjustment strategy, the hazard ratio (95%CI (raised to power 60 to transform from minutes/day to hours/day)) for continuous screen time from the fully adjusted model for total incident CVD was very similar to the one in our study, more specifically 1.06 (1.01–1.11) [16]. Compared to the literature, the current study provides important novel findings by showing independent associations between television viewing time and not only non-fatal CVD, imposing a large burden on society, but also coronary heart disease. Additional analyses for incident cardiac failure and stroke suggested that the strength of association may differ to some extent according to the subtype of CVD outcome examined. Although results for incident cardiac failure and stroke should be interpreted with caution given the limited number of incident cases for both outcomes, these novel findings should encourage future larger prospective studies with longer follow up time to further explore these independent associations for specific types of CVD.

Every hour increase (or decrease) in television viewing time was independently associated with a 6% higher (lower) hazard for total and non-fatal CVD and an 8% higher (lower) hazard for coronary heart disease. Although this effect size is relatively modest, television viewing is a highly prevalent behaviour, both at intra- (e.g. on average nearly three hours/day in Great Britain) and inter-individual levels (e.g. reported by 80% of Great Britain adults) [19]. It therefore provides an important target for health-related behaviour change and cardiovascular health benefits at a population level.

Examining pathways explaining independent effects from television viewing time and (prolonged) sitting per se on cardiovascular health is a relatively unexplored area. A unique and important finding of this study is that a large range of intermediate behavioural and biological characteristics only partially mediated the observed associations. Effect attenuation by these characteristics is supported by results from previous studies. These include independent associations of television viewing time and objectively measured sedentary behaviour with increased energy intake [23], weight gain/incident obesity [2]–[4], dyslipidaemia [5], hypertension [6] and insulin resistance/diabetes [3], [7], [8], which are independent of physical activity and other relevant covariates; partial mediation by triglycerides and HDL cholesterol is additionally supported by findings in rat-models showing acute effects on skeletal muscle lipoprotein lipase which are, both quantitatively and qualitatively, very specific to enforced muscle inactivity conditions (simulating human sedentary behaviour) and largely uninfluenced by exercise conditions [9]. However, none of all potential mediators in the current study could completely account for the associations between television viewing and the outcomes. Further research is required on the (pathophysiological) processes which are specific for sedentary behaviours and not just the effect of absence of physical activity. As television viewing time and these intermediate characteristics were measured at the same time point, we cannot however disentangle direction of causality and therefore true mediation by these variables. Future observational and experimental studies with multiple time-points are warranted to further disentangle these pathways.

An important strength of this study is the large population-based sample of men and women, in which a relatively large proportion developed incident CVD. This enabled us to study three different types of CVD which have not been associated with television viewing time before, and also preliminary examine associations for incident cardiac failure and stroke. Given the sample size, we were able to exclude everyone with baseline history of self-reported and/or diagnosed stroke, myocardial infarction and cancer, and in sensitivity analysis those who developed CVD within the first year of follow-up. This, together with the prospective design and extensive adjustment for confounders, supports a causal inference of the associations found and reduces the possibility of reverse causality. We adjusted for overall PAEE (four different domains), a measure which has previously been shown to predict cardiovascular mortality in this cohort [27], in addition to demonstrated criterion-validity [26]. Furthermore, a comprehensive set of potential mediators were considered, including energy intake and indicators of diet quality, a wide spectrum of objectively measured cardiovascular risk factors and a summary cardio-metabolic risk score. A potential weakness that merits discussion is that television viewing time was self-reported, and only measured at baseline. Therefore, measurement error and misclassification might have resulted in an underestimation of the true associations. Further, although a large range of potential confounders was considered, residual confounding might exist for those that relied on self-report or were categorical. Additionally, although television viewing is the most-prevalent leisure time behaviour [19]–[21] and likely to be representative of leisure time sedentary behaviour [22], other types of sedentary behaviour correlating with television viewing time might also have contributed to residual confounding. Based on recent findings in American adults, showing significant but fair associations between TV viewing time and total sedentary time derived from accelerometry, TV viewing time has been suggested to be a potentially useful indicator for total sedentary time in epidemiological analyses. Of particular relevance to our population sample, rank order correlations in mid- and old-aged adults (Spearman's rho: 40–59 years olds: 0.17 (P<0.001); 60+ year olds: 0.23 (P<0.001) were stronger compared to those in younger adults (20–39 years old: 0.05 (P = 0.04)) [31]. As information on types of video or TV programs watched (e.g. highly charged sports events versus other shows) was not available, we could not examine the possibility of differences in strength of association according to type of entertainment content. The specific mediation effect of eating while watching TV viewing was also not studied as this behaviour was not measured separately, but we did adjust for general dietary variables. Finally, as we excluded participants with baseline chronic health conditions to minimize the possibility of reverse causality and increase internal validity, the associations found are relative to this cohort of relatively healthy middle-aged white adults. Additionally, non-randomly missing data might have biased results, especially in models including dietary intake, which involved exclusion of relatively large proportions of participants. Further prospective large population-based studies in different settings, using repeated measurements, with longer duration of follow-up and preferably using objective measures of physical activity and sedentary time are needed to extend the present results and confirm dose-response effects in different populations.

To conclude, television viewing time, the predominant leisure activity in modern society [19]–[21], was associated with increased risk of total incident CVD, non-fatal CVD and coronary heart disease in healthy middle-aged white adults, independent of total level of physical activity and other relevant confounders. Furthermore, an extensive range of intermediate CVD risk factors only partially mediated these independent associations. These observations suggest the need to consider separate public health recommendations on reducing television viewing time.

Materials and Methods

Study participants

EPIC Norfolk is part of the 10-country collaborative European Prospective Investigation into Cancer and Nutrition (EPIC) study. A cohort of 25,633 residents of Norfolk (UK) aged 45–79 and recruited via general practitioners agreed to participate between 1993–1997. A detailed description of the study design and cohort characteristics has previously been published [24]. Participants were invited for a follow-up assessment between January 1998–October 2000. This follow-up also encompassed introduction of the EPIC physical activity questionnaire (EPAQ2), which is a more comprehensive physical activity instrument also including questions on television viewing as described in more detail below. This follow-up examination, including 15,784 (61.6%) attendees, is the baseline for the current analyses. Complete data for television viewing time, PAEE, education level, smoking status, alcohol consumption, medication use, baseline diabetes status, parental history of CVD and sleep duration were provided by 15,021 participants. Those with diagnosed or self-reported baseline history of stroke (n = 367), myocardial infarction (n = 635), other vascular disease (n = 3) and/or cancer (n = 1,449) were excluded. As a result, 12,608 participants (5,465 men, 7,143 women) were included in the current analyses. The study complies with the Declaration of Helsinki, and was approved by The Norwich District Health Authority Ethics Committee. All participants provided written informed consent.

Measures

Incident CVD events.

All EPIC Norfolk participants are followed up for fatal and non-fatal CVD events. The present study includes follow-up until 31st March 2007, covering 6.9±1.9 (mean±SD) follow-up years. The Office of National Statistics, UK, flagged all individuals for death certification and trained nosologists coded death certificates according to the International Classification of Disease (ICD). Cardiovascular death was defined as ICD 410–448 (ICD 9) or ICD I10–I79 (ICD 10) as underlying cause of death, which comprizes coronary heart disease (410–414 (ICD 9) or I20–I25 (ICD 10)), stroke (430–438 (ICD 9) or I60–I69 (ICD 10)), cardiac failure (428 (ICD 9) or I50 (ICD 10)) and other vascular causes. Cause-specific (same ICD coding) hospital admission was determined via ENCORE (East Norfolk Commission Record, the hospital admissions database kept by the East Norfolk Health Commission) [25], using individuals' unique National Health Service number. Death certificates and ENCORE show high accuracy in correctly identifying incident disease, as previously shown in EPIC Norfolk for incident stroke [25].

The primary outcome was any first CVD event during follow-up, a combined end point defined as the first of any of these events: hospital admission or death because of coronary heart disease, stroke or other vascular disease. Those with non-fatal CVD events were identified as having been admitted to hospital for any of those reasons without dying during follow-up. The subgroup of individuals with coronary heart disease events only included those who had coronary heart disease as the underlying cause of hospital admission and/or death during follow-up.

Television viewing time and PAEE.

By means of the EPAQ2, participants self-reported their physical (in)activity behaviour in different sub-domains (at home, to work, at work and during leisure-time), using the past year as a reference frame [26]. As previously described in detail [27], PAEE (MET*hours/week) was calculated in each of these four (mutually exclusive) sub-domains by multiplying participation (hours/week) by the metabolic cost of each activity (metabolic equivalent (MET)) according to Ainsworth et al. [28]. Total PAEE was then determined by summing energy expenditure of the different sub-domains. Time spent watching television and video (hours/week) was calculated based on responses to four questions about watching before and after six pm at week- and weekend-days. The EPAQ2 scored high for repeatability, both in terms of television viewing time and PAEE [26]. It is valid for ranking individuals, as shown by comparison against minute-by-minute heart rate monitoring and maximal aerobic capacity (VO2max) [26]. In the current study, television viewing time and PAEE were expressed in hours/day and MET*hours/day respectively, for ease of interpretation; they were mutually exclusive (i.e. television viewing was not part of PAEE).

Covariates.

Participants self-reported their education level (low, O level, A level, degree), social class (manual, non manual), smoking status (current, former, never), alcohol consumption (units/week), medication for diabetes (yes, no), medication for hypertension (yes, no), medication for dyslipidaemia (yes, no), medication for depression (yes, no), hormone replacement therapy status (current, former, never), baseline history (yes, no) of myocardial infarction, stroke, cancer and diabetes, and family history of CVD (yes, no) by means of a detailed health and lifestyle questionnaire [24]. Sleep duration (hours/day) was assessed using the EPAQ2 [26]. Diagnosed history of cancer was identified by the ECRiC (Eastern Cancer Registry and Information Centre) and the Office of National Statistics. Total energy intake (kj/day), fruit and vegetable intake (g/day) and saturated fatty acids intake (% total energy intake) were derived from a validated 130-item semi-quantitative food frequency questionnaire [29].

As previously described in detail [24], trained nurses measured height, weight and waist circumference following standardized protocols. Body mass index (BMI) was calculated as weight divided by (height)2 (kg/m2). Systolic and diastolic blood pressures were measured in duplicate using an Accutorr sphygmomanometer (Datascore, UK) after three minutes of sitting. A non-fasting venous blood sample was examined for triglycerides and high density lipoprotein (HDL) cholesterol (mmol/L) using RA1000 (Bayer Diagnostics, UK). Glycated haemoglobin (HbA1c) was measured using Diamat ion exchange HPLC (Bio-Rad Laboratories, UK). A clustered metabolic risk score was constructed to summarize cardio-metabolic risk based on continuously distributed indicators of central obesity (waist circumference), dyslipidaemia (triglycerides and HDL cholesterol), hypertension (systolic and diastolic blood pressure), and hyperglycemia (HbA1c) [30]. These variables were standardized (i.e. z-scores were computed (z = ([value−mean]/SD)), after normalizing (log10) triglycerides and HbA1c. Systolic and diastolic blood pressure z-scores were averaged and the HDL z-scores were inverted. Z-scores were summed and the sum was divided by five to express the cardio-metabolic risk score in SD units. Standardization (z-scores) was performed, stratified by sex (i.e. in males, the male-specific mean and SD values were used for each metabolic variable, in females, the female-specific mean and SD values were used for each metabolic variable) including all participants with complete data for each metabolic variable.

Statistical analysis

Analyses were conducted using SPSS15.0 (SPSS, Inc., Chicago, IL) and STATA10.0 (Stata, Corp., TX). Statistical significance was set at P<0.05. Baseline characteristics were compared by incident CVD status (unpaired t-tests and chi-squared tests) and by television viewing tertiles (lowest: <2.5, middle: 2.5–3.6, highest tertile: >3.6 hours/day; one-way ANOVA and chi-squared tests).

To examine the association between baseline television viewing time (hours/day) and incident CVD Cox proportional hazards regression was used. The proportional hazards assumption was checked by examining Schoenfeld residuals and Kaplan-Meier plots for all three outcome variables. The Schoenfeld residuals did not suggest evidence of deviations from proportionality, which was consistent with observations in the Kaplan-Meier plots. To examine linearity of the association between television viewing time and the outcomes, age-adjusted CVD event rates (95%CI) per 10,000 person-years of follow-up were plotted by 1-hour increments in television viewing time. Additionally, a log likelihood-ratio test examined whether addition of a quadratic term for television viewing time to the adjusted model (Model C) resulted in a statistically significant improvement in model fit. Based on evidence for a linear association, hazard ratios (95%CI) per hour/day increase in television viewing time were estimated. Initial models (Model A) were adjusted for baseline age and gender, and were then further adjusted (Model B) for baseline education level, smoking status, alcohol consumption, antihypertensive medication, lipid-lowering medication, antidepressant medication, diabetes status (based on diabetes history and diabetes medication use), family history of CVD and sleep duration. To examine whether television viewing time was independently associated with incident CVD, we then included baseline PAEE in the model (Model C). Analyses were repeated after additionally excluding participants who experienced any CVD event within the first year of follow-up, to examine the possibility that preclinical disease might influence television habit.

Multiplicative interaction terms were added to Model C to examine effect modification by physical activity level (median split: ≤15.7, >15.7 MET*hours/day), BMI (normal weight: <25, overweight/obese: ≥25 kg/m2), gender, education level (low/O level, A level/degree), age (≤60, >60 yrs of age) and clustered metabolic risk (median split: ≤−0.05, >−0.05). As 6.6% of participants had missing data for clustered metabolic risk, this specific model included 11,776 participants.

The following potential mediators of the association between television viewing time and incident CVD were examined by introducing them into Model C: total energy intake (kJ/day; 21% missing; with and without additional inclusion of fruit and vegetable intake (g/day) and saturated fatty acids intake (21% missing)); BMI, waist circumference, triglycerides (5.8% missing), HDL-cholesterol (5.9% missing), systolic and diastolic blood pressure, HbA1c (4.9% missing) and the clustered metabolic risk score (6.6% missing). Finally, total energy intake, BMI and the clustered metabolic risk score were entered simultaneously (25.9% missing) to Model C. Due to additional missing data for some of these potential mediators, the corresponding models did not include all 12,608 participants.

Acknowledgments

We are tremendously grateful to Prof. Sheila Rodwell, principal investigator and key contributor to the EPIC Norfolk study, who sadly passed away before finalisation of this manuscript. The authors also thank Jian'an Luan for statistical advice and EPIC Norfolk staff and participants for their invaluable contributions.

Author Contributions

Conceived and designed the experiments: KW UE NJW K-TK. Performed the experiments: KW. Analyzed the data: KW SJS. Contributed reagents/materials/analysis tools: AB RL SJS. Wrote the paper: KW. Performed statistics: KW SJS. Interpretation of data: KW SB UE. Responsible for data collection and study design, critical input on all versions and approval of final version of the manuscript: KW SB HB KT-K SJS RL AB NJW UE.

Citation: Wijndaele K, Brage S, Besson H, Khaw K-T, Sharp SJ, Luben R, et al. (2011) Television Viewing and Incident Cardiovascular Disease: Prospective Associations and Mediation Analysis in the EPIC Norfolk Study. PLoS ONE6(5): e20058. https://doi.org/10.1371/journal.pone.0020058

References

1. Matthews CE, Chen KY, Freedson PS, Buchowski MS, Beech BM, et al. (2008) Amount of time spent in sedentary behaviors in the United States, 2003–2004. Am J Epidemiol 167: 875–881.CE MatthewsKY ChenPS FreedsonMS BuchowskiBM Beech2008Amount of time spent in sedentary behaviors in the United States, 2003–2004.Am J Epidemiol167875881

2. Ball K, Brown W, Crawford D (2002) Who does not gain weight? Prevalence and predictors of weight maintenance in young women. Int J Obes 26: 1570–1578.K. BallW. BrownD. Crawford2002Who does not gain weight? Prevalence and predictors of weight maintenance in young women.Int J Obes2615701578

3. Hu FB, Li TY, Colditz GA, Willett WC, Manson JE (2003) Television watching and other sedentary behaviors in relation to risk of obesity and type 2 diabetes mellitus in women. JAMA 289: 1785–1791.FB HuTY LiGA ColditzWC WillettJE Manson2003Television watching and other sedentary behaviors in relation to risk of obesity and type 2 diabetes mellitus in women.JAMA28917851791

4. Blanck HM, McCullough ML, Patel AV, Gillespie C, Calle EE, et al. (2007) Sedentary behavior, recreational physical activity, and 7-year weight gain among postmenopausal U.S. women. Obesity Res 15: 1578–1588.HM BlanckML McCulloughAV PatelC. GillespieEE Calle2007Sedentary behavior, recreational physical activity, and 7-year weight gain among postmenopausal U.S. women.Obesity Res1515781588

5. Fung TT, Hu FB, Yu J, Chu NF, Spiegelman D, et al. (2000) Leisure-time physical activity, television watching, and plasma biomarkers of obesity and cardiovascular disease risk. Am J Epidemiol 152: 1171–1178.TT FungFB HuJ. YuNF ChuD. Spiegelman2000Leisure-time physical activity, television watching, and plasma biomarkers of obesity and cardiovascular disease risk.Am J Epidemiol15211711178

6. Beunza JJ, Martinez-Gonzalez MA, Ebrahim S, Bes-Rastrollo M, Nunez J, et al. (2007) Sedentary behaviors and the risk of incident hypertension - The SUN cohort. Am J Hypert 20: 1156–1162.JJ BeunzaMA Martinez-GonzalezS. EbrahimM. Bes-RastrolloJ. Nunez2007Sedentary behaviors and the risk of incident hypertension - The SUN cohort.Am J Hypert2011561162

7. Helmerhorst HJF, Wijndaele K, Brage S, Wareham NJ, Ekelund U (2009) Objectively measured sedentary time predicts insulin resistance, independent of moderate and vigorous physical activity. Diabetes 58: 1776–1779.HJF HelmerhorstK. WijndaeleS. BrageNJ WarehamU. Ekelund2009Objectively measured sedentary time predicts insulin resistance, independent of moderate and vigorous physical activity.Diabetes5817761779

8. Hu FB, Leitzmann MF, Stampfer MJ, Colditz GA, Willett WC, et al. (2001) Physical activity and television watching in relation to risk for type 2 diabetes mellitus in men. Arch Intern Med 161: 1542–1548.FB HuMF LeitzmannMJ StampferGA ColditzWC Willett2001Physical activity and television watching in relation to risk for type 2 diabetes mellitus in men.Arch Intern Med16115421548

9. Hamilton MT, Healy GN, Dunstan DW, Zderic TW, Owen N (2008) Too little exercise and too much sitting: Inactivity physiology and the need for new recommendations on sedentary behaviour. Curr Cardiov Risk Rep 2: 292–298.MT HamiltonGN HealyDW DunstanTW ZdericN. Owen2008Too little exercise and too much sitting: Inactivity physiology and the need for new recommendations on sedentary behaviour.Curr Cardiov Risk Rep2292298

10. Katzmarzyk PT, Church T, Craig CL, Bouchard C (2009) Sitting time and mortality from all causes, cardiovascular disease and cancer. Med Sci Sports Exerc 41: 998–1005.PT KatzmarzykT. ChurchCL CraigC. Bouchard2009Sitting time and mortality from all causes, cardiovascular disease and cancer.Med Sci Sports Exerc419981005

11. Patel AV, Bernstein L, Deka A, Feigelson HS, Campbell PT, et al. (2010) Leisure time spent sitting in relation to total mortality in a prospective cohort of US adults. Am J Epidemiol 172: 419–429.AV PatelL. BernsteinA. DekaHS FeigelsonPT Campbell2010Leisure time spent sitting in relation to total mortality in a prospective cohort of US adults.Am J Epidemiol172419429

12. Manson JE, Greenland P, LaCroix AZ, Stefanick ML, Mouton CP, et al. (2002) Walking compared with vigorous exercise for the prevention of cardiovascular events in women. New Engl J Med 347: 716–725.JE MansonP. GreenlandAZ LaCroixML StefanickCP Mouton2002Walking compared with vigorous exercise for the prevention of cardiovascular events in women.New Engl J Med347716725

13. Dunstan DW, Barr ELM, Healy GN, Salmon J, Shaw JE, et al. (2010) Television viewing time and mortality: the Ausdiab Study. Circulation 121: 384–391.DW DunstanELM BarrGN HealyJ. SalmonJE Shaw2010Television viewing time and mortality: the Ausdiab Study.Circulation121384391

14. Wijndaele K, Brage S, Besson H, Khaw KT, Sharp SJ, et al. (2010) Television viewing time independently predicts all-cause and cardiovascular mortality: the EPIC Norfolk Study. Int J Epidemiol. K. WijndaeleS. BrageH. BessonKT KhawSJ Sharp2010Television viewing time independently predicts all-cause and cardiovascular mortality: the EPIC Norfolk Study.Int J EpidemiolAdvance access: June 23 2010 (DOI:10.1093/ije/dyq105). Advance access: June 23 2010 (DOI:10.1093/ije/dyq105).

15. Warren TY, Barry V, Hooker SP, Sui X, Church TS, et al. (2010) Sedentary behaviors increase risk of cardiovascular disease mortality in men. Med Sci Sports Exerc 42: 879–885.TY WarrenV. BarrySP HookerX. SuiTS Church2010Sedentary behaviors increase risk of cardiovascular disease mortality in men.Med Sci Sports Exerc42879885

16. Stamatakis E, Hamer M, Dunstan DW (2011) Sceen-based entertainment time, all-cause mortality, and cardiovascular events. Population-based study with ongoing mortality and hospital events follow-up. J Am Coll Cardiol 57: 282–299.E. StamatakisM. HamerDW Dunstan2011Sceen-based entertainment time, all-cause mortality, and cardiovascular events. Population-based study with ongoing mortality and hospital events follow-up.J Am Coll Cardiol57282299

17. Leal J, Luengo-Fernandez R, Gray A, Petersen S, Rayner M (2006) Economic burden of cardiovascular diseases in the enlarged European Union. Eur Heart J 27: 1610–1619.J. LealR. Luengo-FernandezA. GrayS. PetersenM. Rayner2006Economic burden of cardiovascular diseases in the enlarged European Union.Eur Heart J2716101619

18. Tendera M (2006) How much does Europe invest in the treatment of cardiovascular diseases? Eur Heart J 27: 1521–1522.M. Tendera2006How much does Europe invest in the treatment of cardiovascular diseases?Eur Heart J2715211522

19. Office for National Statistics (2006) The time use survey, 2005. How we spend our time. London, UK: Office for National Statistics. Office for National Statistics2006The time use survey, 2005. How we spend our timeLondon, UKOffice for National Statistics

20. Nielsen Media Research (2000) Nielsen report on television. New York, USA: Nielsen Media Research. Nielsen Media Research2000Nielsen report on televisionNew York, USANielsen Media Research

21. Australian Bureau of Statistics (2008) How Australians use their time, 2006. Canberra, Australia: Commomwealth of Australia. Australian Bureau of Statistics2008How Australians use their time, 2006Canberra, AustraliaCommomwealth of Australia

22. Sugiyama T, Healy GN, Dunstan DW, Salmon J, Owen N (2008) Is television viewing time a marker of a broader pattern of sedentary behavior? Ann Behav Med 35: 245–250.T. SugiyamaGN HealyDW DunstanJ. SalmonN. Owen2008Is television viewing time a marker of a broader pattern of sedentary behavior?Ann Behav Med35245250

23. Harris JL, Bargh JA, Brownell KD (2009) Priming effects of television food advertising on eating behavior. Health Psychol 28: 404–413.JL HarrisJA BarghKD Brownell2009Priming effects of television food advertising on eating behavior.Health Psychol28404413

24. Day N, Oakes S, Luben R, Khaw KT, Bingham S, et al. (1999) EPIC-Norfolk: study design and characteristics of the cohort. Br J Cancer 80: 95–103.N. DayS. OakesR. LubenKT KhawS. Bingham1999EPIC-Norfolk: study design and characteristics of the cohort.Br J Cancer8095103

25. Sinha S, Myint PK, Luben RN, Khaw KT (2008) Accuracy of death certification and hospital record linkage for identification of incident stroke. BMC Med Res Meth 8: S. SinhaPK MyintRN LubenKT Khaw2008Accuracy of death certification and hospital record linkage for identification of incident stroke.BMC Med Res Meth8

26. Wareham NJ, Jakes RW, Rennie KL, Mitchell J, Hennings S, et al. (2002) Validity and repeatability of the EPIC-Norfolk Physical Activity Questionnaire. Int J Epidemiol 31: 168–174.NJ WarehamRW JakesKL RennieJ. MitchellS. Hennings2002Validity and repeatability of the EPIC-Norfolk Physical Activity Questionnaire.Int J Epidemiol31168174

27. Besson H, Ekelund U, Brage S, Luben R, Bingham S, et al. (2008) Relationship between subdomains of total physical activity and mortality. Med Sci Sports Exerc 40: 1909–1915.H. BessonU. EkelundS. BrageR. LubenS. Bingham2008Relationship between subdomains of total physical activity and mortality.Med Sci Sports Exerc4019091915

28. Ainsworth BE, Haskell WL, Whitt MC (2000) Compendium of physical activity: an update of activity codes and MET intensities. Med Sci Sports Exerc 32: S498–S516.BE AinsworthWL HaskellMC Whitt2000Compendium of physical activity: an update of activity codes and MET intensities.Med Sci Sports Exerc32S498S516

29. Bingham SA, Gill C, Welch A, Cassidy A, Runswick SA, et al. (1997) Validation of dietary assessment methods in the UK arm of EPIC using weighed records, and 24-hour urinary nitrogen and potassium and serum vitamin C and carotenoids as biomarkers. Int J Epidemiol 26: S137–S151.SA BinghamC. GillA. WelchA. CassidySA Runswick1997Validation of dietary assessment methods in the UK arm of EPIC using weighed records, and 24-hour urinary nitrogen and potassium and serum vitamin C and carotenoids as biomarkers.Int J Epidemiol26S137S151

30. Brage S, Wedderkopp N, Ekelund U, Franks PW, Wareham NJ, et al. (2004) Features of the metabolic syndrome are associated with objectively measured physical activity and fitness in Danish children - the European Youth Heart Study (EYHS). Diabetes Care 27: 2141–2148.S. BrageN. WedderkoppU. EkelundPW FranksNJ Wareham2004Features of the metabolic syndrome are associated with objectively measured physical activity and fitness in Danish children - the European Youth Heart Study (EYHS).Diabetes Care2721412148

31. Clark BK, Healy GN, Winkler EAH, Gardiner PA, Sugiyama T, et al. (2010) Relationship of television time with accelerometer-derived sedentary time: NHANES. BK ClarkGN HealyEAH WinklerPA GardinerT. Sugiyama2010Relationship of television time with accelerometer-derived sedentary time: NHANES.Advance access: October 11 2010 (DOI:10.1249/MSS.0b013e3182019510). Advance access: October 11 2010 (DOI:10.1249/MSS.0b013e3182019510).