Abstract
The recently published Global Lung Function Initiative (GLI) carbon monoxide transfer factor (TLCO) reference equations provide an opportunity to adopt a current, all-age, widely applicable reference set. The aim of this study was to document the effect of changing to GLI from commonly utilised reference equations on the interpretation of TLCO results.
33 863 TLCO results (48% female, 88% Caucasian, n=930 aged <18 years) from clinical pulmonary function laboratories within three Australian teaching hospitals were analysed. The lower limit of normal (LLN) and proportion of patients with a TLCO below this value were calculated using GLI and other commonly used reference equations.
The average TLCO LLN for GLI was similar or lower than the other equations, with the largest difference seen for Crapo equations (median: −1.25, IQR: −1.64, −0.86 mmol·min−1·kPa−1). These differences resulted in altered rates of reduced TLCO for GLI particularly for adults (+1.9% versus Miller to −27.6% versus Crapo), more so than for children (−0.8% versus Kim to −14.2% versus Cotes). For adults, the highest raw agreement for GLI was with Miller equations (94.7%), while for children it was with Kim equations (98.1%). Results were reclassified from abnormal to normal more frequently for younger adults, and for adult females, particularly when moving from Roca to GLI equations (30% of females versus 16% of males).
The adoption of GLI TLCO reference equations in adults will result in altered interpretation depending on the equations previously used and to a greater extent in adult females. The effect on interpretation in children is less significant.
Abstract
Adoption of GLI TLCO reference equations in adults will result in altered interpretation depending on the equations previously used and to a greater extent in adult females. The effect on interpretation in children is less significant. http://bit.ly/3cmRzsY
Introduction
Carbon monoxide transfer factor (TLCO) is a widely used test of respiratory function [1] and plays a vital role in the assessment of gas exchange and the diagnosis and management of various respiratory diseases [2–6]. The accurate interpretation of TLCO relies on the comparison with predicted values that are calculated from reference equations. However, the selection of appropriate reference equations can be problematic with at least 15 sets of equations published in the decade from 1995 to 2004 [7]. The European Respiratory Society (ERS)/American Thoracic Society (ATS) guidelines from 2005 [7] provide little guidance regarding the choice of TLCO reference equations, with the guidelines merely listing the most commonly used in North America and Europe.
In 2017 the Global Lung Function Initiative (GLI) published a new set of reference equations for TLCO [8]. These were derived from a very large normal data set (9710 subjects) compared to previous equations. The data was obtained from multiple international sites, with strict selection criteria. To be included in the analysis, the data needed to be collected after the year 2000, on modern equipment and adhere to strict quality control requirements. This large volume of data was used to create reference equations using the lambda, mu and sigma (LMS) method, which allows modelling of variability and skewedness of data, accounts for the interactive effects of age, height and sex and uses splines to model the non-linear age effects across the life span. The GLI reference equations also span the age range from 5 to 85 years, avoiding the need to use separate paediatric and adult equations, as well as reducing the need to extrapolate reference equations in the elderly.
The strengths of the GLI TLCO reference equations mean that they are most appropriate to use globally, which will lead to more consistent interpretation of TLCO results across centres. However, this transition could impact the interpretation of TLCO results in a clinical respiratory laboratory.
Such differences in the interpretation of spirometry have been well documented when those reference equations were updated [9–13]. In those studies, there were significant differences in the proportion of people classified as abnormal despite reasonable overall agreement across reference equations. To our knowledge, a similar analysis has not been performed for TLCO. The increased complexity in the measurement of TLCO compared with spirometry has the potential to create larger differences across reference equations, and hence a larger impact on TLCO interpretation.
The aim of this study was to document the effect of changing to the GLI TLCO reference equations from commonly utilised older reference equations on the interpretation of TLCO results in a large clinical dataset.
Materials and methods
Data were obtained from clinical pulmonary function laboratory databases at three Australian teaching hospitals (Austin Hospital in Victoria, Concord Repatriation General Hospital in New South Wales and The Royal Children's Hospital in Victoria). All three are large university-affiliated tertiary referral centres involved in the management of a broad range of respiratory diseases. Local Ethics Committee approval was obtained from each of the three hospitals.
A search was conducted on each of the hospital's databases for all TLCO results from 2008 until August 2018. Data were collected for patients aged 5 to 85 years. All testing was performed in accordance with the ERS/ATS guidelines [14]. Any test results not meeting the ERS/ATS guidelines were excluded.
The older reference equations used for comparison with the GLI equations for adults were those from Miller et al. [15], Roca et al. [16], Crapo and Morris [17] and European Community of Coal and Steel (ECCS) [18]. These reference equations were extrapolated to cover the range from 18 to 85 years. A separate analysis was performed with no age extrapolation (covering only the age range specified by each equation). For children, the comparisons with the GLI equations were made with the equations from Cotes et al. [19] and Kim et al. [20]. Data from both Caucasians and non-Caucasians were included in the analysis, with no race adjustment made for non-Caucasians.
The reference equations were used as published, with the exception of Roca, which contain a weight term for TLCO and carbon monoxide transfer coefficient (KCO) for females and alveolar volume (VA) and KCO for males. Due to the small weight range of the reference population, when these equations are applied to obese patients, the predicted values become non-physiological. A common practice to deal with this issue is to use the equations with a limitation on the maximum weight that is used to calculate the predicted values (limited to the maximum weight in the reference population). For example, the maximum weight in the reference population for females was 86 kg; so, for a female whose weight is above this value, a weight of 86 kg is used to calculate the predicted values.
Statistical analysis
The older reference equations were specifically for either adults or children and these data were analysed separately with patients <18 years of age assessed using the equations for children. The GLI equations that span the age range from 5 to 85 years were used for all patients. For each TLCO, KCO and VA result the lower limit of normal (LLN) was calculated using each of the relevant reference equations. A measured value below the LLN was considered abnormally low. The proportion of patients below LLN for each of the relevant reference equations, as well as the number of patients that changed from a normal to abnormal classification (or vice versa) was calculated. The level of raw agreement between each of the older equations compared with the GLI equations was calculated based on the percentage of patients that were classified in the same manner (either both normal or both abnormal) using both equations. The level of agreement was also assessed with the kappa statistic.
Results
Patient characteristics
The patient characteristics are shown in table 1. Data from 33 863 patients were available for analysis after excluding data (approximately 7% not meeting ERS/ATS acceptability criteria). Of these patients, 4028 (11.9%) were non-Caucasian. With the clinical indication for measurement of TLCO being less common in children, only 3% (n=930) of the patients were under 18 years of age. There was an even distribution of both sexes (48% female) and a wide distribution of height and weight in the dataset, with a tendency for the adult patients to be older and in the overweight range, based on BMI.
Summary of patient demographics and carbon monoxide transfer factor (TLCO) results for the entire patient group
Adults
For the adult equations, the median difference (IQR) in the TLCO LLN for GLI compared with the older equations was Crapo: −1.25 (−1.64−0.86), Roca: −1.02 (−1.33−0.67), ECCS: −0.35 (−0.62−0.07), Miller: 0.14 (−0.04–0.30) mmol·min−1·kPa−1. Figure 1 plots the mean LLN for each of the adult reference equations as a function of age, separated by sex. The GLI reference equations tend to produce a slightly lower LLN for a large portion of the age range for both sexes, particularly at a younger age, compared with all other equations. Above the age of approximately 70 years there is a tendency for the GLI equations to produce an LLN which is slightly higher than the LLN from the Miller and ECCS equations (fig. 1).
Mean carbon monoxide transfer factor (TLCO) lower limit of normal (LLN) for each of the adult reference equations as a function of age, separated into (a) females and (b) males.
The proportion of adult males and females with a reduced TLCO (below LLN) using each of the reference equations, as a function of age is illustrated in figure 2. Although the differences in the LLN were not large, they did result in altered rates of reduced TLCO for adults (Miller: 34.2% (n=11 249), GLI: 36.1% (n=11 902), ECCS: 43.2% (n=14 219), Roca: 58.8% (n=19 370), Crapo: 63.7% (n=20 991)). The largest difference in rates of a reduced TLCO for GLI equations occurs when compared with those of Crapo and Roca, particularly for females. There is also a large difference for younger females when comparing the GLI and ECCS equations. The closest agreement with the GLI equations for rates of reduced TLCO is with the Miller equations.
The proportion of (a) females and (b) males with a carbon monoxide transfer factor (TLCO) below the lower limit of normal (LLN) using each reference equations, as a function of age for adults.
The proportion of adult males and females who changed from an abnormally low TLCO result to within the normal range, when moving from each of the older equations to the GLI equations is illustrated in figure 3. Result classification changed (abnormal to normal) more frequently for adult females compared with males (for example reclassification occurred for 30% (n=4741) of females versus 16% (n=2779) of males when moving from Roca to GLI equations).
The proportion of (a) females and (b) males who changed from an abnormally low carbon monoxide transfer factor (TLCO) result to within the normal range, when moving from each of the older equations to the Global Lung Function Initiative (GLI) equations. LLN: lower limit of normal.
The largest proportion of patients changing category were younger rather than older adults. This pattern was most prominent for females when shifting from the Crapo and ECCS equations. As would be expected from the data in figure 2 there were a small number of older adults whose TLCO result changed from within the normal range to abnormally low when moving to the GLI equations from those of Miller and ECCS (supplementary figure S1).
Paediatric
For the paediatric equations, the median difference (IQR) for GLI in the TLCO LLN was Cotes: −0.53 (−0.63−0.39), Kim: 0.00 (−0.07–0.07) mmol·min−1·kPa−1. Figure 4 plots the mean LLN for each of the paediatric reference equations as a function of age, separated by sex. It can be seen that GLI reference equations, produce a slightly lower LLN for the entire age range for both sexes compared with the Cotes equations.
Mean carbon monoxide transfer factor (TLCO) lower limit of normal (LLN) for each of the paediatric reference equations as a function of age, separated into (a) females and (b) males.
Figure 5 demonstrates the proportion of males and females with an abnormally low TLCO using each reference equation, as a function of age for children. The rates of reduced TLCO for children were more similar across equations (GLI: 39.5% (n=367), Kim: 40.3% (n=375), Cotes: 53.7% (n=499)) than for adults. The GLI equations produced lower rates of a reduced TLCO for both sexes compared with those from Cotes (44.3% (n=224) versus 57.9% (n=293) for boys, 33.7% (n=143) versus 48.6% (n=206) for girls). The proportion of boys and girls who changed from an abnormally low TLCO result to within the normal range, when moving from each of the older equations to the GLI equations is illustrated in figure 6.
The proportion of females (a) and males (b) with an abnormally low carbon monoxide transfer factor (TLCO) using each of the different reference equations, as a function of age for children. LLN: lower limit of normal.
The proportion of (a) female and (b) male children who changed from an abnormally low carbon monoxide transfer factor (TLCO) result to within the normal range, when moving from each of the older equations to the Global Lung Function Initiative (GLI) equations. LLN: lower limit of normal.
Agreement between equations
Table 2 shows the level of agreement of the GLI equations with the older equations in identifying abnormality. The overall raw agreement varied from 72.3% with Crapo equations in adults to 98.1% with the Kim equations in children. For adults, the highest level of agreement with the GLI equations for both sexes is with the equations from Miller, whereas the lowest level of agreement is with the Crapo equations (table 2). The level of agreement between the GLI equations and those of Roca, Crapo, and ECCS is lower for females than it is for males. For children, the level of agreement between the GLI equations and those of Cotes is lower than the equations from Kim for both sexes.
Levels of agreement for a reduced carbon monoxide transfer factor (TLCO) (below lower limit of normal) comparing Global Lung Function Initiative equations to each of the older reference equations
In a separate analysis, the limits of agreement were calculated without any age extrapolation of the older equations. The results of this analysis were very similar to those seen with the extrapolated data (supplementary table S1).
Bland–Altman plots comparing the LLN for each of the older equations with the GLI are included in the supplementary material (supplementary figures S2 and S3 and supplementary table S2), showing the median difference (and 95% confidence intervals) between the older equations and GLI separated into males and females. These plots confirm that for adults the largest difference between equations occur when the LLN is larger, which would be seen in younger and taller people.
KCO and VA
The comparison with the GLI equations for adults shows that the LLN for KCO is similar to those of Miller and ECCS, however, they tend to be smaller compared to those of Roca and Crapo, particularly for females (supplementary figures S4 and S5, table S3). The difference was more pronounced at a younger age and higher KCO values. The comparison with the GLI equations for children shows that the LLN for KCO are similar to the Kim equations but the values tend to be smaller compared with the Cotes equations (supplementary figures S6 and S7, table S4). These differences in the LLN result in altered rates of an abnormal KCO with different equations, with the largest differences (and lowest agreements) occurring with the Crapo and Roca equations for adult females and the Cotes equations for children (supplementary figure S8, table S5).
The LLN for VA comparison with the GLI equations yielded mixed results (supplementary figures S9–S12, tables S3 and S4). Values are similar to those of Miller, and for the Kim and ECCS equations for females but not males. Although the median difference was not large for the Roca equations there was a larger spread. No predicted VA is available for the Crapo and Cotes equations. These differences in the LLN result in altered rates of an abnormal VA with different equations, with the largest differences (and lowest agreements) occurring with the ECCS equations for adult males and the Kim equations for male children (supplementary figure S13, table S6).
Discussion
Our analysis shows significant differences in TLCO interpretation using GLI reference equations compared with some commonly used older equations for adults. These differences in abnormality rates are smaller for children. Within these groups, the differences also alter with sex and age. The largest change in abnormality rate will be a reduction of approximately 33% for females if transitioning from the Crapo to the GLI equations. There will be a similar change for females if moving from the Roca equations. There will also be a reduction in abnormality rates (approximately 25%) for younger adult females when transitioning from the ECCS equations to the GLI equations and a similar change for younger adult males when moving from the Roca equations. The Bland–Altman plots comparing the TLCO LLN for older equations with GLI equations (supplementary figures S2 and S3) suggest that for any given TLCO, the change may be variable when switching equations and not easily predicted.
The recent publication of the GLI TLCO reference equations is likely to cause many pulmonary function laboratory directors to re-appraise their choice of reference equations. Changing reference equations is problematic for any clinical laboratory, with characterisation of disease presence and severity potentially altered. Given that our data are from three large hospital laboratories, we feel that they provide an important overview of the effects on interpretation in a representative clinical population.
The observed change in TLCO abnormality rates is not surprising. Similar to our study, a recent letter [21] identified a 1–21% increase in the number of patients who would qualify for clinical trials based on a TLCO cut-off of ≥30% predicted when using GLI compared with older equations, and the best agreement to be with the Miller equations. This strong agreement between the GLI and Miller equations is reassuring given that the Miller equations have previously been shown to be the best at predicting survival in a large group of patients [22]. This suggests that the GLI equations would be appropriate to use.
There are several explanations for the differences in abnormality rates between reference equations. Most of the older reference equations utilised data which were collected prior to 1993, before widespread standardisation of the TLCO test technique. There may have also been discordance in the actual values quoted (e.g. mean versus highest value) and most of the data were collected using equipment which may have had different performance characteristics compared with modern equipment. In addition, the GLI reference equations utilised data after correction for equipment deadspace and test altitude above mean sea level, and used more complex statistical analyses.
Another potential explanation is that the population used to generate the reference equations is not representative of the clinical population. While validation of reference equations is possible by comparing results from local normal subjects, the number required to identify any real differences may be up to 300 [23], which is not realistic for most respiratory function laboratories.
The levels of agreement between GLI equations and the older equations for classification of abnormality vary from 72.3% to 98.1%. These are lower than those seen with similar comparisons for spirometry equations [13]. This would be easily explained by the increased complexity of TLCO measurement, with the requirement to measure additional variables such as inspired and expired gas concentrations.
Overall for adults, adopting GLI will tend to reduce abnormality rates for adults, particularly up to the age of 70 years. This calls to question whether patients were falsely being identified as having a reduced TLCO using older equations. This may well be the case as there are several publications identifying an unexplained reduction in TLCO as a common clinical scenario [24–26]. The adoption of the GLI TLCO reference equations may result in fewer unnecessary clinical investigations following the identification of an isolated reduction in TLCO.
The difference in abnormality rates between the older equations and the GLI equations is most evident in young adults. One potential explanation for this is the use of linear models over the entire adult range in the older equations. The GLI data suggest that there is a plateau in TLCO from the age of approximately 18 to 25 years. Fitting a linear equation to this age range, when it is included with data from older adults may result in an overestimation of predicted values in young adults with the older equations.
In contrast, the GLI equations in both sexes above the age of 70 years produce higher rates of abnormality than the ECCS and Miller equations, with a small proportion of patients moving from within the normal range to abnormally low TLCO (supplementary figure S1). Survival data has suggested that the GLI spirometry equations overestimate predicted values in the elderly, when compared with extrapolated data from other equations [27]. This raises the possibility that the elderly subjects who are able to participate in the GLI normal values study may not be truly representative of the normal elderly population.
The use of KCO and VA in TLCO interpretation remains a controversial topic [7]. The levels of agreement for VA between the older prediction equations and GLI tend to be higher than the levels of agreement for KCO (supplementary tables S5 and S6). The levels of agreement for KCO tend to match well the levels of agreement for TLCO in the adult equations with the exception of the Crapo equations (table 2 and supplementary table S5). For children, there tends to be a difference in the levels of agreement between KCO and TLCO. The variability of agreement levels for these two parameters and the difference compared with the levels of agreement with TLCO is reflective of the fact that the prediction equations are produced independently for each of the parameters.
Although there may be concern about the differences in interpretation compared with older equations, the GLI equations have numerous advantages. First, the GLI equations cover a wider age range than the older equations, reducing the need to extrapolate equations beyond the age in which the data were collected. Secondly, use of the GLI equations eliminate the need to switch from paediatric to adult equations. The issues with switching equations at a certain age are well documented for spirometry [28] and are also relevant for TLCO. As previously mentioned, the other strength of the GLI equations is their scientific validity. All of these advantages suggest that the adoption of the GLI equations is likely to be widespread and rightly so.
The major weakness with the GLI TLCO compared with the GLI spirometry equations is the lack of race-specific equations for non-Caucasians. This poses a significant practical issue for most laboratories; however, this issue is also relevant to the older reference equations which are currently used. The current analysis did not apply an adjustment for non-Caucasians, given that the most commonly used adjustment [29] is based on small numbers, and a single race and sex. Interpreting TLCO results for non-Caucasians remains problematic and the creating race-specific TLCO equations should be a high priority.
There are some limitations to the current analysis. The most obvious is the assumption that the LLN is sensitive enough to separate normal from abnormal. We acknowledge that the certainty of interpretation is reduced, and likelihood of altered classification increases when results are close to the LLN. In these circumstances, the test results are best interpreted with the use of additional information such as other test results, the clinical picture and pre-test probability. Despite this limitation, the LLN is the value that underpins interpretative strategies [7] and we believe that our approach is justified in describing likely effects of changing reference equations. Our analysis did not apply the Roca reference equations exactly as published, which may be considered a limitation. However, our approach of applying a weight limit to the Roca equations [16] prevented extrapolated non-physiological reference values. Without this weight limit, the difference in abnormality rates for females between the Roca and GLI equations would be even larger.
The 2017 GLI TLCO equations provide a unique opportunity to enable accurate prediction of normal values, and their widespread uptake into clinical practice appears likely. Our analysis provides an all-age summary of changes that can be expected, and we believe this information will assist in adoption of these reference equations, facilitating further standardisation in the interpretation of this valuable diagnostic tool.
Supplementary material
Supplementary Material
Please note: supplementary material is not edited by the Editorial Office, and is uploaded as it has been supplied by the author.
Supplementary material ERJ-01905-2019.SUPPLEMENT
Shareable PDF
Supplementary Material
This one-page PDF can be shared freely online.
Shareable PDF ERJ-01905-2019.Shareable
Footnotes
This article has supplementary material available from erj.ersjournals.com
Conflict of interest: D.J. Brazzale has nothing to disclose.
Conflict of interest: L.M. Seccombe has nothing to disclose.
Conflict of interest: L. Welsh has nothing to disclose.
Conflict of interest: C.J. Lanteri has nothing to disclose.
Conflict of interest: C.S. Farah has nothing to disclose.
Conflict of interest: W.R. Ruehland has nothing to disclose.
- Received September 26, 2019.
- Accepted February 20, 2020.
- Copyright ©ERS 2020