Diversity of Mycobacterium tuberculosis and drug resistance in different provinces of Papua New Guinea
© Ley et al.; licensee BioMed Central Ltd. 2014
Received: 28 July 2014
Accepted: 24 November 2014
Published: 5 December 2014
Papua New Guinea (PNG) is a high tuberculosis (TB) burden country of the WHO Western Pacific Region, but so far research on drug resistance (DR) and genotypes of Mycobacterium tuberculosis (M. tuberculosis) was only conducted in few provinces in the country. The aim of the present study was to obtain baseline data on the level of drug resistance and the genotypic diversity of circulating M. tuberculosis in additional provinces and to investigate the differences between three selected sites across PNG.
Genotyping of 147 M. tuberculosis clinical isolates collected in Goroka, Eastern Highlands Province, in Alotau, Milne Bay Province and in Madang, Madang Province revealed three main lineages of M. tuberculosis: Lineage 4 (European-American lineage), Lineage 2 (East-Asian lineage) and Lineage 1 (Indo-Oceanic lineage). All three lineages were detected in all three sites, but the individual lineage compositions varied significantly between sites. In Madang Lineage 4 was the most prevalent lineage (76.6%), whereas in Goroka and Alotau Lineage 2 was dominating (60.5% and 84.4%, respectively) (p < 0.001). Overall, phenotypic drug susceptibility testing showed 10.8% resistance to at least one of the first-line drugs tested. Of all resistant strains (23/212) 30.4% were Streptomycin mono-resistant, 17.4% were Isoniazid mono-resistant and 13% were Rifampicin mono-resistant. Multi-drug resistant (MDR) TB was found in 2.8% of all tested cases (6/212). The highest amount of MDR TB was found in Alotau in Milne Bay Province (4.6%).
A large number of drug resistant TB infections are present in the country and MDR TB has already been detected in all three surveyed regions of PNG, highlighting the importance of monitoring drug resistance and making it a high priority for the National Control Program. Due to the high prevalence of Lineage 2 in Milne Bay Province and given the frequent association of this lineage with drug resistance, monitoring of the latter should especially be scaled up in that province.
Over the last decade, evidence for the impact of the bacterial genetic background on TB infection and disease has strongly increased . A phylogeography of M. tuberculosis based on large sequence polymorphisms and confirmed by multi-locus sequence analysis could furthermore be established, showing an association between specific M. tuberculosis strains and a particular geographic region ,. The different M. tuberculosis strains were grouped into 6 main lineages (a 7th was recently added ). Various studies investigated the impact of different MTBC lineages on the clinical presentation of the disease. Infection with Lineage 2 has for example been found to be associated with faster progression to disease  and drug resistance ,. Lineage 4 on the other hand, has been associated with pulmonary TB rather than extrapulmonary TB . In addition, few studies have analysed the impact of the genetic background of both, humans and bacteria on disease development , and found correlations between different human genetic polymorphisms with specific M. tuberculosis lineages. However, findings from studies investigating lineage specific associations have also reflected within lineage variability ,. Therefore, further discriminatory strain differentiation methods should be considered when investigating M. tuberculosis genetic diversity. Nevertheless, the various findings support the idea of a longstanding host-pathogen co-evolution and the hypothesis that TB spread together with the human out-of-Africa migration ,,. To investigate the M. tuberculosis genetic diversity within and between different populations could therefore give important insights into the dynamics of TB disease and might help to inform national TB programs to develop better control strategies.
Papua New Guinea (PNG) is a high TB burden country with an estimated TB incidence rate of 348/100’000 in 2012 and its proportion of multidrug resistance (MDR) estimated at 4.9% in new cases is higher than the estimated global average of 3.6% ,. Publications on drug resistant TB in PNG are rare, and previous data were mainly derived from patients from Western Province, diagnosed in Australia -. Gilpin et al. reported 25% MDR-TB among patients from Western Province diagnosed between 2000 and 2006 . Only two publications on DR data from other provinces exist: Ballif et al. found that 5.2% of tested isolates of adult patients from Madang were MDR-TB , and a recently published study from Kikori in Gulf Province of PNG reported 9% of suspected MDR-TB (based on Xpert® MTB/RIF) in the investigated population .
PNG harbours a vast human genetic diversity and has been isolated from the outside world for a long time. PNG was populated through several waves of human migrations and people living in the Highlands represent the oldest population from the first migration wave about 50,000 years ago ,. The highlands region has only been ‘discovered’ in the 1930s  and was sparsely populated at the time. Hence it would be expected that evolutionary ‘ancient’ lineages of M. tuberculosis (e.g. Lineage 1 ), proposed to be adapted to low density populations , would be found there, whilst the ‘modern’ lineages such as Lineage 2 and Lineage 4 would be expected at the highly populated coastal regions. However, there is limited information on the TB situation available in PNG. Since various lineages have been found to differ in their prevalence between different regions of the same country (e.g. in Indonesia  or Taiwan ) it is important to also identify the circulating strains in various communities in PNG. To our knowledge, apart from a previous study conducted in Madang province , another study where isolates from patients from Western Province were analysed  and a very recent study from Gulf Province , no other data on the M. tuberculosis population structure from PNG have been published.
The aim of the present study was to obtain baseline data on the level of drug resistance and the genotypic diversity of circulating M. tuberculosis in selected sites of three provinces of PNG (Goroka in Eastern Highlands Province, Alotau in Milne Bay Province and Madang in Madang Province) and to investigate the differences between these sites.
Results and discussion
Study population characteristics
Characteristics of study population for each study site
Drug susceptibility testing result available*
Total TB positive (any type)
Goroka n = 56 n (%)
Alotau n = 66 n (%)
Madang n = 90 n (%)
Total n = 212 n (%)
All sites n = 449 n(%)
Median age [IQR]
Region of origin
History of TB treatment
TB positivity peaked in the age group of 15 to 34 years, which is consistent with country wide data . And 79.1% (355/449) of patients clinically diagnosed with TB had no history of the disease indicating ongoing transmission. However, these data cannot be used to extrapolate for the whole country and further investigations are required to analyse the transmission dynamics of the disease in PNG.
Figure 1 depicts the flow of all isolates collected. Of samples sent for culture 67.2% (225/335) could be recovered showing a positive correlation between initial diagnostic bacterial count and culture success (likelihood ratio χ2 = 140.6, df = 4, p < 0.001). The reduced culture success might be likely due to the long distance transport of samples between laboratories that was experienced with samples from other studies.
One of the Lineage 1 strains from Goroka exhibits a special spoligotyping pattern with only spacers 40–43 being present (sample No. 4 in Figure 3). That spoligotyping pattern has been confirmed with in silico spoligotyping using KvarQ . Per standard definition, this would define that strain as belonging to the Beijing family of Lineage 2 ,. However, SNP-typing assigned this strain to Lineage 1, therefore representing a ‘pseudo-Beijing’ strain with convergent evolution leading to an independent deletion event in the direct repeat locus of MTBC, similar to what has already been described for strains of Lineage 3 .
All Lineage 4 samples belonging to the LAM spoligotyping family were exclusively found in Madang, also reflecting the different M. tuberculosis lineage composition in the spoligotyping data. All Lineage 2 samples belonged to the Beijing family and were subsequently subtyped into three different monophyletic groups which are sublineages defined by the presence or absence of specific regions of difference (RD) . Over all three sites, one sample of 67 Beijing strains belonged to sublineage 1, here defined as having no deletion of RD181, RD150 or RD142. One sample belonged to sublineage 2 with RD181 deleted but RD150 and RD142 being present, but the majority of Beijing strains (63/65), including all Beijing strains from Alotau, belonged to sublineage 3 with RD180 and RD150 deleted. None of the samples had a deletion for RD142. Sublineage 3 is usually rare, hence the high frequency in our sample is surprising. Several studies have shown that sublineage 3 has rapidly spread in Cape Town, South Africa during the last decade, probably because of a founder effect and/or adaptation to the local host-population -. Due to the lack of longitudinal data, it is not possible to draw any conclusion about the time of introduction or the duration of spread of this subtype in PNG. However, except for two outliers, our subtyping data based on RDs implies a single introduction of the Beijing strain and a subsequent clonal expansion through PNG. Starting with this observation, it would be important to follow the prevalence of this lineage in Madang and in the country as a whole to monitor a potential increase of Beijing strains in the future.
The reasons for the observed differences in lineage frequencies in the three study sites surveyed remain speculative. These sites differ in several aspects, e.g. Goroka in the Highlands had been much longer isolated from the outside world for much longer than the two coastal sites. The population in the Highlands are believed to be descendants of the oldest human migration wave that populated the country , however, no statistically significant difference in the prevalence of Lineage 1 in Goroka could be detected.
Univariate and multivariate logistic regressions for the risk of an infection with a Beijing strain of lineage 2
Lineage 2 n = 67 (45.58%)
Non-lineage 2 n = 80 (54.42%)
Univariate regression: associations withM. tuberculosislineage 2
Multivariate regression: associations withM. tuberculosislineage 2
Crude OR (95% CI)
Adj. OR (95% CI)
Study site n (%)
3.5 (1.11 - 11.18)
3.3 (1.0 - 10.8)
0.2 (0.08 - 0.43)
0.2 (0.1 - 0.4)
Sex, n (%)
1 (0.52 - 1.91)
Age - median years [IQR]
1 [0.99 - 1.05]
1.0 [0.98 - 1.1]
12.74 (2.57 -
TB treatment history, n (%)
0.60 (0.19 - 1.86)
Drug resistance (any type)
2.44 (0.85 - 7.02)
2.5 (0.7 - 8.5)
Co-infection with HIV has several times been shown to be associated with infection with strains of the Beijing type ,. We also found a significant association between HIV positivity and infection with a Lineage 2 strain (p = 0.002), but due to sample size multivariate logistic regression could not be performed. However, HIV prevalence in PNG is still comparably low, in particular in Alotau (2.1% in Milne Bay Province (Alotau), 5.2% in Eastern Highlands Province (Goroka) and 2.8% in Madang Province)  and is therefore unlikely to play a major role in the distribution of lineages.
It is important to note that all differences are between the frequencies of Lineage 2 and Lineage 4, whilst Lineage 1 played no major role, although we would have expected more ancient lineages (Lineage 1) in the highlands. Lineage 4 and Lineage 2 are not only the most prevalent M. tuberculosis lineages in PNG, but are also predominating globally . Thus, our findings could support the notion of Hershberg and colleagues  that the lineage distribution between and within countries might become homogenized with increasing migration. In other words, more virulent strains such as Lineage 4 and Lineage 2 might slowly replace the ancient lineages as has already been observed in Cameroon  and might have started in PNG.
Observed phenotypic drug resistance per study site
GKA n = (%)
ALO n = (%)
MAG n = (%)
TOTAL n = (%)
INH (0.1 mg/L)
INH (0.4 mg/L)
STR + INH
STR + RMP
INH + RMP
STR + INH + RMP + EMB
STR + INH + RMP
STR + INH + RMP + ETH
INH + RMP + PZA + ETH
Streptomycin resistance has been shown to be associated with mutations in the rrs, rpsL or gidB gene ,. For two STR mono-resistant samples only a mutation in gidB was detected (sample 1 in Figure 4) or no mutation was found in any of the analysed genes, respectively (sample 2 in Figure 4). All observed mutations in that gene had also been observed in our previous study in Madang. The sole observation of the synonymous mutation A205G and the non-synonymous mutation A92C in Lineage 2 strains, suggests that these mutations might be lineage specific , with no mutations occurring in Lineage 4. In three samples, all belonging to Lineage 4, we found the A10P mutation. In contrast to our previous study where the A10P mutation was absent from 21 pan-susceptible samples, we observed this mutation in one of the pan-susceptible samples (sample 17 in Figure 4), probably suggesting that this mutation plays no role in STR resistance. Whether the mutation V77G found here in one STR resistant strain - and yet only described from PNG - is involved in STR resistance  remains to be determined.
Ninety-five per cent of the rifampicin (RMP) resistance conferring mutations occur in an 81 bp core region of the so called rifampin resistance determining region (RRDR) . We found only one RMP monoresistant sample (sample 9 in Figure 4) which had no mutation in the 849 bp (including the RRDR) of the rpoB gene we amplified, confirming that RMP resistant strains without a typical mutation in the RRDR are not more common in our sample. This is of crucial importance for the PNG DR surveillance being based on Xpert® MTB/RIF (Cepheid)  which determines RMP resistance only through detection of mutations in RRDR . However, 26.1% (6/23) of strains were isoniazid (INH) monoresistant or INH/STR polyresistant. INH resistance is a precursor to MDR-TB and is not detected by the Xpert® MTB/RIF, forming an additional challenge for the control of DR TB that should be addressed by the NTP in the future.
All low level INH resistant samples (resistant to a concentration of 0.1 mg/L INH) of our study showed a mutation in the inhA promoter region, whereas all high INH resistant samples (0.4 mg/L) showed a mutation at codon 315 of katG, including all MDR samples and one of the polyresistant strains with INH and STR resistance (sample 11 in Figure 4). No mutations were detected in the ahpC promoter region for any of the samples.
Two MDR samples had no mutation in 850 bp of katG sequenced, none in the ahpC promoter or the inhA promoter. One of these samples was also resistant to ethambutol (EMB) (sample 16 in Figure 4) but had no mutation in the embB region including codon 306, which is mutated in up to 68% of clinical EMB resistant strains . For that sample the whole genome sequence was available (data not shown) and was used to screen for mutations outside the amplified regions of katG or embB. Screening revealed a 14 bp deletion at position 2156047 to 2156060 (H37Rv reference, GeneBank AL123456) causing a truncation of katG through a frameshift, explaining the phenotypic high INH resistance of that sample. Outside of the sequenced embB region the non-synonymous mutation G406S was found which had already been described in several other studies suggesting its role in EMB resistance ,. For the second MDR sample without a katG mutation, no sequence data was available and the INH resistance conferring mutation could not be determined. It remains unclear whether a deletion in katG or mutations in other genes associated with INH resistance, for example in kasA could be responsible for the INH resistance.
No mutation was found in the amplified regions of pncA for any of the genotyped isolates. However, since for the only sample with phenotypic pyrazinamide resistance (Table 3) no DNA could be obtained, genotyping of that isolate was not possible.
Univariate and multivariate logistic regressions for the risk of drug resistance
Pan-susceptible n = 189 (89.15%)
Resistant n = 23 (10.85%)
Univariate regression: associations with DR
Multivariate regression: associations with DR
Crude OR (95% CI)
Adj. OR (95% CI)
Study site n (%)
1.1 (0.4 - 3.5)
0.9 (0.3 - 2.8)
Sex, n (%)
1.0 (0.4 - 2.5)
Age - median years [IQR]
1.0 [1.0 - 1.1]
1.4 (0.3 - 5.7)
TB treatment history, n (%)
5.5 (1.8 - 16.8)
4.2 (1.1 - 16.7)
2.44 (0.85 - 7.02)
3.4 (1.0 - 11.2)
Ballif et al. previously found a significant association between an infection with a Lineage 2 strain and drug resistance in a previous study conducted in Madang (p < 0.010, OR = 5.2, CI (95%): 1.8 - 15.1) . Also in our current sample set from Madang a significant association between Lineage 2 strains and drug resistance was found, but only after correction for previous TB treatment (p = 0.041, OR = 3.4, CI (95%): 1.0 – 11.2), probably due to the limited sample size (Table 4). To test whether the two sample sets from Madang (Ballif et al. versus Madang samples from current study) differed significantly in the drug resistance data, we compared the results of the drug resistance and Lineage 2 analyses from both sample sets by a χ2 test of ORs, but no significant difference was found (χ2 = 0.747).
Because of the small sample size we also only found a borderline significant association between multiple drug resistance and Lineage 2 infection (p = 0.058; CI (95%) 0.9 – 68.4; OR = 8): there was an 8 times higher risk of being infected with a polyresistant (resistant to more than one drug but not MDR) or MDR strain when infected with an Lineage 2 strain compared to an infection with a strain of a different lineage (in this case Lineage 4 or Lineage 1).
The direct comparison of M tuberculosis population structures from distinct sites in PNG demonstrated a statistically significant difference between subpopulations. In Madang Lineage 4 was the dominating lineage, whereas Lineage 2 was more frequently detected in Alotau and Goroka. Although the reasons for the observed significant differences of the circulating M. tuberculosis strains between study sites are not yet understood, these differences might have a major impact on disease and transmission dynamics in different populations of PNG. Different control strategies for places with a different M. tuberculosis lineage composition are not available yet, i.e. the same control strategies apply for all provinces in PNG, namely to detect cases and treat them accordingly. However, by knowing about the increased prevalence of the Beijing type of M. tuberculosis in Milne Bay Province, and with the known association between this lineage and drug resistance, monitoring of the latter should especially be scaled up in that province, as it could become a hot spot for drug resistance and MDR TB.
Our data show that a significant number of drug resistant TB infections are present across the country and that MDR TB can already be detected in all three surveyed regions of PNG. Nearly all phenotypical resistances were confirmed by sequence analysis.
No inferences can be made from this study for the whole country because of the small sample size and data being derived from only three major towns. Nevertheless, our findings highlight the importance to monitor drug resistance in PNG, and for making it a high priority for the National TB Control Program.
Study sites and patient characteristics
The study was conducted in three different sites across PNG, one site in each region of PNG: in Madang, Madang Province in the Momase Region; Goroka, Eastern Highlands Province, in the Highlands Region, and in Alotau, Milne Bay Province in the Southern Region of PNG. In Madang, patients were enrolled into the study from November 2010 onwards. In Goroka, patient enrolment started in June 2011 and Alotau was added as a study site in July 2011. In all three study sites enrolment was completed in July 2012. Three consecutive sputum samples were collected from adult TB suspect patients (15 years or older) with chronic productive cough who presented at any department of the provincial hospitals (Modilon Hospital, in Madang; Goroka Provincial Hospital in Goroka; Alotau Provincial Hospital and Gurney Health Centre in Alotau). Questionnaire based interviews were conducted to obtain socio-demographic and behavioural information of each patient. TB was diagnosed by either direct smear light microscopy (Ziehl-Neelson staining), fluorescent microscopy (Morse Stain; TB Fluorescent Stain Kit M, Becton, Dickinson and Company, USA), chest X-ray, clinical examination or a combination of these methods. All TB positive study patients were automatically enrolled into the PNG National TB Program (NTP). Therefore, patient management, i.e. treatment and follow up procedures, was carried out according to the NTP guidelines . For a subset of study patients the HIV status could be obtained from the NTP, which recommends HIV testing of TB positive patients.
Sample processing and drug susceptibility testing
Sputa were obtained from all tuberculosis patients with pulmonary involvement who were able to produce sputum. Sputum samples were decontaminated according to Petroff’s method . Subsequently, these samples were inoculated into Mycobacterial Growth Indicator Tubes (BACTEC™ MGIT™ 960 system; BD, Franklin Lakes, NJ, USA) and sent to the Queensland Mycobacterium Reference Lab in Brisbane, Australia, for culture and drug susceptibility testing (DST). DST utilising the BACTEC 960 MGIT system was conducted as described previously .
Genotyping of Mycobacterium tuberculosis
DNA was extracted from culture either by InstaGene Matrix (Bio-Rad, Hercules, CA, USA) following the manufacturer’s protocol, with one bacterial colony resuspended in 1 ml of dH2O as starting material, or by mixing 100 μl of M. tuberculosis inoculated into Dubos broth (prepared tubed medium for cultivation of mycobacteria; Becton, Dickinson and Company, Maryland, USA) with 100 μl distilled water and subsequent heat killing at 90°C for 1 hour. The DNA was then used for molecular analyses.
M. tuberculosis isolates were classified into the main phylogenetic lineages  by a TaqMan real-time PCR assay using single nucleotide polymorphism (SNP) typing . Lineages were further discriminated into families by spoligotyping , using a commercial membrane produced by Ocimum Biosolutions Ltd, India (product IM9702). Information on the shared international type (SIT) and the spoligo family were obtained from SITVIT WEB . All Beijing strains were further sub-classified into monophyletic groups based on the presence or absence of the regions of difference (RDs) RD181, RD150 and RD142 as described by Tsolaki et al. . PCR conditions and primers were used as previously described by Gagneux et al..
Drug resistance genotyping
Of all phenotypically drug resistant strains for which DNA could be obtained, drug resistance associated regions of the following genes were amplified by PCR and sequenced by Macrogen (The Netherlands): katG, rpoB, ahpC (promoter), inhA (promoter), gidB, pncA, gyrA, rrs, rpsL, embB. Additionally, two fully susceptible isolates of patients for whom smear conversion took longer than the usual 2 months were included. Primers and PCR conditions were used as previously described by Ballif et al. with the following modifications; for rpoB a new set of primers was designed: forward primer 5′AYATCGACCACTTCGGYAACC3′, reverse primer 5′TCCTCGATGACGCCGCTTTCT3′ (Y = C/T). PCR was run with an annealing temperature (TA) of 62°C, an elongation time (ET) of 60 seconds and 37 cycles, leading to a product length of 849 bp. For the inhA promoter and katG amplification, primers remained the same as published but the TA was increased from 60°C to 65°C and the cycle number from 35 to 39 for inhA, and from 64°C to 66°C and from 35 to 40 cycles for katG.
Statistical analysis was carried out with Stata 12.1 (StataCorp, College Station, TX, USA). Differences between study sites were assessed by cross-tabulation and significance testing using Fisher’s exact and χ2 testing. P-values <0.05 were considered statistically significant. Univariate logistic regressions were performed to assess associations of drug resistance or Lineage 2 infection with known potential predictors or confounders. Independent variables with a significance level of p < 0.2 in the univariate analysis and a plausible causal link where further analysed in a multivariate logistic regression. Model selection for multivariate regressions was based on the AIC-criterion. The comparison of two different sample sets from the same study site was done with a χ2 test of odds ratios with χ2 > 3.84 considered to show a statistically significant difference (5% level).
Ethical approval for this study was granted by the PNG IMR Institutional Review Board (IRB No. 0913) and the PNG Medical Research Advisory Council (MRAC No. 10.02). The Ethik-Kommission beider Basel (EKBB) has been informed and had approved the study. Written informed consent was obtained from all study participants.
We thank all study participants whose samples were collected and analysed as well as the health authorities of Madang Province, Alotau Province and Goroka Province for approving the study and establishing the contact to the hospitals and health centres. We are indebted to the PNG IMR study teams and the hospital staff in all three sites for sample collection, diagnosis, access to infrastructure and general support. The contribution of the staff of the Queensland Mycobacterium Reference Laboratory is gratefully acknowledged. The study was funded by the Stanley Thomas Johnson Foundation and the Medicor Foundation Liechtenstein.
- Coscolla M, Gagneux S: Does M. tuberculosis genomic diversity explain disease diversity?. Drug Discov Today Dis Mech. 2010, 7: e43-e59. 10.1016/j.ddmec.2010.09.004.PubMed CentralView ArticlePubMedGoogle Scholar
- Hershberg R, Lipatov M, Small PM, Sheffer H, Niemann S, Homolka S, Roach JC, Kremer K, Petrov DA, Feldman MW, Gagneux S: High functional diversity in Mycobacterium tuberculosis driven by genetic drift and human demography. PLoS Biol. 2008, 6: e311-10.1371/journal.pbio.0060311.PubMed CentralView ArticlePubMedGoogle Scholar
- Gagneux S, Deriemer K, Van T, Kato-Maeda M, de Jong BC, Narayanan S, Nicol M, Niemann S, Kremer K, Gutierrez MC, Hilty M, Hopewell PC, Small PM: Variable host-pathogen compatibility in Mycobacterium tuberculosis. Proc Natl Acad Sci U S A. 2006, 103: 2869-2873. 10.1073/pnas.0511240103.PubMed CentralView ArticlePubMedGoogle Scholar
- Firdessa R, Berg S, Hailu E, Schelling E, Gumi B, Erenso G, Gadisa E, Kiros T, Habtamu M, Hussein J, Zinsstag J, Robertson BD, Ameni G, Lohan AJ, Loftus B, Comas I, Gagneux S, Tschopp R, Yamuah L, Hewinson G, Gordon SV, Young DB, Aseffa A: Mycobacterial lineages causing pulmonary and extrapulmonary tuberculosis, Ethiopia. Emerg Infect Dis. 2013, 19: 460-463. 10.3201/eid1903.120256.PubMed CentralView ArticlePubMedGoogle Scholar
- de Jong BC, Hill PC, Aiken A, Awine T, Antonio M, Adetifa IM, Jackson-Sillah DJ, Fox A, Deriemer K, Gagneux S, Borgdorff MW, McAdam KP, Corrah T, Small PM, Adegbola RA: Progression to active tuberculosis, but not transmission, varies by Mycobacterium tuberculosis lineage in The Gambia. J Infect Dis. 2008, 198: 1037-1043. 10.1086/591504.PubMed CentralView ArticlePubMedGoogle Scholar
- Bifani PJ, Mathema B, Liu Z, Moghazeh SL, Shopsin B, Tempalski B, Driscol J, Frothingham R, Musser JM, Alcabes P, Kreiswirth BN: Identification of a W variant outbreak of Mycobacterium tuberculosis via population-based molecular epidemiology. JAMA. 1999, 282: 2321-2327. 10.1001/jama.282.24.2321.View ArticlePubMedGoogle Scholar
- Van Rie A, Warren RM, Beyers N, Gie RP, Classen CN, Richardson M, Sampson SL, Victor TC, van Helden PD: Transmission of a multidrug-resistant Mycobacterium tuberculosis strain resembling “strain W” among noninstitutionalized, human immunodeficiency virus-seronegative patients. J Infect Dis. 1999, 180: 1608-1615. 10.1086/315054.View ArticlePubMedGoogle Scholar
- Caws M, Thwaites G, Dunstan S, Hawn TR, Lan NT, Thuong NT, Stepniewska K, Huyen MN, Bang ND, Loc TH, Gagneux S, van Soolingen D, Kremer K, van der Sande M, Small P, Anh PT, Chinh NT, Quy HT, Duyen NT, Tho DQ, Hieu NT, Torok E, Hien TT, Dung NH, Nhu NT, Duy PM, van Vinh CN, Farrar J: The influence of host and bacterial genotype on the development of disseminated disease with Mycobacterium tuberculosis. PLoS Pathog. 2008, 4: e1000034-10.1371/journal.ppat.1000034.PubMed CentralView ArticlePubMedGoogle Scholar
- van Crevel R, Parwati I, Sahiratmadja E, Marzuki S, Ottenhoff TH, Netea MG, van der Ven A, Nelwan RH, van der Meer JW, Alisjahbana B, van de Vosse E: Infection with Mycobacterium tuberculosis Beijing genotype strains is associated with polymorphisms in SLC11A1/NRAMP1 in Indonesian patients with tuberculosis. J Infect Dis. 2009, 200: 1671-1674. 10.1086/648477.View ArticlePubMedGoogle Scholar
- Kato-Maeda M, Shanley CA, Ackart D, Jarlsberg LG, Shang S, Obregon-Henao A, Harton M, Basaraba RJ, Henao-Tamayo M, Barrozo JC, Rose J, Kawamura LM, Coscolla M, Fofanov VY, Koshinsky H, Gagneux S, Hopewell PC, Ordway DJ, Orme IM: Beijing sublineages of mycobacterium tuberculosis differ in pathogenicity in the guinea Pig. Clin Vaccine Immunol. 2012, 19: 1227-1237. 10.1128/CVI.00250-12.PubMed CentralView ArticlePubMedGoogle Scholar
- Reiling N, Homolka S, Walter K, Brandenburg J, Niwinski L, Ernst M, Herzmann C, Lange C, Diel R, Ehlers S, Niemann S: Clade-specific virulence patterns of Mycobacterium tuberculosis complex strains in human primary macrophages and aerogenically infected mice. MBio. 2013, 4: e00250-13. 10.1128/mBio.00250-13.PubMed CentralView ArticlePubMedGoogle Scholar
- Comas I, Coscolla M, Luo T, Borrell S, Holt KE, Kato-Maeda M, Parkhill J, Malla B, Berg S, Thwaites G, Yeboah-Manu D, Bothamley G, Mei J, Wei L, Bentley S, Harris SR, Niemann S, Diel R, Aseffa A, Gao Q, Young D, Gagneux S: Out-of-Africa migration and Neolithic coexpansion of Mycobacterium tuberculosis with modern humans. Nat Genet. 2013, 45: 1176-1182. 10.1038/ng.2744.PubMed CentralView ArticlePubMedGoogle Scholar
- Gagneux S: Host-pathogen coevolution in human tuberculosis. Philos Trans R Soc Lond B Biol Sci. 2012, 367: 850-859. 10.1098/rstb.2011.0316.PubMed CentralView ArticlePubMedGoogle Scholar
- World Health Organization: Global Tuberculosis Report 2013. 2013. apps.who.int/iris/ bitstream/10665 /91355/1/9789241564656_eng.pdf.Google Scholar
- WHO: Papua New Guinea Tuberculosis Country Profile. WHO; 2013. ., [http://www.who.int/tb/country/en/]
- Lumb R, Bastion I, Carter R, Jelfs P, Keehner T, Sievers A: Tuberculosis in Australia: bacteriologically confirmed cases and drug resistance, 2008 and 2009. A report of the Australian Mycobacterium Reference Laboratory Network. Commun Dis Intell. 2011, 35: 154-161.Google Scholar
- Gilpin CM, Simpson G, Vincent S, O’Brien TP, Knight TA, Globan M, Coulter C, Konstantinos A: Evidence of primary transmission of multidrug-resistant tuberculosis in the Western Province of Papua New Guinea. Med J Aust. 2008, 188: 148-152.PubMedGoogle Scholar
- Simpson G, Coulter C, Weston J, Knight T, Carter R, Vincent S, Robertus L, Konstantinos A: Resistance patterns of multidrug-resistant tuberculosis in Western Province, Papua New Guinea. Int J Tuberc Lung Dis. 2011, 15: 551-552. 10.5588/ijtld.10.0347.View ArticlePubMedGoogle Scholar
- McBryde M: Evalutation of Risks of Tuberculosis in Western Province Papua New Guinea. Australian Aid; 2012. ., [http://aid.dfat.gov.au/countries/pacific/png/Documents/png-tb-evaluation-of-risk.pdf]
- Ballif M, Harino P, Ley S, Coscolla M, Niemann S, Carter R, Borrell S, Siba P, Phuanukoonnon S, Gagneux S, Beck HP:Drug resistance-conferring mutations in Mycobacterium tuberculosis from Madang. Papua New Guinea. BMC Microbiol. 2012, 12: 191-10.1186/1471-2180-12-191.PubMed CentralView ArticlePubMedGoogle Scholar
- Cross GB, Coles K, Nikpour M, Moore OA, Denholm J, McBryde ES, Eisen DP, Warigi B, Carter R, Pandey S, Harino P, Siba P, Coulter C, Mueller I, Phuanukoonnon S, Pellegrini M: TB incidence and characteristics in the remote gulf province of Papua New Guinea: a prospective study. BMC Infect Dis. 2014, 14: 93-10.1186/1471-2334-14-93.PubMed CentralView ArticlePubMedGoogle Scholar
- Main P, Attenborough R, Chelvanayagam G, Bhatia K, Gao X: The peopling of New Guinea: evidence from class I human leukocyte antigen. Hum Biol. 2001, 73: 365-383. 10.1353/hub.2001.0036.View ArticlePubMedGoogle Scholar
- Yoshida M, Ohtsuka R, Nakazawa M, Juji T, Tokunaga K: HLA-DRB1 frequencies of non-Austronesian-speaking Gidra in south New Guinea and their genetic affinities with Oceanian populations. Am J Phys Anthropol. 1995, 96: 177-181. 10.1002/ajpa.1330960206.View ArticlePubMedGoogle Scholar
- Laehy M: The central highlands of New Guinea.Geogr J, 1936, 87.,Google Scholar
- Brosch R, Gordon SV, Marmiesse M, Brodin P, Buchrieser C, Eiglmeier K, Garnier T, Gutierrez C, Hewinson G, Kremer K, Parsons LM, Pym AS, Samper S, van Soolingen D, Cole ST: A new evolutionary scenario for the Mycobacterium tuberculosis complex. Proc Natl Acad Sci U S A. 2002, 99: 3684-3689. 10.1073/pnas.052548299.PubMed CentralView ArticlePubMedGoogle Scholar
- Portevin D, Gagneux S, Comas I, Young D: Human macrophage responses to clinical isolates from the Mycobacterium tuberculosis complex discriminate between ancient and modern lineages. PLoS Pathog. 2011, 7: e1001307-10.1371/journal.ppat.1001307.PubMed CentralView ArticlePubMedGoogle Scholar
- Parwati I, Van CR, Sudiro M, Alisjahbana B, Pakasi T, Kremer K, van der Zanden A, van Soolingen D: Mycobacterium tuberculosis population structures differ significantly on two Indonesian Islands. J Clin Microbiol. 2008, 46: 3639-3645. 10.1128/JCM.00605-08.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen YY, Chang JR, Huang WF, Kuo SC, Yeh JJ, Lee JJ, Jang CS, Sun JR, Chiueh TS, Su IJ, Dou HY: Molecular epidemiology of Mycobacterium tuberculosis in aboriginal peoples of Taiwan, 2006–2011. J Infect. 2014, 68: 332-337. 10.1016/j.jinf.2013.12.004.View ArticlePubMedGoogle Scholar
- Ballif M, Harino P, Ley S, Carter R, Coulter C, Niemann S, Borrell S, Fenner L, Siba P, Phuanukoonnon S, Gagneux S, Beck HP: Genetic diversity of Mycobacterium tuberculosis in Madang, Papua New Guinea. Int J Tuberc Lung Dis. 2012, 16: 1100-1107. 10.5588/ijtld.11.0779.View ArticlePubMedGoogle Scholar
- Hiatt T, Nishikiori N: Epidemiology and control of tuberculosis in the Western Pacific Region: analysis of 2012 case notification data. Western Pac Surveill Response J. 2014, 5: 25-34. 10.5365/wpsar.2014.5.1.013.PubMed CentralView ArticlePubMedGoogle Scholar
- Stucki D, Malla B, Hostettler S, Huna T, Feldmann J, Yeboah-Manu D, Borrell S, Fenner L, Comas I, Coscollà M, Gagneux S: Two new rapid SNP-typing methods for classifying Mycobacterium tuberculosis complex into the main phylogenetic lineages. PLoS One. 2012, 7: e41253-10.1371/journal.pone.0041253.PubMed CentralView ArticlePubMedGoogle Scholar
- Ley SD, Riley I, Beck HP: Tuberculosis in Papua New Guinea: from yesterday until today. Microbes Infect. 2014, 16: 607-614. 10.1016/j.micinf.2014.06.012.View ArticlePubMedGoogle Scholar
- Yen S, Bower JE, Freeman JT, Basu I, O’Toole RF: Phylogenetic lineages of tuberculosis isolates in New Zealand and their association with patient demographics. Int J Tuberc Lung Dis. 2013, 17: 892-897. 10.5588/ijtld.12.0795.View ArticlePubMedGoogle Scholar
- Albanna AS, Reed MB, Kotar KV, Fallow A, McIntosh FA, Behr MA, Menzies D: Reduced transmissibility of East African Indian strains of Mycobacterium tuberculosis. PLoS One. 2011, 6: e25075-10.1371/journal.pone.0025075.PubMed CentralView ArticlePubMedGoogle Scholar
- Steiner A, Stucki D, Coscolla M, Borrell S, Gagneux S: KvarQ: targeted and direct variant calling fastq reads of bacterial genomes. BMC Genomics. 2014, 15: 881-10.1186/1471-2164-15-881.PubMed CentralView ArticlePubMedGoogle Scholar
- Tsolaki AG, Gagneux S, Pym AS, de la Salmoniere YO G, Kreiswirth BN, Van Soolingen D, Small PM: Genomic deletions classify the Beijing/W strains as a distinct genetic lineage of Mycobacterium tuberculosis. J Clin Microbiol. 2005, 43: 3185-3191. 10.1128/JCM.43.7.3185-3191.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- van Soolingen D, Qian L, de Haas PE, Douglas JT, Traore H, Portaels F, Qing HZ, Enkhsaikan D, Nymadawa P, van Embden JD: Predominance of a single genotype of Mycobacterium tuberculosis in countries of east Asia. J Clin Microbiol. 1995, 33: 3234-3238.PubMed CentralPubMedGoogle Scholar
- Fenner L, Malla B, Ninet B, Dubuis O, Stucki D, Borrell S, Huna T, Bodmer T, Egger M, Gagneux S: “Pseudo-Beijing”: evidence for convergent evolution in the direct repeat region of Mycobacterium tuberculosis. PLoS One. 2011, 6: e24737-10.1371/journal.pone.0024737.PubMed CentralView ArticlePubMedGoogle Scholar
- Hanekom M, van der Spuy GD, Streicher E, Ndabambi SL, McEvoy CR, Kidd M, Beyers N, Victor TC, van Helden PD, Warren RM: A recently evolved sublineage of the Mycobacterium tuberculosis Beijing strain family is associated with an increased ability to spread and cause disease. J Clin Microbiol. 2007, 45: 1483-1490. 10.1128/JCM.02191-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Hanekom M, van der Spuy GD, van Pittius NCG, McEvoy CR, Ndabambi SL, Victor TC, Hoal EG, van Helden PD, Warren RM: Evidence that the spread of Mycobacterium tuberculosis strains with the Beijing genotype is human population dependent. J Clin Microbiol. 2007, 45: 2263-2266. 10.1128/JCM.02354-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Cowley D, Govender D, February B, Wolfe M, Steyn L, Evans J, Wilkinson RJ, Nicol MP: Recent and rapid emergence of W-Beijing strains of Mycobacterium tuberculosis in Cape Town, South Africa. Clin Infect Dis. 2008, 47: 1252-1259. 10.1086/592575.View ArticlePubMedGoogle Scholar
- Brites D, Gagneux S: Old and new selective pressures on Mycobacterium tuberculosis. Infect Genet Evol. 2012, 12: 678-685. 10.1016/j.meegid.2011.08.010.PubMed CentralView ArticlePubMedGoogle Scholar
- Caws M, Thwaites G, Stepniewska K, Nguyen TN, Nguyen TH, Nguyen TP, Mai NT, Phan MD, Tran HL, Tran TH, van Soolingen D, Kremer K, Nguyen VV, Nguyen TC, Farrar J: Beijing genotype of Mycobacterium tuberculosis is significantly associated with human immunodeficiency virus infection and multidrug resistance in cases of tuberculous meningitis. J Clin Microbiol. 2006, 44: 3934-3939. 10.1128/JCM.01181-06.PubMed CentralView ArticlePubMedGoogle Scholar
- National Department of Health: Papua New Guinea National Health Plan 2011–2020. 2010. ., [http://www.adi.org.au/wp-content/uploads/2013/07/2011_AidsCouncilPNG_FACT-SHEET.pdf]
- Assam JP, Beng VP, Cho-Ngwa F, Toukam M, Ngoh AA, Kitavi M, Nzuki I, Nyonka JN, Tata E, Tedom JC, Skilton RA, Pelle R, Titanji VP: Mycobacterium tuberculosis is the causative agent of tuberculosis in the southern ecological zones of Cameroon, as shown by genetic analysis. BMC Infect Dis. 2013, 13: 431-10.1186/1471-2334-13-431.PubMed CentralView ArticleGoogle Scholar
- Maddocks I, Anders EM, Dennis E: Donovanosis in Papua New Guinea. Br J Vener Dis. 1976, 52: 190-196.PubMed CentralPubMedGoogle Scholar
- Okamoto S, Tamaru A, Nakajima C, Nishimura K, Tanaka Y, Tokuyama S, Suzuki Y, Ochi K: Loss of a conserved 7-methylguanosine modification in 16S rRNA confers low-level streptomycin resistance in bacteria. Mol Microbiol. 2007, 63: 1096-1106. 10.1111/j.1365-2958.2006.05585.x.View ArticlePubMedGoogle Scholar
- Zhang Y, Yew WW: Mechanisms of drug resistance in Mycobacterium tuberculosis. Int J Tuberc Lung Dis. 2009, 13: 1320-1330.PubMedGoogle Scholar
- Spies FS, da Silva PE, Ribeiro MO, Rossetti ML, Zaha A: Identification of mutations related to streptomycin resistance in clinical isolates of Mycobacterium tuberculosis and possible involvement of efflux mechanism. Antimicrob Agents Chemother. 2008, 52: 2947-2949. 10.1128/AAC.01570-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Ramaswamy S, Musser JM: Molecular genetic basis of antimicrobial agent resistance in Mycobacterium tuberculosis: 1998 update. Tuber Lung Dis. 1998, 79: 3-29. 10.1054/tuld.1998.0002.View ArticlePubMedGoogle Scholar
- National Department of Health, Disease Control Program, National Tuberculosis Program: Papua New Guinea National Tuberculosis Management Protocol. 2011:1–38. ., [http://www.adi.org.au/health-in-png-2/health-reports/#tb]
- Hillemann D, Rusch-Gerdes S, Richter E: Evaluation of the GenoType MTBDRplus assay for rifampin and isoniazid susceptibility testing of Mycobacterium tuberculosis strains and clinical specimens. J Clin Microbiol. 2007, 45: 2635-2640. 10.1128/JCM.00521-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Plinke C, Cox HS, Zarkua N, Karimovich HA, Braker K, Diel R, Rüsch-Gerdes S, Feuerriegel S, Niemann S: embCAB sequence variation among ethambutol-resistant Mycobacterium tuberculosis isolates without embB306 mutation. J Antimicrob Chemother. 2010, 65: 1359-1367. 10.1093/jac/dkq120.View ArticlePubMedGoogle Scholar
- Safi H, Lingaraju S, Amin A, Kim S, Jones M, Holmes M, McNeil M, Peterson SN, Chatterjee D, Fleischmann R, Alland D: Evolution of high-level ethambutol-resistant tuberculosis through interacting mutations in decaprenylphosphoryl-beta-D-arabinose biosynthetic and utilization pathway genes. Nat Genet. 2013, 45: 1190-1197. 10.1038/ng.2743.View ArticlePubMedGoogle Scholar
- Petroff SA: A new and rapid Method for the Isolatin and Cultivation of Tubercle Bacilli directly from the Sputum and Feces. J Exp Med. 1915, 21: 38-42. 10.1084/jem.21.1.38.PubMed CentralView ArticlePubMedGoogle Scholar
- Kamerbeek J, Schouls L, Kolk A, van Agterveld M, van Soolingen D, Kuijper S, Bunschoten A, Molhuizen H, Shaw R, Goyal M, van Embden J: Simultaneous detection and strain differentiation of mycobacterium tuberculosis for diagnosis and epidemiology. J Clin Microbiol. 1997, 35: 907-914.PubMed CentralPubMedGoogle Scholar
- Demay C, Liens B, Burguière T, Hill V, Couvin D, Millet J, Mokrousov I, Sola C, Zozio T, Rastogi N: SITVITWEB – a publicly available international multimarker database for studying Mycobacterium tuberculosis genetic diversity and molecular epidemiology. Infect Genet Evol. 2012, 12: 755-766. 10.1016/j.meegid.2012.02.004.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.