Abstract
Monitoring of chronic lung disease requires reference values of lung function indices, including putative markers of small airway function, spanning a wide age range.
We measured spirometry, transfer factor of the lung for carbon monoxide (TLCO), static lung volume, resistance and ventilation distribution in a healthy population, studying at least 20 subjects per sex and per decade between the ages of 20 and 80 years.
With respect to the Global Lung Function Initiative reference data, our subjects had average z-scores for forced expiratory volume in 1 s (FEV1), forced vital capacity (FVC) and FEV1/FVC of −0.12, 0.04 and −0.32, respectively. Reference equations were obtained which could account for a potential dependence of index variability on age and height. This was done for (but not limited to) indices that are pertinent to asthma and chronic obstructive pulmonary disease studies: forced expired volume in 6 s, forced expiratory flow, TLCO, specific airway conductance, residual volume (RV)/total lung capacity (TLC), and ventilation heterogeneity in acinar and conductive lung zones.
Deterioration in acinar ventilation heterogeneity and lung clearance index with age were more marked beyond 60 years, and conductive ventilation heterogeneity showed the greatest increase in variability with age. The most clinically relevant deviation from published reference values concerned RV/TLC values, which were considerably smaller than American Thoracic Society/European Respiratory Society-endorsed reference values.
Abstract
We present reference equations for median and variability of most lung function indices, in the age range 20–80 years http://ow.ly/SsZFW
Introduction
Studies in patients with asthma or smoking-induced lung disease require appropriate reference values to determine lung function abnormality, especially in the early stages of the disease process. While asthma can arise at any age, chronic obstructive pulmonary disease (COPD) usually develops later in life, yet both chronic conditions lead to follow-up periods of lung function potentially extending through old age. Both conditions are usually marked by a spirometric deficit which can now be set against reference equations that are valid worldwide [1, 2]. However, the characterisation of asthma and smoking-induced lung disease can also benefit from lung function and ventilation distribution indices that are thought to be sensitive to two distinct but sometimes concomitant phenomena: peripheral airway abnormality and abnormal heterogeneity of airway function at all lung depths. Hence, clinical research studies regularly include indices dedicated to represent peripheral and heterogeneous airway function [3]. Abnormality of these lung function indices can sometimes be assessed by comparison with reference values, either obtained for a particular index obtained in one study population [4, 5] or based on reference equations constructed from the summary of previous equations [6]. To date, the most complete set of reference values acquired on a well-defined population (n=179; 18–86 years) is by Roberts et al. [7], since it includes spirometry, plethysmography, transfer factor of the lung for carbon monoxide (TLCO) and single-breath washout indices (phase III slope) on nonsmoking adults.
The main objective of this work was to provide easy-to-use reference equations on all lung function and ventilation indices that have proven to be relevant to asthma, COPD or cystic fibrosis (CF) to some extent in the past. On the lung function side, TLCO and the ratio of residual volume (RV) over total lung capacity (TLC) are particularly relevant to COPD and asthma. On the ventilation distribution side, indices of acinar and conductive (Sacin, Scond) ventilation heterogeneity derived from the normalised phase III slope analysis of the multiple-breath washout (MBW) have been used to characterise bronchial hyperresponsiveness in asthma [8] or asthma control [9], or to pick up different stages of smoking-induced lung disease [10]. In more severe lung disease such as in patients with COPD or patients with CF, the lung clearance index (LCI) [11] and the curvilinearity (Curv) derived from the MBW concentration curve [12] may suffice to quantify degree of heterogeneity. Airway resistance is also seen to be increased in asthma and COPD, with values depending on the method (forced oscillation or plethysmography) [13] or the equipment used [14]. Even if most indices are bound to be partly apparatus- and method-dependent, the dependence of any physiological index on age and/or height will likely be similar across laboratories. Here, we aimed to obtain all the above indices on the same population, in the same lung function laboratory of the university hospital and coached by the same well-trained lung function technician.
Materials and methods
The study protocol was approved by the local UZ Brussel research ethics committee (B14320097554 and B143201112421 pertaining to data covering ≤65 and >65 years age groups, respectively); selected spirometric and MBW indices from our centre have been published previously for the 25–65 years group [15]. Lung function and MBW tests were collected on healthy nonsmoking subjects with <3 pack-years (in case of ex-smokers) in the age range 20–80 years. All participating subjects were Caucasian, were required to have a body mass index <33 kg·m−2 and were defined as healthy through clinical screening, with no medical history of respiratory disease. At least 40 subjects (20 males/20 females) were recruited per decade in the age range 20–80 years, based on the recommended number of subjects described by Quanjer et al. [16]. We adapted the numbers needed to the age range under study here (20–80 years) where age dependency of most indices was expected to be monotonous, in contrast to what occurs in the age range <20 years [1]. However, because we anticipated different variability around the median in different age ranges for at least some indices, we applied generalised additive models for location, scale and shape method proposed by the Global Lung Function Initiative (GLI) [1, 2] to obtain LambdaMuSigma (LMS) model fits. From the resulting parameters L, M and S, predicted values and lower and upper limit of normal (LLN, ULN) can be readily computed; see example in online supplementary material.
All lung function measurements were performed in the same centre (centrally located in Belgium), by the same lung function technician, according to the American Thoracic Society (ATS)/European Respiratory Society (ERS) guidelines for lung function measurements [17] using commercial equipment (Vmax Encore VE models 20c, 22 and 22d; Cardinal Health, Dublin, OH, USA). After standard lung function testing including spirometry, TLCO, static lung volumes and airway resistance by body plethysmography, MBW tests were performed as previously described [15]. Volumes from all above-mentioned tests were reported at body temperature, ambient pressure, saturated with water vapour.
Spirometry
In each subject, three acceptable manoeuvres were performed that satisfied the reproducibility criteria of the ATS/ERS guidelines [18]. First, we selected the largest forced expired vital capacity (FVC), the largest forced expired volume in 1 s (FEV1) and the largest forced peak expiratory flow (PEF) obtained over the three blows. Then, we selected the blow with the largest sum of FEV1 and FVC, from which we extracted forced expired volume in 6 s (FEV6), forced expiratory flow at 75% of FVC (FEF75%) and forced expiratory flow at 25–75% of FVC (FEF25–75%). Finally, derived indices FEV1/FVC and FEV1/FEV6 were considered for developing reference equations.
Carbon monoxide transfer measurement
Two acceptable tests were performed in each subject [19]. Using the single-breath technique and 21% O2/0.3% CO/0.3% C2H2/0.3% CH4 with balance N2 as the inspiratory gas mixture, we obtained TLCO, alveolar volume at total lung capacity (VA) and transfer coefficient of the lung for carbon monoxide (KCO) (i.e. TLCO/VA); all indices were reported as the mean of values obtained from both tests.
Static lung volumes and resistance using body plethysmography
Two acceptable tests were performed in each subject [20]. We retained the measurement of RV determined by whole-body plethysmography (RVpleth), functional residual capacity determined by whole-body plethysmography (FRCpleth), total lung capacity determined by whole-body plethysmography (TLCpleth), airway resistance (Raw) and specific airway conductance (sGaw). All indices were reported as the mean of values obtained from both tests. To ensure extraction of proper FRCpleth values with the body plethysmograph model that we used (Vmax Encore 22d), values labelled in the report as “FRC” and “Vtg” (thoracic gas volume) were examined to determine whether FRC corresponded to the actual FRC level observed on the traces. This had to be done because biological calibrators and tests with two independent plethysmograph boxes had indicated an artefact in FRCpleth (despite proper RVpleth and TLCpleth). With our particular bodybox and its software settings this was due to an unpredictable volume drift upon closing the shutter valve, such that the relative levels corresponding to Vtg and FRC needed readjusting.
MBW test
A bag-in-box and valve system was used, with a re-inspired dead space (i.e. re-inspired volume in subsequent MBW inhalations) of 50 mL and a N2 analyzer (PK-Morgan, Chatham, UK). The MBW test involved 1-L tidal breathing from FRC. Three acceptable tests were performed in each subject. Each MBW test was analysed with a semi-automated break-point analysis [21] from which phase III slope and mean expired concentration were obtained for every breath. FRCmbw was obtained by mass balance including all breaths up to the one corresponding to mean expired concentration just under 1/0.40th of the initial N2 concentration. From three acceptable tests, data were cumulated to obtain plots of mean expired concentration and phase III slope versus lung turnover (TO; cumulative expired volume/FRCmbw for any given test). From the semi-log plot of mean expired concentration (as a percent of initial concentration) versus TO we derived: 1) LCI, as the TO value corresponding to the intercept between the semi-log plot of mean expired concentration and the 1/0.40th initial concentration level; 2) Curv, as unity minus the ratio of two regression slopes, where the regression slope in the denominator is that determined in the range TO=0 and TO=LCI/2 and in the numerator that between TO=LCI/2 and TO=LCI [12]. From the plot of normalised phase III slope (normalised by mean expired concentration) versus TO, we derived indices of acinar and conductive (Sacin, Scond) ventilation heterogeneity, as well as two modified indices (Sacin*, Scond*) for use in patient groups such as CF where ventilation heterogeneity is so marked that Scond becomes invalid [22]. Finally, we determined Fowler dead space volume of the first breath (VDF) in order to assess its potential impact on LCI when exploring the relationship between LCI, which is an overall index of ventilation heterogeneity, and its potential determinants (Sacin, Scond).
Statistical analysis
Simple comparison between selected indices, linear and stepwise multiple regressions were performed using MedCalc (version 14.8.1; MedCalc, Mariakerke, Belgium). For the regression equations, we followed the guidelines and procedures described in Quanjer et al. [1] and Lum et al. [23], using the GAMLSS package (version 4.3-4) in R (version 2.15.2; R Foundation, Vienna, Austria). After Box–Cox–Cole–Green transformation of the data, an LMS model fit was obtained that allows the median and spread of the distribution to vary with height and age (or age squared). We also included sex and tested for interaction between age and sex in the initial estimate of the median. For Curv, which had some zero values, the Box–Cox–Cole–Green transformation required the addition of a very small positive value (+0.001) to the actual Curv measurement value. With the Box–Cox–Cole–Green transformation, the obtained coefficients for M are to be used in a linear expression to obtain the predicted value, whereas the obtained coefficients for S are to be used in an exponential expression to obtain the spread; L corresponds to nu estimated by a constant coefficient. The LLN and ULN (5th and 95th percentile) can then be computed as M×(1–1.645×L×S)(1/L) and M×(1+1.645×L×S)(1/L) [23]. Guided by the Schwarz Bayesian Criterion, the most suitable model with likely contributors was selected for each variable. In all the above analyses, statistical significance was accepted at p<0.05.
Results
The total cohort comprised of 252 subjects, who were homogeneously distributed over the six decades between 20 and 80 years, and over both sexes (figure 1a) and with maximum BMI just under 32 kg·m−2. To characterise our study population with respect to GLI2012, we determined z-scores for FEV1, FVC and FEV1/FVC based on the GLI2012 reference values for Caucasians [2]: mean (±sd) z-scores for FEV1, FVC and FEV1/FVC were −0.12 (±0.90), +0.04 (±0.85) and −0.32 (±0.71), respectively. Table 1 summarises the mean±sd of all variables studied here.
For most variables, reference equations were constructed, resulting in coefficients for sex, height, and age and/or age squared to obtain L, M and S (table 2); an example computation of predicted, ULN and LLN is provided in the online supplementary material. Figure 2 shows the raw data and their LMS fit (median, LLN, ULN) for two lung function indices that are particularly relevant for studies in asthma and COPD: TLCO and RV/TLC. Figure 3 shows the raw data and their LMS fit for three relevant ventilation distribution indices: LCI, Sacin and Scond. In figures 2 and 3, left and right panels are, respectively, female and male data drawn to the same vertical axis; in order to be able to also plot the LMS fit for those indices that depend on height, height was taken into account in each sex group by use of the height–age relationships of figure 1(b). The raw data plots with superimposed LMS fit for all indices of table 1 and their colour-coded comparison to other reference values, wherever available [4–7, 24–30], are in the online supplementary material.
No reference equations were drawn up for FRCpleth and VA due to methodological considerations (see Discussion section), but these indices were compared with FRCmbw and TLCpleth, respectively (figure 4). A Bland–Altman representation of FRCmbw and FRCpleth showed good overall agreement, while VA consistently fell below TLCpleth (by 0.46 L). Multilinear regression of LCI versus its potential contributors (Sacin, Scond, VDF/FRC and FRC) retained the first three indices, which together accounted for 58% of the LCI variability: partial correlation coefficients were +0.67 for Sacin (p<0.0001), +0.13 for Scond (p=0.04) and +0.20 for VDF/FRC (p=0.002).
Discussion
We have developed reference equations, LLN and ULN for indices derived from spirometry, TLCO, static lung volume, resistance and ventilation distribution tests. The lung function indices studied are directly relevant to common respiratory illnesses, including asthma, COPD and CF.
Spirometry
As GLI2012 already offers worldwide valid reference equations [1], we used these to verify that our subjects obtained z-scores close to zero. The coefficients provided in table 2 for the most widely used indices already covered by the GLI, such as FEV1 and FVC, served to verify that the statistical methodology was sound and that the expected explanatory variables were obtained. We then produced reference equations for FEV6 and FEV1/FEV6 which have been suggested for use to screen smokers for restriction and obstruction in the primary care setting [31]. In the online supplementary figures S4 and S5, the LLN values for FEV6 and FEV1/FEV6 are seen to be close to those reported by Hankinson et al. [24]; they almost coincide in the age range 45–65 years, i.e. the age category where screening for smoking-induced lung disease is probably most relevant. The diagnostic value of FEF75% or FEF25–75% in adults has been questioned in recent years, e.g. in terms of their added value with respect to other spirometric indices when classifying patients [32]. However, for completeness we also included FEF75% and FEF25–75% in the online supplementary material. Finally, the PEF (online supplementary figure S6) showed excellent agreement with both Kuster et al. [26] and Hankinson et al. [24], identifying similar independent contributors (sex, age, height) (table 2).
Transfer factor and body plethysmography
TLCO and TLCO/VA show the expected dependence on sex, height and age consistent with fragmented data from the literature [4, 7, 25]. In figure 2 it can be appreciated that TLCO shows a monotonous decrease over the entire age range with variability remaining almost constant throughout. Normal RVpleth/TLCpleth values (figure 2b and online supplementary figure S13) were considerably smaller than those in the most recent ATS/ERS report [27]. The authors rightly pointed out that in the absence of more recent data, reference equations for lung volumes effectively refer to the ERS1983 reference equation for adults constructed from earlier equations resulting from plethysmographic and ventilation lung volume studies, and which potentially included smokers [6]. The RVpleth/TLCpleth obtained here were closer to RV/TLC data obtained in 1991 by Roberts et al. [7] using helium dilution by closed-circuit rebreathing and obtained in nonsmokers.
Finally, VA was consistently lower than TLCpleth and indeed paralleled TLCpleth over the entire range (figure 4), as could be expected [17]. It has been postulated that the ratio of TLC derived from a single- or multiple-dilution test could be a marker of incomplete ventilation distribution [33], potentially making VA/TLC a useful index in its own right. However, using TLCpleth instead of TLC by dilution (multibreath or rebreathing) adds complexity to the interpretation of VA/TLCpleth as a surrogate index, especially in obstructive patients. As a result of manufacturer-implemented corrections in the plethysmographic bodybox and the way the so-called thoracic gas volume (“Vtg”; determined after the valve closure) is related to actual FRC, FRCmbw was chosen over FRCpleth as the preferred method to draw up reference equations for FRC. Nevertheless, when FRCpleth was compared with FRCmbw (figure 4), overall agreement was good, indicating general coherence of the data set. FRCmbw slightly decreased with age, via its dependence on height (which decreases with age; figure 1b), but age was not an independent predictor. This agrees with the very small age dependency obtained from dilution FRC in Roberts et al. [7], which is the same in both sexes. This contrasts with the 10-fold greater age dependence in males versus females in the most used reference equations for FRC [27], probably due to a problem with summarised prior reference equations and the mixture of FRC values derived from ventilation and plethysmographic measurement [34].
To the best of our knowledge, there are no reference equations for airway resistance (Raw) obtained from body plethysmography and we can only compare the slight increase with age that we find in the age range >45 years with that obtained by 5 Hz forced oscillations [5]. The most noteworthy observation of our data set spanning 20–80 years is that there appears to be a minimum in Raw around the age of 45 years (online supplementary figure S14). Specific conductance (sGaw) is seen to remain relatively constant up to 50–55 years, after which it slightly deteriorates (online supplementary figure S15). No published reference values exist for Raw or sGaw covering the 20–80 years range.
MBW test
First, using the raw data, we examined whether there was a link between LCI and Sacin and Scond, as was found to be the case in asthma patients [35]. In the present study, the “disease” driving Sacin and Scond values would be ageing, and indeed both appeared to correlate independently with LCI. An independent contribution from VDF as a fraction of FRC also contributed to an increased LCI. In this case, ageing per se was probably not the driving factor since FRC and VDF/FRC (online supplementary figures S16 and S17) do not show a significant dependence on age (table 2).
The reference equations for LCI, Scond and Sacin (figure 3 and table 2) essentially extend those previously obtained across two laboratories in a narrower age range (25–65 years) [15]. The greatest difference between that study and the present one is that in the age range up to 80 years, LCI and Sacin no longer vary linearly with age, but both increase more rapidly as age increases, resulting in an age-squared dependence. This implies that in the case of smoking-induced lung disease, which often appears later in life, increased acinar ventilation heterogeneity previously reported by our group may have been partly due to ageing [36]. Nevertheless, the average Sacin values observed in the COPD patient groups without emphysema (0.295 L−1; mean age 60 years) and with emphysema (0.444 L−1; mean age 68 years) still largely exceeded the ULN obtained here in the corresponding 60–68 years range [36].
All three MBW indices of figure 3 are also characterised by a greater variability in the advanced age range, which is most marked for Scond (and its modified version Scond*). The increase of average Scond, as well as intersubject variability in Scond, with age is probably due to a greater intersubject heterogeneity of elastic properties of lung units larger than acini. The most likely explanation is loss of elastic recoil and changes in the elastic fibre network that are suggested to occur with ageing [37], bringing about a slight enlargement of respiratory bronchioles and alveolar ducts, without clear destruction of the alveolar walls. Probably, such structural changes are heterogeneously distributed over the lungs and do not occur to the same extent in all subjects. A recent computed tomography study also demonstrated age-dependent air trapping and airway remodelling in normal subjects aged 60–90 years, and the authors identified this to be a heterogeneous process [38].
Limitations of the study
Scaling down the number of subjects required, based on the n=150 per sex group proposed by Quanjer et al. [16] to obtain data in the range 5–90 years, we would have needed to include 15 subjects per decade and sex or 180 subjects in total. For lung volumes, the recommended number by Stocks et al. [27] to obtain adult reference equations is 100 for men and women. One could argue that the recommended numbers may not hold for all other indices, yet the reasoning behind the recommended sample size in Quanjer et al. [16] was based on z-scores, which should in principle apply to all indices measured on normal subjects. Nevertheless, we accept that there is something arbitrary about the choice of a criterion for the number of subjects. Here, we refer to our data set (with a homogeneous distribution over the entire age range and a low number of outliers) and the overall agreement with predictions from other studies wherever available, thus attesting to the adequacy of the number of subjects studied here. Another limitation concerns the equipment used, which mainly pertains to indices derived from MBW. While MBW indices Scond and Sacin are less sensitive to equipment dead space than LCI [15], all three indices are sensitive to the gas used. This implies that the absolute values provided by the reference equations may not be blindly applied to other devices or laboratories elsewhere. However, the observed sex or age dependence of these indices and their variability should be similar, and can be borne in mind when selecting a group of healthy subjects locally in the validation phase of a newly installed MBW device in any given laboratory. A final limitation of the study is that the reference equations cannot be readily used for subjects of ethnicity other than Caucasian or for Caucasian subjects outside the normal length and weight ranges.
In conclusion, we obtained reference equations for most lung function indices currently implemented in a lung function laboratory with a keen interest in pulmonary research. For some variables (Raw, sGaw, VDF/FRC), no such reference data existed; for others (RV/TLC, FRC), existing reference equations were judged inadequate, and for ventilation distribution data, these showed an accelerated deterioration (LCI, Sacin) and increased variability (Scond) beyond 60 years.
Footnotes
This article has supplementary material available from erj.ersjournals.com
Support statement: This study was supported by the Fund for Scientific Research – Flanders. Funding information for this article has been deposited with FundRef.
Conflict of interest: None declared.
- Received May 4, 2015.
- Accepted September 3, 2015.
- Copyright ©ERS 2016