miércoles, 20 de febrero de 2019

Cross-Border Movement of Highly Drug-Resistant Mycobacterium tuberculosis from Papua New Guinea to Australia through Torres Strait Protected Zone, 2010–2015 - Volume 25, Number 3—March 2019 - Emerging Infectious Diseases journal - CDC

Cross-Border Movement of Highly Drug-Resistant <em>Mycobacterium tuberculosis</em> from Papua New Guinea to Australia through Torres Strait Protected Zone, 2010–2015 - Volume 25, Number 3—March 2019 - Emerging Infectious Diseases journal - CDC

Volume 25, Number 3—March 2019


Cross-Border Movement of Highly Drug-Resistant Mycobacterium tuberculosis from Papua New Guinea to Australia through Torres Strait Protected Zone, 2010–2015

Arnold Bainomugisa, Sushil Pandey, Ellen Donnan, Graham Simpson, J’Belle Foster, Evelyn Lavu, Stenard Hiasihri, Emma S. McBryde, Rendi Moke, Steven Vincent, Vitali Sintchenko, Ben J. Marais, Lachlan J.M. Coin, and Christopher CoulterComments to Author 
Author affiliations: Pathology Queensland, Brisbane, Queensland, Australia (A. Bainomugisa, S. Pandey, C. Coulter)The University of Queensland, Brisbane (A. Bainomugisa, L.J.M. Coin)Queensland Department of Health, Brisbane (E. Donnan, C. Coulter)Cairns Tuberculosis Unit, Cairns, Queensland, Australia (G. Simpson, S. Vincent)James Cook University, Townsville, Queensland, Australia (J. Foster, E.S. McBryde)Central Public Health Laboratory, Port Moresby, Papua New Guinea (E. Lavu)Western Province Health Office, Daru, Papua New Guinea (S. Hiasihri)Port Moresby General Hospital, Port Moresby (R. Moke)Westmead Hospital, Sydney, New South Wales, Australia (V. Sintchenko)The University of Sydney, Sydney (V. Sintchenko, B.J. Marais)


In this retrospective study, we used whole-genome sequencing (WGS) to delineate transmission dynamics, characterize drug-resistance markers, and identify risk factors of transmission among Papua New Guinea residents of the Torres Strait Protected Zone (TSPZ) who had tuberculosis diagnoses during 2010–2015. Of 117 isolates collected, we could acquire WGS data for 100; 79 were Beijing sublineage, which was associated with active transmission (odds ratio 6.190, 95% CI 2.221–18.077). Strains were distributed widely throughout the TSPZ. Clustering occurred more often within than between villages (p = 0.0013). Including 4 multidrug-resistant tuberculosis isolates from Australia citizens epidemiologically linked to the TSPZ into the transmission network analysis revealed 2 probable cross-border transmission events. All multidrug-resistant isolates (33/104) belonged to Beijing sublineage and had high-level isoniazid and ethionamide co-resistance; 2 isolates were extensively drug resistant. Including WGS in regional surveillance could improve tuberculosis transmission tracking and control strategies within the TSPZ.
Tuberculosis (TB) is the leading infectious cause of death globally (1). To reduce the burden of TB, many countries committed to achieving a 90% reduction in TB incidence by 2035 as part of the End TB Strategy (2). Australia has already achieved preelimination targets (<10 cases/1 million population) in the nonindigenous, Australia-born population (3). However, achieving TB elimination in Australia remains a daunting challenge, given high population mobility, increased importation of TB cases from high-incidence settings, and cross-border spread from neighboring countries, such as Papua New Guinea (4). In 2013, ≈90% of the TB cases reported in Australia were in persons born overseas (3). The state of Queensland has one of the lowest TB notification rates in Australia (4.0 cases/100,000 population) (5); the Western Province of neighboring country Papua New Guinea has a significantly higher incidence, estimated at 2,901 cases/100,000 population at the provincial capital, Daru (6,7).
Thumbnail of Communities within and boundaries of Torres Strait Protected Zone. Green arrow indicates Saibai Island, Australia; red arrow Mabaduan village, Papua New Guinea; and blue arrow Daru Island, Western Province Provincial Capital, Papua New Guinea. Inset depicts location of Torres Strait Protected Zone between Australia and Papua New Guinea. Operational details on how the 1978 Torres Strait Treaty functions can be accessed online (https://dfat.gov.au/geo/torres-strait/Documents/torres-st
Figure 1. Communities within and boundaries of Torres Strait Protected Zone. Green arrow indicates Saibai Island, Australia; red arrow Mabaduan village, Papua New Guinea; and blue arrow Daru Island, Western Province Provincial Capital,...
Daru General Hospital (Port Moresby, Papua New Guinea) is the main health center that offers TB health services to residents of Western Province, including those in the surrounding areas, such as the Torres Strait Protected Zone (TSPZ). The TSPZ is an area where free bidirectional cross-border movement (without passports or visas) is permitted for purposes of traditional customs and economic activities (Figure 1); the zone was created with the signing of the 1978 Torres Strait Treaty by Papua New Guinea and Australia. This area contains a number of Papua New Guinea villages and 14 Australia island communities (estimated population 1,526) that are part of Queensland (9). Cross-border movement of populations within the TSPZ provides a potential route of entry of Mycobacterium tuberculosis into northern Queensland and its spread elsewhere in Australia.
In 2012, an Australia resident of the TSPZ was reported to have pulmonary TB resistant to streptomycin only. After initial improvement with standard therapy, the patient deteriorated clinically. Rifampin resistance was detected by Xpert MTB/RIF assay (Cepheid, http://www.cepheid.com), and multidrug-resistant (MDR) TB was confirmed after positive culture and drug susceptibility testing. The mycobacterial interspersed repetitive unit 24 (MIRU-24) profile of the new isolate differed from the initial isolate and was characteristic of TB strains in Western Province, indicating a possible reinfection with an MDR TB isolate (10). This incident was the first notified case of MDR TB in a citizen of Australia in the Torres Strait Islands. During outbreak epidemiologic investigations, 3 additional MDR TB diagnoses were made in citizens of Australia who were current or previous residents of the TSPZ. These 4 cases were epidemiologically linked through a network of close contacts. At the time of diagnosis, 2 of the patients were living in different Queensland cities distant from the TSPZ.
We sought to use whole-genome sequencing (WGS) to determine the genomic relationship between these isolates and the MDR TB isolates from Papua New Guinea as proof of principle for transmission of MDR TB through the TSPZ. The discriminatory power of WGS enables the delineation of TB transmission with a higher resolution than conventional genotyping (11,12). The aim of this study was to use WGS and epidemiologic data to determine strain diversity within the TSPZ, characterize geno-resistance markers, and identify potential risk factors for transmission.


Strain Selection
All patients receiving TB diagnoses in Queensland (including Papua New Guinea citizens receiving TB diagnoses in Queensland clinics in the TSPZ) must be reported to the Queensland Department of Health in Brisbane. In this study, we included isolates from Papua New Guinea citizens residing in the TSPZ who received TB diagnoses in the TSPZ during January 1, 2010–December 31, 2015, which we refer to as cross-border TB isolates. We extracted demographic and clinical data from the Queensland Notifiable Conditions System on April 6, 2016. We included 1 patient given a TB diagnosis previously in Papua New Guinea who subsequently sought treatment in the TSPZ for extensively drug-resistant (XDR) TB; we also included the 4 previously mentioned MDR TB cases in residents of Australia because these infections were linked to the TSPZ. During the study period, a small number of drug-susceptible TB cases (n = 14) were reported among Australia citizens residing in the TSPZ, but we did not include these cases in this study.
Drug Susceptibility Testing and Genotyping
We performed M. tuberculosis culture, species identification, genotyping, and drug susceptibility testing at Queensland Mycobacterium Reference Laboratory (Pathology Queensland, Brisbane). MIRU-24 or MIRU-15 results obtained as previously described (13) were available for most isolates. We determined phenotypic susceptibility to first- and second-line drugs as previously described (7). On a case-by-case basis, we used SENSITITRE Microbroth Dilution Method (Trek Diagnostic Systems, https://www.thermofisher.com) to resolve differences between resistance mutations detected and phenotypic drug susceptibility test results.
WGS Analysis
We retrieved M. tuberculosis isolates from −80°C storage, cultured them on Lowenstein–Jensen medium, and extracted isolate DNA using an organic enzymatic method (7). We performed WGS at the Australia Genome Research Facility (Brisbane) using Illumina HiSeq 2000 (https://www.illumina.com) for paired-end reads, checked the quality of reads using FastQC version 0.11.2 (http://www.bioinformatics.babraham.ac.uk/projects/fastqc), and trimmed using trimmomatic version 0.27 (14). We mapped reads to reference genome M. tuberculosis H37Rv (GenBank accession no. NC_000962.3) using BWA-MEM (https://arxiv.org/abs/1303.3997) and used GATK UnifiedGenotyper (15) to call single-nucleotide polymorphisms (SNPs) and small insertion/deletions (indels). We selected the SNPs and small indels with >10 times read depth, 80% allele frequency, and >10-bp difference between neighboring SNPs and indels. Using SnpEff version 4.1 (16), we annotated high-quality SNPs and indels. We included for analysis mutations in all known genes (including regulatory genes) that conferred resistance to TB drugs (Appendix Table 1) and excluded mutations in repetitive regions, such as the PE and PPE gene family regions. We submitted raw reads in the form of FASTQ files to the Sequence Read Archive (project file no. PRJNA401368). Scripts of the raw sequence data are available (https://github.com/arnoldbaino/Daru_scripts).
Phylogenetics and Bayesian Coalescent Analysis
We used data sets from 2 previous independent WGS studies, PRJEB7281 and PRJEB2358 (17), as TB global representatives in phylogenetic analysis. Using concatenated SNP alignment, we constructed a maximum-likelihood phylogenetic tree with RAxML version 7.4.2 (18). We used a general time-reversible model, with rate heterogeneity accommodated by using discrete rate categories (i.e., GTRCAT algorithm), with 1,000 bootstraps and visualized using FigTree version 1.4.2 (http://tree.bio.ed.ac.uk/software/figtree). We performed molecular dating of the Beijing lineage isolates using BEAST version 1.8.2 (19), as previously described (7).
Transmission Assessment
Using ape library in R statistical package (http://cran.r-project.org), we calculated pairwise SNP differences (excluding SNPs in known drug-resistance genes) between isolates. We assessed SNP differences by lineage and compared these differences between different localities to set a threshold as a measure of putative transmission. We constructed genomic clusters using a median-joining network (Network version 5, http://www.fluxus-engineering.com/sharenet.htm). We structured links within genomic clusters assuming the existence of linear transmission among isolates from the same locality and putative transmission among isolates from different localities (if the SNP difference was within the limit of the threshold).
Statistical Analysis and Ethics
We compared patient characteristics using Fisher exact test. We performed univariate and multivariable logistic regression analyses to evaluate associations between patient characteristics and genomic clusters as outcome variables to infer transmission. We used R statistical package for all statistical analyses. This study was approved by the institutional review board of the University of Queensland (Brisbane, Queensland, Australia; approval no. HREC2015000572) and, as required by the Public Health Act of 2005, by Queensland Health (approval no. RD006697).


Thumbnail of Flow diagram of included Mycobacterium tuberculosis isolates from Papua New Guinea citizens residing in Torres Strait Protected Zone, 2010–2015. *Isolates unable to grow or were contaminated. †Included were 4 additional isolates among Queensland residents that were a part of an epidemiologic cluster linked to the Torres Strait Protected Zone. TB, tuberculosis; WGS, whole-genome sequencing.
Figure 2. Flow diagram of included Mycobacterium tuberculosis isolates from Papua New Guinea citizens residing in Torres Strait Protected Zone, 2010–2015. *Isolates unable to grow or were contaminated. †Included were 4 additional isolates...
In total, 134 Papua New Guinea citizens were reported to have TB in the TSPZ during the study period (Figure 2), and 117 (87.3%) had culture-confirmed disease. We acquired and successfully sequenced 100 (85.5%) of 117 isolates from these citizens plus 4 isolates from the Australia citizens who resided or previously resided in the TSPZ. Isolates were sequenced with a mean coverage depth of 74 times (range 42–110 times) and mean coverage breadth of 98.5% (range 96.3%–99.8%). MIRU profiles were available for 98 of 104 patients (Appendix Table 2). Sequence comparisons revealed a median SNP difference of 31 (interquartile range [IQR] 24–38) for isolates with the same MIRU profile and 1,101 (IQR 1,066–1,119) for isolates with different MIRU profiles (Appendix Figure 1). SNP differences between all possible pairs within the data set were bimodal; 2 large peaks represent 2 different lineages. The possibility of cross-contamination among isolates with no SNP differences was negligible because specimens were processed and their DNA sequenced on different days. Of the 104 isolates sequenced, phylogenetic analysis revealed that 83, including the 4 MDR TB isolates from Australia citizens, were part of the modern Beijing sublineage and 21 the Euro-American lineage, which consisted of 6 different sublineages (AppendixTable 3, Figure 2). Although the isolates were widely distributed throughout the TSPZ, 50% originated from a single Papua New Guinea village, Mabadauan (Appendix Figure 3).
Evolutionary analysis of the Beijing sublineage isolates revealed that 59% (49/83) were related to clades previously identified in nearby Daru Island (Appendix Figure 4) (7) and 41% were a part of unique clades. The isolates from the Australia residents were clade C, which we inferred to have emerged in the 1980s.
Thumbnail of Genomic clusters of highly drug-resistant tuberculosis cases among Papua New Guinea and Australia citizens residing or previously residing in the Torres Strait Protected Zone, 2010–2015, by lineage. The top shows the Beijing lineage, and the box at bottom shows the Euro-American lineage. Each circle represents &gt;1 isolate. Larger circles represent &gt;2 isolates with identical sequences; the number of specimens with identical sequences are indicated in parentheses. Details indicat
Figure 3. Genomic clusters of highly drug-resistant tuberculosis cases among Papua New Guinea and Australia citizens residing or previously residing in the Torres Strait Protected Zone, 2010–2015, by lineage. The top shows the...
To determine links between isolates from the same and different localities, especially isolates from villages with high numbers of isolates from different lineages, we deduced a threshold of 8 SNPs to distinguish transmission links (Appendix Figure 5). We found 17 genomic clusters (14 Beijing lineage, 3 Euro-American lineage) constituting 74 isolates (65 Beijing lineage, 9 Euro-American lineage) among PNG citizens (Figure 3). In total, 89% (34/38) of isolates from Mabadauan (identified as Beijing lineage) and 100% (4/4) of isolates from Sigabadaru (characterized as Euro-American lineage) clustered. More cases (36/41) formed genomic clusters in 2011 than in any other year, indicative of enhanced transmission leading up to this time point (Appendix Figure 6). The number of isolates from the same locality that formed genomic clusters (41 Beijing, 5 Euro-American) was significantly higher than the number that formed among isolates from different localities (24 Beijing, 4 Euro-American; p = 0.0013). The median SNP difference was 2 (IQR 1–3) among genomic clusters from the same locality and 4 (IQR 2–6) among genomic clusters from different localities. Of the 4 Australia cases, 3 formed 1 MDR TB cross-border genomic cluster having 1 SNP difference and a 15-month difference in sample collection dates. The fourth isolate was linked to another MDR TB cross-border cluster with no SNP differences and a 5-month difference in sample collection dates, suggesting 2 independent episodes of cross-border transmission into Australia citizens.
Most cross-border TB isolates in Papua New Guinea residents (74%, 74/100) were found in young persons (<35 years of age) (Table 1), and male and female sexes were equally represented. Only 1 (1.4%) of 74 patients tested had HIV co-infection. No major differences in characteristics were detectable between patients infected with the Beijing and Euro-American lineage, except for drug resistance, which was strongly associated with the Beijing lineage (p<0.001).
Univariate analysis revealed that the Beijing lineage (odds ratio [OR] 6.190, 95% CI 2.221–18.077; p<0.0005) and resistance to first- and second-line TB drugs (OR 4.677, 95% CI 1.452–20.943; p = 0.019) were strongly associated with transmission (Table 2). In multivariate logistic analysis, we adjusted for these 2 characteristics, and this analysis showed the Beijing lineage was associated with transmission (adjusted OR 4.484, 95% CI 1.526–13.891). The HIV-infected patient was not part of a cluster.
Phenotypic drug susceptibility testing of all isolates showed that 53 (51.0%) were susceptible to all oral first-line drugs tested (Appendix Table 4); 1 showed streptomycin monoresistance. Drug-resistance patterns observed included 13.5% (14/104) with isoniazid and streptomycin co-resistance 31.7% (33/104) with MDR (1 pre-XDR with fluoroquinolone resistance), and 1.9% (2/104) with XDR. In 2010 and 2011, eleven TB cases were MDR or XDR (Appendix Figure 6). All MDR and XDR isolates, including the isolates from the 4 Australia citizens, belonged to clades C and D of Beijing sublineage
Thumbnail of Correlation of phenotypic and genotypic drug resistance among modern Beijing isolates from Papua New Guinea and Australia citizens residing or previously residing in the Torres Strait Protected Zone, 2010–2015. No resistance mutation detected in 1 kanamycin-resistant strain, despite targeted sequencing of rrs and eis genes. *Pairs of co-occurring mutations observed within the same isolates were embA-embB (C12T) with embB (p.Asp354Ala) and embB (p.Gly406Val) with embB (p.Met306Ile).
Figure 4. Correlation of phenotypic and genotypic drug resistance among modern Beijing isolates from Papua New Guinea and Australia citizens residing or previously residing in the Torres Strait Protected Zone, 2010–2015. No resistance...
Phenotypic and genotypic resistance profiles correlated, except for ethambutol (Figure 4). Of the 35 MDR and XDR isolates, which were also tested for second-line drug susceptibility, all carried the fabG1-inhA mutation (C15T), and 28 (80%) had high-level (0.4 µg/mL) isoniazid resistance. In total, 27 of these 28 had both ndh (ΔG304) and inhA (p.Ile21Val) mutations, and 1 had a katG (p.Trp191Arg) mutation. All the rifampin-resistant isolates had the same rpoB mutation (p.Ser450Leu), and 29 (82.9%) of 35 had compensatory mutations in the rpoC gene (27 p.Val483Gly and 2 p.Trp484Gly). Of note, nearly all isolates with compensatory rpoC mutations (28/29) and isolates with the fabG1-inhA mutation (48/49) were part of a genomic cluster. Of the 2 isolates with the p.Asp94Ala gyrA mutation conferring fluoroquinolone resistance, 1 was XDR; this isolate had pan–second-line injectable resistance associated with an rrs (A1401G) mutation. The second XDR isolate (from a patient previously treated in Papua New Guinea) had a p.Asp89Asn gyrA mutation and phenotypic kanamycin resistance without a detected mutation conferring aminoglycoside resistance. Mutations in genes potentially conferring resistance to newly available agents were not identified.


We characterized M. tuberculosis strain diversity, transmission dynamics, and drug-resistance profiles among cross-border isolates from TB patients residing in the TSPZ. This region is known for its sparse population, unique geography, and ethnic and social diversity, all of which can make disease surveillance difficult. The identification of drug-resistant isolates in every studied year and evidence of cross-border links underscore the risk for cross-border importation of drug-resistant TB with the potential to spread. Studies conducted in other parts of Australia suggest limited local TB transmission, despite a high proportion of imported TB cases (20,21), but those findings mainly reflected occurrences in urban areas. Residents of and visitors to the TSPZ might be vulnerable to community outbreaks; cultural and family relationships were associated with TB transmission in studies of community outbreaks within Australia among indigenous populations (22). We show that most transmission in the TSPZ occurred during local community outbreaks, as demonstrated by the high genomic clustering among patients from the same locality. With the economic growth anticipated in this region (23), interactions between Australia and Papua New Guinea populations could increase and affect disease notification trends.
The Beijing sublineage dominated among isolates, and the Euro-American sublineage isolates were diverse. In a previous report, we attributed isolates from the same Beijing sublineage responsible for a large MDR TB outbreak on Daru Island, which is ≈50 km east of the outer islands of the TSPZ (7). The MDR TB outbreak on Daru Island was unprecedented in scale (24), but no published information described the extent of this outbreak beyond the shores of Daru Island.
In our study, most cases were identified to have come from Mabadauan, indicating that this setting could have a higher migratory rate than other TSPZ villages. The high number of drug-resistant TB cases from Mabadauan suggests this setting could be another hotspot for MDR TB transmission. Control of MDR TB transmission with effective clinical and public health–related interventions in this remote setting is urgently needed. Only 59% of Beijing lineage isolates from the TSPZ were part of clades previously identified on Daru Island (7), illustrating a greater genomic diversity for this lineage. This finding might be an indication of ongoing microevolution, which could further influence transmissibility, acquisition of drug resistance, or severity of disease (25).
Our study identified 2 plausible independent episodes of drug-resistant TB transmission to Australia residents in the TSPZ, a finding that would not have been identified by conventional genotyping techniques. Others have shown the superiority of WGS over conventional genotyping tools for resolving transmission (11,26). Two of the Australia patients with MDR TB were residing in major urban centers in northern Queensland at the time of diagnosis; 1 was living in a residential congregate setting, demonstrating the potential for diffuse community spread. WGS confirmed that the initial isolate noted in the index Australia patient was a lineage 4 isolate with only streptomycin resistance (gidB, Leu79Ser) (data not shown). This isolate was unrelated to the other Euro-American isolates identified and different from the subsequent MDR TB Beijing isolate identified, consistent with exogenous reinfection or endogenous reactivation rather than persistent infection with acquired resistance (27).
In another study conducted in Europe, WGS was used to investigate MDR TB outbreaks among immigrants and traced transmission routes to strains circulating in northern Somalia or Djibouti (28). Although no evidence of transmission to citizens of Europe was noted in that study, our investigation proves the principle of MDR TB cross-border transmission. Our findings highlight the challenge faced by Australia and Papua New Guinea TB control programs to prevent TB transmission through the TSPZ. Over the study period, the number of TB cases among citizens of Australia in the TSPZ was low (18 cases). Most TB notifications on the Australia island communities of the TSPZ were of diagnoses among Papua New Guinea citizens treated in outreach clinics located on the outer Torres Strait Islands close to the Papua New Guinea border. These clinics were closed in 2012, and all Papua New Guinea patients were handed over to the Western Province TB program for treatment, and enhanced TB services were coordinated from Daru Island. This closure and change in services accounts for the sharp fall in notifications in 2012 (Appendix Figure 6).
We observed transmission of a locally evolved Beijing sublineage strain that was associated with resistance to first- and second-line TB drugs. Beijing strains are known to be associated with increased transmissibility (2931), drug resistance, treatment failure, and rapid progression to active disease (32,33).
Although the World Health Organization estimates a high TB incidence among persons living with HIV in Papua New Guinea (44 cases/100,000 population) (1), HIV incidence in our study was low and did not account for the highly successful transmission of Beijing sublineage in this region. In another study in Gulf Province, Papua New Guinea, only 2 HIV-infected persons were identified among 105 TB patients (34).
The high proportion (74%) of TB cases among young adults (<35 years of age) in our study suggests ongoing TB transmission and is consistent with the finding of another study that identified 76.6% of pediatric TB notifications in Queensland were from Papua New Guinea residents in the TSPZ (35). Further assessment is needed of other risk factors, such as medical conditions (e.g., diabetes mellitus, high blood pressure) and socioeconomic status, that could affect TB transmission.
Excluding ethambutol, we observed a near-perfect correlation between resistance mutations and phenotypic in vitro resistance. One XDR TB isolate had resistance to kanamycin (mycobacterial growth indicator tube system critical concentration 2.5 µg/mL and a microbroth dilution MIC of 10 µg/mL), but no known resistance-conferring mutations were documented. Isolates with low-level kanamycin resistance are generally thought to display resistance to only kanamycin among the aminoglycoside drugs (36), but this isolate had a high level of kanamycin resistance. A previous study determined that only 80% of phenotypic kanamycin resistance could be accounted for by known mutations, indicating unexplained mechanisms of kanamycin resistance (37). Although the fabG1-inhA mutation is typically associated with low-level isoniazid resistance, we usually observed high-level isoniazid resistance with this mutation, possibly because of the accompanying inhA(p.Ile21Val) mutation (38). The fabG1-inhA (C15T) mutation does not seem to compromise successful spread and has previously been documented in clustered M. tuberculosis isolates in South Africa (39). One isolate had a fabG1-inhA (C15T) mutation with a rare katG mutation (p.Trp191Arg), which could affect the conformation of the catalase peroxidase protein, considering this mutation is located outside the active binding site of isoniazid (40).
Unexpectedly, we observed most rifampin-resistant isolates had mutations in rpoC. The fitness cost associated with some drug-resistant mutations can be ameliorated by compensatory mutations (41). Most isolates (96.5%) with rpoC compensatory mutations demonstrated genomic clustering, signifying their probable role in MDR or XDR TB transmission. In an assessment of MDR TB isolates in Argentina, no evidence was found to support the role of rpoC mutations in fitness restoration and increased transmission (42). Higher proportions of compensatory mutations in MDR TB outbreaks were reported in Shanghai, China (66%) (32), and Samara, Russia (87%) (11). More studies are needed to assess the role of compensatory mutations in MDR and XDR strain transmission.
This study had a few limitations. The retrospective study design involved the selection of only 1 isolate per patient, and some of the stored isolates were unable to be studied because of failure to recover from culture. Close social contact between patients with genomically linked isolates could not be confirmed. In other studies, a refined breakdown of transmission patterns was shown with a similar approach that included data on social contacts (28,32,43). Although the current scope of the study comprised mostly Papua New Guinea residents, additional information is needed regarding the evolution and transmission dynamics among all residents in the TSPZ, including drug-susceptible isolates in Australia citizens residing in the TSPZ, to assess the implications of transmission to populations on the Australia mainland.
The combination of genomic and epidemiologic data in this study highlighted the wide distribution of a Beijing sublineage strain associated with MDR and XDR TB in the TSPZ that threatens regional TB control through cross-border transmission. Further investigation of the environmental, sociocultural, and clinical factors that facilitate TB transmission in this region is warranted. Supplementation of regional TB surveillance programs with WGS technology could improve control strategies.
Dr. Bainomugisa is a research scientist at the University of Queensland and a part-time scientist at the Queensland Mycobacterium Reference Laboratory (World Health Organization collaborating center) in Brisbane, Queensland, Australia. His primary research interests include bacterial genomics, bioinformatics, and epidemiology of infectious diseases.


We thank the staff of Cairns TB control unit, Torres and Cape TB Control Unit, and Queensland Mycobacterium Reference Laboratory for their assistance. Derek Bensons’ assistance with data analysis is gratefully acknowledged. The dedication and cooperation of Papua New Guinea health officials working in TB control in the South Fly District and Port Moresby should be recognized.
This research was supported by NeCTAR Research Cloud, Queensland Cyber Infrastructure Foundation and University of Queensland’s Research Computing Centre. The NeCTAR Research Cloud is a collaborative Australian research platform supported by the National Collaborative Research Infrastructure Strategy.


  1. World Health Organization. Global tuberculosis report 2017. Geneva: The Organization; 2017.
  2. Uplekar  MWeil  DLonnroth  KJaramillo  ELienhardt  CDias  HMet al.for WHO’s Global TB ProgrammeWHO’s new end TB strategy. Lancet2015;385:1799801DOIPubMed
  3. Toms  CStapledon  RWaring  JDouglas  PNational Tuberculosis Advisory Committee, for the Communicable Diseases Network Australia, and the Australian Mycobacterium Reference Laboratory NetworkTuberculosis notifications in Australia, 2012 and 2013. Commun Dis Intell Q Rep2015;39:E21735.PubMed
  4. Lönnroth  KMor  ZErkens  CBruchfeld  JNathavitharana  RRvan der Werf  MJet al. Tuberculosis in migrants in low-incidence countries: epidemiology and intervention entry points. Int J Tuberc Lung Dis2017;21:62437DOIPubMed
  5. State of Queensland (Queensland Health). Tuberculosis in Queensland, 2016. 2017 Oct [cited 2018 Jun 17]. https://www.health.qld.gov.au/__data/assets/pdf_file/0028/679123/report-tb-qld-2016.pdf
  6. World Health Organization. Let’s kick TB out of PNG. 2018 Feb 13 [cited 2019 Jan 16]. http://www.wpro.who.int/papuanewguinea/mediacentre/releases/20180213-kick-tb-out-png/en/.
  7. Bainomugisa  ALavu  EHiashiri  SMajumdar  SHonjepari  AMoke  Ret al. Multi-clonal evolution of multi-drug-resistant/extensively drug-resistant Mycobacterium tuberculosis in a high-prevalence setting of Papua New Guinea for over three decades. Microb Genom2018;4DOIPubMed
  8. Butler  JRATawake  ASkewes  TTawake  LMcGrath  VIntegrating traditional ecological knowledge and fisheries management in the Torres Strait, Australia: the catalytic role of turtles and dugong as cultural keystone speciesEcol Soc2012;17:34DOI
  9. Australian Bureau of Statistics. Census community profiles. Torres Strait Island (R). Canberra. 2016 [cited 2018 Jun 17]. http://www.censusdata.abs.gov.au/census_services/getproduct/census/2016/communityprofile/LGA36960?opendocument
  10. State of Queensland (Queensland Health). Tuberculosis in Queensland, 1 January 2013–31 December 2014. 2016 May [cited 2018 Jun 17]. https://www.health.qld.gov.au/__data/assets/pdf_file/0033/444759/report-tb-qld-2013-14.pdf
  11. Casali  NNikolayevskyy  VBalabanova  YHarris  SRIgnatyeva  OKontsevaya  Iet al. Evolution and transmission of drug-resistant tuberculosis in a Russian population. Nat Genet2014;46:27986DOIPubMed
  12. Merker  MBlin  CMona  SDuforet-Frebourg  NLecher  SWillery  Eet al. Evolutionary history and global spread of the Mycobacterium tuberculosis Beijing lineage. Nat Genet2015;47:2429DOIPubMed
  13. Supply  PAllix  CLesjean  SCardoso-Oelemann  MRüsch-Gerdes  SWillery  Eet al. Proposal for standardization of optimized mycobacterial interspersed repetitive unit-variable-number tandem repeat typing of Mycobacterium tuberculosis. J Clin Microbiol2006;44:4498510DOIPubMed
  14. Bolger  AMLohse  MUsadel  BTrimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics2014;30:211420DOIPubMed
  15. McKenna  AHanna  MBanks  ESivachenko  ACibulskis  KKernytsky  Aet al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data.Genome Res2010;20:1297303DOIPubMed
  16. Cingolani  PPlatts  AWang  LCoon  MNguyen  TWang  Let al. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly (Austin)2012;6:8092DOIPubMed
  17. Guerra-Assunção  JACrampin  ACHouben  RMMzembe  TMallard  KColl  Fet al. Large-scale whole genome sequencing of M. tuberculosis provides insights into transmission in a high prevalence area. eLife2015;4:e05166DOIPubMed
  18. Stamatakis  AHoover  PRougemont  JA rapid bootstrap algorithm for the RAxML Web servers. Syst Biol2008;57:75871DOIPubMed
  19. Drummond  AJSuchard  MAXie  DRambaut  ABayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol2012;29:196973DOIPubMed
  20. Gurjav  UOuthred  ACJelfs  PMcCallum  NWang  QHill-Cawthorne  GAet al. Whole genome sequencing demonstrates limited transmission within identified Mycobacterium tuberculosisclusters in New South Wales, Australia. PLoS One2016;11:e0163612DOIPubMed
  21. Gurjav  UJelfs  PHill-Cawthorne  GAMarais  BJSintchenko  VGenotype heterogeneity of Mycobacterium tuberculosis within geospatial hotspots suggests foci of imported infection in Sydney, Australia. Infect Genet Evol2016;40:34651DOIPubMed
  22. Outhred  ACHolmes  NSadsad  RMartinez  EJelfs  PHill-Cawthorne  GAet al. Identifying likely transmission pathways within a 10-year community outbreak of tuberculosis by high-depth whole genome sequencing. PLoS One2016;11:e0150550DOIPubMed
  23. Torres Strait Regional Authority. Torres Strait Regional economic investment strategy, 2015–2018. 2015 Sep [cited 2018 Jun 17]. http://www.tsra.gov.au/__data/assets/pdf_file/0007/17755/SQUIZ-Torres-Strait-Regional-Economic-Investment-Strategy.pdf
  24. Furin  JCox  HOutbreak of multidrug-resistant tuberculosis on Daru Island. Lancet Respir Med2016;4:3479DOIPubMed
  25. Mokrousov  IInsights into the origin, emergence, and current spread of a successful Russian clone of Mycobacterium tuberculosis. Clin Microbiol Rev2013;26:34260DOIPubMed
  26. Merker  MKohl  TARoetzer  ATruebe  LRichter  ERüsch-Gerdes  Set al. Whole genome sequencing reveals complex evolution patterns of multidrug-resistant Mycobacterium tuberculosisBeijing strains in patients. PLoS One2013;8:e82551DOIPubMed
  27. Hatherell  HAColijn  CStagg  HRJackson  CWinter  JRAbubakar  IInterpreting whole genome sequencing for investigating tuberculosis transmission: a systematic review. BMC Med2016;14:21DOIPubMed
  28. Walker  TMMerker  MKnoblauch  AMHelbling  PSchoch  ODvan der Werf  MJet al.MDR-TB Cluster ConsortiumA cluster of multidrug-resistant Mycobacterium tuberculosis among patients arriving in Europe from the Horn of Africa: a molecular epidemiological study. Lancet Infect Dis2018;18:43140DOIPubMed
  29. Buu  TNvan Soolingen  DHuyen  MNLan  NTQuy  HTTiemersma  EWet al. Increased transmission of Mycobacterium tuberculosis Beijing genotype strains associated with resistance to streptomycin: a population-based study. PLoS One2012;7:e42323DOIPubMed
  30. Cowley  DGovender  DFebruary  BWolfe  MSteyn  LEvans  Jet al. Recent and rapid emergence of W-Beijing strains of Mycobacterium tuberculosis in Cape Town, South Africa. Clin Infect Dis2008;47:12529DOIPubMed
  31. Yang  CLuo  TSun  GQiao  KSun  GDeRiemer  Ket al. Mycobacterium tuberculosis Beijing strains favor transmission but not drug resistance in China. Clin Infect Dis2012;55:117987DOIPubMed
  32. Yang  CLuo  TShen  XWu  JGan  MXu  Pet al. Transmission of multidrug-resistant Mycobacterium tuberculosis in Shanghai, China: a retrospective observational study using whole-genome sequencing and epidemiological investigation. Lancet Infect Dis2017;17:27584DOIPubMed
  33. Parwati  Ivan Crevel  Rvan Soolingen  DPossible underlying mechanisms for successful emergence of the Mycobacterium tuberculosis Beijing genotype strains. Lancet Infect Dis2010;10:10311DOIPubMed
  34. Cross  GBColes  KNikpour  MMoore  OADenholm  JMcBryde  ESet al. TB incidence and characteristics in the remote gulf province of Papua New Guinea: a prospective study. BMC Infect Dis2014;14:93DOIPubMed
  35. Donnan  EJCoulter  CSimpson  GClark  JNourse  CPaediatric tuberculosis in Queensland, Australia: overrepresentation of cross-border and Indigenous children. Int J Tuberc Lung Dis2017;21:2639DOIPubMed
  36. Reeves  AZCampbell  PJSultana  RMalik  SMurray  MPlikaytis  BBet al. Aminoglycoside cross-resistance in Mycobacterium tuberculosis due to mutations in the 5′ untranslated region of whiB7. Antimicrob Agents Chemother2013;57:185765DOIPubMed
  37. Zaunbrecher  MASikes  RD JrMetchock  BShinnick  TMPosey  JEOverexpression of the chromosomally encoded aminoglycoside acetyltransferase eis confers kanamycin resistance in Mycobacterium tuberculosis. Proc Natl Acad Sci U S A2009;106:200049DOIPubMed
  38. Zhang  YYew  WWMechanisms of drug resistance in Mycobacterium tuberculosis: update 2015. Int J Tuberc Lung Dis2015;19:127689DOIPubMed
  39. Johnson  RWarren  RMvan der Spuy  GDGey van Pittius  NCTheron  DStreicher  EMet al. Drug-resistant tuberculosis epidemic in the Western Cape driven by a virulent Beijing genotype strain. Int J Tuberc Lung Dis2010;14:11921.PubMed
  40. Bertrand  TEady  NAJones  JNJesmin Nagy  JMJamart-Grégoire  Bet al. Crystal structure of Mycobacterium tuberculosis catalase-peroxidase. J Biol Chem2004;279:389919DOIPubMed
  41. Comas  IBorrell  SRoetzer  ARose  GMalla  BKato-Maeda  Met al. Whole-genome sequencing of rifampicin-resistant Mycobacterium tuberculosis strains identifies compensatory mutations in RNA polymerase genes. Nat Genet2011;44:10610DOIPubMed
  42. Eldholm  VMonteserin  JRieux  ALopez  BSobkowiak  BRitacco  Vet al. Four decades of transmission of a multidrug-resistant Mycobacterium tuberculosis outbreak strain. Nat Commun2015;6:7119DOIPubMed
  43. Walker  TMIp  CLHarrell  RHEvans  JTKapatai  GDedicoat  MJet al. Whole-genome sequencing to delineate Mycobacterium tuberculosis outbreaks: a retrospective observational study.Lancet Infect Dis2013;13:13746DOIPubMed



Cite This Article

DOI: 10.3201/eid2503.181003

No hay comentarios:

Publicar un comentario