Next Article in Journal
Analyzing the Interplay between COVID-19 Viral Load, Inflammatory Markers, and Lymphocyte Subpopulations on the Development of Long COVID
Previous Article in Journal
Co-Circulation of SARS-CoV-2 and Other Respiratory Pathogens in Upper and Lower Respiratory Tracts during Influenza Season 2022–2023 in Lazio Region
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Microorganisms
Volume 11
Issue 9
10.3390/microorganisms11092240
microorganisms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genetic Variability of the Internal Transcribed Spacer and Pyruvate:Ferredoxin Oxidoreductase Partial Gene of Trichomonas vaginalis from Female Patients

by
Fernando Martinez-Hernandez
1,
Fabiola Sanchez-Aguillon
2,
Joel Martinez-Ocaña
1,
Nelly Raquel Gonzalez-Arenas
1,
Mirza Romero-Valdovinos
3,
Eduardo Lopez-Escamilla
1,
Pablo Maravilla
1,* and
Guiehdani Villalobos
1,4,*
1
Departamento de Ecologia de Agentes Patogenos, Hospital General “Dr. Manuel Gea Gonzalez”, Mexico City 14080, Mexico
2
Laboratorio de Investigación del Departamento de Biologia Molecular e Histocompatibilidad, Hospital General “Dr. Manuel Gea Gonzalez”, Mexico City 14080, Mexico
3
Laboratorio de Patogenos Emergentes, Departamento de Biologia Molecular e Histocompatibilidad, Hospital General “Dr. Manuel Gea Gonzalez”, Mexico City 14080, Mexico
4
Laboratorio de Biologia Molecular del Departamento de Produccion Agricola y Animal, Universidad Autonoma Metropolitana, Mexico City 04960, Mexico
*
Authors to whom correspondence should be addressed.
Microorganisms 2023, 11(9), 2240; https://doi.org/10.3390/microorganisms11092240
Submission received: 10 August 2023 / Revised: 31 August 2023 / Accepted: 4 September 2023 / Published: 5 September 2023
(This article belongs to the Section Parasitology)
Browse Figures
Review Reports Versions Notes

Abstract

:
In the present study, we evaluated the genetic variability of the internal transcribed spacer (ITS) region and the pyruvate:ferredoxin oxidoreductase (pfor) A gene of Trichomonas vaginalis from female patients and its possible implications in the host–parasite relationship. Phylogenetic and genetics of populations analyses were performed by analyzing sequences of the ITS region and partial pfor A gene of clinical samples with T. vaginalis, as previously documented. Alignments of protein sequences and prediction of three-dimensional structure were also performed. Although no correlation between the main clinical characteristics of the samples and the results of phylogeny was found, a median-joining analysis of ITS haplotypes showed two main clusters. Also, pfor A, due to its phylogenetic divergence, could be used as a marker to confirm the genus and species of trichomonads. Alignment of protein sequences and prediction of three-dimensional structure showed that PFOR A had a highly conserved structure with two synonymous mutations in the PFOR domain, substituting a V for a G or a S for a P. Our results suggest that the role of genetic variability of PFOR and ITS may not be significant in the symptomatology of this pathogen; however, their utility as genus and species markers in trichomonads is promising.

1. Introduction

Trichomoniasis by Trichomonas vaginalis is the most common nonviral sexually transmitted infection (STI) in humans worldwide; in 2016, the WHO estimated 156 million new infections [1]. T. vaginalis infection in females is symptomatic in approximately 50% of cases, and approximately 30% of patients with asymptomatic cases develop some symptoms in the 6-month period postinfection [2]. A wide variety of clinical symptoms have been described in women, such as urethral discharge, itching, abdominal pain and dysuria, which can lead to premature birth and low birth weight in newborns. In contrast, infection of the male genitourinary tract is generally asymptomatic, although mild urethritis, epididymitis and prostatitis can occur [3,4].
Under both aerobic and anaerobic conditions, T. vaginalis uses carbohydrates as its main energy source via fermentative metabolism. Carbohydrate metabolism occurs in the cytoplasm and hydrogenosome, which is analogous to mitochondria and is the site of fermentative oxidation of pyruvate [4]. Anaerobic energy metabolism is accomplished through catalyzed reactions by several key enzymes, such as pyruvate:ferredoxin oxidoreductase (PFOR) [5]. In addition, the draft of the T. vaginalis genome sequence revealed pfor-like genes such as pfor A, pfor BI, pfor BII and pfor E [6]; these genes are available in the recent released T. vaginalis genome sequence (GenBank project number PRJNA885811). In different organisms, anaerobic energy enzymes are known to function in different subcellular compartments in such a way that they perform divergent functions as a consequence of their localization and environmental conditions [7,8], as observed for PFOR in T. vaginalis [9], which is located on the membranes of hydrogenosomes and is responsible for the decarboxylation of pyruvate to acetyl-CoA [10]. Although the catalytic mechanism of PFOR is not fully understood, it is assumed that Fe-S groups are involved in the catalytic mechanism of the enzyme. In bacteria, preserved motifs of Fe-S centers with cysteines (Cys), which coordinate the Fe-S centers present in the enzyme, have been described [5]. In addition, this enzyme is an example of a surface-associated cell-binding protein that has enzyme activity and is involved in cytoadherence [11]. Parasite adherence to vaginal epithelial cells is fundamental for initiating and maintaining Trichomonas infection and for survival in the host [12]. PFOR is involved in proliferation, adhesion to host cells and abscess formation, so it could influence the pathogenicity of this parasite [9,12,13]. Interestingly, this enzyme is present in protozoan parasites but is absent in humans [14]. It is important to highlight that, to date, PFOR studies have been focused on their characterization and function depending on their localization, but whether pfor genes present genetic variations and are associated with pathogenicity or the development of specific symptoms have not been described.
On the other hand, the internal transcribed spacers (ITSs), ITS1 and ITS2, are spacer regions divided by the 5.8S rRNA gene (ITS1+5.8S+ITS2); these loci evolve much faster than coding regions because substitutions occurring in spacers may be considered neutral mutations without any constraints. Thus, there is a large amount of information related to the usefulness of ITS, mainly clarifying problematic taxa, such as cryptic species [15]. Although there are some reports about the usefulness of ITS1-5.8S rRNA-ITS2 as a marker for differentiation among T. vaginalis isolates, information obtained using this locus is still scarce [16,17,18,19]. Therefore, the objective of the present study was to evaluate the genetic variability of the ITS1+5.8S+ITS2 region and the pyruvate:ferredoxin oxidoreductase (pfor) A gene of T. vaginalis from female patients and its possible implications in the host–parasite relationship.

2. Materials and Methods

2.1. Ethical Approval and Sampling

The Ethics and Research Committees at Hospital General “Dr. Manuel Gea Gonzalez” (HGMGG) approved this study, with reference number 12-75-2015. Recovered samples, used primarily for pathogen diagnosis of 28 women with trichomoniasis who attended the gynecology service of the HGMGG between 2015 and 2019, were analyzed [20]; thus, vaginal swab samples confirmed by microscopy with T. vaginalis were placed in sterile plastic tubes with 0.9% physiological saline solution and kept at −70 °C until DNA isolation.

2.2. PCR Sequencing

DNA was extracted from approximately 500 µL of the physiological saline solution in which vaginal swabs were immersed using the Puregene DNA Purification System (Gentra Systems, Hilden, Germany) according to the manufacturer’s protocol. In samples with unclear microscopic observations, the presence of T. vaginalis was confirmed by analyzing the ITS1+5.8S+ITS2 region using primers and PCR conditions that were previously described [21].
To analyze the pfor gene, specific primers for variant A of different T. vaginalis sequences were designed at the laboratory based on information available in GenBank (sequences AY661465, U16822 and XM_001582310). A suitable region of ~852 bp was chosen for amplification with the forward primer 5′ TCAAGGTCCACCTCTTCCGCCC 3′ and the reverse primer 5′ GCATGGTCGACCATGTCCCAG 3′. This region includes the important sites of the enzyme, such as the domain II and III of the ferredoxin oxidoreductase site, which are all located between amino acids 305 and 603 of the protein. PCR amplifications were carried out in a final volume of 25 μL containing each primer at 10 pmol, 1X PCR buffer (8 mM Tris-HCl, pH 8; 20 mM KCl), 1.5 mM MgCl2, 0.5 mM dNTPs and 2 U of Taq DNA Polymerase (Promega, Madison, WI, USA). Approximately 500 ng of DNA was used as a template to amplify genomic sequences. The amplification conditions were 1 cycle at 94 °C for 5 min; 35 cycles including denaturation, annealing and extension steps at 94 °C for 30 s, 54 °C for 1 min and 72 °C for 30 s, respectively; and a final extension step at 72 °C for 7 min. The presence of amplicons was observed by staining with ethidium bromide followed by electrophoresis in a 1.5% agarose gel; in most samples, amplification assays were performed at least twice, both forward and reverse, to achieve the quality and concentration required for sequencing, then the amplicons were purified using an AxyPrep PCR clean-up kit (Axigen Biosciences, Union City, CA, USA) and sequenced on both strands in an ABI prism sequencer (Applied Biosystems, Foster City, CA, USA) using the Sanger method, and chromatograms were obtained.

2.3. Phylogenetic and Genetic Variation Analysis

The chromatograms for each sequence were evaluated with the Mesquite V2.5 software using the Chromaseq package [22] that partially makes use of the Phred and Phrap algorithms for base calling, assigning quality values to each called base and assembling contigs [23], because more than one sequence was obtained for most of the samples, and from them obtains a consensus. All sequences were subjected to a BLAST search of the GenBank database and submitted with the access numbers OR493214, OR069659-OR069683 for pfor A partial sequences and OR005496-OR005511 for the ITS1-5.8S rRNA-ITS2 region; multiple alignments were performed using CLUSTAL W and MUSCLE programs with manual adjustment in the MEGA program [24]. The molecular evolution model used to build the phylogenetic tree was the Hasegawa–Kishino–Yano model with a gamma distribution and invariant sites (HKY+G+I) for ITS sequences and the general time-reversible model with a gamma distribution and invariant sites (GTR+G+I) for the pfor A partial gene, which were selected using MODELTEST software version 3.7 [25]. The algorithms used to carry out the phylogenetic analyses were maximum likelihood (ML), maximum parsimony (MP) and the Bayesian algorithm (BA). ML and MP were performed with 1000 bootstrap replicates under the Tamura–Nei model of evolution, and the algorithms were implemented in MEGA7. BA reconstructions were performed using the MrBayes software version 3.2.7 [26], while the analysis was performed for two million generations with sampling trees every 100 generations including two independent assays with four strands each, considering codon splits. Trees with scores below the burn-in phase were discarded, and those remainders were collected and used to build majority consensus trees. A median-joining network analysis was performed using NETWORK 4.611 with sequences of T. vaginalis downloaded from GenBank, which showed high identity with the samples obtained in the present study according to the BLAST search. Haplotype networks were established under default settings and assumptions; sequences of the ITS1+5.8S+ITS2 region sequences of the ITS1+5.8S+ITS2 region for BA reconstructions were MK172845-7, KM095107; KX459478-79, 86; KX977518-20 and KC215389-90 from T. gallinae, T. tenax, T. brixi and T. stableri, respectively, using T. stableri as the outgroup, while for the pfor A tree, sequences of T. tenax and T. gallinae (OCTD01001480 and MRSU01008911, respectively) were included as the outgroup. In addition, some T. vaginalis sequences were downloaded from GenBank as controls.
A genetic diversity analysis among populations was performed using DnaSPv4 [27]; specifically, expected heterozygosity (Hd), nucleotide diversity (π), haplotype polymorphism (θ), Tajima’s D and Fu and Li’s D were calculated. The meaning of these keys are as follows: π, average number of nucleotide differences among all possible pairs of sequences in the sample; θ, proportion of nucleotide sites that can be predicted to be polymorphic from this region of the genome—both indices are used to measure polymorphisms ranging from 0 to 1. Finally, negative values for Tajima’s D and Fu and Li’s D tests allow us to test for past population expansion or neutrality of mutations; to this end, negative values suggest a recent expansion process or an effect of purifying selection [28,29].
To search for an association between the phylogenetic analysis of both markers (ITS and pfor A) and the main clinical characteristics of the samples, since the small sample size and frequency in at least one cell was less than 5, Fisher’s exact tests for two tails were performed, and odds ratios and 95% confidence intervals were also calculated using Epi-Info6 software v6.04 (Centers for Disease Control and Prevention, Atlanta, GA, USA).

2.4. Alignment of Protein Sequences and Prediction of Three-Dimensional Structure

The sequences obtained with the primers specific to the A variant of the pfor gene were translated to amino acid sequences and aligned with those for PFOR of T. vaginalis Q4KY23, Q27088 and A0A8U0WP41 deposited in the UniProt database. This analysis was performed using the ClustalW program. The I-TASSER [30] server was used to model the 3D structure of T. vaginalis PFOR variant A. This server is under active development with the intention of providing the most accurate predictions of protein structure and function using state-of-the-art algorithms. After analysis, the models with the highest confidence scores (C-scores) were selected. Three-dimensional structures were analyzed and visualized using PyMOL 2.5.2 [31].

3. Results

3.1. Characteristics of the Study Population

Twenty-eight vaginal swab samples of women with a mean age of 34 ± 13.7 years were analyzed; the main gynecological and microbiological findings for the women were previously reported [20]. Briefly, most of the patients described genital burning/itching (35.7%), 21.4% had cervicovaginitis and 17.8% exhibited vaginal fetid discharge and 42.8% were pregnant. Interestingly, in six samples (21.4%), coinfections with Blastocystis ST1-3 were found. DNA extractions and PCRs were performed for each sample from two to four times to reach the concentration and quality recommended for sequencing. Unfortunately, since the DNA obtained from the samples was used for a previous study [20], it was not possible to obtain amplicons of ITS1+5.8S+ITS2 and pfor of sufficient quality and concentration for sequencing in all cases. All sequences were subjected to a BLAST search against the GenBank database (OR493214, OR069659-OR069683 to pfor A partial sequences and OR005496-OR005511 for the ITS1-5.8S rRNA-ITS2 region).

3.2. Phylogenetic and Genetic Variation Findings

Sixteen samples sequenced for the ITS1+5.8S+ITS2 region showed two haplotypes; the median-joining network tree for this locus (Figure 1) included the haplotypes obtained here as a main cluster (subgroup 1), near the divergent cluster of T. tenax. Interestingly, a second cluster of T. vaginalis haplotypes was identified (subgroup 2), close to the T. gallinae cluster.
The phylogenetic tree grouped all samples analyzed with the T. vaginalis clade with strong support (97/99/1.00); although there were no strong support values, three subclades were also observed. In contrast, other trichomonads, such as T. gallinae, T. tenax, T. stableri and T. brixi, were clearly separated (Figure 2).
Sequences for a partial pfor A gene of 26 samples were obtained, and some haplotypes were identified; no correlations between the main clinical characteristics of the samples and the results of phylogenetic analysis for pfor A were found (Appendix A, Table A1). The median-joining network tree for the partial pfor A gene (Figure 3) showed a wide dispersion of haplotypes, without specific clustering. All pfor A haplotypes of T. vaginalis samples were included in a specific clade (100/100/1.00), separated from T. gallinae and T. tenax; although three small subclades were distinguished, the bootstrap support values were below 70 (Figure 4).
Genetic polymorphism indices of both studied genes were calculated, and Table 1 summarizes the results. The expected heterozygosity values for both markers (ITS or pfor A) ranged from 0.49 to 0.79; in contrast, for the π and θ indices, very low values (~0.002) were obtained. Regarding Tajima’s D and Fu and Li’s D, most of the values were negative but without statistical significance, except for those of the ITS1+5.8S+ITS2 haplotypes, which were statistically significant. Interestingly, when ITS subgroups were independently analyzed, subgroup 2 showed higher variability (π > 0.005) than subgroup 1; although Tajima’s D and Fu and Li’s D tests showed negative values for both subgroups, only subgroup 2 values had significant support.

3.3. Alignment of pfor A and Prediction of Its Three-Dimensional Structure

Regarding the protein deduction analysis, the nucleotide sequences of pfor A were translated, yielding sequences of approximately 284 amino acids. Protein alignments of such sequences showed four highly conserved cysteine residues. In addition, two mutations were found, one in the OR069675 sample of Val by Gly-137 and another in Q4KY23 of Ser by Pro-463 (Appendix B, Figure A1). The predicted three-dimensional structure of the protein (Figure 5A) shows the domain of pyruvate ferredoxin oxidoreductase in cyan color, consisting of six β-strand structures, exhibiting mutations in red and four α-helix structures. The consensus sequence Val and Ser (Figure 5B), and mutations to Gly and Pro (Figure 5C,D, respectively) can also be observed.

4. Discussion

Several factors have been documented that can modify the virulence of T. vaginalis, such as transcriptomics and proteomics of the adhesion process [32], the presence of the double-stranded RNA virus known as T. vaginalis virus (TVV) [33], interactions with human-associated bacteria such as Mycoplasma hominis [34] infection and the promotion of vaginal microbiota imbalance by increasing the proportions of Parvimonas, Sneathia and other anaerobes [35]; however, association among the symptomatology and the genetic variability of ITS and pfor A gene of T. vaginalis clinical isolates has not been documented.
PFOR is a metabolic enzyme described in several anaerobic microbial eukaryotes, such as Entamoeba histolytica [36], Giardia duodenalis [37] and T. vaginalis [9,38]. In trichomonads, it has been highlighted for possessing alternative, nonenzymatic functions [39]; thus, PFOR is also associated with the process of pathogenicity in T. vaginalis, such as adherence to host cells and proliferation [9,39]. In the present study, no correlations between the main clinical characteristics of the samples and the results of phylogenetic analysis for pfor A were found. As is often the case in many studies, phylogenetic inferences do not show clear agreement between the evolution of the pathogen and the development of the disease [40].
Although in the present study the ITS1+5.8S+ITS2 loci were originally used as a tool for the molecular identification of T. vaginalis, the sequences obtained and their subsequent analysis provided interesting information about the population genetic structure of this pathogen. The phylogenetic tree grouped our samples within the T. vaginalis main cluster; although subclustering for the main cluster were observed, the low bootstrap value (below 70) suggest that this subclustering should be interpreted with caution; in contrast, the median-joining network allowed us to sort all sequences into two subgroups: one with less variability and from which T. gallinae diverged and another undergoing an expansion event and that is close to the T. tenax cluster. On the other hand, an interesting global population genetics study using a panel of microsatellite markers and SNPs comprising isolates from different parts of the world, including 11 samples from Mexico, was performed [33]. In that study, it was shown that T. vaginalis is a genetically diverse parasite with a unique population structure, exhibiting two types (type 1 and type 2) present in equal proportions worldwide; therefore, our results for the ITS1+5.8S+ITS2 loci are consistent with those presented by Conrad et al. [33].
For pfor A sequences of T. vaginalis, a robust phylogenetic tree was generated that clustered all sequences in a clear and specific cluster, away from the T. gallinae and T. tenax clusters; thus, our results suggest that this gene could eventually be used as a marker to confirm genus and species identities among trichomonads. However, because there are few sequences available for this gene in the Trichomonas genus, further studies are required to confirm this finding. When two different loci are submitted to a phylogenetic analysis, it is common to observe dissimilarities through inference [40]. A study focused on the genetic variability of pfor and the small-subunit ribosomal RNA (18S rRNA) from Blastocystis sp. also found similarities between the phylogenies of the two markers [41]. In another study, during the sequencing analysis of the Fe-hydrogenase gene of T. gallinae isolates, the median joining network analysis showed a distribution of haplotypes similar to the median joining network inference for pfor halplotypes of T. vaginalis obtained in the present study; both Fe-hydrogenase and PFOR are enzymes placed in the hydrogenosomes of Trichomonas spp. [42]
Regarding the genetic variability of ITS and pfor A, both sequences showed very low values (~0.002) for the π and θ indices. In addition, Tajima’s D and Fu and Li’s D tests showed negative values with statistical significance only for ITS1+5.8S+ITS2 haplotypes. A study on the identification of genetic variants of T. vaginalis in China, analyzing fragments of the 18S rRNA, found that π values ranged from 0.002 to 0.005 among different isolates, with a significant negative value for Fu’s test [29]. Thus, our results are in concordance with those of Mao and Liu [43] and with those of Conrad et al. [33] because they support the finding that, as noted above, T. vaginalis is a genetically diverse parasite with a unique population structure, exhibiting two types (types 1 and 2), while the Mexican samples were grouped exclusively into a single type, which could be undergoing a process of expansion.
On the other hand, since the diversity and phylogeny of pfor A may allow us to infer the function of the protein or may be related to characteristics of the population analyzed (i.e., in Blastocystis sp., pfor did not allow discrimination between subtypes of parasites), the phylogeny of this enzyme may provide information about the biochemical function instead of the relationship of the group [41].
The alignment of the amino acid sequences of the samples obtained in this study with the PFOR sequences of T. vaginalis Q4KY23, Q27088 and A0A8U0WP41 showed four highly conserved cysteines; interestingly, these cysteines are far from each other, i.e., CX90-CX60-CX87-CX, and there is no similarity to any cysteine-rich motif reported in the literature. Since our results showed that Gly and Pro mutations in the PFOR sequence are not located near the cysteine-rich sites mentioned above, it is likely that these mutations have no impact on the PFOR activity of T. vaginalis; in addition, some mutations in nitroreductases (enzymes involved in the processing of metronidazole) of T. vaginalis have been identified, e.g., mutations from Val to Asp-51, Val to Ala-153, Val to Ala-159, Val to Ala-117, Val to Ala-132 and Ser to Gly-29, and only the Val mutation to Ala-59 has been associated with resistance to metronidazole [44]. Thus, with the design of the present study, it is not possible to clarify the role of the PFOR mutations.
Until now, the heterogeneity of clinical manifestations of trichomoniasis associated with differences in phenotypic and genetic expression of T. vaginalis has not been clear [16,17,43]. The population or isolates of T. vaginalis show significantly different levels of genetic variation; thus, the proper characterization of these populations is vital to avoid generalizing the behavior of each trichomonad, as observed for the Mexican population. Understanding these relationships is essential for future strategies to control or prevent T. vaginalis in each region.
Since the present study was based on sequence analysis of a fragment of the pfor A gene and ITS marker from clinical isolates of T. vaginalis recovered from 26 patients, the small sample size constitutes a potential bias in the analysis of their genetic variability; however, in both markers it was possible to yield robust phylogenetic inferences that allows identification among different clusters and haplotypes as well as a probable population expansion process in this parasite. In addition, the development of future studies using next-generation sequencing (NGS) for a more extensive search and accurate results applied to a case-control design would help to clarify the etiological role of these markers in clinical isolates of T. vaginalis and include the study of samples from other populations and male patients.

5. Conclusions

Our results suggest that the role of genetic variability of pfor and ITS may not be significant in the symptomatology of this pathogen; however, their utility as genus and species markers in trichomonads is promising.

Author Contributions

Conceptualization, F.M.-H., F.S.-A. and G.V.; methodology, F.M.-H., F.S.-A., J.M.-O., N.R.G.-A. and E.L.-E.; software, F.M.-H.; J.M.-O. and N.R.G.-A.; validation, F.M.-H., M.R.-V. and P.M.; formal analysis, F.M.-H., M.R.-V., P.M. and G.V.; investigation, F.M.-H., F.S.-A., E.L.-E., P.M. and G.V.; resources, F.M.-H. and P.M.; data curation, F.S.-A., E.L.-E. and N.R.G.-A.; writing—original draft preparation, F.M.-H., P.M. and G.V.; supervision, F.M.-H., M.R.-V., P.M. and G.V.; funding acquisition, F.M.-H., M.R.-V. and P.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Secretaria de Salud, Mexico, Hospital General “Dr. Manuel Gea Gonzalez”, Grant #E022.

Institutional Review Board Statement

The Research Ethics and Research Committees of the General Hospital “Dr. Manuel Gea Gonzalez” (HGMGG) approved the study with reference number 12-75-2015 (see reference [22]).

Informed Consent Statement

Signed, informed consent was obtained from all participants (see reference [22]).

Data Availability Statement

All relevant data are within the article. The sequences data were submitted to the GenBank database under the accession numbers OR005496-OR005511 for the ITS1-5.8S rRNA-ITS2 region and OR069659-OR069683 to pfor locus.

Acknowledgments

Laura Marquez and Nelly Lopez, LaNaBio, Instituto de Biologia UNAM, are acknowledged for their assistance during the sequencing technique.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Association analysis between main clinical status of patients with T. vaginalis and subclustering for pfor A.
Table A1. Association analysis between main clinical status of patients with T. vaginalis and subclustering for pfor A.
SymptomsNumber of PatientsGenBank Accession Numberspfor A Subclustering
(Clusters according to Figure 4)
A n = 11p Value * A vs. B+C+DB n = 3p Value * B vs. A+C+DC n = 4p Value * C vs. A+B+DD n = 8p Value * D vs. A+B+C
Genital burning/itching9OR069659 OR069661 OR069663 OR069665 OR069667 OR069670 OR069673 OR069679 OR06966420.21620.26730.10320.675
Cervicovaginitis6OR069661 OR069663 OR069673 OR069664 OR069682 OR06967021.00001.00030.027 **10.627
Vaginal discharge with unpleasant smell5OR069659 OR069662 OR069671 OR069681 OR06968321.00020.08400.55511.000
Myomatosis1OR06966710.42301.00001.00001.000
Dysuria1OR06966801.00001.00010.15301.000
Dyspareunia1OR06967301.00001.00010.15301.000
Gestational diabetes3OR069663 OR069673 OR06967800.23001.00020.56511.000
Overweight/obesity15OR069659 OR069662 OR069667 OR069668 OR069669 OR069670 OR069673 OR069674 OR069676 OR069680 OR069677 OR069678 OR069679 OR069681 OR06966461. 00030.23821.00030.2183
Pregnancy15OR069661 OR069660 OR069663 OR069665 OR069666 OR069671 OR069673 OR069674 OR069676 OR069675 OR069677 OR069678 OR069664 OR069682 OR06968361.00030.23821.00051.000
* p value by Fisher’s exact test; ** for this item, odds ratio (95% confidence interval) = 19 (1.07–709.5), due to the wide range of the confidence interval, its significance must be interpreted with care.

Appendix B

Microorganisms 11 02240 g0a1
Figure A1. Alignment of the sequences obtained in this study (translated into amino acids) with the PFOR sequences of T. vaginalis Q4KY23, Q27088 and A0A8U0WP41 deposited in the UniProt database. Only the segment of the alignment corresponding to amplification with the primers of this study is included. The arrows in the opposite directions indicate the forward and reverse primers. Asterisks represent identities. Highly preserved cysteine residues are indicated by a black triangle. Gly and Pro mutations are indicated by gray shading. Alignment was obtained with the multisequence alignment program CLUSTAL W (version 1.5).
Figure A1. Alignment of the sequences obtained in this study (translated into amino acids) with the PFOR sequences of T. vaginalis Q4KY23, Q27088 and A0A8U0WP41 deposited in the UniProt database. Only the segment of the alignment corresponding to amplification with the primers of this study is included. The arrows in the opposite directions indicate the forward and reverse primers. Asterisks represent identities. Highly preserved cysteine residues are indicated by a black triangle. Gly and Pro mutations are indicated by gray shading. Alignment was obtained with the multisequence alignment program CLUSTAL W (version 1.5).
Microorganisms 11 02240 g0a1

References

  1. Rowley, J.; Vander, H.S.; Korenromp, E.; Low, N.; Unemo, M.; Abu-Raddad, L.J.; Chico, R.M.; Smolak, A.; Newman, L.; Gottlieb, S.; et al. Chlamydia, gonorrhoea, trichomoniasis and syphilis: Global prevalence and incidence estimates, 2016. Bull. World Health Organ. 2019, 97, 548–562P. [Google Scholar] [CrossRef]
  2. Schwebke, J.R.; Burguess, D. Trichomoniasis. Clin. Microbiol. Rev. 2004, 17, 794–803. [Google Scholar] [CrossRef]
  3. Kissinger, P. Trichomonas vaginalis: A review of epidemiologic, clinical and treatment issues. BMC Infect. Dis. 2015, 15, 307. [Google Scholar] [CrossRef]
  4. Petrin, D.; Delgaty, K.; Bhatt, R.; Garber, G. Clinical and microbiological aspects of Trichomonas vaginalis. Clin. Microbiol. Rev. 1998, 11, 300–317. [Google Scholar] [CrossRef]
  5. Van der Giezen, M.; Tovar, J.; Clark, C.G. Mitochondrion-derived organelles in protest and fungi. Int. Rev. Cytol. 2005, 244, 175–225. [Google Scholar]
  6. Carlton, J.M.; Hirt, R.P.; Silva, J.C.; Delcher, A.L.; Schatz, M.; Zhao, Q.; Wortman, J.R.; Bidwell, S.L.; Alsmark, U.C.; Besteiro, S.; et al. Draft genome sequence of the sexually transmitted pathogen Trichomonas vaginalis. Science 2007, 315, 207–212. [Google Scholar] [CrossRef]
  7. Jeffery, C.J. Mass spectrometry and the search for moonlighting proteins. Mass Spectrom. Rev. 2005, 24, 772–782. [Google Scholar] [CrossRef]
  8. Hug, L.A.; Stechmann, A.; Roger, A. Phylogenetic distribution and histories of proteins involved in anaerobic pyruvate metabolism in eukaryotes. Mol. Biol. Evol. 2009, 27, 311–324. [Google Scholar] [CrossRef]
  9. Meza-Cervantes, P.; Gonzalez-Robles, A.; Cardenas-Guerra, R.E.; Ortega-Lopez, J.; Saavedra, E.; Pineda, E.; Arroyo, R. Pyruvate: Ferredoxin oxidoreductase (PFO) is a surface-associated cell-binding protein in Trichomonas vaginalis and is involved in trichomonal adherence to host cell. Microbiology 2011, 157, 3469–3482. [Google Scholar] [CrossRef]
  10. Hrdy, I.; Müller, M. Primary structure and eubacterial relationship of the pyruvate: Ferredoxin oxidoreductase of the amitochondriate eukaryote Trichomonas vaginalis. J. Mol. Evol. 1995, 41, 388–396. [Google Scholar] [CrossRef]
  11. Alderete, J.F.; Garza, G.E. Specific nature of Trichomonas vaginalis parasitism of host cell surfaces. Infect. Immun. 1985, 50, 701–707. [Google Scholar] [CrossRef]
  12. Song, H.O. Influence of 120 kDa Pyruvate: Ferredoxin Oxidoreductase on Pathogenicity of Trichomonas vaginalis. Korean J. Parasitol. 2016, 54, 71–74. [Google Scholar] [CrossRef]
  13. Moreno-Brito, V.; Yañez-Gomez, C.; Meza-Cervantes, P.; Avila-Gonzalez, L.; Rodriguez, M.A.; Ortega-Lopez, J.; Gonzalez-Robles, A.; Arroyo, R. A Trichomonas vaginalis 120 KDa protein with identity to hydrogenosome pyruvate: Ferrodoxin oxidoreductase is a surface adhesin induce by iron. Cell Microbiol. 2005, 7, 245–258. [Google Scholar] [CrossRef]
  14. Upcroft, J.A.; Delgadillo-Correa, M.G.; Dunne, R.L.; Sturm, A.W.; Johnson, P.J.; Upcroft, P. Genotyping Trichomonas vaginalis. Int. J. Parasitol. 2006, 36, 821–828. [Google Scholar] [CrossRef]
  15. Villalobos, G.; Orozco-Mosqueda, G.E.; Lopez-Perez, M.; Lopez-Escamilla, E.; Córdoba-Aguilar, A.; Rangel-Gamboa, L.; Olivo-Diaz, A.; Romero-Valdovinos, M.; Maravilla, P.; Martinez-Hernandez, F. Suitability of internal transcribed spacers (ITS) as markers for the population genetic structure of Blastocystis spp. Parasit. Vectors 2014, 7, 461. [Google Scholar] [CrossRef]
  16. Snipes, L.J.; Gamard, P.M.; Narcisi, E.M.; Beard, C.B.; Lehmann, T.; Secor, W.E. Molecular epidemiology of metronidazole resistance in a population of Trichomonas vaginalis clinical isolates. J. Clin. Microbiol. 2000, 38, 3004–3009. [Google Scholar] [CrossRef]
  17. Rivera, W.L.; Ong, V.A.; Masalunga, M.C. Molecular characterization of Trichomonas vaginalis isolates from the Philippines. Parasitol. Res. 2009, 106, 105–110. [Google Scholar] [CrossRef]
  18. Ibáñez-Escribano, A.; Nogal-Ruiz, J.J.; Arán, V.J.; Escario, J.A.; Gómez-Barrio, A.; Alderete, J.F. Determination of internal transcribed spacer regions (ITS) in Trichomonas vaginalis isolates and differentiation among Trichomonas species. Parasitol. Int. 2014, 63, 427–431. [Google Scholar] [CrossRef]
  19. Ertabaklar, H.; Ertuğ, S.; Çalışkan, S.Ö.; Malatyalı, E.; Bozdoğan, B. Use of internal transcribed spacer sequence polymorphisms as a method for Trichomonas vaginalis genotyping. Turkiye Parazitol. Derg. 2018, 42, 6–10. [Google Scholar] [CrossRef]
  20. Villalobos, G.; Sanchez-Aguillon, F.; Carmona-Maldonado, M.V.; Gonzalez-Arenas, N.R.; Lopez-Escamilla, E.; Hernandez-Castro, R.; Romero-Valdovinos, M.; Martinez-Flores, W.A.; Ramirez-Hinojosa, J.P.; Maravilla, P.; et al. Unexpected presence of Blastocystis subtype 1-3 DNA in human vaginal and sperm samples coinfected with Trichomonas vaginalis. Korean J. Parasitol. 2022, 60, 195–200. [Google Scholar] [CrossRef]
  21. Lopez-Escamilla, E.; Sanchez-Aguillon, F.; Alatorre-Fernandez, C.P.; Aguilar-Zapata, D.; Arroyo-Escalante, S.; Arellano, T.; Moncada-Barron, D.; Romero-Valdovinos, M.; Martinez-Hernandez, F.; Rodriguez-Zulueta, P.; et al. New Tetratrichomonas species in two patients with pleural empyema. J. Clin. Microbiol. 2013, 51, 3143–3146. [Google Scholar] [CrossRef]
  22. Maddison, W.P.; Maddison, D.R. Mesquite: A Modular System for Evolutionary Analysis. (Version 3.81). Available online: http://www.mesquiteproject.org/ (accessed on 30 July 2023).
  23. Ewing, B.; Hillier, L.; Wendl, M.C.; Green, P. Base-calling of automated sequencer traces using Phred. I. Accuracy assessment. Genome Res. 1998, 8, 175–185. [Google Scholar] [CrossRef]
  24. Tamura, K.; Peterson, D.; Peterson, N.; Stecher, G.; Nei, M.; Kumar, S. MEGA5: Molecular Evolutionary Genetics Analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol. Biol. Evol. 2011, 28, 2731–2739. [Google Scholar] [CrossRef]
  25. Posada, D.; Crandall, K.A. Modeltest: Testing the model of DNA substitution. Bioinformatics 1998, 14, 817–818. [Google Scholar] [CrossRef]
  26. Ronquis, F.; Huelsenbeck, J.P. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 2003, 19, 1572–1574. [Google Scholar] [CrossRef]
  27. Rozas, J.; Sanchez-DelBarrio, J.C.; Messeguer, X.; Rozas, R. DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics 2003, 19, 2496–2497. [Google Scholar] [CrossRef]
  28. Hedrick, P.W. Measures of Genetics Variation. In Genetics of Populations, 2nd ed.; Jones and Bartlett Publishers: Sudbury, MA, USA, 1999; pp. 47–87. [Google Scholar]
  29. Fu, Y.-X.; Li, W.-H. Statistical test of neutrality of mutations. Genetics 1993, 133, 693–709. [Google Scholar] [CrossRef]
  30. Yang, J.; Yan, R.; Roy, A.; Xu, D.; Poisson, J.; Zhang, Y. The I-TASSER Suite: Protein structure and function prediction. Nat. Methods 2015, 12, 7–8. [Google Scholar] [CrossRef]
  31. Schrodinger, L.L.C. The PyMOL Molecular Graphics System, Version 1.8; Schrodinger, LLC: New York, NY, USA, 2015. [Google Scholar]
  32. Hirt, R.P. Trichomonas vaginalis virulence factors: An integrative overview. Sex. Transm. Infect. 2013, 89, 439–443. [Google Scholar] [CrossRef]
  33. Conrad, M.D.; Gorman, A.W.; Schillinger, J.A.; Fiori, P.R.; Arroyo, R.; Malla, N.; Dubey, M.L.; Gonzalez, J.; Blank, S.; Secor, W.E.; et al. Extensive genetic diversity, unique population structure and evidence of genetic exchange in the sexually transmitted parasite Trichomonas vaginalis. PLoS Neglected Trop. Dis. 2012, 6, e1573. [Google Scholar] [CrossRef]
  34. Dessì, D.; Rappelli, P.; Diaz, N.; Cappuccinelli, P.; Fiori, P.L. Mycoplasma hominis and Trichomonas vaginalis: A unique case of symbiotic relationship between two obligate human parasites. Front. Biosci. 2006, 11, 2028–2034. [Google Scholar]
  35. Brotman, R.M.; Bradford, L.L.; Conrad, M.; Gajer, P.; Ault, K.; Peralta, L.; Forney, L.J.; Carlton, J.N.; Abdo, Z.; Ravel, J. Association between Trichomonas vaginalis and vaginal bacterial community composition among reproductive-age women. Sex. Transm. Dis. 2012, 39, 807–812. [Google Scholar] [CrossRef]
  36. Pineda, E.; Encalada, R.; Rodríguez-Zavala, J.S.; Olivos-García, A.; Moreno-Sánchez, R.; Saavedra, E. Pyruvate: Ferredoxin oxidoreductase and bifunctional aldehyde-alcohol dehydrogenase are essential for energy metabolism under oxidative stress in Entamoeba histolytica. FEBS J. 2010, 277, 3382–3395. [Google Scholar] [CrossRef]
  37. Townson, S.M.; Upcroft, J.A.; Upcroft, P. Characterization and purification of pyruvate: Ferredoxin oxidoreductase from Giardia duodenalis. Mol. Biochem. Parasitol. 1996, 79, 183–193. [Google Scholar] [CrossRef]
  38. Steinbüchel, A.; Müller, M. Anaerobic pyruvate metabolism of Tritrichomonas foetus and Trichomonas vaginalis hydrogenosomes. Mol. Biochem. Parasitol. 1986, 20, 57–65. [Google Scholar] [CrossRef]
  39. Alderete, J.F.; Millsap, K.W.; Lehker, M.W.; Benchimol, M. Enzymes of microbial pathogens and Trichomonas vaginalis: Molecular mimicry and functional diversity. Cell. Microbiol. 2001, 3, 359–370. [Google Scholar] [CrossRef]
  40. Baum, D.A.; Smith, S.D. Tree Thinking: An Introduction to Phylogenetic Biology; Roberts and Company Publishers: Greenwood Village, CO, USA, 2013; pp. 56–78. [Google Scholar]
  41. Alarcon-Valdes, P.; Villalobos, G.; Martinez-Flores, W.A.; Lopez-Escamilla, E.; Gonzalez-Arenas, N.R.; Romero-Valdovinos, M.; Martinez-Hernandez, F.; Santillan-Benitez, J.G.; Maravilla, P. Can the pyruvate: Ferredoxin oxidoreductase (PFOR) gene be used as an additional marker to discriminate among Blastocystis strains or subtypes? Parasit. Vectors 2018, 11, 564. [Google Scholar] [CrossRef]
  42. Alrefaei, A.F.; Albeshr, M.F.; Alharbi, S.N.; Alrefaei, A.F.; Almutairi, M.H.; Almutairi, B.O.; Nader, J.L.; Manoharadas, S. Molecular characterization of the Fe-hydrogenase gene marker in Trichomonas gallinae isolated from birds in Riyadh, Saudi Arabia. Parasitol. Int. 2021, 8, 102263. [Google Scholar] [CrossRef]
  43. Mao, M.; Liu, H.L. Genetic diversity of Trichomonas vaginalis clinical isolates from Henan province in central China. Pathog. Glob. Health 2015, 109, 242–246. [Google Scholar] [CrossRef]
  44. Paulish-Miller, T.E.; Augostini, P.; Schuyler, J.A.; Smith, W.L.; Mordechai, E.; Adelson, M.E.; Gygax, S.E.; Secor, W.E.; Hilbert, D.W. Trichomonas vaginalis metronidazole resistance is associated with single nucleotide polymorphisms in the nitroreductase genes ntr4Tv and ntr6Tv. Antimicrob. Agents Chemother. 2014, 58, 2938–2943. [Google Scholar] [CrossRef]
Microorganisms 11 02240 g001
Figure 1. Haplotype network tree generated using ITS1+5.8S+ITS2 region sequences of T. vaginalis from Mexico and GenBank; the numbers on branches refer to mutational changes. Haplotypes for T. gallinae are shown in blue, for T. tenax in red, for Mexican samples obtained here in black and downloaded from GenBank in yellow. Two subgroups (1 and 2, respectively) can be differentiated, one close to the species T. gallinae and another to T. tenax; both subgroups are shown by two large circles. Numbers in branches are mutational changes; sizes of circles are proportional to haplotype frequencies.
Figure 1. Haplotype network tree generated using ITS1+5.8S+ITS2 region sequences of T. vaginalis from Mexico and GenBank; the numbers on branches refer to mutational changes. Haplotypes for T. gallinae are shown in blue, for T. tenax in red, for Mexican samples obtained here in black and downloaded from GenBank in yellow. Two subgroups (1 and 2, respectively) can be differentiated, one close to the species T. gallinae and another to T. tenax; both subgroups are shown by two large circles. Numbers in branches are mutational changes; sizes of circles are proportional to haplotype frequencies.
Microorganisms 11 02240 g001
Microorganisms 11 02240 g002
Figure 2. Phylogenetic tree based on ITS1+5.8S+ITS2 region sequences of T. vaginalis from Mexico and GenBank. The numbers at the nodes indicate bootstrap support values and Bayesian posterior probabilities, respectively, given by ML, MP and Bayesian analyses. The sequences obtained in this study are shown in bold. Bar = estimated number of substitutions per site.
Figure 2. Phylogenetic tree based on ITS1+5.8S+ITS2 region sequences of T. vaginalis from Mexico and GenBank. The numbers at the nodes indicate bootstrap support values and Bayesian posterior probabilities, respectively, given by ML, MP and Bayesian analyses. The sequences obtained in this study are shown in bold. Bar = estimated number of substitutions per site.
Microorganisms 11 02240 g002
Microorganisms 11 02240 g003
Figure 3. Haplotype network tree generated using pfor A sequences of T. vaginalis from Mexico and GenBank; the numbers on branches refer to mutational changes. Haplotypes for T. gallinae are shown in blue and T. tenax in red. Those of Mexican samples obtained here are shown in black, and haplotypes downloaded from GenBank are shown in yellow. Numbers in branches are mutational changes; sizes of circles are proportional to haplotype frequencies.
Figure 3. Haplotype network tree generated using pfor A sequences of T. vaginalis from Mexico and GenBank; the numbers on branches refer to mutational changes. Haplotypes for T. gallinae are shown in blue and T. tenax in red. Those of Mexican samples obtained here are shown in black, and haplotypes downloaded from GenBank are shown in yellow. Numbers in branches are mutational changes; sizes of circles are proportional to haplotype frequencies.
Microorganisms 11 02240 g003
Microorganisms 11 02240 g004
Figure 4. Phylogenetic tree based on pfor A sequences of T. vaginalis from Mexico and GenBank. The numbers at the nodes indicate bootstrap support values and Bayesian posterior probabilities, respectively, given by ML, MP and Bayesian analyses. The sequences obtained in this study are shown in bold. Bar = estimated number of substitutions per site.
Figure 4. Phylogenetic tree based on pfor A sequences of T. vaginalis from Mexico and GenBank. The numbers at the nodes indicate bootstrap support values and Bayesian posterior probabilities, respectively, given by ML, MP and Bayesian analyses. The sequences obtained in this study are shown in bold. Bar = estimated number of substitutions per site.
Microorganisms 11 02240 g004
Microorganisms 11 02240 g005
Figure 5. Prediction of three-dimensional structure for PFOR A of T. vaginalis. In (A), the whole three-dimensional structure is shown. The image shows the domain of pyruvate ferredoxin oxide reductase in cyan. In addition, mutated amino acids are shown in red. Panel (B) provides a zoomed-in view of the consensus sequence with Val and Ser. In (C), the mutation of Val to Gly is observed, and, in (D), the mutation of Ser to Pro is shown.
Figure 5. Prediction of three-dimensional structure for PFOR A of T. vaginalis. In (A), the whole three-dimensional structure is shown. The image shows the domain of pyruvate ferredoxin oxide reductase in cyan. In addition, mutated amino acids are shown in red. Panel (B) provides a zoomed-in view of the consensus sequence with Val and Ser. In (C), the mutation of Val to Gly is observed, and, in (D), the mutation of Ser to Pro is shown.
Microorganisms 11 02240 g005
Table 1. Genetic polymorphism indices between different Trichomonas vaginalis sequences.
Table 1. Genetic polymorphism indices between different Trichomonas vaginalis sequences.
PopulationnExpected Heterozygosity (Hd)Nucleotide Diversity (π)Haplotype Polymorphism (θ)Tajima’s D
(p)		Fu and Li’s D
(p)
Whole pfor A *320.7900.00190.0033−1.267
(p > 0.10)		−1.540
(p > 0.10)
Mexican samples pfor A260.8150.00210.0031−1.041
(p > 0.10)		−1.069
(p > 0.10)
Whole ITS ‡680.4930.00270.0087−1.942
(p < 0.05)		−3.253
(p < 0.05)
Mexican samples ITS160.3250.00110.00100.155
(p > 0.10)		0.688
(p > 0.10)
ITS type 1280.6480.00420.0087−1.723
(p > 0.10)		−1.997
(p > 0.10)
ITS type 2 +420.4870.00580.0022−2.443
(p < 0.05)		−2.708
(p < 0.05)
* Five sequences of pfor A gene were downloaded from GenBank, and they were added as references. ‡ Fifty-two sequences of the ITS1+5.8S+ITS2 region were downloaded from GenBank, and they were added as references. + Twenty-six sequences of the ITS1+5.8S+ITS2 region were downloaded from GenBank, and they were added as references. Characters in bold indicate statistically significant values
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Martinez-Hernandez, F.; Sanchez-Aguillon, F.; Martinez-Ocaña, J.; Gonzalez-Arenas, N.R.; Romero-Valdovinos, M.; Lopez-Escamilla, E.; Maravilla, P.; Villalobos, G. Genetic Variability of the Internal Transcribed Spacer and Pyruvate:Ferredoxin Oxidoreductase Partial Gene of Trichomonas vaginalis from Female Patients. Microorganisms 2023, 11, 2240. https://doi.org/10.3390/microorganisms11092240

AMA Style

Martinez-Hernandez F, Sanchez-Aguillon F, Martinez-Ocaña J, Gonzalez-Arenas NR, Romero-Valdovinos M, Lopez-Escamilla E, Maravilla P, Villalobos G. Genetic Variability of the Internal Transcribed Spacer and Pyruvate:Ferredoxin Oxidoreductase Partial Gene of Trichomonas vaginalis from Female Patients. Microorganisms. 2023; 11(9):2240. https://doi.org/10.3390/microorganisms11092240

Chicago/Turabian Style

Martinez-Hernandez, Fernando, Fabiola Sanchez-Aguillon, Joel Martinez-Ocaña, Nelly Raquel Gonzalez-Arenas, Mirza Romero-Valdovinos, Eduardo Lopez-Escamilla, Pablo Maravilla, and Guiehdani Villalobos. 2023. "Genetic Variability of the Internal Transcribed Spacer and Pyruvate:Ferredoxin Oxidoreductase Partial Gene of Trichomonas vaginalis from Female Patients" Microorganisms 11, no. 9: 2240. https://doi.org/10.3390/microorganisms11092240

APA Style

Martinez-Hernandez, F., Sanchez-Aguillon, F., Martinez-Ocaña, J., Gonzalez-Arenas, N. R., Romero-Valdovinos, M., Lopez-Escamilla, E., Maravilla, P., & Villalobos, G. (2023). Genetic Variability of the Internal Transcribed Spacer and Pyruvate:Ferredoxin Oxidoreductase Partial Gene of Trichomonas vaginalis from Female Patients. Microorganisms, 11(9), 2240. https://doi.org/10.3390/microorganisms11092240

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop