Molecular epidemiological study of multidrug-resistant tuberculosis isolated from sputum samples in Eastern Cape, South Africa

Molecular epidemiological study of multidrug-resistant tuberculosis isolated from sputum samples in Eastern Cape, South Africa

Journal Pre-proof Molecular epidemiological study of multidrug-resistant tuberculosis isolated from sputum samples in Eastern Cape, South Africa Nolw...

610KB Sizes 0 Downloads 24 Views

Journal Pre-proof Molecular epidemiological study of multidrug-resistant tuberculosis isolated from sputum samples in Eastern Cape, South Africa

Nolwazi Londiwe Bhembe, Ezekiel Green PII:

S1567-1348(20)30014-9

DOI:

https://doi.org/10.1016/j.meegid.2020.104182

Reference:

MEEGID 104182

To appear in:

Infection, Genetics and Evolution

Received date:

12 June 2019

Revised date:

31 December 2019

Accepted date:

6 January 2020

Please cite this article as: N.L. Bhembe and E. Green, Molecular epidemiological study of multidrug-resistant tuberculosis isolated from sputum samples in Eastern Cape, South Africa, Infection, Genetics and Evolution(2019), https://doi.org/10.1016/ j.meegid.2020.104182

This is a PDF file of an article that has undergone enhancements after acceptance, such as the addition of a cover page and metadata, and formatting for readability, but it is not yet the definitive version of record. This version will undergo additional copyediting, typesetting and review before it is published in its final form, but we are providing this version to give early visibility of the article. Please note that, during the production process, errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.

© 2019 Published by Elsevier.

Journal Pre-proof

Molecular epidemiological study of multidrug-resistant tuberculosis isolated from sputum samples in Eastern Cape, South Africa

Nolwazi Londiwe Bhembe 1* and Ezekiel Green1

Department of Biotechnology and Food Technology, Faculty of Science, University of

f

1

pr

e-

*To whom correspondence should be addressed;

oo

Johannesburg, Doornfontein, 2028, South Africa

Dr Nolwazi Londiwe Bhembe, Department of Biotechnology and Food Technology, Faculty

Pr

of Science, University of Johannesburg, Doornfontein, 2028, South Africa.

al

Tel: +276032441675.

Jo u

rn

E-mail: [email protected]

1

Journal Pre-proof ABSTRACT Drug-resistant tuberculosis prevalence is still a global challenge. Making it imperative to examine the molecular epidemiology of drug resistant tuberculosis. Molecular epidemiology methods can evaluate transmission patterns and risk factors, ascertain transmission cases of multidrug-resistant tuberculosis (MDR-TB) and furthermore determine transmission patterns in a human populace. This work focuses on MDR-TB isolates in distinguishing them into several species and genotyping the MDR-TB isolates, mainly for epidemiological studies

oo

f

using the genomic regions of difference and the spoligotyping techniques. A total of 184 deoxyribonucleic acid isolated from sputum samples that showed resistance against the two

pr

major first-line anti-tuberculosis drugs (Rifampicin and Isoniazid) were examined. The

e-

deoxyribonucleic acid samples were amplified with primers specific for each flanking region

Pr

of the genomic regions of difference for the identification of different MTBC species. Isolates were further characterized into different lineages using the spoligotyping commercial kit. The

al

M. tuberculosis species was detected in 83.7% (154/184) of the deoxyribonucleic acid

rn

isolates, followed by the M. caprae in 8.7% (16/184) and the least detected species was the

Jo u

M. africanum in 2.2% (4/184). Nineteen spoligotype international types (SITs) were identified in this study. The pre-existing shared types were from 94.6% (174/184) isolates with 1.1% (2/184) isolates recognized as orphans and 4.3% (8/184) isolates were not found in the SITVIT database. The predominant family (spoligotype) was the Beijing with 67.4% (124/184) strains. This study gives a general overview of drug resistant strains and the circulating strains in the Eastern Cape, South Africa and it shows that the common Mycobacteria in the province is the Beijing strain.

Keywords: tuberculosis; resistant strains; genetic variability; spoligotyping; drug-resistance

2

Journal Pre-proof 1. Introduction Tuberculosis (TB) is still an encumbrance global and South Africa is one of the countries that has the highest burden (TB Statistics South Africa, 2019). Approximately 322,000 cases of active TB have been reported in South Africa in 2017 (World Health Organization, 2018a). A diminution was observed in reported TB cases, from 400,000 to 300,000 in 2014 (World Health Organization, 2015). However, there has been a noticeable increase of multidrugresistant TB (MDR-TB) leading to South Africa being part of the countries with the highest

oo

f

MDR-TB incidences globally (World Health Organization, 2015). The highest burden of rifampicin-resistant TB (RR-TB) has been reported in South Africa, as well as MDR-TB

pr

(World Health Organization, 2016). Multidrug-resistant TB is tuberculosis resistant to the

e-

two principal first line anti-TB drugs (isoniazid and rifampicin) (World Health Organization,

Pr

2018). Approximately 490,000 people in 2016 became ill with MDR-TB globally, and about 110,000 people were resistant to rifampicin (World Health Organization, 2017). Moreover, it

al

is anticipated that about 6.2% cases of the 110,000 reported as MDR-TB were XDR-TB cases

rn

(World Health Organization, 2018b).

Jo u

In 2014, about 35,000 multidrug-resistant TB cases were recorded in South Africa (Eyewitness News, 2017). Yet again, there were 19,000 incidences of MDR-TB reported in 2016 from South Africa (World Health Organization, 2017). Tuberculosis in general is still one of the leading, vague health problems in this country, including the Eastern Cape Province with the situation worsening irrespective of the implemented TB control programs such as the direct observed treatment strategy (DOTS) and the vaccine Bacille CalmetteGuerin (BCG) (Karim et al, 2009). The province has the fourth highest TB prevalence with approximately

662/100,000

population

infected

Organization, 2015). 3

with

TB

annually

(World

Health

Journal Pre-proof Different members of the Mycobacterium tuberculosis complex (MTBC) (Forrellad et al, 2013) instigate tuberculosis. The MTBC members are M. tuberculosis, M. africanum, M. bovis, M. bovis BCG vaccine strain, M. canettii, M. caprae, M. pinnipedii, M. microti and M. mungi infecting animals as well as humans (Smith et al, 2006; Young et al, 2008; Alexander et al, 2010). Of these entire complex members, the M. tuberculosis is the most infamous member of the MTBC, causing TB in humans, resulting in an annual report of about 1.5 million deaths (Raviglione et al, 2016). The MTBC encompasses of closely interrelated

oo

f

subspecies that have 99.9% genome sequence identity (Comas et al, 2013). The genomic region of difference is reported to be involved in host interaction and pathogenicity, this

e-

pr

determines the host spectrum of individual MTBC species (Ru et al, 2017). Laboratory diagnosis of human TB relies on acid-fast staining, culture and microscopy

Pr

mainly (Warren et al, 2006). These traditional diagnostic methods have limitations that include sensitivity, specificity and they are time consuming (Parkash et al, 2009). The latest

al

tool used in South Africa (GeneXpert/RIF) can diagnose TB from sputum specimens

rn

identifying the susceptibility or rifampicin resistance in two hours (Rachow et al, 2011).

Jo u

However, this tool cannot differentiate between the different MTBC member and it cannot delineate the strain lineages. This is where the genotyping techniques come into play because of their ability to map disseminating strains in a specific region. Studies on molecular epidemiology of MDR-TB are vital because they offer cutting-edge methods to study spreading dynamic forces and evolutionary genetics of the pathogen with supremacy on TB control actions (Kamerbeek et al, 1997). There are several molecular methods that are used to review the molecular epidemiology of TB. These methods include Insertion Sequence 6110 (IS6110) restriction fragment length polymorphism (RFLP) typing,

spoligotyping and

mycobacterial interspersed repetitive unit variable number of tandem DNA repeats (MIRU4

Journal Pre-proof VNTR) (Gupta et al, 2014). In this study there are two tools used, the genomic regions of difference and spoligotyping. We consider these tools imperative for accurate diagnosis of TB and public health surveillance, especially in the Eastern Cape Province, where data on molecular epidemiology of MDR-TB is still scanty. Materials and Methods 2.1.

Study setting and Ethics authorization

oo

f

The Eastern Cape Province has the third highest population, with approximately 7 million people (Census, 2011) and it has the second worst burden of about 49% of the population

pr

living below the poverty mark in rural areas (Bradshaw et al, 2000). This province is

e-

bordered by Lesotho and the Free State in the north, KwaZulu Natal in the northeast and the

Pr

Northern and Western Cape in the west (Bradshaw et al, 2000). The University of Fort Hare Research Ethics Committee (UREC) permitted this current study (issued ethical clearance

Sample collection

rn

2.2.

al

certificate: REC-270710-028-RA Level 01).

Jo u

Samples were collected as earlier described in our previous publication (Bhembe et al, 2014); three thousand eight hundred and ten sputum specimens were collected from patients with tuberculosis symptoms in different infirmaries and clinics in the Eastern Cape Province over a period of 24 months (January 2012 to December 2013). Succinctly, 4.8% (184/3810) DNA samples were used in this present study. A written informed consent form of participants was attained and the confidentiality of the patient’s identities was secured. The Biodata of the patients, such as the gender and age group (0-14, 15-29, 30-44, 45-59, 60 and above) was also collected. The patients were summarized according to age groups of 15 year intervals into 5 groups (0-14 young people, 15-29 youth, 30-44 adults, 45-59 middle age, 60 and above 5

Journal Pre-proof which is the elderly). This grouping was achieved through the use of pivot table by age range, using the Microsoft office excel to perform the analysis and conception of data and information. The samples were transported to the Microbiology laboratory at the University of Fort Hare for further analyses. The personal identifiers of all patients used in this study are imperceptible to safeguard the identity of the patients. 2.3.

Mycobacterial culture

f

In the Molecular Pathogenicity and Molecular Epidemiology Research Group (MPMERG)

oo

laboratory, samples were decontaminated and inoculated on Lowenstein-Jensen (LJ) slants

pr

that were incubated for approximately eight weeks at 37°C. Mycobacteria isolation and

e-

identification for each isolate was carried out by the hospitals and clinics using the acid-fast

2.4.

Pr

staining (AFB) technique.

Antibiotic susceptibility testing

al

An antibiotic susceptibility test was executed in the MPMERG laboratory at the University of

rn

Johannesburg. Briefly, Mycobacterium colony growth from the LJ slants was transferred into

Jo u

0.85% normal saline that was in screw-cap tubes. This method was executed using the agardilution proportion method as outlined by Balows et al. (1991) for isoniazid and rifampicin with some modifications as mentioned earlier in our published work (Bhembe et al, 2014). With all sets of experimentations, the M. tuberculosis (H37Rv) strain was used as a positive control. 2.5.

Deoxyribonucleic acid (DNA) isolation

To isolate DNA, the AFB positive isolates were heat-killed by incubating the cells in a water bath for an hour at the temperature of 80o C following the protocol defined by Berg et al.

6

Journal Pre-proof (2009). A volume of 200 µl fractions of the supernatant was aliquated into sterile 1.5 ml Eppendorf tube for each isolate and kept in a -80o C freezer. 2.6.

Molecular detection of MDR-TB

One hundred and eighty-four (184) DNA isolates were confirmed to be MTBC members with the Seeplex® MTB Nested Ace detection assay (Seegene Inc, Korea) following the manufacturer’s instructions. The constituted polymerase chain reaction (PCR) mixture was

f

amplified using the MyCyclerT M thermal cycle (Bio-Rad, Cape Town, South Africa) and the

oo

amplification products were viewed with the Alliance 4.7 transilluminator (UVITEC Limited,

pr

Cambridge, UK). The DNA isolates were further amplified through a multiplex PCR with

e-

four primer sets per reaction. For this study, 3 primers per primer set (RD1, RD4, RD9 and RD12) as described by Warren et al. (2006) with some modifications were used. Briefly, each

Pr

PCR reaction contained 2 μl DNA template, 5 μl Q-buffer, 2.5 μl of 10X buffer, 2 μl of 25

al

mM MgCl2 , 4 μl of the 10 mM dNTPs, 0.5 μl of each primer (50 pmol/μl) 0.125 μl HotStarTaq DNA polymerase (Fermentas, Thermo Scientific, US) and sterile nuclease free

rn

H2 O was added to a total volume of 25 μl. The H37Rv (ATCC 27294), strain was used as a

Jo u

positive control and sterile H2 O was used as a negative control. Amplification was attained using the MyCyclerT M thermal cycler (BioRad, Cape Town, South Africa). The cycling conditions comprised of 1 cycle at 95o C for 10 minutes (Taq activation cycle), followed by 1 cycle denaturation at 95o C for 15 minutes, 45 cycles of denaturation at 94 o C for 1 minute, primer annealing at 62o C for 1 minute and extension at 72o C for 1 minute followed by 1 cycle of the final extension at 72o C for 10 minutes. The amplified products were electrophoretically segregated on 3.0% agarose gel containing 5 μl Ethidium bromide at 70V for 4 hours. The size of amplicons was estimated by comparison with a 100 base pair (bp) DNA ladder (Thermo Scientific). Images were visualized with a 4.7 7

Journal Pre-proof Alliance

transilluminator

(UVITEC

Limited,

Cambridge,

UK) and

the

results were

interpreted as described by Warren et al. (2006). 2.7.

Sub-lineage classification of MDR-TB isolates and data analysis

The isolates were further characterized into different strains using the spoligotyping commercial kit (Isogene Life Science B.V., Utrecht, The Netherlands) following the protocol as described by Kamerbeek et al. (1997). The genotypes were articulated in binary code and

the MIRU-VNTRplus (http://www.miru-vntrplus.org) (Weniger et al, 2010) and

oo

into

f

octal formats in Microsoft Excel spreadsheets and all the spoligotyping data was deposited

pr

SpolDB4.0 database (Brudey et al, 2006). The SPOTCLUST (a combination of models) by

Pr

their spacer oligonucleotide typing patterns. 3. Results

Isolation, identification and susceptibility profiles

al

3.1.

e-

Vitol et al. (2006) was used to determine the strain families of the MDR-TB isolates based on

rn

The 4.8% (184/3810) of the sputum samples investigated in this present study was resistant to

Jo u

both the two major first line anti-TB drugs (rifampicin and isoniazid). Multidrug-resistant TB was identified in all 184 samples that were targeted during sample collection. This study explicates the different MTBC members and the lineages detected in different age groups of the study population. We observed that the M. tuberculosis was the most detected species with most patients within the age group range of 15 to 29 years, 35.3% (65/184) followed by those in the age group 30 to 44 years 31.5% (58/184) (Table 1). 3.2.

Differentiation of MDR-TB isolates using RDs

The different amplification products matched the expected base pairs, which designate the presence and absence of the genomic regions of difference (RD1, RD4, RD9 and RD12). It 8

Journal Pre-proof was possible to differentiate between Mycobacterium tuberculosis, which showed four different band sizes (146 bp, 172 bp, 235 bp and 369 bp) and the M. canettii species showed three different band sizes (146 bp, 172 bp and 235 bp). However, it was again possible to differentiate between the M. caprae, which showed two bands (146 bp and 172 bp) and M. africanum species that showed three bands (146 bp, 172 bp and 369 bp) using the primer sets reported by Warren et al. (2006). However, these species were not the only ones detected;

Delineation of MTBC isolates

oo

3.3.

f

the M. bovis species was also identified in this study as shown in Table 2.

pr

The results show the possibility of differentiating between the different Mycobacterium

e-

tuberculosis complex members M. tuberculosis, M. caprae, M. bovis, M. africanum and M. canettii. The M. tuberculosis species were detected in 83.7% (154/184) of the DNA isolates

Pr

followed by M. caprae in 8.7% (16/184) and the least detected species was M. africanum in

al

2.2% (4/184) (Table 2). Nevertheless, in this study, we were able to detect species that are known to be primarily hosted by humans (M. tuberculosis, M. africanum and M. canettii) and

rn

species known to be primarily hosted by animals such as goats (M. caprae) and cattle (M.

3.4.

Jo u

bovis) from human sputum specimens (Table 2). Distribution of the different genotypes

Nineteen spoligotype international types (SITs) were identified in this study. The pre-existing shared types were from 94.6% (174/184) isolates, with 1.1% (2/184) isolates recognized as orphans and 4.3% (8/184) isolates were not found in the SITVIT database. However, the TBinsight SPOTCLUST (SpolDB3-based model) was able to identify the families of the isolates that were not found in the database as well as the orphans. A total of 17 families were identified: with the Beijing 67.4% (124/184) as the most prevalent family and the S, BOV3

9

Journal Pre-proof and AFRI were the least identified families with only one sample in each family 0.5% (1/184) (Table 3). 4. Discussion The epidemiology of TB has indicated that TB spread from one person/ animal to the next through the air by nuclei droplets produced when the infected person with either pulmonary or laryngeal TB, sneezes or coughs (National tuberculosis Management Guidelines 2014).

f

This mode of transmission is more prevalent in a certain social state of affairs (for instance,

oo

food insecurity and malnutrition, geographic and cultural hindrances to health maintenance

pr

access and underprivileged housing and environmental conditions) present in South Africa

e-

(Hargreaves et al. 2011).

investigation

and

control

Pr

Information on molecular epidemiology of MDR-TB is valuable as it supports TB program

of

countries.

However,

MDR-TB

molecular

al

epidemiological data of the Eastern Cape Province is still meagre. For better understanding of

rn

the moving and intensifying clusters of MDR-TB within the province, we have characterized

Jo u

the mycobacteria isolated from TB patients using the genomic region of difference and the spoligotyping tools. In this study, the circulating strains in the province were known and their drug sensitivity patterns defined.

The genomic region of difference using RD1, RD4, RD9 and RD12 designated the studied strains as M. tuberculosis, M. bovis, M. caprae, M. canettii and M. africanum (Table1). The M. tuberculosis species proved to be the most dominant that infected most of the TB patients. However, most of the patients infected with MDR-TB were between the age groups of 30 to 44 years 42.4% (78/184) and 15 to 29 years 38% (70/184). This could be due to the economical activeness of people at this ages. Yet again, in this study more females than males are reported to be infected with M. tuberculosis, which is different from what has been 10

Journal Pre-proof reported by other studies, a higher male ratio as compared to females infected with M. tuberculosis (Murray, 1991; Hudelson, 1996). However, this could also be due to selection bias where more female samples than those of males as reported before (Bhembe et al, 2014) were collected from the clinics and hospitals. We noted that species known to be primary hosted by animals were also detected from the samples gathered from human sputum. This shows the possibility of transmission of MTBC

f

from animals to humans and vice versa. Nevertheless, this is not easy to learn using the

oo

spoligotyping tool; so, it needs further discriminatory techniques, which will be utilized for

pr

future studies on these isolates. Tuberculosis is endemic in South Africa and identifying the

e-

circulating strains within populations is essential and is important for understanding the epidemiology of the infection in a region, along with supporting the initiation of better

Pr

control programs that will avoid the expansion of MDR-TB at both local and global level. In 2016, the World Health Organization (World Health Organization, 2016) reported about

al

19,613 patients detected with rifampicin-resistant TB in South Africa. This number indicated

rn

a substantial growth from the 2011 reported prevalence of about 13,000 rifampicin-resistant

Jo u

cases in this country (World Health Organization, 2016). It is not surprising to find 83.7% of M. tuberculosis in this study samples because this species is primarily hosted by humans. However, the prevalence of TB by this strain proves to be high, which is a problem in the public health domain. There are several reports in South Africa in addition to other countries of the M. tuberculosis Beijing strain reported amongst cases of drug-resistant MDR-TB (Glynn et al, 2002; Stavrum et al, 2009). Previous studies reported in the Eastern and Western Cape Provinces postulate that there is an association of MDR-TB and the Beijing family strains (Bifani et al, 2002; Johnson et al, 2006; Cliota et al,

11

Journal Pre-proof 2012; Klopper et al, 2013). However, we cannot prove the association of the Beijing strain and MDR-TB. Results from the SITVIT database and spoligotyping international typing (SIT) number indicates that the Beijing strain is the most prevalent shared type in this study. The strain diversity detected is generally constant with studies reported from other parts of South Africa (Chiota et al, 2007; Mlambeo et al, 2008; Marais et al, 2013). The LAM genotype is the

f

second predominant strain in this study. This is not startling because this genotype has been

oo

reported to be prevalent in the Eastern and Western Cape Provinces of South Africa (Stavrum

pr

et al, 2009), which could be one of the reasons of the prevalence of this genotype reported in

e-

this study. This could suggest a continuing spread of this strain in the province.

Pr

The X lineage was also detected in this study (Table 3); this lineage is prevalent in the Northern and Western Cape Provinces (Stavrum et al, 2009). This could mean the strain has

al

been transmitted from these provinces, because these provinces borders the Eastern Cape

rn

Province; the Western Cape Province border this province from the west and the Northern

Jo u

Cape Province border this province in the north-west of the escarpment. This study reports 2.7% (5/184) of the BOV (M. bovis subsp. bovis respectively). The primary host of M. bovis is cattle, but a variety of domestic animals, wild animals along with humans may possibly be infected with this species (Ayele et al, 2004; Gibson et al, 2004). Most human cases of bovine TB are observed in countries where bovine TB is not controlled, mostly affecting young people subsequently from either drinking or handling contaminated milk (Cosivi et al, 1998). Data on human TB instigated by the M. bovis species in both industrialized and unindustrialized countries is deficient. A study by Elsabban et al. (1992) showed the 12

Journal Pre-proof proportion of sputum-positive TB patients infected with M. bovis to be 0.4%, 5.4% and 6.4%, while Nafeh et al. (1992) reported that nine out of twenty indiscriminately chosen patients with TB peritonitis were infected with M. bovis and the rest with M. tuberculosis. The M. bovis subsp. M. caprae species was also found in 1.1% (2/184) of the isolates. The species is a causative agent of TB in caprine herds (Duet et al, 2008). This is not surprising as the Eastern Cape Province has the largest numbers of goats in South Africa (Statistics South

f

Africa, 2014). This pathogen was isolated from cattle in Central Europe (Pavlik et al, 2002)

oo

and pigs (Kongpetchsatit et al, 2006). There are less than 1% reported TB cases caused by

pr

this species in industrialized countries, but other areas or population might experience a

e-

higher prevalence (Amafu 2006; Rodwell et al, 2008). A study by Bhembe et al. (2017) also reported on this species isolated from slaughtered cattle lymph nodes in the Eastern Cape

Pr

Province. These results may possibly suggest transmission of this pathogen from livestock to humans. The connection between humans and animals has indicated that the mode of

al

contagion, such as the use of unpasteurized milk and the proximity between mammals occurs

rn

between man and animals (McDaniel et al, 2014).

Jo u

This study also reports on the AFRI lineage 2.2% (4/184) belonging to the M. africanum species which was detected from the isolates. This lineage has been reported as a prevailing M. bovis clonal complex (Egbe et al, 2017). The M. africanum (subtype I) has traits that resemble those of M. bovis and is widespread in West Africa (Mostowy et al, 2004). The possible route of transmission of M. africanum may perhaps be through the movement of persons from one region to another. Family 33 is also reported in this study with the prevalence of 2.2% (4/184) reported to be belonging to the Euro-American lineage (MagugaPhasha et al, 2017). It is not the first time this family is reported in South Africa, it has been

13

Journal Pre-proof reported in the Limpopo Province (Maguga-Phasha et al, 2017) and the Eastern Cape Province (Bhembe et al, 2017). The T family has also been reported in different provinces of South Africa, which includes the Free State (van Dijk et al, 2016) and the Eastern Cape (Stavrum et al, 2009); which was described as the most predominant family in the province. This present study reports a prevalence of 1.1% (2/184) of the T family. This could mean there is continuing transmission

f

of this strain, although it is at a marginal rate. The least detected family in this study is the S

oo

family 0.5% (1/184); this is not surprising, because the family has been reported with limited

pr

dissemination through African countries (Mbugi et al, 2016).

e-

The SpolDB4 database also identified the Manu 1.1% (2/184) and Canettii 2.7% (5/184)

Pr

families. The Canettii family is associated with the M. canettii species. The M. canettii species is a sporadic, smooth, and a variant of the M. tuberculosis commonly isolated in

al

Africa from humans (Fabre et al, 2010). The species has been reported worldwide and the

Jo u

5. Conclusions

rn

majority of reported cases are from patients who have either travelled or lived in Africa.

This study validates that the Beijing strain is the prevalent strain disseminating in the Eastern Cape Province. The results of this study denote a reference point of the population structure in the province and may serve, as a drive for future molecular studies. Genotyping of multidrug-resistant Mycobacterium tuberculosis strains is highly imperative for multidrugresistant tuberculosis control as it allows the detection of alleged occurrences and transmission chains. This study evidences that multidrug-resistant TB is instigated by diverse species of Mycobacterium causing the infection in humans, which together form the Mycobacterium tuberculosis Complex. 14

Journal Pre-proof Ethical considerations The experiment was designed and conducted in agreement with the ethical guidelines highlighted by the National Committee for Research Ethics in Science and Technology of South Africa and all sampling and analytical procedures were followed as documented in the University of Fort Hare’s health and safety guidelines. Acknowledgements

f

The authors extend their gratitude to the National Research Foundation (NRF) of South

oo

Africa (grant number: 76255) for a subvention, the University of Fort Hare and the

pr

University of Johannesburg for letting us use their amenities.

e-

Conflict of interest

Pr

No conflict of interest declared. References

al

1. Alexander, K.A., Laver, P.N., Michel, A.L., Williams, M., van Helden, P.D., Gey van

rn

Pittius, N.C., 2010. Novel Mycobacterium tuberculosis complex pathogen, M. mungi. Emerging Infectious Diseases 16, 1296-1299.

Jo u

2. Amanfu, W., 2006. The situation of tuberculosis and tuberculosis control in animals of economic interest. Tuberculosis 86, 330-335. 3. Ayele, W.Y., Neill, S.D., Zinsstag, J., Weiss, M.G., Pavlik, I., 2004. Bovine tuberculosis: an old disease but a new threat to Africa. International Journal of Tuberculosis and Lung Diseases 8, 924-937. 4. Balows, A., Hausler, W.J., Herrmann Jr., K.L., Isenberg, H.D., Shadomy, H.J., 1991. Manual of Clinical Microbiology, fifth ed. Washington, DC: American Society for Microbiology. 5. Berg, S., Firdessa, R., Habtamu, M., Gadisa, E., Mengistu, A., Yamuah, L., Ameni, G., Vordermeier, M., Robertson, B.D., Smith, N.H., Engers, H., Young, D., 15

Journal Pre-proof Hewinson, G.R., Aseffa, A., Gordon, S.V., (2009). The burden of Mycobacterial disease in Ethiopian cattle: Implications for public health. PLoS ONE 4(4), e5068. 6. Bhembe, N.L., Jaja, I.F., Nwodo, U.U., Okoh, A.I., Green, E., 2017. Prevalence of tuberculous lymphadenitis in slaughtered cattle in Eastern Cape, South Africa. International Journal Infectious Diseases 61(2017), 27-37. 7. Bhembe, N.L., Nwodo, U.U., Govender, S., Hayes, C., Ndip, R.N., Okoh, A.I., Green E.,

2014.

Molecular

detection

and

characterization

of

resistant

genes

in

Mycobacterium tuberculosis complex from DNA isolated from tuberculosis patients

oo

f

in the Eastern Cape Province South Africa. BMC Infectious Diseases 4(14), 479. Doi: 10.1186/1471-2334-14-479. P.J.,

Mathema,

B.,

Kurepina,

N.E., Kreiswirth, B.N., 2002. Global

pr

8. Bifani,

dissemination of the Mycobacterium tuberculosis W-Beijing family strains. Trends in

e-

Microbiology 10, 45-52.

S.A., Allix,

C., Aristimuño, A., Cheong,

L., Arora,

S., Diel,

J., Baumanis,

R., Ellermeier,

al

P., Cataldi,

Pr

9. Brudey, K., Driscoll, J.R., Rigouts, L., Prodinger, W.M., Gori, A., Al-Hajoj,

M., Ferdinand,

S., Garcia

de

Viedma,

V., Binder,

C., Evans,

D., Garzelli,

L., Cafrune,

J.T., Fauville-Dufaux,

C., Gazzola,

L., Gomes,

rn

H.M., Guttierez, M.C., Hawkey, P.M., van Helden, P.D., Kadival, G.V., Kreiswirth, B.N., Kremer,

K., Kubin,

M., Kulkarni,

S.P., Liens,

B., Lillebaek,

T., Ho,

Jo u

M.L., Martin, C., Martin, C., Mokrousov, I., Narvskaïa, O., Ngeow, Y.F., Naumann, L., Niemann,

S., Parwati,

V., Rasolonavalona, S., Shemyakin, D., Streicher,

T., Rossetti,

I.G., Singh,

E.M., Suffys,

I., Rahim,

Z., Rasolofo-Razanamparany,

M.L., Rüsch-Gerdes,

U.B., Somoskovi, P.N., Tortoli,

S., Sajduda,

A., Skuce,

E., Tracevska,

A., Samper,

R.A., van

Soolingen,

T., Vincent,

V., Victor,

T.C., Warren, R.M., Yap, S.F., Zaman, K., Portaels, F., Rastogi, N., Sola, C., 2006. Mycobacterium

tuberculosis

complex

genetic

diversity:

mining

the

fourth

international spoligotyping database (SpolDB4) for classification, population genetics and epidemiology. BMC Microbiology 6, 23.

16

Journal Pre-proof 10. Census,

2011.

Provincial

profile:

Eastern

Cape,

2011.

http://www.statssa.gov.za/publications/Report-03-01-71/Report-03-01-712011.pdf (accessed 23 February 2018). 11. Chihota, V., Apers, L., Mungofa, S., Kasongo, W., Nyoni, I.M., Tembwe, R., Mbulo, G., Tembo, M., Streicher, E.M., van der Spuy, G.D., Victor, T.C., van Helden, P., Warren, R.M., 2007. Predominance of a single genotype of Mycobacterium tuberculosis in regions of southern Africa. International Journal of Tuberculosis and Lung Diseases 11, 311-318.

oo

f

12. Chihota, V.N., Muller, B., Mlambo, C.K., Pillay, M., Tait, M., Streicher, E.M., Marais, E., van der Spuy, G.D., Hanekom, M., Coetzee, G., Trollip, A., Hayes, C.,

pr

Bosman, M.E., Gey van Pittius, N.C., Victor, T.C., van Helden, P.D., Warren, R.M., 2012. Population structure of multi- and extensively drug-resistant Mycobacterium

e-

tuberculosis strains in South Africa. Journal of Clinical Microbiology 50, 995-1002.

Pr

13. Comas, I., Coscolla, M., Luo, T., Borell, S., Holt, K.E., Kato-Maeda, M., 2013. Out of Africa migration and Neolithic coexpansion of Mycobacterium tuberculosis with

al

modern humans. Nature Genetics 45, 1176-1182. 10.1038/ng.2744. 14. Cosivi, O., Grange, J.M., Daborn, C.J., Raviglione, M.C., Fujikura, T., Cousins, D.,

rn

Robinson, R.A., Huchzermeyer, H.F., de Kantor, I., Meslin, F.X., 1998. Zoonotic Mycobacterium

bovis

in developing countries. Emerging

Jo u

tuberculosis due to

Infectious Diseases 4, 59-70. 15. Duet, E.L., Domingos, M., Amado, A., Botelho, A., 2008. Spoligotyping diversity of Mycobacterium caprae animal isolates. Veterinary Microbiology 130, 415-421. 16. Egbe, N.F., Muwonge, A., Ndip, L., Kelly R.F., Sander, M.., Tanya, V., Ngu Ngwa, V., Handel1, I.G., Novak, A., Ngandalo, R., Mazeri, S., Morgan K.L., Asuquo A., Bronsvoort, B.M., 2017. Molecular epidemiology of Mycobacterium bovis in Cameroon, Scientific Reports - Nature 7, 4652. Doi:10.1038/s41598-017-04230-6. 17. Elsabban, M.S., Lofty, O., Awad, W.M., Soufi, H.S., Mikhail, D.G., Hammam, H., 1992. Bovine tuberculosis and its extent of spread as a source of infection to man and animals in Arab Republic of Egypt. In: Proceedings of the International Union 17

Journal Pre-proof Against Tuberculosis and Lung Disease Conference on Animal Tuberculosis in Africa and the Middle East. Cairo, Egypt. 18. Eyewitness News, 2017. Concern over sharp increase in drug-resistant TB cases in South

Africa.

http://ewn.co.za/2017/02/06/concern-over-sharp-increase-in-drug-

resistant-tb-cases-in-sa (accessed 16 December 2018). 19. Fabre, M., Hauck, Y., Soler, C., Koeck, J.L., Van Ingen, J., van Sloongen, D., 2010. Molecular characteristics of “Mycobacterium canettii’’ the smooth Mycobacterium

f

tuberculosis bacilli. Journal of Infection, Genetics and Evolution 10, 1165-1173.

oo

20. Forrellad, M.A., Klepp, L.I., Gioffré, A., Sabio, J., García, Y., Morbidoni, H.R., de la Paz Santangelo, M., Cataldi, A.A., Bigi, F., 2013. Virulence factors of the

pr

Mycobacterium tuberculosis complex. Virulence 4(1), 3-66.

e-

21. Gibson, G.R., Probert, H.M., Van loo, J., Rastall, R.A., Roberfroid, M.B., 2004. Dietary modulation of the human colonic microbiota: updating the concept of

Pr

prebiotics. Nutrition Research Reviews 17, 259-275. 22. Glynn, J.R., Whiteley, J., Bifani, P.J., Kremer, K., van Soolingen, D., 2002.

al

Worldwide occurrence of Beijing/W strains of Mycobacterium tuberculosis: a

rn

systematic review. Emerging Infectious Diseases 8, 843-849. 23. Gupta, A., Kulkarni, S., Rastogi, N., Anupurba, S., 2014. A study of Mycobacterium

Jo u

tuberculosis genotypic diversity & drug resistance mutations in Varanasi, north India. Indian Journal of Medical Research 139(6), 892-902. 24. Hargreaves, J.R., Boccia, D.,

Evans, C.A., Adato, M., Petticrew, M., Porter, J.D.H.,

2011. The Social Determinants of Tuberculosis: From Evidence to Action. American Journal of Public Health 101(4), 654-662. 25. Hudelson, R., 1996. Gender differentials in tuberculosis: the role of socio-economic and cultural factors. Tuberculosis and Lung Diseases 77, 391-400. 26. Johnson, R., Warren, R., Strauss, O.J., Jordaan, A.M., Falmer, A.A., Beyers, N., Schaaf, H.S., Murray, M., Cloete, K., van Helden, P.D., Victor, T.C., 2006. An outbreak of drug-resistant tuberculosis caused by a Beijing strain in the Western 18

Journal Pre-proof Cape, South Africa. International Journal of Tuberculosis and Lung Diseases 10(12), 1412-1414. 27. Kamerbeek, J., Schouls, L., Kolk, A., van Agterveld, M., van Soolingen, D., Kuijper, S., 1997. Simultaneous detection and strain differentiation of Mycobacterium tuberculosis for diagnosis and epidemiology. Journal of Clinical Microbiology. 35, 907-14. 28. Karim, S.A., Churchyard, G., Karim, Q.A., Lawn, S.D., 2009. HIV infection and tuberculosis in South Africa: an urgent need to escalate the public health response.

oo

f

Lancet 374, 921-933.

29. Klopper, M., Warren, R.M., Hayes, C., Gey van Pittius, N.C., Streicher, E.M., Muller,

pr

B., Sirgel, F.A., Chabula-Nxiweni, M., Hoosain, E., Coetzee, G., David van Helden, P., Victor, T.C., Trollip, A.P., 2013. Emergence and spread of extensively and totally

e-

drug-resistant tuberculosis, South Africa. Emerging Infectious Diseases Journal 19,

30. Kongpetchsatit,

O.,

Pr

449-455.

Phatihattakorn,

W.,

Mahakunkijcharoen,

Y.,

Eampokalarp,

al

Boonyasopun, J., Ramasoota, P., 2006. Mutation in the rpoB gene of the rifampicin resistant M. avium Complex strains from Thailand. Southeast Asian Journal of

rn

Tropical Medicine and Public Health 37(suppl 3), 165-173.

Jo u

31. Maguga-Phasha, N.T.C., Munyai, N.S., Mashinya, F., Makgatho, M.E., Mbajiorgu, E.F., Genetic diversity and distribution of Mycobacterium tuberculosis genotypes in Limpopo,

South

Africa.

BMC

Infectious

Diseases

(2017)

17:764.

Doi:

10.1186/s12879-017-2881-z. 32. Marais, B.J., Mlambo, C.K., Rastogi, N., Zozio, T., Duse, A.G., Victor, T.C., Marais, E., Warren, R.M., 2013. Epidemic Spread of Multidrug-Resistant Tuberculosis in Johannesburg. South Africa Journal of Clinical Microbiology 51(6), 1818-1825; Doi: 10.1128/JCM.00200-13. 33. Mbugi, E.V., Katale, B.Z., Streicher, E.M., Keyyu, J.D., Kendall, S.L., Dockrell, H.M., Michel, A.L., Rweyemamu, M.M., Warren, R.M., Matee, M.I., van Helden, P.D., Couvin D., Rastogi, N., 2016. Mapping of Mycobacterium tuberculosis complex 19

Journal Pre-proof genetic diversity profiles in Tanzania and Other African Countries. PLoS ONE 11(5), e0154571. Doi:10.1371/ journal.pone.0154571. 34. McDaniel, C.J., Cardwell, D.M., Moeller Jr, R.B., Gray, G.C., 2014. Humans and Cattle: A Review of Bovine Zoonoses. Vector Borne Zoonotic Diseases 14(1): 1–19. Doi: 10.1089/vbz.2012.1164. 35. Mlambo, C,K., Warren, R.M., Poswa, X., Victor, T.C., Duse, A.G., Marais, E., 2008. Genotypic diversity of extensively drug-resistant tuberculosis (XDR-TB) in South

f

Africa. International Journal of Tuberculosis and Lung Diseases 12, 99-104.

oo

36. Mostowy, S., Onipede, A., Gagneux, S., Niemann, S., Kremer, K., Desmond, E.P., Behr, M., 2004. Genomic analysis distinguishes Mycobacterium africanum. Journal of

pr

Clinical Microbiology 42(8), 3594-3599.

e-

37. Murray, C.J., 1991. Social, economic and operational research on tuberculosis: recent studies and some priority questions. Bulletin of the International Union against

Pr

Tuberculosis and Lung Disease 66, 149-156.

38. Nafeh, M.A., Medhat, A., Abdul-Hameed, A.G., Ahmad, Y.A., Rashwan, N.M.,

al

Strickland, G.T., 1992. Tuberculous peritonitis in Egypt: the value of laparoscopy in

rn

diagnosis. American Journal of Tropical Medicine and Hygiene 47, 470-477.

National

Jo u

39. National Strategic Plan on HIV, STIs and TB 2012-2016, 2017. South African AIDS

Council,

2017

www.gov.za/documents/national-strategic-plan-hiv-stis-and-tb-2012-2016

(accessed

26 September 2019).

40. National Tuberculosis Management Guidelines, 2014. http://www.tbonline.info/media/uploads/documents/ntcp_adult_tb- guidelines27.5.2014.pdf (accessed 23 April 2018). 41. Parkash, O., Singh B.P., Pai, M., 2009. Regions of Differences Encoded Antigens as Targets for Immuno diagnosis of Tuberculosis in Humans.

Scandinavian Journal of

Immunology 70(4), 345-357. 42. Pavlik, I., Dvorska, L., Bartos, M., Parmova, I., Melicharek, I., Jesenska, A., Havelkova, M., Slosarek, M., Putova, I., Martin, G., Erler, W., Kremer, K., van 20

Journal Pre-proof Soolingen, D., 2002. Molecular epidemiology of bovine tuberculosis Czech Republic and Slovakia in the period 1965-2001 studied by Spoligotyping. Veterinary Medicine 47, 181-194. 43. Rachow A, Zumla A, Heinrich N, Rojas-Ponce G, Mtafya B, Reither K, Ntinginya, E.N., O'Grady, J., Huggett, J., Dheda, K., Boehme, C., Perkins, M., Saathoff, E., Hoelscher, M., 2011. Rapid and accurate detection of Mycobacterium tuberculosis in sputum samples by Cepheid Xpert MTB/RIF assay a clinical validation study. PLoS One 6(6), 20458. Doi: 10.1371/journal.pone.0020458.

oo

f

44. Raviglione, M., Sulis, G., 2016. Tuberculosis 2015: burden, challenges and strategy for control and elimination. Infectious Disease Reports 8, 6570.

pr

45. Rodwell, T.C., Moore, M., Moser, K.S., Brodine, S.K., Strathdee, S.A., 2008. Tuberculosis from Mycobacterium bovis in binational communities, United States.

e-

Emerging Infectious Diseases 14, 909-916.

Pr

46. Ru, H., Liu, X., Lin, C., Yang J., Chen, F., Sun, R., Zhang, l., Liu, J., 2017. The impact of genome region of difference 4 (RD4) of Mycobacterial virulence and BCG Frontiers

cell

infectious

microbiology

7,

239.

Doi

al

efficacy.

rn

10.3389/fcimb.201700239.

47. Smith, N.H., Gordon, S.V., de la Rua-Domenech, R., Clifton-Hadley, R.S., Hewinson,

Jo u

R.G., 2006. Bottlenecks and broomsticks: the molecular evolution of Mycobacterium bovis. Nature Review Microbiology 4, 670-681. 48. Statistics

South

Africa,

2014.

http://beta2.statssa.gov.za/?p=1447

(Accessed

1

February 2014).

49. Stavrum, R., Mphahlele, M., Ovreås, K., Muthivhi, P.T., Fourie, B., Weyer, K., Grewal, H.M., 2009. High diversity of Mycobacterium tuberculosis genotypes in South Africa and preponderance of mixed infections among ST53 isolates. Journal of Clinical Microbiology 47(6), 1848-1856. 50. TB

Statistics

South

Africa,

2019.

National,

incidence

and

https://tbfacts.org/tb-statistics-south-africa/ (accessed 29 September 2019).

21

provincial

Journal Pre-proof 51. van Dijk, A.V.S., Makhoahle, P.M., Rigouts, L., Kamaldeen Baba, K., 2016. Diverse molecular genotypes of Mycobacterium tuberculosis complex isolates circulating in the Free State, South Africa. International Journal of Microbiology 2016, 6572165, 7. 52. Vitol, I., Driscoll, J., Kreiswirth, B., Kurepina, N., Bennett, K.P., 2006. Identifying Mycobacterium tuberculosis Complex strain families using spoligotypes. Infections, Genetics and Evolution 6(6), 491-504. 53. Warren, R.M., Gey van Pittius, N.C., Barnard, M., Hesseling, A., Engelke, E., de Kock, M., Gutierrez, M.C., Chege, G.K., Victor, T.C., Hoal, E.G., van Helden, P.D.,

oo

f

2006. Differentiation of Mycobacterium tuberculosis complex by PCR amplification of genomic regions of difference. International Journal of Tuberculosis and Lung

pr

Diseases 10, 818-822.

54. Weniger, T., Krawczyk, J., Supply, P., Niemann, S., Harmsen, D., 2010. web

tool for polyphasic genotyping of Mycobacterium

e-

MIRUVNTRplus: a

Pr

tuberculosis complex bacteria. Nucleic Acids Research 38, W326-331. 55. World Health Organization, 2018a. What is multidrug-resistant tuberculosis (MDR-

56. World

Health

rn

September 2019).

al

TB) and how do we control it? http://www.who.int/features/qa/79/en/ (accessed 23

Organization,

2018b.

Global

Tuberculosis

Report.

Jo u

www.who.int/tb/publications/global_report/en/ (accessed 26 September 2019). 57. World Health Organization, 2016. Global Tuberculosis Report. Contract No.: WHO/HTM/

TB/2016.13.

https://apps.who.int/medicinedocs/en/d/Js23098en/

(accessed 21 March 2018). 58. World

Health

Organization,

2015.

Global

Tuberculosis

Report.

http://apps.who.int/iris/bitstream/10665/191102/1/9789241565059_eng.pdf?ua=1 (accessed 12 February 2018). 59. Young, D.B., Perkins, M.D., Duncan, K., Barry, C.E., 2008. Confronting the scientific obstacles to global control of tuberculosis. Journal of Clinical Investigations 118, 1255-1265.

22

Journal Pre-proof

Tables Table 1 Demographic profile of the study population defined using genomic regions of difference Age (years)

Gender of patients with species detected M. bovis

M. caprae

M. canettii

M. africanum

Male

Female

Male

Female

Male

Total

1

0

0

0

0

5

0

pr

1

1

0

0

70

e-

Male

Female

Male

Female

0-14

3

1

0

0

0

15-29

38

27

0

0

3

30-44

36

22

0

3

6

5

2

1

2

1

78

45-59

6

17

1

1

1

0

0

0

1

0

27

60+

3

1

0

0

0

0

0

0

0

4

Total

86

68

1

10

6

3

2

3

1

184

Pr

Female

oo

f

M. tuberculosis

al

0

Jo u

rn

4

23

Journal Pre-proof Table 2 Segregation of MTBC into species as identified following genomic regions of difference Number of DNA

isolated from humans

isolates

Humans

154

83.7%

M. caprae

Goats

16

8.7%

M. bovis

Cattle

5

2.7%

M. africanum

Humans

M. canettii

Humans

oo

M. tuberculosis

f

Primary host, but

Jo u

rn

al

Pr

e-

Total

24

Percentage

4

2.2%

5

2.7%

184

100%

pr

Species identified

Journal Pre-proof Table 3 Prevalence of the two major first line anti-TB drugs of the detected genotypes in the Eastern Cape Province of South Africa Family

East-Asian Euro-American Indo-Oceanic Euro-American Euro-American Euro-American Euro-American Euro-American Euro-American Indo-Oceanic Euro-American Euro-American Euro-American Euro-American M. bovis M. bovis M. africanum M. bovis M. canettii Euro-American Euro-American Euro-American Euro-American Euro-American Euro-American Euro-American Euro-American Euro-American Euro-American

Beijing LAM4 X1 X3 LAM3 T1 S LAM9 LAM4 MANU2 X3 X3 X3-variant1 LAM2 BOV1 BOV3 AFRI Bov_4 Caprea Canettii X3 LAM9 Family33 T1 Family33 X3 Family33 Family33 LAM3 LAM9

e-

pr

oo

f

SpolDB4 ID

rn

al

Pr

Octal code (spoligotype description) 000000000003771 777777607760731 777777377670371 300076777760771 776177607760771 777777777760771 776377777760771 777777607760771 777767607760731 777777607763771 700076777760771 700076777760671 700036777760771 677717607760771 676773777777600 640013777777600 770777777777671 277773777377600 000000000101000 500046637740661 777567607760771 777357757763771 557347637740661 557347607743661 100046637740771 776147607763771 300076777763771 756177607760771 557347607740661

Jo u

Shared International Type (SIT) 1 60 1329 2286 33 53 34 42 1750 1247 92 70 91 2271 482 479 181 647 592 ORPHAN ORPHAN Not in SITVIT Not in SITVIT Not in SITVIT Not in SITVIT Not in SITVIT Not in SITVIT Not in SITVIT Not in SITVIT Total

Prevalence no (%) 124 (67.4) 11 (6) 1 (0.5) 6 (3.3) 5 (2.7) 1 (0.5) 1 (0.5) 1 (0.5) 1 (0.5) 2 (1.1) 2 (1.1) 1 (0.5) 1 (0.5) 1 (0.5) 4 (2.2) 1 (0.5) 4 (2.2) 2 (1.1) 5 (2.7) 1 (0.5) 1 (0.5) 1 (0.5) 1 (0.5) 1 (0.5) 1 (0.5) 1 (0.5) 1 (0.5) 1 (0.5) 1 (0.5) 184 (100)

25

Journal Pre-proof Credit author statement

Nolwazi Londiwe Bhembe: Data curation, Investigation, Formal analysis, Validation, Visualization, Writing - Original draft preparation.

Ezekiel Green: Conceptualization,

Jo u

rn

al

Pr

e-

pr

oo

f

Methodology, Software, Resources, Supervision.

26

Journal Pre-proof Highlights 

This study gives a general overview of multidrug-resistant strains and circulating strains present in the Eastern Cape, South Africa.



Diverse spoligotype patterns were detected in this study, which delineates a high diversity of Mycobacterium tuberculosis strains disseminating in the Eastern Cape Province.

f

The Beijing strain was the predominant strain among the multidrug-resistant

rn

al

Pr

e-

pr

oo

tuberculosis isolates.

Jo u



27