1. Introduction
Global public health continues to be threatened by tuberculosis (TB), an infectious disease caused by
Mycobacterium tuberculosis (
M. tuberculosis) and ranked among the top 10 causes of mortality globally [
1]. In 2020, the estimated number of incident cases of TB stands at 9.9 million out of which Africa had a 25% share. South Africa is one of the 30 countries with a high TB burden and has the eighth highest TB incidence globally having more than 500 cases per 100,000 population, which towers above the global average incidence of 127 cases per 100,000 population [
2].
In South Africa, TB continues to be a disease of major importance and remained the leading cause of death for three consecutive years that is from 2016 to 2018 [
3]. In 2019 alone, an estimated 360,000 South Africans became ill with TB and 58,000 people were estimated to have died from the disease [
2]. The COVID-19 pandemic's effects have undone years of progress made in reducing the number of TB deaths worldwide, with the predicted number of deaths in 2020 returning to the level of 2017 [
4]. The eradication of TB by 2035, a strategic goal of World Health Organization (WHO) cannot be actualized unless the emergence of resistance in TB is seriously addressed and controlled [
5].
Drug-resistant TB (DR-TB) has emerged as a major risk to global TB control. Different mutations in genes such as
rpoB,
katG,
inhA,
pncA,
embB,
rpsL, gyrA,
ethA and
rrs have been identified conferring resistance to TB first-line drugs, second-line drugs or injectables and fluoroquinolones [
6,
7]. Mutations in codon 507 to 533 regions of the
rpoB gene, called the rifampicin-resistance determining region (RRDR) are majorly responsible for rifampicin resistance while mutations in the
katG and
inhA genes are associated with INH resistance. The
katG and
inhA mutations give rise to high-level and low-level INH resistance respectively. Although mutations in both
gyrA and
gyrB genes are responsible for fluoroquinolone resistance,
gyrA accounts for 60%-70% of all mutations [
8,
9,
10]. Recent research has revealed that different mutations in
M. tuberculosis can confer varying levels of phenotypic resistance to anti-TB medications [
6,
7,
11]. Consequently, the aggregation of mutations at several positions has a comprehensive effect on drug resistance [
12].
In rural areas of the Eastern Cape, drug resistance and gene mutations remain major barriers to effective control and management of TB. However, there was no report on the frequency of gene mutations associated with resistance to rifampicin (RIF) and isoniazid (INH). Hence, in this study, we report the prevalence of mutations in drug resistance genes (rpoB, katG and inhA) and identified the strains and lineages of DR TB strains.
4. Discussion
Spontaneous chromosomal mutations in particular locations of the bacterial genome is generally thought to be the cause of resistance to RIF and INH which several studies have proven [
9]. It is uncertain how
M. tuberculosis strains are distributed, how much has been transmitted recently and in the past, and how DR strains transmit in rural Eastern Cape and according to our knowledge, this is the first study detailing the distribution of drug-resistant genes, mutation sites and genotypes.
Clinicians are concerned about MDR-TB produced by mutations in
M. tuberculosis, which are mostly caused by the
rpoB,
katG, and
inhA genes [
18]. RIF resistance is typically thought to as a marker for MDR-TB. Hence, screening mutations in candidate genes constitute the foremost significant step to making a definite diagnosis in drug-resistant strains. The 81 bp core region of the
rpoB gene's nucleotide sequences was examined for mutations. The prevalence of mutations in the
rpoB gene in this study was higher than in the other genes (
katG and
inhA genes). Analysis of the RIF-associated mutations revealed a prevalence of 65.8% of the
rpoB gene. This is comparable to the study done by Otchere et al., with a prevalence of 52% [
19]. Contrarily, the prevalence of mutation reported in studies [
9,
20,
21] was higher than that found in our investigation, at 93.5%, 94.9%, and 91.2%, respectively, in the
katG gene. Our analysis found that the S531L codon was the site of the majority of rifampin resistance-causing mutations in the
rpoB gene. This finding is consistent with earlier research [
9,
20,
22,
23] and may be related to the propagation of a common clone. This codon's mutation is known to be a hot area for
rpoB gene mutations in M. tuberculosis, and it has also been observed in other South African provinces [
24] showing that these mutations are prevalent in the country. This high frequency of occurrence may be due to the low fitness cost associated with
rpoB S531L [
24] and have been associated with major MDR-TB outbreaks. [
23]. Low frequency of mutation was observed in codon 526 at 3.8% in this study but higher in Uganda at 12.5% [
25]; China at 14.9% [
20]; Brazil at 9.9%, [
26] indicating that frequency varies with different geographic locations. Codons 526D and 531 co-occurred at a lesser rate in this study (2.8%) than they did in Iran (23.9%) [
27].
Due to its early significant bactericidal activity, INH is a first-line TB medication and a very important medicine for the treatment of TB. For doctors, finding mutations in the
katG or
inhA promoter region is crucial since it predicts the degree of INH resistance and helps them choose the best course of treatment [
28]. Given the high degree of INH medication resistance, the catchment regions of the clinics in this study area need to be watched for any changes that may arise during patient TB treatment. In this investigation, both the
katG S315T gene mutation and very resistant INH strains were seen. The majority of mutations occur at codon 315, which is present in 30% to 90% of INH-resistant bacteria [
16]. This claim was corroborated by the study's findings, which showed a correlation between elevated levels of INH medication resistance and S315T mutations. Several nations, such as Zambia and Brazil, have observed a tendency of codon S315T mutation in the
katG gene [
9,
29]. Mutations in
katG occurred only in codon agc/acc S315Tb in this study but the study of Jagielski et al. [
30] found that mutations occurred in eleven other codons. The prevalence of
katG S315T varied according to the geographic region: Sub-Saharan Africa (94.9%) [
9], West Africa (64%) [
31], Southeast Asia (29.3%) [
32], and the United States (38%) [
33] while the global frequency of
katG S315 is estimated at 64.2% [
34]. In addition to the
katG gene, the
inhA, fabG1, and oxyR-ahpC genes are also associated with
M. tuberculosis INH resistance. It has been discovered that 20–42% of INH-resistant bacteria carry mutations in the
inhA promoter region [
33,
34]. Previous studies have shown that polymorphisms in the promoter region of the
inhA gene cause low-level resistance to INH, which ranges from 8% to 43% [
35]. In this study, low-level resistance ratio was 31.2%, close to the high limit of the reported range. Other percentages of
inhA mutations have been found in various parts of South Africa, including Kwazulu Natal (27.5%) [
36] and Free State (13.4%) [
8]. On the other hand, INH-resistant strains were found in Zambia and Ethiopia at rates of 0.8% and 2.0%, respectively, of mutations in the
inhA promoter region [
9,
37]. Contrarily, Lempens et al. [
38], reported a lower percentage of occurrence of C-15T gene mutation with a high level of INH resistance. This suggests that mutation of the
inhA gene does not always indicate a low resistance level. Mutations in the
inhA gene also give cross-resistance to ethionamide (ETH), a second-line medication used in MDR therapy, and are therefore thought of as a surrogate marker for early diagnosis of ETH resistance [
8,
39,
40]. This is because the two drugs share the same target of action. The ETH was once a component of the treatment plan for MDR-TB in South Africa, according to the National Department of Health (NDoH) in that country. In the presence of
inhA mutations, employing ETH to treat MDR-TB would not have been successful due to this cross-resistance [
8]. Consequently, in the clinical management of MDR-TB cases displaying
inhA mutations, ETH must be excluded from the regimen. The majority of our C-15TB isolates (88.5%) included the most prevalent
inhA gene mutation, supporting the findings of Seifert et al. [
41], who stated that the most prevalent
inhA gene mutation was frequently seen in C-15TB.
The clinical and molecular characteristics of the
M. tuberculosis strains vary in different areas. This was observed in this study in consonance with the findings of Liu et al. [
12]. Investigating the evolutionary lineages of M. tuberculosis can benefit from the range of mutations. The prevalence of
rpoB,
katG, and
inhA mutations in various regions of Mthatha may aid in determining whether to standardize treatment plans or provide tailored care in each region where these mutations have been discovered. There is an indication of
M. tuberculosis strains that are constantly mutating as we observed combined mutations (
Table 2). This data can be used in the development of new anti-TB drugs.
Heteroresistance in the study area increased with time, by the third year of the study period, the heteroresistance rate was almost tripled that of the previous year (
Table 1). The
rpoB and
katG combination had the highest number of heteroresistant isolates followed by
rpoB and
inhA combination. In
rpoB and
katG combination, the mutation regions
rpoB S315L and
katG 531ST had the highest number of isolates (
Table 2). Under the selective pressure of inadequate anti-TB medication, separation into susceptible and resistant organisms most likely elucidates heteroresistance caused by infection with single strains. Several reports have documented the development of resistance as a result of insufficient treatment [
42]. It has also been proven that mixed-strain infections caused by heteroresistant bacteria might have a negative effect on treatment outcomes. Treatment of such instances with first-line anti-tuberculosis medicine may select for and increase the drug-resistant strain in the host because heteroresistance makes it possible for the drug-resistant strain to go undetected [
43]. The rate of 17.9% of heteroresistance in this study is similar to the finding by Rinder et al. [
44], who reported a rate of 17%. Other studies have reported significantly lower rates [
45]. In our study one strain of either the Beijing, LAM, or X genotype induced heteroresistance, indicating that the division of a single strain into susceptible and resistant organisms is most likely the main underlying mechanism.
The
M. tuberculosis population in this study area was genetically diverse. From the 441 clinical isolates, 23 spoligotypes were observed and classified into major
M. tuberculosis lineages; lineages 1, 2, 3, and 4 as shown in Table 3. Beijing and Euro-American (LAM, T, S, and X) strains dominate the population structure of rifampicin-resistant tuberculosis (RRTB) isolates in South Africa, which can be explained by the historical movement of strains as South Africa was located in a geographically central position on the historical trade route between East and West for hundreds of years [
46].
One of the most widespread genotypes of
M. tuberculosis found globally is the Beijing family, often known as lineage 2. It is usually linked to immune evasion and antibiotic resistance, which promotes rapid bacterial replication, spread, and transmission [
47]. The Beijing family, which is more transmissible than other families [
48] were prevalent in this study (42%). This lineage has been detected in studies reported from other parts of South Africa including Limpopo, Western Cape, and Mpumalanga [
49,
50]. According to Said et al. [
49] Beijing family is predominant in the Eastern Cape followed by LAM. According to the study of Chihota et al. [
51] based on the review of the repository and databases, South African
M. tuberculosis strains revealed the greatest diversity and greatest abundance of Beijing families. Furthermore, the association between HLA-B27 and host-pathogen compatibility has accounted for the success of the Beijing lineage in South Africa [
51]. Given the growing concern over the prevalence of Beijing strains and their success in evolving to fit into various human groups, suitable measures should be implemented for public health surveillance. The knowledge of the lineages circulating in the study area will help in understanding the drivers of drug resistance and their impact on treatment outcomes and management of TB transmission.
Beijing lineage was documented for the first time in East Asian nations, according to van Soolingen et al. [
52], with a particular spoligotype pattern distinguished by the inclusion of the last 10 spacers (spacers 34– 43). The Beijing strains were historically introduced to South Africa, not from their primary origin (China), but from its secondary origin (Indonesia), according to historical evidence corroborated by genetic data [
46]. Pokam et al. [
53] argued that the dominance of the Beijing family is related to the current influx of Asian population into the African region, as well as increased trade relations of Africans returning from business trips to China, which led to the spread in the continent. Various theories have been proposed to support the introduction of the Beijing family into Africa. It is important to actively carry out surveillance of the Beijing family to verify its heightened transmission and understand its importance in the management plan of TB in this area.
In our investigation, the LAM lineage was the next most common (18.8%). This is not unexpected given that this genotype has been found to be widespread in the provinces of the Eastern Cape and the Free State [
54,
55]. This suggests there is a continuous TB transmission strain still ongoing throughout the province. Nonetheless, the LAM genotype predominates in Gauteng, Northern Cape, and KwaZulu Natal [
49] and some other neighboring countries in Southern Africa, including Zambia and Zimbabwe [
56,
57]; but with least prevalence in Western Cape [
49] The LAM strain is predominant in KwaZulu Natal, which is a neighboring Province to Eastern Cape; this might eventually bring it to par with the predominant Beijing strain thereby causing more havoc. Hence, continuous surveillance of the genetic diversity must be carried out to profile these strains. Concerning delineation of the LAM sub-lineages, five were found in our study namely LAM 3, 4, 5, 9 and LAM II-ZWE with LAM 3 being predominant at 74.7%. On the other hand, another study in South Africa reported that out of the 12 sub-lineages delineated globally, six sub-lineages were found of which LAM4 was the most predominant [
54]. LAM11-ZWE has been reported to be a dominant subfamily in Zambia and Zimbabwe [
56] whose origin has been traced to Portugal [
51]. Two strains from the LAM 3 sub-lineage were labeled as orphans because they had no matches in the SITVIT2 database. Compared to the results obtained by Maguga-Phasha et al., LAM 3 in our investigation, which corresponds to SIT33 and SIT 719, represented 8.2% and 3.9%, respectively [
50] in Limpopo had LAM 3 (7.0%) corresponding to SIT 33 only. LAM 1, LAM 2 and LAM 6 were not reported in this study but reported elsewhere [
29].