Abstract
There is little literature about the mortality associated with bronchiectasis. The aim of the present study was to investigate factors affecting mortality in patients with bronchiectasis.
In total, 91 patients were examined for aetiology, pulmonary function tests, high-resolution computed tomography, sputum microbiology and quality of life scores and were then followed over 13 yrs.
Overall, 29.7% of the patients died. On multivariate analysis, age, St George’s Respiratory Questionnaire activity score, Pseudomonas aeruginosa infection, total lung capacity (TLC), residual volume/TLC and the transfer factor coefficient were all independently associated with mortality.
In patients with moderate to severe bronchiectasis, mortality is associated with a degree of restrictive and obstructive disease, poor gas transfer and chronic pseudomonas infection. These features should guide future research into disease progression, and identify those patients needing intensive treatment.
Bronchiectasis is a chronic respiratory condition characterised by the abnormal dilatation of the bronchial lumen 1. Damaged airways predispose the patient to recurrent bacterial infections, which in more severe cases become chronic. This leads to a host inflammatory response that causes tissue damage, leading to a vicious circle of host-mediated, bacteria-stimulated disease progression 2. However, why in the majority of patients such progression is slow, whereas in others it occurs much more quickly is poorly understood.
Investigation of bronchiectasis aims to identify possible aetiologies and predisposing factors that may be treatable 3. Management then centres on physiotherapy and antibiotics are used promptly in infective exacerbations. There is a lack of good randomised controlled trials for evidence-based management in bronchiectasis, and many of the management programmes are extrapolated from the treatment of other conditions such as pneumonia and cystic fibrosis 4.
A clearer picture of the natural history and mortality in bronchiectasis would be of great value, both in allowing more accurate prognostic information and in determining the specific characteristics of the disease that confer an increased risk of progression and mortality. This information will form the basis of future studies assessing the impact of therapies in bronchiectasis. Neutrophilic inflammation is thought to be a key component of the pathogenesis of chronic obstructive pulmonary disease (COPD), a condition in which infective exacerbations have been linked to disease progression 5, 6 and in which bronchiectasis is quite commonly seen in severe cases on computed tomography (CT) scans 7. A better understanding of disease mechanisms in bronchiectasis might also improve COPD research.
In 1994, 111 patients with bronchiectasis were studied to validate the St George’s Respiratory Questionnaire (SGRQ) as a tool in assessing the health status of patients with bronchiectasis 8. In the current study, 14 yrs on, we have reassessed these patients and carried out an analysis of the factors in the original assessment that predict mortality. The purpose of the present study was to provide prognostic information at baseline.
METHODS
In the original study in 1994, 120 sequential patients with clinically diagnosed bronchiectasis were approached to take part in this outpatient study in the Host Defence Unit of the Royal Brompton Hospital (London, UK), a tertiary referral centre. Of these, 111 agreed to take part in the study. The exclusion factors were bronchiectasis not being the predominant pathology and cystic fibrosis; however, no patients were in these groups 8. Each patient provided signed consent and the protocols were approved by the Royal Brompton Hospital Ethics Committee. Each patient had full host defence investigations, which have been described elsewhere 9, including sweat testing for cystic fibrosis, two sputum samples that were also tested for nontuberculous mycobacteria, and a high-resolution CT scan. Comprehensive lung function tests with body plethysmography were performed on two separate occasions 6 months apart at entry to the study.
The outcome of these patients was determined as of March 2007, 13 yrs after the initial recruitment. Patients not followed up in our institution were tracked via general practitioners, other hospitals and the Office for National Statistics, and death certificates were obtained for all patients that died. The CT scans were reassessed and on this basis the patient number was reduced further to only include patients who also had bronchiectasis according to CT criteria.
CT images of 1–1.5 mm thickness at 10 mm intervals were scored for bronchiectasis 10, 11 by two blinded, independent radiologists. The scoring system used was modified from the system of Bhalla et al. 12 and has been shown to be associated with low inter-observer variation 13. Bronchiectasis extent and bronchial wall thickness were scored in each lobe on a scale of 0 to 3 (3: most severe changes). Bronchial dilatation and mucus plugging were graded from 0 to 2. A mosaic attenuation pattern (reflecting small airways obliteration) and emphysema were scored to the nearest 5%. The average score of the two observers, as a percentage of the maximum score possible to account for patients with lobectomies, was used for analysis.
Statistical analysis
Kaplan–Meier curves were used to illustrate survival data. Variables were examined against survival using Cox’s proportional hazards model: in multivariate models, a stepwise approach was used to discard variables not independently linked to mortality (p>0.05). Hazard ratios are given with confidence intervals (CIs) for significant (p<0.05) and marginal (0.05<p<0.15) trends.
Stepwise multiple linear regression was used to identify independent CT determinants of lung function impairment. Lung function indices shown to have an independent impact on mortality (total lung capacity (TLC), residual volume (RV)/TLC and transfer coefficient of the lung for carbon monoxide (KCO)) were used as the dependent variables in separate models. The validity of the parametric assumptions of linear regression was confirmed with testing for heteroscedasticity and RV/TLC was subjected to zero skewness logarithmic transformation on the basis of this test.
Statistical analyses were performed using STATA (version 4; StataCorp LP, College Station, TX, USA).
RESULTS
Patient demographics and survival
There were 111 patients that consented to the original SGRQ study in 1994. Review of CTs found no bronchiectasis in 12 patients and these were not included in further analysis. These patients had a clinical presentation compatible with bronchiectasis, but they did not satisfy the CT criteria. The data set was incomplete in a further eight patients and they were also excluded, leaving 91 patients in the study (92% of the eligible patients; fig. 1⇓). Patient demographics are shown in table 1⇓.
Patients were followed up from recruitment in 1994, until March 2007. During the 13-yr follow-up, 29.7% of the bronchiectasis patients died. Using life expectancy data from the Office for National Statistics, the expected death rate of a 52-yr-old (mean age in our study) is 14.7% for males and 8.9% females over a similar 13-yr period 14. Accounting for the mortality in our study, the mean±sd length of follow-up was 126.9±38.0 months. All the 91 patients were successfully accounted for at the end of the study. The cause of death was bronchiectasis, respiratory infection or respiratory failure in 19 (70.4%) out of 27 cases. Haemoptysis was the terminal event in one of these patients. Two patients each died of renal failure and colon cancer, and one patient each died of heart failure, cerebrovascular accident, liver metastasis and pulmonary embolism. Post mortem examinations were performed in three patients and the cause of death was classified as suppurative bronchitis, bronchiectasis and heart failure. Of those who died during the study, the median age was 60 yrs. A Kaplan–Meier survival plot is shown in figure 2⇓.
The aetiology of the bronchiectasis was as defined at recruitment and was similar to previous studies 9, 15. Idiopathic aetiology accounted for more cases (56.0%) than a recent study (34.6%) 9. Nevertheless, the same broad trends were observed, with idiopathic and post-infective causes of bronchiectasis being the most common aetiologies. The majority of the patients in the current study were never-smokers (76.9%) and the proportion of emphysema on the high-resolution CT scans was also low (mean score <4%). There was a range of severity scores (table 1⇑), suggesting that the present study population covered a spectrum of disease activity; however, broadly the cohort described a population with moderate to severe bronchiectasis. Bronchiectasis affected a total of 304 (56.8%) out of a possible 535 lobes examined on CT scan. Spirometry showed that airflow obstruction was the predominant finding in this group of bronchiectasis patients, with an average forced expiratory volume in 1 s (FEV1) of 1.85 L·min−1 (65.8% predicted) and forced vital capacity of 3.05 L (88.7% pred), which is in accordance with other studies 15–17.
Mortality risk factors
All the measured variables were included in univariate analysis, as shown in table 2⇓. A relative rate (RR) >1 describes a positive correlation between the parameter and mortality, whereas RR <1 suggests a negative correlation. There were large ranges of lung function, CT and health status variables that had a strong relationship to mortality in the current analysis. There were particularly strong associations for bronchiectasis extent, bronchiectasis severity and wall thickness on CT assessment; FEV1 and transfer factor of the lung for carbon monoxide on lung function measurement; and dyspnoea and activity score on the health status assessment. The RRs and CIs for the CT and lung function criteria appear much lower than, for example, the existence of pseudomonas. This is a function of them being continuous as opposed to categorical variables. For lung function tests, a risk ratio reflects a change in risk per unit lung function. In other words, a risk ratio of 0.97 for FEV1 means that for every 1% change in FEV1 there is a 3% change in risk. A difference of 30% FEV1 consequently confers a relative rate of 2.71 (95% CI 1.65–4.45).
Several of these variables will represent different ways of measuring the same effect. The multivariate analysis demonstrated the variables with independent effects on mortality. Complex CT scoring systems are not practical for routine clinical practice and would not be available in most centres, and thus these variables were excluded from the multivariate analysis. Increased age, male sex, Pseudomonas aeruginosa colonisation and higher SGRQ activity scores (describing lower activity levels), in addition to higher RV/TLC and lower TLC, and KCO lung function measurements were independent predictors of mortality (table 3⇓).
Although not used in the main multivariate model, the CT variables were assessed together in a separate multivariate analysis to find the most important CT factors that had strong independent effects on mortality. Increased wall thickness (RR 1.06, 95% CI 1.03–1.10; p<0.005) and emphysema (RR 1.04, 95% CI 1.00–1.07; p = 0.027), despite the low prevalence of this latter feature, were the strongest predictors of mortality in bronchiectasis from the CT findings.
Seven patients had lobectomies and four of these died. The mean±sd TLC in the lobectomy group was 4.63±0.91 L, which was 80.43±10.97% pred, compared with 6.00±1.42 L (105.04±15.55% pred) in the other patients. However, the lung function parameters that predicted mortality in the survival analysis were unchanged when these seven patients were removed.
There was no difference in mortality rates among the different aetiologies of bronchiectasis. As the number in several of the groups was very small, we separated them into idiopathic and known causes (p = 0.57 Cox proportional hazard) and have illustrated mortality differences using Kaplan–Meier curves (fig. 2⇑).
Functional–morphological correlations
The lung function measures that independently predict survival in bronchiectasis are TLC, RV/TLC and KCO, as shown previously. The morphological correlates of these were examined by separate linear regression models with these three lung function measures as the dependent variables and the CT variables as cofactors. A regression coefficient (RC) >0 describes a positive correlation between the parameter measurements, whereas RC <0 suggests a negative correlation. TLC was positively correlated to mosaicism on the CT scan, RV/TLC was positively correlated with wall thickness and mosaicism, and KCO was negatively correlated with emphysema (table 4⇓).
Longitudinal data
The purpose of the current study was to construct prognostic indices from baseline data. We did not routinely repeat these tests during the course of the study. Nevertheless, the majority of patients did have several repeat lung function and sputum investigations during follow-up. There was an average decrease in KCO per year of follow-up of 0.03 (95% CI 0.01–0.04). This was significantly greater in those who died compared with the surviving patients (0.05 (0.005–0.09) versus 0.02 (0.004–0.03), respectively; p = 0.02, Mann–Whitney). The longitudinal changes in TLC (0.02 (-0.02–0.06) yr−1) and RV/TLC (0.1% (-0.62–0.82%) yr−1) were insignificant, with no difference between the patients that survived or died. Repeat sputum cultures were also performed in 17 out of the 20 patients with P. aeruginosa. Of these, P. aeruginosa was a consistent finding in 11 out of 17 patients with intermittent cultures in six out of 17 patients.
DISCUSSION
Prior to the present study, the outcome in adult bronchiectasis had not been well described and there has been a lack of long-term prospective studies. The current study is unique in having a very well defined bronchiectatic population, all of whom had extensive measurements of disease at baseline. Furthermore, all patients were followed prospectively until death or the end of the study 13 yrs later. The extensive nature of the testing in the present study and the length of follow-up mean that it is ideally placed to comment on features of bronchiectasis patients that influence mortality.
The average age of death in patients with bronchiectasis was found to be <55 yrs in a study in the 1960s 18, and most patients succumbed to the disease before the age of 40 yrs in a report of 400 bronchiectatics in 1940 19. Since this time, however, the condition, aetiologies and treatment have altered with the advent of vaccination programmes and the increased use of antimicrobials for childhood infections and tuberculosis. More recent studies have reported a variation in mortality rates in bronchiectasis, with the worst being a 4-yr survival of 58% 16, a retrospective study showing a 75% survival at 8.8 yrs 20, and the best being an 81% survival at 14 yrs 17. The survival rates during our study were broadly comparable to the latter two studies, being 91% at 4 yrs, 83.5% at 8.8 yrs and 68.3% at the end of 12.3 yrs (fig. 2⇑). The primary cause of death in our patients was respiratory, suggesting that patients die from, not just with, bronchiectasis.
The current study looked specifically at the factors associated with mortality rate in bronchiectasis over a 13-yr period. There are few previous studies assessing this, but a poor initial ventilatory capacity 17, male sex 20, and a secondary diagnosis such as asthma or COPD 20 have all been suggested as inferring a worse prognosis in bronchiectasis. A more recent report prospectively followed 98 patients with bronchiectasis over a 4-yr period in Turkey. Risk factors for mortality in this study were only looked at in a univariate analysis, and included age, low body mass index, hypoxia, hypercapnia, radiographic extent, Medical Research Council dyspnoea scale and lack of vaccinations 16. The results of our study have demonstrated that the factors with an independent impact on mortality are age, SGRQ activity score, P. aeruginosa colonisation, and TLC, RV/TLC and KCO lung function measurements.
Among a larger number of variables examined against mortality, most were not independently predictive and were discarded using a stepwise methodology, to develop a final model that contained only six covariates. Failure to discard nonsignificant variables would have resulted in confounding by collinearity (as, for example, between pulmonary function measures of airflow obstruction), as well as a degree of model over-fitting. However, it should also be acknowledged that in multivariate analysis, the retention of one collinear variable (quantifying, for example, airflow obstruction), as opposed to another, is sometimes a very close call. The multivariate modelling strategy used is predominantly for enabling mechanistic conclusions as opposed to predicting individual patient outcomes. The results should be interpreted as identifying the broad determinants of mortality (which includes the severity of airflow obstruction), in the search for pathogenetic insights. The model is not advanced as a multivariate clinical index for immediate use. Nevertheless, the stronger univariate relationships we report are of immediate clinical relevance.
Formal CT scoring, as performed in our study, is unrealistic in routine clinical practice 21 and was not included in final multivariate modelling. This approach was further validated by the finding that CT features did not reach statistical significance for an independent effect on mortality when re-entered into the modelling equations. The CT variables are nevertheless very useful in order to appreciate the pathogenetic correlates of the functional indices. The CT scan correlates of the lung function parameters with an independent impact on mortality, TLC, RV/TLC and KCO were mosaicism, wall thickness and mosaicism, and emphysema, respectively.
Of the variables shown to have an independent effect on mortality in bronchiectasis, P. aeruginosa has previously been shown to infect patients with more extensive disease and severe airflow obstruction 22. Two studies have shown that infection leads to an increased progression of disease 23 with a more rapid decline in lung function parameters 24. The cause and effect relationship is less clear, as this pathogen may instead be a marker of disease progression as opposed to the cause for progression 25. The independent effect of infection with P. aeruginosa on mortality shown in the present study suggests that it is likely to impact on survival more than just being a marker of severity. Infection control measures, such as clinic segregation, to avoid patient exposure to the bacterium, should be investigated. Furthermore, treatment regimes designed to eradicate P. aeruginosa, when it occurs for the first time, have shown some benefit in delaying the decline of lung function in patients with cystic fibrosis 26, 27, and may similarly impact on the progression of bronchiectasis 4. This should be investigated in randomised controlled studies.
Poor lung function measurements have previously been shown to be associated with mortality in bronchiectasis 17. Our study has shown an increased mortality in patients with a high RV/TLC ratio. This suggests the importance of obstruction in the presence of restriction in this condition. In the multivariate model, RV/TLC was more important than the other variables that measured obstruction, suggesting that obstruction had the greatest effect on mortality in the presence of some restriction. The importance of a component of restrictive lung disease in bronchiectasis is also confirmed by the correlation between mortality and a lower TLC in the present study. This relationship was still present when those patients who had lobectomies were removed. The final lung function measurement with an independent effect on mortality in this study is KCO. This demonstrates a vascular or emphysematous component that impacts on mortality in bronchiectasis. The low prevalence of emphysema in CT scans suggests the former as the more likely. There appear to be separate and independent effects of the three individual components of functional impairment on lung function testing, with obstruction, restriction and impairment of gas transfer all separately implicated in mortality. These components are likely to represent separate pathogenic components of bronchiectasis that may be important in helping to understand the impact of disease, in addition to providing specific targets for future treatments. To elucidate the potential morphological correlates of these lung function measures, they were correlated with CT findings.
The measure of obstruction, RV/TLC, was correlated to mosaicism and wall thickness on CT scans. The correlation of functional obstruction to wall thickness does not differentiate between secretions and ablative bronchiolitis as the cause; however, the relationship with mosaicism favours the latter as a morphological explanation of the obstructive functional element in bronchiectasis. This is in agreement with a previous study correlating CT findings with FEV1, which is another measure of airways obstruction 28. The authors concluded that airways obstruction in bronchiectasis was primarily related to intrinsic disease of small and medium airways, and not to the large airway abnormalities, airway mucus plugging or emphysema. Small airways disease is prominent in histological studies of bronchiectasis, with obstructive and inflammatory bronchiolitis leading to the features of airways obstruction and decreased radiographic attenuation 29.
In contrast, CT features could not explain the restrictive element, with TLC showing only a positive correlation with mosaicism. Possible pathological correlates of lung restriction in bronchiectasis include parenchymal scarring secondary to previous infections, pleural disease, peribronchial fibrosis and atelectasis. Peribronchial fibrosis has been demonstrated in histological studies of bronchiectasis 30, and this together with atelectasis distal to the ablated airways may not be detectable on CT scanning 28.
The contribution of poor gas exchange to mortality may be secondary to emphysema or to a pulmonary vascular component. The KCO was positively correlated to emphysema and not mosaicism, suggesting emphysema as the morphological correlate of low transfer factor. However, CT measurements of the pulmonary vasculature are poor, and the possible importance of pulmonary hypertension in bronchiectasis is further highlighted by an earlier study in which cor pulmonale was present in 37% of patients who died from bronchiectasis 17. From a management point of view, the importance of assessing pulmonary artery pressures in bronchiectasis patients with low transfer factors should be emphasised, as these patients may be amenable to some of the newer anti-hypertensive agents.
In addition to helping understand the pathogenesis of the condition and providing additional targets for treatment, knowledge of the lung function factors associated with mortality may provide useful correlates for mortality against which future treatments can be tested, and it may be used in defining subgroups of patients into prognostic categories. The severity of bronchiectasis is presently defined by the extent and grade of large airway dilatation. These features are not associated with prognosis and hence a paradigm involving lung function may be more appropriate. The current study was not designed to investigate the rate of change of parameters over time and hence routine repeat measurements were not performed. The longitudinal investigations we have described were taken at different times and for different indications. With these limitations, a decline in KCO was demonstrated over time, predominantly in the patients that died, and this may be helpful in providing a basis for assessing the stability of patients in the long term.
In summary, the present study has suggested the importance of lung function measurements in addition to age, P. aeruginosa infection and health questionnaire activity scores in the prognosis of bronchiectasis. It is important to note that the conclusions refer to a tertiary centre and are subject to the limitation of a predominantly retrospective analysis. Nevertheless, knowledge of the results from the current study should change investigation protocols and define patients that need more intensive management plans, such as antibiotics or anti-inflammatory agents. The results may also help to define useful subgroups for the assessment of treatment regimes in much needed randomised controlled trials, to help advance bronchiectasis treatment towards a stronger evidence base.
Statement of interest
None declared.
Acknowledgments
The current authors would like to acknowledge P. Jones (St George’s Hospital, London, UK) and C. Wilson (no affiliation; previously Royal Brompton Hospital, London) for their major contribution to the initial investigation of the cohort of patients, and their helpful comments on the present work.
- Received January 8, 2009.
- Accepted March 31, 2009.
- © ERS Journals Ltd