Aims: To examine the association between self reported physical workload and low back pain (LBP) in younger twins. To investigate whether genetic factors interact with physical workload in relation to LBP.
Methods: A twin control study was performed within a population based twin register using 1910 complete monozygotic (MZ) and same sexed dizygotic (DZ) twin pairs aged 25–42 and discordant for LBP. LBP in the affected twins was divided into two groups: “LBP for ≤30 days during the past year”, and “LBP for >30 days during the past year”. Physical workload was divided into four categories: “sitting”, “sitting/walking”, “light physical”, and “heavy physical”. Data were analysed in a matched design using conditional logistic regression. MZ and DZ twins were analysed separately and together in order to determine possible genetic influences in relation to physical workload and LBP.
Results: Statistically significant graded relations were found for increasing workload and LBP of longer duration but not for LBP of shorter duration (≤30 days during the past year). In both LBP groups the “sitting” and “sitting/walking” groups were not statistically different. MZ and DZ twins did not differ significantly with respect to LBP in the various workload groups.
Conclusions: There is evidence for a dose-response relation between physical workload and LBP of longer duration. Attention to clinically relevant subgroups based on duration, for example, is necessary in epidemiological studies dealing with LBP. Physical workload might be more important than genetic factors in LBP.
- risk factor
- twin study
- DZ, dizygotic
- LBP, low back pain
- MZ, monozygotic
Statistics from Altmetric.com
Low back pain (LBP) is an enormous public health problem in developed countries. In the United Kingdom alone the combined direct and indirect costs associated with LBP are estimated to be over £12 billion annually.1 This is more than the estimated costs of lower respiratory tract infections, Alzheimer’s disease, stroke, diabetes, multiple sclerosis, and epilepsy combined.1 Considerable research efforts have been devoted to determining potential risk factors in order to identify possible preventive measures. Physical workload, a suspected risk factor, has been investigated extensively in relation to LBP over the past three decades. We are now beginning to see a pattern of heavier physical work as being hazardous for the low back and sedentary work as having either a neutral or protective effect in relation to LBP.2,3 However, the evidence is not completely consistent and many studies are generally considered to be of low methodological quality.2–4
Identification of different varieties or subgroups of LBP was listed as a primary research priority by the Second International Forum for Primary Care Research on Low Back Pain in 1998.5 In 2001 Leboeuf-Yde and Manniche proposed that probably embedded in the term “low back pain” were several, yet unidentified, subentities, some of which might be identified through their selective association with risk factors6—that is, excessive body weight is associated with recurrent or longer lasting LBP, but not with shorter or fewer attacks.7 It is possible that inconsistencies between studies dealing with physical workplace exposures in relation to LBP might be caused by lack of attention to such subentities.
Twin studies provide a special research methodology which, as well as assessing the relative contributions of genetic and environmental factors to a disease, may also be used to examine environmental risk factors in twin pairs discordant for that disease trait.8 This so called “twin control method” can provide information about the way genes and relevant environmental factors interact.9 We performed such a twin control study nested within a population based twin study. Our aims were to examine the pattern of association between self reported physical work load and LBP in workers aged 25–42, and to investigate whether the associations differ depending on the duration of LBP. Furthermore we wanted to investigate whether genetic factors influence this relation.
MATERIALS AND METHODS
This study was based on the “young cohort” of the Danish Twin Register established in 1991. In 1994, a questionnaire survey dealing with a large number of health, lifestyle, and occupational factors was performed. The subsequent database comprises 29 433 twin individuals (11 992 complete pairs) born between 1 January 1953 and 31 December 1982. The establishment and characteristics of this cohort have been described in detail elsewhere.10
The present study was restricted to birth cohorts 1953–68 aged 25–42 at the time of data collection. Of the 13 592 twins (6796 pairs) born 1953–68, 8776 (4388 pairs) were either monozygotic (MZ) or same sexed dizygotic (DZ) pairs. For the analysis of risk factors for LBP, the informative twins were the twin pairs discordant for the disease trait. In this sample 1910 complete pairs (3820 individuals) were discordant with respect to LBP status, meaning that one twin reported LBP within the past year and the co-twin did not (table 1).
LBP was assessed by asking participants how many days they had experienced trouble (pain, tenderness, stiffness, or other trouble) in the lowest part of their back during the preceding year. The question was accompanied by a drawing showing a person from behind with the low back area clearly shaded. The questions relating to the low back were all from the Standardized Nordic Questionnaire.11 Possible answers were then grouped into three categories: “no low back pain”; “LBP short”, that is, LBP for ≤30 days during the past year; and “LBP long”, that is, LBP for >30 days during the past year.
Physical workload was assessed by asking: “How will you describe your work or daily occupation?”. Possible answers were: “sitting”, “sitting and walking”, “light physical work”, and “heavy physical work”.
Occupations for this entire twin cohort were coded according to the DISCO-88 classification system developed by the International Labour Organisation.12,13 Table 2 shows the proportions of twins from discordant pairs in the 10 occupational groups, along with self rated physical work load (table 2).
Twin control method
The twin control method, or the co-twin control method, is a classic matched case-control design that has several advantages over other matched pair designs. First, twin pairs are genetically matched—that is, MZ twins share 100% of their genes and DZ twins share roughly 50%. Second, twins are matched closely for a large number of known and unknown environmental factors: they share prenatal environment and they usually share early life environment such as location, socioeconomic background, nutrition, leisure time activities, and possible toxic exposures. Third, twins are of the exact same age and, in this study, also of the same gender.
In a twin-control study twin pairs discordant for a trait or disease are thus used as matched pairs. Analysis is then performed in a way analogous to a case-control study, with the affected twin being treated as the case and the unaffected co-twin as the control.
Zygosity was determined based on answers to specific questions relating to similarity and mistaken identity. This is a well established and valid method in large twin populations, and comparison with laboratory methods (serological markers) has shown less than 5% misclassification.14
Outcome variables and potential risk indicators (physical workload) were tabulated and checked for logical errors. Incomplete pairs were excluded from further analysis, and because both LBP and physical workload might be gender dependent, opposite sexed DZ twin pairs were also excluded. Paired analyses were then performed using conditional logistic regression, the statistical method of choice in matched designs with categorical data and with more than one level of exposure.9,15 Regression models for the subgroups “LBP short” and “LBP long” were constructed, and odds ratios (OR) with 95% confidence intervals (CI) were estimated. This was done for MZ and DZ twins separately and combined. Subsequently, to accurately estimate the size of the difference between the light and heavy workload groups, ORs for the difference between the two were computed using the post-estimation command “lincom”. To gauge the appropriateness of each model individually, p values for the ORs were computed based on the Wald test, using physical workload as a categorical/continuous variable coded as 1, 2, 3, 4 (table 3); p < 0.05 was considered significant. All analyses were performed using the Stata statistical software package, release 6.0.
Table 1 describes the study sample with respect to gender, age, LBP, and physical workload. The response rate, after one reminder, was 86.2%.
Table 2 shows proportions of twins in 10 occupational groups according to the DISCO-88 classification and the distribution of physical workload in the various groups. Twins in the missing category include persons supplying ambiguous or no specific occupational information.
In the paired analysis (table 3, fig 1) a statistically significant positive gradient, or dose-response like relation, was found between physical workload and LBP of more than 30 days duration during the past year, with ORs ranging from 1.07 to 2.17 for MZ and DZ twins combined. For this group both differences between the light and heavy physical group and the Wald test remained significant. Conversely, for LBP of shorter duration no significant difference between light and heavy work was found; in fact ORs and 95% CIs were practically identical for the two exposures.
For both “LBP short” and “LBP long”, no significant differences were found between the sitting only and the sitting/walking groups, indicating that completely sedentary work and partially sedentary work exhibit an equally neutral effect in relation to LBP.
No significant differences between MZ and DZ twins were found in the two groups, even though in the “LBP long” group for DZ twins the difference between light and heavy physical work reached statistical significance. The ORs, however, were practically identical, and with roughly twice the sample size (528 versus 298), the narrower CIs are hardly surprising and probably do not indicate a real difference.
In this study using 3820 MZ and same sexed DZ twin individuals (1910 complete pairs) we found a significant positive gradient, or “dose-response” relation, between self reported physical workload and LBP of more than 30 days duration, but not for LBP of shorter duration. Interestingly, the Wald test showed a significant trend between increasing workload and LBP of shorter duration; however the difference between the light and heavy groups was far from significant. Thus we conclude that increasing physical workload is a risk factor for LBP of longer duration but probably not for LBP of shorter duration, implying different causal mechanisms for the two subgroups. Therefore treating LBP as a non-specific entity without attention to clinically relevant subgroups, based on duration, for example, is probably inappropriate in epidemiological studies.16
We found no difference with respect to LBP between the group reporting to have only sitting work and the group reporting their work to be both sitting and walking around. This might seem surprising since it has been a common belief that either too little or too much load on musculoskeletal tissues, in this case the lumbar spine, might provoke pain and disability, whereas moderate loads tend to be protective.17 In the light of the present findings and two recent systematic reviews,2,3 we question this assumption, at least in relation to LBP.
We further found the OR for MZ and DZ twins to be very similar for all four categories of workload (table 3). However, ORs for MZ and DZ twins were not completely identical (particularly in the “LBP short” group), indicating a possible interaction between genetic and environmental factors. For “LBP long”, however, the ORs and 95% CIs were similar for MZ and DZ twins, indicating that the influence from physical workload on LBP of longer duration might be greater than the genetic influence. Consequently, it will be necessary in future studies on genes and LBP to include physical workload in the statistical models in order not to overestimate the genetic contribution to this condition. This recommendation is in line with findings by Sambrooke and colleagues,18 Matsui and colleagues,19 and Battié and colleagues20 who, using twins to investigate lumbar disc degeneration (but not pain), also found that heaviness of work influenced the heritability estimates significantly. However, to properly address the issue of the relative importance of genetic versus environmental factors in relation to LBP, path analysis using specialised computer software should be used, and this was not the purpose of the present study.
The results of this study should be interpreted in light of several potential limitations. First, twins have to be representative of the normal singleton population for these results to be valid on a larger scale. Indeed the representativeness of twin studies has been questionned.21 However twins, in spite of an average lower birth weight, have the same prevalence of many adult diseases, including diabetes mellitus22 and thyroid disease23; they have the same fecund ability as ordinary siblings,24 and the same ischaemic heart disease mortality25 and general mortality rate as the general population.26 Furthermore the one year period prevalence (46%) of LBP in this twin sample is comparable to the one year period prevalence (44–54%) found in other studies from the Nordic countries using population based samples.27 Therefore twins are considered to be representative of the general population,8 and we found no reason to challenge this assumption.
Second, the definition of LBP has to be valid. To determine validity an accepted gold standard definition of the disease or condition in question is required. Unfortunately, for LBP such a gold standard does not exist at present. This makes the issues of common acceptance and reproducibility important. The LBP definition used in this study is identical to the definition used in the Standardised Nordic Questionnaire.11 This is a validated questionnaire and widely used in epidemiological LBP research. The Nordic Questionnaire has been criticised for not being specific in distinguishing between different LBP syndromes.4 However, in spite of possible shortcomings, using a commonly accepted and validated questionnaire enables comparison between studies and is normally considered a sign of quality. Also, using a test-retest procedure, Biering-Sorensen and Hilden found the one year recall of LBP using these exact questions to be acceptable.28
Third, questionnaire data on physical workload have been criticised for not offering accurate information on exposure.4 However, questionnaires do classify groups with heterogeneous occupational tasks,29 and questionnaire scores have been shown to be significantly related to current30 and previous31 work status. Furthermore, although efforts have been made to better characterise the variability of exposures between and within occupational groups, no consensus has emerged regarding how to precisely measure exposure to physical loads, or how to apply these measurements in real life.32 Furthermore, detailed measurements on an individual basis in large samples, such as in this database, are both difficult and very costly if possible at all. The information on occupational group and self reported physical workload presented in table 2 adds credibility to our exposure information given detailled knowledge of this classification system and the specific occupations included in each group.13
A statistically significant graded relation exists between physical workload and low back pain of longer duration.
People sitting down all the time do not experience more low back pain than people who both sit and move around.
Attention to clinically relevant subgroups is important in epidemiological studies dealing with low back pain.
Physical workload might more important than genetic factors in low back pain.
Programmes aiming at reducing low back pain in the workplace should target heavier physical work and not sedentary occupations.
The present study was nested in a large general health survey. Thus LBP was not singled out to the participants as the aim of the questionnaire or any subsequent study, and bias arising from attention to any one health problem, in this case LBP, was probably avoided. In addition the high response rate (86.2%) ensures representativeness of this sample. Our results, using a research design not used previously to investigate workload as a risk factor for LBP, confirms the results of two recent systematic reviews in this area.2,3 Furthermore it is in line with the results of Bengtsson and Thorson, who, in the only previous study dealing specifically with LBP and workload in twins, found the proportion of both MZ and DZ twins with “back pain” significantly higher in the “heavy work” group compared to the “light work” group.33
Longitudinal studies are needed to study the long term effects of various workloads in this cohort, and to better study interactions between genetic and environmental factors, such as physical workload or other occupational exposures, and LBP. Some workers, for example, are able to tolerate heavy physical work for many years, whereas others have very limited tolerance and move on to other types of jobs, the so called “healthy worker” effect. Much of the variation in response to seemingly identical exposures could possibly be explained by interactions between the environment, in this case physical workload, and the genetic make-up of the individual. Studies aiming at identifying such interactions are currently being undertaken at the Danish Twin Register.
The authors greatfully acknowledge professor Werner Vach, Dept of Demographics and Biostatistics, University of Southern Denmark, for statistical advice. The principal author was assisted by funding from The Foundation for Chiropractic Research and Postgraduate Education.