Reference Values of Pulmonary Function Tests for Rural Canadians
Chandima P Karunanayake1*, James A Dosman1, Louise Hagel1, Donna C Rennie1,2, Joshua A Lawson1,3, Punam Pahwa1,4 and Saskatchewan Rural Health Study Group
1Canadian Centre for Health and Safety in Agriculture, University of Saskatchewan, Canada
2College of Nursing, University of Saskatchewan, 104 Clinic Place, University of Saskatchewan, Canada
3Department of Medicine, University of Saskatchewan, Royal University Hospital, Canada
4Department of Community Health & Epidemiology, College of Medicine, University of Saskatchewan, Canada
*Corresponding author: Dr. Chandima Karunanayake, Canadian Centre for Health and Safety in Agriculture, 104 Clinic Place, University of Saskatchewan, Saskatoon, Saskatchewan, Canada, Tel: 1- 306- 966-1647, Fax: 1-306-966-8799, E-mail: email@example.com
Int J Respir Pulm Med, IJRPM-2-021, (Volume 2, Issue 2), Original Research Article; ISSN: 2378-3516
Received: June 6, 2015 | Accepted: June 19, 2015 | Published: June 22, 2015
Citation: Karunanayake CP, Dosman JA, Hagel L, Rennie DC, Lawson JA, et al. (2015) Reference Values of Pulmonary Function Tests for Rural Canadians. Int J Respir Pulm Med 2:021. 10.23937/2378-3516/1410021
Copyright: © 2015 Karunanayake CP, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Background: Pulmonary function tests are used routinely to evaluate patients with respiratory diseases and those who are at risk of developing respiratory diseases. Lung function prediction equations are generally developed for urban populations and whether these differ on rural population is not well studied. There is limited information on prediction equations for rural populations.
Objective: The purpose of the present study was to derive prediction equations for commonly used pulmonary function tests in a rural Canadian Caucasian population and to examine whether results of tests were comparable with published equations currently in use.
Methods: The Saskatchewan Rural Health Study was a prospective cohort study conducted in 2010-11. Clinical assessments were conducted on 1675 adults and of those, 1609 had satisfactory pulmonary function tests. Ninety two healthy, asymptomatic lifetime nonsmoking rural Caucasian adults were selected to derive rural Canadian predicted equations. The predictive models were compared with three published models (Crapo, Canadian, and Global Lung Function Initiative 2012) for selected tests. In addition to age and height, other factors such as weight and abdominal girth were considered.
Results: Comparisons with published references showed good agreement with Canadian equations and present study for males but less agreement for females. Weight or abdominal girth did not improve the predictive model significantly.
Conclusion: Although there were slight differences, the prediction equations for rural Caucasians were very much similar to the three published equations, of these two were primarily based on urban subjects (Crapo and Canadian). Therefore, all three published studies can be recommended for the prediction of reference values for rural Canadian Caucasians.
Rural, Pulmonary function tests, Reference values, Spirometry
PFTs:: Pulmonary Function Tests, ATS: American Thoracic Society, SRHS: Saskatchewan Rural Health Study, SPT: Skin Prick Testing, FEV1: Forced Expired Volume in one Second, FVC: Forced Vital Capacity, FEF25-75: Maximum Mid-Expiratory Flow Rate, SEE: Standard Error of the Estimate, GLI: Global Lung Function Initiative
Pulmonary function tests (PFTs) are used routinely to evaluate patients with respiratory diseases and those who are at risk of developing respiratory diseases [1,2]. Normal lung function values and ranges are calculated according to variables such as age, height and weight, and should be based on measurement techniques in accordance with guidelines published by the American Thoracic Society (ATS) [3,4]. The common practice in recent studies [5,6] for determining normal lung function values has been to define the reference population as being that portion of the total population compared exclusively of healthy, lifetime nonsmokers. Some studies discuss other standards for population differences such as race or ethnicity [1,5,7-9], region or country [2,6,10,11], and urban/rural settings [1,12]. In addition, there may be differences in lung function between populations due to varying occupational and environmental exposures and need to be considered when estimating prediction equations [13,14]. Due to these reasons healthy lifelong nonsmokers in a rural environment may be have different lung function results compared to urban counterparts.
Lung function prediction equations are generally developed for urban populations and whether these differ on rural population is not well studied. There is limited information on prediction equations for rural populations. A recent study conducted in Saskatchewan consists of rural Canadian Caucasian population who are mainly exposed to farming and other rural exposures. In our previous studies [15,16], we have used Crapo equations (based on urban Caucasian population) to predict lung function values of rural populations. In this report we would like to investigate whether or not it is appropriate to use prediction equations (such as Crapo) to predict lung function of rural populations by comparing the predicted lung function values based on our study population with those obtained from other established prediction equations [1,5,6].
Study design for adult baseline survey
The Saskatchewan Rural Health Study (SRHS) was a prospective rural cohort and at present, the baseline survey information has been collected. The rural population was defined as consisting of those persons living in towns and municipalities outside the commuting zone of larger urban centres with a population of 10,000 or more . In 2006, there were 44,329 farm enterprises in Saskatchewan encompassing 111,600 household residents . The baseline survey for adults consisted of three stages: (1) recruitment of populations from rural municipalities and small towns; (2) administration of a mail out household questionnaire; (3) clinical assessment that involved anthropometric measures, lung function measurements, and allergy testing on sub-population. Details of the study design are provided elsewhere . In brief, Dillman's method, which involves a series of mail contacts with all prospective participants, was utilized to recruit study participants [19,20]. After excluding ineligible households information on variables of demographics, lifestyle factors, occupational exposures, and socio-economic status based on the Population Health Framework [21,22] was collected by self-administered questionnaires. The study population comprised 8261 individuals (males and females, 18 years of age or older) from 4624 households living in 32 rural municipalities and 15 rural towns in the study area. The study was approved by the Biomedical Research Ethics Board of the University of Saskatchewan, Canada. The SRHS was conducted with the understanding and informed consents of the all participants.
The final question on the baseline questionnaire was 'Would you be willing to be contacted about having breathing and/or allergy tests at a nearby location?' Those who responded positively to this question were sent a letter of invitation to participate in a clinical assessment. Research nurses trained in spirometry and allergy assessment and located in each study quadrant telephoned each household of consenting participants to arrange a time and a place (usually no greater than 60 kilometers from their residence) for this clinical assessment. Mobile clinics were set up in small towns located in the study area. The clinical testing consisted of lung function testing with spirometry, allergy skin prick testing (SPT), blood pressure, and anthropometric measurements (height, weight, abdominal girth). The protocol is described below.
Measurements and selection criteria
Blood pressure was measured as recommended by the Canadian Coalition for High Blood Pressure . Spirometry was completed using the forced expiratory maneuver. The measures of forced expired volume in one second (FEV1), forced vital capacity (FVC), FEV1/FVC ratio, and maximum mid-expiratory flow rate (FEF25-75) were obtained via Sensormedics (Anaheim, CA) dry rolling seal volume displacement spirometers [16,24,25]. Measurements were taken according to standards of the American Thoracic Society [3,4]. There were 3209 persons who consented on the questionnaire to having further testing. Of those, 2863 individuals aged 18-75 years were contacted of which 1675 persons completed clinical testing and 1609 had satisfactory pulmonary function tests. A technically unsatisfactory pulmonary function test is defined as if excessive hesitation of false start, coughing during the first second of maneuver, glottis closure, a leak in the system or around the mouthpiece, and an obstructed mouthpiece [3,4]. For this analysis, 92 rural Caucasian normal healthy subjects were selected (46 males and 46 females). These were subjects those who met the following criteria defined as "normal": a lifetime non-smoker; free of any of following non-respiratory chronic conditions: diabetes, heart disease, heart attack, hardening of the arteries, high blood pressure; free of any of following respiratory/chest illnesses: attack of bronchitis, pneumonia, hay fever, sinus trouble, chronic bronchitis, emphysema, chronic obstructive pulmonary disease, sleep apnea, asthma, cystic fibrosis, tuberculosis and other chest illnesses; and no respiratory symptoms (cough, phlegm or wheeze) or shortness of breath.
Regression models for each PFT, with age, height, weight and abdominal girth as predictive covariates were obtained for each sex. The goodness of fit of the models was determined by their coefficient of determination (R2) and standard error of the estimate (SEE or residual SD) and analysis of residuals. The transformation of the dependent variable into logarithms was tested and did not improve the predictive ability of the model substantially. Assumption of normality of the dependent variables was tested using the Shapiro-Wilk test. If the p-value of the Shapiro-Wilk test was greater than 0.05, the data was considered normal . Predicted values of pulmonary function tests were calculated using Canadian , Crapo  and Global Lung Function Initiative 2012 (GLI2012)  equations and were compared graphically.
The mean, standard deviation (SD) and range of age, height, weight and abdominal girth are given in Table 1. The subjects ranged from 20 to 69 years. Values of R2 and SEE from the predicted equations for the present study (SRHS) and from two previous studies (Canadian study  and Crapo ) are shown in Table 2 (for males) and Table 3 (for females). For the SRHS, both weight and abdominal girth were not significant predictors for all four PFTs for males. For females, both weight and abdominal girth were significant predictors for FVC and FEV1, and abdominal girth was significant predictor for FEV1/FVC ratio. But, there was no significant difference in model fit when abdominal girth was added to the model for females instead of weight. Overall, there was no significant difference in model fit adding covariates weight or abdominal girth to the model for females with predictive covariates age and height. Therefore, the models with predictive covariates age and height were used for comparison purposes. For the SRHS regression models with age and height as predictive covariates, R2 was greater for FVC and FEV1 and R2 was lower for FEF25-75 and FEV1/FVC ratio for males compared to the Canadian  and Crapo  predicted equations. For females, R2 values for the SRHS prediction equations were greater for FVC, FEV1 and FEV1/FVC ratio compared to Canadian study prediction equations. For the SRHS prediction equations, SEE was lower than the value for all PFTs compared to predicted equations of other two studies for both males and females. We were unable to compare GLI2012 equations as different methodology was used for the computation of prediction equations. Hence, there were no available published values of R2 and SEE from GLI2012 equations for Caucasians.
Table 1: Physical Characteristics of normal subjects in SRHS View Table 1
Table 2: Pulmonary function test prediction equations for rural males View Table 2
Table 3: Pulmonary function test prediction equations for rural females View Table 3
Figures 1-3 depict the graphical comparison of predicted pulmonary function values based on the four different prediction equations, namely, SRHS, Canadian , Crapo , and GLI2012 . As indicated in Figure 1 (Panel B), predicted values of FEV1, when corrected for height and age appear to be similar for ages 20 to 70 years old men 175cm in height based on the SRHS and the Canadian study prediction equations . The predicted values of FVC for males based on Crapo  and GLI2012  equations are lower than those based on equations of the Canadian  and the SRHS (Figure 1,Panel A). For ages 20-45 years old, 175cm in height men, values for FEV1/FVC ratio based on SRHS equations are lower than those of the three other studies (Figure 1,Panel C). Values of FVC and FEV1 were similar for the Canadian  and the SRHS prediction equations for heights 155 to 190 cm men 45 years of old age when corrected for height and age (Figure 2,Panel A & B). The FEV1/FVC ratio are higher using Crapo equations  than predicted values obtained using the Canadian, GLI2012 and the SRHS prediction equations for heights 155 to 185 cm and 45 years of old men (Figure 2,Panel C).
Figure 1: Predicted Forced Vital Capacity (FVC) (A), Predicted Forced Expiratory Volume in 1 s (FEV1) (B) and FEV1/FVC ratio (C) by age for 20 to 70 years old men 175cm in height and Predicted Forced Vital Capacity (FVC) (D), Predicted Forced Expiratory Volume in 1 s(FEV1) (E) and FEV1/FVC ratio (F) by age for 27 to 70 years old women 165 cm in height using the present study(SRHS), Canadian, Crapo and GLI2012 equations. View Figure 1
Figure 2: Predicted Forced Vital Capacity (FVC) (A), Predicted Forced Expiratory Volume in 1 s(FEV1) (B) and FEV1/FVC ratio (C) by height for 155cm to 190cm men 45 years old age and Predicted Forced Vital Capacity (FVC) (D), Predicted Forced Expiratory Volume in 1 s(FEV1) (E) and FEV1/FVC ratio (F) by height for 150cm to 175cm women 45 years old age using the present study(SRHS), Canadian, Crapo and GLI2012 equations. View Figure 2
Figure 3: Predicted Forced Expiratory Flow during the middle half of the Forced Vital Capacity (FEF25-75) for a male 175cm in height for 20 to 70 years (A), Predicted forced expiratory flow during the middle half of the forced vital capacity (FEF25-75) for a female 165cm in height for 27 to 70 years (B) using the present study (SRHS), Crapo and GLI2012 equations. View Figure 3
For FVC, we observed slightly higher values using the SRHS prediction equations compared to the predictions from other three studies for ages 27 to 55 and 165cm tall female (Figure 1,Panel D) and for heights 150 to 175cm and 45 years old female (Figure 2,Panel D). For FEV1, we observed slightly higher values using the SRHS prediction equations compared to the values predicted from other three studies for ages 27 to 45 years and 165cm tall female and for heights 150 to 165cm and 45 years old female (Figure 1,Panel E and Figure 2,Panel E respectively). Also we observed slightly lower values for FEV1 using prediction equations of the SRHS compared to Canadian study  for ages 50 to 70 years and 165cm tall female (Figure 1,Panel E). The SRHS predicted values for FEV1/FVC ratio were lower than those predicted using prediction equations of the Canadian , GLI2012  and Crapo  for females with ages 27 to 70 (Figure 1,Panel F) and heights 157 to 175cm (Figure 2,Panel F).
For the Canadian study , the prediction equations for FEF25-75 was not available, hence the comparison of predicted FEF25-75 values with the SRHS predicted equations was not possible.
However, the predicted FEF25-75 values using SRHS equations are higher than those using prediction equations for the GLI2012  and Crapo study  for males of 175cm aged from 20 to 70 years (Figure 3,Panel A). The predicted FEF25-75 values were not much different using the SRHS and GLI2012  prediction equations for females of 165cm aged from 27 to 70 years. The predicted FEF25-75 values using SRHS equations are lower than those using prediction equations of the GLI2014  and Crapo study  for females of 165cm aged from 27 to 70 years (Figure 3,Panel B). However there were differences of FEF25-75 predicted values using the predicted equations of three studies (Crapo, GLI2012 and SRHS) for 45 year old male or female with heights ranges from 155 to 190cm and 150 to 175cm, respectively.
We have modeled equations for several PFTs based on data available for normal healthy, rural Caucasian lifetime nonsmokers from Saskatchewan, Canada. We found that, although there were slight differences, the predictions of FVC and FEV1 were very much similar to those based on the Canadian study  equations for males, which mostly included urban subjects. Therefore there is no difference in PFTs for males with respect to urban or rural setting. For females, the predictions of FVC and FEV1 based on the SRHS equations did not show very good agreement with predicted values obtained using Canadian , Crapo  and GLI2012  equations. For both males and females we did not see a very good agreement of FEV1/FVC ratio predictions using the prediction equations of the present study and other three studies [1,5,6]. There was a very good agreement for females of FEF25-75 predicted values of the SRHS equations with GLI2012  equations and we did not observe this in males.
We thought that rural Caucasians might have values for pulmonary function that would be affected by the high rates of obesity in rural people . However that did not seem to be true. We did not observe significant improvement of the fit of the predictive models when covariates weight and abdominal girth were included in the models. Supporting this idea, a recent study reported that there was no effect from obesity on spirometry tests .
Caucasians were selected in all four studies to derive prediction equations and were believed to have a similar ethnic origin. Crapo's model was generated in the United States in Caucasian urban populations and the Canadian study, mostly included urban subjects from Canada. The GLI study included Caucasians from Brazil, Australia, Canada, USA, Israel, Mexican Americans, Chile, Mexico, Uruguay, Venezuela, Algeria, Tunisia and many European countries. All four studies complied with ATS recommendations and techniques. Differences in prediction equations based on these studies could be due to several reasons, such as, sample size variations among the four studies and rural versus urban populations. It can be assumed that the rural populations tend to be more homogeneous than urban populations in terms of environmental exposures. However, these findings on differences in prediction equations warrant further investigation.
A strength of this study was that the selected subjects lived in a rural setting and no urban residents were included in the SRHS. This helped us to establish prediction equations for rural Canadian Caucasians. One of the limitations of the present study is the relatively small sample size with 46 males and 46 females. According to Schneider et al.  for linear regression model in general, the number of observations must be at least 20 times greater than the number of variables under study. As our sample sizes are greater than 40 for males and females, the multiple linear regression models with age and height are reasonable for prediction of PFTs. However, due to small samples results from models with three variables ([age, height, weight] and [age, height, abdominal girth]) need to be interpreted with caution. Another potential limitation was that the selection criteria depended on participants consent to taking part in clinical assessment and then screening of self-reported responses of the subjects, which can lead to biased results . In the present study subjects were screened as normal healthy subjects if they were asymptomatic to respiratory diseases and symptoms and lifetime nonsmokers using the responses from self-reported questionnaire. In the other three studies participants were highly screened with a clinical examination and laboratory tests [1,5,6]. Another limitation was the length of rural residency was unknown. Currently, the participants lived in rural area, but do not know how long they were there.
In summary, we compared our models for selected PFTs with Canadian, Crapo and GLI2012 prediction equations, and found them to be fairly similar overall, although there were some differences for both males and females. Therefore, the prediction equations from all three published studies [1,5,6] can be recommended for the prediction of reference values for rural Canadian Caucasians.
The Saskatchewan Rural Health Study Team consists of: James Dosman, MD (Designated Principal Investigator, University of Saskatchewan, Saskatoon, SK Canada); Dr. Punam Pahwa, PhD (Co-principal Investigator, University of Saskatchewan, Saskatoon SK Canada); Dr. John Gordon, PhD (Co-principal Investigator, University of Saskatchewan, Saskatoon SK Canada); Yue Chen, PhD (University of Ottawa, Ottawa Canada); Roland Dyck, MD (University of Saskatchewan, Saskatoon SK Canada); Louise Hagel (Project Manager, University of Saskatchewan Saskatoon SK Canada); Bonnie Janzen, PhD (University of Saskatchewan, Saskatoon SK Canada); Chandima Karunanayake, PhD (University of Saskatchewan, Saskatoon SK Canada); Shelley Kirychuk, PhD (University of Saskatchewan, Saskatoon SK Canada); Niels Koehncke, MD (University of Saskatchewan, Saskatoon SK Canada); Joshua Lawson, PhD, (University of Saskatchewan, Saskatoon SK Canada); William Pickett, PhD (Queen's University, Kingston ON Canada); Roger Pitbaldo, PhD (Professor Emeritus, Laurentian University, Sudbury ON Canada); Donna Rennie, RN, PhD, (University of Saskatchewan, Saskatoon SK Canada); Ambikaipakan Senthilselvan, PhD (University of Alberta, Edmonton, AB, Canada). This study was funded by a grant from the Canadian Institutes of Health Research "Saskatchewan Rural Health Study", Canadian Institutes of Health Research MOP-187209-POP-CCAA-11829. We are grateful for the contributions of the rural municipality administrators and the community leaders of the towns included in the study that facilitated access to the study populations and to all of participants who donated their time to complete and return the survey.
Crapo RO, Morris AH, Gardner RM (1981) Reference spirometric values using techniques and equipment that meet ATS recommendations. Am Rev Respir Dis 123: 659-664.
Gore CJ, Crockett AJ, Pederson DG, Booth ML, Bauman A, et al. (1995) Spirometric standards for healthy adult lifetime nonsmokers in Australia. Eur Respir J 8: 773-782.
Miller MR Crapo R, Hankinson J, Brusasco V, Burgos F, et al. (2005) General considerations for lung function testing. Eur Respir J 26: 153-161.
Miller MR, Hankinson J, Brusasco V, Burgos F, Casaburi R, et al. (2005) Standardisation of spirometry. Eur Respir J 26: 319-338.
Quanjer PH, Stanojevic S, Cole TJ, Baur X, Hall GL, et al. (2012) Multi-ethnic reference values for spirometry for the 3-95-yr age range: the global lung function 2012 equations. Eur Respir J 40: 1324-1343.
Gutierrez C, Ghezzo RH, Abboud RT, Cosio MG, Dill JR, et al. (2004) Reference values of pulmonary function tests for Canadian Caucasians. Can Respir J 11: 414-424.
Withers RT, Bourdon PC, Crockett A (1989) Spirometric standards for healthy male lifetime nonsmokers. Hum Biol 61: 327-342.
Withers RT, Lemmey AB (1989) Lung volume and spirometric standards for healthy female lifetime nonsmokers. Hum Biol 61: 343-368.
Gibson J, Gallagher H, Johansen A, Webster I (1979) Lung function in an Australian population. I. Spirometric standards for nonsmoking adults. Med J Aust 1: 292-295.
Roberts CM, MacRae KD, Winning AJ, Adams L, Seed WA (1991) Reference values and prediction equations for normal lung function in a non-smoking white urban population. Thorax 46: 643-650.
Duong M, Islam S, Rangarajan S, Teo K, O'Byrne PM, et al. (2013) Global differences in lung function by region (PURE): an international, community-based prospective study. Lancet Respir Med 1: 599-609.
Morris JF, Koski A, Johnson LC (1971) Spirometric standards for healthy nonsmoking adults. Am Rev Respir Dis 103: 57-67.
Schnabel E, Chen CM, Koch B, Karrasch S, Jörres RA, et al. (2010) Regional differences in prediction models of lung function in Germany. Respir Res 11: 40.
Pellegrino R, Viegi G, Brusasco V, Crapo RO, Burgos F, et al. (2005) Interpretative strategies for lung function tests. Eur Respir J 26: 948-968.
Pahwa P, Karunanayake CP, Rennie DC, Chen Y, Schwartz DA, et al. (2009) Association of the TLR4 Asp299Gly polymorphism with lung function in relation to body mass index. BMC Pulm Med 9: 46.
Chénard L, Senthilselvan A, Grover VK, Kirychuk SP, Lawson JA, et al. (2007) Lung function and farm size predict healthy worker effect in swine farmers. Chest 131: 245-254.
du Plessis V, Beshiri R, Bollman RD, Clemenson H (2004) Definitions of "Rural".Agriculture and Rural Working Paper Series, Working Paper No. 61.Catalogue no. 21-601-MIE- No. 061. Ottawa, Canada: Agricultural Division, Statistics Canada.
Statistics Canada (2006) Census of Agriculture. Ottawa, ON: Minister of Industry.
Pahwa P, Karunanayake CP, Hagel L, Janzen B, Pickett W, et al. (2012) The Saskatchewan Rural Health Study: An Application of a Population Health Framework to Understand Respiratory Health Outcomes. BMC Res Notes 5: 400.
Dillman DA (2000) Mail and internet surveys: the tailored design method. (2nd edn). New York, NY: Wiley & Sons Inc, USA.
Health Canada (1994) Strategies for Population Health: Investing in the Health of Canadians, Cat. No. H39-316/1994E. Ottawa, On: Health Canada, Minister of Supply and Services Canada.
Pickett W, Day L, Hagel L, Brison RJ, Marlenga B, et al. (2008) The Saskatchewan Farm Injury Cohort: rationale and methodology. Public Health Rep 123: 567-575.
Abbott D, Campbell N, Carruthers-Czyzewski P, Chockalingam A, David M, et al. (1994) Guidelines for measurement of blood pressure, follow-up, and lifestyle counselling. Canadian Coalition for High Blood Pressure Prevention and Control. Can J Public Health 85 Suppl 2: S29-S43.
Senthilselvan A, Chénard L, Ulmer K, Gibson-Burlinguette N, Leuschen C, et al. (2007) Excess respiratory symptoms in full-time male and female workers in large-scale swine operations. Chest 131: 1197-1204.
Senthilselvan A, Dosman JA, Pahwa P (1991) Respiratory symptoms and alterations in pulmonary tests in poultry producers. Am Rev Resp Dis 143: A100.
Shapiro SS, Wilk MB (1965) An Analysis of Variance Test for Normality (Complete Samples). Biometrika 52: 591-611.
Shields M, Tjepkema M (2006) Regional differences in obesity. Health Rep 17: 61-67.
Al Ghobain M (2012) The effect of obesity on spirometry tests among healthy non-smoking adults. BMC Pulm Med 12: 10.
Schneider A, Hommel G, Blettner M (2010) Linear regression analysis: part 14 of a series on evaluation of scientific publications. Dtsch Arztebl Int 107: 776-782.
Pahwa P, Karunanayake C, Hagel L, Janzen B, Rennie D, et al. (2012) Self-selection bias in an epidemiological study of respiratory health of a rural population. J Agromedicine 17: 316-325.