Skip to main content
Download PDF

Association between fatty acid metabolism gene mutations and Mycobacterium tuberculosis transmission revealed by whole genome sequencing

BMC Microbiology volume 23, Article number: 379 (2023) Cite this article

A Correction to this article was published on 03 January 2024

This article has been updated

Abstract

Background

Fatty acid metabolism greatly promotes the virulence and pathogenicity of Mycobacterium tuberculosis (M.tb). However, the regulatory mechanism of fatty acid metabolism in M.tb remains to be elucidated, and limited evidence about the effects of gene mutations in fatty acid metabolism on the transmission of M.tb was reported.

Results

Overall, a total of 3193 M.tb isolates were included in the study, of which 1596 (50%) were genomic clustered isolates. Most of the tuberculosis isolates belonged to lineage2(n = 2744,85.93%), followed by lineage4(n = 439,13.75%) and lineage3(n = 10,0.31%).Regression results showed that the mutations of gca (136,605, 317G > C, Arg106Pro; OR, 22.144; 95% CI, 2.591-189.272), ogt(1,477,346, 286G > C ,Gly96Arg; OR, 3.893; 95%CI, 1.432–10.583), and rpsA (1,834,776, 1235 C > T, Ala412Val; OR, 3.674; 95% CI, 1.217–11.091) were significantly associated with clustering; mutations in gca and rpsA were also significantly associated with clustering of lineage2. Mutation in arsA(3,001,498, 885 C > G, Thr295Thr; OR, 6.278; 95% CI, 2.508–15.711) was significantly associated with cross-regional clusters. We also found that 20 mutation sites were positively correlated with cluster size, while 11 fatty acid mutation sites were negatively correlated with cluster size.

Conclusion

Our research results suggested that mutations in genes related to fatty acid metabolism were related to the transmission of M.tb. This research could help in the future control of the transmission of M.tb.

Peer Review reports

Background

Tuberculosis (TB) is a highly contagious infectious disease caused by Mycobacterium tuberculosis(M.tb) that primarily affects the lungs and spreads through the respiratory tract [1, 2]. It is estimated that 25% of the global population is infected with M.tb. Despite the large global efforts at curbing the spread of M.tb complex strains, 10.6 million new patients develop TB every year [3,4,5]. China still has the second-highest number of TB infections globally [6]. Therefore, having a thorough and comprehensive understanding of the transmission mechanisms of M.tb is of great significance in the prevention and treatment of TB.

The metabolism of fatty acids is critical to the survival of M.tb within the host. M.tb utilizes diverse lipids as major carbon and energy source during infection. Fatty acids are degraded via beta-oxidation to generate reduced power and energy [7,8,9]. At the same time, fatty acids play a crucial role in the composition of the cell wall of M.tb [10]. Furthermore, it is worth noting that M.tb utilizes fatty acids to produce essential metabolic intermediates closely related to its virulence [11]. M.tb is capable of incorporating fatty acids into phospholipids or utilizing them as a source of carbon for energy storage through their conversion into triglycerides. This conversion process has been linked to the promotion of drug resistance in M.tb [12, 13]. In response to hypoxia, M.tb within macrophages loaded with lipids undergoes a process of accumulating neutral lipids, which results in the loss of acid resistance and the development of antibiotic resistance [14]. Some virulence genes can facilitate the spread of M.tb. Fatty acid metabolism plays a significant role in enhancing the virulence and pathogenicity of M.tb. Nonetheless, the exact regulatory mechanism of fatty acid metabolism in M.tb is still unclear, and there is limited research on how mutations in fatty acid metabolism genes affect the transmission of M.tb. Therefore, further investigation is necessary to gain a better understanding of these aspects.

Whole genome sequencing (WGS) is a reliable tool for studying M.tb transmission. In this study, WGS was used to study the influence of fatty acid metabolism gene mutations on the transmission of M.tb in China. Specifically, the genomic cluster was used to represent the transmission of M.tb [15].

Results

Sample description

The M.tb isolates were classified according to the seven geographical regions of China. The vast majority of M.tb isolates from Eastern China (66.8%), Southern China (15.4%) and Central China (4.5%), as shown in Fig. 1. The analysis revealed that the majority of M.tb isolates belonged to lineage2 (n = 2744, 85.93%), followed by lineage4 (n = 439, 13.75%), and a smaller number of isolates belonged to lineage3 (n = 10, 0.31%). Most of the isolates belonged to sub-lineage2.2, while there were fewer isolates belonging to sub-lineage4.4 and sub-lineage4.5. The M.tb isolates were clustered into 499 groups, with sizes ranging from 2 to 108 isolates. Of these clusters, those containing 2 isolates of M.tb were defined as small clusters, those containing 3–5 isolates were defined as medium clusters, and those containing 6 or more isolates were defined as large clusters. There were 86 cross-regional clusters, ranging in size from 2 to 6 regions, as shown in Table 1. The phylogenetic tree of M.tb isolates was constructed as described in Fig. 2.

Fig. 1
figure 1

Distribution of 3197 isolates of Mycobacterium tuberculosis in seven natural geographical regions of China

Full size image
Table 1 Characteristics of Mycobacterium tuberculosis in China
Full size table
Fig. 2
figure 2

Phylogenetic tree for the 3197 Mycobacterium tuberculosis isolates from China

Full size image

The effect of mutations in fatty acid metabolism genes on clustering

After excluding positions with a mutation frequency lower than 0.01, we analyzed 73 mutation positions. During the comparison between clustered and non-clustered isolates, we observed significant differences (P < 0.05) for 43 mutation positions in fatty acid metabolism genes in the univariate analysis, as detailed in Supplement Table 1. Following univariate analysis, 73 mutation positions variables were selected for multivariate regression. To correct for possible confounding factors, we used the lineage and geographical location of M.tb as covariates in addition to the 73 mutation positions in fatty acid metabolism genes. Finally, five mutation positions of fatty acid metabolism genes with significant influence on clustering were determined (P < 0.05), as shown in Table 2. Among these, three mutation positions of fatty acid metabolism genes were identified as risk factors for clustering. The mutations were gca (136,605, 317G > C, Arg106Pro; OR, 22.144; 95% CI, 2.591-189.272), ogt(1,477,346, 286G > C ,Gly96Arg; OR, 3.893; 95%CI, 1.432–10.583), and rpsA (1,834,776, 1235 C > T, Ala412Val; OR, 3.674; 95% CI, 1.217–11.091), respectively.

Table 2 Analysis of the effect of fatty acid metabolism gene mutations on clustering
Full size table

Effects of mutations in fatty acid metabolism genes on clustering of lineage2

Positions with mutation in genes involved in fatty acid metabolism in lineage2 frequency less than 0.01 were removed. We analyzed 55 mutation positions, and 18 fatty acid metabolism gene mutation positions showed statistically significant differences between clustered and non-clustered isolates (P < 0.05), as detailed in Supplementary Table 2. Following univariate analysis, 55 mutation positions of fatty acid metabolism genes were analyzed by multivariate regression. In order to correct the possible confounding factors, we used the geographical location of M.tb as a covariate in addition to the 55 mutation positions in fatty acid metabolism genes. The results showed that mutations in six fatty acid metabolism gene positions were significantly associated with the clustered isolates of lineage2 (P < 0.05), see Table 3. Among these three mutation positions were identified as risk factors for clustering, including ogt (1,477,346, 286G > C, Gly96Arg; OR, 3.952; 95% CI, 1.453–10.749), rpsA (1,834,776, 1235 C > T, Ala412Val; OR,3.636; 95% CI,1.204–10.982), and gca (136,605, 317G > C, Arg106Pro; OR, 22.789; 95% CI, 2.669-194.569).

Table 3 Analysis of the effect of fatty acid metabolism gene mutations on clustering of lineage 2
Full size table

Effects of mutations in fatty acid metabolism genes on clustering of lineage4

The analysis focused on 33 fatty acid metabolism gene mutation positions in lineage4, which were selected by excluding those with frequencies lower than 0.01. In the comparison between clustered and non-clustered isolates of lineage4, the difference in the mutation of two fatty acid metabolism gene positions was statistically significant (P < 0.05). Results can be found in Supplementary Table 3. Following univariate analysis, 33 fatty acid metabolism gene mutation positions were included in a multivariate regression analysis. However, we included the geographical location of M.tb as a covariate in our analysis in order to control for possible confounding effects. The results showed that there was no risk factor for the clustered isolates of lineage4, see Table 4.

Table 4 Analysis of the effect of fatty acid metabolism gene mutations on clustering of lineage 4
Full size table

The effect of mutations in fatty acid metabolism genes on the cross-regional transmission of M.tb

After screening out the positions with clustering mutation frequency less than 0.01, 61 mutation positions of fatty acid metabolism genes were analyzed. In comparison between the cross-regional and non-cross-regional clusters, 26 fatty acid metabolism gene mutation positions showed significant differences (P < 0.05), as detailed in Supplementary Table 4. Following univariate analysis, 61 mutation positions were included in multiple regression analysis, and we also included the lineage as covariate to correct for potential confounding factors. The results showed that five mutation positions of fatty acid metabolism genes had a significant influence on regional factors (P < 0.05), see Table 5. Among these, mutation position of arsA(3,001,498) was identified as cross-regional risk factors (885 C > G, Thr295Thr; OR, 6.278; 95% CI, 2.508–15.711). Notably, the arsA was synonymous mutations.

Table 5 Analysis of the influence of fatty acid metabolism gene mutations on cross-regional
Full size table

Effects of mutations in fatty acid metabolism genes on cluster size of M.tb

A total of 61 mutation positions of fatty acid metabolism genes were analyzed. The results showed that 31 mutation positions were significantly associated with cluster size (P < 0.05). Among these, 20 mutation positions were found to be positively related to cluster size. Notably, seven of these mutation positions were synonymous, including fgd1 (491,742, 960T > C, Phe320Phe), fadB (957,117, 825T > C, Asp275Asp), fadH (1,306,259, 1968T > C, Ala656Ala), rpsA (1,834,177, 636 A > C, Arg212Arg), fadD15 (2,449,629, 1470G > A, Gln490Gln), fas (2,847,281, 2052T > C, Asp684Asp), and agpS (3,476,350, 612 C > T, Ser204Ser). For further details refer to Fig. 3.

Fig. 3
figure 3

Correlation analysis of fatty acid metabolism gene mutation positions and clusters

Full size image

Discussion

Fatty acid metabolism plays a crucial role in the growth of M.tb. To investigate the impact mutations of fatty acid metabolism gene mutations on the spread of TB in China, we analyzed 3107 isolates of M.tb and 83 fatty acid metabolism genes. In China, most of the M.tb isolates belonged to lineage2 (Beijing lineage), followed by lineage 4 (European lineage), and lineage3 (South Asia lineage).Most of the clustered isolates (n = 1463,91.67%) also belonged to lineage2, which indicated that the main isolates of transmission belonged to lineage2 in China.

Based on our findings, we observed a missense mutation (317G > C, Arg106Pro) at position 136,605 of gca (Rv0112), and another missense mutation (1235 C > T, Ala412Val) at position 1,834,776 of rpsA (Rv1630). These mutations have been associated with increased risk of transmission of M.tb, particularly within lineage 2, and are also correlated with cluster size. Although some functions of gca remain unclear, they may be associated with the transport of the M.tb cell membrane and the synthesis of the cell wall, both of which play critical roles in the pathogenesis of TB. Further research is needed to fully understand the mechanism by which this mutation promotes transmission. RpsA (Rv1630) is the largest 30 S protein in the ribosome and plays a crucial role in translation. Mutations or deletions of rpsA can have a significant impact on the growth and metabolism of M.tb [16,17,18]. A missense mutation (c.1235 C > T p.Ala412Val) has been identified at position 1,834,776 of rpsA. This mutation promotes the spread of TB isolates and lineage2 isolates and is associated with cluster size. Interestingly, both the Beijing isolate of M.tb and multidrug-resistant isolates exhibit two non-synonymous single nucleotide polymorphisms in the ogt gene [19,20,21]. The researchers hypothesized that these mutations in ogt (Rv1316c) may contribute to the successful global distribution of these isolates, which is consistent with our findings. Our results revealed a missense mutation (286 g > C, Gly 96 arg) at position 1,477,346 of the ogt gene. The ogt gene encodes an enzyme called N-acetylglucosamine (O-GlcNAc) transferase, which is a glycosyltransferase responsible for catalyzing the addition of O-GlcNAc modification onto specific serine or threonine residues of proteins. O-GlcNAc transferase may play a role in regulating M.tb growth, adaptability, and pathogenicity by modifying and affecting key M.tb proteins. This missense mutation potentially promotes the transmission of M.tb isolates, including lineage 2 isolates, and could have implications for M.tb metabolism, cell wall synthesis, drug resistance, and other characteristics [22, 23]. In our study, we did not find any mutations in fatty acid metabolism genes that had an impact on the transmission of lineage 4. This could be partially attributed to the fact that our sample size may have been insufficient to capture rare lineage 4 strains or related mutations, and a larger sample set might provide more accurate results.

A synonymous mutation at position 3,001,498 of arsA (Rv2684) (885 C > G, Thr295Thr) has been found to affect the transmission of isolates across different regions. The expression of arsA allows M.tb to adapt to different environments within the host’s body. Specifically, arsA helps the bacterium to evade the host immune response [2, 24].

In addition, our results confirmed that both synonymous and non-synonymous mutations can affect the transmission of M.tb, indicating that synonymous mutations in fatty acid metabolism of M.tb are not all neutral mutations, which is consistent with the result that synonymous mutations in yeast genes studied by Xukang Shen are mostly strong non-neutral mutations [25].

Conclusion

The results of this study suggest that mutations in fatty acid metabolism genes may increase the transmission risk of M.tb, which highlights the need for further investigation into the effects of these mutations on M.tb control and dissemination. These findings provide valuable insights into the therapy of TB.

Method

Sample collection

A total of 1550 M.tb culture-positive cases were collected from two medical institutions from 2011 to 2018 in China: Shandong Public Health Clinical Research Center (SPHCC) and Weifang Respiratory Clinical Hospital (WRCH). All samples were collected anonymously and informed consent was not required. Our research was approved by the Ethics Committee of Shandong Provincial Hospital, which is affiliated with Shandong First Medical University.

DNA extraction and sequencing

Genomic DNA from 1447 isolates was extracted with Cetyltrimethylammonium Bromide (CTAB) and underwent quality control (QC). The Illumina HiSeq 4000 system was used to sequence the genomes [26], and the sequence data were deposited in the National Center for Biotechnology Information (NCBI) under BioProject PRJNA1002108. In addition, 1755 isolates of M.tb from 23 provinces, 4 municipalities, and 5 autonomous regions in China were included in this study [27,28,29,30,31,32,33,34]. See Supplementary Tables 5–6 for the sample number. A total of 3202 genomes were analyzed, and M.tb H37Rv was used as the reference genome sequence.

Single nucleotide polymorphism (SNP) analysis

To map the sequencing reads to the standard isolate H37Rv, the BWA Mem algorithm (version 0.7.17-r1188) was used. We only included samples with a coverage rate of 98% or higher and a minimum depth of at least 20% [35]. Variant calling was performed using Samclip (version 0.4.0) and SAMtools (version 1.15), and the resulting variants were further filtered by Free Bayes (version 1.3.2) and Bcftools (version 1.15.1). We excluded Single nucleotide polymorphisms (SNPs) located in repeat regions, such as polymorphic GC-rich sequences (PE/PPE genes) and direct repeat SNPs, as well as repeat bases identified by Tandem Repeat Finder (version 4.09) and RepeatMask (version 4.1.2-P1) [36, 37]. Finally, the SNP was annotated with SnpEff v 4.1 l, and the result was obtained with Python programming language [38].

Prediction of drug resistance

To identify drug resistance mutations, we compared known indels and SNPs using TBProfiler (version 2.8.12) and the tuberculosis database (TBDB) [39, 40]. We then searched for genotypic markers of drug resistance mutations in both first-line drugs (such as isoniazid, rifampicin, pyrazinamide, ethambutol, and streptomycin) and second-line drugs (such as ethionamide, quinolones, amikacin, capreomycin, and kanamycin), using a set of genetic polymorphisms. Mutations that were not correlated with phenotypic drug resistance were excluded as markers of genetic drug resistance [41]. For more information about the mutations detected as molecular resistance predictions in 3202 isolates, please refer to Supplementary Table 7.

Phylogenetic analysis

The isolates were divided into different lineages according to Coll et al. [42] (Supplementary Tables 5–6). The maximum likelihood phylogenetic tree construction was performed using IQ-TREE (verdion1.6.12) using the JC nucleotide substitution model, the gamma model of rate heterogeneity, and 100 bootstrap replicates [43]. M carneti CIPT140010059 was considered an outlier, and five isolates belonging to two lineages were excluded. The phylogenetic tree was visualized by iTOL (https://itol.Embl.De/). However, isolates of lineage1were excluded from further analysis because of their small number. Therefore, a total of 3193 isolates were included in the final analysis.

Propagation analysis

Cluster analysis was used to study the effect of fatty acid metabolism gene mutation on the transmission of M.tb. Clustering was defined as a group of isolates with less than 10 SNPs among each other (see Supplementary Table 8). To study the regional variations, the geographical location of the isolates in China was divided into seven natural regions. Then, the clusters were classified as cross-regional clusters or non-cross-regional clusters. The cross-regional cluster means that the strains in the cluster come from two or more different regions.

Acquisition of fatty acid metabolic genes

According to the NCBI database, a total of 83 fatty acid metabolism genes were obtained. Mutations in genes involved in fatty acid metabolism were done by bcftools (version 1.15.1) with an included filter parameter ‘FMT/GT="1/1” && QUAL > = 100 && FMT/DP > = 10 && (FMT/AO)/(FMT/DP) > = 0’. The results were shown in Supplementary Table 9.

Statistical analysis

The data are presented as a number (percent). The positions with mutation frequency < 0.01 in fatty acid metabolism genes were excluded from the analysis [44]. SPSS version 26 was used for statistical analysis. The comparison of categorical variables was done using the Pearson’s chi-square test or Fisher exact test as appropriate between clustered and non-clustered, as well as cross-regional and non-cross-regional clusters. Variables with univariate analysis were included in the binary logistic regression model for multivariate analysis. To analyze the effect of fatty acid metabolism gene mutations on cluster size, the rank correlation analysis of Spearman was carried out by using R version 4.1.0. All reported statistical tests were 2-sided, and P values < 0.05 were considered statistically significant.

Data Availability

The whole genome sequences have been submitted to the NCBI under the accession.

number PRJNA1002108.

Change history

Abbreviations

CTAB:

Cetyltrimethylammonium Bromide

M.tb:

Mycobacterium tuberculosis

QC:

Quality control

SNP:

Single nucleotide polymorphism

SNPs:

Single nucleotide polymorphisms

SPHCC:

Shandong Public Health Clinical Research Center

TB:

Tuberculosis

TBDB:

Tuberculosis database

WGS:

Whole genome sequencing

WRCH:

Weifang Respiratory Clinical Hospital

References

  1. Rylance J, Pai M, Lienhardt C, Garner P. Priorities for Tuberculosis research: a systematic review. Lancet Infect Dis. 2010;10:886–92.

    Article  PubMed  Google Scholar 

  2. Furin J, Cox H, Pai M, Tuberculosis. The Lancet. 2019;393:1642–56.

    Article  Google Scholar 

  3. Linh NN, Viney K, Gegia M, Falzon D, Glaziou P, Floyd K, et al. World Health Organization treatment outcome definitions for Tuberculosis: 2021 update. Eur Respir J. 2021;58:2100804.

    Article  PubMed  Google Scholar 

  4. Atamna-Mawassi H, HuBerman-Samuel M, Hershcovitz S, Karny-Epstein N, Kola A, Cortés LEL, et al. Interventions to reduce Infections caused by multidrug resistant Enterobacteriaceae (MDR-E): a systematic review and meta-analysis. J Infect. 2021;83:156–66.

    Article  PubMed  CAS  Google Scholar 

  5. Van Rie A, Enarson D. XDR Tuberculosis: an indicator of public-health negligence. Lancet. 2006;368:1554–6.

    Article  PubMed  Google Scholar 

  6. World Health Organization. Global Tuberculosis report 2016. Geneva: World Health Organization; 2016.

    Google Scholar 

  7. Singh V, Mani I, Chaudhary DK, Somvanshi P. The β-ketoacyl-ACP synthase from Mycobacterium tuberculosis as potential drug targets. Curr Med Chem. 2011;18:1318–24.

    Article  PubMed  CAS  Google Scholar 

  8. Nazarova EV, Montague CR, Huang L, La T, Russell D, VanderVen BC. The genetic requirements of fatty acid import by Mycobacterium tuberculosis within macrophages. Elife. 2019;8:e43621.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Laval T, Chaumont L, Demangel C. Not too fat to fight: the emerging role of macrophage fatty acid metabolism in immunity to Mycobacterium tuberculosis. Immunol Rev. 2021;301:84–97.

    Article  PubMed  CAS  Google Scholar 

  10. Kinsella RJ, Fitzpatrick DA, Creevey CJ, McInerney JO. Fatty acid biosynthesis in Mycobacterium tuberculosis: lateral gene transfer, adaptive evolution, and gene duplication. Proc Natl Acad Sci USA. 2003;100:10320–5.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  11. Salaemae W, Azhar A, Booker GW, Polyak SW. Biotin biosynthesis in Mycobacterium tuberculosis: physiology, biochemistry and molecular intervention. Protein Cell. 2011;2:691–5.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  12. Grimes KD, Aldrich CC. A high-throughput screening fluorescence polarization assay for fatty acid adenylating enzymes in Mycobacterium tuberculosis. Anal Biochem. 2011;417:264–73.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  13. Sacco E, Covarrubias AS, O’Hare HM, Carroll P, Eynard N, Jones TA, et al. The missing piece of the type II fatty acid synthase system from Mycobacterium tuberculosis. Proc Natl Acad Sci U S A. 2007;104:14628–33.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  14. McKinney JD, Höner zu Bentrup K, Muñoz-Elías EJ, Miczak A, Chen B, Chan WT, et al. Persistence of Mycobacterium tuberculosis in macrophages and mice requires the glyoxylate shunt enzyme isocitrate lyase. Nature. 2000;406:735–8.

    Article  PubMed  CAS  Google Scholar 

  15. Alonso Rodríguez N, Martínez Lirola M, Chaves F, lñigo J, Herranz M, Ritacco V, et al. Differences in the robustness of clusters involving the Mycobacterium tuberculosis strains most frequently isolated from immigrant cases in Madrid. Clin Microbiol Infect. 2010;16:1544–54.

    Article  PubMed  Google Scholar 

  16. Akhmetova A, Kozhamkulov U, Bismilda V, Chingissova L, Abildaev T, Dymova M, et al. Mutations in the < I > pnc A and < I > rps a genes among 77 < I > Mycobacterium tuberculosis isolates in Kazakhstan. int j Tuberc lung dis. 2015;19:179–84.

    Article  PubMed  CAS  Google Scholar 

  17. Vallejos-Sánchez K, Lopez JM, Antiparra R, Toscano E, Saavedra H, Kirwan DE, et al. Mycobacterium tuberculosis ribosomal protein S1 (RpsA) and variants with truncated C-terminal end show absence of interaction with pyrazinoic acid. Sci Rep. 2020;10:8356.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Singh A, Somvanshi P, Grover A. Pyrazinamide drug resistance in RpsA mutant (∆438A) of Mycobacterium tuberculosis: Dynamics of essential motions and free-energy landscape analysis. J Cell Biochem. 2018. https://0-doi-org.brum.beds.ac.uk/10.1002/jcb.28013

    Article  PubMed  PubMed Central  Google Scholar 

  19. Ebrahimi-Rad M, Bifani P, Martin C, Kremer K, Samper S, Rauzier J, et al. Mutations in putative mutator genes of Mycobacterium tuberculosis strains of the W-Beijing family. Emerg Infect Dis. 2003;9:838–45.

    Article  PubMed  Google Scholar 

  20. Olano J, López B, Reyes A, del Pilar Lemos M, Correa N, Del Portillo P, et al. Mutations in DNA repair genes are associated with the Haarlem lineage of Mycobacterium tuberculosis independently of their antibiotic resistance. Tuberculosis. 2007;87:502–8.

    Article  PubMed  CAS  Google Scholar 

  21. Mestre O, Luo T, Dos Vultos T, Kremer K, Murray A, Namouchi A, et al. Phylogeny of Mycobacterium tuberculosis Beijing strains constructed from polymorphisms in genes involved in DNA replication, recombination and repair. PLoS ONE. 2011;6:e16020.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  22. Miggiano R, Perugino G, Ciaramella M, Serpe M, Rejman D, Páv O, et al. Crystal structure of Mycobacterium tuberculosis O 6-methylguanine-DNA methyltransferase protein clusters assembled on to damaged DNA. Biochem J. 2016;473:123–33.

    Article  PubMed  CAS  Google Scholar 

  23. Rad ME, Bifani P, Martin C, Kremer K, Samper S, Rauzier J et al. Mutations in putative mutator genes of Mycobacterium tuberculosis strains of the W-Beijing Family. Emerg Infect Dis. 2003;9.

  24. Hu K, Jordan AT, Zhang S, Dhabaria A, Kovach A, Rangel MV et al. Characterization of guided entry of tail-anchored proteins 3 homologues in Mycobacterium tuberculosis. J Bacteriol. 2019;201.

  25. Shen X, Song S, Li C, Zhang J. Synonymous mutations in representative yeast genes are mostly strongly non-neutral. Nature. 2022;606:725–31.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  26. Foox J, Tighe SW, Nicolet CM, Zook JM, Byrska-Bishop M, Clarke WE, et al. Performance assessment of DNA sequencing platforms in the ABRF Next-Generation sequencing study. Nat Biotechnol. 2021;39:1129–40.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  27. Chen X, He G, Wang S, Lin S, Chen J, Zhang W. Evaluation of whole-genome sequence method to Diagnose Resistance of 13 anti-tuberculosis Drugs and characterize resistance genes in clinical Multi-drug Resistance Mycobacterium tuberculosis isolates from China. Front Microbiol. 2019;10:1741.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Yang C, Luo T, Shen X, Wu J, Gan M, Xu P, et al. Transmission of multidrug-resistant Mycobacterium tuberculosis in Shanghai, China: a retrospective observational study using whole-genome sequencing and epidemiological investigation. Lancet Infect Dis. 2017;17:275–84.

    Article  PubMed  CAS  Google Scholar 

  29. Zhang H, Li D, Zhao L, Fleming J, Lin N, Wang T, et al. Genome sequencing of 161 Mycobacterium tuberculosis isolates from China identifies genes and intergenic regions associated with drug resistance. Nat Genet. 2013;45:1255–60.

    Article  PubMed  CAS  Google Scholar 

  30. Hicks ND, Yang J, Zhang X, Zhao B, Grad YH, Liu L, et al. Clinically prevalent mutations in Mycobacterium tuberculosis alter propionate metabolism and mediate multidrug tolerance. Nat Microbiol. 2018;3:1032–42.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  31. Liu Q, Ma A, Wei L, Pang Y, Wu B, Luo T, et al. China’s Tuberculosis epidemic stems from historical expansion of four strains of Mycobacterium tuberculosis. Nat Ecol Evol. 2018;2:1982–92.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Huang H, Ding N, Yang T, Li C, Jia X, Wang G, et al. Cross-sectional whole-genome sequencing and epidemiological study of Multidrug-resistant Mycobacterium tuberculosis in China. Clin Infect Dis. 2019;69:405–13.

    Article  PubMed  CAS  Google Scholar 

  33. Luo T, Comas I, Luo D, Lu B, Wu J, Wei L, et al. Southern East Asian origin and coexpansion of Mycobacterium tuberculosis Beijing family with Han Chinese. Proc Natl Acad Sci USA. 2015;112:8136–41.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  34. Jiang Q, Liu Q, Ji L, Li J, Zeng Y, Meng L, et al. Citywide transmission of Multidrug-resistant Tuberculosis under China’s Rapid Urbanization: a Retrospective Population-based genomic spatial epidemiological study. Clin Infect Dis. 2020;71:142–51.

    Article  PubMed  Google Scholar 

  35. Jung Y, Han D, BWA-MEME. BWA-MEM emulated with a machine learning approach. Bioinformatics. 2022;:btac137.

  36. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25:2078–9.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Liu F, Zhang Y, Zhang L, Li Z, Fang Q, Gao R, et al. Systematic comparative analysis of single-nucleotide variant detection methods from single-cell RNA sequencing data. Genome Biol. 2019;20:242.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  38. Cingolani P, Platts A, Wang LL, Coon M, Nguyen T, Wang L, et al. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly (Austin). 2012;6:80–92.

    Article  PubMed  CAS  Google Scholar 

  39. Tayyab N, Zaman G, Satti L, Ikram A, Gardezi AH, Khadim MT. Direct susceptibility testing on MGIT 960 TB system: a Rapid Method for detection of drug resistant Tuberculosis. J Coll Physicians Surg Pak. 2018;28:590–3.

    Article  PubMed  Google Scholar 

  40. Phelan JE, O’Sullivan DM, Machado D, Ramos J, Oppong YEA, Campino S, et al. Integrating informatics tools and portable sequencing technology for rapid detection of resistance to anti-tuberculous Drugs. Genome Med. 2019;11:41.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Banerjee R, Allen J, Westenhouse J, Oh P, Elms W, Desmond E, et al. Extensively drug-resistant Tuberculosis in California, 1993–2006. Clin Infect Dis. 2008;47:450–7.

    Article  PubMed  Google Scholar 

  42. Coll F, McNerney R, Guerra-Assunção JA, Glynn JR, Perdigão J, Viveiros M, et al. A robust SNP barcode for typing Mycobacterium tuberculosis complex strains. Nat Commun. 2014;5:4812.

    Article  PubMed  CAS  Google Scholar 

  43. Nguyen L-T, Schmidt HA, von Haeseler A, Minh BQ. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol. 2015;32:268–74.

    Article  PubMed  CAS  Google Scholar 

  44. Farhat MR, Freschi L, Calderon R, Ioerger T, Snyder M, Meehan CJ, et al. GWAS for quantitative resistance phenotypes in Mycobacterium tuberculosis reveals resistance genes and regulatory regions. Nat Commun. 2019;10:2128.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

Download references

Acknowledgements

We thank Shandong Provincial Hospital, Shandong Provincial Chest Hospital, 13 municipal-level and 21 county-level local health departments for drug susceptibility data, demographic, and clinical data.

Funding

This work was supported by the Department of Science & Technology of Shandong Province (CN) (No.2007GG30002033; No.2017GSF218052) and the Jinan Science and Technology Bureau (CN) (No.201704100). The funding body/bodies did not provide any assistance in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Author information

Authors and Affiliations

  1. Deartment of Chinese Medicine Integrated with Western Medicine, College of Traditional Chinese Medicine, Shandong University of Traditional Chinese Medicine, 16369 Jingshi Road, Lixia District, Jinan, 250355, Shandong, People’s Republic of China

    Yameng Li, Tingting Wang, Yingying Li, Yawei Hou & Huaichen Li

  2. Artificial Intelligence Institute, Qilu University of Technology (Shandong Academy of Sciences), Jinan, 250011, Shandong, People’s Republic of China

    Xianglong Kong

  3. Department of Respiratory and Critical Care Medicine, The Third Affiliated Hospital of Shandong First Medical University (Affiliated Hospital of Shandong Academy of Medical Sciences), Jinan, 250031, Shandong, People’s Republic of China

    Yifan Li

  4. Department of Respiratory and Critical Care Medicine, Shandong Provincial Hospital, Shandong University, Shandong Provincial Hospital Affiliated to Shandong First Medical University, Jingwuweiqi Road, Huaiyin District, Jinan, 250021, Shandong, People’s Republic of China

    Ningning Tao, Qiqi An, Yao Liu & Huaichen Li

  5. Shandong First Medical University & Shandong Academy of Medical Sciences, Jinan, 250117, Shandong, People’s Republic of China

    Xuehan Zhu, Qilin Han & Yuzhen Zhang

Authors
  1. Yameng Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Xianglong Kong
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Yifan Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Ningning Tao
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Tingting Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Yingying Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Yawei Hou
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Xuehan Zhu
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Qilin Han
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Yuzhen Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  11. Qiqi An
    View author publications

    You can also search for this author in PubMed Google Scholar

  12. Yao Liu
    View author publications

    You can also search for this author in PubMed Google Scholar

  13. Huaichen Li
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

HCL, YL, and YML participated in the study design. YL, HCL, YML, XLK, and YFL performed data collection and statistical analyses. YL, NNT, TTW, YYL, and YWH helped draft the manuscript. XHZ, QLH, YZZ, and QQA overviewed and supervised the project. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Yao Liu or Huaichen Li.

Ethics declarations

Ethics approval and consent to participate

This study complies with the Declaration of Helsinki, and was approved by the Ethics Committee of Shandong Provincial Hospital, affiliated with Shandong University (SPH) and the Ethics Committee of Shandong Provincial Chest Hospital (SPCH), which waived informed patient consent because all patient records and information were anonymized and deidentified before the analysis.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1: Supplement Table 1

Correlation analysis of fatty acid metabolism gene mutations between clustered and non-clustered isolates

Supplementary Material 2: Supplement Table 2

Correlation analysis of fatty acid metabolism gene mutations between clustered and non-clustered isolates of lineage2

Supplementary Material 3: Supplement Table 3

Correlation analysis of fatty acid metabolism gene mutations between clustered and non-clustered isolates of lineage4

Supplementary Material 4: Supplement Table 4

Correlation analysis of fatty acid metabolism gene mutations between cross-regional and non-cross-regional clusters

Supplementary Material 5: Supplement Table 5

Information about 1447 isolates of Mycobacterium tuberculosis

Supplementary Material 6: Supplement Table 6

Information about 1755 isolates of Mycobacterium tuberculosis

Supplementary Material 7: Supplement Table 7

Drug resistance of 3202 isolates of Mycobacterium tuberculosis

Supplementary Material 8: Supplement Table 8

Information on clustering of 3193 isolates of Mycobacterium tuberculosis

Supplementary Material 9: Supplement Table 9

Mutations of fatty acid metabolism genes in 3193 isolates of Mycobacterium tuberculosis

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Li, Y., Kong, X., Li, Y. et al. Association between fatty acid metabolism gene mutations and Mycobacterium tuberculosis transmission revealed by whole genome sequencing. BMC Microbiol 23, 379 (2023). https://0-doi-org.brum.beds.ac.uk/10.1186/s12866-023-03072-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12866-023-03072-9

Keywords

BMC Microbiology

ISSN: 1471-2180

Contact us