You are here

Article Text

Article menu

Download PDFPDF

Coronary artery disease
Original research article
Peak flow rate and death due to coronary heart disease: 30-year results from the Northwick Park Heart cohort study
  1. Tim C Clayton1,
  2. Tom W Meade2,
  3. Elizabeth L Turner3,4 and
  4. Bianca L De Stavola1
  1. 1Department of Medical Statistics, London School of Hygiene and Tropical Medicine, London, UK
  2. 2Department of Non-Communicable Disease Epidemiology, London School of Hygiene and Tropical Medicine, London, UK
  3. 3Department of Biostatistics and Bioinformatics, Duke University School of Medicine, Durham, USA
  4. 4Duke Global Health Institute, Duke University, Durham, USA
  1. Correspondence to Mr Tim C Clayton; tim.clayton{at}lshtm.ac.uk

Abstract

Objective Numerous studies have reported that chronic obstructive pulmonary disease or impaired lung function are associated with later coronary heart disease (CHD). However, it is unclear if lung function is an independent risk factor, as many of these studies have included only limited measures of other factors associated with CHD.

Methods In total 2167 men of all ages in the first Northwick Park Heart Study were followed for a median of 30 years. Cox regression models were used to assess the relationship between peak flow rate (PFR) and CHD mortality adjusted for potential confounders measured at baseline. Analyses allowed for missing data, and secondary analyses for repeat measures on some men and competing risks of CHD death.

Results There were 254 CHD deaths with some evidence of an association between PFR and CHD mortality. The adjusted HRs (95% CIs) from the lowest to the highest of four PFR quartiles were 1.53 (1.04 to 2.25), <430 L/min; 1.43 (0.99 to 2.08), 430 – <490 L/min; and 1.31 (0.93 to 1.86), 490 – <550 L/min; compared with the reference group of ≥550 L/min (trend test p=0.025). Other associations with CHD mortality were observed for systolic blood pressure (p<0.0001), body mass index (p=0.0002), smoking status (p=0.015), blood cholesterol (p=0.005), plasma fibrinogen (p=0.001) and high-risk ECG (p=0.021). There were no strong associations for factors V and VIII or platelet count.

Conclusions After allowing for a range of other risk factors associated with CHD, there was only limited evidence of a relation between PFR and CHD mortality.

This is an Open Access article distributed in accordance with the terms of the Creative Commons Attribution (CC BY 4.0) license, which permits others to distribute, remix, adapt and build upon this work, for commercial use, provided the original work is properly cited. See: http://creativecommons.org/licenses/by/4.0/

https://doi.org/10.1136/openhrt-2014-000164

Statistics from Altmetric.com

    Request Permissions

    If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

    KEY MESSAGES

    What is already known about this subject?

    • Impaired lung function is thought to be a risk factor for coronary heart disease.

    What does this study add?

    • While there is evidence of an association between respiratory function and coronary heart disease it is unlikely to be a strong risk factor.

    How might this impact on clinical practice?

    • It is important to deal with impaired respiratory function but it is unlikely to impact strongly on the risk of coronary heart disease to the extent suggested by other studies which indicated a doubling in the risk whereas our study, after adjusting for a range of risk factors, indicates an increased risk of approximately 40%.

    Introduction

    There have been many reports that chronic obstructive pulmonary disease (COPD) or measures found in impaired pulmonary function such as peak flow rate (PFR) or forced expiratory volume (FEV) are associated with coronary heart disease (CHD) events.1–16 Consequently, impaired respiratory function is generally considered as a risk factor for CHD. Results similar to those for CHD have been reported for stroke.17 ,18 Associations between respiratory function with the development of carotid atherosclerotic plaques19 and with central arterial stiffness have also been described. This report uses PFR as an index of respiratory function to assess the relationship with CHD mortality. The study's first advantage is that the results presented include data on some variables in the coagulation system influencing the thrombotic component of CHD that have not previously been reported in relation to PFR and COPD. Second, analyses of results have included methods to allow for imputation of missing values and also consideration of competing risks of non-CHD deaths. Third, repeat measurements have enabled assessments of the effects of changes over time using time-updated modelling techniques. Finally, participants have been followed up for a median duration of 30 years, so that the value or otherwise of data on PFR and CHD mortality over a longer time period than in other studies has been possible.

    Methods

    This report is based on data from the first Northwick Park Heart Study (NPHS-1). Full details have been reported elsewhere.20

    Participants

    Between 1972 and 1978, NPHS-I recruited men of all ages working at three industrial or occupational sites in North West London. The response to invitations to participate was about 80%. Methods for the clinical examination, laboratory measurements and clotting factor assays and follow-up have been described before.20 Follow-up for non-fatal events ended in 1986. However, follow-up for fatal events, with which this report is concerned, continued until June 2006, with notifications of deaths and their causes having been provided throughout the study by the Office for National Statistics. Blind assessments of causes of death have been made as before,20 based on information from general practitioners, hospitals and coroners.

    As in previous reports of men in the NPHS-1 study, analyses are confined to white participants who had not previously had myocardial infarcts (but including 78 men with a history of definite or possible angina).20 Most men (1703/2167; 78%) had repeat blood samples taken at follow-up examinations between 1978 and 1984, at a median of 6.5 years after entry to the study. Comparison of baseline values for those who were or were not available for follow-up (excluding the 81 men who died within 6.5 years of entry to the study in order to include only those potentially available for follow-up) showed the former to be on average 4 years older, to use less alcohol and to have higher mean body mass index (BMI), blood cholesterol, red blood cell count, haemoglobin and packed cell volumes than those not followed up. Those not followed had lower age-adjusted all-cause mortality and similar CHD mortality compared with those followed up. Bearing these differences in mind, it has been possible to establish prospective data on the impact of PFR on cardiac mortality over time, using time-updated models to take account of the repeat measurements made.

    Measurements

    With participants standing up, PFR was measured three times using the original Wright peak flow metre. The clinic nurse gave instructions on what was required, with a demonstration of the appropriate respiratory effort. With a brief rest interval between each of the three measures, the highest of these was recorded. Besides PFR itself and age and gender, characteristics recorded20 included smoking status, alcohol use, height and weight and hence BMI, blood pressure (average of three readings), social class, total cholesterol and triglycerides (fasting), and automated Coulter full blood count (details shown or noted in tables). Biological activities of clotting factors V, VII and VIII were assayed,21–23 and fibrinogen was measured by a clot weight method.24 Of the clotting factors, only results for factor VII and fibrinogen are shown, as those for V and VIII were not associated with CHD, although they were included in the statistical models to allow as fully as possible for any confounding effect on the relationship between PFR and CHD. ECG findings were coded according to the Minnesota code,25 and classified as indicating higher or lower risk of CHD according to the groupings established, and successfully validated, for the Medical Research Council (MRC) Hypertension Trial.26

    NPHS-1 started in 1972, which was before the development of ethics committees in the UK. It was therefore not submitted to ethics committees. Those approached about the study were given a full verbal explanation of the reasons for it, its nature and what it would involve, and agreed to take part.

    Statistical analyses

    PFR was considered in four categories: (1) <430 L/min (2) 430–<490 L/min (3) 490–<550 L/min and (4) ≥550 L/min. These categories were defined to include approximately a quarter of the CHD deaths in each. PFR was categorised in this way as there was no observed association between PFR and CHD mortality at higher levels of PFR. The characteristics of the patients at baseline and also for those who were followed up have been tabulated, and CHD death rates calculated. The primary analysis has focused on the association of PFR measured at baseline with CHD death, using Cox regression models. Unadjusted HRs for PFR were estimated in order to compare with adjusted HRs, and Kaplan-Meier curves showing the cumulative percentage of CHD death by PFR groups have been presented. All other analyses have been adjusted for age by setting the timescale in the Cox models to be current age, while also including age at baseline to allow that baseline characteristics may have been measured several years previously.27 A Cox regression model was developed for the primary analysis to allow for the confounding effect of other potential risk factors collected at baseline. Covariates of interest were included irrespective of the strength of evidence for their association with CHD mortality to more fully allow for any potential confounding due to these covariates. Data were missing for some variables, particularly for some of the laboratory measurements, although the numbers were low for baseline data (<5% for any one variable), being higher for follow-up data (<12%). In order to adjust for these missing data, multiple imputation techniques with chained equations were used, assuming missingness is at random, with 100 imputations using all covariates in the model including the outcome variable of CHD mortality and time in the study.28 Analyses were also conducted using a complete-case analysis (ie, those with values for all variables) and by imputing the median value for missing covariates to assess whether results were consistent with the findings from the primary analysis.

    In order to assess whether the impact of PFR changed according to time in the study or with age, multivariable Cox models were fitted with interaction terms between PFR and (1) time in follow-up (≤15 or >15 years) or (2) current age (≤75 or >75 years) after first splitting the data according to the relevant timescales.

    A secondary analysis was also conducted using information for those in whom a follow-up visit was conducted (acknowledging, as already indicated, that those followed up may not have been truly representative of the initial population studied). Time-updated Cox models accounting for repeat measurements of variables for those in whom a follow-up visit was conducted were used to consider age-adjusted and more fully adjusted HRs.

    Finally, an analysis was conducted to consider whether there was any impact of competing risks due to deaths from non-cardiac causes using the method of Fine and Gray29 for competing risks regression.

    All analyses were conducted using Stata V.13.0 (StataCorp. 2013. Stata Statistical Software: Release 13. College Station, Texas, USA: StataCorp LP).

    Results

    In total, data for 2167 men were included in the analyses, and the men were followed for a median of 30 years. Table 1 shows that there were 254 cardiac deaths (rate 4.58/1000 person-years). In addition, 720 men died of non-cardiac causes over the follow-up period. Table 1 also shows that, as expected, advancing age is strongly associated with increasing HRs for CHD death. The unadjusted HRs increase with decreasing PFR (figure 1), although after adjustment for age the HRs move towards 1 (table 2). Table 2 gives results from the multivariable model and shows that after adjusting for age and other risk factors measured at baseline, including BMI and smoking, the impact of PFR on CHD is reduced further. The hazard for the lowest PFR quartile is 54% more than for the highest quartile, although the evidence for an association, though suggestive, is not strong (trend test p=0.025). There was an increased risk of CHD mortality with smoking, BMI, systolic blood pressure, cholesterol, fibrinogen and an abnormal ECG considered high risk by the Minnesota coding. NPHS-1 has not previously included results based on ECG data, and our results now show that ‘high-risk’ abnormalities are associated with increased CHD mortality, HR 1.58 (95% CI 1.11 to 2.27). The association of factor VII with CHD mortality is of marginal statistical significance (p=0.058). Other factors included in this analysis but with little evidence for an association with cardiac death were triglycerides, haemoglobin, red blood cell count, white blood cell count or packed cell volume (data not shown though referred to in table 2). Neither alcohol use nor social class appeared to be independent risk factors. There was no evidence that the association between PFR and CHD mortality changed over time of follow-up (interaction p value 0.79) or with age (interaction p=0.24).

    View this table:
    Table 1

    Characteristics and cardiac mortality of individuals at baseline and follow-up

    View this table:
    Table 2

    Multivariable associations with cardiac mortality using baseline characteristics

    Figure 1
    Figure 1

    Kaplan-Meier curve of cumulative percentage of cardiac mortality by peak flow rate category (adjusted HRs (95% CI) compared with ≥550 L/min: (<430 L/min) 1.53 (1.04 to 2.25); (430 – <490 L/min) 1.43 (0.99 to 2.08); (490 – <550 L/min) 1.31 (0.93 to 1.86)).

    In the secondary analysis using updated values for those in whom a follow-up visit was available, table 3 shows that the association of PFR with CHD mortality was less than for those with baseline values (trend test p=0.50), though the hazard was 15% higher in the lowest quartile of PFR than in the highest quartile (≥550 L/min). There was no evidence for an association of factor VII with CHD mortality, whereas the impact of a ‘high-risk’ ECG abnormality was stronger than when considering baseline ECG alone.

    View this table:
    Table 3

    Multivariable associations with cardiac mortality using time-updated model

    Further analyses considered complete-case analyses only, and also by imputation of the median values of missing variables. Results were largely consistent with those shown. Finally, an analysis was undertaken allowing for the potential competing risk of non-cardiac mortality and results were again consistent with the results presented (results not shown).

    Discussion

    The prospective studies that have included measures of impaired lung function such as PFR, FEV or FEV1 and that have shown an association with the incidence of fatal CHD have come from many settings, often the USA,4 ,5 ,7 ,10 but also from the UK and other European countries,2 ,6 ,8 ,12 ,14 China,11 ,15 Australia and New Zealand,9 and India.13 Some have collected data from more than one country.1 ,3 Many studies have typically reported relative risks (RRs) of 2 or less for the association between measures of impaired lung function and CHD, though in two North American studies,7 risks up to 4 were recorded. A systematic review of the relationship between FEV1 and cardiovascular mortality reported a pooled RR of 1.75 comparing those in the lowest quintile to the highest quintile of FEV1.16 In NPHS-1, the apparent effect of PFR on CHD, showing a HR difference of 53% between the lowest and highest quartiles, was modest using baseline characteristics (ie, analysis based on all participants) and had largely disappeared based on those who had follow-up examinations.

    There do not appear to be large or consistent differences between FEV, FEV1 and PFR that would affect the findings of epidemiological studies. A study of older participants in east Boston5 also used only PFR, and gave much the same results as ours and as those for other studies (differences may be important in assessing lung function in clinical settings, such as acute asthma attacks in children.)

    Data for men who had previously had major CHD events were not included in the analyses. However, abnormal ECGs often reflecting early pathological changes of CHD were, not surprisingly, significantly associated with CHD mortality. Omission of ECG data from other studies may thus well account for the generally stronger associations of PFR and other measures of impaired lung function with CHD than in NPHS–1, since these other studies have not been in a position to allow for the early pathology of CHD and the correlation of PFR with ECG findings. Smoking, cholesterol and systolic blood pressure were other independent risk factors associated with CHD death. We are aware of only one other study that has included measurement of plasma fibrinogen, which increases with the level of smoking but remains a strong independent risk factor for CHD. However, this study,30 CARDIA, dealt mainly with the association of fibrinogen with risk factors including COPD that contribute to CHD, rather than prospective associations with CHD itself.

    In the present study, the association of plasma fibrinogen with fatal CHD was highly statistically significant in the multivariable analysis (table 2), which takes account of correlations between different risk factors. An explanation might be the contribution of fibrinogen to thrombosis in CHD, or as a marker of the degree of an inflammatory process in COPD and CHD. The trend was similar but reduced in the time-updated model (table 3), although there was little evidence for an association. There was a suggestion that factor VII may be associated with CHD in the multivariable model using baseline data, and may thus have contributed to the thrombotic element in CHD, though the finding was at a marginal level and there was no association in the time-updated model. There was no evidence that factors V and VIII and platelet count (as other components of the haemostatic system) were associated with fatal CHD. The Whitehall Cohort II study concluded that socioeconomic differences were not important in explaining the association between lung function and mortality.12

    Virtually all prospective studies are missing data for some variables. Long-term follow-up periods almost inevitably mean that their data will be subject to competing risks of non-cardiac mortality. However, few allow for these points. We have used well-tested and accepted methods of imputation to allow for missing values in our data, and the results were similar compared with the analyses not using these methods. Allowance for competing risks of non-CHD mortality also made little difference.

    Treatment of COPD appears to reduce the onset of cardiovascular disease.31 However, Lange et al6 concluded that “although impaired ventilatory function is a significant predictor of death from myocardial infarction and other cardiovascular diseases, it should not be regarded as a genuine risk factor for ischaemic heart disease.” We would agree with this view, in that the inclusion of a wide range of characteristics associated with CHD, including fibrinogen and abnormal ECG findings, suggests that measures such as PFR are not strongly and independently associated with fatal CHD, and that the view that they usefully predict CHD should be modified accordingly.

    References

    1. ↵
      1. Friedman GD,
      2. Klatsky AL,
      3. Siegelaub AB
      . Lung function and risk of myocardial infarction and sudden cardiac death. N Engl J Med 1976;294:1071–5.
      OpenUrlCrossRefPubMedWeb of Science
    2. ↵
      1. Menotti A,
      2. Mariotti S,
      3. Seccareccia S,
      4. et al
      . The 25 year estimated probability of death from some specific causes as a function of twelve risk factors in middle aged men. Eur J Epidemiol 1988;4:60–7.
      OpenUrlCrossRefPubMedWeb of Science
    3. ↵
      1. Ebi-Kryston KL,
      2. Hawthorne VM,
      3. Rose G,
      4. et al
      . Breathlessness, chronic bronchitis and reduced pulmonary function as predictors of cardiovascular disease mortality among men in England, Scotland and the United States. Int J Epidemiol 1989;18:84–8.
      OpenUrlAbstract/FREE Full Text
    4. ↵
      1. Marcus EB,
      2. Curb JD,
      3. MacLean CJ,
      4. et al
      . Pulmonary function as a predictor of coronary heart disease. Am J Epidemiol 1989;129:97–104.
      OpenUrlAbstract/FREE Full Text
    5. ↵
      1. Cook NR,
      2. Evans DA,
      3. Scherr PA,
      4. et al
      . Peak expiratory flow rate and 5-year mortality in an elderly population. Am J Epidemiol 1991;133:784–94.
      OpenUrlAbstract/FREE Full Text
    6. ↵
      1. Lange P,
      2. Nyboe J,
      3. Jensen G,
      4. et al
      . Ventilatory function impairment and risk of cardiovascular death and of fatal or non-fatal myocardial infarction. Eur Respir J 1991;4:1080–7.
      OpenUrlAbstract/FREE Full Text
    7. ↵
      1. Schroeder EB,
      2. Welch VL,
      3. Couper D,
      4. et al
      . Lung function and incident coronary heart disease: the Atherosclerosis Risk in Communities Study. Am J Epidemiol 2003;158:1171–81.
      OpenUrlAbstract/FREE Full Text
    8. ↵
      1. Engstrom G,
      2. Hedblad B,
      3. Janzon L
      . Reduced lung function predicts increased fatality in future cardiac events. A population-based study. J Intern Med 2006;260:560–7.
      OpenUrlCrossRefPubMedWeb of Science
    9. ↵
      1. Young RP,
      2. Hopkins R,
      3. Eaton TE
      . Forced expiratory volume in one second: not just a lung function test but a marker of premature death from all causes. Eur Respir J 2007;30:616–22.
      OpenUrlAbstract/FREE Full Text
    10. ↵
      1. Lee HM,
      2. Le H,
      3. Lee BT,
      4. et al
      . Forced vital capacity paired with Framingham Risk Score for prediction of all-cause mortality. Eur Respir J 2010;36:1002–6.
      OpenUrlAbstract/FREE Full Text
    11. ↵
      1. Yang L,
      2. Zhou M,
      3. Smith M,
      4. et al
      . Body mass index and chronic obstructive pulmonary disease-related mortality: a nationally representative prospective study of 220,000 men in China. Int J Epidemiol 2010;39:1027–36.
      OpenUrlAbstract/FREE Full Text
    12. ↵
      1. Sabia S,
      2. Shipley M,
      3. Elbaz A,
      4. et al
      . Why does lung function predict mortality? Results from the Whitehall II Cohort Study. Am J Epidemiol 2010;172:1415–23.
      OpenUrlAbstract/FREE Full Text
    13. ↵
      1. Hebert JR,
      2. Pednekar MS,
      3. Gupta PC
      . Forced expiratory volume predicts all-cause and cancer mortality in Mumbai, India: results from a population-based cohort study. Int J Epidemiol 2010;39:1619–27.
      OpenUrlAbstract/FREE Full Text
    14. ↵
      1. Wannamethee SG,
      2. Shaper AG,
      3. Rumley A,
      4. et al
      . Lung function and risk of type 2 diabetes and fatal and nonfatal major coronary heart disease events: possible associations with inflammation. Diabetes Care 2010;33:1990–6.
      OpenUrlAbstract/FREE Full Text
    15. ↵
      1. Smith M,
      2. Zhou M,
      3. Wang L,
      4. et al
      . Peak flow as a predictor of cause-specific mortality in China: results from a 15-year prospective study of ∼170,000 men. Int J Epidemiol 2013;42:803–15.
      OpenUrlAbstract/FREE Full Text
    16. ↵
      1. Sin DD,
      2. Wu L,
      3. Man SF
      . The relationship between reduced lung function and cardiovascular mortality: a population-based study and a systematic review of the literature. Chest 2005;127:1952–9.
      OpenUrlCrossRefPubMedWeb of Science
    17. ↵
      1. Wannamethee SG,
      2. Shaper AG,
      3. Ebrahim S
      . Respiratory function and risk of stroke. Stroke 1995;26:2004–10.
      OpenUrlAbstract/FREE Full Text
    18. ↵
      1. Hozawa A,
      2. Billings JL,
      3. Shahar E,
      4. et al
      . Lung function and ischemic stroke incidence: the Atherosclerosis Risk in Communities study. Chest 2006;130:1642–9.
      OpenUrlCrossRefPubMedWeb of Science
    19. ↵
      1. Zureik M,
      2. Kauffmann F,
      3. Touboul PJ,
      4. et al
      . Association between peak expiratory flow and the development of carotid atherosclerotic plaques. Arch Intern Med 2001;161:1669–76.
      OpenUrlCrossRefPubMedWeb of Science
    20. ↵
      1. Meade TW,
      2. Brozovic M,
      3. Chakrabarti RR,
      4. et al
      . Haemostatic function and ischaemic heart disease: principal results of the Northwick Park Heart Study. Lancet 1986;2:533–7.
      OpenUrlPubMedWeb of Science
    21. ↵
      1. Brozovic M,
      2. Stirling Y,
      3. Harricks C,
      4. et al
      . Factor VII in an industrial population. Br J Haematol 1974;28:381–91.
      OpenUrlCrossRefPubMedWeb of Science
    22. ↵
      1. Brozovic M,
      2. Chakrabarti R,
      3. Stirling Y,
      4. et al
      . Factor V in an industrial population. Br J Haematol 1976;33:543–50.
      OpenUrlCrossRefPubMed
    23. ↵
      1. Meade TW,
      2. Cooper JA,
      3. Stirling Y,
      4. et al
      . Factor VIII, ABO blood group and the incidence of ischaemic heart disease. Br J Haematol 1994;88:601–7.
      OpenUrlCrossRefPubMedWeb of Science
    24. ↵
      1. Fearnley GR
      . Measurement of spontaneous fibrinolytic activity. J Clin Pathol 1964;17:307–9.
      OpenUrlFREE Full Text
    25. ↵
      1. Rose GA,
      2. Blackburn H
      . Cardiovascular survey methods. East Afr Med J 1969;46:220–7.
      OpenUrlPubMed
    26. ↵
      1. Miall WE,
      2. Greenberg G
      . Mild hypertension—is there pressure to treat? An account of the MRC trial: Press Syndicate of the University of Cambridge, 1987.
    27. ↵
      1. Sperrin M,
      2. Buchan I
      . Modelling time to event with observations made at arbitrary times. Stat Med 2013;32:99–109.
      OpenUrlCrossRefPubMed
    28. ↵
      1. White IR,
      2. Royston P,
      3. Wood AM
      . Multiple imputation using chained equations: issues and guidance for practice. Stat Med 2011;30:377–99.
      OpenUrlCrossRefPubMed
    29. ↵
      1. Fine JP,
      2. Gray RJ
      . A proportional hazards model for the subdistribution of a competing risk. J Am Stat Assoc 1999;94:496.
      OpenUrlCrossRefWeb of Science
    30. ↵
      1. Thyagarajan B,
      2. Jacobs DR,
      3. Apostol GG,
      4. et al
      . Plasma fibrinogen and lung function: the CARDIA Study. Int J Epidemiol 2006;35:1001–8.
      OpenUrlAbstract/FREE Full Text
    31. ↵
      1. Stone IS,
      2. Barnes NC,
      3. Petersen SE
      . Chronic obstructive pulmonary disease: a modifiable risk factor for cardiovascular disease? Heart 2012;98:1055–62.
      OpenUrlAbstract/FREE Full Text

    Footnotes

    • The work was carried out at the London School of Hygiene and Tropical Medicine.

    • Contributors TWM initiated the NPHS-1 and proposed the analyses of the associations of PFR with death due to CHD. TWM also wrote the paper along with TCC. TCC also carried out the main statistical analyses assisted by ELT, using files previously prepared by BLDS. BLDS also provided advice on the data and statistical analyses.

    • Funding The NPHS-1 Study was funded by the UK Medical Research Council. When TWM moved to the London School of Hygiene and Tropical Medicine (LSHTM) on his retirement from the MRC, he transferred the NPHS data to LSHTM, and the other authors have contributed unpaid time from their appointments at LSHTM.

    • Competing interests TCC, as principal author, has had full access to the data and has been responsible for the decision to submit for publication.

    • Ethics approval NPHS-1 started in 1972, which was before the development of ethics committees in the UK. It was therefore not submitted to ethics committees. Those approached about the study were given a full verbal explanation of the reasons for it, its nature and what it would involve, and agreed to take part.

    • Provenance and peer review Not commissioned; externally peer reviewed.

    • Data sharing statement Data from NPHS-1 has been published extensively and this provides the long-term outcomes from this study.