To the Editor:
Described as a “template for success” by Hanekom et al. [1] the Beijing genotype of Mycobacterium tuberculosis has been associated with hypervirulence, drug resistance, evasion of the bacille Calmette–Guérin (BCG) vaccine and differential immunoregulation [1]. The genotype is itself diverse; accordingly the fact that specific traits have been associated with Beijing only in certain settings may be explained by the variation in the subtypes that predominate in each.
We sought to explore the role that the Beijing genotype plays in the epidemiology of multidrug-resistant (MDR) tuberculosis (TB) in two very different epidemiological settings: England and Peru. In the UK nearly three quarters of TB disease occurs amongst migrants, which may result from latent M. tuberculosis infections acquired from high-burden countries of origin. Thus the distribution of lineages in the UK reflects a composite of the circulating strains in other nations. Whilst increasing levels of MDR-TB in the UK have triggered understandable concern, only a small fraction of cases of TB are MDR (in 2011 this was 1.0%) [2]. In contrast, the Peruvian TB epidemic is relatively “home-grown” and the prevalence of drug resistance is greater (in 2011 5.3% were MDR-TB) [2]. Our objectives were thus: 1) identify any association between the Beijing genotype and MDR in these two settings, and 2) to examine if trends in Beijing might explain changes in the prevalence of MDR-TB over time.
Four datasets were used during this study: two English studies and two Peruvian studies.
The English datasets were taken from Public Health England's Enhanced Tuberculosis Surveillance System and Mycobacterial Surveillance Network and were referred to as England 1 and England 2. England 1 consisted of 8 859 individuals diagnosed and 24-loci MIRU-VNTR (mycobacterial interspersed repetitive units variable number tandem repeats) strain-typed between 2010 and 2011. Of these, 4 661 had at least first-line drug sensitivity data, a minimal 12-loci predefined profile [3] and were pulmonary. England 2 consisted of 189 MDR-TB cases from 2004 to 2007, which had undergone 24-loci typing. 117 were pulmonary and has at least the minimal standard 12-loci profile. 24-loci typing was not available for non-MDR cases during this period. Susceptibility to rifampicin, isoniazid, ethambutol, pyrazinamide, streptomycin, amikacin, kanamycin, capreomycin, azithromycin, clarithromycin, ciprofloxacin, ofloxacin, moxifloxacin, prothionamide, ethionamide, cycloserine, 4-aminosalicylic acid, linezolid, clofazamine and rifabutin had been tested using the proportion or resistance ratio methods.
The Peruvian datasets were referred to as Peru 1 and Peru 2. Peru 1 consisted of 323 sputum-positive cases of pulmonary TB, recruited from clinical facilities in Lima between 2004 and 2006 and extracted from Taype et al. [4]. 229 had been tested for sensitivity to rifampicin, isoniazid, ethambutol and streptomycin, of which 101 (44.1%) were MDR; we note that this figure is much higher than the overall proportion of MDR in this setting [2], suggesting that drug susceptibility testing was prioritised for those at highest risk of MDR, consistent with previous Peruvian testing policies. All cases had a complete standard 12-loci profile. Peru 2 consisted of 230 cases of pulmonary MDR-TB identified in Lima between 1996 and 2004 [5, 6]. Isolates underwent drug susceptibility testing to rifampicin, isoniazid, ethambutol, pyrazinamide, streptomycin, kanamycin, capreomycin, ethionamide, 4-aminosalicylic acid and levofloxacin. One individual lacking a complete standard 12-loci profile was excluded, leaving 229.
MIRU-VNTRplus was utilised to designate lineages for all four datasets to ensure comparability [7]. Substantial differences in the lineage and sublineage profiles of English (England 1, 4661 individuals) and Peruvian (Peru 1, 229 individuals) pulmonary TB cases were apparent (Chi-squared test of independence p-value <0.005). The vast majority of Peruvian cases were of the Euro-American lineage (212 (92.6%) out of 229), with Latin American-Mediterranean (LAM) the majority sublineage (102 (48.1%) out of 212). The remainder of cases were nearly all Beijing (15 (6.6%) out of 229). In contrast, English cases comprised of a variety of lineages consistent with their multinational origins: Euro-American 2412 (51.8%) out of 4661, Delhi-CAS (central-Asia) 1335 (28.6%) out of 4661, East African Indian (EAI) 429 (9.2%) out of 4661, Beijing 392 (8.3%) out of 4661. The Euro-American sublineages observed were also diverse.
Restricting our analysis to MDR cases revealed very different profiles between Peru and England (p-value <0.005). Peruvian cases (Peru 2, 229 individuals) were dominated by the Euro-American lineage (215 (93.9%) out of 229), of which 121 (56.3%) out of 215 was the LAM sublineage. Regression modelling of Peru 1 to explore the relationship between MDR and M. tuberculosis lineage showed no association (Beijing crude OR (95% CI) 1.10 (0.38–3.14), Wald p-value 0.86; Euro-American baseline). In stark contrast to the overall distribution of lineages, among English MDR cases (England 1) the proportion Beijing increased to 36 (47.4%) out of 76. Regression modelling revealed a strong positive association between Beijing and MDR among English pulmonary cases, (crude OR 11.51 (6.65–19.95), p<0.005; adjusted OR 9.53 (4.87–18.65), p<0.005; Euro-American baseline and England 1 used for both models; models excluded non-M. tuberculosis cases and those with missing data). Age, sex, a composite social risk factor variable (homelessness, imprisonment, drug/alcohol abuse), country of birth, ethnic group and previous diagnosis were adjusted for in the latter; this could not be done for Peru 1 as these variables were not available.
We also tested whether Beijing was enriched among strains with resistance to a greater number of drugs (fig. 1a). Pre-MDR-TB (pre-MDR) was defined as resistance to either, but not both, isoniazid and rifampicin. Pre-extensively drug resistant (pre-XDR) cases were MDR plus resistant to either a fluoroquinolone or a second-line injectable. Beijing was enriched from pre-MDR to XDR cases in the English data (Pearson's correlation coefficient (ptrend) Beijing versus non-Beijing p <0.005; England 1). This could be specifically observed among non-UK born, but not UK born, cases (data not shown). No trend was seen in the Peruvian data (ptrend = 0.76; Peru 1 and Peru 2). Further analysis to determine whether the association has changed over time indicated a recent increase in the proportion of cases caused by Beijing (ptrend<0.005; England 1 and England 2) (fig. 1b).
a) The Beijing genotype is enriched with increasing drug resistance in English (England 1) (ptrend<0.005 Bejing versus non-Beijing), but not Peruvian (Peru 1 and Peru 2) pulmonary tuberculosis (TB) cases. b) Lineages of the English pulmonary multidrug resistant (MDR) TB cases (England 1 and England 2), presented by year of diagnosis. The proportion of cases associated with the Beijing genotype may be increasing with time in the latter years. CAS: central Asian; M.: Mycobacterium; EAI: East-African Indian; XDR: extensively drug resistant. #: 24-loci MIRU-VNTR typing was not undertaken between 2008 and 2009.
Our analysis suggests that the Beijing genotype plays a different role in the current epidemiology of MDR-TB in England and Peru. In England, the observed association likely reflects the rising importation of resistant strains, rather than the development of resistance within the UK. In Peru there is likely to be a similarly high prevalence of MDR among all genotypes, potentially explaining why we did not observe an association. Interestingly, a recent paper by Iwamoto et al. [8] suggests that the overall prevalence of Beijing has increased in Peru over the last decade, but the proportion of Beijing strains that are MDR has decreased.
We utilised available datasets from studies undertaken at different times with variable designs, thus some interpretative caution is required. For example, Peru 2 recruited from households with more than one pulmonary MDR case, and thus could have oversampled clustered cases. However, the similarity of the lineage profiles of Peru 2 with Peru 1, which was not collected through such means, suggests this effect is limited. English cases were tested for resistance to a wider range of drugs than the Peruvian cases, thus the “no known resistance” category is not strictly comparable between the two datasets. Nevertheless, only 1.3% of the England 1 dataset was resistant to drugs not tested within the Peru 1 dataset.
Molecular epidemiological surveillance and addressing and preventing resistance have been highlighted as vital for TB elimination by the European Respiratory Society [9, 10]. We demonstrate within this study how the tools of the former could be vital to aid the latter: within Peru the MDR epidemic could worsen should fitter drug resistant Beijing strains evolve or enter the country and within the UK where the lineage profile is highly dynamic.
Acknowledgments
We acknowledge our utilisation of the published data of Taype et al. [4]. Our study made use of the MIRU-VNTRplus database website (www.miru-vntrplus.org) developed by D. Harmsen, S. Niemann, P. Supply and T. Weniger.
Footnotes
Support statement: H.R. Stagg was supported by a European Respiratory Society Fellowship STRTF 24-2011. The strain typing of the English datasets was funded by Public Health England's Strain Typing Project. I. Abubakar was funded through an NIHR Senior Research Fellowship NIHR-2011-04-001.
Conflict of interest: Disclosures can be found alongside the online version of this article at www.erj.ersjournals.com
- Received July 3, 2013.
- Accepted August 31, 2013.
- ©ERS 2014