Skip to main content

Advertisement

Diversity of Mycobacterium tuberculosis and drug resistance in different provinces of Papua New Guinea

Article metrics

Abstract

Background

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.

Results

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%).

Conclusion

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.

Background

Over the last decade, evidence for the impact of the bacterial genetic background on TB infection and disease has strongly increased [1]. 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 [2],[3]. The different M. tuberculosis strains were grouped into 6 main lineages (a 7th was recently added [4]). 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 [5] and drug resistance [6],[7]. Lineage 4 on the other hand, has been associated with pulmonary TB rather than extrapulmonary TB [8]. In addition, few studies have analysed the impact of the genetic background of both, humans and bacteria on disease development [8],[9] 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 [10],[11]. 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 [2],[12],[13]. 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% [14],[15]. Publications on drug resistant TB in PNG are rare, and previous data were mainly derived from patients from Western Province, diagnosed in Australia [16]-[19]. Gilpin et al. reported 25% MDR-TB among patients from Western Province diagnosed between 2000 and 2006 [17]. 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 [20], 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 [21].

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 [22],[23]. The highlands region has only been ‘discovered’ in the 1930s [24] and was sparsely populated at the time. Hence it would be expected that evolutionary ‘ancient’ lineages of M. tuberculosis (e.g. Lineage 1 [25]), proposed to be adapted to low density populations [26], 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 [27] or Taiwan [28]) 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 [29], another study where isolates from patients from Western Province were analysed [17] and a very recent study from Gulf Province [21], 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

A total of 449 patients enrolled into our study were diagnosed with TB of any type. From 396 of these patients sputum samples could be collected and of these 335 (74.6%) samples were available for study purposes (see Figure 1). From 212 (63.3%) of the available samples M. tuberculosis were successfully grown in culture and drug susceptibility testing (DST) was performed. Details of population characteristics per study site are described in Table 1.

Figure 1
figure1

Flow chart of cohort sample size development from screening to genotyping.

Table 1 Characteristics of study population for each study site

TB positivity peaked in the age group of 15 to 34 years, which is consistent with country wide data [30]. 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.

Genotyping

Genotyping of MTBC could successfully be conducted of 147/212 (69%) samples. With a TaqMan Real-time PCR assay [31] three of seven worldwide reported MTBC lineages [3],[12], were detected in all three study sites. Overall, Lineage 4 was the most prevalent lineage with 75/147 isolates (51.0%), followed by Lineage 2 with 67/147 isolates (45.6%), whilst Lineage 1 was rare with 5/147 isolates (3.4%). This composition of circulating lineages was similar as previously described for Madang [29] and other parts of PNG [17],[21],[32] and reflects the reduced MTBC diversity within PNG compared to other countries of the region, e.g. Indonesia [27] or New Zealand [33] where Lineage 3 strains were also found. When lineage compositions (shown in Figure 2) were compared between the study sites, statistically significant differences were observed between all three sites (Fisher’s exact test (p < 0.001)). Lineage 1 was generally rare in all three sites. Lineage 4 was the most abundant lineage in Madang (76.6%) similar to what has previously been described [29]. In contrast, in Alotau Lineage 2 was the dominant lineage (84.4%). In Goroka, a trend towards a higher prevalence of Lineage 2 (60.5%) was found but was not as high as in Alotau. However, it cannot be ruled out that the low proportion of Lineage 1 in our cohort is due to the strain genetic background itself, e.g. exhibiting reduced transmissibility [34] and a lower growth rate in macrophages compared to Lineage 2 and Lineage 4 [11] potentially leading to a decreased culture recovery of Lineage 1.

Figure 2
figure2

Mycobacterium tuberculosis lineage composition for each study site. Black bars = lineage 1 (Indo-Oceanic); dashed bars = lineage 2 (East-Asian); white bars = lineage 4 (Euro-American).

We further performed spoligotyping, and the spoligotyping pattern with the corresponding families and frequencies are shown in Figure 3. Thirty different spoligotyping patterns belonged to nine different families, including 14 orphans with 12 different spoligotyping patterns. Orphans were strains with no matching entry in the SITVIT database and are therefore considered as undefined. Whether these strains represent PNG specific M. tuberculosis strains remains to be confirmed using other methods.

Figure 3
figure3

Spoligotypes for the observed lineages (n = 145). SIT = shared international type.

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 [35]. Per standard definition, this would define that strain as belonging to the Beijing family of Lineage 2 [36],[37]. 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 [38].

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) [36]. 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 [39]-[41]. 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 [22], however, no statistically significant difference in the prevalence of Lineage 1 in Goroka could be detected.

Other factors might also influence the lineage distribution. Goroka and Madang are connected to the country through the highlands highway, whilst Madang and Alotau have direct access to the island of PNG and beyond. Similarly, the host genetic background and environmental or circumstantial factors such as co-infections or age could well influence the M. tuberculosis lineage distribution. We therefore performed univariate and multivariate logistic regressions to test whether infection with a Lineage 2 strain (equals an infection with Beijing strain) was associated with other factors such as DR, HIV status, age and gender (Table 2). Samples not belonging to Lineage 2 were pooled for the analysis, consisting of five samples of Lineage 1 and 75 samples of Lineage 4. Univariate logistic regression confirmed the differences in lineage distribution between the study sites. These differences remained significant also when correcting for possible confounders in a multivariate regression.

Table 2 Univariate and multivariate logistic regressions for the risk of an infection with a Beijing strain of lineage 2

Co-infection with HIV has several times been shown to be associated with infection with strains of the Beijing type [42],[43]. 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) [44] 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 [3]. Thus, our findings could support the notion of Hershberg and colleagues [2] 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 [45] and might have started in PNG.

Drug resistance

In total 23/212 (10.8%) samples were resistant to at least one of the drugs tested. The details of the DR patterns and their frequency in each study site are shown in Table 3. With 6.6% (14/212) the proportion of isolates resistant (mono-resistant or poly-resistant) to Streptomycin (STR) was the highest in our cohort. High STR resistance in PNG is well known and has been found in previous studies [18]-[20]. In the past STR has been used frequently as single drug to treat TB but also urinary tract and Klebsiella infections [46] which might have led to drug resistance. Drug resistance of any type and MDR frequencies were highest in Alotau with 4.6% of MDR cases, but monoresistance was observed more often in Goroka (8.9%) and Madang (6.7%) compared to Alotau (4.6%). However, none of these differences were statistically significant (p = 0.960).

Table 3 Observed phenotypic drug resistance per study site

In order to determine drug resistance mechanisms of 16 phenotypically drug resistant isolates for which DNA could be obtained (16/23), we determined mutations in ten genes known to be associated with resistance. Two pan-susceptible strains of patients with late sputum conversion were also included. Sequence data were analysed and mutations observed are shown in Figure 4.

Figure 4
figure4

Overview of genes/gene regions sequenced and mutations observed. No mutations were found in ahpC, pncA and gyrA. Black squares: mutation detected at indicated position; white squares: no mutation found at indicated position; dotted squares: no mutation found at any of the indicated positions, which are the most commonly known sites of mutation for the detected phenotypic drug resistance patterns. All positions indicate the amino acid change at the codon position, except for rrs gene where nucleotide position and change is indicated.

Streptomycin resistance has been shown to be associated with mutations in the rrs, rpsL or gidB gene [47],[48]. 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 [20],[49] 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 [20] 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) [50]. 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) [51] which determines RMP resistance only through detection of mutations in RRDR [52]. 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 [48]. 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 [53],[54]. 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[50] 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.

Possible associations with known risk factors for DR were tested using univariate logistic regressions (Table 4). Risk of being infected with a DR strain was 5.5 times higher for patients with a history of TB treatment and 2.4 times higher for patients infected with an Lineage 2 strain although the latter was not statistically significant (p = 0.097). To correct for possible confounders multivariate logistic regression analysis was conducted and the adjusted odds ratio (OR) for risk of being infected with a DR strain with a history of TB treatment decreased from 5.5 to 4.2 but remained significant (p = 0.040). This association highlights the importance of constant access to treatment and compliance for the control of DR TB. Strengthening the DOTS strategy and increasing awareness of TB in the population is crucial also in PNG and needs to be maintained at a high level in order not to delay diagnosis and to prevent the possible spread of TB.

Table 4 Univariate and multivariate logistic regressions for the risk of drug resistance

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) [29]. 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).

Conclusions

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.

Methods

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 [51]. 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 [55]. 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 [29].

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 [12] by a TaqMan real-time PCR assay using single nucleotide polymorphism (SNP) typing [31]. Lineages were further discriminated into families by spoligotyping [56], 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 [57]. 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. [36]. PCR conditions and primers were used as previously described by Gagneux et al.[3].

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.[20] 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

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

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.

References

  1. 1.

    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.

  2. 2.

    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.

  3. 3.

    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.

  4. 4.

    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.

  5. 5.

    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.

  6. 6.

    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.

  7. 7.

    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.

  8. 8.

    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.

  9. 9.

    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.

  10. 10.

    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.

  11. 11.

    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.

  12. 12.

    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.

  13. 13.

    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.

  14. 14.

    World Health Organization: Global Tuberculosis Report 2013. 2013. apps.who.int/iris/ bitstream/10665 /91355/1/9789241564656_eng.pdf.

  15. 15.

    WHO: Papua New Guinea Tuberculosis Country Profile. WHO; 2013. ., [http://www.who.int/tb/country/en/]

  16. 16.

    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.

  17. 17.

    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.

  18. 18.

    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.

  19. 19.

    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]

  20. 20.

    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.

  21. 21.

    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.

  22. 22.

    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.

  23. 23.

    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.

  24. 24.

    Laehy M: The central highlands of New Guinea.Geogr J, 1936, 87.,

  25. 25.

    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.

  26. 26.

    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.

  27. 27.

    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.

  28. 28.

    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.

  29. 29.

    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.

  30. 30.

    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.

  31. 31.

    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.

  32. 32.

    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.

  33. 33.

    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.

  34. 34.

    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.

  35. 35.

    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.

  36. 36.

    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.

  37. 37.

    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.

  38. 38.

    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.

  39. 39.

    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.

  40. 40.

    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.

  41. 41.

    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.

  42. 42.

    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.

  43. 43.

    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.

  44. 44.

    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]

  45. 45.

    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.

  46. 46.

    Maddocks I, Anders EM, Dennis E: Donovanosis in Papua New Guinea. Br J Vener Dis. 1976, 52: 190-196.

  47. 47.

    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.

  48. 48.

    Zhang Y, Yew WW: Mechanisms of drug resistance in Mycobacterium tuberculosis. Int J Tuberc Lung Dis. 2009, 13: 1320-1330.

  49. 49.

    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.

  50. 50.

    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.

  51. 51.

    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]

  52. 52.

    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.

  53. 53.

    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.

  54. 54.

    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.

  55. 55.

    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.

  56. 56.

    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.

  57. 57.

    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.

Download references

Acknowledgements

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.

Author information

Correspondence to Hans-Peter Beck.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contribution

SDL: contributed to the study design, implemented and coordinated the study, participated in the sample collection and processing, carried out the molecular analysis, conducted data analysis and drafted the manuscript. PH: contributed to study design and coordinated the study in one of the study sites and was involved in sample collection and responsible for patient management. KV and RK were responsible for patient management. CC, SP and RC supervised culturing and phenotypic DST and helped to draft the paper. JF carried out culturing at Swiss TPH and was involved in the molecular analysis. MB contributed to the study design and was involved in data analysis. SPM: was involved in conceiving the study, contributed to its design and helped to draft the manuscript. PS: contributed to the study design, SG: participated in the design of the study, coordinated the molecular work and helped to draft the paper. HPB: conceived the study and contributed to its design and coordination and acquired the funding for the study and helped to draft the manuscript. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Keywords

  • Tuberculosis
  • Papua New Guinea
  • Genotyping