Abstract
After an experimental neonatal screening program for cystic fibrosis (CF) from 1973–1979, a follow-up study took place from 1980–1997. Patients were treated at specialized centres (C) or at local hospitals (non-C). Aims of the study were: 1) to determine whether the previously reported benefits from screening persisted with time and after adjustment for confounding variables; and 2) to investigate whether centre treatment was associated with improved prognosis of CF patients.
Prognosis of patients detected by screening (S; n=24) was compared with patients detected clinically, born during (non-S; n=29) and after the screening programme (post-S; n=39). In addition, prognosis was compared between 45 C and 47 non-C patients. Multivariable regression analysis was used to compare survival and mixed-effects model regression analysis was used to compare clinical outcome between patients. The analyses included the variables screening, centre treatment, sex, meconium ileus and genotype.
S patients had a significantly smaller decline in forced expiratory volume in one second (FEV1) (difference +2.74% predicted) and significantly lower immunoglobulin-G (IgG) levels (difference −473.69 mg·dL−1) than non-S patients until 12 yrs of age. At 12 yrs of age, vital capacity was significantly higher in S patients than in non-S patients (difference +362.79 mL). Survival seemed to be best for S patients compared to both non-S and post-S patients. Post-S patients were significantly heavier (difference in sd weight +0.77), had a significantly smaller decline in FEV1 (difference +2.80% pred) and lower IgG levels (difference −453.04 mg·dL−1) than non-S patients until 12 yrs of age. C patients had a significantly improved survival (relative risk (RR) 0.18, 95% confidence interval 0.05–0.57) than non-C patients.
Early diagnosis through neonatal screening leads to better preservation of lung function in the long term in cystic fibrosis patients. Management of cystic fibrosis patients in specialized centres improves survival.
This study was supported by ZorgOnderzoekNederland (ZON) grants 28.566, 28.566-1 and 28.566-2.
Pulmonary disease is the major cause of morbidity and mortality in cystic fibrosis (CF). Shortly after birth, the lower airways of most CF patients become colonized with bacteria which leads to lung infection and inflammation 1. Bacterial lower respiratory tract infections in CF are associated with significant pulmonary deterioration but are amenable to antibiotic therapy at an early stage, resulting in clearance of pathogens and reduction of inflammation 2. Therefore, early diagnosis may favourably influence prognosis of patients with CF as a result of early intervention, which may prevent or postpone early irreversible damage of the airways.
Any proof of the medical benefit of an early CF diagnosis through neonatal screening requires studies with many years of follow-up, given the relatively long life span of patients with CF. Consequently, patients with CF may have been denied an effective form of intervention for a long time while awaiting conclusive results of randomized, controlled trials of neonatal screening for CF 3–5. The effectiveness of neonatal screening for CF can also be assessed by nonrandomized, controlled studies with a long follow-up period, focusing on important clinical outcome measures such as pulmonary status 6.
In the present 17-yr follow-up study of patients with CF, earlier observations were extended over a longer period and a greater number of patients. In contrast to previous reports 7, 8, data in this study were analysed on an intention-to-screen basis and analyses were adjusted for possible confounding variables, such as centre treatment, sex, meconium ileus, and genotype. CF care is probably best provided in specialized centres because of the all-round management, although few data exist to prove this assumption.
The aims of this study were: 1) to determine whether the previously reported benefits from neonatal screening 7, 8 persist with a longer duration of follow-up and after adjustment for several potential biases; and 2) to investigate whether centralized treatment was associated with improved prognosis of patients with CF.
Methods
Study design
The design of this study has been described in detail elsewhere 8. In summary, a clinical follow-up study of patients with CF was performed from April 1980–April 1997, after an experimental neonatal screening programme for CF. Screening for CF was carried out in a geographically defined area in the north of the Netherlands from March 1973–March 1979, using the determination of the albumin content of meconium 9. Of all neonates born in this area, 45% were involved in the screening programme. The organizations caring for mother and child decided whether to participate in the screening programme based on their willingness to perform a screening test in newborns. These organizations operate quite independently from hospital care in the Netherlands. Outcome of patients identified by screening (S) was compared with two cohorts of patients diagnosed clinically: one cohort born during the screening programme (non-S) and one cohort born in the 6 yrs after (post-S) the screening programme (table 1⇓). Post-S patients were added to increase the number of patients in this follow-up study. Previous results of this follow-up study have been described elsewhere 7, 8. From 1980–1990, clinical data were collected prospectively. Data from 1990–1997 were retrospectively collected.
After confirmation of diagnosis, patients were treated either at or in close cooperation with a CF centre (C) or solely at a local hospital (non-C). In CF centres, teams including paediatric pulmonologists and gastroenterologists, dieticians, physiotherapists, CF nurses, psychologists, microbiologists, and geneticists, all experienced in treating patients with CF, provide all-round management of CF patients based on predefined protocols with regular follow-up. In addition, CF centres perform research on different aspects of CF. In local hospitals, CF care is provided by less disciplines without specialized CF clinics, and fewer CF patients are treated at local hospitals.
Allocation of treatment was left to the decision of the general paediatrician taking care of the patients at the time of diagnosis. The parents of the patients sometimes requested referral to a CF centre. CF patients who were seen once or twice a year at a CF centre and treated at a local hospital were considered to be treated in close cooperation with a specialized centre. The University Hospital Vrije Universiteit Amsterdam and the University Hospital Groningen ethics committees approved this study and all subjects gave informed consent.
Patients
In the present study, 92 patients were identified (table 1⇑). Before the start of the follow-up, 10 patients died: four S patients, five non-S patients and one post-S patient. Seventy per cent of these patients were treated solely at local hospitals. Three post-S patients treated at a specialized centre were excluded from the study because of comorbidity, including congenital heart disease, spastic tetraplegia, and severe bowel dysfunction leading to early death. During follow-up, three non-C treated patients were lost to follow-up: two S patients and one post-S patient. After 1990, 12 CF patients could be added to the present analysis. Seven patients were identified with the help of the Dutch Cystic Fibrosis Registry (four non-S and three post-S patients, 57% were treated at a CF centre) and five deceased patients (all post-S) were identified with the help of Statistics Netherlands in Voorburg. Data were retrospectively collected from the first seven patients. It was assumed that the last five patients received care at local hospitals since they were not known in the CF centres in this study. No clinical data were available from these five patients.
Diagnosis of CF in all patients was confirmed by a positive sweat test (sweat chloride concentration >60 mmol·L−1) performed by quantitative pilocarpine iontophoresis. Furthermore, in most patients deoxyribonucleic acid (DNA)-analysis was performed to determine the mutations involved. CF genotypes were determined by the Dept of Medical Genetics, University of Groningen. No data on mutation analysis were available from the patients who died before 1990.
Forty-five (49%) of the 92 patients were treated at or in close cooperation with a specialized centre. In the clinical follow-up, 72 patients could be included (table 1⇑) of whom 39 (54%) received centralized treatment (table 2⇓).
Outcome measures
Survival
Patient survival was calculated as the period from birth to the end of the observation period on 1 April 1997. Patients with meconium ileus were excluded from survival analysis, they have an increased risk of dying in the first month of life, and most patients with meconium ileus were referred to a CF centre. Patients who were lost to follow-up and patients who died before the start of the study were included in the analysis until the date they were last known to be alive.
Clinical outcome
Clinical data were obtained every 6 months. For nutritional status, anthropometric variables were converted to sd scores for height and weight, by comparison with the 50th percentile for mean height and weight of Dutch children of the same age and sex 10. Patients were classified as having pancreatic insufficiency or sufficiency according to the results of the faecal fat excretion as assessed by the Van de Kamer method and/or the presence or absence of symptoms of malabsorption responding to pancreatic enzyme supplementation. Lung function was measured using a wet 10 L spirometer (Lode Spirograph Type DL; Lode Instruments, Groningen, the Netherlands), or a pneumotachograph (Jaeger, Masterlab, Wuerzburg, Germany), taking the best of three manoeuvres. Inspiratory vital capacity (VC) and forced expiratory volume in one second (FEV1), expressed as a percentage of reference values based on patient height, age and sex 11, were used as parameters.
Immunoglobulins were determined by standard methods. Immunoglobulin-G (IgG) levels were used as an indicator of the severity of chronic lung infection and inflammation.
Statistical analysis
Patients with a false-negative screening test were included in the S group in order to perform an intention-to-screen analysis, in contrast to the previously reported analysis of this follow-up study 8. The Chi-squared test was used to compare patient groups with respect to categorical variables. The main study interest was to examine the influence of screening on prognosis of patients with CF. Because of differences in patient management, which may also influence outcome of CF patients, a separate analysis of the effect of centre treatment on prognosis of patients was performed. Moreover, analyses were adjusted for differences in sex and meconium ileus between patients.
Survival
Multivariable regression analysis, based on the Cox's proportional hazard model, 12 was used to compare estimated survival between patients including the predictor variables screening, centre treatment, age at diagnosis, and sex.
Clinical outcome
Clinical outcome between patients was compared using linear mixed-effects (LME) model regression analysis 13. The LME analysis implemented in S-Plus 14 accounts for irregularly spaced longitudinal data and different follow-up times. Screening, centre treatment, sex, and meconium ileus were used as predictor variables. The analysis of VC was also adjusted for differences in height between patients. Interaction terms with age were included to test differences in change over time of the outcome variable for various subgroups. Separate models for each clinical outcome measure, including sd scores for height and weight, FEV1 % predicted, VC, and IgG, were fitted.
At the end of the follow-up in 1997, the eldest patients were 24 yrs of age (S and non-S patients) and the youngest patients were 12 yrs of age (post-S patients). The analysis of the clinical follow-up data was, therefore, divided into two periods: clinical outcome until, and above, 12 yrs of age. In the analysis of outcome until 12 yrs of age, the intercept or starting point was taken at the age of the first available lung function data, mostly around 6 yrs.
To explore a possible influence of genotype on clinical outcome, separate LME analyses including the variable genotype were performed.
The data are presented as differences in intercepts and slopes between patient groups. In all analyses, the significance level was taken at p<0.05.
Results
During the follow-up, 24 screened and 68 nonscreened CF patients were identified (table 1⇑). Frequency of CF at birth was similar in the S and non-S groups, but was lower in the post-S group. Age at diagnosis and number of patients with meconium ileus was highest in the post-S group.
Survival
Analysis of survival was performed in 74 patients, excluding 18 patients with meconium ileus or comorbidity.
Influence of screening
After adjustment for age at diagnosis, sex, and centre treatment, the S group seemed to have a better survival throughout follow-up than the non-S and post-S groups (fig. 1⇓). The relative risk (RR) of dying for the S group was 0.43 (95% confidence interval (CI) 0.13–1.38, p=0.16) in comparison with the non-S group, and 0.35 (95% CI 0.09–1.41, p=0.14) in comparison with the post-S group. As compared with the non-S group, the RR of dying for the post-S group was 1.23 (95% CI 0.42–3.60, p=0.70).
Influence of centre treatment
After adjustment for screening, age at diagnosis, and sex, the C group had a significantly better survival during total follow-up than the non-C group (fig. 2⇓). The RR of dying for patients from the C group compared with patients from the non-C group was 0.18 (95% CI 0.05–0.57, p<0.01).
Clinical outcome
Longitudinal clinical evaluations were performed in 72 patients: 18 S, 24 non-S, and 30 post-S patients (table 2⇑). From the 72 patients, 39 (54%) received centre care. Patients with meconium ileus were included in this analysis. Characteristics of patients participating in the longitudinal analysis are compared in table 2⇑. The number of patients treated at a specialized centre was significantly higher in the S group than in the non-S group. The post-S group received more centre care than the S group.
Observed values throughout total follow-up from FEV1 % pred in the S, non-S and post-S patients are shown in figure 3⇓. In addition, observed values from sd weight and FEV1 % pred in the C and non-C patients are presented in figures 4 and 5⇓⇓. In these figures, no adjustment is made for the influence of the several predictor variables used in the LME analyses, such as screening, centre treatment, sex, and meconium ileus.
Influence of screening
Analysis of outcome until 12 yrs of age
Until the age of 12 yrs, S patients showed a tendency towards better sd scores for height and weight than non-S patients. Difference in sd weight at the age of 6 yrs between S and non-S groups was 0.52, p=0.06 (table 3⇓). As a result of a higher increase of height over time in the non-S group, differences in height at the end of this period between these groups were much smaller. No significant differences in nutritional status until the age of 12 yrs were found between the S and post-S groups. However, in comparison with the non-S group, the post-S group was significantly heavier at the age of 6 yrs (difference in intercept of sd weight: 0.77, p<0.01; table 3⇓).
At the age of 6 yrs, the S group seemed to have better lung function parameters than the non-S and post-S groups. The post-S group showed the worst lung function at this age.
As expected, all patients showed a decline in FEV1 % pred throughout the observation period. However, the change in FEV1 over time was significantly different between the S and non-S groups, favouring the S group (difference in slope 2.74% pred, p<0.05), and between the post-S and non-S groups, favouring the post-S group (difference in slope 2.80% pred, p<0.05) (table 3⇑). The change in FEV1 over time in the S group was not different from that in the post-S group. The S and post-S groups tended to have a higher increase in VC over time than the non-S group.
In this analysis, the S group had significantly lower IgG levels than the non-S group (difference in intercept −473.69 mg·dL−1, p<0.01; table 3⇑). Levels of IgG were similar in patients from the S and post-S groups, but the post-S group had significantly lower IgG levels than the non-S group (difference in intercept −453.04 mg·dL−1, p<0.01; table 3⇑).
Analysis of outcome above 12 yrs of age
After the age of 18, final height is expected to be achieved among patients in the follow-up. Therefore, the analysis of nutritional status after the age of 12 yrs was performed until the age of 18 yrs. Differences in height and weight between the S and non-S groups were smaller in patients >12 yrs of age than <12 yrs of age and were statistically not significant. In addition, no significant differences in nutritional status between either the S and post-S groups or the post-S and non-S groups were found.
The S group tended to have overall better lung function parameters than the non-S and post-S groups. Average values of VC at the age of 12 yrs in the S group were significantly higher than in the non-S group (difference in intercept 362.79 mL, p=0.03; table 4⇓). In addition, a tendency towards higher values of FEV1 % pred was observed in the S group. Differences in the two lung function parameters between the S and post-S groups were not significant, but the S group tended to have a higher increase of VC over time after the age of 12 yrs than the post-S group. At the age of 12 yrs, patients from the post-S group tended to have better lung function parameters than patients from the non-S group, but with increasing age the changes in lung function over time were similar.
Differences in the IgG levels between the S and non-S groups in patients <12 yrs old, in favour of the screened patients, persisted after the age of 12 yrs but were no longer significant. The IgG levels throughout this observation period were similar in the S and post-S groups. Although the IgG values appeared to be lower in the post-S group, differences in IgG levels between post-S and non-S groups were not significant in this analysis.
Influence of centre treatment
Analysis of outcome until 12 yrs of age
Observed values of weight were higher in the C group than in the non-C group (fig. 4⇑). Similar results were observed for sd height. Nevertheless, in the LME analysis, the increase of both nutritional parameters over time was similar in both groups, and differences in intercept and slope from sd scores for height and weight were statistically not significant. Until the age of 12 yrs, FEV1 declined in both groups over time (fig. 5⇑). This decline was slightly, but not significantly, less in the C group. A greater increase in VC over time was found in the C group, but differences in slopes of VC between C and non-C groups were statistically not significant.
IgG levels at the age of 6 yrs were found to be lower in the C group than in the non-C group. Until the age of 12 yrs, both groups showed an increase in IgG over time. This increase was slightly, but not significantly, higher in the C group.
Analysis of outcome above 12 yrs of age
Although centre-treated patients seemed heavier and taller than the noncentre-treated patients at the age of 12 yrs (fig. 4⇑), differences in height and weight did not reach significance. The increase of both height and weight over time was similar in both groups.
After 12 years of age, the decline in FEV1 % pred and the increase in VC over time tended to be in favour of the C group. Figure 5⇑ shows that after 12 yrs of age FEV1 continues to decline in the non-C group, whereas average observed values of FEV1 stay rather constant over time in the in the C group In addition, the observed values of VC were much higher in the C group than in the non-C group. Nevertheless, in the LME analysis, the differences in intercept and slope of FEV1 and VC were not significant for centre-treated and noncentre-treated patients.
Due to the slightly higher increase in IgG over time before the age of 12 yrs, the C group showed, in patients >12 yrs of age, a tendency towards higher IgG levels than the non-C group. However, the differences in both intercept and slope of IgG between the two groups were statistically not significant.
Influence of genotype on the linear mixed-effects analyses of clinical outcome of patients with cystic fibrosis
Data on genotype were available from 55 patients. Genotype did not significantly affect the influence of screening in the LME analyses of clinical outcome. Also, the influence of centre treatment on clinical outcome was not influenced by differences in genotype between patients, except for the parameter IgG. The LME analysis of outcome above 12 yrs of age, including the variable genotype, showed a difference in IgG between the C and non-C groups, with respect to the intercept of −0.11 mg·dL−1 and the slope of +8.73 mg·dL−1. In contrast, the same LME analysis without the variable genotype showed a difference between the C and non-C groups in intercept from IgG of +99.47 mg·dL−1 and a difference in slope from IgG of −13.94 mg·dL−1. However, this difference was, like the other differences, statistically not significant.
Discussion
Neonatal screening
This study demonstrates the long-term pulmonary benefits of neonatal screening for CF. S patients showed a smaller decline in FEV1 % pred, mainly before 12 yrs of age, together with higher levels of VC than non-S patients. In addition, S patients had significantly lower IgG levels until 12 yrs of age than non-S patients. A better preservation of lung function in the S group was found, despite an overall better clinical condition in the non-S group in patients >12 yrs of age, due to a greater loss of severe patients who died at a young age in the non-S group. Improved lung function associated with neonatal screening for CF has also been reported in a recent study from Australia 6. S patients tended to have a longer survival into adulthood than non-S patients (fig. 1⇑).
Although the S group appeared to be heavier and taller than the non-S group, differences in nutritional status were not significant. A greater catch-up growth after diagnosis, as well as an overall better clinical condition after the age of 12 yrs in the non-S group, may be possible explanations for the similar nutritional status. Early nutritional benefits of neonatal screening have also been found in the randomized, controlled trial from Wisconsin 3. Significantly, beneficial effects of neonatal screening on clinical outcome and survival could not be found in the comparison between S patients and the younger post-S patients. However, this comparison was biased for the following reasons. First, CF patients born later are likely to have a better prognosis than CF patients born earlier 15, 16. In the present study, the post-S group had a significantly better clinical outcome until 12 yrs of age than the non-S group. After 1980, CF patients received better pancreatic enzyme replacement therapy and were no longer treated with a low fat diet. Despite better nutritional management from time of diagnosis, the outcome of the post-S group was not significantly different from that of the S group. Moreover, S patients seemed to have a better preservation of lung function than post-S patients. These results indicate that early intervention through neonatal screening can have long-term pulmonary benefits that are superior to the recognized beneficial effects of a good nutritional status on lung function 17. Secondly, despite a longer follow-up than in a previous report 8, the frequency of CF at birth in this group was relatively low and it is suspected that the post-S group was (still) incomplete. It is speculated that in this group, patients died at a young age prior to diagnosis of CF because of the high age at diagnosis and the proportionally high number of patients with meconium ileus. Incomplete ascertainment is likely to result in biased comparisons of clinical outcome and survival.
Centre treatment
Treatment at a specialized centre improved long-term survival of patients with CF (fig. 2⇑). The RR of dying for the C group compared with the non-C group was 0.18 (95% CI 0.05–0.57, p<0.01). Studies from Denmark 18 and Australia 19 have also shown longer survival in CF patients receiving care at a specialized centre. Over time, treatment was improving, as more and more patients were being treated at specialized centres. This may have resulted in a biased comparison between C and non-C patients. However, the differences in clinical outcome and survival between C and non-C groups were adjusted for the number of patients born during and after the screening programme in each group.
Although not statistically significant, centre-treated patients seemed to be heavier and taller and tended to have a better preservation of lung function in this 17-yr follow-up study than patients without centralized care. Since differences in nutritional status and lung function between treatment groups were in the same direction, they are, therefore, likely to be real and not due to chance. Moreover, differences in nutritional status and lung function did not change significantly after adjustment for genotype.
IgG levels were slightly, but not significantly, higher in the C group than in the non-C group after the age of 12 yrs. Serum immunoglobulin levels often rise after chronic colonization with Pseudomonas auruginosa in CF 20. The percentage of patients with pseudomonas colonization was also slightly, but not significantly, higher among C patients: 57.9% in the C group and 51.6% in the non-C group, p=0.60. The observed differences in IgG above 12 yrs of age may be caused by differences in genotype between patients. Severe genotypes are often associated with a higher colonization rate with P. aeruginosa than mild genotypes 21, which may lead to higher IgG levels. The C group contained more patients with a severe mutation than the non-C group. After adjustment for genotype, similar IgG levels above 12 yrs of age in both groups were found.
The results of the comparison of clinical outcome between C and non-C patients may have been affected by differences in survival. Before the age of 12 yrs, none of the patients from the C group died, in contrast with the non-C group in whom 17 out of 47 patients (36%) died. Overall, clinical condition of non-C patients was better in the analysis in patients >12 yrs old than <12 yrs old, probably due to the absence of patients with a severe clinical expression who died at a young age. The selection bias introduced by survival of the fittest may be an explanation for the nonsignificant differences found in the comparison of long-term clinical outcome between the C and non-C groups.
Although this study was not randomized in design, the allocation of patient groups was nonselective. Complete ascertainment of the S and non-S groups is probable as they were comparable in size and frequency of CF at birth (table 1⇑). Age at diagnosis in the non-S group did not differ significantly from that of patients in the national Dutch Cystic Fibrosis Registry 22. However, the number of patients who received centre care was significantly greater in the S group than in the non-S group (table 2⇑). Centre treatment was associated with an improved survival of patients with CF in the present study (fig. 2⇑). Therefore, the analyses of survival and clinical outcome were adjusted for centre treatment. In addition, comparison of clinical outcome was adjusted for the possible confounding variables sex and meconium ileus. Furthermore, the analyses did not change significantly after adjustment for genotype 23. Therefore, differences in outcome between screened and nonscreened patients cannot be explained by differences in centre treatment, sex, meconium ileus or genotype, and selection bias is highly improbable.
The LME model regression analysis applied in this study to compare clinical outcome between patients, is considered as the statistical technique of choice in clinical trials involving the longitudinal analysis of lung function data 24. The advantage of this methodology is that estimates of average changes can be computed even when individual subjects have had variable length of follow-up and timing of observations. Moreover, covariables can be incorporated in the analysis. Selection bias can be prevented with this type of analysis.
Although differences in prognosis between S and non-S patients and C and non-C patients observed in the present study were in the same direction in favour of the S and C groups, only a few were statistically significant. However, the statistical power of the present study was too low, which may lead to false-negative results. Recently, it has been calculated that studies of this kind need ≥66 patients in each patient group to detect a 10% difference in an outcome measure such as FEV1 25.
The use of the meconium test led to five patients from the S group being undetected. These patients were included in the S group in the present report, in order to perform an intention-to-screen analysis. However, present screening tests for CF possess sensitivities of ≥97% 26, 27.
Conclusions
The results of the present longitudinal study show that neonatal screening for cystic fibrosis lead to pulmonary benefits in adulthood. Moreover, screened patients probably have a longer survival than non-screened patients. The statistically significant differences in this small sized study although few in number, suggest that while selection bias is improbable early diagnosis by screening and early treatment leads to substantial clinical benefit for patients with cystic fibrosis. A striking finding was the tendency towards better lung function parameters in screened patients born earlier, compared to the nonscreened patients born later, despite better nutritional management in the younger cohort. Centralized treatment resulted in an improved long-term survival of patients with cystic fibrosis. The present findings imply that management of all patients with cystic fibrosis should be carried out at, or in close cooperation with, a specialized centre, immediately after a diagnosis is made by neonatal screening.
Acknowledgments
The authors thank G. Slingerland for helpful assistance in the collection of data and G. te Meerman for helpful discussions on the statistical analysis. The help of the Dutch Cystic Fibrosis registration and Statistics Netherlands in Voorburg for finding all patients from the area is gratefully acknowledged. The authors are grateful to all patients, paediatricians, and pulmonologists without whom this study would not have been possible.
- Received September 12, 1999.
- Accepted April 17, 2001.
- © ERS Journals Ltd