- Open Access
Proteomic signatures of metronidazole-resistant Trichomonas vaginalis reveal novel proteins associated with drug resistance
Parasites & Vectors volume 13, Article number: 274 (2020)
Trichomoniasis is the most common non-viral sexually transmitted disease caused by the protozoan parasite Trichomonas vaginalis. Metronidazole (MTZ) is a widely used drug for the treatment of trichomoniasis; however, increased resistance of the parasite to MTZ has emerged as a highly problematic public health issue.
We conducted iTRAQ-based analysis to profile the proteomes of MTZ-sensitive (MTZ-S) and MTZ-resistant (MTZ-R) parasites. STRING and gene set enrichment analysis (GESA) were utilized to explore the protein-protein interaction networks and enriched pathways of the differentially expressed proteins, respectively. Proteins potentially related to MTZ resistance were selected for functional validation.
A total of 3123 proteins were identified from the MTZ-S and MTZ-R proteomes in response to drug treatment. Among the identified proteins, 304 proteins were differentially expressed in the MTZ-R proteome, including 228 upregulated and 76 downregulated proteins. GSEA showed that the amino acid-related metabolism, including arginine, proline, alanine, aspartate, and glutamate are the most upregulated pathways in the MTZ-R proteome, whereas oxidative phosphorylation is the most downregulated pathway. Ten proteins categorized into the gene set of oxidative phosphorylation were ATP synthase subunit-related proteins. Drug resistance was further examined in MTZ-S parasites pretreated with the ATP synthase inhibitors oligomycin and bafilomycin A1, showing enhanced MTZ resistance and potential roles of ATP synthase in drug susceptibility.
We provide novel insights into previously unidentified proteins associated with MTZ resistance, paving the way for future development of new drugs against MTZ-refractory trichomoniasis.
Trichomonas vaginalis is a unicellular, flagellated eukaryote that causes the most widespread non-viral sexually transmitted infection (STI), with more than 275 million cases reported annually worldwide . Although men are often asymptomatic carriers of T. vaginalis infection, dysuria, discharge and an increased risk of prostate cancer have been reported . Infected women develop vaginitis, urethritis, and cervicitis, potentially leading to adverse pregnancy outcomes, such as infertility, preterm delivery, and low-birth-weight infants . Trichomoniasis has also been linked to an increased risk of human immunodeficiency virus (HIV) transmission  and cervical cancer . Metronidazole (MTZ) and other 5-nitroimidazoles are the only two drugs currently approved by the FDA for the treatment of trichomoniasis, but MTZ-resistant (MTZ-R) strains are on the rise [6, 7]. It has been shown that 4.3% of T. vaginalis isolates display MTZ resistance in the USA . Higher oral MTZ doses can sometimes cure refractory trichomoniasis but appear to be poorly tolerated . Additionally, the teratogenic effect of MTZ on animal models is well documented [10,11,12]. Hence, it is essential to discover alternative chemotherapeutic agents against MTZ-R T. vaginalis.
Several mechanisms have been proposed to elucidate how pathogens resist to MTZ treatment, such as altered reduction efficiency , drug inactivation , reduced drug uptake , active drug efflux , and increased DNA damage repair . Previous work showed that deficient pyruvate:ferredoxin oxidoreductase (PFOR) activity in T. vaginalis leads to low levels of anaerobic resistance to MTZ . Another report indicated that MTZ-R trichomonads have decreased ferredoxin (Fd) levels . A recent transcriptomic study in Entamoeba histolytica demonstrated that MTZ resistance is associated with specific transcriptional changes . These findings suggest that various mechanisms are involved in MTZ resistance in different protozoan parasites. It is worthwhile to advance our knowledge of the mechanism of MTZ resistance in T. vaginalis by using multi-omics approaches, such as next-generation RNA sequencing (RNA-seq) or quantitative proteomic analysis, which is helpful to identify the novel players regulating MTZ resistance.
Although the draft genome sequence of T. vaginalis has been reported , the progress in studying the biology of trichomonads by using proteomic technologies is slow. The initial proteomic study utilized 2-dimentional gel electrophoresis (2-DE) combined with matrix-assisted laser desorption ionization time-of-flight (MALDI-TOF) mass spectrometric analysis to establish the proteome reference map of T. vaginalis . Subsequent studies focused on the profiling of subcellular proteomes, such as the hydrogenosome  and surface proteomes . Comparative proteomic analyses were conducted to identify the potentially important proteins involved in the adherence of trophozoites to fibronectin , iron availability , and the transformation of trophozoites to pseudocysts . Additionally, differentially expressed proteins were identified from bovine and feline Tritrichomonas foetus genotypes, providing evidence of host-specific adaptation . A recent comparative proteomic analysis of the membrane proteins has been carried out in another trichomonad Trichomonas gallinae, showing distinct proteins associated with virulence .
As an increased resistance to MTZ has emerged as a highly problematic public health issue in several protozoan infections, we aim to unveil the molecular events of MTZ resistance in T. vaginalis using the isobaric tags for relative and absolute quantitation (iTRAQ)-based quantitative proteomic approach. The findings of this study will advance our understanding of the mechanism of MTZ resistance in the amitochondriate organisms and pave the way for development of alternative drugs against MTZ-refractory trichomoniasis.
Trichomonas vaginalis strains and culture conditions
The T. vaginalis MTZ-S strain (ATCC 30236) and MTZ-R strain (ATCC 50143) were maintained in YIS medium , pH 5.8, containing 10% heat-inactivated horse serum and 1% glucose at 37 °C. Growth of the parasites was monitored by using trypan blue exclusion hemocytometer counts.
MTZ susceptibility assay
To determine the concentration and incubation time of MTZ for proteomic analysis, the cell density of both MTZ-S and MTZ-R parasites was monitored every 4 h after treatment with different concentrations of MTZ (5, 10, 20 μM) (Sigma-Aldrich, Saint Louis, MO, USA) compared with that of the sterile distilled water (SDW)-treated control.
Sample preparation and iTRAQ labeling
Whole cell lysates extracted from the MTZ-S and MTZ-R parasites treated with MTZ (20 μM for 8 h) or SDW were harvested and resuspended in 100 µl RapiGest SF surfactant (Waters Corporation, MA, USA). The protein concentrations in samples were determined using a BCA protein assay kit (Pierce, Thermo Fisher Scientific, Illinois, USA). The protein samples (10 µg) were separated with 10% sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) (Future Scientific Co., Taoyuan, Taiwan) followed by Coomassie blue staining using standard procedures. Each sample (20 µg) was first reduced with 5 mM tris-(2-carboxyethyl)-phosphine (TCEP; Sigma-Aldrich) at 60 °C for 30 min, followed by cysteine-blocking with 10 mM methyl methanethiosulfonate (MMTS; Sigma-Aldrich) at 25 °C for 30 min, and digested at 37 °C for 16 h by trypsin (1 µg; Promega, Madison, USA) in a solution containing 150 mM TEABC. The peptides were then labeled with iTRAQ reagent (ABSciex, CA, USA) according to the manufacturer’s protocol. After incubation at room temperature for 1 h, the six peptide mixtures were pooled, dried by vacuum centrifugation, and store at − 80 °C until use.
LC-MS/MS analysis for protein identification
The dried peptide mixtures were reconstituted with 30% acetonitrile/0.1% formic acid for analysis using a LTQ-Orbitrap ELITE mass spectrometer (Thermo Fisher Scientific) as previously described . The MS raw data files were analyzed by Proteome Discoverer software (version 1.4; Thermo Fisher Scientific) including the reporter ions quantifier node for iTRAQ quantification. The MS/MS spectra was searched against the TrichDB-33_TvaginalisG3 sequence database (97,471 entries) using the Mascot search engine (version 2.5; Matrix Science, London, UK). For peptide identification, 10 ppm mass tolerance was permitted for intact peptide masses, and 0.05 Da for HCD fragment ions with an allowance for two missed cleavages made from the trypsin or semi-trypsin digestion: oxidized methionine, acetyl (protein N-terminal), iTRAQ (N-terminal), and iTRAQ (lysine) as variable modifications; and Methylthio (cysteine) as the fixed modifications. The peptide-spectrum match (PSM) was then filtered based on high confidence and Mascot search engine rank 1 of peptide identification to ensure an overall false discovery rate below 0.01. Proteins with a single peptide hit were removed. The identified proteins with fold changes higher than 2-fold or lower than 0.5-fold in the MTZ-R proteome upon MTZ treatment were considered differentially expressed in this study.
Proteomic data analysis
Differential protein expression analysis in the MTZ-R proteome in response to drug treatment compared with that of MTZ-S was performed using Gene Set Enrichment Analysis (GSEA) . Trichomonas vaginalis-specific databases with functional annotation of all identified proteins based on the Kyoto Encyclopedia of Genes and Genomes (KEGG)  and Gene Ontology (GO)  were incorporated into the GSEA software as reference databases. Protein-protein interaction networks were analyzed by STRING . The STRING database collects and integrates known and predicted protein-protein association data for a large number of organisms, including T. vaginalis. A protein identifier was inserted in STRING and the network view summarizes the network of predicted associations for a particular group of proteins.
RNA extraction, cDNA synthesis and quantitative PCR (qPCR)
Total RNA was extracted from MTZ-S and MTZ-R parasites treated with MTZ (20 μM) or SDW for 8 h using the SV total RNA isolation system (Premega). Reverse transcription (RT) was carried out using the SuperScriptTMIII First-Strand synthesis system (Thermo Fisher Scientific). qPCR was performed as previously described with minor modifications . Briefly, the qPCR was performed by using SYBR Green qPCR Master Mix (Bio-Rad, CA, USA) on a ViiA7 real-time PCR system (Thermo Fisher Scientific). 60S rRNA was used as an internal control for normalization of gene expression in all experimental groups. Primer pairs used in this study are listed in Additional file 1: Table S1.
ATP synthase inhibition assay
To verify the possible role of ATP synthase in MTZ resistance, MTZ-S and MTZ-R trophozoites were pretreated with different concentrations of oligomycin (1 μM, 5μM and 10 μM) (Sigma-Aldrich, Saint Louis, MO, USA) or bafilomycin A1 (100 nM and 500 nM) (Sigma-Aldrich, Saint Louis, MO, USA) for 1 h, and the susceptibility to MTZ was monitored as compared to the DMSO-treated control.
Quantitative data were expressed as mean ± SD of three independent experiments unless otherwise indicated. A Student’s t-test (two-tailed) was used to evaluate the significant differences between groups. P < 0.05 was considered statistically significant.
The effect of MTZ on the growth of MTZ-S and MTZ-R T. vaginalis
The workflow started from the extraction of proteins of MTZ-S and MTZ-R parasites treated with or without MTZ to quantitative proteomic analysis, followed by in silico pathway analysis and functional validation, is shown in Fig. 1. To determine the concentration and incubation period of MTZ that could differentiate the growth patterns between MTZ-S and MTZ-R parasites, the cell density of trophozoites treated with different concentrations of MTZ (5, 10 and 20 μM) was monitored for 24 h (Fig. 2a, b). Treatment with 5 μM MTZ maintained the survival of MTZ-S parasites at least within 24 h, whereas MTZ-R parasites can still replicate after 16 h of the same treatment. Treatment with 10 μM MTZ inhibited the growth of MTZ-R parasites and induced cell death of MTZ-S parasites after 20 h of treatment. The inhibitory effect of 20 μM MTZ treatment on the growth of MTZ-R parasites was similar to that of 10 μM MTZ treatment, but this dosage can lead to rapid cell death of MTZ-S parasites after 8 h of treatment. Since 8 h of 20 μM MTZ treatment was able to differentiate the growth patterns between MTZ-S and MTZ-R parasites, we used this condition for subsequent proteomic analysis.
Quantitative proteomic analysis of MTZ-S and MTZ-R isolates
To identify the differentially expressed proteins associated with MTZ resistance, whole cell lysates extracted from MTZ-S and MTZ-R trophozoites treated with MTZ (20 μM for 8 h) were subjected to quantitative proteomic profiling using iTRAQ combined with LC-MS/MS analysis (Fig. 2c). Proteomic profiling of MTZ-S and MTZ-R trophozoites without drug treatment was also analyzed, representing strain-specific protein expression patterns. A total of 3123 and 3125 proteins were identified from both isolates treated with or without MTZ, respectively (Fig. 2d, Additional file 2: Table S2). Among the identified proteins, we obtained 304 differentially expressed proteins (9.7% of the identified MTZ-R/MTZ-S proteomes) in the MTZ-R proteome upon MTZ treatment, including 76 downregulated and 228 upregulated proteins. Additionally, we obtained 244 differentially expressed proteins (7.8% of the identified MTZ-R/MTZ-S proteomes) in the MTZ-R proteome without MTZ exposure, including 124 downregulated and 120 upregulated proteins. To validate the proteomic data, we determined the mRNA expression levels of 4 genes with 2-fold upregulation or downregulation in the MTZ-R proteome upon MTZ treatment (Additional file 1: Table S1). qPCR analysis confirmed the expression patterns of these transcripts, which were in agreement with the protein expression patterns (Additional file 3: Figure S1).
Protein-protein interaction network analysis of the most differentially expressed proteins in the MTZ-R proteome upon MTZ treatment
Among the 304 differentially expressed proteins between the MTZ-R and MTZ-S proteomes in response to MTZ treatment, 21 proteins (0.7% of the identified MTZ-R/MTZ-S proteomes) showed more than a 5-fold increase in MTZ-R parasites (Table 1), whereas 13 proteins (0.4% of the MTZ-R/MTZ-S proteomes) showed more than a 5-fold decrease in MTZ-R parasites (Table 2). The most upregulated proteins in MTZ-R parasites treated with MTZ were ankyrin 2,3/unc44 (TVAG_465440, 29.25-fold), receptor expression-enhancing protein (TVAG_232140, 11.46-fold), and histidinol-phosphate aminotransferase (TVAG_108400, 9.76-fold). On the other hand, the most downregulated proteins in MTZ-R parasites treated with MTZ were glucose kinase (TVAG_442070, 0.09-fold), iron-sulfur flavoprotein (TVAG_370510, 0.11-fold), and NAD(P)H dehydrogenase (TVAG_311580, 0.13-fold). Protein-protein interaction networks of these significantly dysregulated proteins in the MTZ-R proteome upon drug treatment were further constructed by STRING. The enrichment tests for functional associations have been incorporated in STRING, such as GO, KEGG, Pfam and InterPro. Ankyrin 2,3/unc44 was predicted to interact with several kinases (Fig. 3a). Receptor expression-enhancing protein was shown to be associated with a cluster of proteins with hydrolase activity and GDP-binding site on the ER membrane compartment (Fig. 3b). Histidinol-phosphate aminotransferase and glucose kinase were linked to metabolic pathways (Fig. 3c, d). It is noteworthy to mention that iron-sulfur flavoprotein and its interactive proteins with the iron hydrogenase subunit have also been identified to be downregulated in other MTZ-R T. vaginalis stains  (Fig. 3e). NAD(P)H dehydrogenase was associated with the proteins with 2F-2S iron sulfur cluster binding domain (Fig. 3f).
Differentially expressed proteins involved in the hydrogenosomal metabolism of the MTZ-R proteome upon MTZ treatment
It has been shown that laboratory-generated MTZ resistance is associated with downregulation of specific hydrogenosomal enzymes that reduce MTZ, such as PFOR and Fd, but the clinical resistant strains showed no decrease in PFOR or Fd transcription [37, 38]. We thus analyzed the differentially expressed proteins involved in hydrogenosomal metabolism  between the MTZ-R and MTZ-S proteomes in the presence or absence of MTZ (Fig. 4, Additional file 4: Table S3). Two of five PFOR proteins (TVAG_230580, TVAG_242960) were downregulated in the MTZ-R proteome upon drug treatment. Among the four identified Fd proteins, Fd1 (TVAG_003900, 0.635-fold) was the most downregulated protein in the MTZ-proteome compared with that of MTZ-S in response to drug treatment. Particularly noteworthy is that all identified succinate thiokinase (STK) (TVAG_183500, TVAG_144730, TVAG_259190, TVAG_047890 and TVAG_318670), which catalyze ATP synthesis in the hydrogenosome, were significantly downregulated in the MTZ-R proteome in response to drug treatment. This suggests that the hydrogenosomal energy metabolism was suppressed in MTZ-R parasites upon drug treatment.
GSEA revealed altered metabolic signatures in the MTZ-R proteome after drug treatment
To globally assess the differences between the MTZ-R and MTZ-S proteomes upon drug treatment, we utilized Gene Set Enrichment Analysis (GSEA) to identify the enriched protein sets, both positive and negative regulation, according to KEGG (Additional file 5: Table S4, Additional file 6: Table S5) and GO (Additional file 7: Table S6, Additional file 8: Table S7) functional annotations. GSEA for KEGG pathway mapping showed that the most upregulated pathways in the MTZ-R proteome upon drug treatment were arginine and proline metabolism (enrichment score, ES = 0.66), alanine, aspartate, and glutamate metabolism (ES = 0.65), and ribosome biogenesis in eukaryotes (ES = 0.44) (Fig. 5a), whereas the most downregulated pathways were oxidative phosphorylation (ES = − 0.89), citrate cycle (ES = − 0.58), and phagosome (ES = − 0.45) (Fig. 5b). GSEA based on GO categories revealed that the most upregulated pathways in the MTZ-R proteome upon drug treatment were actin binding (ES = 0.63), cysteine type peptidase activity (ES = 0.61), and response to heat (ES = 0.58) (Fig. 5c), whereas the most downregulated protein sets belonged to protein disulfide oxidoreductase activity (ES = − 0.52), glycolysis (ES = − 0.43), and threonine type endopeptidase activity (ES = − 0.43) (Fig. 5d). Based on the GSEA data enriched from the MTZ-R proteome after drug treatment, it seems that there was a metabolic reprogramming toward reduced glycolysis but increased amino acid metabolism in MTZ-R parasites.
Several aminotransferases were upregulated in the MTZ-R proteome upon MTZ exposure
We further analyzed the individual dysregulated proteins in these enriched protein sets. Specifically, several aminotransferases categorized into arginine, proline, aspartate, and glutamate metabolism, such as glucosamine-fructose-6-phosphate aminotransferase (TVAG_478320 and TVAG_011780), aspartate aminotransferase (TVAG_431780, TVAG_430210 and TVAG_419720), alanine aminotransferase (TVAG_088220), and glutamate dehydrogenase (GDH) (TVAG_025910) were upregulated and statistically enriched in the MTZ-R proteome after drug treatment. On the other hand, several proteins involved in glycolysis that link to hydrogenosomal energy metabolism, such as glyceraldehyde 3-phosphate dehydrogenase (TVAG_193800), phosphoglycerate mutase (TVAG_209020), and PFOR (TVAG_230580 and TVAG_242960), were downregulated in the MTZ-R proteome after drug treatment. These results suggest that specific reactions catalyzed by these aminotransferases are activated in MTZ-R parasites in response to MTZ; however, it requires further investigation whether this metabolic reprogramming is an outcome or a cause of MTZ resistance.
Potential roles of ATP synthase subunit-related proteins in MTZ susceptibility
Oxidative phosphorylation was the most downregulated pathway in the MTZ-R proteome in response to drug treatment. It is noteworthy that most proteins categorized into the oxidative phosphorylation protein set were ATP synthase subunit-related proteins (Table 3), suggesting their potential roles in MTZ resistance. In the T. vaginalis genome, the genes encoding the subunits of putative V-type ATPase (A to H, a, c and d) have been identified , but their biological functions have not yet been investigated. Given that all the subunits of ATPase were significantly downregulated in the MTZ-R proteome in response to MTZ, we postulate that these ATP synthase-related proteins may sensitize trophozoites to drug treatment. We thus examined whether the ATP synthase inhibitors oligomycin and bafilomycin A1 (BafA1) are able to enhance the resistance of MTZ-S parasites to MTZ. The effects of MTZ on the growth of MTZ-S and MTZ-R parasites pretreated with different concentrations of oligomycin (1, 5 and 10 μM) (Fig. 6a, b) and BafA1 (100 nM and 500 nM) (Fig. 6c, d) were monitored. Notably, MTZ-S parasites pretreated with oligomycin and BafA1 enhanced the resistance to MTZ compared with the DMSO-treated control. However, the effects of ATP synthase inhibitors on the growth of MTZ-R following MTZ treatment were not observed. Together, these results suggest that inhibition of ATP synthase-related proteins may play a crucial role in MTZ resistance of T. vaginalis.
MTZ resistance has been observed in 4.3–9.6% of clinical isolates of T. vaginalis in the USA [8, 39]; however, the underlying mechanism is far from understood. In the current study, we sought to advance our understanding of the global changes between the MTZ-R and MTZ-S proteomes in response to drug treatment, unraveling novel molecules or pathways associated with MTZ resistance. To our knowledge, this is the first quantitative proteomic analysis of MTZ resistance in trichomonads, showing many previously uncharacterized proteins with a potential impact on MTZ resistance.
A recent study demonstrated that the intracellular iron status is associated with early progression of MTZ resistance in T. vaginalis, but no significant changes in transcription of PFOR and Fd was observed between the MTZ-S and MTZ-R strains . Additionally, no marked difference in PFOR mRNA levels is detected between the HM-1 and HM-1 derived MTZ-R strains of E. histolytica . Our results showed inconsistent protein expression patterns of different isoforms of PFOR and Fd in the MTZ-R proteome. These findings suggest that multiple molecules or pathways are involved in MTZ resistance and that the underlying mechanism is regulated by different vaginal microenvironments. We found that all hydrogenosomal succinate STK family proteins are downregulated in the MTZ-R proteome after MTZ treatment, which is not consistent with the previous finding that shows almost the same protein expression levels of α-subunits of STK in the MTZ-S strain and its drug-resistant derivatives . Hence, it is likely that various mechanisms of MTZ resistance exist in different isolates of T. vaginalis, which is similar to the recent study showing that the majority of differentially expressed proteins are unique to each of the three MTZ-R strains of Giardia lamblia . Additionally, a transcriptomic study on MTZ resistance in tens of T. vaginalis strains  also demonstrates that their TrxR gene expression results are opposite to the expression pattern shown in previous studies .
Ankyrin 2,3/unc44 is the most upregulated protein in the MTZ-R proteome upon drug treatment, with about 30-fold higher expression compared to that in the MTZ-S proteome. It has been demonstrated that several bacterial and viral ankyrin repeat (ANK)-containing proteins play a role in host-pathogen interaction and the evolution of infectious diseases , but the biological significance of the ANK family proteins have not yet been well characterized in protists. A previous study reported that a gene homologous to ANKRD1 (ankyrin repeat domain 1[cardiac muscle]) is associated with the platinum sensitivity in ovarian cancer cell lines, and decreasing ANKRD1 expression is a potential strategy to sensitize tumors to platinum-based drugs . Additionally, ankyrin binding to the multidrug transporter MDR1 (P-glycoprotein) results in the efflux of chemotherapeutic drugs in both human breast and ovarian tumor cells . We have also identified two multidrug resistance proteins (TVAG_542450 and TVAG_542460) in the MTZ-R proteome after drug treatment (Additional file 2: Table S2). It is worth clarifying whether the identified MDR proteins are involved in MTZ resistance in T. vaginalis.
We found that several aminotransferases were upregulated in the MTZ-R proteome in response to MTZ. NAD-specific GDH, one of the aminotransferases that catalyzes the oxidative deamination of glutamate to α-ketoglutarate, has been shown to be upregulated as a part of the stress response in the MTZ-R strain of E. histolytica . Our previous transcriptomic analysis in T. vaginalis also demonstrated the overexpression of aminotransferases upon glucose restriction, especially for GDH, which may serve as an energy compensation mechanism caused by glucose starvation . Since we observed that the gene set of glycolysis was statistically downregulated and all the hydrogenosomal succinate STK family proteins that catalyze ATP synthesis were repressed in the MTZ-R proteome upon drug treatment, it is possible that the overexpression of aminotransferases also plays a compensatory role for energy production and parasite survival. It remains to be determined whether downregulation of hydrogenosomal STK in MTZ-R parasites upon MTZ exposure involves in the regulation of drug resistance in T. vaginalis.
Using GSEA, we have demonstrated that the oxidative phosphorylation protein set that contains ten ATP synthase subunit-related proteins was the most downregulated pathway in the MTZ-R proteome upon drug treatment. The drug resistance of MTZ-S parasites has been proved to be enhanced after pretreatment with two ATP synthase inhibitors oligomycin and bafilomycin A1, suggesting the potential association of these putative ATP synthases with drug susceptibility in T. vaginalis. ATP synthase is located on the membranes of mitochondria, bacteria, and chloroplast thylakoids as well as on the surfaces of various cell types, such as endothelial cells [47, 48], keratinocytes , and adipocytes . It is noteworthy that downregulation of ATP synthase components has been reported in the majority of cancers [51, 52] and shown to be associated with chemotherapy resistance . A recent study mechanistically demonstrated that loss of the ATP synthase subunit ATP5H is strongly linked to multimodal cancer therapy resistance and poor survival in cancer patients, and the mechanism is mediated by ROS accumulation and HIF-1 α stabilization that confers to tumor cells a stem-like and invasive phenotype . A previous study in Trypanosoma brucei indicated that the V-ATPase subunit depleted strains display more than a 100-fold increase in the half-maximal effective concentration for the trypanocidal drug isometamidium chloride (ISM) . Additionally, bafilomycin A1, a specific inhibitor of the V-ATPase that binds to the c subunit of V0, strongly antagonized ISM action. These findings conclude that V-ATPase sensitizes T. brucei to ISM. Trichomonas vaginalis lacks conventional mitochondria and instead contains divergent mitochondrial-related organelles called hydrogenosomes; however, the aforementioned ATP synthase subunit-related proteins have not been identified in the hydrogenosomal proteome [23, 26]. Some of the ATPase-related proteins have been identified by surface proteomic analysis of T. vaginalis, but their biological significance has not yet been verified. It is worthwhile to study which subunit of ATP synthase plays the most crucial role in MTZ resistance by gene-specific silencing. On the other hand, it will be also feasible to overexpress the specific ATP synthase in MTZ-R parasites and verify their susceptibility to MTZ. Together, our results suggest that loss of ATP synthase subunit-related proteins may be associated with MTZ resistance in T. vaginalis, with the potential to overcome drug resistance in trichomonads.
Collectively, to the best of our knowledge, this study revealed for the first time the global proteomic changes between the MTZ-R and MTZ-S strains following drug treatment, enhancing our understanding of the molecular events potentially involved in MTZ resistance. We have highlighted the significantly enriched pathways as well as the specific proteins associated with MTZ resistance. Our results demonstrate that several aminotransferases categorized into amino acid metabolism were upregulated, whereas several proteins involved in glycolysis and hydrogenosomal energy metabolism were downregulated in the MTZ-R proteome in response to drug treatment. Specifically, ATP synthase subunit-related proteins classified into the oxidative phosphorylation gene set were the most downregulated family proteins in the MTZ-R proteome, and functional validation by ATP synthase inhibition assay supported the potential involvement of these proteins in MTZ susceptibility. We provide novel insights into the MTZ resistance mechanism at the proteome level, paving the way for future development of new drug targets against MTZ-refractory trichomoniasis.
Availability of data and materials
Data supporting the conclusions of this article are included within the article and its additional files. The datasets used and/or analyzed during the present study will be made available by the corresponding author upon reasonable request. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via PRIDE partner repository with the dataset identifier PXD018522.
sodium dodecyl sulfate polyacrylamide gel electrophoresis
WHO. Global incidence and prevalence of selected curable sexually transmitted infection. Geneva: World Health Organization; 2012.
Sutcliffe S, Alderete JF, Till C, Goodman PJ, Hsing AW, Zenilman JM, et al. Trichomonosis and subsequent risk of prostate cancer in the Prostate Cancer Prevention Trial. Int J Cancer. 2009;124:2082–7.
Leitsch D. Recent advances in the Trichomonas vaginalis field. F1000Res. 2016;5:162.
Sorvillo F, Kerndt P. Trichomonas vaginalis and amplification of HIV-1 transmission. Lancet. 1998;351:213–4.
Yang S, Zhao W, Wang H, Wang Y, Li J, Wu X. Trichomonas vaginalis infection-associated risk of cervical cancer: a meta-analysis. Eur J Obstet Gynecol Reprod Biol. 2018;228:166–73.
Schmid G, Narcisi E, Mosure D, Secor WE, Higgins J, Moreno H. Prevalence of metronidazole-resistant Trichomonas vaginalis in a gynecology clinic. J Reprod Med. 2001;46:545–9.
Crowell AL, Sanders-Lewis KA, Secor WE. In vitro metronidazole and tinidazole activities against metronidazole-resistant strains of Trichomonas vaginalis. Antimicrob Agents Chemother. 2003;47:1407–9.
Kirkcaldy RD, Augostini P, Asbel LE, Bernstein KT, Kerani RP, Mettenbrink CJ, et al. Trichomonas vaginalis antimicrobial drug resistance in 6 US cities, STD Surveillance Network, 2009–2010. Emerg Infect Dis. 2012;18:939–43.
Sobel JD, Nyirjesy P, Brown W. Tinidazole therapy for metronidazole-resistant vaginal trichomoniasis. Clin Infect Dis. 2001;33:1341–6.
Cantu JM, Garcia-Cruz D. Midline facial defect as a teratogenic effect of metronidazole. Birth Defects Orig Artic Ser. 1982;18:85–8.
Burtin P, Taddio A, Ariburnu O, Einarson TR, Koren G. Safety of metronidazole in pregnancy: a meta-analysis. Am J Obstet Gynecol. 1995;172:525–9.
Menegola E, Broccia ML, Di Renzo F, Massa V, Giavini E. Craniofacial and axial skeletal defects induced by the fungicide triadimefon in the mouse. Birth Defects Res B Dev Reprod Toxicol. 2005;74:185–95.
Leiros HK, Kozielski-Stuhrmann S, Kapp U, Terradot L, Leonard GA, McSweeney SM. Structural basis of 5-nitroimidazole antibiotic resistance: the crystal structure of NimA from Deinococcus radiodurans. J Biol Chem. 2004;279:55840–9.
Ralph ED, Clarke DA. Inactivation of metronidazole by anaerobic and aerobic bacteria. Antimicrob Agents Chemother. 1978;14:377–83.
Lacey SL, Moss SF, Taylor GW. Metronidazole uptake by sensitive and resistant isolates of Helicobacter pylori. J Antimicrob Chemother. 1993;32:393–400.
Pumbwe L, Chang A, Smith RL, Wexler HM. BmeRABC5 is a multidrug efflux system that can confer metronidazole resistance in Bacteroides fragilis. Microb Drug Resist. 2007;13:96–101.
Land KM, Johnson PJ. Molecular basis of metronidazole resistance in pathogenic bacteria and protozoa. Drug Resist Updat. 1999;2:289–94.
Rasoloson D, Vanacova S, Tomkova E, Razga J, Hrdy I, Tachezy J, et al. Mechanisms of in vitro development of resistance to metronidazole in Trichomonas vaginalis. Microbiology. 2002;148:2467–77.
Yarlett N, Yarlett NC, Lloyd D. Ferredoxin-dependent reduction of nitroimidazole derivatives in drug-resistant and susceptible strains of Trichomonas vaginalis. Biochem Pharmacol. 1986;35:1703–8.
Penuliar GM, Nakada-Tsukui K, Nozaki T. Phenotypic and transcriptional profiling in Entamoeba histolytica reveal costs to fitness and adaptive responses associated with metronidazole resistance. Front Microbiol. 2015;6:354.
Carlton JM, Hirt RP, Silva JC, Delcher AL, Schatz M, Zhao Q, et al. Draft genome sequence of the sexually transmitted pathogen Trichomonas vaginalis. Science. 2007;315:207–12.
Huang KY, Chien KY, Lin YC, Hsu WM, Fong IK, Huang PJ, et al. A proteome reference map of Trichomonas vaginalis. Parasitol Res. 2009;104:927–33.
Schneider RE, Brown MT, Shiflett AM, Dyall SD, Hayes RD, Xie Y, et al. The Trichomonas vaginalis hydrogenosome proteome is highly reduced relative to mitochondria, yet complex compared with mitosomes. Int J Parasitol. 2011;41:1421–34.
de Miguel N, Lustig G, Twu O, Chattopadhyay A, Wohlschlegel JA, Johnson PJ. Proteome analysis of the surface of Trichomonas vaginalis reveals novel proteins and strain-dependent differential expression. Mol Cell Proteomics. 2010;9:1554–66.
Huang KY, Huang PJ, Ku FM, Lin R, Alderete JF, Tang P. Comparative transcriptomic and proteomic analyses of Trichomonas vaginalis following adherence to fibronectin. Infect Immun. 2012;80:3900–11.
Beltran NC, Horvathova L, Jedelsky PL, Sedinova M, Rada P, Marcincikova M, et al. Iron-induced changes in the proteome of Trichomonas vaginalis hydrogenosomes. PLoS One. 2013;8:e65148.
Dias-Lopes G, Wisniewski JR, de Souza NP, Vidal VE, Padron G, Britto C, et al. In-depth quantitative proteomic analysis of trophozoites and pseudocysts of Trichomonas vaginalis. J Proteome Res. 2018;17:3704–18.
Stroud LJ, Slapeta J, Padula MP, Druery D, Tsiotsioras G, Coorssen JR, et al. Comparative proteomic analysis of two pathogenic Tritrichomonas foetus genotypes: there is more to the proteome than meets the eye. Int J Parasitol. 2017;47:203–13.
Martinez-Herrero MDC, Garijo-Toledo MM, Gonzalez F, Bilic I, Liebhart D, Ganas P, et al. Membrane associated proteins of two Trichomonas gallinae clones vary with the virulence. PLoS ONE. 2019;14:e0224032.
Diamond LS, Clark CG, Cunnick CC. YI-S, a casein-free medium for axenic cultivation of Entamoeba histolytica, related Entamoeba, Giardia intestinalis and Trichomonas vaginalis. J Eukaryot Microbiol. 1995;42(3):277–8.
Michalski A, Damoc E, Lange O, Denisov E, Nolting D, Muller M, et al. Ultra high resolution linear ion trap orbitrap mass spectrometer (orbitrap elite) facilitates top down LC MS/MS and versatile peptide fragmentation modes. Mol Cell Proteomics. 2012;11(O111):013698.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005;102:15545–50.
Wixon J, Kell D. The Kyoto Encyclopedia of genes and genomes—KEGG. Yeast. 2000;17:48–55.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000;25:2–9.
Szklarczyk D, Morris JH, Cook H, Kuhn M, Wyder S, Simonovic M, et al. The STRING database in 2017: quality-controlled protein-protein association networks, made broadly accessible. Nucleic Acids Res. 2017;45:D362–8.
Bradic M, Warring SD, Tooley GE, Scheid P, Secor WE, Land KM, et al. Genetic indicators of drug resistance in the highly repetitive genome of Trichomonas vaginalis. Genome Biol Evol. 2017;9:1658–72.
Mead JR, Fernadez M, Romagnoli PA, Secor WE. Use of Trichomonas vaginalis clinical isolates to evaluate correlation of gene expression and metronidazole resistance. J Parasitol. 2006;92:196–9.
Wright JM, Webb RI, O’Donoghue P, Upcroft P, Upcroft JA. Hydrogenosomes of laboratory-induced metronidazole-resistant Trichomonas vaginalis lines are downsized while those from clinically metronidazole-resistant isolates are not. J Eukaryot Microbiol. 2010;57:171–6.
Schwebke JR, Barrientes FJ. Prevalence of Trichomonas vaginalis isolates with resistance to metronidazole and tinidazole. Antimicrob Agents Chemother. 2006;50:4209–10.
Argaez-Correa W, Alvarez-Sanchez ME, Arana-Argaez VE, Ramirez-Camacho MA, Novelo-Castilla JS, Coral-Martinez TI, et al. The role of iron status in the early progression of metronidazole resistance in Trichomonas vaginalis under microaerophilic conditions. J Eukaryot Microbiol. 2018;66:309–15.
Emery SJ, Baker L, Ansell BRE, Mirzaei M, Haynes PA, McConville MJ, et al. Differential protein expression and post-translational modifications in metronidazole-resistant Giardia duodenalis. Gigascience. 2018;7:giy024.
Leitsch D, Kolarich D, Binder M, Stadlmann J, Altmann F, Duchene M. Trichomonas vaginalis: metronidazole and other nitroimidazole drugs are reduced by the flavin enzyme thioredoxin reductase and disrupt the cellular redox system. Implications for nitroimidazole toxicity and resistance. Mol Microbiol. 2009;72:518–36.
Al-Khodor S, Price CT, Kalia A, Abu Kwaik Y. Functional diversity of ankyrin repeats in microbial proteins. Trends Microbiol. 2010;18:132–9.
Scurr LL, Guminski AD, Chiew YE, Balleine RL, Sharma R, Lei Y, et al. Ankyrin repeat domain 1, ANKRD1, a novel determinant of cisplatin sensitivity expressed in ovarian cancer. Clin Cancer Res. 2008;14:6924–32.
Bourguignon LY, Peyrollier K, Xia W, Gilad E. Hyaluronan-CD44 interaction activates stem cell marker Nanog, Stat-3-mediated MDR1 gene expression, and ankyrin-regulated multidrug efflux in breast and ovarian tumor cells. J Biol Chem. 2008;283:17635–51.
Huang KY, Chen YY, Fang YK, Cheng WH, Cheng CC, Chen YC, et al. Adaptive responses to glucose restriction enhance cell survival, antioxidant capability, and autophagy of the protozoan parasite Trichomonas vaginalis. Biochim Biophys Acta. 2014;1840:53–64.
Moser TL, Kenan DJ, Ashley TA, Roy JA, Goodman MD, Misra UK, et al. Endothelial cell surface F1-F0 ATP synthase is active in ATP synthesis and is inhibited by angiostatin. Proc Natl Acad Sci USA. 2001;98:6656–61.
Moser TL, Stack MS, Asplin I, Enghild JJ, Hojrup P, Everitt L, et al. Angiostatin binds ATP synthase on the surface of human endothelial cells. Proc Natl Acad Sci USA. 1999;96:2811–6.
Burrell HE, Wlodarski B, Foster BJ, Buckley KA, Sharpe GR, Quayle JM, et al. Human keratinocytes release ATP and utilize three mechanisms for nucleotide interconversion at the cell surface. J Biol Chem. 2005;280:29667–76.
Kim BW, Choo HJ, Lee JW, Kim JH, Ko YG. Extracellular ATP is generated by ATP synthase complex in adipocyte lipid rafts. Exp Mol Med. 2004;36:476–85.
Cuezva JM, Krajewska M, de Heredia ML, Krajewski S, Santamaria G, Kim H, et al. The bioenergetic signature of cancer: a marker of tumor progression. Cancer Res. 2002;62:6674–81.
Isidoro A, Martinez M, Fernandez PL, Ortega AD, Santamaria G, Chamorro M, et al. Alteration of the bioenergetic phenotype of mitochondria is a hallmark of breast, gastric, lung and oesophageal cancer. Biochem J. 2004;378:17–20.
Shin YK, Yoo BC, Chang HJ, Jeon E, Hong SH, Jung MS, et al. Down-regulation of mitochondrial F1F0-ATP synthase in human colon cancer cells with induced 5-fluorouracil resistance. Cancer Res. 2005;65:3162–70.
Song KH, Kim JH, Lee YH, Bae HC, Lee HJ, Woo SR, et al. Mitochondrial reprogramming via ATP5H loss promotes multimodal cancer therapy resistance. J Clin Invest. 2018;128:4098–114.
Baker N, Hamilton G, Wilkes JM, Hutchinson S, Barrett MP, Horn D. Vacuolar ATPase depletion affects mitochondrial ATPase function, kinetoplast dependency, and drug sensitivity in trypanosomes. Proc Natl Acad Sci USA. 2015;112:9112–7.
We would like to thank the “Molecular Medicine Research Center, Chang Gung University” from The Featured Areas Research Center Program within the framework of the Higher Education Sprout Project by the Ministry of Education (MOE) in Taiwan.
This study was supported by grants from the Ministry of Science and Technology, Taiwan (MOST 107-2320-B-016-008-MY3) to KYH and (MOST 107-2320-B-182-028-MY3) to LJC, Tri-Service General Hospital Songshan Branch, Taiwan (10810 and 10818) to HAL and JYW, Tri-Service General Hospital, Taiwan (TSGH-C108-213) to RMC, Chang Gung Memorial Hospital, Linkou, Taiwan (CMRPD1H0251) to JSY.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 3: Figure S1. Validation of the proteomic data by qPCR analysis. Four genes with 2-fold upregulation (receptor expression-enhancing protein and histidinol-phosphate aminotransferase) (a) and downregulation (glucose kinase and iron-sulfur flavoprotein) (b) in the MTZ-R proteome upon MTZ treatment compared with the MTZ-S proteome were confirmed by qPCR analysis. *P < 0.05, ** P < 0.01, *** P < 0.001.
About this article
Cite this article
Lin, H., Chu, L.J., Huang, P. et al. Proteomic signatures of metronidazole-resistant Trichomonas vaginalis reveal novel proteins associated with drug resistance. Parasites Vectors 13, 274 (2020). https://doi.org/10.1186/s13071-020-04148-5
- Trichomonas vaginalis
- Metronidazole resistance