|Year : 2018 | Volume
| Issue : 3 | Page : 278-286
Genetic diversity of Mycobacterium tuberculosis in south coastal Karnataka, India, using spoligotyping
Kiran Chawla1, Ajay Kumar1, Vishnu Prasad Shenoy1, Devendra Singh Chauhan2, Pragya Sharma2
1 Department of Microbiology, Kasturba Medical College, MAHE, Manipal, India
2 Department of Microbiology & Molecular Biology, ICMR-National JALMA Institute for Leprosy & other Mycobacterial Diseases, Agra, India
|Date of Submission||23-Dec-2016|
|Date of Web Publication||18-Jun-2018|
Dr Kiran Chawla
Department of Microbiology, Kasturba Medical College, MAHE, Manipal 576 104, Karnataka
Source of Support: None, Conflict of Interest: None
| Abstract|| |
Background & objectives: Despite high occurrence of tuberculosis in India very little information is available about the genetic diversity of Mycobacterium tuberculosis isolates prevailing in coastal Karnataka, India. Thus, the present study was undertaken to explore the genetic biodiversity of M. tuberculosis isolates prevailing in south coastal region of Karnataka (Udupi District), India.
Methods: A total of 111 Mycobacterial isolates were cultured in Lowenstein Jensen (LJ) medium and after obtaining growth, DNA was extracted and spoligotyping was performed. SITVIT WEB database was used to locate families of spoligotypes.
Results: On analyzing the hybridization results of all 111 isolates on SITVIT WEB database 57 (51.35%) isolates were clustered into 11 Spoligotype International Types (SIT). The largest cluster of 14 (12.61%) isolates was SIT-48 (EAI1-SOM), followed by SIT-1942 (CAS1-Delhi) with 11 isolates (9.9%) and SIT-11 with seven (6.30%). Moreover, 23 isolates (20.72%) had unique spoligotypes and 31 (27.92%) were orphans. Spotclust analysis revealed that majority (67%) of orphan isolates were variants of CAS (37%) and EAI-5 (34%).
Interpretation & conclusions: The present study revealed high biodiversity among the circulating isolates of M. tuberculosis in this region with the presence of mixed genotypes earlier reported from north and south India along with certain new genotypes with unique SITs. The study highlights the need for further longitudinal studies to explore the genetic diversity and to understand the transmission dynamics of prevailing isolates.
Keywords: Genetic diversity - genotyping - Mycobacterium tuberculosis - spoligotyping
|How to cite this article:|
Chawla K, Kumar A, Shenoy VP, Chauhan DS, Sharma P. Genetic diversity of Mycobacterium tuberculosis in south coastal Karnataka, India, using spoligotyping. Indian J Med Res 2018;147:278-86
|How to cite this URL:|
Chawla K, Kumar A, Shenoy VP, Chauhan DS, Sharma P. Genetic diversity of Mycobacterium tuberculosis in south coastal Karnataka, India, using spoligotyping. Indian J Med Res [serial online] 2018 [cited 2019 Jul 18];147:278-86. Available from: http://www.ijmr.org.in/text.asp?2018/147/3/278/234608
In 2016, an estimated 10.4 million people developed tuberculosis and 1.67 million died from the disease globally, highlighting the devastation caused by Mycobacterium tuberculosis (MTB). India reported the largest number of cases (25%) of the global total . Comparing the share in global burden, the knowledge about the genetic diversity of strains prevailing in India is inadequate. Study on genetic diversity of MTB is important for effectively controlling the emerging drug-resistant strains. Further, it also enables us to understand the transmission dynamics and differentiate between reinfection and reactivation in patients with relapse. Application of molecular methods such as spoligotyping, restriction fragment length polymorphism (RFLP) and Mycobacterial interspersed repetitive units-Variable number of tandem repeats (MIRU-VNTRs) is useful for understanding the prevailing MTB genotypes and their transmission within the population. Over the years, spacer oligonucleotide genotyping (spoligotyping) has emerged as the most widely used method for MTB genotyping after IS6110-based fingerprinting, which requires large quantities of good quality DNA and has a limitation in isolates having low copies of IS6110 elements . Spoligotyping is based on DNA polymorphism in direct repeat (DR) locus of MTB. The latter contains multiple 36 base pair DRs which are well conserved and are interspersed with non-repetitive spacer sequences which are 34-41 base pairs long . On comparing DR regions of several isolates of MTB, it was observed that the deletions and insertions in DRs occur, but the order of the spacers remains same. Thus, variation in the sequence of spacers is used to determine genetic diversity among isolates by targeting the DR locus for in vitro DNA amplification. Spoligotyping is a simple, rapid and sensitive test, and can be performed directly from the clinical samples even if the bacteria are not viable or are in tissues in paraffin-embedded blocks .
The population of Udupi district in Karnataka State of India comprises a section of migratory people from other parts of Karnataka. It will be interesting to know about the circulating genotypes of MTB in this geographical area. Thus the present study was aimed to characterize the MTB isolates from pulmonary tuberculosis patients using spoligotyping and to identify the predominant MTB clades prevalent in Udupi district of south coastal Karnataka.
| Material & Methods|| |
A cross-sectional epidemiological study was carried out from January to August 2016 in the department of Microbiology, Kasturba Medical College (KMC), Manipal, Karnataka, in association with National JALMA Institute for Leprosy and other Mycobacterial Diseases (ICMR), Agra, after obtaining the approval of Institutional Research and Ethical Committee of KMC, Manipal.
Sample processing & DNA extraction:The stock isolates of MTB obtained from sputum of smear-positive pulmonary tuberculosis patients collected from designated microscopy centres of three taluks (Udupi, Kundapura and Karkala) in a stratified sampling technique from Udupi district from September 2011 to August 2014 in our previous study , were used. For spoligotyping, the stock isolates were subcultured on Lowenstein–Jensen (LJ) medium (HiMedia, Mumbai) and incubated for four to eight weeks at 37°C. After obtaining the growth on LJ medium, DNA extraction was performed using a commercially available kit from Qiagen (Hilden, Germany) as per manufacturer's instructions.
Spoligotyping: Spoligotyping was carried out by amplifying the whole DR region using the commercially available kit (Mapmygenome, Hyderabad) according to a standardized method using the designated primers' pairs of DRa and DRb . Master cycler gradient 5331 (Eppendorf, Germany) was used for DNA amplification. The amplified PCR products were hybridized with nitrocellulose membrane having covalently linked 43 spacer oligonucleotides following the manufacturer's instructions. MTB H37Rv and M. bovis were used as positive control and double-distilled water as a negative control. The hybridized fragments were identified using enhanced chemiluminescence system (GE Healthcare, UK). The spoligotypes were initially reported as 43 digits binary representation of 43 spacers, one (1) was scored for positive hybridization and zero (0) for negative hybridization. The binary codes were converted into octal codes. The octal codes were analyzed using the SITVIT WEB (http://www.pasteur-guadeloupe.fr:8081/SITVIT_ONLINE/) database to locate the families of spoligotypes . A cluster was defined as two or more isolates with identical spoligotype patterns. However, if the spoligotype pattern in the SITVIT WEB database corresponded to only one isolate, then the pattern was called 'unique'. The spoligotypes that did not match any pattern in the SITVIT WEB database were defined as 'orphan'. The orphans were further analyzed by 'Spotclust', which assigns families based on SpolDB3.0 (http://tbinsight.cs.rpi.edu/run_spotclust.html). The evolutionary relationship of the clinical isolates was analysed based on spoligotyping data by calculating NJ tree dendrogram [Figure 1] which was based on Jaccard Index and minimum spanning tree [Figure 2] using MIRU-VNTR plus (http://www.miru-vntrplus.org/MIRU/index.faces),,. Discriminatory power of spoligotyping was determined by calculating Hunter Gaston Discriminatory Index (HGDI). Association between demographic characteristics of the patients with the occurrence of clustered and non-clustered MTBisolates was examined using univariate odds ratio.
|Figure 1: A spoligotyping-based dendogram generated through N-J tree analysis using Mycobacterial Interspersed Repetitive Units-Variable Number of Tandem Repeats plus database showing genetic diversity of 111 clinical isolates from Udupi district, Karnataka, India.|
Click here to view
|Figure 2: A minimum spanning tree based on phylogenetic lineages. The structure of the tree is represented by branches (continuous, dashed and dotted lines) and circles representing individual pattern. The length of the branches represents the distance between patterns while the complexity of the lines (continuous, grey dashed and grey dotted) denotes the number of allele/spacer changes between two patterns: solid lines, 1 or 2 or more changes (thicker ones indicate a single change, while the thinner one indicates 2 changes); grey dashed lines represent 3 changes; and grey dotted lines represent 4 or more changes. The size of the circle is proportional to the total number of isolates in our study, illustrating unique isolates (smaller nodes) versus clustered isolates (bigger nodes). The colour of the circles indicates the lineage to which the specific pattern belongs.|
Click here to view
| Results|| |
A total of 111 clinical isolates of MTB obtained from sputum samples of pulmonary tuberculosis patients were spoligotyped and analyzed. On SITVIT WEB analysis, a total of 57 (51.35%) isolates were clustered into 11 spoligotype international types (SITs) with SIT-48 (EAI1-SOM) being the largest cluster of 14 (12.61%) isolates, followed by SIT-1942 (CAS1-Delhi) and SIT-11 (EAI3-IND) with 11 (9.9 %) and seven (6.3%) isolates, respectively. Further 23 isolates (20.72%) had unique spoligotypes and 31 (27.92%) were orphans [Table 1]. The most prevalent clades prevailing in this region were EAI1-SOM (15.31%), CAS1-Delhi (12%) and EAI-5 (10%). Beijing (4%) and Harlem (3%) families, which are emerging worldwide and of major concern due to their association with drug resistance, were found less prevalent in this region [Figure 3]. MANU 4.5 per cent (5/111), X1 1.8 per cent (2/111), T2 1.8 per cent (2/111), T1 5.4 per cent (6/111), T5 0.9 per cent (1/111), H1 0.9 per cent (1/111), H3 1.8 per cent (2/111), EAI2-Manilla 0.9 per cent (1/111) and LAM-6 0.9 per cent (1/111) were the other clades reported in this study [Table 1]. We further analyzed 31 orphans isolates by Spotclust (SPOLDB3 based) to identify their most probable families. Spotclust analyses revealed that majority (67%) of orphans in this study were variants of clades CAS (37%) and EAI-5 (34%) with 11 and 10 clinical isolates, respectively. One clinical isolate each was found to be the variant of EAI1-SOM, EAI3, T2 and Family 33. Further two clinical isolates were found to be variants of Family 35, whereas three clinical isolates were found to be variants of T1. The distribution of clades in three taluks of Udupi district is shown in [Figure 4]. Clinical and demographic characteristics of patients in relation to clustering of MTBisolates are shown in [Table 2]. Two out (1.80%) of 111 patients had multidrug-resistant tuberculosis (MDR-TB) which shared SIT-48 (EAI-SOM) and SIT-1(Beijing) spoligotypes. Further, four (3.60%) isolates were monoresistant to isoniazid, 27 (24.32%) were monoresistant to streptomycin, seven (6.30%) were polydrug-resistant other than MDR and 71 (63.96%) were sensitive to all first-line drugs. The detailed description of drug resistance in isolates with respect to lineage is given in [Table 3]. Thirteen (11.71%) patients were seropositive for HIV (mean age=41.54, SD=8.10, range: 30-57 yr), SIT-48 was shared by two isolates, one isolate each belonged to SIT-236, SIT-126, SIT-100, SIT-288, SIT-11, SIT-119, SIT-490, SIT-1268, SIT-52 and orphan. Discriminatory power of spoligotyping for 111 clinical isolates of Udupi district based on HGDI was found to be 0.93.
|Figure 3: A pie chart showing the distribution of MTB clades prevalent in Udupi district on the basis of SITVIT WEB analysis.|
Click here to view
|Figure 4: Distribution of clades in three taluks of Udupi district, Karnataka, India.|
Click here to view
|Table 1: Spoligotypes shared by clinical isolates of Mycobacterium tuberculosis evaluated by SITVIT WEB analysis (based on database search August 14, 2016)|
Click here to view
|Table 2: Clinical and demographic characteristics of patients in relation to clustering of Mycobacterium tuberculosis isolates|
Click here to view
|Table 3: Drug susceptibility profiles of clinical isolates with respect to lineage|
Click here to view
| Discussion|| |
The largest cluster in the present study belonged to EAI1-SOM (SIT-48) shared by 14 (12.61%) clinical isolates. EAI1-SOM has earlier been also reported in variable prevalence from south India ,,. East African Indian lineages originated in Guinea-Bissau (Africa) and believed to be transmitted to other continents through dispersion of modern humans from that region . EAI1-SOM, a sublineage of EAI originated in Somalia, has been reported from India, with high prevalence in southern regions, whereas prevalence in northern regions of India is reported variably low ,. SIT-1942 (CAS1-Delhi) was the second largest cluster in the present study. CAS1-Delhi lineage supposed to be originated from India was later disseminated to regions such as Saudi Arabia, Kenya, South Africa, Malaysia, Myanmar, Australia, the USA and parts of Europe through frequent migration ,. Beijing, the East Asian lineage originated in China, is frequently reported to be associated with drug resistance and ubiquitously prevalent worldwide . Low prevalence of Beijing has been reported (2.4-14.9%) in most studies in India except Mumbai (18.8%) and Assam (35.4%). It correlated with the presence of high drug resistance in these geographical regions. In the present study, it was found in about four per cent (4/111) of MTB isolates. One of these isolates, belonging to Beijing lineage, was detected as MDR. The low prevalence of Beijing lineage might be the possible reason for lesser number of multidrug resistance cases in Udupi district. MANU, an ancient lineage believed to be the probable ancestral of CAS and EAI with worldwide prevalence , was found to be less prevalent (4.5%) in the present study. It was represented by SIT-100, which was initially reported from Delhi . Some of the unique SITs found in this investigation, such as SIT-349, SIT-298, SIT-46, SIT-52, SIT-119, SIT-64 and SIT-1952, have been reported earlier from north India ,,,,,, and a few such as SIT-256, SIT-138 and SIT-1268 have been reported from south India ,,,. Spoligotypes of 31 clinical isolates from the present work did not match with the spoligo patterns available in SITVIT WEB database, suggesting the high evolutionary pressure on the isolates in this region which resulted in emergence of newer spoligotypes. Clustering of isolates was not found significantly associated with any of demographic variables and drug resistance patterns due to limited sample size. IS6110-RFLP analysis was considered as the gold standard for genotyping of MTB . However, MIRU-VNTR typing has now been considered as the new molecular reference tool for being simpler, faster, less DNA demanding and with a similar resolution than RFLP typing . We preferred spoligotyping as MIRU-VNTR typing was an expensive technique for limited resource laboratory like ours.
Although this was one of the first genotyping studies on clinical isolates of MTB in Udupi district of Karnataka, limited number of isolates and the inability to perform MIRU-VNTR typing were the major limitations.
In conclusion, our study showed the presence of highly diverse patterns of spoligotypes in Udupi district. High number of orphans (28%) in this study highlights the presence of a large number of genotypes which are yet to be identified. Longitudinal studies need to be done in this region to understand the transmission dynamics of MTB using different genotyping methods.
| Acknowledgment|| |
The authors are thankful to Shri Vijay Dua and Shrimati Rajani K, for help in experimental work.
Financial support & sponsorship: None.
Conflicts of Interest: None.
| References|| |
Desikan S, Narayanan S. Genetic markers, genotyping methods & next generation sequencing in Mycobacterium tuberculosis. Indian J Med Res
Jagielski T, Minias A, van Ingen J, Rastogi N, Brzostek A, Żaczek A, et al
. Methodological and clinical aspects of the molecular epidemiology of Mycobacterium tuberculosis
and other mycobacteria. Clin Microbiol Rev
Kremer K, van Soolingen D, Frothingham R, Haas WH, Hermans PW, Martín C, et al.
Comparison of methods based on different molecular epidemiological markers for typing of Mycobacterium tuberculosis
complex strains: Interlaboratory study of discriminatory power and reproducibility. J Clin Microbiol
Rao P, Chawla K, Shenoy VP, Mukhopadhyay C, Brahmavar V, Kamath A, et al.
Study of drug resistance in pulmonary tuberculosis cases in South Coastal Karnataka. J Epidemiol Glob Health
Kamerbeek J, Schouls L, Kolk A, van Agterveld M, van Soolingen D, Kuijper S, et al.
Simultaneous detection and strain differentiation of Mycobacterium tuberculosis
for diagnosis and epidemiology. J Clin Microbiol
Demay C, Liens B, Burguière T, Hill V, Couvin D, Millet J, et al.
SITVITWEB - A publicly available international multimarker database for studying Mycobacterium tuberculosis
genetic diversity and molecular epidemiology. Infect Genet Evol
Vitol I, Driscoll VJ, Kreiswirth B, Kurepina N, Bennett KP. Identifying Mycobacterium tuberculosis
complex strain families using spoligotypes. Infect Genet Evol
Jaccard P. New research on floral distribution. Bull Soc Vaudoise Sci Nat
Allix-Béguec C, Harmsen D, Weniger T, Supply P, Niemann S. Evaluation and strategy for use of MIRU-VNTR plus
, a multifunctional database for online analysis of genotyping data and phylogenetic identification of Mycobacterium tuberculosis
complex isolates. J Clin Microbiol
Weniger T, Krawczyk J, Supply P, Niemann S, Harmsen D. MIRU-VNTR plus
: A web tool for polyphasic genotyping of Mycobacterium tuberculosis
complex bacteria. Nucleic Acids Res
Hunter PR, Gaston MA. Numerical index of the discriminatory ability of typing systems: An application of Simpson's index of diversity. J Clin Microbiol
Singh J, Sankar MM, Kumar P, Couvin D, Rastogi N, Singh S, et al.
Genetic diversity and drug susceptibility profile of Mycobacterium tuberculosis
isolated from different regions of India. J Infect
Joseph BV, Soman S, Radhakrishnan I, Hill V, Dhanasooraj D, Kumar RA, et al.
Molecular epidemiology of Mycobacterium tuberculosis
isolates from Kerala, India using IS6110-RFLP, spoligotyping and MIRU-VNTRs. Infect Genet Evol
Kandhakumari G, Stephen S, Sivakumar S, Narayanan S. Spoligotype patterns of Mycobacterium tuberculosis
isolated from extra pulmonary tuberculosis patients in Puducherry, India. Indian J Med Microbiol
Liu H, Prugnolle F, Manica A, Balloux F. A geographically explicit genetic model of worldwide human-settlement history. Am J Hum Genet
Singh UB, Arora J, Suresh N, Pant H, Rana T, Sola C, et al.
Genetic biodiversity of Mycobacterium tuberculosis
isolates from patients with pulmonary tuberculosis in India. Infect Genet Evol
Singh UB, Suresh N, Bhanu NV, Arora J, Pant H, Sinha S, et al.
Predominant tuberculosis spoligotypes, Delhi, India. Emerg Infect Dis
Gupta A, Kulkarni S, Rastogi N, Anupurba S. A study of Mycobacterium tuberculosis
genotypic diversity & drug resistance mutations in Varanasi, North India. Indian J Med Res
Arora J, Bhalla M, Sidiq Z, Lal P, Behera D, Rastogi N, et al
. Predominance of Beijing genotype in extensively drug resistant Mycobacterium tuberculosis
isolates from a tertiary care hospital in New Delhi, India. Int J Mycobacteriol
Chatterjee A, D'Souza D, Vira T, Bamne A, Ambe GT, Nicol MP, et al.
Strains of Mycobacterium tuberculosis
from Western Maharashtra, India, exhibit a high degree of diversity and strain-specific associations with drug resistance, cavitary disease, and treatment failure. J Clin Microbiol
Devi KR, Bhutia R, Bhowmick S, Mukherjee K, Mahanta J, Narain K, et al.
Genetic diversity of Mycobacterium tuberculosis
isolates from Assam, India: Dominance of Beijing family and discovery of two new clades related to CAS1_Delhi and EAI family based on spoligotyping and MIRU-VNTR typing. PLoS One
Brudey K, Driscoll JR, Rigouts L, Prodinger WM, Gori A, Al-Hajoj SA, et al
. Mycobacterium tuberculosis
complex genetic diversity: Mining the fourth international spoligotyping database (SpolDB4) for classification, population genetics and epidemiology. BMC Microbiol
Sankar MM, Singh J, Diana SCA, Singh S. Molecular characterization of Mycobacterium tuberculosis
isolates from North Indian patients with extrapulmonary tuberculosis. Tuberculosis (Edinb)
Mathuria JP, Sharma P, Prakash P, Samaria JK, Katoch VM, Anupurba S, et al.
Role of spoligotyping and IS6110-RFLP in assessing genetic diversity of Mycobacterium tuberculosis
in India. Infect Genet Evol
Sharma P, Chauhan DS, Upadhyay P, Faujdar J, Lavania M, Sachan S, et al.
Molecular typing of Mycobacterium tuberculosis
isolates from a rural area of Kanpur by spoligotyping and mycobacterial interspersed repetitive units (MIRUs) typing. Infect Genet Evol
Shanmugam S, Selvakumar N, Narayanan S. Drug resistance among different genotypes of Mycobacterium tuberculosis
isolated from patients from Tiruvallur, South India. Infect Genet Evol
Thomas SK, Iravatham CC, Moni BH, Kumar A, Archana BV, Majid M, et al.
Modern and ancestral genotypes of Mycobacterium tuberculosis
from Andhra Pradesh, India. PLoS One
Braden CR, Crawford JT, Schable BA. Quality assessment of Mycobacterium tuberculosis
genotyping in a large laboratory network. Emerg Infect Dis
Supply P, Allix C, Lesjean S, Cardoso-Oelemann M, Rüsch-Gerdes S, Willery E, et al
. Proposal for standardization of optimized mycobacterial interspersed repetitive unit-variable-number tandem repeat typing of Mycobacterium tuberculosis
. J Clin Microbiol
[Figure 1], [Figure 2], [Figure 3], [Figure 4]
[Table 1], [Table 2], [Table 3]