Tuberculosis elimination in the Netherlands.

This study assessed progress towards tuberculosis (TB) elimination in the Netherlands by using DNA finger-printing. Mycobacterium tuberculosis strains were defined as new if the IS6110 restriction fragment length polymorphism pattern had not been observed in any other patient during the previous 2 years. Other cases were defined as clustered and attributed to recent transmission. In the period 1995-2002, the incidence of TB with new strains was stable among non-Dutch residents and declined among the Dutch. However, the decline among the Dutch was restricted to those >65 years of age. Moreover, the average number of secondary cases per new strain did not change significantly over time. We conclude that the decline of TB in the Netherlands over the past decade was mainly the result of a cohort effect: older birth cohorts with high infection prevalence were replaced by those with lower infection prevalence. Under current epidemiologic conditions and control efforts, TB may not be eliminated.

**********

Reported rates of tuberculosis (TB) in the Netherlands in 2003 were 3.5 per 100,000 among Dutch residents and 125 per 100,000 among the non-Dutch. The non-Dutch are formally defined as those without a Dutch passport, but in practice they include mostly foreign-born persons. During the past 10 years, TB reports in the Netherlands declined among the Dutch (from 693 in 1993 to 531 in 2002) and remained approximately stable among the non-Dutch (at an average of 892 per year). To what extent these patterns were attributable to changes in TB transmission and to what extent to changes in the introduction of new strains from abroad or from reactivation of latent infection are unclear.

According to mathematical models, eliminating TB as a public health problem, defined as a prevalence of Mycobacterium tuberculosis infection of <1%, has been predicted to occur among the Dutch by 2030 if no further changes in transmission occur (1,2). However, if elimination is defined as an incidence of smear-positive disease of <1 case per million population per year (1), TB is not expected to be eliminated by that date because of transmission from foreign-source cases (2). The proportion of all Dutch TB cases in the period 1993-1998 attributable to recent transmission from a non-Dutch source case has been estimated at 12% to 20% with DNA-fingerprinting data (3). With time, an increasing proportion of TB cases among the Dutch were expected to be attributable to recent transmission from non-Dutch source cases (2). However, direct observations with fingerprinting results have not yet been used to evaluate these model predictions.

Immigration patterns in the Netherlands have varied in the past decade: large numbers of persons from countries with high TB endemicity have sought asylum in the early 1990s. In recent years, these numbers became smaller after stricter immigration laws were passed. Shifts in countries of origin of immigrants have also occurred, and some of these countries had much higher TB rates than others. Therefore, the introduction of new strains from abroad may be expected to have varied over time. Control measures, in contrast, have shown little change over the study period. TB screening is obligatory at entry and is offered every 6 months for 2 years on a voluntary basis. No routine screening for and treatment of latent infection exist for immigrants.

This study attempted to determine, by DNA-fingerprinting of M. tuberculosis isolates, to what extent TB trends from 1995 to 2002 were determined by changes in the introduction of new strains and by changes in ongoing transmission. We also describe the trend of TB transmission from non-Dutch source patients to the Dutch population. The combined evidence is used to assess the prospects for eliminating TB in the Netherlands.

Methods

Patient and treatment data since 1993 were available in the Netherlands Tuberculosis Register, an anonymous case register maintained by the KNCV Tuberculosis Foundation. Reporting to the register is voluntary, but cross-matching with mandatory reports to the ministry of health on all patients who have started TB treatment suggests >99% completeness. The register includes data on demographic characteristics, clinical features, risk groups, treatment given, and treatment outcome. From January 1993 to December 2002, a total of 15,331 TB patients were registered, including those without bacteriologic confirmation.

Over the same period, 10,356 first M. tuberculosis isolates of TB patients were subjected to standard IS6110 restriction fragment length polymorphism (RFLP) analysis (4). Subtyping with the polymorphic GC-rich sequence probe was carried out for strains with <5 IS6110 copies. IS6110 RFLP patterns were analyzed by using the Bionumerics software, version 3.5 for Windows (Applied Maths, Sint-Martens-Latem, Belgium).

Information from the 2 databases was combined; sex, date of birth, postal area code, and year of diagnosis were used as identifiers. A perfect match was obtained for 7,529 (73%) isolates and a near-perfect match for 981 (9%). Both groups were included, yielding a total study size of 8,510 (82%) culture-positive patients. Mismatches may be due to administrative errors in a database, unreliable date of birth (e.g., for some immigrant groups), or postal area code (e.g., homeless), and the exclusion of persons with identical identifiers.

TB can occur soon after primary infection or reinfection (recent transmission) or as the result of endogenous reactivation of latent infection (5). The cut-off point for separating recent from remote transmission is arbitrary: some researchers used 5 years (5-7), others 2 years (8), and others 1 year (9). We estimated the percentage of cases with identical RFLP patterns occurring within a given period after each culture-positive case with Kaplan-Meier survival analysis, as suggested by Jasmer et al. (9). The Kaplan-Meier estimate of the probability that a patient was followed by another with an identical fingerprint was 46.2% for the total study period and 33.7% for a 2-year period. Thus, of all cases followed by a patient with an identical fingerprint within 10 years, 73% were followed within 2 years. Using this information, we defined strains as new if the RFLP pattern had not been observed in another patient during the previous 2 years. Other strains were attributed to ongoing transmission. During the first 2 study years (1993-1994), judging whether strains were new was not possible; therefore, data from these 2 years were used to define new strains from 1995 onwards but were otherwise excluded from the analysis.

The observation period in which secondary cases can be observed is longer for strains introduced earlier in the study period than for those introduced later. To obtain an unbiased estimate of the trend of the number of secondary cases generated by source cases, secondary cases arising >2 years after a new strain was introduced were excluded. Thus, all patients were assigned to 1 of the following 3 mutually exclusive categories: case with a new strain, secondary case within 2 years of the introduction of a new strain, and secondary case >2 years after the introduction of a new strain. To assess the trend of transmission between Dutch and non-Dutch persons, secondary cases were attributed to a source case-patient, defined as the patient from whom the new strain was first isolated (10).

Population data by year, age group, sex, and (Dutch/non-Dutch) nationality were obtained from Statistics Netherlands (available from http://statline. cbs.nl/StatWeb) and used as denominators for incidence rates. Relative risks of TB by year of diagnosis, age, sex, and Dutch or non-Dutch nationality were determined separately for new strains and secondary cases with Poisson regression. Risk factors for the average number of secondary cases per new strain were also identified with Poisson regression (3).

Results

Of the 8,510 TB patients with known RFLP results in the period 1993-2002, 1,580 were found in 1993 to 1994, and 6,930 in 1995 to 2002. Of the latter, 4,594 (66%) had new strains, 1,198 (17%) had secondary cases within 2 years of the introduction of a new strain, and 1,138 (16%) had secondary cases >2 years after a new strain was introduced.

The incidence of TB with new strains was on average 52 per 100,000 among the non-Dutch and 1.4 per 100,000 among the Dutch. The incidence declined over the study period among the Dutch (rate ratio per year 0.96, 95% confidence interval [CI] 0.94-0.98) and was stable among the non-Dutch (rate ratio per year 1.02, 95% CI 1.00-1.03, p = 0.06) (Figure 1). The incidence of all cases attributed to recent transmission, regardless of the duration of the cluster, was 23 per 100,000 among the non-Dutch and 0.9 per 100,000 among the Dutch. The incidence declined among the Dutch (rate ratio per year 0.97, 95% CI 0.95-1.00, p = 0.03) but not among the non-Dutch (rate ratio per year 0.99, 95% CI 0.97-1.02) (Figure 1).

[FIGURE 1 OMITTED]

Reduction of TB incidence with new strains among the Dutch was restricted to those [greater than or equal to] 65 years of age (Figure 2). Incidence was stable at 0.85/100,000 in the age group <65 years (rate ratio per year 1.0, 95% CI 0.97-1.03), declined from 3.5 to 2.2 in those 65-74 years of age (rate ratio per year 0.91, 95% CI 0.86-0.95), and declined from 9.4 to 4.8 in those [greater than or equal to] 75 years of age (rate ratio per year 0.92, 95% CI 0.87-0.95). The incidence rate in the [greater than or equal to] 75-year age group declined more rapidly than the number of cases in that age group, since the population in that age group increased from 847,000 in 1995 to 979,000 in 2002.

Of the 4,594 patients with new strains in the period 1995-2002, a total of 3,459 were found in the period 1995-2000 and could be followed up for 2 years. Of the 1,318 Dutch patients with new strains, 182 (14%) generated secondary cases at an average of 1.7 cases per cluster (1.2 Dutch and 0.5 non-Dutch) (Table 1). The average number of secondary cases generated was 0.23 per new strain and declined steeply with the age of the source case-patient (rate ratio per age group 0.74, 95% CI 0.70-0.78) (Table 1). The average number of secondary cases generated did not depend on the sex of the source patient (p > 0.5). The average number of secondary case-patients per new strain did not differ significantly over time (rate ratio per year 0.96, 95% CI 0.89-1.02).

Of the 2,141 non-Dutch patients with new strains, 283 (13%) generated secondary cases at an average of 1.9 cases per cluster (0.5 Dutch and 1.4 non-Dutch) (Table 2). The average number of secondary cases generated was 0.25 overall, declined with age of the source patient (rate ratio per age group 0.86, 95% CI 0.80-0.92), and was lower for female than male source patients (rate ratio 0.68, 95% CI 0.57-0.81) (Table 2). The average number of secondary cases per new strain over time did not change (rate ratio per year 0.97, 95% CI 0.93-1.02).

In clusters starting in the period 1995-2000, an increasing proportion of Dutch secondary TB cases was attributable to a non-Dutch source case as time progressed ([chi square] trend 4.49, p = 0.03) (Table 3). This trend was observed not only among those cases arising within 2 years of the start of the cluster (Table 3) but also among all Dutch secondary cases, regardless of cluster duration ([[chi square].sub.trend] 42, p < 0.001, data not shown). The proportion of Dutch secondary cases attributable to a non-Dutch source case declined steeply with age, both among all Dutch secondary case-patients ([[chi square].sub.trend] 41, p < 0.001) and among those arising within 2 years of the start of the cluster ([[chi square].sub.trend] 27, p < 0.001) (Table 3). The proportion of cases attributed to a non-Dutch source patient was not associated with sex of the Dutch secondary case-patient (Table 3).

Discussion

This study suggests that the declining TB incidence among the Dutch in the Netherlands during the past decade has been achieved under stable control conditions. Among the Dutch, the incidence of TB attributable to new strains declined, particularly among the elderly. The incidence of TB cases due to recent transmission declined as well, a result of fewer new strains being introduced. The average number of secondary cases per new strain did not change significantly. An overall reduction in incidence of clustered cases among the U.S.-born population was also observed in San Francisco (9) and New York (11) and was attributed to improved control. In the Netherlands, we do not attribute the decline to improved control but to a cohort effect. Research may determine to what extent the reported declines in San Francisco and New York could be explained by a reduction in number of secondary cases per newly introduced strain and whether a cohort effect played a role in those settings.

In industrialized countries the annual risk for M. tuberculosis infection has declined steeply over the past century (12); as a result, compared to younger persons, older persons were exposed to much higher risks for infection in their youth (13). Thus, the prevalence of infection increases sharply with age. The risk for TB due to reactivation of latent infection therefore increases with age as well. Within the older age groups, this risk is now declining with each calendar year as earlier birth cohorts leave and more recent birth cohorts with lower infection prevalence enter the age group.

The incidence of culture-positive TB with new strains among Dutch persons <65 years of age was stable in the past decade, at 0.85 per 100,000 population in our matched dataset, and thus [approximately equal to] 1/100,000 or 10 per million if failure to match is taken into account. Of these new case-patients with a known sputum smear result, 63% had smear-positive TB. If elimination of TB as a public health problem is defined as achieving an incidence of new smear-positive TB cases of <1 per million (1), elimination is unlikely to be achieved under current epidemiologic conditions and control efforts.

Our study confirms previous predictions from a mathematical model about the increasing importance of transmission from immigrants to the Dutch population (2). The number of cases observed among the Dutch is best explained by immigrant scenarios 1 and 2 in the modeling study (2), which assume that a Dutch TB patient is 8 times more likely than a non-Dutch TB patient to infect a Dutch person.

Over time, the decline of TB incidence among elderly Dutch will become less important as the birth cohorts with a high prevalence of infection are replaced with cohorts with much lower infection rates. Contact with highly TB-endemic countries through immigrants and international travel, on the other hand, is becoming increasingly important as a determinant of TB trends in the Netherlands. This finding was shown in this study by the increasing proportion over time of Dutch patients with secondary cases attributed to a non-Dutch source patient. This finding suggests the need for further reorientation of the focus of TB control within the Netherlands towards immigrants and their contacts and reemphasizes the importance of global TB control for achieving TB elimination in countries with low incidence of this disease (14).

The separation of TB patients into those with new strains, attributed to reactivation or acquisition abroad, and secondary cases attributed to recent transmission is likely to be imperfect for the following reasons. Some strains identified as new may have represented ongoing transmission in the presence of strain evolution. Some strains attributed to ongoing transmission may represent remote transmission, particularly among the elderly (15). In this national database, epidemiologic confirmation of linkage between patients was far from complete (16). However, in a recent, more detailed study in Amsterdam, most clustered patients were found to have epidemiologic links (17). Missing data as a result of incomplete matching may have contributed to a slight underestimate of the observed clustering percentage and of the number of secondary cases per source case (18,19). However, since the matching percentage was not associated with calender year (data not shown), this underestimate should not affect the trend estimates. Some source cases may have been misclassified, in particular in large clusters. These sources of misclassification are expected to reduce the observed difference between cases with new strains and those attributed to recent transmission but do not invalidate the main conclusion that TB incidence among the Dutch was reduced mainly because of fewer reactivation cases among persons >65 years of age.

We conclude that the decline of TB in the Netherlands during the past decade was mainly the result of a cohort effect: older birth cohorts with high infection prevalence were replaced by those with lower infection prevalence. Contact through immigrants and international travel with countries with high TB incidence increasingly determines TB trends in the Netherlands and will prevent achieving TB elimination under current conditions. Global TB control is required to achieve TB elimination in countries with a low incidence of this disease.

Table 1. Tuberculosis cases with new strains among the Dutch,
1995-2000, and their secondary cases within 2 years

                                                      Other cases in
                                                      these clusters
                                                      within 2 years
                                                       of first case

                                      Cases with
                                     new strains
                      Cases with    being first of
                      new strains   cluster, n (%)   Dutch   Non-Dutch

Year
  1995                    229          34 (15)        52         17
  1996                    255          36 (14)        40         19
  1997                    227          29 (13)        32          8
  1998                    193          31 (16)        34         16
  1999                    220          29 (13)        31         10
  2000                    194          23 (12)        27         17
Age group
  <25                     102          32 (31)        43         16
  25-34                   154          41 (27)        43         30
  35-44                   140          28 (20)        37         10
  45-54                   115          11 (10)        14          6
  55-64                   150          19 (13)        22          9
  65-74                   219          22 (10)        28          8
  [greater than
    or equal to] 75       438           29 (7)        29          8
Sex
  Male                    740          103 (14)       121        44
  Female                  578           79 (14)       95         43
Total                    1,318         182 (14)       216        87

Table 2. Tuberculosis cases with new strains among the Dutch,
1995-2000, and their secondary cases within 2 years

                                                      Other cases in
                                                      these clusters
                                                      within 2 years
                                                       of first case

                                      Cases with
                                     new strains
                      Cases with    being first of
                      new strains   cluster, n (%)   Dutch   Non-Dutch

Year
  1995                    346          39 (11)        19        68
  1996                    334          46 (14)        21        70
  1997                    361          46 (13)        19        52
  1998                    346          43 (12)        20        72
  1999                    383          63 (16)        29        77
  2000                    371          46 (12)        19        61
Age group
  <25                     580          101 (17)       51        151
  25-34                   783          98 (13)        24        124
  35-44                   370          42 (11)        33        78
  45-54                   171          18 (11)         5        22
  55-64                   116          13 (11)         9        15
  65-74                   88            6 (7)          2         4
  [greater than
    or equal to] 75       33            5 (15)         3         6
Sex
  Male                   1,226         171 (14)       82        264
  Female                  915          112 (12)       45        136
Total                    2,141         283 (13)       127       400

Table 3. Dutch tuberculosis cases attributed to recent transmission
and diagnosed within 2 years of the start of clusters

                            First case of cluster

                      Dutch          Non-Dutch
                              (% non-Dutch first case)

Year
  1995                 16              3 (16)
  1996                 40             19 (32)
  1997                 42             22 (34)
  1998                 32             18 (36)
  1999                 34             27 (44)
  2000                 27             23 (46)
  2001                 19              9 (32)
  2002                  6              6 (50)
Age group
  <25                  29             43 (60)
  25-34                46             32 (41)
  35-44                32             16 (33)
  45-54                27             12 (31)
  55-64                24              9 (27)
  65-74                34             10 (23)
  [greater than
    or equal to] 75    24              5 (17)
Sex
  Male                 124            80 (39)
  Female               92             47 (34)
Total                  216            127 (37)

References

(1.) Styblo K. The elimination of tuberculosis in the Netherlands. Bull Int Union Tuberc Lung Dis. 1990;65:49-55.

(2.) Wolleswinkel-van den Bosch JH, Nagelkerke N J, Broekmans JF, Borgdorff MW. The impact of immigration on the elimination of tuberculosis in The Netherlands: a model-based approach. Int J Tuberc Lung Dis. 2002;6:130-6.

(3.) Borgdorff MW, Nagelkerke N J, De Haas PE, Van Soolingen D. Transmission of tuberculosis depending on the age and sex of source cases. Am J Epidemiol. 2001;154:934-43.

(4.) Van Embden JDA, Cave MD, Crawford JT, Dale JW, Eisenach KD, Gicquel B, et al. Strain identification of Mycobacterium tuberculosis by DNA fingerprinting: recommendations for a standardized methodology. J Clin Microbiol. 1993;31:406-9.

(5.) Vynnycke E, Fine PEM. The natural history of tuberculosis: the implications of age-dependent risks of disease and the role of reinfection. Epidemiol Infect. 1997;119:183-201.

(6.) Rieder HL. Epidemiologic basis of tuberculosis control. Paris: International Union against Tuberculosis and Lung Disease; 1999. p. 66.

(7.) Sutherland I. Recent studies in the epidemiology of tuberculosis, based on the risk of being infected with tubercle bacilli. Adv Tuberc Res. 1976;19:1-63.

(8.) Cronin WA, Golub JE, Lathan M J, Mukasa LN, Hooper N, Razeq JH, et al. Molecular epidemiology of tuberculosis in a low- to moderate-incidence state: are contact investigations enough? Emerg Infect Dis. 2002;8:1271-9.

(9.) Jasmer RM, Hahn JA, Small PM, Daley CL, Behr MA, Moss AR, et al. A molecular epidemiologic analysis of tuberculosis trends in San Francisco, 1991-1997. Ann Intern Med. 1999;130:971-8.

(10.) Behr MA, Warren SA, Salamon H, Hopewell PC, Ponce de Leon A, Daley CL, et al. Transmission of Mycobacterium tuberculosis from patients smear-negative for acid-fast bacilli. Lancet. 1999;353:444-9.

(11.) Geng E, Kreiswirth B, Driver C, Li J, Burzynski J, DellaLatta P, et al. Changes in the transmission of tuberculosis in New York City from 1990 to 1999. N Engl J Med. 2002;349:1453 8.

(12.) Styblo K, Meijer J, Sutherland I. The transmission of tuberculosis bacilli. Its trend in a human population. Bull Int Union Tuberc. 1969;42:5-104.

(13.) Frost WH. The age selection of mortality from tuberculosis in successive decades. Am J Hyg. 1939;30:91-6.

(14.) Bloom BR. Tuberculosis-the global view. N Engl J Med. 2002;346:1434-5.

(15.) Lillebaek T, Dirksen A, Vynnycky E, Baess I, Thomsen VO, Andersen AB. Stability of DNA patterns and evidence of Mycobacterium tuberculosis reactivation occurring decades after the initial infection. J Infect Dis. 2003;188:1032-9.

(16.) Lambregts-van Weezenbeek CS, Sebek MM, van Gerven PJ, van Gerven P J, de Vries G, Verver S, et al. Tuberculosis contact investigation and DNA fingerprint surveillance in The Netherlands: 6 years' experience with nation-wide cluster feedback and cluster monitoring. Int J Tuberc Lung Dis. 2003;7(Suppl 3): 463-70.

(17.) van Deutekom H, Hoijng SP, de Haas PE, Langendam MW, Horsman A, van Soolingen D, et al. Clustered tuberculosis cases: do they represent recent transmission and can they be detected earlier? Am J Respir Crit Care Med. 2004; 169:806-10.

(18.) Glynn JR, Vynnycky E, Fine PE. Influence of sampling on estimates of clustering and recent transmission of Mycobacterium tuberculosis derived from DNA fingerprinting techniques. Am J Epidemiol. 1999;149:366-71.

(19.) Van Soolingen D, Borgdorff MW, De Haas PE, Sebek MM, Veen J, Dessens M, et al. Molecular epidemiology of tuberculosis in The Netherlands: a nationwide study from 1993 through 1997. J Infect Dis. 1999;180:726-36.

Martien W. Borgdorff, * ([dagger]) Marieke J van der Werf, * Petra E.W. de Haas, ([double dagger]) Kristin Kremer, ([double dagger]) and Dick van Soolingen ([double dagger])

* KNCV Tuberculosis Foundation, The Hague, the Netherlands; ([dagger]) University of Amsterdam, the Netherlands; and National Institute of Public Health and the Environment, Bilthoven, the Netherlands

Dr. Borgdorff is professor of International Health at the University of Amsterdam and executive director of the KNCV Tuberculosis Foundation. His areas of interest include the molecular epidemiology of TB and other communicable diseases.

Address for correspondence: Martien W. Borgdorff, KNCV Tuberculosis Foundation, PO Box 146, 2501 CC The Hague, the Netherlands; fax: [phone]; email: [e-mail]