Whole genome sequencing analysis to evaluate the influence of T2DM on polymorphisms associated with drug resistance in M. tuberculosis

doi:10.21203/rs.3.rs-1451730/v1

Background

Type 2 diabetes mellitus (T2DM) has been associated with treatment failure, and the development of drug resistance in tuberculosis (TB). Also, whole-genome sequencing has provided a better understanding and allowed the growth of knowledge about polymorphisms in genes associated with drug resistance. Considering the above, this study analyzes genome sequences to evaluate the influence of type 2 diabetes mellitus in the development of mutations related to tuberculosis drug resistance. M. tuberculosis isolates from individuals with (n=74), and without (n=74) type 2 diabetes mellitus was recovered from online repositories, and further analyzed.

Results

The results showed the presence of 591 SNPs with similar proportions between diabetics, and non-diabetics individuals (48.2% vs. 51.7%), but with no significant relationship. A greater number of mutations associated with rifampicin resistance was observed in the T2DM-TB individuals (19.3% vs. 12.5%), and the exclusive presence of rpoBQ432L, rpoBQ432P, rpoBS441L, and rpoBH445L variants. While these variants are not private to T2DM-TB cases they are globally rare highlighting a potential role of T2DM. The phylogenetic analysis showed 12 sublineages, being 4.1.1.3, and 4.1.2.1 the most prevalent in T2DM-TB individuals but not differing from those most prevalent in their geographic location. Four clonal complexes were found, however, no significant relationship with T2DM was observed. Samples size and potential sampling biases prevented us to look for significant associations.

Conclusions

The occurrence of globally rare rifampicin variants identified only in isolates from individuals with T2DM could be due to the hyperglycemic environment within the host. Therefore, further studies about the dynamics of SNPs’ generation associated with antibiotic resistance in patients with diabetes mellitus are necessary.

diabetes

tuberculosis

drug resistance

polymorphisms

Tuberculosis (TB) maintains a severe impact on global public health; according to the World Health Organization (WHO), in 2020 it caused 1.5 million deaths, and 10 million people sickened. One of the factors associated with the disease´s high mortality is the development of drug resistance. In 2018, approximately 500,000 individuals developed resistance to rifampicin, and from this patient’s group, 82% developed multi-drug resistance[1].

Studies have observed an increased risk of treatment failure, and the development of drug resistance in people with type 2 diabetes mellitus (T2DM). However, the impact of this condition on the development of drug resistance remains unclear[2–4]. A possible explanation may lie in the altered pharmacokinetics of T2DM-TB cases leading to heterogenous drug exposure during treatment[3].

Whole-genome sequencing (WGS) has opened the possibility to increase our understanding of drug resistance mechanisms. Particularly, the identification of the presence of polymorphisms in target genes associated with the performance of anti-tuberculosis drugs will aid in the development of the predictive diagnosis of drug resistance in TB[5].

Considering all the above, the present study aims to characterize the polymorphisms associated with drug resistance, and to determine the influence of diabetes in the occurrence of these mutations by analyzing whole genome sequences of M. tuberculosis (Mtb) from individuals with TB, and T2DM-TB binomial.

2.1 Patient characteristics

From nearly 1 million Mtb WGS currently available in the repositories analyzed, only 827 genomes had clinical information related to the resistance profile and confirmation of the patient with or without T2DM, and with L4 in the metadata. From this group, only 74 genomes from individuals with T2DM had the required sequencing depth, and coverage to be included in the study. An additional set of 74 WGS from individuals with drug resistance and without T2DM were retrieved and randomly included to match both study groups. Individuals carrying the 148 Mtb strains were mostly male 84 (57.4%), with a mean age of 47 ± 13 years at the time of sputum sample collection.

Patients came from nine countries, mainly Georgia with 36 participants (24%), Mexico, 34 (22.9%), Moldova, 20 (13.5%), and Belarus with 15 (10.1%). According to the type of resistance, four groups of isolates were identified: 40 were monoresistant (MR) (27%), 18 polyresistant (POL) (12.1%), 62 multidrug-resistant (MDR) (41.8%), and 28 extensively drug resistant (XDR) (18.9%). (Table 1)

Table 1

Epidemiological characteristics of the sample
	T2DM isolates^a	TB isolates^b	Total
	(n = 74)	(n = 74)	(n = 148)
Sex
Male	44(59%)	40 (54%)	84(57.4%)
Female	30 (41%)	34 (45%)	64(42.6%)
Age
Mean	50	42	47
(± SD)	(± 10)	(± 15)	(± 13)
Country
Georgia	8 (10.8%)	28 (37%)	36 (24.3%)
Mexico	23 (31%)	11 (14.8%)	34 (23%)
Moldova	6 (8.1%)	14 (18.9%)	20 (13.5%)
Belarus	4 (5.4%)	11 (14.8%)	15 (10.1%
Romania	2 (2.7%)	10 (13.5%)	12 (8.1%)
Peru	11 (14.8%)	0 (0.0%)	11 (7.4%)
Indonesia	9 (12.1%)	0 (0.0%)	9 (6.1%)
Spain	8 (10.8%)	0 (0.0%)	8 (5.4%)
Azerbaijan	3 (4%)	0 (0.0%)	3 (2%)
Resistance type^c
MR	20 (27%)	20 (27%)	40 (27%)
POL	9 (12.1%)	9 (12.1%)	18 (12.1%)
MDR	31 (48.8%)	31 (48.8%)	62 (41.8%)
XDR	14 (18.9%)	14 (18.9%)	28 (18.9%)
Drug-resistance
Isoniazid	58 (67.5%)	55 (74.3%)	113 (76.3%)
Rifampicin	50 (67.5%)	53 (71.6%)	103 (69.5%)
Streptomycin	21 (28.3%)	42 (56.7%)	63 (42.5%)
Ethambutol	27 (36.4%)	37 (50%)	64 (43.2%)
Pyrazinamide	27 (36.4%)	24 (32.4%)	51 (34.4%)
Amikacin	19 (25.6%)	30 (40.5%)	49 (33.1%)
Fluoroquinolones	16 (21.6%)	17 (22.9%)	33 (22.3%)
^aT2DM: Host with diabetes mellitus type 2 and tuberculosis. ^bTB Isolates: Host without diabetes mellitus type 2 and tuberculosis. ^cMR: mono-resistant tuberculosis, POL: poly-resistant tuberculosis, MDR: multidrug-resistant tuberculosis, XDR: extensively drug resistant tuberculosis

2.2 Characterization of variants associated with resistance

Genotypic resistance in the isolates was observed for several drugs, with the highest proportion of resistance to isoniazid (INH) with 113 isolates (76.3%), and rifampicin (RIF) including 103 strains (69.5%), followed by ethambutol (EMB) 64 (43.2%), streptomycin (STR) 63 (42.5%), and pyrazinamide (PZA) in 51 strains (34.4%). Genotypic resistance to second-line drugs had a lower representation with only 49 sequences resistant to amikacin (AMK) (33.1%), and 33 (22.3%) to fluoroquinolones (FQ). Among these drugs, only STR showed a significant difference in the number of WGS with this type of resistance with 21 (28.3%) isolates in the T2DM group vs 42 (56.7%) in the TB group (p = .00048).

During this study were identified 431 resistance-related SNPs in the diabetic, and non-diabetic groups; 207 (48%) vs 224 (52%) respectively. Only 23 (12 [2.7%] vs 11 [2.5%]) were classified as non-fixed SNPs. The SNPs mean was 2.79 vs 3.02 with a 2.7 vs 2.8 standard deviation, with no significant association in the distribution between groups (p = 0.4413). Regarding the distribution of polymorphisms by type of resistance, the sequences from individuals with diabetes presented a higher frequency in isolates with monoresistance, 20 (7%) vs 21 (6.9%), and polyresistant, 25 (8.8%) vs 14 (4.6%). In contrast, lower frequencies of polymorphisms were identified in isolates classified as MDR, 136 (47.7%) vs 159 (52.0%), and XDR, 104 (36.5%) vs 112 (36.6%); however, non-significant differences were identified.

Fifty-six high confidence resistance-related variants were identified, from which, 16 had a proportion greater than 1% in the dataset. Genotypic resistance was given by genomic variants associated with anti-tuberculosis drugs, and mainly concentrated in genes associated with this property. Resistance to isoniazid was observed in 76 isolates (26.6%) from the T2DM group, vs 84 on the isolates (27.4%) with only TB, rifampicin 67 (23.51%) vs 66 (21.5%), ethambutol 33 (11.58%) vs 45 (14.71%), and streptomycin 26 (9.12%) vs 43 (14.0%).

Three SNPs were observed with a frequency greater than 15% in both groups; katG S315T was found in 91 sequences (61.4%), and fabG1-inhA in 62 (41.8%); these were the only genes associated with resistance to isoniazid observed in the dataset. These mutations were followed by rpoB S450L in 58 isolates (39.1%). (Table 2) Nevertheless, no differences on these mutations were observed in the groups.

Table 2

Allele frequency of SNPs associated with drug resistance^a
Antibiotic	Gene	Aminoacid	Frequency	Frequency	Frequency
resistance		change	T2DM^b	TB^c	Total
			(n = 74)	(n = 74)	(n = 148)
			n(%)	n(%)	n(%)
Isoniazid
	katG	S315T	44 (59.4%)	47 (63.5%)	91 (61.4%)
	fabG1-inhA	c15t	28 (37.8%)	34 (45.9%)	62 (41.8%)
Rifampicin
	rpoB	S450L	28 (37.8%)	30 (40.5%)	58 (39.1%)
	rpoB	H445D	5 (6.7%)	11 (14.8%)	16 (10.8%)
	rpoB	D435Y	2 (2.7%)	6 (8.1%)	8 (5.4%)
	rpoC	V483G	8 (10.8%)	10 (13.5%)	18 (12.1%)
	rpoC	V483A	3 (4.0%)	2 (2.7%)	5 (3.38%)
Streptomycin
	rpsL	K88R	7 (9.4%)	15 (20.2%)	22 (14.8%)
	rpsL	K43R	4 (5.4%)	14 (18.9%)	18 (12.1%)
Ethambutol
	embB	M306I	12 (16.2%)	9 (12.1%)	21 (14.1%)
	embB	Q497R	5 (6.7%)	12 (16.2%)	17 (11.4%)
	embB	S297A	4 (5.4%)	9 (12.1%)	13 (8.7%)
Pyrazinamide
	pncA	D49G	5 (6.7%)	6 (8.1%)	11 (7.4%)
	pncA	L120P	7 (9.4%)	2 (2.7%)	9 (6.8%)
Fluoroquinolone
	gyrA	A90V	4 (5.4%)	6 (8.1%)	10 (6.7%)
	gyrA	D94G	4 (5.4%)	3 (4.0%)	8 (5.4%)
^aSNPs with a proportion greater than 3% in the sample (n = 148), ^bT2DM: Host with diabetes mellitus type 2 and tuberculosis
^cTB: Host without diabetes mellitus type 2 and tuberculosis

In only two SNPs was detected a difference in a greater proportion than 10% between diabetic, and non-diabetic groups; these were, rpsL K88R (9.4% vs 20.2%), and rpsL K43R (5.41% vs 18.9%), both related with resistance to STR. Surprisingly, we found a higher proportion of specific variants related to rifampicin resistance in people with T2DM, 17 (19.3%) vs 11 (12.5%), respectively; four mutations were identified exclusively in isolates from T2DM group; rpoB H445L (observed in four isolates) and rpoB Q432L, rpoB Q432P and rpoB S441L (detected in only one isolate each). Likewise, it was found that in 20 isolates carrying the rpoB S450L variant also had the compensatory mutation rpoC V483A/V483G. Among these, eight isolates (10.8%) from patients with T2DM carried two mutations, two rpoC V483A, and six rpoC V483G, meanwhile, in 12 strains (16.2%) from the non-diabetic group these mutations were also observed, 2 rpoC V483A, and 10 rpoC V483G.

2.3 Phylogenetic characterization

Phylogenetic analysis of the isolates showed that 148 isolates were included in 12 sublineages of L4, being the most frequent 4.3.3(LAM9) with 31 sequences (21%), and 4.10 (PGG3) with 29 (20%), followed by 4.1.2.1 (Haarlem) with 25 (16%), and 4.2.1 (Ural) with 24 (16%). The LAM9 (22%), Ural (20.9%), and Haarlem (19%) sublineages showed the highest frequency of MDR isolates, on the other hand, LAM9 concentrated the highest proportion of XDR tuberculosis (43%).

Regarding the phylogenetic distribution based on diabetes comorbidity, differences in the prevalence of some sublineages were observed. Sublineages 4.1.1.3 (X3/X1), 13% vs 3%, and 4.1.2.1 (Haarlem), 24% vs 9%, were predominantly found in isolates from patients with T2DM, respectively. By comparison, the sublineages 4.3.3 (LAM9), 28% vs. 14%, 4.10 (PGG3), 26% vs. 14%, and 4.2.1 (Ural), 23% vs. 9%, were found in a higher proportion within non-diabetic individuals.

According to the drug resistance profile, 40% of the MR isolates belonged to the 4.10 (PGG3) sublineage, and 33.3% of the POL isolates were found in the 4.1.2.1 (Haarlem) sublineage. The highest proportion of MDR sequences was classified into three sublineages: 19.3% in 4.1.2.1 (Haarlem), 20.9% in 4.2.1 (Ural), and 22.5% in 4.3.3 (LAM9). XDR isolates were mainly classified in the 4.2.1 (Ural) sublineage 21.4%, and 4.3.3 (LAM9) 42.8%. (Table 3). Non-significant differences were observed with drug resistance pattern, and sublineage association in the groups.

Phylogenetic analysis in the 148 isolates, with a 12 SNPs cut-off point, identified 36 isolates (24%) grouped in four clonal complexes (C): C1, with sublineage 4.3.3 (LAM9), composed of five isolates, all of them from the TB group, and from Belarus; C2, with lineage 4.2.1 (Ural), which includes 17 sequences, six T2DM vs 11 TB, 88% originating from Moldova; C3, with sublineage 4.1.1.3 (X3/X1), composed of nine isolates, six T2DM vs three TB) all from Mexico and C4 with sublineage 4.10 (PGG3), which includes five sequences from Georgia.

Table 3. Mtb L4 sublineages identified by WGS

		T2DM isolates^a	TB isolates^b	Total
Sublineage	Code	(n=74)	(n=74)	(n=148)

PGG3	4.10	10 (13.5%)	19 (25.6%)	29 (19.5%)
X	4.1.1	3 (4.0%)	2 (2.7%)	5 (3.3%)
X2	4.1.1.1	0 (0.0%)	1 (1.3%)	1(.06)
X3/X1	4.1.1.3	10 (13.5%)	2 (2.7%)	12 (8.1%)
T1;H1	4.1.2	1 (1.3%)	2 (2.7%)	3 (2%)
Haarlem	4.1.2.1	18 (24.3%)	7 (9.4%)	25 (16.8%)
Ural	4.2.1	7 (9.4%)	17 (22.9%)	24 (16.2%)
LAM3	4.3.2	5 (6.7%)	0 (0.0%)	5 (3.3%)
LAM9	4.3.3	10 (13.5%)	21 (28.3%)	31 (20.9%)
LAM1	4.3.4.1	3 (4.0%)	1 (1.3%)	4 (2.7%)
LAM11	4.3.4.2	1 (1.3%)	4 (5.3%)	5 (3.3%)
S	4.4.1	1 (1.3%)	0 (0.0%)	1(.06)

^aT2DM: Host with diabetes mellitus type 2 and tuberculosis.

^bTB Isolates: Host without diabetes mellitus type 2 and tuberculosis

A pattern of six mutations was observed in all C1 isolates: katG S315T, rss (position 1472359), embA (position 4243221), embB Q497R, pncA D49G, and rpoB H445D. Similarly, C2 included 14 sequences (82.3%), and shared the katG S315T, fabG1-inhA, rpoB S450L, and rpsL K88R mutation pattern, all these isolates came from Moldova. However, no relationship between diabetes comorbidity, and the presence of these variant patterns were identified. (Fig. 1)

The occurrence of rifampicin resistance-associated variants identified only in patients with T2DM had a heterogeneous phylogenetic distribution suggesting that they are not driven by geographically prevalent strains. The variants observed in a single isolate were rpoB S441L, and rpoB Q432P, belonging to sublineage 4.10 (PGG3), whereas rpoB Q432L, was observed in another isolate, and was found in sublineage 4.3.4.2 (LAM11). On the other hand, the rpoB H445L variant was observed in four isolates from three different sublineages 4.1.2 (T1;H1), 4.3.3 (LAM9), and 4.3.4.2 (LAM11).

Drug resistance generation in Mtb is caused by chromosomal mutations resulting in clones that proliferate, and subsequently fixate on the host[6]. The mutation rate of Mtb is low with an estimated 0.3–0.5 SNPS changes per year[7, 8], but this condition seems to be affected by the internal environment of the host, and by the pharmacokinetic variability of the plasma concentration of anti-tuberculosis drugs[9]. In this study, genome sequences of Mtb from individuals with drug resistant tuberculosis, and with or without T2DM, were analyzed to determine the impact of T2DM on the generation of drug resistance, given that diabetes increases the risk of developing drug resistance, and multidrug resistance[2–4, 10], also modifies the host’s internal environment by promoting the generation of reactive oxygen species in mycobacteria[11], which could lead to changes in the mutation rate in these subjects.

The polymorphism analysis showed that the most frequent variants were associated with the first-line drugs, isoniazid, and rifampicin. Mutations in the katG S315T and fabG1-inhA genes showed the highest proportions in both groups of individuals, this is consistent with previous studies, [12, 13] and could be due to the fact that these variants are the first to appear in the genome of resistant strains from different parts of the world[14].

In the T2DM group four specific rpoB gene variants were found: H445L, S441L, Q432P, and Q432L. These variants have been reported in studies where comorbidities are not specified[15–18]. According to the recent Catalogue of mutations related with drug resistance in Mycobacterium tuberculosis complex of WHO, these mutations have a low prevalence (1.17%, 0.26%, 0.29% and 0.21% respectively), and have been classified as SNPs associated with phenotypic resistance against rifampicin[19]. These variants were found in four different sublineages: 4.10 (PGG3), 4.3.4.2 (LAM11), 4.1.2 (T1; H1) and 4.3.3 (LAM9). Only two isolates with the H445L variant were found in a clonal complex (C2), which could indicate that remaining mutations were not acquired through clonal transmission but were the product of processes specific of the host, possibly by consequence to the complications in the drug treatment in T2DM individuals, and partially explaining the increased risk for the development of drug resistance [2–4]. Also with relation to RIF resistance, the rpoB Q432L, and rpoB Q432P mutations have been identified in isolates with a high minimum inhibitory concentration (MIC), while rpoB S441L, and rpoB H445L showed medium, and low MIC respectively.[17]. Likewise, pharmacokinetics studies have found that individuals with T2DM had 25% lower average plasma concentration, and an increase in time to reach the maximum concentration of rifampicin[3]. The rpoB Q432L and rpoB S441L variants, have been associated with a moderate fitness cost in in vitro studies[20], the previously mentioned changes in drug distribution could favor the proliferation of bacterial population carrying these mutations.

As explained above, the specific occurrence of these mutations in the T2DM group could be due to poor glycemic control that alters tissue perfusion, and rifampicin pharmacokinetics; however, the sample size in which we found these variants prevented us from looking for associations between the groups. Further studies with a larger number of individuals are needed to determine whether certain variants in the rifampicin resistance determining region in subjects with diabetes could be due to internal environmental conditions resulting from a hyperglycemic condition.

The phylogenetic analysis showed heterogeneity in the 12 sublineages found in both groups of individuals. The 4.1.1.3 (X3/X1), and 4.1.2.1 (Haarlem) sublineages included the higher number of isolates coming from individuals with T2DM. In particular, the largest difference in the isolates’ proportion between the groups was observed in the X3/X1 lineage, this could be linked to the country of origin of these isolates, since 83% came from Mexico. A similar situation happened with Haarlem sublineage, where 66% of the isolates from this sublineage, came from Mexico, and Indonesia. A total of 33 isolates were grouped into four clonal complexes, of which only C3, was composed of sublineage 4.1.1.3 (X3/X1) and showing a higher proportion of isolates from individuals with T2DM (six TB-T2DM vs. three TB). The high proportion, and differences observed in the distribution of isolates between diabetic, and non-diabetic individuals in the lineages from these countries could be explained with the fact that Mexico, and Indonesia are among the top 3 places in the worldwide prevalence of T2DM[21].

It has been described that the phylogenetic sublineages of M. tuberculosis complex are associated with specific geographic regions[22, 23]; therefore, the fact that diabetes was found with higher proportions in certain lineages could be a consequence of the disease burden in the specific countries; moreover, although the relationship between lineage, and glycemic level has not been characterized, a possible association between dysglycemia, and the distribution of Haarlem strains has been recently described[24]. To establish a clearer relationship, are required studies incorporating a larger number of WSGs from individuals with T2DM from different lineages.

The strength of this study lies in the fact that, to the best of our knowledge, there are no studies using WSG to analyze the impact of T2DM on the emergence of SNPs associated with drug resistance in L4 isolates the largest Mtb lineage distributed globally. Furthermore, the main limitation was the lack of information about comorbidities in genomic sequences of M. tuberculosis isolates available in repositories. From the 827 drug resistant isolates initially identified, only 74 fulfilled the criteria to be included in the TB-T2DM group, this stands out the need to include information related to comorbidity in the metadata available in the genome repositories. This would be of great help to improve further analysis of genomes, also to evaluate the specific effect of the host in the development of drug resistant in TB.

Finally, the lack of differences in the polymorphisms between the T2DM-TB, and TB groups could be due to the fact that these sequences come from isolates that had already been diagnosed with a resistance profile, as evidenced by the high proportion of fixed SNPs identified in the sample[25], which hinders the identification of the time of appearance of the resistance mutations, and the generation of these changes as consequence of the internal environment caused by type 2 diabetes mellitus in the host.

With this study it was possible to provide preliminary information about the impact of T2DM on drug resistance in whole genome sequences of M. tuberculosis, the results showed no significant association among the number of polymorphisms associated with resistance, sublineages, and the presence of T2DM. However, four variants in rpoB were observed only in individuals with T2DM. Undoubtedly, further studies about the dynamics of SNPs’ generation associated with antibiotic resistance in patients with diabetes mellitus are necessary, the results will have an important effect on the establishment of new therapeutic, attention, and control models for tuberculosis in the context of T2DM.

5.1 Sample selection

A search for whole genome sequences (WGSs) of Mtb from patients with drug resistance, and with or without T2DM was carried out in public databases; GenBank of the National Center for Biotechnology Information, the European Nucleotide Archive, and TB Portals of the National Institute of Health. In addition, 8 Mtb sequences were kindly provided by Dr. Iñaki Comas Espadas from Centro de Investigaciones Biomédicas de Valencia, Spain (manuscript in preparation). Only the sequences that incorporate the following information were included in the study: presence/absence of a diagnosis of diabetes in the patient, phenotypic and genotypic profile of drug resistance, coverage greater than 99%, depth greater than 100X, and belonging to lineage 4 (L4). Furthermore, those sequences were grouped according to host characteristics into two groups: 1) Host with TB drug resistant, and diabetes mellitus type 2 (T2DM), and 2) Host with TB drug resistant, and without diabetes mellitus type 2.

Accession numbers of genomes selected are available in the Additional file 1.

5.2 Bioinformatic analysis

The bioinformatic analysis was carried out using a previously validated pipeline (http://tgu.ibv.csic.es/?page_id=1794)[26, 27]. The workflow was composed by four sections: 1) Read cleaning and contamination removal, done by means of Trimmomatic v0.36, and Kraken v0.10.5; 2) Mapping, and coverage calculation, with BWA v0.7.12, and Bedtools v2.26; 3) Variant (SNPs and indels) calling was carried out by using Samtools v1.3.1, VarScan v2.3.7, and GATK v 4.0.2.1 and, 4) Variant filtering performed by using a customed Python code.

Genomic variants that were present in at least 20 reads, and at ≥ 90% of frequency within each isolate were used to detect phylogenetic mutations and confirm the belonging to L4 and sublineage inclusion. Variants with at least 10 reads at ≥ 10% to ≤ 90% frequency were called no fixed-SNP, and used to detect antibiotic resistance, and to confirm or predict the drug resistant profile.

5.3 Phylogenetic analysis

The phylogenetic tree was inferred from the resulting maximum likelihood SNP alignment in RAxML v8.2.4, by applying the General Time Reversible model of nucleotide substitution with the gamma distribution (GTRGAMMA) as used in other studies[26]. A 12 SNPs cut-off point was used to determine clonal complexes as recommended in previous studies[28, 29]. Phylogenetic tree visualization was done with iTOL v5 (https://itol. embl.de)[30].

5.4 Statistical analysis

The data collected were summarized using descriptive statistics, to identify the differences in polymorphisms between the groups, the chi-square test, and spearman correlation were performed using SPSS v.25 software.

T2DM: Type 2 diabetes mellitus; TB: Mycobacterium tuberculosis; WGS: Whole-genome sequencing; SNP: Single nucleotide polymorphism; MR: Monoresistant; XDR: Extensively drug resistant; POL: Polyresistant; MDR: multidrug-resistant; INH: Isoniazid; RIF: Rifampicin; EMB: Ethambutaol; STR: Streptomycin; PZA: Pyrazinamide; AMK: Amikacin; FQ: Fluoroquinolones; L4: Lineage 4; MIC: minimum inhibitory concentration.

Supplementary Information

The online version contains supplementary material available at: Additional file 1: Accession numbers of genomes and databases.

Acknowledgements

Computational infrastructure and support were provided by Biomedical Institute of Valencia IBV-CSIC, Valencia. We are grateful to Secretaría de Salud de Veracruz for providing some of the samples used in this study.

Authors’ contributions

RZ and GA-BH designed and coordinated the study, I-CM and CF-MM processed the genome dataset, GA-BH and DE-PM performed data analyses, IC contributed to data interpretation. RZ, GA-BH and DE-PM wrote the manuscript. All authors read and approved final manuscript.

Funding

Gustavo Bermúdez was funded by CONACyT, fellow number 710466. Roberto Zenteno was funded by CONACYT: "Fondo sectorial de investigación para la educación” A1-S-22956.

Availability of data and materials

The datasets supporting the conclusions of this article are available in the BioSample and ENA repositories. The accession numbers can be found in the Additional file 1.

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Author details

1. Health Sciences Doctoral Program. Health Sciences Institute. University of Veracruz, Mexico. 2. Biomedical Institute of Valencia IBV-CSIC, Valencia, Spain. 3. CIBER of Epidemiology and Public Health, Madrid, Spain. 4. Public Health Institute, University of Veracruz, Mexico. 5. Multidisciplinary network of tuberculosis research, Mexico. www.remitb.org6.

Global tuberculosis report 2021: World Health Organization; 2021.
Kornfeld H, Sahukar SB, Procter-Gray E, Kumar NP, West K, Kane K, Natarajan M, Li W, Babu S, Viswanathan V: Impact of diabetes and low body mass index on tuberculosis treatment outcomes. Clinical Infectious Diseases 2020, 71(9):e392-e398.
Alfarisi O, Mave V, Gaikwad S, Sahasrabudhe T, Ramachandran G, Kumar H, Gupte N, Kulkarni V, Deshmukh S, Atre S: Effect of diabetes mellitus on the pharmacokinetics and pharmacodynamics of tuberculosis treatment. Antimicrobial agents and chemotherapy 2018, 62(11):e01383-01318.
Yoon YS, Jung J-W, Jeon EJ, Seo H, Ryu YJ, Yim J-J, Kim YH, Lee B-H, Park YB, Lee BJ: The effect of diabetes control status on treatment response in pulmonary tuberculosis: a prospective study. Thorax 2017, 72(3):263–270.
Schleusener V, Köser CU, Beckert P, Niemann S, Feuerriegel S: Mycobacterium tuberculosis resistance prediction and lineage classification from genome sequencing: comparison of automated analysis tools. Scientific reports 2017, 7(1):1–9.
Al-Saeedi M, Al-Hajoj S: Diversity and evolution of drug resistance mechanisms in Mycobacterium tuberculosis. Infection and drug resistance 2017, 10:333.
Hakamata M, Takihara H, Iwamoto T, Tamaru A, Hashimoto A, Tanaka T, Kaboso SA, Gebretsadik G, Ilinov A, Yokoyama A: Higher genome mutation rates of Beijing lineage of Mycobacterium tuberculosis during human infection. Scientific reports 2020, 10(1):1–11.
Takiff HE, Feo O: Clinical value of whole-genome sequencing of Mycobacterium tuberculosis. The Lancet Infectious Diseases 2015, 15(9):1077–1090.
Allue Guardia A, Garcia JI, Torrelles JB: Evolution of drug-resistant Mycobacterium tuberculosis strains and their adaptation to the human lung environment. Frontiers in Microbiology 2021, 12:137.
Perez-Navarro LM, Restrepo BI, Fuentes-Dominguez FJ, Duggirala R, Morales-Romero J, López-Alvarenga JC, Comas I, Zenteno-Cuevas R: The effect size of type 2 diabetes mellitus on tuberculosis drug resistance and adverse treatment outcomes. Tuberculosis 2017, 103:83–91.
Ferlita S, Yegiazaryan A, Noori N, Lal G, Nguyen T, To K, Venketaraman V: Type 2 diabetes mellitus and altered immune system leading to susceptibility to pathogens, especially Mycobacterium tuberculosis. Journal of clinical medicine 2019, 8(12):2219.
Lempens P, Meehan CJ, Vandelannoote K, Fissette K, de Rijk P, Van Deun A, Rigouts L, de Jong BC: Isoniazid resistance levels of Mycobacterium tuberculosis can largely be predicted by high-confidence resistance-conferring mutations. Scientific reports 2018, 8(1):1–9.
Charoenpak R, Santimaleeworagun W, Suwanpimolkul G, Manosuthi W, Kongsanan P, Petsong S, Puttilerpong C: Association Between the Phenotype and Genotype of Isoniazid Resistance Among Mycobacterium tuberculosis Isolates in Thailand. Infection and drug resistance 2020, 13:627.
Manson AL, Cohen KA, Abeel T, Desjardins CA, Armstrong DT, Barry CE, Brand J, Chapman SB, Cho S-N, Gabrielian A: Genomic analysis of globally diverse Mycobacterium tuberculosis strains provides insights into the emergence and spread of multidrug resistance. Nature genetics 2017, 49(3):395–402.
Ng KC, Meehan CJ, Torrea G, Goeminne L, Diels M, Rigouts L, de Jong BC, André E: Potential application of digitally linked tuberculosis diagnostics for real-time surveillance of drug-resistant tuberculosis transmission: validation and analysis of test results. JMIR medical informatics 2018, 6(1):e12.
Bhembe NL, Nwodo UU, Govender S, Hayes C, Ndip RN, Okoh AI, Green E: Molecular detection and characterization of resistant genes in Mycobacterium tuberculosis complex from DNA isolated from tuberculosis patients in the Eastern Cape province South Africa. BMC infectious diseases 2014, 14(1):1–9.
Li M-c, Lu J, Lu Y, Xiao T-y, Liu H-c, Lin S-q, Xu D, Li G-l, Zhao X-q, Liu Z-g: rpoB Mutations and Effects on Rifampin Resistance in Mycobacterium tuberculosis. Infection and Drug Resistance 2021, 14:4119.
Miotto P, Cabibbe AM, Borroni E, Degano M, Cirillo DM: Role of disputed mutations in the rpoB gene in interpretation of automated liquid MGIT culture results for rifampin susceptibility testing of Mycobacterium tuberculosis. Journal of clinical microbiology 2018, 56(5):e01599-01517.
Catalogue of mutations in Mycobacterium tuberculosis complex and their association with drug resistance. 2021.
Gagneux S: Fitness cost of drug resistance in Mycobacterium tuberculosis. Clinical Microbiology and Infection 2009, 15:66–68.
Lin X, Xu Y, Pan X, Xu J, Ding Y, Sun X, Song X, Ren Y, Shan P-F: Global, regional, and national burden and trend of diabetes in 195 countries and territories: an analysis from 1990 to 2025. Scientific reports 2020, 10(1):1–11.
Chihota VN, Niehaus A, Streicher EM, Wang X, Sampson SL, Mason P, Källenius G, Mfinanga SG, Pillay M, Klopper M: Geospatial distribution of Mycobacterium tuberculosis genotypes in Africa. PLoS One 2018, 13(8):e0200632.
Brynildsrud OB, Pepperell CS, Suffys P, Grandjean L, Monteserin J, Debech N, Bohlin J, Alfsnes K, Pettersson JO-H, Kirkeleite I: Global expansion of Mycobacterium tuberculosis lineage 4 shaped by colonial migration and local adaptation. Science Advances 2018, 4(10):eaat5869.
Lopez K, Arriaga MB, Aliaga JG, Barreda NN, Sanabria OM, Huang C-C, Zhang Z, García-de-la-Guarda R, Lecca L, Calçada Carvalho AC: Dysglycemia is associated with Mycobacterium tuberculosis lineages in tuberculosis patients of North Lima—Peru. PloS one 2021, 16(1):e0243184.
Séraphin MN, Norman A, Rasmussen EM, Gerace AM, Chiribau CB, Rowlinson M-C, Lillebaek T, Lauzardo M: Direct transmission of within-host Mycobacterium tuberculosis diversity to secondary cases can lead to variable between-host heterogeneity without de novo mutation: a genomic investigation. EBioMedicine 2019, 47:293–300.
Xu Y, Cancino-Muñoz I, Torres-Puente M, Villamayor LM, Borrás R, Borrás-Máñez M, Bosque M, Camarena JJ, Colomer-Roig E, Colomina J: High-resolution mapping of tuberculosis transmission: Whole genome sequencing and phylogenetic modelling of a cohort from Valencia Region, Spain. PLoS medicine 2019, 16(10):e1002961.
Comas I, Coscolla M, Luo T, Borrell S, Holt KE, Kato-Maeda M, Parkhill J, Malla B, Berg S, Thwaites G: Out-of-Africa migration and Neolithic coexpansion of Mycobacterium tuberculosis with modern humans. Nature genetics 2013, 45(10):1176–1182.
Malm S, Linguissi LSG, Tekwu EM, Vouvoungui JC, Kohl TA, Beckert P, Sidibe A, Rüsch-Gerdes S, Madzou-Laboum IK, Kwedi S: New Mycobacterium tuberculosis complex sublineage, Brazzaville, Congo. Emerging infectious diseases 2017, 23(3):423.
Jajou R, Kohl TA, Walker T, Norman A, Cirillo DM, Tagliani E, Niemann S, de Neeling A, Lillebaek T, Anthony RM: Towards standardisation: comparison of five whole genome sequencing (WGS) analysis pipelines for detection of epidemiologically linked tuberculosis cases. Eurosurveillance 2019, 24(50):1900130.
Letunic I, Bork P: Interactive Tree Of Life (iTOL) v5: an online tool for phylogenetic tree display and annotation. Nucleic acids research 2021, 49(W1):W293-W296.

No competing interests reported.

Additionalfile11.xlsx

Whole genome sequencing analysis to evaluate the influence of T2DM on polymorphisms associated with drug resistance in M. tuberculosis

Status:

Version 1

Abstract

Figures

1. Introduction

2. Results

3. Discussion

4. Conclusions

5. Methods

Abbreviations

Declarations

References

Additional Declarations

Supplementary Files

Status:

Version 1