A Transcriptomic Biomarker Predicting Linezolid-Associated Neuropathy During Treatment of Drug-Resistant Tuberculosis

Authors

Nika Zielinski^1,2,3, Dragos Baiceanu^4,5, Antonela Dragomir^4,5,6, Jan Heyckendorf⁷, Elmira Ibraim^4,5, Niklas Köhler^1,2,3, Christoph Leschczyk^2,8, Cristina Popa^4,5, Andrea Rachow^9,10,11, Jens Sachsenweger¹², Patricia Sanchez Carballo^1,2,3, Dagmar Schaub^1,2,3, Hajo Zeeb^13,14, Begna Tulu^1,2,3, Andrew R. DiNardo^15,16, Christoph Lange^1,2,3,15, Maja Reimann^1,2,3

Affiliated Institutions

¹Division of Clinical Infectious Diseases, Research Center Borstel, Borstel, Germany
²German Center for Infection Research (DZIF) Partner Site Hamburg-Lübeck-Borstel-Riems, Borstel, Germany
³Respiratory Medicine and International Health, University of Lübeck, Lübeck, Germany
⁴Marius Nasta Institute of Pneumophtiziology (MNI), Bucharest, Romania
⁵Eastern-European Study Site of DZIF in MNI, Bucharest, Romania
⁶UMF Carol Davila, Bucharest, Romania
⁷Clinic for Internal Medicine I, Leibniz Lung Clinic, University Hospital Schleswig-Holstein (UKSH) Campus Kiel, Kiel, Germany
⁸Division of Cellular Microbiology, Research Center Borstel, Borstel, Germany
⁹Division of Infectious Diseases and Tropical Medicine, Medical Centre of the University of Munich (LMU), Munich, Germany
¹⁰German Centre for Infection Research (DZIF), Partner Site Munich, Munich, Germany
¹¹Unit Global Health, Helmholtz Zentrum München, German Research Centre for Environmental Health (HMGU), Neuherberg, Germany
¹²Department of Pneumology, Asklepios Clinic Hamburg-Harburg, Hamburg, Germany
13Department of Prevention and Evaluation, Leibniz Institute for Prevention Research and Epidemiology – BIPS, Bremen, Germany
¹⁴Faculty of Human and Health Sciences, University of Bremen, Bremen, Germany
¹⁵Baylor College of Medicine and Texas Children’s Hospital, Global TB Program, Houston, Texas
¹⁶Radboud University Medical Center, Internal Medicine, Nijmegen, Netherlands

Corresponding Author

Dr. Maja Reimann
Email: [email protected]

SUGGESTED CITATION

Zielinski N, Baiceanu D, Dragomir A, Heyckendorf J, Ibraim E, Köhler N, Leschczyk C, Popa C, Rachow A, Sachsenweger J, Sanchez Carballo P, Schaub D, Zeeb H, Tulu B, DiNardo AR, Lange C, Reimann M. A Transcriptomic Biomarker Predicting Linezolid-Associated Neuropathy During Treatment of Drug-Resistant Tuberculosis. Pathogens and Immunity. 2024;9(2):25–42. doi: 10.20411/pai.v9i2.705

DOI

10.20411/pai.v9i2.705

Abstract

Background: Neuropathic adverse events occur frequently in linezolid-containing regimens, some of which remain irreversible after drug discontinuation.

Objective: We aimed to identify and validate a host RNA-based biomarker that can predict linezolid-associated neuropathy before multidrug-resistant/rifampicin-resistant tuberculosis (MDR/RR-TB) treatment initiation and to identify genes and pathways that are associated with linezolid-associated neuropathy.

Methods: Adult patients initiating MDR/RR-TB treatment including linezolid were prospectively enrolled in 3 independent cohorts in Germany. Clinical data and whole blood RNA for transcriptomic analysis were collected. The primary outcome was linezolid-associated optic and/or peripheral neuropathy. A random forest algorithm was used for biomarker identification. The biomarker was validated in an additional fourth cohort of patients with MDR/RR-TB from Romania.

Results: A total of 52 patients from the 3 identification cohorts received linezolid treatment. Of those, 24 (46.2%) developed peripheral and/or optic neuropathies during linezolid treatment. The majority (59.3%) of the episodes were of moderate (grade 2) severity. In total, the expression of 1,479 genes differed significantly at baseline of treatment. Suprabasin (SBSN) was identified as a potential biomarker to predict linezolid-associated neuropathy. In the validation cohort, 10 of 42 (23.8%) patients developed grade ≥3 neuropathies. The area under the curve for the biomarker algorithm prediction of grade ≥3 neuropathies was 0.63 (poor; 95% confidence interval: 0.42 – 0.84).

Conclusions: We identified and preliminarily validated a potential clinical biomarker to predict linezolid-associated neuropathies before the initiation of MDR/RR-TB therapy. Larger studies of the SBSN biomarker in more diverse populations are warranted.

Keywords

adverse events; linezolid; neurotoxicity; precision medicine; SBSN; tuberculosis; multidrug resistance

Introduction

According to the latest report by the World Health Organization (WHO), 10.6 million people developed tuberculosis (TB) in 2022. Second to COVID-19, TB has been the leading cause of death by an infectious disease worldwide in 2022 [1], and the emergence of drug-resistant TB is challenging control and successful treatment of this disease [2]. Based on the availability of novel drugs and treatment regimens, the treatment of drug-resistant TB has recently been revised substantially [3]. The one-armed Nix-TB Trial showed that a 3-drug regimen including bedaquiline, pretomanid, and linezolid (BPaL) administered over 6 months resulted in 90% treatment success in patients who had either failed multidrug-resistant tuberculosis (MDR-TB) treatment or had extensive drug-resistant tuberculosis (XDR-TB) [4]. However, a high rate of linezolid-associated adverse events was observed on a 1,200 mg dose administered once daily over 6 months, with 80.7% of patients experiencing polyneuropathy and 47.7% of patients experiencing bone marrow toxicity with thrombocytopenia and/or anemia.

Results from an individual patient database including 9,178 patients with multidrug-resistant/rifampicin-resistant tuberculosis (MDR/RR-TB) indicate that linezolid is the most toxic of the second-line anti-TB medicines and that 14.1% (95% confidence interval [95%CI]: 9.9% – 19.6%) of patients receiving linezolid had to discontinue the drug due to the occurrence of adverse events [5]. Neuropathies associated with linezolid can remain irreversible if the drug is not discontinued in time, and no effective therapy has been established yet for patients who are affected [6]. With the aim to maintain the high efficacy of the BPaL regimen but to reduce the linezolid-attributed toxicity, investigators of the ZeNix-Trial evaluated standard 6-month dosages of bedaquiline and pretomanid with variable dosages and durations of linezolid over 2 and 6 months with 600 mg and 1,200 mg in 4 parallel study arms.

The best benefit-risk profile of linezolid was observed in the arm with 600 mg linezolid over 6 months. Here, 24.4% of patients still experienced neuropathies compared to 37.8% patients in the arm with 1,200 mg linezolid over 6 months [7]. In the TB-PRACTECAL study, standard dosages of bedaquiline, pretomanid, and moxifloxacin were administered with 600 mg linezolid over 16 weeks followed by 300 mg linezolid over 8 weeks with 9.3% of patients developing neuropathy [8]. As a result of these trials, the WHO has issued novel management guidelines for the treatment of drug-resistant TB in 2022 prioritizing the 6-month BPaLM treatment regimen with a linezolid dosage of 600 mg daily [3].

So far, the mechanism of action for linezolid-associated neuropathies is yet to be fully understood. Linezolid was found to inhibit autophagy flux [9] and mitochondrial protein synthesis [10], which could have an effect on the peripheral nerves by disrupting their normal function. It is important to note that a dose reduction of linezolid resulted in a decrease of attributed adverse events [7]. However, adverse events still occurred in parts of the population even at a low dose of linezolid, and it was also found that linezolid trough concentrations correlate with mitochondrial toxicity-related adverse events [11]. Adverse events, even when mild, can lead to decreased therapy adherence and consequently decreased cure rates [12]. Given this situation, it would be highly desirable to be able to predict which patients are likely to develop neuropathies when treated with linezolid. Presently, there are no biomarkers available that allow a prediction of linezolid-related adverse events in patients initiating treatment for drug-resistant TB.

Transcriptomic signatures are promising biomarkers for the prognosis of the development of TB [13, 14], diagnosis of active TB [15], and for treatment monitoring [16]. We aimed to identify and validate a biomarker model that predicts the occurrence of peripheral and/or optic neuropathy in patients receiving linezolid as part of an MDR/RR-TB treatment regimen before treatment initiation.

Methods

Study Design and Participants

This study comprises a subgroup of patients recruited as part of an observational cohort study to evaluate various endpoints and clinical questions before, during, and after treatment of sensitive and MDR-TB with regard to transcriptome profiles for potential biomarker identification. Adult patients with culture-confirmed pulmonary MDR/RR-TB and (pre-) XDR-TB without HIV infection from 4 multicenter prospective cohort studies who received 600 mg linezolid once daily as part of a treatment regimen were included in this study (Figure 1). Three of those cohorts were enrolled in Germany: The German Identification Cohort (GIC; enrollment period: March 01, 2013 to September 30, 2016 in 5 clinical centers), the German Validation Cohort (GVC; enrollment period: March 01, 2015 to April 30, 2018 in 7 clinical centers), and the Second German Validation Cohort (SGVC; enrollment period: October 01, 2018 to July 31, 2021 at the Research Center Borstel). The fourth cohort (Eastern European Study Side Cohort [EESSC]) was enrolled between May 01, 2019 and March 31, 2023 at the Marius-Nasta Institute in Bucharest.

Sputum samples were genotypically tested for rifampicin resistance using Xpert MTB/RIF (Cepheid, Sunnyvale) and tested for second-line drug susceptibility using GenoType MTB-DRplus and MTBDRsl (Hain Life Sciences, Nehren). Phenotypic resistance testing of bacteria was additionally performed. The treatment regimen for MDR/RR-TB was individualized based on the results of genotypic and phenotypic testing of the pathogen and the current treatment guidelines [17–22]. For the entire inpatient stay as well as the follow-up period, information on the therapy regimen including antibiotics used, therapy duration (in days), therapy interruptions, and the occurrence of adverse events was captured. Clinical data of the patients, such as age, sex, time to culture positivity (TTP; MGIT liquid culture), and time to culture conversion (TCC), were also collected. At the start of treatment, all patients were tested for HIV infection.

The onset of linezolid-associated neuropathy was dated to the time when the first report of symptoms was noted in the patient’s medical record or when linezolid was discontinued in response to clinical signs of neurotoxicity. Optic neuropathy was evaluated by an ophthalmologist if symptoms were stated, using fundoscopy visual accuracy testing and visual field testing. There was no standardized testing for peripheral neuropathies. Symptoms that have led to treatment discontinuation were stocking and/or glove type of sensory loss to touch and/or temperature misperception for peripheral neuropathies, and visual impairment for optic neuropathies. Severity grading was retrospectively performed by physicians using the endTB severity grading scale, version 5 [23]. Peripheral and optic neuropathies were considered together as a binary outcome variable in most analyses.

Whole blood transcriptomic RNA was collected with PAXgene tubes (Qiagen®, Venlo). In the identification cohort, PAXgene tubes were collected before the start of treatment, after 14 days of treatment, at the time of smear conversion, at the time of culture conversion, and after 6, 10, 15, and 20 months of treatment. If the treatment was longer than 20 months, another PAXgene tube was collected at the actual end of treatment. In the validation cohort, PAXgene tubes were collected before the start of treatment, after 13 days, after 8 weeks, and after 4, 6, 10, 15, and 20 months of treatment. In this study, transcriptomic RNA before the start of treatment was used only. The processing of the microarray data for the transcriptomic analyses has been described elsewhere in detail [16].

Statistical Analysis

Patient characteristics were analyzed using absolute and relative frequencies as well as mean with standard deviation or median with interquartile range (IQR) depending on distribution. Differences between the group of patients with and without linezolid-associated neuropathies were identified through t-test, Wilcoxon rank sum test, or Fisher’s exact test. The significance level was set at a P-value <0.05 for all analyses. The cumulative probability for linezolid-associated neuropathies during the course of treatment was calculated using a Kaplan-Meier estimator.

For the analysis of differentially expressed genes and pathways, over-representation analysis was performed by using genes with P-value <0.05 that were analyzed using the molecular signature of database for Hallmark, Reactome, KEGG, and Wikipathways. Ingenuity pathway analysis was used to calculate pathway z scores to clarify directionality.

Data from GIC, GVC, and SGVC were combined as the identification cohort for the development of the biomarker algorithm. The preselection of genes out of 44,000 z score standardized transcripts was done by using a LASSO regression with a binary classification for linezolid-associated neuropathy as the dependent variable. After setting up the first random forest algorithm with the randomForest package including the preselected genes from the LASSO regression, features with a mean decrease of the Gini coefficient <3 were removed. The random forest algorithm was externally validated in the EESSC cohort with the receiver operator characteristic (ROC) curve and corresponding area under the curve (AUC) including the 95%CI. Severe linezolid-associated neuropathies with a severity grade ≥3 were considered for the validation of the biomarker algorithm.

All statistical analyses were performed with R software (version 4.3.0) and GraphPad Prism (version 10.2.1).

Ethics

The study was initially approved by the Ethics Committee of the University of Lübeck, Germany (AZ 12-233) and then approved by the local Ethic Committees of all participating centers in Germany. Study approval for the Romanian cohort was granted by the Marius Nasta Institute (3181/25.03.2015; Bucharest, Romania).

Results

In total, 52 patients were eligible for inclusion in the identification cohort (Figure 1). The median age of patients in the identification cohort was 35.5 years (IQR 26.0 – 42.0 years); 75.0% were male, and half of the patients had MDR/RR-TB (51.9%). The median TTP was 15.0 days (IQR 11.0 – 20.0 days) at baseline, the median TCC was 53.5 days (IQR 24.0 – 77.0 days), and the median linezolid therapy duration was 242.0 days (IQR 145.0 – 532.5 days). Almost three- quarters (73.1%) of the patients completed treatment with a definition of cure according to the Tuberculosis Network European Trials Group (TBNET). Six patients could not be followed up to the end of study. Patients with and without linezolid-associated neuropathies did not differ significantly with regard to socio-demographic characteristics, comorbidities, and markers for disease severity (Table 1). In the identification cohort, 24 (46.2%) patients developed linezolid-associated neuropathies. Of these, 20 patients had sensory peripheral neuropathies only, 1 patient had optic neuropathy only, and 3 patients had both sensory and optic neuropathy. Accordingly, a total of 27 events of linezolid-associated neuropathies occurred, of which 10 (37.0%) events could be assigned to severity grade 1, 16 (59.3%) events to severity grade 2, 1 (3.7%) event to severity grade 3, and 0 to severity grade 4. In all cases, therapy with linezolid was discontinued at the onset of symptoms. Linezolid-associated neuropathies occurred after a median of 198.5 days (IQR 173.0 – 330.5 days). The cumulative probability increased continuously over the course of therapy and reached over 50.0% after 381 days of linezolid treatment duration in the identification cohort (Figure 2). The validation cohort consisted of 42 patients (Figure 1), the majority of whom were male (29 patients, 69.0%), similar to the identification cohort (P-value = 0.644). While there were some age differences between the cohorts, with a higher median age in the validation than in the identification cohort (49.0 years [IQR 36.0 – 55.0] vs 35.5 years [IQR 26.0 – 42.0]; P-value <0.001), none of these characteristics differed significantly between the group of patients with and without linezolid-associated neuropathy within the validation cohort (Table 1). In the validation cohort, 17 patients (40.5%) developed linezolid-associated neuropathies, 12 patients had sensory peripheral neuropathy only, and 5 patients had a combined form of both sensory peripheral and optic neuropathy. Severity was classified on an individual patient level and not on event level in the validation cohort. Here, 3 patients were assigned to severity grade 1, 4 patients to severity grade 2, 6 patients to severity grade 3, and 4 patients to severity grade 4. Linezolid was discontinued in 6 patients.

In the identification cohort, 1,479 genes were statistically significant differentially expressed before treatment initiation when comparing the group of patients with and without linezolid-associated neuropathy. Two clusters were found with regard to the statistically significant differentially expressed genes that correlate with occurrence of linezolid-associated neuropathy (Figure 3).

Patients with linezolid-associated neuropathies showed a higher expression of genes associated with cellular proliferation, autophagy, Wnt signaling, p53, and immune activation. In contrast, patients without linezolid-associated neuropathy had an increase in genes associated with NRF2-mediated oxidative stress response, the ATP-dependent DNA damage response, tyrosine kinase receptor signaling, MAPK signaling, PI3K-AKT signaling, and the Il-1 pathway (Figure 4A, Supplementary Table 1). Specific genes increased in patients without linezolid-associated neuropathy included ATG4C, BAX, and IGF1, while patients on linezolid who suffered from neuropathy had an increase in AURKB, CASP3, CASP8, CDK1, CDKN1B, CEBPB, NADSYN1, NRF1, GPX7, and SBSN (Figure 4B, Figure 4C, Supplementary Table 2).

The final biomarker algorithm consists of one gene (SBSN). SBSN expression was significantly lower in patients with linezolid-associated neuropathy compared to patients without linezolid-associated neuropathy (median -0.53 vs. 0.38; P-value <0.001, Figure 4C). The biomarker random forest algorithm was trained in the entire identification cohort.

The discrimination for predicting severe linezolid-associated neuropathies using the biomarker algorithm in the validation cohort showed an AUC of 0.628 (95%CI 0.415 – 0.841) as displayed in Figure 5. This is considered to be poor accuracy.

Out of the 42 patients starting linezolid treatment, 32 did not develop linezolid-associated neuropathies. The biomarker algorithm correctly predicted this outcome in 21 patients corresponding to a specificity of 65.6% (95%CI 46.8% – 81.4%). Ten patients developed severe linezolid-associated neuropathies, and the biomarker algorithm predicted this correctly in 6 patients, corresponding to a sensitivity of 60.0% (95%CI 26.2% – 87.8%). Changes in sensitivity and specificity with different cut-off values as well as positive and negative predictive values for the biomarker algorithm can be found in Supplementary Table 3.

Table 1. Patient Characteristics in the Identification and Validation Cohort

(a) Three missing values. (b) Eight missing values. (c) One missing value.
SD = Standard deviation; MDR/RR-TB = Multidrug-resistant/rifampicin-resistant tuberculosis; XDR-TB = Extensive drug-resistant tuberculosis; TTP = Time to culture positivity; IQR = Interquartile range; TCC = Time to culture conversion; TBNET = Tuberculosis Network European Trials Group.

Figure 1. Flow chart for study inclusion and reason for exclusion in identification and validation cohort. GIC = German Identification Cohort; GVC = German Validation Cohort, SCVC = Second German Validation Cohort; EESSC = Eastern European Study Side Cohort; MDR/RR-TB = Multidrug-resistant/rifampicin-resistant tuberculosis; XDR-TB = Extensive drug-resistant tuberculosis

Figure 2. Kaplan-Meier curve with 95% confidence interval for the cumulative probability of linezolid-associated neuropathies during the course of linezolid treatment. The entire identification cohort (n=52) from 0 – 600 days is considered. Censored patients are shown with a cross in the curve.

Figure 3. Heat map of 1,479 statistically significant differentially expressed genes prior to treatment initiation. The heat map shows the genes on patient-level and z score normalized. The in rows and columns are hierarchically clustered. Upregulated genes are colored red and downregulated genes are colored blue. The annotation highlights patients with (orange; AE = linezolid-associated neuropathy) and patients without linezolid-associated neuropathy (green; NOAE = no linezolid-associated neuropathy) across the top of the heat map.

Figure 4. Selected statistically significant differentially expressed pathways and genes in the group of patients with linezolid-associated neuropathy (AE) or without linezolid-associated neuropathy (NOAE). Figure 4A. Over-representation analysis using the molecular signature of database for Hallmark, Reactome, KEGG, and Wikipathways. Figure 4B. Violin plots of z score standardized gene expression levels grouped by biological function. Figure 4C. Violin plot of the SBSN z score standardized gene expression levels.

Figure 5. Receiver-Operator Characteristic (ROC) for evaluation of the prediction power of the biomarker algorithm in the validation population. Shown is the ROC curve with the area under the curve (AUC) for predicting linezolid-associated neuropathies with grade 3 or higher in the validation cohort (AUC = 0.628 [95%CI: 0.415 – 0.841]).

Discussion

In 4 prospective cohorts of patients undergoing treatment for drug-resistant TB with a linezolid-containing regimen, we identified and validated a transcriptomic biomarker model to predict the occurrence of neuropathies before treatment initiation. The biomarker model with SBSN predicts the occurrence of neuropathies with an accuracy of AUC = 0.628. This is, to our knowledge, the first biomarker that has the potential to predict adverse events prior to initiation of second-line anti-TB therapy, although the prediction to discriminate patients who develop neuropathies from those who do not following the initiation of a linezolid-containing regimen was still poor with this biomarker.

SBSN is a gene that is located in the q13 region of chromosome 19 and codes for the protein suprabasin. Suprabasin was first described for its function in epidermal differentiation [24]. Studies have linked suprabasin expression to human diseases such as malignant tumors [25, 26] and autoimmune diseases [27, 28]. Even though its exact function in the development of neuropathies is not established, recent findings have shown that suprabasin expression in the nervous system could play a role in neuronal function, development, and regeneration, as well as axonal growth [29, 30]. Axonal growth, as well as neuronal survival and plasticity, may be influenced through its effect on signaling pathways (AKT, WNT/β-catenin, and/or p38MAPK). These pathways are known to have an impact on several cellular processes including migration, growth, apoptosis, and immune resistance [31]. Thus, dysfunction or dysregulation of suprabasin in the nervous system could potentially contribute to the development or progression of neuropathic conditions. Suprabasin is secreted in exosomes that regulate NF-κB via NEMO [32]. Should suprabasin be activating NF-κB to induce neurotoxicity, this mechanism could possibly provide a therapeutic biomarker or target for prevention. This finding should be validated in future studies.

Comparable studies are lacking so far. In one study, expression levels of 2 genes (MKI67 and SLC22A8) were associated with neuropathies in rats receiving linezolid [33]. MKI67, a gene involved in cellular proliferation, was not increased in our study. However, our study similarly found that patients suffering from linezolid-associated neuropathy had an increase in genes related to cellular proliferation. SLC22A8 was also not increased in our study. A possible explanation could be the limited transferability from animal models to the human organism and the use of pre-treatment gene expression data in our study, whereby linezolid-based changes in gene expression after linezolid exposure were not represented. The identification of genes for the biomarker algorithm followed an exploratory, statistical approach, which offered the advantage that previously unknown genes could be identified in this context. To date, there are no studies that have emphasized SBSN in the context of drug-associated neuropathies.

The biomarker algorithm was developed in a multicenter German cohort with closely monitored clinical data and microarray transcriptomic data from standardized blood collection and laboratory preparation. The incidence of linezolid-associated neuropathies overall and the time to sensory peripheral neuropathies is comparable to other recently published cohort studies that used an individualized treatment regimen [6, 34, 35]. However, the incidence of linezolid-associated neuropathies was lower in the ZeNix-Trial (24.4%) and TB-PRACTECAL (9.3%) where the treatment duration of linezolid was shorter in the BPaL(M) regimen [7, 8]. The prediction power of the biomarker algorithm was evaluated in an additional fourth, non-German external cohort with an AUC of 0.628 (95%CI 0.415 – 0.841). If the biomarker algorithm had been applied before treatment initiation, 65.6% of linezolid-associated polyneuropathies could have been prevented.

However, in 4 out of 10 cases, the biomarker algorithm did not make the correct prediction. Based on our results, using the algorithm in clinical practice would have meant that linezolid be withheld in 34.4% of patients who did not develop a neuropathy. Close monitoring of symptoms in affected patients could therefore also be considered as an alternative recommendation. One aspect that is not considered in this biomarker algorithm but was found to be associated with the development of linezolid-associated neuropathies, are linezolid trough levels. While one study did not find statistically significant differences between trough levels of >2 mg/L and ≤2 mg/L and the incidence of linezolid-associated neuropathies [34], others found in their study that in patients with MDR-TB, the incidence of peripheral neuropathies was significantly higher in patients with a maximum trough level of >2 mg/L compared to ≤2 mg/L (46.2% vs 26.6%, P-value = 0.02) [6]. The combination of the biomarker algorithm with therapeutic drug monitoring could possibly fine-tune the prediction of neuropathies.

A limitation of this study is that the identified biomarker algorithm only focuses on the prediction of linezolid-associated neuropathies but does not consider other severe and potential life-threatening adverse drug reactions such as myelosuppression or lactic acidosis, which can play an important role in clinical practice. The biomarker algorithm has only been trained and preliminary validated in adult patients with pulmonary TB treated with 600 mg of linezolid once daily and not in people living with HIV. Further validation should be done in larger studies including, for example, children, elderly patients, and people living with HIV, in the BPaL(M) regimen as well as in patients from different areas of the world. In addition, an interplay of individual predisposing factors, consisting of physical and genetic characteristics, comorbidities, and lifestyle and environmental influences, can already lead to nerve damage independent of the intake of linezolid, which can occur spontaneously at the same time as taking the drug or only be unmasked by taking the drug [36]. This could limit the prediction power of the biomarker algorithm. The sample size in the identification and validation cohort is rather small. However, we would assume that the risk of overfitting is small because we used a random forest algorithm that has been shown to be suitable for dealing with OMICS data and offers the advantage that it is relatively stable against overfitting due to ensemble learning [37]. As the prediction accuracy of the SBSN biomarker remains poor at present, future studies should explore it further, including the measurement of levels from circulating blood to verify these findings.

Conclusion

In conclusion, we identified and preliminarily validated a potential clinical biomarker to predict linezolid-associated neuropathies before the initiation of MDR/RR-TB therapy. Linezolid-associated neuropathies can be therapy limiting and considerably restrict the quality of life of affected patients. Early detection of linezolid-associated neuropathies can bring a decisive advantage for the prognosis of the neuropathy. Biomarker-based prediction of adverse events opens a novel avenue to improve the management of patients at risk for adverse events in the future.

Acknowledgments

We thank Jessica Hofmeister, Franziska Daduna, Sandra Nyenhuis, Frauke Koops, and Lasse Möller for technical assistance and the DZIF study group for data contribution.

Funding

This work was supported by the German Center of Infection Research under grant agreement TTU-TB 02.709.

Potential Conflicts of Interest

The authors have no conflicts of interest to declare.

Supplementary Data

Supplementary materials are available at the Pathogens and Immunity website. Supplementary data may be provided by the authors to benefit the reader. Supplementary data are not copyedited and are the sole responsibility of the authors. Questions or comments related to supplementary materials should be addressed to the corresponding author.

Supplementary Tables

References

1. World Health Organization. Global tuberculosis report 2023. Geneva: World Health Organization; 2023.

2. Lange C, Dheda K, Chesov D, Mandalakas AM, Udwadia Z, Horsburgh CR, Jr. Management of drug-resistant tuberculosis. Lancet. 2019;394(10202):953-66. doi: 10.1016/S0140-6736(19)31882-3. PubMed PMID: 31526739.

3. World Health Organization. WHO consolidated guidelines on tuberculosis: module 4: treatment: drug-resistant tuberculosis treatment. 2022 update ed. Geneva: World Health Organization; 2022.

4. Conradie F, Diacon AH, Ngubane N, Howell P, Everitt D, Crook AM, Mendel CM, Egizi E, Moreira J, Timm J, McHugh TD, Wills GH, Bateson A, Hunt R, Van Niekerk C, Li M, Olugbosi M, Spigelman M, Nix TBTT. Treatment of Highly Drug-Resistant Pulmonary Tuberculosis. N Engl J Med. 2020;382(10):893-902. doi: 10.1056/NEJMoa1901814. PubMed PMID: 32130813; PMCID: PMC6955640.

5. Lan Z, Ahmad N, Baghaei P, Barkane L, Benedetti A, Brode SK, Brust JCM, Campbell JR, Chang VWL, Falzon D, Guglielmetti L, Isaakidis P, Kempker RR, Kipiani M, Kuksa L, Lange C, Laniado-Laborin R, Nahid P, Rodrigues D, Singla R, Udwadia ZF, Menzies D, Collaborative Group for the Meta-Analysis of Individual Patient Data in MDRTBt. Drug-associated adverse events in the treatment of multidrug-resistant tuberculosis: an individual patient data meta-analysis. Lancet Respir Med. 2020;8(4):383-94. doi: 10.1016/S2213-2600(20)30047-3. PubMed PMID: 32192585; PMCID: PMC7384398.

6. Eimer J, Frechet-Jachym M, Le Du D, Caumes E, El-Helali N, Marigot-Outtandy D, Mechai F, Peytavin G, Pourcher V, Rioux C, Yazdanpanah Y, Robert J, Guglielmetti L, group L. Association Between Increased Linezolid Plasma Concentrations and the Development of Severe Toxicity in Multidrug-Resistant Tuberculosis Treatment. Clin Infect Dis. 2023;76(3):e947-e56. doi: 10.1093/cid/ciac485. PubMed PMID: 35717636.

7. Conradie F, Bagdasaryan TR, Borisov S, Howell P, Mikiashvili L, Ngubane N, Samoilova A, Skornykova S, Tudor E, Variava E, Yablonskiy P, Everitt D, Wills GH, Sun E, Olugbosi M, Egizi E, Li M, Holsta A, Timm J, Bateson A, Crook AM, Fabiane SM, Hunt R, McHugh TD, Tweed CD, Foraida S, Mendel CM, Spigelman M, ZeNix Trial T. Bedaquiline-Pretomanid-Linezolid Regimens for Drug-Resistant Tuberculosis. N Engl J Med. 2022;387(9):810-23. doi: 10.1056/NEJMoa2119430. PubMed PMID: 36053506; PMCID: PMC9490302.

8. Nyang’wa BT, Berry C, Kazounis E, Motta I, Parpieva N, Tigay Z, Solodovnikova V, Liverko I, Moodliar R, Dodd M, Ngubane N, Rassool M, McHugh TD, Spigelman M, Moore DAJ, Ritmeijer K, du Cros P, Fielding K, Collaborators T-PS. A 24-Week, All-Oral Regimen for Rifampin-Resistant Tuberculosis. N Engl J Med. 2022;387(25):2331-43. doi: 10.1056/NEJMoa2117166. PubMed PMID: 36546625.

9. Yuan Y, Li J, Chen Y, Cai Q, Xu Y, Lin L, Lang Y, Guo S, Zhang R, Cai X. Mechanism underlying linezolid-induced peripheral neuropathy in multidrug-resistant tuberculosis. Front Pharmacol. 2022;13:946058. doi: 10.3389/fphar.2022.946058. PubMed PMID: 36160387; PMCID: PMC9500448.

10. De Vriese AS, Coster RV, Smet J, Seneca S, Lovering A, Van Haute LL, Vanopdenbosch LJ, Martin JJ, Groote CC, Vandecasteele S, Boelaert JR. Linezolid-induced inhibition of mitochondrial protein synthesis. Clin Infect Dis. 2006;42(8):1111-7. doi: 10.1086/501356. PubMed PMID: 16575728.

11. Song T, Lee M, Jeon HS, Park Y, Dodd LE, Dartois V, Follman D, Wang J, Cai Y, Goldfeder LC, Olivier KN, Xie Y, Via LE, Cho SN, Barry CE, 3rd, Chen RY. Linezolid Trough Concentrations Correlate with Mitochondrial Toxicity-Related Adverse Events in the Treatment of Chronic Extensively Drug-Resistant Tuberculosis. EBioMedicine. 2015;2(11):1627-33. doi: 10.1016/j.ebiom.2015.09.051. PubMed PMID: 26870788; PMCID: PMC4740314.

12. Lange C, Aarnoutse RE, Alffenaar JWC, Bothamley G, Brinkmann F, Costa J, Chesov D, van Crevel R, Dedicoat M, Dominguez J, Duarte R, Grobbel HP, Gunther G, Guglielmetti L, Heyckendorf J, Kay AW, Kirakosyan O, Kirk O, Koczulla RA, Kudriashov GG, Kuksa L, van Leth F, Magis-Escurra C, Mandalakas AM, Molina-Moya B, Peloquin CA, Reimann M, Rumetshofer R, Schaaf HS, Schon T, Tiberi S, Valda J, Yablonskii PK, Dheda K. Management of patients with multidrug-resistant tuberculosis. Int J Tuberc Lung Dis. 2019;23(6):645-62. doi: 10.5588/ijtld.18.0622. PubMed PMID: 31315696.

13. Zak DE, Penn-Nicholson A, Scriba TJ, Thompson E, Suliman S, Amon LM, Mahomed H, Erasmus M, Whatney W, Hussey GD, Abrahams D, Kafaar F, Hawkridge T, Verver S, Hughes EJ, Ota M, Sutherland J, Howe R, Dockrell HM, Boom WH, Thiel B, Ottenhoff THM, Mayanja-Kizza H, Crampin AC, Downing K, Hatherill M, Valvo J, Shankar S, Parida SK, Kaufmann SHE, Walzl G, Aderem A, Hanekom WA, Acs, groups GCcs. A blood RNA signature for tuberculosis disease risk: a prospective cohort study. Lancet. 2016;387(10035):2312-22. doi: 10.1016/S0140-6736(15)01316-1. PubMed PMID: 27017310; PMCID: PMC5392204.

14. Scriba TJ, Fiore-Gartland A, Penn-Nicholson A, Mulenga H, Kimbung Mbandi S, Borate B, Mendelsohn SC, Hadley K, Hikuam C, Kaskar M, Musvosvi M, Bilek N, Self S, Sumner T, White RG, Erasmus M, Jaxa L, Raphela R, Innes C, Brumskine W, Hiemstra A, Malherbe ST, Hassan-Moosa R, Tameris M, Walzl G, Naidoo K, Churchyard G, Hatherill M, Team C-S. Biomarker-guided tuberculosis preventive therapy (CORTIS): a randomised controlled trial. Lancet Infect Dis. 2021;21(3):354-65. doi: 10.1016/S1473-3099(20)30914-2. PubMed PMID: 33508224; PMCID: PMC7907670.

15. Heyckendorf J, Reimann M, Marwitz S, Lange C, group D-Tcs. Pathogen-free diagnosis of tuberculosis. Lancet Infect Dis. 2021;21(8):1066. doi: 10.1016/S1473-3099(21)00337-6. PubMed PMID: 34331877.

16. Heyckendorf J, Marwitz S, Reimann M, Avsar K, DiNardo AR, Gunther G, Hoelscher M, Ibraim E, Kalsdorf B, Kaufmann SHE, Kontsevaya I, van Leth F, Mandalakas AM, Maurer FP, Muller M, Nitschkowski D, Olaru ID, Popa C, Rachow A, Rolling T, Rybniker J, Salzer HJF, Sanchez-Carballo P, Schuhmann M, Schaub D, Spinu V, Suarez I, Terhalle E, Unnewehr M, Weiner J, 3rd, Goldmann T, Lange C. Prediction of anti-tuberculosis treatment duration based on a 22-gene transcriptomic model. Eur Respir J. 2021;58(3). doi: 10.1183/13993003.03492-2020. PubMed PMID: 33574078.

17. Schaberg T, Bauer T, Castell S, Dalhoff K, Detjen A, Diel R, Greinert U, Hauer B, Lange C, Magdorf K, Loddenkemper R. [Recommendations for therapy, chemoprevention and chemoprophylaxis of tuberculosis in adults and children. German Central Committee against Tuberculosis (DZK), German Respiratory Society (DGP)]. Pneumologie. 2012;66(3):133-71. doi: 10.1055/s-0031-1291619. PubMed PMID: 22328186.

18. Schaberg T, Bauer T, Brinkmann F, Diel R, Feiterna-Sperling C, Haas W, Hartmann P, Hauer B, Heyckendorf J, Lange C, Nienhaus A, Otto-Knapp R, Priwitzer M, Richter E, Rumetshofer R, Schenkel K, Schoch OD, Schonfeld N, Stahlmann R. [Tuberculosis Guideline for Adults - Guideline for Diagnosis and Treatment of Tuberculosis including LTBI Testing and Treatment of the German Central Committee (DZK) and the German Respiratory Society (DGP)]. Pneumologie. 2017;71(6):325-97. doi: 10.1055/s-0043-105954. PubMed PMID: 28651293.

19. World Health Organization. Companion handbook to the WHO guidelines for the programmatic management of drug-resistant tuberculosis. Geneva: World Health Organization; 2014.

20. World Health Organization. WHO treatment guidelines for drug-resistant tuberculosis. 2016 update ed. Geneva: World Health Organization; 2016.

21. World Health Organization. WHO consolidated guidelines on drug-resistant tuberculosis treatment. Geneva: World Health Organization; 2019.

22. World Health Organization. WHO operational handbook on tuberculosis: module 4: treatment: drug-resistant tuberculosis treatment. Geneva: World Health Organization; 2020.

23. endTB Consortium. TB Severity Grading Scale Version 5. 2016.

24. Park GT, Lim SE, Jang SI, Morasso MI. Suprabasin, a novel epidermal differentiation marker and potential cornified envelope precursor. J Biol Chem. 2002;277(47):45195-202. doi: 10.1074/jbc.M205380200. PubMed PMID: 12228223; PMCID: PMC1283087.

25. Zhu J, Wu G, Li Q, Gong H, Song J, Cao L, Wu S, Song L, Jiang L. Overexpression of Suprabasin is Associated with Proliferation and Tumorigenicity of Esophageal Squamous Cell Carcinoma. Sci Rep. 2016;6:21549. doi: 10.1038/srep21549. PubMed PMID: 26899563; PMCID: PMC4761926.

26. Shao C, Tan M, Bishop JA, Liu J, Bai W, Gaykalova DA, Ogawa T, Vikani AR, Agrawal Y, Li RJ, Kim MS, Westra WH, Sidransky D, Califano JA, Ha PK. Suprabasin is hypomethylated and associated with metastasis in salivary adenoid cystic carcinoma. PLoS One. 2012;7(11):e48582. doi: 10.1371/journal.pone.0048582. PubMed PMID: 23144906; PMCID: PMC3492451.

27. Ichinose K, Ohyama K, Furukawa K, Higuchi O, Mukaino A, Satoh K, Nakane S, Shimizu T, Umeda M, Fukui S, Nishino A, Nakajima H, Koga T, Kawashiri SY, Iwamoto N, Tamai M, Nakamura H, Origuchi T, Yoshida M, Kuroda N, Kawakami A. Novel anti-suprabasin antibodies may contribute to the pathogenesis of neuropsychiatric systemic lupus erythematosus. Clin Immunol. 2018;193:123-30. doi: 10.1016/j.clim.2017.11.006. PubMed PMID: 29162406.

28. Aoshima M, Phadungsaksawasdi P, Nakazawa S, Iwasaki M, Sakabe JI, Umayahara T, Yatagai T, Ikeya S, Shimauchi T, Tokura Y. Decreased expression of suprabasin induces aberrant differentiation and apoptosis of epidermal keratinocytes: Possible role for atopic dermatitis. J Dermatol Sci. 2019;95(3):107-12. doi: 10.1016/j.jdermsci.2019.07.009. PubMed PMID: 31399284.

29. Sjostedt E, Zhong W, Fagerberg L, Karlsson M, Mitsios N, Adori C, Oksvold P, Edfors F, Limiszewska A, Hikmet F, Huang J, Du Y, Lin L, Dong Z, Yang L, Liu X, Jiang H, Xu X, Wang J, Yang H, Bolund L, Mardinoglu A, Zhang C, von Feilitzen K, Lindskog C, Ponten F, Luo Y, Hokfelt T, Uhlen M, Mulder J. An atlas of the protein-coding genes in the human, pig, and mouse brain. Science. 2020;367(6482). doi: 10.1126/science.aay5947. PubMed PMID: 32139519.

30. Kim SW, Kim KT. Expression of Genes Involved in Axon Guidance: How Much Have We Learned? Int J Mol Sci. 2020;21(10). doi: 10.3390/ijms21103566. PubMed PMID: 32443632; PMCID: PMC7278939.

31. Pribyl M, Hodny Z, Kubikova I. Suprabasin-A Review. Genes (Basel). 2021;12(1). doi: 10.3390/genes12010108. PubMed PMID: 33477529; PMCID: PMC7831088.

32. Chen H, Chen X, Zhang Z, Bao W, Gao Z, Li D, Xie X, Zhou P, Yang C, Zhou Z, Pan J, Kuang X, Tang R, Feng Z, Zhou L, Zhu D, Yang J, Wang L, Huang H, Tang D, Liu J, Jiang L. Extracellular vesicles-transferred SBSN drives glioma aggressiveness by activating NF-kappaB via ANXA1-dependent ubiquitination of NEMO. Oncogene. 2022;41(49):5253-65. doi: 10.1038/s41388-022-02520-6. PubMed PMID: 36316443.

33. Uehara T, Kondo C, Yamate J, Torii M, Maruyama T. A toxicogenomic approach for identifying biomarkers for myelosuppressive anemia in rats. Toxicology. 2011;282(3):139-45. doi: 10.1016/j.tox.2011.01.027. PubMed PMID: 21296123.

34. Jaspard M, Butel N, El Helali N, Marigot-Outtandy D, Guillot H, Peytavin G, Veziris N, Bodaghi B, Flandre P, Petitjean G, Caumes E, Pourcher V. Linezolid-Associated Neurologic Adverse Events in Patients with Multidrug-Resistant Tuberculosis, France. Emerg Infect Dis. 2020;26(8):1792-800. doi: 10.3201/eid2608.191499. PubMed PMID: 32687026; PMCID: PMC7392460.

35. Zhang P, Li W, Liu M, Zhan S, Zhang H, Deng G, Chen X. Linezolid-Associated Neuropathy in Patients with MDR/XDR Tuberculosis in Shenzhen, China. Infect Drug Resist. 2022;15:2617-24. doi: 10.2147/IDR.S365371. PubMed PMID: 35634579; PMCID: PMC9139335.

36. Head KA. Peripheral neuropathy: pathogenic mechanisms and alternative therapies. Altern Med Rev. 2006;11(4):294-329. PubMed PMID: 17176168.

37. Couronne R, Probst P, Boulesteix AL. Random forest versus logistic regression: a large-scale benchmark experiment. BMC Bioinformatics. 2018;19(1):270. doi: 10.1186/s12859-018-2264-5. PubMed PMID: 30016950; PMCID: PMC6050737.

Footnotes

Submitted April 1, 2024 | Accepted May 28, 2024 | Published June 25, 2024

Copyright

Copyright © 2024 The Authors. This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International License.