Skip to main content

Phenotypic and genotypic drug sensitivity profiles of Mycobacterium tuberculosis infection and associated factors in northeastern Ethiopia

Abstract

Background

Tuberculosis is a devastating and a deadly disease despite the novel advances in its diagnostic tools and drug therapy. Drug resistant Mycobacterium contributes a great share to tuberculosis mortality. Status of drug resistance and patients’ awareness toward the disease is unknown in northeastern Ethiopia. Thus, the aim of this study was to determine the phenotypic and genotypic drug sensitivity patterns and associated factors in Oromia Special Zone and Dessie Town, northeastern Ethiopia.

Methods

In a cross-sectional study, 384 smear positive tuberculosis cases were recruited and Löwenstein-Jensen culture was done. The performance of GenoTypic MTBDRplus assay using the conventional BACTEC MGIT 960 as a “gold standard” was determined. Drug resistant strains were identified using spoligotyping. Pearson Chi-square test was used to determine the association of drug sensitivity test and tuberculosis type, lineages, dominant strains and clustering of the isolates.

Results

The 384 smear positive Mycobacterium samples were cultured on LJ media of which 29.2% (112/384) as culture positive. A fair agreement was found between MTBDRplus assay and the conventional MGIT test in detecting the Mycobacterium tuberculosis with sensitivity, specificity, positive and negative predictive value of 94.2, 30.2, 68.4 and 76.5%, respectively. Among LJ culture positive samples 95 of them gave valid result for MTBDRplus assay and 16.8% (16/95) as drug resistant. Similarly, MGIT subculture was made for the 112 isolates and 69 of them gave positive result with 15.9% (11/69) as drug resistant. Cohen’s kappa value showed almost a perfect agreement between the two testing methods in detecting rifampicin (sensitivity 100% and specificity 98.3%) and multi-drug resistance (sensitivity 83.3% and specificity 100%). Spoligotyping identified 76.5% (13/17) of the drug resistant isolates as Euro-American and family 33 as the predominant family. Significant association was observed between drug resistant isolates and the dominant strains (χ2: 34.861; p = 0.040) of the Mycobacterium.

Conclusion

Higher magnitude of drug resistance was found in the study area. The GenoTypic MDRTBplus assay had an acceptable drug sensitivity testing performance.

Peer Review reports

Background

Ethiopia is one of the 30 top global tuberculosis (TB) burden countries. Antimicrobial susceptibility testing is vital in prescribing an effective drug regime for TB patient, especially in areas where drug resistance incidence is high like Ethiopia. Although the global mortality rate of TB is reducing by 3% annually, the threat of its drug resistance is increasing [1]. The cumulative effects of treatment interruption like lack of awareness about the nature of bacteria, shortage and lack of WHOs recommended diagnostic tools, and prolonged drug consumption period for treatment increases the risk. Unless the new drug regimens are used, multidrug resistant (MDR) TB treatment could take more time causing longer absenteeism from work or even loss of employment, social isolation, and long-term socioeconomic and psychological effects [2].

In 2018, an estimated 3.4% of the global TB cases were new drug resistant (MDR/RR) and 18% were among previously treated cases. In Ethiopia, the estimated incidence of MDR/RR is 0.71 and 16% among new and previously treated cases, respectively [3]. In addition, the estimated cases of MDR-RR had a magnitude of 484,000 which was about 10% downward from the best estimate published by WHO in its 2018 global TB report. Of these estimated cases, about 44.2% (214,000) deaths were due to MDR/RR-TB which was also a downward revision to the best estimates. In fact, the global notified cases rather than estimates of MDR/RR-TB were 186,772 up from 160,772 cases in 2017, and 156,071 cases were enrolled in treatment which was also up from 139,114 in 2017. The number of people enrolled in treatment in the year was equivalent to only 32% of the estimated incidence of the 484,000 cases [3, 4].

The Genotype MTBDRplus assay diagnostic tools were used and correctly identified Rifampicin (RIF), Isoniazid (INH) and MDR-TB resistant isolates with best efficacy of sensitivity and specificity as shown in a study done at Cameroon [5]. The “End TB Strategy” calls early diagnosis of TB and prompt treatment for all persons of all ages with any form of drug susceptible or drug resistant TB. WHO further defines a universal access to those drug sensitivity testing technologies of both LPAs and Mycobacterium growth indicator tubes (MGIT). For proper treatment of the drug resistant tuberculosis on time a lesser time taking testing methods (GeneXpert, LPAs and MGITs) are recommended than the gold standard LJ culture method [6].

Spoligotyping is a widely used DNA fingerprinting methods that allows simultaneous detection and molecular typing of M. tuberculosis complex. It has a meticulous value in population based studies with low cost to define the phylogeographic specificity of the circulating clades/families of tubercle bacilli. It is the DR region used in individual M. tuberculosis strains and in different members of the M. tuberculosis complex to identify and align the spoligotype patterns. This imply that spoligotyping is recommended as the best preliminary screening test for M. tb isolates that circulates in the societies [7].

In Ethiopia, although great expansion of health institutes are built and primary attention is given to TB control and management, drug sensitivity testing is limited. Thus, the aim of this study was to assess the performance of both phenotypic and genotypic drug sensitivity test, status of drug resistance and the circulating resistant strains in Oromia Special Zone and Dessie Town administration in northeastern Ethiopia.

Methods

The study area and sample collection sites

The study was conducted in Oromia Special Zone (OSZ) and South Wollo Zone (SWZ) of the Amhara Regional State, northeastern Ethiopia. Kemise and Bati Town health centers were the data collection sites from OSZ where as Dessie Referral Hospital (DRH), Bikat Higher Diagnostic Laboratories (BHDL), Dessie Health Center (DHC) and Boru Meda Hospital (BMH) were the sample collection sites from South Wollo. Tuberculosis cases confirmed by the health personnel who fulfilled the inclusion criteria were included in the study. The samples (sputum and Fine needle aspirates (FNAs)) were collected on the spot from consenting participants 18 years and older of all TB cases as of April 2015 to January 2017. Those with severe TB who were unable to provide their sputum specimens were excluded from the study.

Study design and laboratory processing

An institution-based cross-sectional study design was used among all forms of TB cases. Dry, translucent, leak-proof 50 ml capacities of falcon tubes were used to collect a minimum of 3-5 ml sputum sample and labeled using indelible labeling markers. The samples were collected from all those smear-positive TB confirmed participants of both previously untreated (new) and previously treated (retreatment) cases.

For bacteriological examination, direct microscopic examination was used after Ziehl-Neelsen staining technique performance at sample collection sites. A portion of Positive samples were kept at a range of − 10 to − 20 °C in the refrigerator of the health institutes. The FNAs from patients suspected for extra-pulmonary tuberculosis (EPTB) was collected using a 21-gauge needle attached to a 10 ml syringe with maximum care and safety by an experienced pathologist [8]. Afterwards, the Ziehl-Neelsen smear technique was performed by the same pathologist to check its positivity. A suction of about 1 ml samples was collected from positive patients for this research purpose and preserved in sterile and tightly closed nunc tubes with phosphate buffer saline of equal amount at pH 7.2 and kept in the same refrigerator [9]. Finally, both sputum and FNA samples were transported using a cold chain of 4 °C to TB laboratory of Aklilu Lemma Institute of Pathobiology (ALIPB) in Addis Ababa where the samples were kept in a deep freezer of – 80 °C until culturing was done [10].

A 37.2 g stock of selective Lowenstein-Jensen (LJ) media was suspended in 600 ml distilled water supplemented with 12 ml of glycerol. The glycerol was replaced by 10 g of sodium pyruvate for the second culturing media. Then frequent agitation was made for a minute by boiling. The final solution was cooled and sterilized in autoclave at 121 °C for 15 min. For the preparation of homogenized whole eggs, fresh hen eggs were washed with running water and sterilized by 70% ethanol before breaking to mix by shaking under aseptic condition. One liter of this homogenized egg was added to the prepared autoclaved media solution, shacked and sieved with sterile sieve. The completely prepared medium of 6-8 ml volume was dispensed into screw-capped culture tubes and their lids were securely fastened. Finally, the tubes were slanted at 15–200 angles in inspissator and heated at 85 °C to solidify/coagulate for 45 min. The prepared media was left for 48 h before culturing the bacteria to test their contamination and those poor quality media were identified and discarded. Finally, the Mycobacterium culture was performed [11].

All the inoculated LJ slants in culture tubes were incubated aerobically at 37 °C in a slanted position and contamination was checked daily for the first week. The inoculated media was positioned in an upright position starting from the second week and colony formation of the bacteria was observed every week for eight consecutive weeks. The grown M. tb was heat killed and freezed whereas the weakly grown colonies were sub-cultured. Both the heat killed and freezed isolates were kept in a deep freezer of -80 °C until drug sensitivity test was done using both phenotypically and genotypically.

The molecular GeneXpert and line probe assay

A total of 95 sputum samples were subjected to the Xpert test and the test was performed according to the manufacturer’s instruction in DRH. The automated readout reports the presence of Mycobacterium with a detection of rifampicin (RIF) resistance and recorded [12]. The GenoType MTBDRplus version 2 LPA was used for identification of M. tuberculosis complex and its resistance to RIF, isoniazid (INH) or both using the heat killed DNA extracts. The assay was based on the DNA-STRIP technology and the whole procedure involves a multiplex PCR amplification with biotinylated primers, and a reverse hybridization. This assay detected for the absence and/or presence of wild type (WT) and/or mutant (MUT) DNA sequences with in specific region of three genes: the rpoB gene (coding for the β-subunit of the RNA polymeraze), for the identification of RIF resistance; the katG gene (coding for the catalase peroxidase), for high level INH resistance; and the promoter region of the inhA gene (coding for the NADH enoylACP reductase), for low-level INH resistance. The procedure of the test was performed based on the manufacturer’s instructions. Missing of bands in rpoB probes showed resistance to RIF whereas missing of katG and/or inhA indicates INH resistance. On the other hand, the mutation probes were also considered as resistant when the bands are as strong as or stronger than the existing AC bands. Absence of the wild type (WT) band is usually accompanied by the presence of mutant (MUT), which indicates resistance. In rare cases, missing of WT band(s) without a corresponding MUT band might be observed which was considered as due to “uncommon mutations” in the probe region. Presence of both WT and MUT bands in the same stripe might be an indication for the presence of hetero-resistance or mixed infection [13].

Phenotypic BD BACTEC™ MGIT™ 960 system (SIRE test)

Sub-culturing was performed for 112 LJ-culture positive freezed isolates using BACTEC MGIT 960 instrument (Becton Dickinson, Baltimore, MD, USA) for a maximum of 42 days from initial incubation date.. The drug sensitivity test to streptomycin (S), isoniazid (I), rifampicin (R) and ethambutol (E) was done for all of the successfully recovered sub-culture positive isolates following the protocols [14]. To control tubes, 0.5 ml of Growth Control working solution was added and labeled as drug free Mycobacterium growth indicator tubes (MGIT), and the other drugs containing tubes were labeled S, I, R and E. The result was interpreted as the growth unit (GU) of the growth control (GC) reached at a minimum value of 400. At the time the GU value of the GC was 400 or more and if the GU value of the drug tube was less than 100, the test result was reported as “susceptible,” while if the GU value of the drug tube was 100 or more the result was interpreted as “resistant”. The GU values of both the DST sets were retrieved and recorded. In case the GU value of the control did not reach 400 within 21 days, the instrument indicated an X200 error, indicating insufficient growth. On the other hand, if the GU reached 400 earlier than day 4, the instrument gave an X400 error, indicating contamination or over inoculation and in such cases the test was repeated [15]. In all cases of the study, the methods were carried out in accordance with relevant guidelines and regulations.

Data analysis

The recorded data was checked for completeness and consistency, and then entered into Microsoft Excel 2007 spreadsheets. The data was then exported to IBM SPSS Statistics for Windows, Version 25.0. (Armonk, NY: IBM Corp., USA) program for analysis. Descriptive statistics was used to determine frequency and percentage. In spoligotyping, the reference data base (SpolDB4) available online through https://www.pasteur-guadeloupe.fr:8081/SITVITDemo/ was used to assign the shared international spoligotype numbers (SIT) to known profiles for the drug resistant isolates and if not available the new patterns were considered as ‘orphan’ types. An online tool Run TB-Lineage with a website of http://tbinsight.cs.rpi.edu/run_tb_lineage.html was also used to identify family/clade, lineages and sublineages of the isolates. Those identical spoligotypes patterns with identical DNA genotypes were considered as a cluster and these clustered strains were identified as an indicator for the recent transmission.

Sensitivity, specificity, positive and negative predictive values were calculated to compare the performances between drug sensitivity tests. Agreement between the tests was assessed by Cohen’s Kappa statistics. The Kappa value was interpreted with values < 0 as indicating no agreement and 0–0.20 as slight, 0.21–0.40 as fair, 0.41–0.60 as moderate, 0.61–0.80 as substantial, and 0.81–1 as almost perfect agreement [16]. The Pearson Chi-square (χ2) was used to determine the association of drug sensitivity test and TB type, lineages, dominant strains and clustering of the isolates. The results were considered as statistically significant when the p-value was less than 5%.

Results

Socio-demographic characteristics and tuberculosis infections

A total of 384 TB cases took part in the study with majority of the participants being males (55.5%). Ages of the study participants ranged from 18 to 75 years and the median age was 30 (IQR = 15). The prevalence of TB cases was highest (67.0%) in the 18–37 years age group. There was no significant difference among males (mean = 34.5, SD = 12.3) and females (mean = 32.7, SD = 11.4).

GeneXpert MTB/RIF, GenoType MTBDRplus assay and BACTEC MGIT 960 test result

The overall LJ-culture positive sample was 29.2% (112/384). Among the collected samples, 95 were identified by rapid molecular diagnostic GeneXpert at DRH while the rest were detected by smear microscopic test (Fig. 1). Of the Xpert MTB/RIF tested samples, two were identified as RIF resistant. Twenty two of the 95 GeneXpert samples were LJ-culture positive and used for further drug sensitivity test using the GenoType MTBDRplus assay and the conventional BACTEC MGIT 960 System. Both methods confirmed the two RIF resistant samples by Xpert as drug resistant and in addition they also detect one more sample as multidrug resistant.

Fig. 1
figure1

A flowchart outlining sample collection and the number of resistant isolates using both phenotypic and genotypic drug sensitivity test

The sub-culture test by MGIT gave 61.6% (69/112) as positive and 64 of these sub-culture positive samples as from new TB cases and the rest five as from retreatment patients. Sixty five samples were correctly identified as M. tb and 30 of the MGIT sub-culture negative samples were detected as M. tb positive by the GenoType MTBDRplus assay. The remaining 13 of the 112 LJ-culture positive samples were negative to MGIT sub-culture and also not detected as M. tb by LPA. Of the 112 isolates tested by line probe assay 15.2% (17/112) have no TUB band. The translated sensitivity, specificity, positive predictive value and negative predictive value were 94.2, 30.2, 68.4 and 76.5%, respectively. There was a fair agreement (Kappa = 0.276; P < 0.001) between the two methods (MGIT and LPA) in detecting the Mycobacterium (Table 1) when MGIT was used as the reference standard in the analysis [17].

Table 1 M. tuberculosis detection rate by line probe assay and Mycobacterium growth indicator tube using the 112 LJ-culture positive samples

Sensitivity test done by LPA detected 16.8% (16/95) as drug resistant of which 18.8% (3/16) were from retreatment cases. The four first line anti-TB drugs (STM, INH, RIF and ETM) by conventional BACTEC MGIT 960 also identified comparable proportion 15.9% (11/69) as drug resistant with 8 and 3 samples as from the new and retreatment ones, respectively. More than half of the sub-culture positive samples by MIGIT 56.5% (39/69) were from Dessie referral hospital and all isolates from BMH were detected as drug resistant by either MGIT, LPA or both.

Drug resistance was detected in 8 common isolates of MGIT and LPA test. In addition, the LPA identified 2 isolates as susceptible which were resistant by MGIT and also no TUB band was observed in LPA strip for one isolate. Six of the resistant isolates detected by LPA were not tested for their sensitivity by MGIT due to their negative sub-culture result and 2 of the isolates resistant by LPA were susceptible to MGIT.

Greater proportion 8.4% (8/95) of the LPA tested isolates were multidrug resistance than mono-resistant. Of these, 5.3% (5/95) was among new patients and 3.2% (3/95) of the MDR as among retreatments. Rifampicin and isoniazid mono-resistance accounted for 6.3% (6/95) and 2.1% (2/95), respectively. The LPA test detects 57.1% (8/14) of rifampicin resistant isolates had multidrug resistant TB. On the other hand, all of the sub-culture positive samples by BACTEC MGIT 960 underwent 276 tests and 15.9% (11/69) of the samples were detected as resistant to any of the four anti-TB drugs. More resistance was detected 5.8% (4/69) to the four first-line antimicrobial drugs (INH + RIF + STM + EMB) by the MGIT test than their mono-resistance (INH =1.4%, RIF =1.4%, STM = 1.4%, EMB = 1.4%) (Fig. 2). The frequency of resistance to isoniazidniazd, rifampicin, streptomycin and ethambutol was in a comparable proportion 11.6% (8/69), 10.1% (7/69), 10.1% (7/69) and 8.7% (6/69); in that order. Greater proportion of resistance was detected in females 22.6% (7/31) than males 14.1% (9/64) by LPA test and also comparable proportion of drug resistance was detected in females 20% (4/20) than males 14.3% (7/49) by MGIT.

Fig. 2
figure2

Drug resistance pattern by genotypic line probe assay (LPA) and phenotypic Mycobacterium growth indicator tube (MGIT) from northeastern Ethiopia, April 2015 to January 2017. INH Isoniazid, RIF Rifampicin, STM Streptomycin, EMB Ethambutol, LPA line probe assay, MGIT Mycobacterium growth indicator tubes

As compared to BACTEC MGIT 960, the sensitivity and specificity of the GenoType MTBDRplus assay for the detection of RIF-resistant M. tb isolates was 100 and 98.3%, respectively. Similarly, the sensitivity and specificity for detection of INH-resistance was 75 and 98.2%, and for MDR resistance 83.3 and 100%, respectively. Kappa value showed that there is almost a perfect agreement between the two methods in detecting RIF and MDR where as there is a substantial agreement for the detection of isoniazid (Table 2).

Table 2 Performance of GenoTypic MTBDRplus assay for detection of RIF, INH and MDR resistance in comparison to phenotypic BACTEC MGIT 960 System

In comparison to MGIT, the drug sensitivity of LPA at katG gene was found as 75% detection and that of inhA promoter region was 12.5%. The sensitivity of LPA to both katG and inhA was also 12.5%. This finding also revealed that the katG, inhA, both katG and inhA results of LPA test had almost perfect agreement, of susceptibility to the MGIT test (Table 3).

Table 3 Detection of isoniazid (katG, inhA and both katG and inhA) resistance by LPA with its resistance by MGIT testing method

Resistance and mutation patterns of rifampicin and isoniazid using the GenoType MTBDRplus assay

For the line probe assay, single probe resistance was detected in 25.0% (4/16) of the isolates and two or more probe resistance was detected in the remaining drug resistant isolates. In addition, single probe mutation was detected in 12.5% (2/16) of the drug resistant isolates. In 10 of the RIF resistant isolates, the WT8 was missed with an additional miss of WT3, WT4 and WT6 probe in the rpoB gene of the three isolates and no mutation was detected. There were also an omission of the rpoBWT3 and WT7 probe for 3 RIF resistant isolates with unknown amino acid change for all the 13 RIF resistant isolates in the rpoB gene region. Concerning isonizid resistance, the resistance was further identified using katG and inhA gene region. Of the 12 INH-resistant isolates, katG resistance was detected in 66.7% (8/12) of the isolates by omitting WT probe at codon 315 and mutations were detected in two of the isolates at Codon S315T1. Mutation was identified in five isolates of the inhA gene region at WT1 (− 15/− 16) and an additional mutation to one isolate at WT2 (− 8). Addition of specific inhA mutations were not observed in place of the omitted probes at the inhA promoter region of the gene.

Genetic diversity for the drug resistant Mycobacterium tuberculosis

Of the major lineages identified using spoligotyping 76.5% (13/17) were Euro-American followed by Indo-Oceanic 17.6% (3/17) lineages. Most of the drug resistant isolates 76.5% (13/17) had no SIT number indicating that they are orphans. Family 33 is also the predominant family identified in this study (Table 4).

Table 4 Spoligotyping result of drug resistant isolates from northeastern Ethiopia

The association between any drug resistance and TB type showed a varied proportion but with no statistically significant difference among pulmonary and extra-pulmonary TB cases (χ2: 0.233; p = 0.629). A patient with previous TB treatment and drug-resistance linked to new cases infected with the same strain (Table 5).

Table 5 Association between any drug sensitivity patterns and TB type, major lineage, sub-lineage and dominant strains

Discussion

GenoTypic and phenotypic drug sensitivity test

The microbiological drug sensitivity test using GeneXpert is crucial to treat the RIF resistant (a surrogate marker for MDR-TB) patients on the spot than the other time taking testing methods. The additional detection of drug resistant isolate using MTBDRplus assay and BACTEC MGIT 960 test was due to single detection of resistance by the Xpert (only RIF resistant cases). The GenoTypic LPA can detect an additional INH resistance where as the phenotypic MGIT can detect even a more additional drugs such as STM and EMB to RIF and INH. This could enhance the chance of detecting more resistant isolates. The lower proportion of MGIT sub-culture positivity detection 61.6% (69/112) for DST test from LJ-culture positive Mycobacterium sample than the heat killed LPA isolates 84.8% (95/112) might be due to the existence of bacterial DNA in the later one. In the case of freezed live Mycobacterium, its survival could be interrupted through time until it was sub-cultured using the MGIT instrument. The higher sensitivity result 94.2% of LPA and its lower specificity value 30.2% is comparable with similar settings done in Kenya where its sensitivity and specificity was 99.2 and 26.9%, respectively [17]. This showed that LPA has a good performance in detecting the true positivity of Mycobacterium and in a fair agreement (Kappa = 0.276; P < 0.001) with the BACTEC MGIT 960 performance.

The overall genotypic drug sensitivity test in this study detected 16.8% (16/95) as drug resistant to either of RIF and/or INH. This is higher than 7.7% (10/130) of similar studies done by using LPA in Arsi South central Ethiopia [18], 8.7% (14/161) and 13.3% (37/279) reported from central Ethiopia [19]. In contrast, our finding is in less proportion than a study report from northwest Ethiopia [20] and southwest Ethiopia [21] having a rate of 18.0% (20/111) and 39.3% (44/112), respectively. Studies from Chad [22] and Punjab state of India [23] also reported greater proportion of resistance as 23.4% (73/312) and 58.4% (163/279), respectively. The possible explanation for such variations of drug resistant proportion could be due to difference in sample size, study sites and subjects where the samples were collected.

The rate of MDR in this study was 8.4% (8/95) by genotypic MTBDRplus assay which is greater than 1.1% (3/279) [19], 1.8% (2/111) [20], 3.1% (5/161) [24] and 6.7% (11/165 [25] similar studies. On the contrary, our MDR finding is less than, 17.5% (33/189) [26], 25.8% (72/279) [23] and 27.7% (31/112) [21]. Regarding the detection of RIF 6.3% (6/95) and INH 2.1% (2/95) mono-resistance, the current prevalence in line with a number of studies [27,28,29] that rifampicin mono-resistant is greater than isoniazid mono-resistant with varying degrees of proportion. Such RIF resistance is an important implication for higher risk of multi-drug resistance as it is a surrogate marker of MDR [26]. Nevertheless, there are also studies [19, 20] that show isoniazid mono-resistance as the leading mono-resistant than RIF.

The 15.9% (11/69) drug resistant detection rate using BACTEC MGIT 960 in this study was at greater proportion than similar studies [30, 31]. On the contrary, less proportion of drug resistances were reported in this study than other research reports by using the same detection system [32, 33] which might be due to different study areas and sample size. Most of the isolates detected as resistant by LPA were also resistant by MGIT. In 2 of the samples which were detected as resistant by MGIT but not by LPA could be due to its extra drug sensitivity detection test. All the MDR resistant samples by MGIT in this study were not restricted to isoniazid and rifampicin, but had an additional resistance to EMB, STM or both. This implied sensitivity detection by MGIT posses an extra benefit to detect other first line TB drugs.

The high sensitivity of the MTBDRplus assay (100%) in detecting RIF resistance in this study is the same with similar study settings in Southern [34] and northwestern [35] Ethiopia, but with little variation to the isolates from central Ethiopia [19]. The specificity of the assay (98.3%) in detecting RIF resistance is comparable to 99.2% [34], 99.8% [19] and 100% [35] in detecting RIF resistant isolates. There is a great variation of MTBDRplus assay sensitivity to INH in our study 75% and a study from the Southern Ethiopia 33.3% [34]. The sensitivity of the assay to INH was 82.7 in a study of isolates from central [19] and 91.7% in isolates collected from northwestern part of the country [35]. The susceptibility of MTBDRplus to INH in this finding 98.2% corresponds to similar reports (97.2–100%) in the country. Although the assay’s sensitivity result to detect MDR have some discrepancy, the susceptibility result accounts 100% which coincides exactly with other studies reported from the country [19, 34, 35].

The study also found that there was an excellent agreement between BD BACTEC MGIT 960 and MTBDRplus assay in detecting RIF and MDR with a Kappa value of 0.925 and 0.901, respectively. However, it was found as a good agreement (Kappa value, 0.774) between the two testing methods in detecting INH. Similarly, it was also reported as an excellent agreement in detecting MDR (Kappa value, 1) but with good (0.663) and moderate (0.494) agreement in detecting RIF and INH resistance, respectively [34]. Such difference between the two findings might be due to the difference in number of drug resistant isolates. Further, the agreement between Culture MGIT and LPA in detecting katG, inhA and both (katG and inhA) resistance was almost perfect and agrees with a study from Kenya [17].

Frequency of gene mutations associated with rifampicin and isoniazid resistant isolates by GenoType MTBDRplus assay

The wild-type rpoB probe hybridization band pattern showed an omission of the bands at WT3, WT4, WT6, WT7 and WT8 probes. Of these, the miss of WT8 probe predominates and accounted for 71.4% (10/14) among the RIF resistant isolates in a small region of amino acids located between position numbers 530–533 of the rpoB gene. Omission of all the rpoB gene probes was without any gain of probes in the MUT region and depicted as “unknown” mutation. Such lack of binding of a WT probe without simultaneous binding of a mutant probe is likely caused by the presence of a resistant mutation. Similarly, missing of probes without any gain of the corresponding MUT was reported from northwestern Ethiopia [20]. A report from central Ethiopia also in line with our finding that no gain of any mutant band was identified for the non-hybridized WT8 probe [19]. Similar findings were also reported for the absence of MUT band from other countries [36, 37]. Despite such reports, there were research findings where gain of MUT probes was identified in a part of the rpoB gene [21].

Greater frequency of resistance 8/16 (50%) to INH occurred due to mutation of the katG gene, whereas lower frequency of resistance 5/16 (31.3%) was caused by the mutations in the promoter region of the inhA gene. Greater frequencies of the katG gene WT omission at codon 315 was also reported from other studies in Ethiopia that agrees with this finding [19,20,21]. A miss of WT probe at Codon 315 in the katG gene without the presence of specific MUT band accounted 75% (6/8) in this study. The remaining 25% of the strains had mutations in the katG gene at codon 315 with amino acid change of S315T1 (AGC → ACC). The inhA promoter region showed a miss of − 15/− 16 codon that was detected in 5/16 (31.3%); and 1/16 (6.3%) at − 8 gene region without conferring any addition of specific mutational band in both cases. This mutation also agrees with the finding from central Ethiopia and also considered as ‘unknown’ mutation [19].

Genetic diversity of the drug resistant Mycobacterium tuberculosis

The higher proportion of Euro-American lineages among drug resistance strains could be due to its higher prevalence 63.4% (59/93) under the interpretable spoligotyping results of the LPA done isolates in this study. This in line with findings that reported greater proportion of the same most frequently reported major lineage from other parts of Ethiopia [20, 38]. The clustered strains of the study also showed higher frequency and a statistically significant association with any anti-TB drug resistance than the unique ones revealing an implication for the highest risk of drug resistance among recent TB transmission.

Conclusions

Resistance of M. tb to common drugs showed agreeable patterns in both conventional phenotypic BD BACTEC MGIT 960 and the GenoTypic MTBDRplus assay test. Both testing methods highlighted greater proportion of drug resistant isolates among newly treated patients than retreatment ones. Rifampicin mono-resistance was accounted for greater proportion of resistance than INH mono-resistance in the study and the majority of the RIF resistant isolates were also found as multi-drug resistant. Euro-American major lineages as the predominant lineage with most of the drug resistant isolates as orphans. Thus, further consideration should be given toward drug sensitivity testing methods in a large number of LJ-positive samples and to identify the circulating strains for better prevention and control programs which might realizes the targeted “End TB strategy”.

Limitations

In this study, the less number of culture positivity from large smear positive samples due to various reasons and lack of second line drug sensitivity test was considered as its limitation.

Availability of data and materials

All the data sets on which our conclusions relayed on were presented in the main section of this manuscript.

Abbreviations

AC:

Amplification control

ALIPB:

Aklilu Lemma Institute of Pathobiology

BHDL:

Biqat higher diagnostic laboratory

BMH:

Boru Meda hospital

BTHC:

Bati Town health centre

CBN:

Conformal Bayesian network

CRHC:

Chefa Robit health center

DHC:

Dessie health center

DNA:

Deoxyribonucleic acid

DOT:

Directly observed therapy

DRH:

Dessie referral hospital

DST:

Drug sensitivity test

EA:

Euro-American

EAI:

East-African-Indian

EMB(E):

Ethambutol

EPTB:

Extra-pulmonary tuberculosis

FNA:

Fine needle aspirate

GC:

Growth control

GU:

Growth unit

HIV:

Human Immunodeficiency Virus

INH(I):

Isoniazid

IO:

Indo-Oceanic

IQR:

inter-quartile range

KBBN:

Knowledge based Bayesian network.

LJ:

Löwenstein Jensen

LPA:

Line probe assay

M. tb :

Mycobacterium tuberculosis

MDR:

Multi-drug resistance

MGIT:

Mycobacterium growth indicator tubes

MUT:

Mutant

NPV:

Negative predictive value

PCR:

Polymeraze chain reaction

PPV:

Positive predictive value

PZA(Z):

Pyrazinamide

RIF(R):

Rifampicin

RNA:

Ribonucleic acid

RR:

Rifampisin resistance

SD:

Standard deviation

SIT:

Spoligotype International Type

Spoligotyping:

Spacer Oligonucleotidetyping

SPSS:

Statistical Package for the Social Sciences

STM(S):

Streptomycin

TB:

Tuberculosis

WHO:

World Health Organization

WT:

Wild type

XDR:

Extensively drug-resistant

References

  1. 1.

    WHO. Global tuberculosis report 2017. 2017. doi: ISBN 9789241565516.

  2. 2.

    Morris MD, Quezada L, Bhat P, Moser K, Smith J, Perez H, et al. Social, economic, and psychological impacts of MDR-TB treatment in Tijuana, Mexico: a patient's perspective. Int J Tuberc Lung Dis. 2013;17(7):954–60.

    CAS  Article  Google Scholar 

  3. 3.

    WHO. Global tuberculosis report 2019. 2019. doi: ISBN 9789241565714.

  4. 4.

    WHO. Global tuberculosis report 2018. 2018. doi: ISBN 9789241565646.

  5. 5.

    Abanda NN, Djieugoué JY, Lim E, Pefura-Yone EW, Mbacham WF, Vernet G, et al. Diagnostic accuracy and usefulness of the Genotype MTBDRplus assay in diagnosing multidrug-resistant tuberculosis in Cameroon: a cross-sectional study. BMC Infect Dis. 2017;17:379.

    Article  Google Scholar 

  6. 6.

    WHO. Molecular assays intended as initial tests for the diagnosis of pulmonary and extra-pulmonary TB and rifampicin resistance in adults and children: rapid communication. Policy update. Geneva: World Health Organization; 2020.

    Google Scholar 

  7. 7.

    Gori A, Bandera A, Marchetti G, Esposti A, Lidia Catozzi L, Nardi G, et al. Spoligotyping and Mycobacterium tuberculosis. Emerg Infect Dis. 2005;11(8):1242–8.

    CAS  Article  Google Scholar 

  8. 8.

    Birhanu T, Belay M, Ameni G. Molecular characterization of Mycobacterium tuberculosis complex isolated from Tuberculous lymphadenitis patients at Dessie private hospitals, Northern Ethiopia. JBAH. 2014;4(15):125–32.

    Google Scholar 

  9. 9.

    Berg S, Schelling E, Hailu E, Firdessa R, Gumi B, Erenso G, et al. Investigation of the high rates of extrapulmonary tuberculosis in Ethiopia reveals no single driving factor and minimal evidence for zoonotic transmission of Mycobacterium bovis infection. BMC Infect Dis. 2015;15:112.

    Article  Google Scholar 

  10. 10.

    Wright C. Fine-needle aspiration biopsy of lymph nodes. CME. 2012;30(2):56–60.

    Google Scholar 

  11. 11.

    Tripathi K, Tripathi PC, Nema S, Shrivastava AK, Dwiwedi K, Dhanvijay AK. Modified Petroff’s method: an excellent simplified decontamination technique in comparison with Petroff’s method. Int J Recent Trends Sci Technol. 2014;10(3):2249–59.

    Google Scholar 

  12. 12.

    Boehme CC, Nabeta P, Hillemann D, Nicol MP, Shenai S, Krapp F, et al. Rapid molecular detection of tuberculosis and rifampin resistance. N Engl J Med. 2010;363:1005–15.

    CAS  Article  Google Scholar 

  13. 13.

    GenoType MTBDRplus - GHDonline. https://www.ghdonline.org/uploads/MTBDRplusV2_0212_304A-02-02.pdf.

  14. 14.

    Aktas AE, Yigit N, Ayyildiz A, Bastopcu A. Comparison of the Mycobacterium growth Indicator tube method and the method of proportion for drug susceptibility testing of Mycobacterium tuberculosis. Eurasian J Med. 2014;46:96–101.

    Article  Google Scholar 

  15. 15.

    Siddiqi S, Ahmed A, Asif S, Behera D, Javaid M, Jani J, et al. Direct Drug Susceptibility Testing of Mycobacterium tuberculosis for Rapid Detection of Multidrug Resistance Using the Bactec MGIT 960 System: a Multicenter Study. J Clin Microbiol. 2012:435–40.

  16. 16.

    McHugh ML. Interrater reliability: the kappa statistic. Biochem Med. 2012;22(3):276–82.

    Article  Google Scholar 

  17. 17.

    Aricha SA, Kingwara L, Mwirigi NW, Chaba L, Kiptai T, Wahogo J, et al. Comparison of GeneXpert and line probe assay for detection of Mycobacterium tuberculosis and rifampicin mono resistance at the National Tuberculosis Reference Laboratory, Kenya. BMC Infect Dis. 2019;19:852.

    CAS  Article  Google Scholar 

  18. 18.

    Tamrat F. Molecular epidemiology and drug sensitivity of Mycobacterium tuberculosis isolates among new pulmonary tuberculosis patients in Arsi, Oromia region, Ethiopia. MSc Thesis. Ethiopia: Addis Ababa University; 2016.

    Google Scholar 

  19. 19.

    Omer Z, Mekonnen Y, Worku A, Zewde A, Medhin G, Mohammed T, et al. Evaluation of the GenoType MTBDRplus assay for detection of rifampicin- and isoniazid-resistant Mycobacterium tuberculosis isolates in Central Ethiopia. Int J Mycobacteriol. 2017;5:475–81.

    Article  Google Scholar 

  20. 20.

    Alelign A, Zewude A, Mohammed T, Tolosa S, Ameni G, Petros B. Molecular detection of Mycobacterium tuberculosis sensitivity to rifampicin and isoniazid in South Gondar zone, Northwest Ethiopia. BMC Infect Dis. 2019;19:343.

    Article  Google Scholar 

  21. 21.

    Tadesse M, Aragaw D, Dimah B, Efa F, Abdella K, Kebede W, et al. Drug resistance-conferring mutations in Mycobacterium tuberculosis from pulmonary tuberculosis patients in Southwest Ethiopia. Inter J Mycobacteriol. 2016;5:185–91.

    Article  Google Scholar 

  22. 22.

    Diallo A, Ossoga GW, Daneau G, Lo S, Ngandolo R, Djaibé CD, et al. Emergence and clonal transmission of multi-drug-resistant tuberculosis among patients in Chad. BMC Infect Dis. 2017;17:579.

    Article  Google Scholar 

  23. 23.

    Rufai SB, Kumar P, Singh A, Prajapati S, Balooni V, Singh S. Comparison of Xpert MTB/RIF with line probe assay for detection of rifampin-Monoresistant Mycobacterium tuberculosis. J Clin Microbiol. 2014;52(6):1846–52.

    CAS  Article  Google Scholar 

  24. 24.

    Bekele S, Derese Y, Hailu H, Mihret A, Dagne K, Yamuah Y, et al. Line-probe assay and molecular typing reveal a potential drug resistant clone of Mycobacterium tuberculosis in Ethiopia. Trop Dis Travel Med Vaccines. 2018;4:15.

    Article  Google Scholar 

  25. 25.

    Daniel TM. The history of tuberculosis, Historical review. Respir Med. 2006;100:1862–70.

    Article  Google Scholar 

  26. 26.

    Aung WW, Ei PW, Nyunt WW, Swe TL, Win TL, Kim KJ, et al. Phenotypic and genotypic analysis of anti-tuberculosis drug resistance in Mycobacterium tuberculosis isolates in Myanmar. Ann Lab Med. 2015;35:494–9.

    CAS  Article  Google Scholar 

  27. 27.

    Mehari K, Asmelash T, Hailekiros H, Wubayehu T, Godefay H, Araya T, et al. Prevalence and Factors Associated with Multidrug-Resistant Tuberculosis (MDR-TB) among Presumptive MDR-TB Patients in Tigray Region, Northern Ethiopia. Can J Infect Dis Med Microbiol. 2019; doi.org/10.1155/2019/2923549.

  28. 28.

    Desikan P, Panwalkar N, Mirza SB, Chaturvedi A, Ansari K, Varathe R, et al. Line probe assay for detection of Mycobacterium tuberculosis complex: an experience from Central India. Indian J Med Res. 2017;145:70–3.

    Article  Google Scholar 

  29. 29.

    Adane K, Ameni G, Bekele S, Abebe M, Aseffa A. Prevalence and drug resistance profile of Mycobacterium tuberculosis isolated from pulmonary tuberculosis patients attending two public hospitals in east Gojjam zone, Northwest Ethiopia. BMC Public Health. 2015;15:572.

    Article  Google Scholar 

  30. 30.

    Damena D, Tolosa S, Hailemariam M, Zewude A, Worku A, Mekonnen B, Mohammed T, et al. Genetic diversity and drug susceptibility profiles of Mycobacterium tuberculosis obtained from Saint Peter’s TB specialized hospital, Ethiopia. PLoS One. 2019;14(6):e0218545.

    CAS  Article  Google Scholar 

  31. 31.

    Maru MH, Mariam S, Airgecho T, Gadissa E, Aseffa A. Prevalence of Tuberculosis, Drug Susceptibility Testing, and Genotyping of Mycobacterial Isolates from Pulmonary Tuberculosis Patients in Dessie, Ethiopia. Tuberc Res Treat. 2015; https://doi.org/10.1155/2015/215015.

  32. 32.

    Mekonnen A, Merker M, Collins JM, Addise D, Aseffa A, Petros B, et al. Molecular epidemiology and drug resistance patterns of Mycobacterium tuberculosis complex isolates from university students and the local community in eastern Ethiopia. PLoS One. 2018;13(9):e0198054.

    Article  Google Scholar 

  33. 33.

    Singh UB, Pandey P, Mehta G, Bhatnagar AK, Mohan A, Goyal V, et al. Genotypic, phenotypic and clinical validation of GeneXpert in extra-pulmonary and pulmonary tuberculosis in India. PLoS One. 2016;11(2):e0149258.

    Article  Google Scholar 

  34. 34.

    Wondale B, Medhin G, Abebe G, Tolosa S, Mohammed T, Teklu T, et al. Phenotypic and genotypic drug sensitivity of Mycobacterium tuberculosis complex isolated from south Omo zone, Southern Ethiopia. Infect Drug Resist. 2018;11:1581–9.

    CAS  Article  Google Scholar 

  35. 35.

    Tessema B, Beer J, Emmrich F, Sack U, Rodloff AC. Analysis of gene mutations associated with isoniazid, rifampicin and ethambutol resistance among Mycobacterium tuberculosis isolates from Ethiopia. BMC Infect Dis. 2012;12:37.

    CAS  Article  Google Scholar 

  36. 36.

    Addo KK, Addo SO, Mensah GI, Mosi L, Bonsu FA. Genotyping and drug susceptibility testing of mycobacterial isolates from population-based tuberculosis prevalence survey in Ghana. BMC Infect Dis. 2017;17:743.

    Article  Google Scholar 

  37. 37.

    Mohan N, Chandrasekhar PB, Padmaja IJ, Raizada N, Rao PS, Kumar BS. Genotype MTBDRplus line probe assay for rapid and direct detection of rifampicin and isoniazid resistance in Mycobacterium tuberculosis complex from sputum samples. J NTR Univ Health Sci. 2014;3(1):23–7.

    Article  Google Scholar 

  38. 38.

    Tulu B, Ameni G. Spoligotyping based genetic diversity of Mycobacterium tuberculosis in Ethiopia: a systematic review. BMC Infect Dis. 2018;18:140.

    Article  Google Scholar 

Download references

Acknowledgements

We thank the zonal and districts health bureau of Oromia Special Zone and South Wollo Zone in Amhara region who gave us permission to do the research under their health institutes. All health institutes of the sample collection sites are duly acknowledged for their primary identification and diagnosis of TB cases. We are thankful to EPHI TB laboratory and all laboratory workers for their material and technical support to do both the conventional and molecular drug sensitivity test. All study participants are also acknowledged for being part of this study. At last but not least, we thank Mr. Samuel Bekele for his technical assistance to culture our samples.

Funding

This work was supported by funds from AAU thematic research project A/C 0162230106072100101. In addition, it was also supported in part by the NIH/ Fogarty International Center Global Infectious Diseases grant D43TW009127 through H3 Africa.

Author information

Affiliations

Authors

Contributions

FG, BE and GA: conceived the idea and designed the study. FG, GA, BE and YM: Consistent follow up and coordinated the sample collection. FG: undertook the field work, drafting and writing of the manuscript. FG, BY and MA: contributed a lot in the laboratory processing of drug sensitivity test. FG and BG did data analysis and interpretation. BE, BG, GA and YM reviewed and edited the manuscript. All authors read and approved the final version of the manuscript.

Corresponding author

Correspondence to Fikru Gashaw.

Ethics declarations

Ethics approval and consent to participate

The study was carried out after obtaining ethical clearance from the Institute Review Board of College of Natural and Computational Sciences, Addis Ababa University (Ref. No. CNSDO/392/07/15) and later renewed in July, 2019 (CNSDO/668/11/2019). Letter of permission was also obtained from the health bureaus of the study sites. Prior to conducting the study, objectives of the study were clarified to the study participants and written informed consent was obtained. Subject confidentiality, any special data security requirements were maintained and ensured to the participants. Positive TB case participants were referred to the DOTs program for possible treatment and follow up.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no any competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Gashaw, F., Erko, B., Mekonnen, Y. et al. Phenotypic and genotypic drug sensitivity profiles of Mycobacterium tuberculosis infection and associated factors in northeastern Ethiopia. BMC Infect Dis 21, 261 (2021). https://doi.org/10.1186/s12879-021-05961-8

Download citation

Keywords

  • Mycobacterium tuberculosis
  • Drug resistance
  • BACTEC MGIT 960
  • MTBDRplus assay
  • Northeastern Ethiopia