The educational gradient in cardiovascular risk factors: impact of shared family factors in 228,346 Norwegian siblings
BMC Public Health volume 17, Article number: 281 (2017)
Various indicators of childhood socioeconomic position have been related to cardiovascular disease (CVD) risk in adulthood. We investigated the impact of shared family factors on the educational gradient in midlife CVD risk factors by assessing within sibling similarities in the gradient using a discordant sibling design.
Norwegian health survey data (1980–2003) was linked to educational and generational data. Participants with a full sibling in the health surveys (228,346 individuals in 98,046 sibships) were included. Associations between attained educational level (7–9 years, 10–11 years, 12 years, 13–16 years, or >16 years) and CVD risk factor levels in the study population was compared with the corresponding associations within siblings.
Educational gradients in risk factors were attenuated when factors shared by siblings was taken into account: A one category lower educational level was associated with 0.7 (95% confidence interval 0.6 to 0.8) mm Hg higher systolic blood pressure (27% attenuation), 0.4 (0.4 to 0.5) mmHg higher diastolic blood pressure (30%), 1.0 (1.0 to 1.1) more beats per minute higher heart rate (21%), 0.07 (0.06 to 0.07) mmol/l higher serum total cholesterol (32%), 0.2 (0.2 to 0.2) higher smoking level (5 categories) (30%), 0.15 (0.13 to 0.17) kg/m2 higher BMI (43%), and 0.2 (0.2 to 0.2) cm lower height (52%). Attenuation increased with shorter age-difference between siblings.
About one third of the educational gradients in modifiable CVD risk factors may be explained by factors that siblings share. This implies that childhood environment is important for the prevention of CVD.
In high income countries, cardiovascular disease (CVD) risk factors are generally more prevalent among adults with lower socioeconomic position (SEP) [1, 2]. A substantial proportion of the inverse associations between SEP and CVD are driven by the higher prevalence of CVD risk factors by lower SEP [3, 4]. This socioeconomic gradient in CVD represent a potential for disease prevention, and it is of great policy interest to determine which periods during the life course are important for development of the cardiovascular disease risk factors; elevated blood pressure, disadvantageous serum lipid profile, obesity, tobacco use and physical inactivity [5, 6].
Childhood SEP, usually indicated by parental level of education or occupational class, is associated with adult CVD risk factors in several populations [6–10]. Proposed mechanisms for this association include both biological factors such as poor maternal nutrition and health, intrauterine growth retardation, poor growth in early childhood, obesity in adolescence and repeated childhood infections, and include also environmental factors such as health behaviour and psychosocial factors in the family and surroundings [5, 6, 10]. CVD risk factors might track from childhood to adulthood , but are also modifiable in adult age. Body height, in contrast, is a stable trait through adult life. In this context, sibling similarity in the educational gradient in height might serve as a strong indicator of shared family environment.
Indicators of childhood and adulthood SEP are likely to be highly correlated, and it is difficult to assess their independent effect on CVD risk [12, 13]. An alternative and potentially robust approach is to use a sibling comparison, where a model of the association between attained education and CVD risk factors within siblings can be compared with a cohort model not using the sibling approach [10, 12]. Although this approach is not directly comparable with a randomized controlled trial, it makes it possible to ask the question; what would the socioeconomic gradient in cardiovascular risk factors be if everyone experienced similar family environment in childhood. A sibling comparison will control for a number of family factors (genetic and environmental) as a way to quantify the impact of early life on associations between exposures and outcomes in adulthood [10, 12, 14–16]. Studies using this method have reported that part of the educational gradient in adult all-cause and CVD mortality was explained by factors shared by siblings [14–16].
Except for one study by Lawlor et al. analysing the within sibling-effect of education on adult body mass index (BMI) , no other studies have to our knowledge investigated the within sibling-effect of education on CVD risk factors. According to that study, the major part of the inverse association between educational attainment and adult BMI may be explained by family factors shared by siblings .
We investigated the educational gradient in CVD risk factors between siblings in a large sample from the Norwegian general population, being in their early 40s. We aimed to assess the impact of factors in early life shared by siblings on educational gradient in adulthood CVD risk factors (blood pressure, total cholesterol, heart rate, BMI, smoking).
Population based health survey participants from the Counties studies (1974-88) , the Age 40 Program (1985–1999) and CONOR (1994–2003) [19, 20], born after 1940, turning 40–45 years at the year of their health survey, and who had at least one full sibling among the health survey participants were selected for this sibling design. Participants with missing information for either/both of the parents, missing educational data and participants without any full siblings in the health surveys were excluded (Fig. 1).
Health survey data were linked to national educational data, to the Norwegian Population Registry, and to the Norwegian Family Based Life Course (NFLC) study  using the unique national personal identification number.
The index person’s mother, father and siblings were identified from the NFLC study . Parental identification has been proven to be reliable for index persons born from 1940 and onwards . Full siblings were defined as persons registered with the same mother and father in the NFLC study. Twins were included.
Education was registered in the National Educational Database and reported in National Population and Housing Censuses every 10th year from 1970–2001. A person’s highest attained educational level was classified as up to 9 years, [or 7 years in the 1960s], 10-11 years, 12 years (or vocational education with corresponding academic level), 13–16 years, and more than 16 years. Participants with no registered education were excluded.
Cardiovascular risk factors
In all health surveys self-assessed questionnaires, clinical measures and non-fasting blood sampling were collected [18–20]. Smoking status and cigarette pack years were collapsed into a graded variable; 1) never smoker, 2) past smoker with < 20 pack years, 3) past smoker with ≥ 20 pack years, 4) current smoker with < 20 pack years, 5) current smoker with ≥ 20 pack years. Self-reported treatment with blood pressure lowering medications was recorded. Blood pressure was initially measured manually by sphygmomanometers (Ercameter, ERKA, Bad Tölz, Germany), and the second of two measurements defined systolic and diastolic blood pressure. Later, the average of the last two available automatic oscillometric measures (Dinamap, Criticon, Tampa, USA)  defined blood pressure. Heart rate, as a proxy for physical activity , was recorded during the automatic blood pressure measurements which were performed after 2 min of rest . Height and weight were measured and BMI (kg/m2) was calculated. Non-fasting serum total cholesterol were initially measured by non-enzymatic, and later enzymatic method, and the non-enzymatic values were converted by a correction factor . Systematic COronary Risk Evaluation (SCORE) risk prediction score of 10-year CVD mortality was calculated based on our variables; age, sex, systolic blood pressure, total cholesterol and current daily smoking. We applied the algorithm presented in Appendix A in the paper by Conroy et al. and chose the coefficients for CVD risk in a high risk population .
Descriptive statistics are presented as mean ± standard deviation (SD) and as counts of persons (%). The Cuzick test for trend was used. Multilevel linear regressions of the exposure education was fitted in separate models for each of the outcomes; systolic blood pressure, diastolic blood pressure, heart rate, total cholesterol, BMI, height, smoking category and the SCORE risk score. The analyses were adjusted for sex, age at examination and examination year. The cohort estimate assesses the association between all individuals regardless of sibling similarities. In the within sibship estimate, each individual together with their sibling(s) made up a sibship cluster. The individual-specific factors that are constant in siblings (shared genes or familial environment) are omitted in the within sibships estimate. In addition, the within sibships term control for confounding from unobserved family-level factors. The dissimilarity of the association between education and the CVD risk factor for the cohort and the within sibships estimates was tested by using the Hausman specification test .
Evaluation of the dissimilarity between the cohort and within sibship estimates helps interpret the role of shared family factors. If the between- and within-associations are equal, it indicates that the unobserved family-specific factors in the within sibships analyses are not important. A weaker within sibships association than cohort association may indicate that these unmeasured family factors confound parts of the association. A stronger within sibships association may indicate that these unmeasured factors have obscured the association . One example could be; among 4 siblings with different attained educational level, all overweight with a slightly lower BMI for the ones with higher educational level, there is a weaker within sibships educational gradient in BMI than among unrelated individuals in a cohort. If numerous sibships had such weaker gradients, the attenuation of the educational BMI gradient from the cohort analyses to the within sibships analyses would suggest that unmeasured family factors confound the association between educational level and BMI.
Variability expressed as SD in CVD risk factors between individuals in the cohort and within the sibships was assessed. The regression beta coefficient (β) represents the number of units more disadvantageous CVD risk factor per one lower educational level. The percentage change from the βCohort to the βWithin sibships was expressed as: % = ((βCohort – βWithin sibships) / βCohort)*100. In a sub-population of sibships with only two siblings, we stratified the population by age-difference between siblings.
Assumptions were examined in the cohort models by standard linear regression, for which diagnostic statistical tests are available. Potential deviations from linearity or non-constant variance were examined in plots of residuals against predicted values. Robustness was examined by plotting delta beta values against the participants’ reference numbers to look for points with high influence, and by plotting leverage-versus-squared-residual plots. We examined whether the cohort models gave negative predictions.
Sensitivity analyses were performed: We examined interaction between education and sex in the cohort models. We stratified the analyses by sex which meant that only same-sex siblings contributed. Also, the analyses were performed in a sub-population including only sibships with discordant educational level. The statistical package STATA version 13 was used.
Of the 457,142 health survey participants born in 1940 and forth, and examined the year they turned 40–45 years, 91% had information on parental status from the NFLC study, while 39,667 individuals (9%) were excluded due to missing parental status (Fig. 1). The educational level was lower for those with missing parental status (21.7% basic education) than for those with registered parental status (16.3% basic education, p < 0.001). Of the participants with registered parental status 0.4% had unregistered length of education and were excluded from the analyses (Fig. 1). These excluded participants had similar mean percentage points SCORE risk score (0.59 ± 0.62 SD) as the included participants with registered length of education (0.62 ± 0.70), p = 0.496. Of the participants with registered parental status and registered length of education 45% did not have a full sibling in the health surveys and were excluded from the analyses (Fig. 1). For these excluded participants a lower proportion had basic education (15.1%) than for the included participants (17.4%), p < 0.001.
In all, 228,346 participants had a full sibling that had participated in one of the health surveys, and they defined the study population (Fig. 1).
The study population included 117,818 (52%) women and 110,528 (48%) men who were born from 1940–1963 and examined in 1980–2003 at median age 41 years. Only 1.2% of the study cohort had missing values for any of the CVD risk factors.
The study population of 228,346 participants made up 98,046 sibships that included up to 9 siblings. Of the sibships, 73% consisted of two siblings, 20% consisted of three siblings, and 6% consisted of four or more siblings. In all 164, 518 participants (72%) were parts of sibling groups with differing educational level between two or more siblings.
More disadvantageous CVD risk factor levels were found by lower levels of education in men and women (Table 1). Smaller variation and thus more similarity in CVD risk factors were found within sibships than between all individuals in the cohort (Table 2).
Cohort and within sibling analyses
The cohort and within sibships analyses differed significantly by Hausman specification test by p < 0.001 for all CVD risk factors in Table 3. In the within sibships analyses the educational gradients were attenuated; by one third for blood pressure, total cholesterol, smoking category and SCORE risk score, by two fifths for BMI, by one fifth for heart rate, and by one half for the opposite gradient in adult height (Table 3).
In a subsample of sibships with only two siblings, analyses stratified for age-difference within a sibling pair indicated that a relatively larger proportion of the educational gradient in most of the CVD risk factors was explained by sibling similarities between siblings forin the sibling pairs born closer in time than for the sibling pairs with larger age-span (Table 4).
We found no important deviations from linearity, nor any important non-constant variance of residuals, nor any points with undue high influence. Negative predictions were found only for SCORE risk prediction score, and with no more than 7% negative predictions.
Sex-stratified analyses included 73,023 sisters in 33,602 sibships and 65,400 brothers in 30,209 sibships, and showed a similar pattern as the main results. We did not formally test the sex differences in the estimates; however the educational gradients for blood pressure, total cholesterol and BMI in women were numerically stronger than in men (Additional file 1: Tables A and B). Analyses performed in the sub-population of only sibhips discordant for educational attainment showed similar pattern as in the main analyses (Additional file 1: Table C vs. Table 3).
In this study on average one third of the inverse educational gradients in blood pressure, serum total cholesterol, smoking, heart rate and BMI were explained by family factors shared by siblings. Furthermore, these educational gradients were explained to a greater extent with closer age-span between siblings.
The main strength of this study is the large sample size with clinical measurements, sex-stratified analyses, and the consistency of our finding across a wide range of CVD risk factors. The narrow age range of 40-45 years makes the risk factor levels comparable between siblings. The strict definition of siblings (with same mother and father) increases the plausibility of similar childhood environment and proportion of common genes. Also, the divorce rate was below 0.4% per year during 1940–1970 . A potential weakness of the within sibships analysis is that it conditions on sibships discordant for education and CVD risk factors. This implies a selection of sibships that might differ in non-shared causes of the educational level attained and of the measured CVD risk factor level . Non-shared environment include systematic elements such as birth order and birth year that might introduce non-shared confounding . Our results were not altered when adjusting for birth order (results not shown), stratifying the analyses by sex, or restricting the analyses to sibships discordant for educational level (Additional file 1). Non-systematic non-shared environment that we did not take into account (season of birth, sibling-sibling interactions, differential parental treatment and peer groups) might also have contributed to non-shared confounding . Our analyses were strengthened by including adult height, which is a composite measure of genetic disposition, health and nutrition during the growing years , is potentially modifiable during early life up to puberty where most siblings share environment, and is unlikely to be causally affected by attained education unlike the other CVD risk factors examined in this paper. Accordingly, shared family factors had the strongest impact on the educational gradient in height, of all the CVD risk factor gradients (Table 3). Our sibling design has thus the potential to capture unobserved shared confounding, as the educational gradient in height that is not explained by shared family factors might result from residual confounding from non-shared factors in early life, childhood and youth.
Our findings are consistent with the one previous study by Lawlor et al , comparing the associations between educational attainment and BMI within siblings of the same family and between individuals form different families in a smaller study comprising 5 467 individuals. Here the effect of education on adult BMI between individuals from different families was attenuated to the null in within sibling pair analyses.
The sibling similarities in the educational gradients in the modifiable midlife CVD risk factors, blood pressure, serum total cholesterol, smoking, heart rate and BMI, that we identify (Table 3), can be interpreted as a result of both the environment that siblings share; parents’ health behaviour, parents’ socioeconomic position, housing, neighbourhood and schooling, in addition to genes that full siblings share (approximately 50% ). These shared family similarities were weakened by a larger age-span between siblings (Table 4). This is in agreement with siblings with longer age–span sharing the same proportion of genes, but to a lesser extent their childhood environment than sibling pairs born closer in time , and suggests that family environment plays an important role in the socioeconomic gradient in CVD.
Around one third of the association between lower educational attainment and less favourable risk factors seem to be explained by shared family factors. The contribution increases with shorter age-difference between siblings. This suggests that there is substantial scope for prevention of CVD that starts early and that childhood environment matters.
Body mass index
- NFLC study:
The Norwegian family based life course study
Systematic COronary risk evaluation
Beta regression coefficient
Stringhini S, Dugravot A, Shipley M, Goldberg M, Zins M, Kivimaki M, Marmot M, Sabia S, Singh-Manoux A. Health behaviours, socioeconomic status, and mortality: further analyses of the British Whitehall II and the French GAZEL prospective cohorts. PLoS Med. 2011;8(2), e1000419.
Kanjilal S, Gregg EW, Cheng YJ, Zhang P, Nelson DE, Mensah G, Beckles GL. Socioeconomic status and trends in disparities in 4 major risk factors for cardiovascular disease among US adults, 1971-2002. Arch Intern Med. 2006;166(21):2348–55.
Stringhini S, Sabia S, Shipley M, Brunner E, Nabi H, Kivimaki M, Singh-Manoux A. Association of socioeconomic position with health behaviors and mortality. JAMA. 2010;303(12):1159–66.
Laaksonen M, Talala K, Martelin T, Rahkonen O, Roos E, Helakorpi S, Laatikainen T, Prattala R. Health behaviours as explanations for educational level differences in cardiovascular and all-cause mortality: a follow-up of 60 000 men and women over 23 years. Eur J Public Health. 2008;18(1):38–43.
Harper S, Lynch J, Smith GD. Social determinants and the decline of cardiovascular diseases: understanding the links. Annu Rev Public Health. 2011;32:39–69.
Ben-Shlomo Y, Kuhn D. A life course approach to chronic disease epidemiology: conceptual models, empirical challenges and interdisciplinary perspectives. Int J Epidemiol. 2002;31(2):285–93.
Galobardes B, Lynch JW, Smith GD. Is the association between childhood socioeconomic circumstances and cause-specific mortality established? Update of a systematic review. J Epidemiol Community Health. 2008;62(5):387–90.
Lamont D, Parker L, White M, Unwin N, Bennett S, Cohen M, Richardson D, Dickinson HO, Adamson A, Alberti K. Risk of cardiovascular disease measured by carotid intima-media thickness at age 49-51: lifecourse study. BMJ. 2000;320(7230):273–8.
Packard CJ, Bezlyak V, McLean JS, Batty GD, Ford I, Burns H, Cavanagh J, Deans KA, Henderson M, McGinty A. Early life socioeconomic adversity is associated in adult life with chronic inflammation, carotid atherosclerosis, poorer lung function and decreased cognitive performance: a cross-sectional, population-based study. BMC Public Health. 2011;11(1):42.
Smith GD. Epidemiology, epigenetics and the ‘Gloomy Prospect’: embracing randomness in population health research and practice. International J Epidemiol. 2011;40(3):537–62.
Mattsson N, Ronnemaa T, Juonala M, Viikari JS, Raitakari OT. Childhood predictors of the metabolic syndrome in adulthood. The Cardiovascular Risk in Young Finns Study. Ann Med. 2008;40(7):542–52.
Gilman SE, Loucks EB. Invited commentary: does the childhood environment influence the association between every x and every y in adulthood? Am J Epidemiol. 2012;176(8):684–8.
Galobardes B, Smith GD, Lynch JW. Systematic review of the influence of childhood socioeconomic circumstances on risk for cardiovascular disease in adulthood. Ann Epidemiol. 2006;16(2):91–104.
Naess O, Hoff DA, Lawlor D, Mortensen LH. Education and adult cause-specific mortality—examining the impact of family factors shared by 871 367 Norwegian siblings. Int J Epidemiol. 2012;1–9.
Søndergaard G, Mortensen LH, Andersen A-MN, Andersen PK, Dalton SO, Madsen M, Osler M. Does shared family background influence the impact of educational differences on early mortality? Am J Epidemiol. 2012;176(8):675–83.
Madsen M, Andersen PK, Gerster M, Andersen AM, Christensen K, Osler M. Are the educational differences in incidence of cardiovascular disease explained by underlying familial factors? A twin study. Soc Sci Med. 2014;118:182–90.
Lawlor DA, Clark H, Davey Smith G, Leon DA. Childhood intelligence, educational attainment and adult body mass index: findings from a prospective cohort and within sibling-pairs analysis. Int J Obes. 2006;30(12):1758–65.
Bjartveit K, Foss OP, Gjervig T, Lund-Larsen PG. The cardiovascular disease study in Norwegian counties. Background and organization. Acta Med Scand Suppl. 1979;634:1–70.
Naess O, Sogaard AJ, Arnesen E, Beckstrom AC, Bjertness E, Engeland A, Hjort PF, Holmen J, Magnus P, Njolstad I. Cohort profile: cohort of Norway (CONOR). Int J Epidemiol. 2008;37(3):481–5.
Bjartveit K, Stensvold I, Lund-Larsen P, Gjervig T, Krüger O, Urdal P. Cardiovascular screenings in Norwegian counties. Background and implementation. Status of risk pattern during the period 1986-90 among persons aged 40-42 years in 14 counties. Tidsskr Nor Laegeforen. 1991;111(17):2063.
Graff-Iversen S. External female hormones, serum lipids and mortality - Population studies of Norwegian women. Oslo: University of Oslo; 2005.
Thelle DS, Selmer R, Gjesdal K, Sakshaug S, Jugessur A, Graff-Iversen S, Tverdal A, Nystad W. Resting heart rate and physical activity as risk factors for lone atrial fibrillation: a prospective study of 309 540 men and women. Heart. 2013;99(23):1755–60.
Naess O, Hoff DA. The Norwegian Family Based Life Course (NFLC) study: data structure and potential for public health research. Int J Public Health. 2013;58(1):57–64.
Lund-Larsen P. Blood pressure measured with a sphygmomanometer and with Dinamap under field conditions–a comparison. Nor Epidemiol. 1997;7(2):235–41.
Jensen MT, Suadicani P, Hein HO, Gyntelberg F. Elevated resting heart rate, physical fitness and all-cause mortality: a 16-year follow-up in the Copenhagen Male Study. Heart. 2013;99(12):882–7.
Foss OP, Urdal P. Cholesterol for more than 25 years: Could the results be compared throughout all this time? Nor Epidemiol. 2003;13(1):85–8.
Conroy R, Pyörälä K, Fitzgerald A, Sans S, Menotti A, De Backer G, De Bacquer D, Ducimetiere P, Jousilahti P, Keil U. Estimation of ten-year risk of fatal cardiovascular disease in Europe: the SCORE project. Eur Heart J. 2003;24(11):987–1003.
Rabe-Hesketh S, Skrondal A. Multilevel and lognitudinal modeling using Stata. Texas: Stata press; 2008.
Mamelund S-E, Brunborg H, Noack T. Divorce in Norway 1886-1995 by Calendar Year and Marriage Cohort. Oslo-Kongsvinger: Statistics Norway; 1997.
Frisell T, Öberg S, Kuja-Halkola R, Sjölander A. Sibling comparison designs: bias from non-shared confounders and measurement error. Epidemiology. 2012;23(5):713–20.
Gunnell D. Can adult anthropometry be used as a 'biomarker' for prenatal and childhood exposures? Int J Epidemiol. 2002;31(2):390–4.
We acknowledge senior researcher Randi Selmer at the Norwegian Institute of Public Health for sharing knowledge on calculation of the SCORE risk score.
The interpretation and reporting of these data is the sole responsibility of the authors, and no endorsement by the Norwegian Cause of Death Registry is intended nor should be inferred.
This work was supported by the Norwegian Research Council [Grant number 213788].
Availability of data and materials
The dataset analyzed during the current study is not available due to restrictions in the permission given by the Regional Ethics Committee South-East.
ØN conceived the study, acquired the data and handled supervision. IA, LHM, SG-I, HS and ØN participated in the design of the study. IA, LHM, HS, MKRK and ØN participated in the analyses and interpretation of the data. IA drafted the manuscript. LHM, SG-I, HS, MKRK and ØN made critical revisions of the manuscript for key intellectual content, and all authors have read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and conset to participate
The data linkage was approved by the Regional Ethics Committee South-East in May 25th 2012, reference number 2012/872, and the study complies with the Declaration of Helsinki.
Written informed consent was obtained from participants in the Age 40 Program and CONOR. Participants in the Counties studies gave written permission for their results to be sent to their physician. Permission to be absolved from this professional secrecy has been granted and concession to handle this personal health information has been given by the authorities.
Table A Cardiovascular risk factor levels according to level of education in the cohort and within sibships in men. Table B Cardiovascular risk factor levels according to level of education in the cohort and within sibships in women. Table C Cardiovascular risk factor levels according to level of education in the cohort and within sibships in sibships discordant for educational level. (DOCX 46 kb)
About this article
Cite this article
Ariansen, I., Mortensen, L.H., Graff-Iversen, S. et al. The educational gradient in cardiovascular risk factors: impact of shared family factors in 228,346 Norwegian siblings. BMC Public Health 17, 281 (2017). https://doi.org/10.1186/s12889-017-4123-0