Novel Potential Diagnostic Serum Biomarkers of Metabolomics in Osteoarticular Tuberculosis Patients: A Preliminary Study
1 Medical School of Chinese People’s Liberation Army (PLA), Beijing, China
2 Department of Clinical Laboratory Medicine, The First Medical Center, Chinese People’s Liberation Army (PLA) General Hospital, Beijing, China
3 Department of Clinical Laboratory Medicine, The Eighth Medical Center, Chinese People’s Liberation Army (PLA) General Hospital, Beijing, China
*Correspondence: Chengbin Wang, wangcbin301@163.com; Hong Lei, leihong_hospital@126.comAbstract
Osteoarticular tuberculosis is one of the extrapulmonary tuberculosis, which is mainly caused by direct infection of Mycobacterium tuberculosis or secondary infection of tuberculosis in other parts. Due to the low specificity of the current detection method, it is leading to a high misdiagnosis rate and subsequently affecting the follow-up treatment and prognosis. Metabolomics is mainly used to study the changes of the body’s metabolites in different states, so it can serve as an important means in the discovery of disease-related metabolic biomarkers and the corresponding mechanism research. Liquid chromatography tandem mass spectrometry (LC-MS/MS) was used to detect and analyze metabolites in the serum with osteoarticular tuberculosis patients, disease controls, and healthy controls to find novel metabolic biomarkers that could be used in the diagnosis of osteoarticular tuberculosis. Our results showed that 68 differential metabolites (p<0.05, fold change>1.0) were obtained in osteoarticular tuberculosis serum after statistical analysis. Then, through the evaluation of diagnostic efficacy, PC[o-16:1(9Z)/18:0], PC[20:4(8Z,11Z,14Z,17Z)/18:0], PC[18:0/22:5(4Z,7Z,10Z,13Z,16Z)], SM(d18:1/20:0), and SM[d18:1/18:1(11Z)] were found as potential biomarkers with high diagnostic efficacy. Using bioinformatics analysis, we further found that these metabolites share many lipid metabolic signaling pathways, such as choline metabolism, sphingolipid signaling, retrograde endocannabinoid signaling, and sphingolipid and glycerophospholipid metabolism; these results suggest that lipid metabolism plays an important role in the pathological process of tuberculosis. This study can provide certain reference value for the study of metabolic biomarkers of osteoarticular tuberculosis and the mechanism of lipid metabolism in osteoarticular tuberculosis and even other tuberculosis diseases.
Introduction
Tuberculosis is a chronic infectious disease that is caused by Mycobacterium tuberculosis (M.tb) and takes the respiratory tract as the main route of transmission (Walzl et al., 2018). According to the Global Tuberculosis Report by the World Health Organization in 2021, 9.9 million new cases and 1.5 million deaths were caused by tuberculosis globally in 2020 (World Health Organization, 2021). Generally, tuberculosis is divided into two types: pulmonary tuberculosis (PTB) and extrapulmonary tuberculosis (EPTB). For the diagnosis of tuberculosis, clinical laboratory tests play an important role, although with the long culture time and the low positive rate, the culture of Mycobacterium tuberculosis is still the gold standard (Darton et al., 2017); the other common tests such as smear microscopy, purified protein derivative (PPD) test, interferon-gamma (IFN-γ) release assay (IGRA), M.tb nucleic acid test, and Xpert/MTB system are also widely used (Fan et al., 2012; Norbis et al., 2014; Fan et al., 2018).
Osteoarticular tuberculosis is a kind of extrapulmonary tuberculosis whose M.tb directly infects the bone and joint tissue or spreads to the bone and joint tissue from other parts, which accounts for about 3%–5% of the total incidence of EPTB (Norbis et al., 2014). The most common site for osteoarticular tuberculosis is spine, especially the thoracic and lumbar spine, and the pathological changes are mostly bone destruction and tuberculous granulation tissue formation. At the same time, tuberculous abscesses that do not show related inflammation are easily formed beside the bone, and in severe cases, sinus tracts may even be formed (Johansen et al., 2015). The early clinical manifestations of osteoarticular tuberculosis are atypical, and the specificity of laboratory tests and imaging examinations is low, resulting in a high rate of missed diagnosis for osteoarticular tuberculosis. Most patients are misdiagnosed during the advanced stage that affected the treatment and prognosis. Hence, there is an urgent need for accurate diagnosis methods of osteoarticular tuberculosis (Yi et al., 2018; Vinhaes et al., 2019).
Metabolomics is mainly used to study the changes of the body’s metabolites in different states (Zhang et al., 2018); the main methods of metabolomics contained nuclear magnetic resonance (NMR) and mass spectrometry (MS) (Zhou and Yin, 2019). Numerous metabolites have been successfully discovered as biomarkers for the diagnosis of various diseases. What is more, the metabolites which are searched by metabolomics can be elucidating the pathological or functional mechanisms by bioinformatics analysis (Wheelock et al., 2013). In the aspect of tuberculosis, there were many metabolomics studies of PTB, including secondary pulmonary tuberculosis (Luies et al., 2017), active pulmonary tuberculosis (Cho et al., 2020), and drug-resistant pulmonary tuberculosis (Tuyiringire et al., 2018; Chen et al., 2020); however, for osteoarticular tuberculosis, because of the low incidence of this disease, the related metabolomics studies are rare.
According to this situation, we used LC-MS/MS to detect serum metabolites in osteoarticular tuberculosis patients (TB Group), osteoarthritis patients (DC Group, including rheumatoid arthritis and ankylosing spondylitis), and healthy controls (HC Group) (Dudka et al., 2021), aiming to find potential markers for early and accurate diagnosis of osteoarticular tuberculosis. We further compared and analyzed serum metabolites of the above three groups and selected differential metabolites by the bioinformatics method to analyze the corresponding pathogenesis (Zhang et al., 2019).
Materials and Methods
Study Cohort
From November 2018 to November 2019, 30 serum samples of diagnosed osteoarticular tuberculosis patients (TB group), 30 serum samples of disease control containing rheumatoid arthritis patients and ankylosing spondylitis patients (DC group), and 30 serum samples of healthy control (HC group) were collected from the 1st and the 8th Medical Center of Chinese PLA General Hospital. The diagnosis of the TB group was based on the following criteria: (a) positive nucleic acid test of M.tb; (b) medical image (X-ray, CT scan, etc.) findings showed specific features of TB infection; (c) positive pathology diagnosis of TB in bone or joint specimens; and (d) effective response to antituberculosis treatments. The diagnosis of the DC group was based on the following criteria: (a) positive inflammatory protein test; (b) specific clinical manifestations (joint stiffness in the morning, bending change, etc.); (c) medical image (X-ray, CT scan, etc.) findings showed specific features; and (d) effective response to hormone treatments (glucocorticoid, etc.). The HC group included adults without any disease clinical diagnosis, and all the tests were negative or normal. Patients with any diagnosis of cancer, metabolic disease, autoimmunity disease, immunodeficiency disease, and other pathogen infections were excluded from this study. Patients who have other organs tuberculosis (e.g., pulmonary tuberculosis or other extrapulmonary tuberculosis) were excluded as well.
This research was carried out in strict accordance with the declaration of Helsinki and approved by the Ethics Committee of Chinese PLA General Hospital. All participants signed an informed consent and gave their permission to use their blood samples for this study.
For each patient, 5 ml peripheral blood was drawn under a vacuum vessel containing separation gel in the morning before any treatments. After blood coagulation and 2,370g centrifugation, the serum was divided into several EP tubes and stored in a -80°C refrigerator for subsequent metabolomics analysis.
Statistical Analysis
The chi-square test was used for the analysis of characteristics of the study participants, and the Kruskal–Wallis H test was used to determine the differences between groups. The data of metabolomics were normalized using Progenesis QI data analysis software (Nonlinear Dynamics, Newcastle, UK). SIMCA 14.1 software was used to analyze the metabolites using the orthogonal partial least squares (OPLS) model. MetaboAnalyst (http://www.metaboanalyst.ca/) was used to analyze the related pathways of specific metabolites. R2X (the interpretability of the model for the categorical variable X) was obtained after cross-validation; R2Y (the interpretability of the model for the categorical variable Y) and Q2 (predictability of the model) were obtained after cross-validation to judge the validity of the model. The Variable Importance in the Projection (VIP) value and the p value of t-test were used for evaluating the difference metabolites between groups. Final results were shown with scatter plots, trend charts, and the receiver operator characteristic curve (ROC curve), and analyses of the AUC, sensitivity, and specificity of each different metabolites were made by GraphPad Prism 6.0 software. The establishment of diagnostic models which contained metabolic biomarkers was made by MedCalc software, including several different metabolites combination and its statistical analysis; meanwhile, logistic regression and ROC curve analysis were used for the establishment of diagnostic models.
Results
Basic Data Preprocessing
There was no statistical difference in the age and gender between TB, DC, and HC groups (p>0.05). The positive rates of clinical laboratory tests and medical imaging features are shown in Table 1 . After LC-MS/MS analysis and peak alignment, picking, and normalization of raw data, metabolites were obtained, while the QC results of three patterns showed good reproducibility, which indicated that the results were credible ( Supplementary Files ).
| TB Group (n = 30) | DC Group (n = 30) | HC Group (n = 30) | p value | |
|---|---|---|---|---|
| Age (median, IQR) | 46 (33–61) | 53.5 (31.5–63.5) | 48.5 (41–55) | >0.05 |
| Gender (male/female) | 16/14 | 14/16 | 17/13 | >0.05 |
| Xpert test positive no. (%) | 17 (56.67%) | / | / | / |
| TB nucleic acid test positive no. (%) | 13 (43.33%) | / | / | / |
| TB antibody test positive no. (%) | 9 (30.0%) | / | / | / |
| IGRA positive no. (%) | 25 (83.3%) | / | / | / |
| Cultivate positive no. (%) | 2 (6.67%) | / | / | / |
| Imaging features positive no. (%) | 24 (80.00%) | / | / | / |
| Pathology positive no. (%) | 17 (56.67%) | / | / | / |
KEGG Enrichment Analysis
For the differential metabolites screened out in the previous analysis, we conducted a bioinformatics analysis on the differential metabolites. We imported these metabolites one by one into the KEGG database for signal pathway analysis and performed further enrichment and statistical analysis of these results. After KEGG enrichment analysis, the main metabolic pathways in these differential metabolites are necroptosis, choline metabolism, sphingolipid signaling, retrograde endocannabinoid signaling, sphingolipid metabolism, and glycerophospholipid metabolism ( Figure 4 ). According to the results of bioinformatics, it can be seen that the main signal pathways are concentrated in cell metabolism and lipid metabolism, which is in good agreement with the currently known pathogenic mechanisms of Mycobacterium tuberculosis.
Discussion
Tuberculosis is caused by Mycobacterium tuberculosis, one of the most widespread infectious diseases worldwide (Procop, 2016), which is always divided into two types: pulmonary tuberculosis (PTB) and extrapulmonary tuberculosis (EPTB). Because of its lower incidence rate, EPTB does not attract enough attention compared to PTB. Osteoarticular tuberculosis is one of EPTB, whose M.tb directly infects the bone and joint tissue or spreads to the bone and joint tissue from other parts, and it is hard to distinguish with rheumatic immune arthritis and bone tumor at its early stage based on current clinical laboratory and clinical imaging methods (Fan et al., 2018). Metabolomics is one of the omics proposed together with proteomics and transcriptomics in recent years (Dutta et al., 2020). It is defined as the study of the complete set of metabolites inside cells, tissues, organs, and biological fluids. It represents a major and rapidly evolving component of systems biology—a new integrative approach to deciphering the complexity of biological systems (Preez et al., 2017). The metabolite changes associated with the specific phenotype being investigated may be classified as characteristics of the perturbation, which, in the context of investigating a disease, could be used toward better disease characterization, diagnostics, treatment, and other clinical applications. The methods of metabolomics include nuclear magnetic resonance (NMR) and mass spectrometry (MS) (Laíns et al., 2019). NMR is appropriate for the detection of all hydrogen-containing compounds, since it determines the magnetic resonance of nuclei in a molecule, and it is considered an unbiased, robust, reproducible, non-destructive, and selective analytical platform, which requires almost no sample pretreatment. On the other side, however, NMR has a low sensitivity and is short of an analyte separation element (Shin et al., 2011). MS is defined as the process of forming gaseous ions, with or without fragmentation, which are then characterized by their m/z ratios and respective relative abundances. Direct MS infusion is a high-throughput method, requiring short time for each sample analysis, and has been applied successfully in metabolomics studies, but it is not preferred for the analyses of complex biological samples such as blood and urine due to matrix interference (Schoeman et al., 2012). For TB, there are also more metabolomics studies focusing on colony culture, sputum specimens, blood specimens, urine specimens, tissue specimens, etc. (Parida and Kaufmann, 2010)
Zhou et al. (2015) and Albors-Vaquer et al. (2020) used NMR to detect differences in the expression of metabolites in the serum of PTB, healthy adults, lung-related benign lesions, and lung cancers; the results showed that some amino acids have been changed among these groups, such as alanine, lysine, glutamate, glutamine, ketone bodies, lactate, and pyruvate. Mendes Rˆego et al. (Rêgo et al., 2021) used MS to detect serum metabolites among drug-sensitive tuberculosis and drug-resistant tuberculosis; there were also some amino acids that were changed such as isoleucine, proline, hercynite, betaine, and pantothenic acid. On another aspect, Vrieling et al. (2019), Luo et al. (2020) and Chen et al. (2021) used MS to detect serum metabolites among PTB and healthy control, tuberculous pleuritis and malignancy, and PTB and PTB with type 2 diabetes, respectively. The results of these studies also showed some changes in amino acids, phospholipids, sphingolipids, etc., which are similar to other studies. According to other studies, amino acid levels such as alanine, lysine, glutamate, and glutamine in patients with tuberculosis infection have decreased significantly, which may be related to the uptake of glutamate, glutamate, and alanine by Mycobacterium tuberculosis for their corresponding life activities (Harth and Horwitz, 2003; Agapova et al., 2019). The body’s immune process for tuberculosis is mainly the activation of T cells, which also changes the law of glucose metabolism. The disease process affects the metabolism of related immune cells and thus affects the secretion of cytokines such as IFN-γ, further leading to a decrease in immune effect (Lande et al., 2003).
According to the characteristics of Mycobacterium tuberculosis, in addition to the cell membrane and peptidoglycan of ordinary bacteria, there are also a large number of lipids and carbohydrates on its surface. These components have strong biological activity on eukaryotic cells, thus suggesting that it has a strong relationship with pathogenicity. For the mechanism of tuberculosis, lipid metabolism and lipid effector molecules play a vital role, such as regulating the production of cytokines, scavenging oxygen free radicals, and producing granulation inflammation and mitochondrial toxicity (Walpole et al., 2018). These related pathogenic components have also been reflected in other studies. Our study focuses on the blood specimens of osteoarticular tuberculosis, compared with PTB; this type of research is currently rarely reported. The results of our study showed that compared with healthy adults and rheumatoid arthritis patients, patients with osteoarticular tuberculosis also have significant changes in serum amino acid and lipid metabolism. Most amino acids such as glycine, L-arginine, and L-alanine were downregulated in osteoarticular tuberculosis patients’ serum; in contrast, most lipids such as phosphatidylcholine, sphingomyelin, and phosphatidylethanolamine were upregulated in osteoarticular tuberculosis patients’ serum.
Further, our study selected five potential serum metabolite biomarkers, namely, PC(o-16:1(9Z)/18:0), PC(20:4(8Z,11Z,14Z,17Z)/18:0), PC(18:0/22:5(4Z,7Z,10Z,13Z,16Z)), SM(d18:1/20:0), and SM(d18:1/18:1(11Z)), which belong to phosphatidylcholine (PC) and sphingomyelin (SM). On the aspect of KEGG enrichments, necroptosis, choline metabolism, sphingolipid signaling, retrograde endocannabinoid signaling, sphingolipid metabolism, and glycerophospholipid metabolism are maybe the main pathways. In most cases, these metabolites are related to lipid accumulation and obesity in the body (Wahl et al., 2012; Reinehr et al., 2015; Hellmuth et al., 2016). At the same time, some other pathological processes that cause enhanced lipid metabolism can also significantly increase these metabolites; for example, Wu et al. (2018) have reported that the sphingomyelinase/ceramide system, which has shown several times to be a crucial factor in the internalization, processing, and killing of diverse pathogens, also modulates the pro-inflammatory response and the state of mycobacteria in macrophages, which highlights the important role of lipid metabolism in the pathogenic mechanism of pathogens. For PTB, sphingosine-1 phosphate (S1P) and ceramide are central molecules and are decisive for sphingolipid signaling; otherwise, they are about the secretion of interferon (IFN)-γ during the course of infection and infiltration of pulmonary CD11b+ macrophages and expression of S-1P receptor-3 (S-1PR3) in the lungs during the course of infection (Braverman et al., 2016; Nadella et al., 2019). Takenami et al. (2018) have a research about the IgM and total IgG antibody response to cardiolipin (CL), phosphatidylcholine (PTC), phosphatidylethanolamine (PE), phosphatidylinositol (PI), and sulfatide (SL-I) as biosignatures that can be used to diagnose PTB and its applicability for monitoring the efficacy of antituberculosis treatment; the antibody concentrations of PTB patients were significantly higher than those of healthy control, which also indicates that lipids play a vital role in the process of Mycobacterium tuberculosis infection. Also, fatty acids can stimulate the activation of dormant Mycobacterium tuberculosis in liquid medium (Shleeva et al., 2013). Eicosanoids, lipid mediators derived from arachidonic acid, have been associated with the modulation of the host response to Mycobacterium tuberculosis infection (Bafica et al., 2005; Tobin et al., 2012). Moreover, it also reported increased eicosanoid ratios in plasma in tuberculosis patients compared to healthy control (Mayer-Barber and Sher, 2015). However, related research about the lipid in osteoarticular tuberculosis or even other extrapulmonary tuberculosis is rare; the results of our research, whether from the differential metabolites screened out or the signal pathways obtained from bioinformatics, are closely related to the results of pulmonary tuberculosis-related research (Huang et al., 2020; Han et al., 2021). It provides more reference value for future research on diagnostic biomarkers from the perspective of metabolomics and provides some hints for the establishment of subsequent auxiliary diagnostic methods for these differential metabolites.
In the future study, we will expand the scope of clinical serum samples, from the single-center collection to the multicenter collection in different regions, in order to minimize the limitations of the experimental results. At the same time, serum samples of different tuberculosis diseases such as pulmonary tuberculosis can be added as follow-up controls to determine whether the relevant markers are widely related to tuberculosis, and then extended to not only the auxiliary diagnosis for osteoarticular tuberculosis but also the establishment of a new method of laboratory examination for tuberculosis for reference.
Data Availability Statement
The data presented in the study are deposited in the MetaboLights repository, accession number MTBLS4187 with URL www.ebi.ac.uk/metabolights/MTBLS4187.
Ethics Statement
The studies involving human participants were reviewed and approved by the Ethics Committee of Chinese PLA General Hospital. The patients/participants provided their written informed consent to participate in this study.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
This study did not receive any specific grants from funding agencies in the public, commercial, or not-for-profit sectors. We sincerely thank the staff of Beijing Omics Biotechnology Co., Ltd., for their help with LC-MS/MS analysis.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2022.827528/full#supplementary-material