Your new experience awaits. Try the new design now and help us make it even better

BRIEF RESEARCH REPORT article

Front. Public Health, 29 September 2025

Sec. Infectious Diseases: Epidemiology and Prevention

Volume 13 - 2025 | https://doi.org/10.3389/fpubh.2025.1657825

This article is part of the Research TopicWorld TB Day 2025: Yes! We Can End TB: Commit, Invest, DeliverView all 9 articles

Population structure and antibiotic resistance profiles of Mycobacterium tuberculosis isolates from Ibadan, Nigeria (2019–2020): a pilot study to improve affordable molecular diagnostic tools

Elsy Carvajal&#x;Elsy Carvajal1Samantha Escandn&#x;Samantha Escandón1Pelumi Daniel Adewole&#x;Pelumi Daniel Adewole2Bernardo Castro-RodríguezBernardo Castro-Rodríguez1ngel Sebastin Rodríguez-PazmioÁngel Sebastián Rodríguez-Pazmiño1Solon Alberto Orlando,Solon Alberto Orlando3,4Alexandra NarvaezAlexandra Narvaez5Olumuyiwa Samuel AlabiOlumuyiwa Samuel Alabi2Miguel Angel Garcia-Bereguiain
Miguel Angel Garcia-Bereguiain1*
  • 1One Health Research Group, Universidad de Las Américas, Quito, Ecuador
  • 2Department of Pharmaceutical Microbiology, Faculty of Pharmacy, University of Ibadan, Ibadan, Nigeria
  • 3Instituto Nacional de Salud Pública e Investigación, Guayaquil, Ecuador
  • 4Universidad Ecotec, Guayaquil, Ecuador
  • 5Universidad Espíritu Santo, Guayaquil, Ecuador

Nigeria ranks as the sixth country globally and the first in Africa with the highest burden of tuberculosis (TB) infection. The emergence and spread of multidrug-resistant TB (MDR-TB) strains have posed significant challenges to effective disease management in the country. In this study, 55 Mycobacterium tuberculosis (MTB) isolates from patients attending a hospital in Ibadan city (Nigeria) were selected. MTB isolates were analyzed using PCR amplification of gene fragments associated with antibiotic resistance, followed by Sanger sequencing and bioinformatics analysis. Additionally, MIRU-VNTR genotyping was performed to address population structure and transmission dynamics. Results show an association between mutations in the rpoB, inhA and gyrA genes and phenotypic resistance to rifampicin, isoniazide and fluoroquinolones in a significant percentage of the MTB isolates. However, an extended panel of genes would enable a better characterization of antibiotic resistance. The population structure of MTB in Ibadan, as determined by using MIRU-VNTR, revealed that 96.1% of the strains belong to lineage 4, distributed in the following sublineages: Uganda I (47.1%), LAM (21.6%), Cameroon (17.6%), and Ghana (9.8%). Meanwhile, 3.9% of the strains correspond to lineage 5 (L5), West African-1 sub-lineage. The population structure was very heterogeneous and no active transmission clusters were detected. Overall, this pilot study demonstrated the utility of cost-effective molecular tools in enhancing TB surveillance and control programs in settings where whole-genome sequencing (WGS) is still an economical challenge.

Introduction

Nigeria is one of the most densely populated countries in West Africa, where significant social and economic inequalities contribute to inadequate healthcare coverage and inequitable access to medical services (1). Data indicate that less than 3% of Nigerians are enrolled in the National Health Insurance Scheme (NHIS), leaving 97% of the population, including vulnerable and disadvantaged groups, without access to this benefit (2). According to the World Health Organization (WHO), Nigeria ranks as a high-burden country for TB, placing sixth globally and first in Africa. In 2023, the WHO reported a total of 499,000 TB incidence cases, representing a 16% increase compared to the 418,000 cases recorded in 2015 (3). The increase in cases has exacerbated the TB crisis in the country, due to the emergence and spread of drug-resistant (DR-TB) and multidrug-resistant (MDR-TB) strains, posing a major public health challenge (4).

MDR-TB is characterized by resistance to isoniazid and rifampicin, the two most potent first-line anti-TB drugs. In 2023, Nigeria reported an estimated 9,400 cases of MDR-TB or rifampicin-resistant tuberculosis (RR-TB) (3). Alarmingly, 57% of individuals with bacteriologically confirmed TB that year exhibited resistance to both rifampicin and fluoroquinolones, meeting the criteria for pre-XDR-TB or XDR-TB (5). These figures highlight the critical severity of MDR-TB in Nigeria, emphasizing the urgent need for more effective strategies in diagnosis, treatment, and prevention. The situation is further exacerbated by poor adherence to anti-TB programs, widespread self-medication, and challenges in completing treatment due to the high cost and frequent shortages of anti-TB drugs (6).

The alarming rise of MDR-TB and XDR-TB in Nigeria is driven by the accumulation of spontaneous mutations in Mycobacterium tuberculosis (MTB), which lead to significant phenotypic changes, particularly in antibiotic resistance. Mutations in key genes such as rpoB, katG, and gyrA play a crucial role in the bacterium’s ability to evade antimicrobial treatments. For instance, mutations in rpoB alter RNA polymerase, reducing its affinity for rifampicin and leading to resistance. Similarly, mutations in katG prevent the activation of isoniazid, one of the most effective anti-TB drugs. Meanwhile, mutations in the promoter region of inhA reduce the drug’s ability to inhibit the enoyl-ACP reductase enzyme, further contributing to resistance to isoniazid. Additionally, mutations in gyrA, which encode a subunit of DNA gyrase, confer resistance to fluoroquinolones, an important second-line treatment. Beyond conferring drug resistance, these genetic adaptations enhance the bacterium’s survival and persistence in diverse environments, further complicating TB control and eradication efforts (7, 8). According to the literature, the following regions in the rpoB, katG, and gyrA genes are considered hotspots for mutations associated with antibiotics resistance, particularly within the rifampin-resistance determining region (RR-DR) between codons 507 and 533 (911). In katG, mutations (particularly at codon 315) can lead to isoniazid resistance by impairing the activation of the drug (10, 11). The gyrA gene, which encodes the A subunit of DNA gyrase, frequently exhibits mutations in its quinolone-resistance determining region, resulting in fluoroquinolone resistance (9). These genetic changes highlight the complexity of resistance mechanisms and emphasize the importance of monitoring these loci in clinical and diagnostic settings (11). These genetic changes not only confer resistance but also influence the bacterium’s adaptability and survival under selective pressures, complicating eradication efforts. Monitoring these genetic loci is crucial for the early detection and effective management of drug-resistant tuberculosis (11). Given the increasing prevalence of MDR and XDR-TB strains, understanding the molecular basis of resistance is imperative for the development of innovative diagnostic tools and treatment strategies (11).

Molecular epidemiology utilizes molecular techniques to analyze the population structure and transmission dynamics of MTB in a particular setting. Currently, various molecular and sequencing methodologies are available for epidemiological studies, including single-nucleotide polymorphism (SNPs) analysis, mycobacterial interspersed repetitive unit variable number tandem repeats (MIRU-VNTR), and whole genome sequencing (WGS) (12). Most of these tools are reproducible and provide a more comprehensive view of epidemiological studies, highlighting WGS as the molecular characterization technique with the most significant potential for discrimination in the diagnosis of M. tuberculosis (13). MIRU-VNTR is structured with 40-100-bp repetitive sequences organized in direct tandem repeats across various loci within the chromosome of MTB H37Rv (with 41 loci present). Functioning akin to minisatellites, these sequences exhibit 12 loci MIRUs that demonstrate variations in tandem repeat copy numbers and sequence variations between repeats (14). The method involves PCR amplification followed by gel electrophoresis, making it a robust and straightforward tool for genotyping. MIRU-VNTR stands out due to its simplicity high reproducibility, cost-effectiveness, and ability to identify species within the Mycobacterium tuberculosis complex (MTBC) and their genotypes, making it particularly accessible for developing countries and facilitating easy integration into a global database for strain comparison worldwide1 (15).

This study aimed to analyze the population structure and mutation profiles of the rpoB, inhA, katG, and gyrA genes in MTB isolates from patients attending a hospital facility in Ibadan, Nigeria, during 2019 and 2020.

Materials and methods

Study population and setting

The study was conducted at the Government Chest Hospital in Ibadan, Oyo State, Nigeria, between December 2019 and October 2020, and involved MTB culture from pulmonary tuberculosis patients. The inclusion criteria were that participants had to consent to be included in the study. No personal information about the patients was collected for this study, and further analysis was focused on MTB cultures isolated from these patients.

MTB culture and drug susceptibility testing

MTB culture and drug susceptibility (first- and second-line drugs) testing from sputum samples was carried out as described elsewhere (11, 16). A total number of 72 MTB isolates were collected during the study period and were sent from Ibadan to Universidad de Las Américas (Ecuador) for further analysis (Table 1).

Table 1
www.frontiersin.org

Table 1. Drug susceptibility testing for MTB isolates included in the study.

DNA extraction

Genomic DNA from MTB isolates was extracted using the boiling method (17). Briefly, MTB cultures were resuspended in TE buffer (10 mM Tris–HCl, 1 mM EDTA, pH 8.0), boiled at 95 °C for 45 min, and centrifuged at 10,000 g for 5 min. Supernatant was used directly for molecular assays. Of the 72 samples received from Ibadan, good quality DNA was recovered for 51 samples. Of those 51 samples, 35 MTB strains that were phenotypically resistant were included in further analysis for PCR and Sanger sequencing.

PCR and sanger sequencing for rpoB, inhA, katG, and gyrA genes

Gene segments associated with antibiotic resistance were amplified using the primers listed in Supplementary Table 1.

For the rpoB gene (RR-DR), the PCR mixture contained DNA, 7.5 μL of 2X Gotaq Green Master Mix (Promega), 0.15 μL of each primer (0.13 μM), and 6.1 μL distilled water for a final volume of 15 μL. The cycling program consisted of 2 min at 94 °C, followed by 30 cycles of 45 s at 94 °C, 1 min at 65 °C, and 60 s at 72 °C, with a final elongation at 72 °C for 5 min during one cycle.

For the inhA and katG genes, the PCR mixture consisted of 7.5 μL of 2x Gotaq Green Master Mix (Promega), 1.5 μL of each primer (1 μM), 3.5 μL distilled water, in a final volume of 15 μL. The cycling program consisted of 2 min at 94 °C, followed by 30 cycles of 45 s at 94 °C, 1 min at 65 °C, 1 min at 72 °C, with a final elongation at 72 °C for 5 min during one cycle.

For the gyrA gene, the PCR mixture contained 7.5 μL of 2X Gotaq Green Master Mix (Promega), 0.15 μL of each primer (0.4 μM), 5.3 μL distilled water, in a final volume of 15 μL. The cycling program was set at 5 min at 94 °C, followed by 35 cycles of 45 s at 94 °C, 30 s at 55 °C, and 50 s at 72 °C, with a final elongation at 72 °C for 10 min during one cycle.

PCR products were amplified in 2% agarose gels with a 100 bp ladder (ABM) to confirm amplification visually.

Sanger sequencing was performed at Laboratorios de Investigación of the Universidad de Las Américas (Quito, Ecuador), using the DNA Sequencing by Capillary Electrophoresis protocol described by ThermoFisher (18). The genetic analysis and the sequence alignment were performed using the software Geneious (version 2023.0.2), with the MTB H37Rv genome as the reference (accession number: NC_000962.3), downloaded from the NCBI Genome database.2

MIRU-VNTR genotyping of MTB isolates

The MIRU-VNTR genotyping was conducted through PCR amplification of 12 loci followed by visualization in 2% agarose gels, as described previously. Genotyping data of the isolates were analyzed using the MIRU-VNTRplus web application available on: www.miru-vntrplus.org. Sublineage identification was performed by similarity search using MIRU-VNTR information through categorical distance measure (MIRU-VNTR weight: 1). The neighbor-joining tree (NJT) and calculation of Minimum Spanning Tree (MST) were also performed using MIRU-VNTR information.

Results

Mutation analysis of rpoB, inhA, katG, and gyrA genes

The drug susceptibility testing results for the 51 MTB isolates with good DNA quality are detailed in Figure 1. There were 35 MTB isolates resistant to rifampicin, seven isolates resistant to isoniazid and five isolates resistant to fluoroquinolones.

Figure 1
A categorical NJ tree diagram depicts genetic relationships among different lineages, labeled as Uganda I, Cameroon, Ghana, LAM, and West African 1, with specific identifiers for each branch. The tree illustrates cluster similarities, with a numerical matrix on the right showing associated genetic data. A color-coded legend maps lineages to colors: blue for Uganda I, red for Cameroon, yellow for Ghana, green for LAM, and light blue for West African 1. A scale bar represents genetic distance.

Figure 1. The population structure of 51 MTB strains from Ibadan Nigeria included in the study. The Neighbor-joining tree was constructed on 12 –loci MIRU-VNTR with default parameters. Each color represents a lineage. The label on the samples contains the resistance of each strain. The analysis was performed using the MIRU-VNTRplus web application.

Sanger sequencing results of the rpoB gene for the 35 MTB strains resistant to rifampicin are showed in Table 2. Twelve MTB strains (34.3%) were discarded due to poor sequencing quality. For the remaining 23 MTB strains, 10 (43.5%) had mutations in the RR-DR region of the rpoB gene. Mutations in codon 531 were found in 7 strains, while mutations in codons 513, 526 and 533 were found in one single MTB isolate each (Table 2). In 13 MTB strains (56.5%) that were phenotypically rifampicin-resistant, there were no mutations in the RD-DR region.

Table 2
www.frontiersin.org

Table 2. Mutations analysis of rpoB, katG, inhA, and gyrA for the MTB isolates included in the study.

Sanger sequencing results of the katG gene and inhA promoter for the 7 MTB strains resistant to isoniazid. While no mutations were found in katG, 2 MTB strains (28.6%) had mutations in inhA (Table 2). In 5 MTB strains (71.4%) that were phenotypically isoniazid-resistant, there were no mutations found in either in katG or inhA.

Sanger sequencing results of the rpoB gene for the 5 MTB strains resistant to fluoroquinolones are shown in Table 2. One MTB strain (20%) was discarded due to poor sequencing quality. For the remaining 4 MTB strains, all of these (100%) had mutations in gyrA. The four MTB strains had mutations in codon 95, and one of them also had a mutation in codon 94 (Table 2).

The population structure of MTB isolates from Ibadan (Nigeria) included in this study

The population structure for the 51 MTB isolates included in the genotypic analysis is shown in the neighbor-joining tree (NJT) in Figure 1. Lineage identification using similarity search in MIRU-VNTRplus showed that 96.1% (49/51) of the strains belong to lineage 4 (L4), distributed as follows: Uganda I (24/51, 47.1%), LAM (11/51, 21.6%), Cameroon (9/51, 17.6%), and Ghana (5/51, 9.8%). Meanwhile, 3.9% (2/51) of the strains correspond to lineage 5 (L5), specifically the West African-1 sub-lineage. The population structure was very heterogeneous, and no active transmission clusters were detected through 12-MIRU-VNTR typing (Figure 1).

The minimum spanning tree (MST) was done with a maximum difference of two loci for the detection of clonal complexes. This analysis shows 7 well-defined clonal complexes (Figure 2). Complex 1 consists of 18 MTB strains. Within this complex, there is a predominance of the Uganda I sublineage, and one strain has been identified as belonging to the LAM sublineage. Complexes 2 to 7 are made up of 2 to 3 strains. And 20 strains were singletons.

Figure 2
Network diagram showing complexes labeled One through Six, with nodes connected by arrows. Nodes are color-coded to represent different lineages: red (Cameroon), green (Ghana), yellow (LAM), blue (Uganda), and purple (West African 1). A legend explains the colors. Complex One is central with branching pathways.

Figure 2. Minimum spanning tree (MST), based on 12-loci MIRU-VNTR analysis for the 51 MTB strains for Ibadan (Nigeria) included in this study. Genotypic lineages are distinguished by various colored dotted circles. The maximum locus difference between clonal complexes (CCs) is 2. The analysis was performed using the MIRU-VNTRplus web application.

In Supplementary Table 2, the drug susceptibility profiles across the different MTB sublineages are detailed. The drug resistance MTB strains are concentrated in lineage 4, especially in strains from the sublineage Uganda I.

Discussion

Nigeria has the highest tuberculosis infection rate in Africa, with 268 deaths reported daily. In 2020, the year this study’s samples were collected, the WHO recorded 138,591 TB cases in the country (a 15% increase from the previous year) (19). Despite this alarming situation, Nigeria has one of the world’s lowest TB detection rates, with only 24% of cases reported to the National Tuberculosis and Leprosy Control Program, further hindering efforts to combat the disease (20).

Our results for the mutation analysis of drug resistance genes are similar to previous report from Nigeria. A previous study from Ibadan reported a rifampicin resistance prevalence of 23.3%, closely aligning with the 28.6% observed in our study (21). Notably, missense mutations were identified, with the most prominent occurring at codon 531, as we also found (22). This mutation confers high-level rifampicin resistance by binding to the RNA polymerase subunit and inhibiting DNA-directed RNA synthesis in M. tuberculosis (22, 23). Moreover, mutations at codons 513, 526, and 533 found in our study have also been previously reported in Nigeria (24) and are associated with low-level rifampicin resistance (25).

Analysis of a gyrA gene fragment for fluoroquinolone resistance revealed a 60% incidence of the Ser95Thr mutation at codon 95 and a 20% occurrence of the Asp94Gly mutation at codon 94. Previous studies worldwide, including Nigeria, have already shown that mutations at codons 94 and 95 are the most common (2628).

The C (−15 T) mutation in the inhA promoter, detected in 20% of strains, is a well-established marker for MDR-TB. Along with three other mutations (−16A > G, −8 T > A, and −8 T > C), it contributes to resistance by upregulating inhA expression (27). A similar study conducted in Southwestern Nigeria (28) has already reported that 89.7% of samples with phenotypic isoniazid resistance carried at least one known resistance-associated mutation. In contrast, our study did not identify any mutations within the amplified region of katG, suggesting that resistance may stem from mutations outside the analyzed section (29). Notably, up to 20% of isoniazid-resistance mutations occur outside katG, often in the ahpC gene (30, 31). Additionally, a significant proportion of clinical isolates exhibit resistance that cannot be fully explained by mutations in the analyzed genes, suggesting the involvement of alternative mechanisms such as efflux pump overexpression (32). M. tuberculosis contains multiple efflux pumps capable of extruding different classes of antibiotics. These multidrug resistance efflux pumps and transporters may be overexpressed due to mutations in regulatory regions, or their expression can be induced by antibiotic exposure through interactions with bacterial regulatory systems (33, 34). Furthermore, during M. tuberculosis infection, certain bacterial subpopulations can develop phenotypic tolerance to antimycobacterial drugs without acquiring genetic mutations (35). Another plausible cause is the presence of mutations located outside the routinely screened regions, such as outside the RRDR in the rpoB gene, which have been associated with low-level rifampicin resistance and may escape detection by standard molecular assays (e.g., I491 > F and V170 > F) (36). Moreover, heteroresistance, the coexistence of resistant and susceptible bacterial subpopulations within the same sample, has been increasingly recognized and may significantly contribute to phenotypic–genotypic discordance (37).

MIRU-VNTR analysis in this study identified L4 as the most prevalent lineage (96.1%), encompassing sub-lineages such as LAM, Ghana, Uganda I, and Cameroon. Previous studies conducted in various Nigerian states (8, 38), have also reported L4 as one of the most common lineages. However, the distribution of sublineages within L4 differs between different locations and study periods in Nigeria. For instance, a previous study from 2012 in Ibadan did not detect Uganda I, the most frequent sublineage in our study (39). However, another study done in Lagos reported a prevalence of 11.5% for Uganda I (30). On the other hand, the Cameroon sublineage has a lower prevalence in our study compared to previous reports from Nigeria (39, 40). Overall, these results highlight the significant variability of MTB sublineages over recent years and across different areas of Nigeria, likely fueled either by active transmission in certain settings and also by TB relapse in a highly exposed population across the country, as it has been reported in other settings (4143). In fact, the lack of active transmission clusters within our study population at a low-resolution level, like 12 loci MIRU-VNTR supports the hypothesis of a very high diversity of MTB strains within Ibadan and the role of TB relapse in the regional epidemiology of TB in this location (4143). Nevertheless, the diversity of sublineages and the lack of clustering may also be influenced by the limitations of the MIRU-VNTR technique itself, as it is not a perfect method for genotyping and is affected by homoplasy events in these molecular markers (12, 44). Therefore, MIRU-VNTR should be complemented with other methods, such as SNP analysis, spoligotyping, or WGS, to enhance analytical accuracy and improve lineage designation when possible.

The high burden of multidrug-resistant tuberculosis (MDR-TB) in developing countries such as Nigeria reflects systemic weaknesses in both diagnostic capacity and access to effective treatment compared to developed countries. In many low-income settings, limited availability of rapid resistance testing tools, such as GeneXpert MTB/RIF, delays appropriate therapy and contributes to the continued transmission of resistant strains (45). Common resistance-conferring mutations in rpoB, katG, and inhA genes, which are associated with resistance to rifampicin and isoniazid, are prevalent in Nigerian isolates, further complicating treatment and reflecting local selective pressures resulting from inconsistent therapy and informal antibiotic use.

The population structure in Nigeria, characterized by widespread and heterogeneous latent TB infection, contrasts with the localized, risk-group-driven patterns seen in high-income countries. This facilitates the sustained transmission and reactivation of diverse, often resistant MTB lineages (46). As a result, control strategies that focus primarily on interrupting transmission chains may be insufficient. Effective TB control in Nigeria requires an integrated approach that combines molecular surveillance, mutation-informed treatment protocols, and interventions targeting social determinants such as poverty and malnutrition (47). This highlights the need for comprehensive public health interventions that go beyond biomedical solutions and target the social conditions underpinning TB vulnerability.

This study emphasizes the implementation of cost-effective molecular approaches for the surveillance of drug-resistant M. tuberculosis, particularly resistance to first-line drugs, in resource-constrained settings. Overall, our results for mutation analysis of key antibiotic resistance genes support its use as an affordable method to address drug resistance analysis of MTB in the absence of WGS. However, a larger panel of target genes should be considered for a more sensitive detection of resistant strains. Locally, in Nigeria, where TB burden is high and public health budgets are limited, the adoption of in-house PCR protocols and targeted sequencing offers a practical solution to improve early detection of resistance and guide timely treatment decisions provides a feasible solution for improving early detection, enabling better tracking of transmission dynamics and drug resistance patterns, and informing timely treatment decisions (45, 48). Regionally, across Africa, where similar challenges in laboratory infrastructure persist, these scalable tools can support the identification of outbreaks, monitoring of treatment failures, and refinement of national control strategies (49). Globally, the deployment of affordable molecular genotyping contributes to enhanced surveillance and a better understanding of resistance evolution, supporting the WHO End TB Strategy by enabling data-driven responses in high-burden regions. Thus, the approaches proposed in this research offer a viable path forward for TB control at multiple levels of intervention (50).

One of the main limitations of this study is the relatively small sample size (n = 51), obtained from a single center, which may restrict the representativeness of the findings. Therefore, the results should be interpreted with caution and cannot be generalized to the broader population. Moreover, the analysis was limited to four resistance-associated genes, which may not fully capture the genetic basis of phenotypic resistance observed in some isolates. Future studies should include additional loci such as embB, pncA, rpsL, and rrs, in order to provide a more comprehensive understanding of drug resistance in the studied population. Also, whole genome sequencing should be considered in future studies for a deeper characterization of MTB isolates.

Conclusion

Continuous epidemiological monitoring of Mycobacterium tuberculosis in Nigeria is essential, particularly given the country’s dense population, challenging socioeconomic conditions, and high levels of human mobility. Ongoing surveillance is vital for tracking the prevalence, distribution, and transmission dynamics of tuberculosis. It supports the identification of outbreaks, the design of more effective prevention and treatment strategies, the detection of drug resistance, and the evaluation of intervention effectiveness. Moreover, given the high burden of TB in Nigeria, there is an urgent need for further studies to understand better and address the challenges posed by the disease. In this sense, implementing affordable methods for monitoring drug resistance mutations monitoring, and MTB population dynamics would help to improve TB surveillance and control in a context where WGS is still an economical challenge.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.

Ethics statement

The studies involving humans were approved by the Oyo State Ministry of Health Ethics Committe (code AD13/479/1548). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

EC: Investigation, Data curation, Writing – review & editing, Writing – original draft, Methodology, Formal analysis. SE: Formal analysis, Writing – original draft, Data curation, Investigation, Writing – review & editing, Methodology. PA: Writing – review & editing, Conceptualization, Investigation, Methodology, Resources, Funding acquisition. BC-R: Methodology, Validation, Investigation, Supervision, Writing – review & editing. ÁR-P: Methodology, Formal analysis, Investigation, Writing – review & editing, Validation. SO: Resources, Methodology, Investigation, Writing – review & editing. AN: Resources, Funding acquisition, Writing – review & editing, Methodology, Investigation. OA: Writing – review & editing, Resources, Validation, Supervision, Investigation, Conceptualization. MG-B: Conceptualization, Data curation, Funding acquisition, Validation, Project administration, Resources, Formal analysis, Methodology, Writing – review & editing, Supervision, Investigation, Writing – original draft.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This study was partially funded by Universidad de Las Américas (MED.MGB.23.01).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.

Generative AI statement

The authors declare that no Gen AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpubh.2025.1657825/full#supplementary-material

Footnotes

References

1. Okeke, C, Ezenwaka, U, Ekenna, A, Onyedinma, C, and Onwujekwe, O. Analysing the progress in service delivery towards achieving universal health coverage in Nigeria: a scoping review. BMC Health Serv Res. (2023) 23:1094. doi: 10.1186/s12913-023-10090-w

PubMed Abstract | Crossref Full Text | Google Scholar

2. Onasanya, AA. Increasing health insurance enrolment in the informal economic sector. J Glob Health. 10:010329. doi: 10.7189/jogh.10.010329

PubMed Abstract | Crossref Full Text | Google Scholar

3. WHO Data. Available online at: https://www.who.int/teams/global-programme-on-tuberculosis-and-lung-health/data (Accesed September 05, 2025).

Google Scholar

4. WHO Global tuberculosis report 2021. Disponible en. Available online at: https://www.who.int/publications/i/item/9789240037021 (Accesed September 05, 2025).

Google Scholar

6. Lessem, E, Cox, H, Daniels, C, Furin, J, McKenna, L, Mitnick, CD, et al. Access to new medications for the treatment of drug-resistant tuberculosis: patient, provider and community perspectives. Int J Infect Dis. (2015) 32:56–60. doi: 10.1016/j.ijid.2014.12.012

PubMed Abstract | Crossref Full Text | Google Scholar

7. Thumamo Pokam, BD, Yeboah-Manu, D, Teyim, PM, Guemdjom, PW, Wabo, B, Fankep, ABD, et al. A pilot study on the genetic diversity of Mycobacterium tuberculosis complex strains from tuberculosis patients in the Littoral region of Cameroon. J Clin Tuberc Other Mycobact Dis. (2020) 21:100182. doi: 10.1016/j.jctube.2020.100182

PubMed Abstract | Crossref Full Text | Google Scholar

8. Makpu, SS, Otokunefor, K, and Otokunefor, TV. Genetic diversity of Mycobacterium tuberculosis complex isolates circulating in north Central Nigeria. Dordrecht, Netherlands: En Atlantis Press (2023).

Google Scholar

9. Chen, J, Chen, Z, Li, Y, Xia, W, Chen, X, Chen, T, et al. Characterization of gyrA and gyrB mutations and fluoroquinolone resistance in Mycobacterium tuberculosis clinical isolates from Hubei Province, China. Braz J Infect Dis. (2012) 16:136–41. doi: 10.1590/S1413-86702012000200005

PubMed Abstract | Crossref Full Text | Google Scholar

10. Cho, EH, Bae, HK, Kang, SK, and Lee, EH. Detection of isoniazid and rifampicin resistance by sequencing of katG, inhA, and rpoB genes in Korea. Korean J Lab Med. (2009) 29:455–60. doi: 10.3343/kjlm.2009.29.5.455

PubMed Abstract | Crossref Full Text | Google Scholar

11. Franco-Sotomayor, G, Garzon-Chavez, D, Leon-Benitez, M, de Waard, JH, and Garcia-Bereguiain, MA. A first insight into the katG and rpoB gene mutations of multidrug-resistant Mycobacterium tuberculosis strains from Ecuador. Microb Drug Resist. (2019) 25:524–7. doi: 10.1089/mdr.2018.0203

PubMed Abstract | Crossref Full Text | Google Scholar

12. Shi, J, Zheng, D, Zhu, Y, Ma, X, Wang, S, Li, H, et al. Role of MIRU-VNTR and spoligotyping in assessing the genetic diversity of Mycobacterium tuberculosis in Henan Province, China. BMC Infect Dis. (2018) 18:447. doi: 10.1186/s12879-018-3351-y

Crossref Full Text | Google Scholar

13. Acharya, B, Acharya, A, Gautam, S, Ghimire, SP, Mishra, G, Parajuli, N, et al. Advances in diagnosis of tuberculosis: an update into molecular diagnosis of Mycobacterium tuberculosis. Mol Biol Rep. (2020) 47:4065–75. doi: 10.1007/s11033-020-05413-7

PubMed Abstract | Crossref Full Text | Google Scholar

14. Supply, P, Mazars, E, Lesjean, S, Vincent, V, Gicquel, B, and Locht, C. Variable human minisatellite-like regions in the Mycobacterium tuberculosis genome. Mol Microbiol. (2000) 36:762–71. doi: 10.1046/j.1365-2958.2000.01905.x

PubMed Abstract | Crossref Full Text | Google Scholar

15. Rovina, N, Karabela, S, Constantoulakis, P, Michou, V, Konstantinou, K, Sgountzos, V, et al. MIRU-VNTR typing of drug-resistant tuberculosis isolates in Greece. Ther Adv Respir Dis. (2011) 5:229–36. doi: 10.1177/1753465811402120

PubMed Abstract | Crossref Full Text | Google Scholar

16. Franco-Sotomayor, G, Rivera-Olivero, IA, Leon-Benitez, M, Uruchima-Campoverde, SE, Cardenas-Franco, G, and Perdomo-Castro, ME. Fast, simple, and cheap: the Kudoh-Ogawa swab method as an alternative to the Petroff–Lowenstein-Jensen method for culturing of Mycobacterium tuberculosis. J Clin Microbiol. (2020) 58:e01424–19. doi: 10.1128/JCM.01424-19

Crossref Full Text | Google Scholar

17. Garzon-Chavez, D, Zurita, J, Mora-Pinargote, C, Franco-Sotomayor, G, Leon-Benitez, M, Granda-Pardo, JC, et al. Prevalence, drug resistance, and genotypic diversity of the Mycobacterium tuberculosis Beijing family in Ecuador. Microb Drug Resist. (2019) 25:931–7. doi: 10.1089/mdr.2018.0429

PubMed Abstract | Crossref Full Text | Google Scholar

19. WHO Nigeria records gain in fight against tuberculosis amidst coronavirus pandemic. (2025). Available online at: https://www.afro.who.int/news/nigeria-records-gain-fight-against-tuberculosis-amidst-coronavirus-pandemic (Accesed September 05, 2025).

Google Scholar

20. WHO Nigeria re-commits to eradicating Tuberculosis and controlling Noncommunicable Diseases by 2030. (2025) Available online at: https://www.afro.who.int/news/nigeria-re-commits-eradicating-tuberculosis-and-controlling-noncommunicable-diseases-2030 (Accesed September 05, 2025).

Google Scholar

21. Falodun, OI, Cadmus, IS, and Fagade, OE. Patterns of drug resistance in Mycobacterium tuberculosis from tuberculosis patients in Ibadan Nigeria. Micro Med. (2021) 9:8–17. doi: 10.5281/zenodo.4288506

Crossref Full Text | Google Scholar

22. Ahmad, S, Mokaddas, E, and Fares, E. Characterization of rpoB mutations in rifampin-resistant clinical Mycobacterium tuberculosis isolates from Kuwait and Dubai. Diagn Microbiol Infect Dis. (2002) 44:245–52. doi: 10.1016/S0732-8893(02)00457-1

PubMed Abstract | Crossref Full Text | Google Scholar

23. Zaw, MT, Emran, NA, and Lin, Z. Mutations inside rifampicin-resistance determining region of rpoB gene associated with rifampicin-resistance in Mycobacterium tuberculosis. J Infect Public Health. (2018) 11:605–10. doi: 10.1016/j.jiph.2018.04.005

PubMed Abstract | Crossref Full Text | Google Scholar

24. Ohiengbomwan, OT, Oguzie, J, Eromon, P, Kayode, AT, Afolabi, TS, and Komolafe, IO. Gene mutation patterns of Mycobacterium tuberculosis complex and associated factors among suspected multidrug-resistant tuberculosis patients in Osun state, south-west, Nigeria. Sci Afr. (2023) 22:e01968. doi: 10.1016/j.sciaf.2023.e01968

Crossref Full Text | Google Scholar

25. Rufai, SB, Kumar, P, Singh, A, Prajapati, S, Balooni, V, and Singh, S. Comparison of Xpert MTB/RIF with line probe assay for detection of rifampin-monoresistant Mycobacterium tuberculosis. J Clin Microbiol. (2014) 52:1846–52. doi: 10.1128/JCM.03005-13

PubMed Abstract | Crossref Full Text | Google Scholar

26. Alangaden, GJ, Manavathu, EK, Vakulenko, SB, Zvonok, NM, and Lerner, SA. Characterization of fluoroquinolone-resistant mutant strains of Mycobacterium tuberculosis selected in the laboratory and isolated from patients. Antimicrob Agents Chemother. (1995) 39:1700–3. doi: 10.1128/AAC.39.8.1700

PubMed Abstract | Crossref Full Text | Google Scholar

27. WHO Catalogue of mutations in Mycobacterium tuberculosis complex and their association with drug resistance. Available online at: https://www.who.int/publications/i/item/9789240082410 (Accesed September 05, 2025).

Google Scholar

28. Senghore, M, Otu, J, Witney, A, Gehre, F, Doughty, EL, Kay, GL, et al. Whole-genome sequencing illuminates the evolution and spread of multidrug-resistant tuberculosis in Southwest Nigeria. PLoS One. 12:e0184510. doi: 10.1371/journal.pone.0184510

Crossref Full Text | Google Scholar

29. Vilchèze, C, and Jacobs, WR. Resistance to isoniazid and ethionamide in Mycobacterium tuberculosis: genes, mutations, and causalities. Microbiol Spectr. (2014) 2:MGM2-0014-2013. doi: 10.1128/microbiolspec.MGM2-0014-2013

PubMed Abstract | Crossref Full Text | Google Scholar

30. Jagielski, T, Bakuła, Z, Roeske, K, Kamiński, M, Napiórkowska, A, Augustynowicz-Kopeć, E, et al. Detection of mutations associated with isoniazid resistance in multidrug-resistant Mycobacterium tuberculosis clinical isolates. J Antimicrob Chemother. (2014) 69:2369–75. doi: 10.1093/jac/dku161

PubMed Abstract | Crossref Full Text | Google Scholar

31. Kang, JY, Hur, J, Kim, S, Jeon, S, Lee, J, Kim, YJ, et al. Clinical implications of discrepant results between genotypic MTBDRplus and phenotypic Löwenstein-Jensen method for isoniazid or rifampicin drug susceptibility tests in tuberculosis patients. J Thorac Dis. (2019) 11:400–9. doi: 10.21037/jtd.2019.01.58

Crossref Full Text | Google Scholar

32. Laws, M, Jin, P, and Rahman, KM. Efflux pumps in mycobacterium tuberculosis and their inhibition to tackle antimicrobial resistance. Trends Microbiol. (2022) 30:57–68. doi: 10.1016/j.tim.2021.05.001

Crossref Full Text | Google Scholar

33. Louw, GE, Warren, RM, Gey van Pittius, NC, McEvoy, CRE, Van Helden, PD, and Victor, TC. A balancing act: efflux/influx in mycobacterial drug resistance. Antimicrob Agents Chemother. (2009) 53:3181–9. doi: 10.1128/AAC.01577-08

PubMed Abstract | Crossref Full Text | Google Scholar

34. Tenover, FC. Mechanisms of antimicrobial resistance in bacteria. Am J Infect Control. (2006) 34:S3–S10. doi: 10.1016/j.ajic.2006.05.219

Crossref Full Text | Google Scholar

35. Kester, JC, and Fortune, SM. Persisters and beyond: mechanisms of phenotypic drug resistance and drug tolerance in bacteria. Crit Rev Biochem Mol Biol. (2014) 49:91–101. doi: 10.3109/10409238.2013.869543

PubMed Abstract | Crossref Full Text | Google Scholar

36. Qadir, M, Faryal, R, Khan, MT, Khan, SA, Zhang, S, Li, W, et al. Phenotype versus genotype discordant rifampicin susceptibility testing in tuberculosis: implications for a diagnostic accuracy. Microbiol Spectr. 12:e0163123. doi: 10.1128/spectrum.01631-23

Crossref Full Text | Google Scholar

37. de Assis Figueredo, LJ, de Almeida, IN, Augusto, CJ, Soares, VM, Suffys, PN, Carvalho, W d S, et al. Characterization of Mycobacterium tuberculosis heteroresistance by genotyping. Int J Mycobacteriol. (2020) 9:368–72. doi: 10.4103/ijmy.ijmy_132_20

Crossref Full Text | Google Scholar

38. Bakuła, Z, Wuyep, VB, Bartocha, Ł, Vyazovaya, A, Ikeh, EI, Bielecki, J, et al. Molecular snapshot of drug-resistant Mycobacterium tuberculosis strains from the plateau state, Nigeria. PLoS One. (2022) 17:e0266837. doi: 10.1371/journal.pone.0266837

PubMed Abstract | Crossref Full Text | Google Scholar

39. Lawson, L, Zhang, J, Gomgnimbou, MK, Abdurrahman, ST, Moullec, SL, Mohamed, F, et al. A molecular epidemiological and genetic diversity study of tuberculosis in Ibadan, Nnewi and Abuja, Nigeria. PLoS One. (2012) 7:e38409. doi: 10.1371/journal.pone.0038409

PubMed Abstract | Crossref Full Text | Google Scholar

40. Noorizhab, MNF, Zainal Abidin, N, Teh, LK, Tang, TH, Onyejepu, N, Kunle-Ope, C, et al. Exploration of the diversity of multi-drug resistant Mycobacterium tuberculosis complex in Lagos, Nigeria using WGS: distribution of lineages, drug resistance patterns and genetic mutations. Tuberculosis. (2023) 140:102343. doi: 10.1016/j.tube.2023.102343

PubMed Abstract | Crossref Full Text | Google Scholar

41. Garzon-Chavez, D, Garcia-Bereguiain, MA, Mora-Pinargote, C, Granda-Pardo, JC, Leon-Benitez, M, Franco-Sotomayor, G, et al. Population structure and genetic diversity of Mycobacterium tuberculosis in Ecuador. Sci Rep. (2020) 10:6237. doi: 10.1038/s41598-020-62824-z

PubMed Abstract | Crossref Full Text | Google Scholar

42. Castro-Rodriguez, B, León-Ordóñez, K, Franco-Sotomayor, G, Benítez-Medina, JM, Jiménez-Pizarro, N, Cárdenas-Franco, G, et al. Population structure of Mycobacterium tuberculosis in El Oro: a first insight into Ecuador-Peru tuberculosis transmission. J Infect Public Health. (2024) 17:527–34. doi: 10.1016/j.jiph.2024.01.015

PubMed Abstract | Crossref Full Text | Google Scholar

43. Castro-Rodriguez, B, Espinoza-Andrade, S, Franco-Sotomayor, G, Benítez-Medina, JM, Jiménez-Pizarro, N, Cárdenas-Franco, C, et al. A first insight into tuberculosis transmission at the border of Ecuador and Colombia: a retrospective study of the population structure of Mycobacterium tuberculosis in Esmeraldas province. Front Public Health. (2024) 12:1343350. doi: 10.3389/fpubh.2024.1343350

Crossref Full Text | Google Scholar

44. Allix-Béguec, C, Harmsen, D, Weniger, T, Supply, P, and Niemann, S. Evaluation and strategy for use of MIRU-VNTRplus, a multifunctional database for online analysis of genotyping data and phylogenetic identification of Mycobacterium tuberculosis complex isolates. J Clin Microbiol. (2008) 46:2692–9. doi: 10.1128/JCM.00540-08

PubMed Abstract | Crossref Full Text | Google Scholar

45. Adepoju, P. Nigeria’s widening tuberculosis gap. Lancet Infect Dis. (2020) 20:29. doi: 10.1016/S1473-3099(19)30712-1

PubMed Abstract | Crossref Full Text | Google Scholar

46. Getahun, H, Matteelli, A, Abubakar, I, Aziz, MA, Baddeley, A, Barreira, D, et al. Management of latent Mycobacterium tuberculosis infection: WHO guidelines for low tuberculosis burden countries. Eur Respir J. (2015) 46:1563–76. doi: 10.1183/13993003.01245-2015

PubMed Abstract | Crossref Full Text | Google Scholar

47. Lönnroth, K, Jaramillo, E, Williams, BG, Dye, C, and Raviglione, M. Drivers of tuberculosis epidemics: the role of risk factors and social determinants. Soc Sci Med. (1982) 68:2240–6. doi: 10.1016/j.socscimed.2009.03.041

PubMed Abstract | Crossref Full Text | Google Scholar

48. Merker, M, Kohl, TA, Niemann, S, and Supply, P. The evolution of strain typing in the Mycobacterium tuberculosis complex. Adv Exp Med Biol. (2017) 1019:43–78. doi: 10.1007/978-3-319-64371-7_3

PubMed Abstract | Crossref Full Text | Google Scholar

49. Ogbuabor, DC, and Onwujekwe, OE. Governance of tuberculosis control programme in Nigeria. Infect Dis Poverty. (2019) 8:45. doi: 10.1186/s40249-019-0556-2

PubMed Abstract | Crossref Full Text | Google Scholar

50. WHO Implementing the end TB strategy: the essentials, (2022) Available online at: https://www.who.int/publications/i/item/9789240065093 (Accesed September 05, 2025).

Google Scholar

Keywords: Nigeria, Mycobacterium tuberculosis, antibiotic resistance, MIRU-VNTR, sublineages

Citation: Carvajal E, Escandón S, Adewole PD, Castro-Rodríguez B, Rodríguez-Pazmiño ÁS, Orlando SA, Narvaez A, Alabi OS and Garcia-Bereguiain MA (2025) Population structure and antibiotic resistance profiles of Mycobacterium tuberculosis isolates from Ibadan, Nigeria (2019–2020): a pilot study to improve affordable molecular diagnostic tools. Front. Public Health. 13:1657825. doi: 10.3389/fpubh.2025.1657825

Received: 01 July 2025; Accepted: 10 September 2025;
Published: 29 September 2025.

Edited by:

Adwoa Asante-Poku, University of Ghana, Ghana

Reviewed by:

Hai xia Wang, The Second Hospital of Shandong University, China
Onaiwu Tola Ohiengbomwan, Redeemer's University, Nigeria

Copyright © 2025 Carvajal, Escandón, Adewole, Castro-Rodríguez, Rodríguez-Pazmiño, Orlando, Narvaez, Alabi and Garcia-Bereguiain. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Miguel Angel Garcia-Bereguiain, bWFnYmVyZWd1aWFpbkBnbWFpbC5jb20=

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.