Identification and Sequencing of Candida krusei Aconitate Hydratase Gene Using Rapid Amplification of cDNA Ends Method and Phylogenetic Analysis

authors:

avatar Roohollah Fateh 1 , 2 , avatar Farideh Zaini 2 , * , avatar Parivash Kordbacheh 2 , avatar Mehraban Falahati 3 , avatar Sasan Rezaie 2 , avatar Roshanak Daie Ghazvini 2 , avatar Nahid Borhani 4 , avatar Mahin Safara 2 , avatar Azam Fattahi 2 , avatar Ali Kanani 2 , avatar Shirin Farahyar 3 , avatar Manzar Bolhassani 5 , avatar Mansour Heidari 5 , 6

Department of Microbiology and Immunology, Faculty of Medicine, Qom University of Medical Sciences, Qom, IR Iran
Department of Medical Mycology and Parasitology, School of Public Health, Tehran University of Medical Sciences, Tehran, IR Iran
Department of Medical Mycology and Parasitology, Faculty of Medicine, Iran University of Medical Sciences, Tehran, IR Iran
Department of Parasitology, Pasteur Institute of Iran, Tehran, IR Iran
Department of Medical Genetics, Tehran University of Medical Sciences, Tehran, IR Iran
Experimental Research Center, Tehran University of Medical Sciences, Tehran, IR Iran

how to cite: Fateh R, Zaini F, Kordbacheh P, Falahati M, Rezaie S, et al. Identification and Sequencing of Candida krusei Aconitate Hydratase Gene Using Rapid Amplification of cDNA Ends Method and Phylogenetic Analysis. Jundishapur J Microbiol. 2015;8(11):e25218. https://doi.org/10.5812/jjm.25218.

Abstract

Background:

The production and development of an effective fungicidal drug requires the identification of an essential fungal protein as a drug target. Aconitase (ACO) is a mitochondrial protein that plays a vital role in tricarboxylic acid (TCA) cycle and thus production of energy within the cell.

Objectives:

The current study aimed to sequence Candida krusei ACO gene and determine any amino acid residue differences between human and fungal aconitases to obtain selective inhibition.

Materials and Methods:

Candida krusei (ATCC: 6258) aconitase gene was determined by 5’Rapid Amplification of cDNA Ends (RACE) method and degenerate Polymerase Chain Reaction (PCR) and analyzed using bioinformatics softwares.

Results:

One thousand-four hundred-nineteen nucleotide of C. krusei aconitase gene were clarified and submitted in Genbank as a partial sequence and then taxonomic location of C. krusei was determined by nucleotide and amino acid sequences of this gene. The comparison of nucleotide and amino acid sequences of Candida species ACO genes showed that C. krusei possessed characteristic sequences. No significant differences were observed between C. krusei and human aconitases within the active site amino acid residues.

Conclusions:

Results of the current study indicated that aconitase was not a suitable target to design new anti-fungal drugs that selectively block this enzyme.

1. Background

Candida species are the most common cause of opportunistic fungal infections and Candida albicans is the most predominant fungal opportunistic pathogen (1). However the relative prevalence and incidence of the Candida species depends on the geographical location, patient population, and clinical settings (2-4). In recent years, an increase is observed in invasive mycosis, particularly Candida infections (5-9). The increase is accompanied by a shift towards a higher proportion of species other than C. albicans, which are frequently resistant to fluconazole. C. krusei is one of the species that inherited resistance to fluconazole and is mainly isolated from immunocompromised patients such as patients with hematologic malignancies (10, 11).

Widespread use of fluconazole to prevent fungal infections in patients with human immunodeficiency virus infection contributes to a significant increase in C. krusei infection, especially due to the high incidence of resistance to this drug (12). Many of classes of anti-fungal drugs are currently available to treat fungal infections. These classes include polyene, azoles, allymines and echinocandins (13). The azole compounds, especially fluconazole are the most frequently used anti-fungal agents (14). Although the azole drugs have good pharmacokinetic properties and low toxicity, they are fungistatic and inhibit only the growth of the fungal cells and do not kill them. A critical problem is the resistance of fungal pathogens to antifungal drugs especially azole drugs. Currently, C. krusei is considered as a resistant to the azole drugs (15). The mortality rate due to fungal infections remains very high because of the limitations of effective available drugs (16).

The production and development of an effective fungicidal drug requires the identification of an essential fungal protein as drug target (17). Aconitate hydratase (ACO) is a mitochondrial protein that plays a vital role in tricarboxylic acid (TCA) cycle and thus production of energy within the cell (18). There are two forms of the protein aconitase (aconitate hydratase) in eukaryotes such as fungi: a mitochondrial and a cytoplasmic form. Mitochondrial aconitase convert citrate to isocitrate via a cis-aconitate intermediate in the TCA cycle. The cytoplasmic aconitase protein, in addition to enzymatic activity, is also an iron regulatory protein (IRP) (19). Targeting the fungal aconitase selectively would lead to the rapid death of fungal cells with little or no toxicity to the human host. It is necessary to obtain selective inhibition to determine any amino acid residue differences between human and fungal aconitases. These amino acids could then be used as a target to selectively bind and inhibit the function of fungal aconitase, thus produce the basis to innovate a new antifungal drug class.

Several studies commonly indicated that C. krusei virulence was less than that of C. albicans in terms of its adherence to epithelial cells and prosthetic surfaces, proteolytic activity and production of phospholipases. Furthermore, it seems that structural and metabolic features of C. krusei are significantly different from those of other medically common Candida spp. and display different behavioral patterns to host defenses, adding credence to the belief that it should be reassigned taxonomically (12, 20, 21). rapid amplification of cDNA ends (RACE) is a method to amplify cDNA fragment from an mRNA template with a known internal site and unknown sequences at either the 3’ or 5’-end of the mRNA (22). This procedure is defined by other researchers as one-sided PCR (23) or anchored PCR (24).

2. Objectives

The current study employed 5’ RACE method and degenerate PCR to partial lysequence the C. krusei aconitase gene. The current study aimed to identify C. krusei taxonomic location based on ACO gene sequence, determine any structural differences between C. krusei, common pathogenic Candida spp. and human aconitase sequences, and produce the basis to create a new antifungal drug to block the C. krusei aconitase function.

3. Materials and Methods

3.1. Strains

The current study used the standard strain of C. krusei (ATCC: 6258). The yeast was cultured on yeast extract (5g/L, Baltimore Biological Laboratory, USA)/ peptone (10g/L, Merck, Germany)/dextrose (20g/L, Merck, Germany) agar (20g/L, Merck, Germany) (YEPD) plate, incubated at 37°C for 48 hours, and then transferred into yeast extract peptone dextrose broth at 37°C for 24 hours.

3.2. 5’ Rapid Amplification of cDNA Ends (5’ RACE) Method

The used 5’ RACE method in this study was introduced by Xiang Zong Shi and Jarvis to identify GC rich sequences (25).

3.2.1. RNA Extraction

Total RNA was extracted from the exponential phase of yeast using the RNeasy Protect Mini Kit (Qiagen, Germany). For mechanical disruption, the yeast cells were sonicated (Hieschler, Germany; cycle 0.6; Amplitude 70%) with acid-washed glass beads (0.4 - 0.5 mm diameter). The quantity and quality of the extracted RNA was measured using the Nanodrop 1000 spectrophotometer (Thermo Scientific, USA).

3.2.2. cDNA Synthesis

For complementary DNA (cDNA) synthesis in a total volume of 20 μL, 2 μg of total RNA from organism was heated at 65°C for 7 minutes, followed by cooling on ice. A master mixture contained 4 μL of 5× reverse transcriptase (RT) buffer (Fermentas, Canada); 1 μL of gene-specific primer (5 ACO1R1; Table 1) (10 pmol/μL); dNTP Mix, 1 μL (10 mM); Ribolock, 1 μL(20 U) (Fermentas, Canada); and diethylpyrocarbonate (DEPC) treated water. Four hundred units (2 μL) of Moloney murine leukemia virus (M-MuLV) reverse transcriptase enzyme (Fermentas, Canada) were added. The RT temperature was set to 42°C for 1 hour and 70°C for 7 minutes. Subsequently, the reaction was treated with an equal volume of 0.6 N NaOH for 20 minutes at 65°C to degrade the mRNA and then the cDNA was precipitated by adjusting the solution to 0.5 M NaCl and adding absolute ethanol (2.5 times the total volume). After incubation at −80°C for 2 hours, the cDNA was separated using centrifugation for 30 minutes at 13,000 × g, washed with 70% ethanol, and dissolved in 10 μL of a commercial ligation-adaptor provided as follows:

Table 1.

Oligonucleotides Used in the Study

NameTmSequence (5’ to 3’)Reference
UPM Long77°CCTAATACGACTCACTATAGGGCAAGCAGTGGTATCAACGCAGAGT(23)
UPM Short53°CCTAATACGACTCACTATAGGGC(23)
Sense Adaptor (SA)71°CAAGCAGTGGTATCAACGCAGAGTGNNNNN(23)
Antisense Adaptor (ASA)62°Cp-ACTCTGCGTTGATACCACTGCTT (5’-phosphorylated)(23)
5ACO1R157°CGTTCAAGTCAATGGCTCTGGIn this study
5ACO1R252°CAGTGGACAGTAGATGGAGTGGIn this study
5ACO1R358°CTGGGCGGTAGCATCTTGACIn this study
3ACO1F54°CACCCAAGATATTGAAAGAGGIn this study
3ACO1R58°CAATCTTGGTTCCAAAGCAGCIn this study
M13 F54°CGTAAAACGACGGCCAGUniversal
M13 R56°CCAGGAAACAGCTATGACUniversal

3.2.3. Adaptor Ligation

The ligation adaptor was provided by mixing two commercial oligonucleotides, named sense adaptor and antisense adaptor (Table 1), in an annealing buffer comprised of 10 mM Tris-HCl (pH. 8.0) and 50 mM NaCl. The mixture was placed in a boiling water bath, which was then allowed to slowly cool to room temperature, and the resulting adaptor was added to the first-strand cDNA, as explained above. Subsequently, a ligation reaction was done in a total volume of 10 μL, using 3 μL of the cDNA-adaptor mixture, 10X T4 DNA ligase buffer, 1.0 μL; and T4 DNA ligase, 1.0 μL. This mixture was placed overnight at room temperature and then the adaptor-ligated cDNA was diluted fivefold with ddH2O and used as the template for subsequent PCRs.

3.2.4. PCR amplification

Primary PCR mixtures included 12.5 μL of Taq DNA Polymerase Master Mix RED (Ampliqon), 2.0 μL of the diluted adaptor ligated first-strand cDNA described above, 2 μL of a gene-specific primer (5ACO1R2; Table 1), 2 μL of a primer complementary to the 3’ end of the adaptor sequence (UPM long; Table 1) and 6.5 μL of ddH2O in a total volume of 25 μL. Touchdown cycling condition was used for the primary amplification, with an initial denaturation step for 5 minutes at 95°C, five PCR cycles each at two different annealing temperatures (63°C and 56°C), and then 30 PCR cycles with an annealing temperature of 50°C. Each PCR cycle included a 40 seconds denaturation step at 95°C, a 40 seconds annealing step at the temperatures indicated above, and a 90 seconds extension step at 72°C and a final extension at 72°C for 10 minutes.

Mixture of nested secondary PCR included 12.5 μL of Taq DNA Polymerase Master Mix RED (Ampliqon), 1.0 μL of the 1:5 diluted PCR product described above, 1 μL of a gene-specific primer (5ACO1R3; Table 1), 1 μL of a primer complementary to the UPM long sequence (UPM short; Table 1) and 9.5 μL of ddH2O in a total volume of 25 μL. Touchdown cycling condition was used for the secondary amplification, with an initial denaturation step for 5 minutes at 95°C, five PCR cycles each at two different annealing temperatures (63°C and 56°C), and then 30 PCR cycles with an annealing temperature of 53°C. Each PCR cycle included a 40 seconds denaturation step at 95°C, a 40 seconds annealing step at the temperatures indicated above, and a 90 seconds extension step at 72°C and a final extension at 72°C for 10 minutes.

3.3. Multiple Sequence Alignment and Degenerate PCR

The 3ACO1F primer was designed from the fragment obtained by 5’ RACE method and degenerate PCR primer 3ACO1R was designed by the alignment and comparison of previously published ACO sequences from five Candida species; C. tropicalis, C. albicans, C. dubliniensis, C. parapsilosis, and C. orthopsilosis. All sequences were obtained from the GenBank (accession numbers XM_002548708, XM_711132, XM_002422275, HE605206, and XM_003869506 respectively). The PCR reaction mixture included 12.5 μL of Taq DNA Polymerase Master Mix RED (Ampliqon), 1 μL (10 pmole/μL) of each forward (3ACO1F) and reverse (3ACO1R) primers, 1 μL C. krusei genomic DNA and 9.5 μL ddH2O in a total volume of 25 μL. The cycling parameters were 95°C for 5 minutes: 30 cycles of denaturation for 40 seconds at 95°C, annealing for 40 seconds at 60°C, an extension for 90 seconds at 72°C, with a final extension of 10 minutes.

3.4. Cloning and Sequencing of PCR Products

All PCR products were cloned using a TOPO TA-cloning kit (Invitrogen, USA) following the manufacturer’s instructions. Products were ligated into the TA vector and then transformed into chemically competent cells using heat shock at 42°C. Cells were then mixed with 250 μL of SOC medium, incubated at 37°C for one hour, and then 50 μL of the cell culture was plated onto X-gal/IPTG/Amp/LB plates. Following overnight incubation at 37°C, a single white colony (~2 mm diameter) was picked from the plate and inoculated into 5 mL of Luria Bertani (LB) broth. After overnight incubation at 37°C, 1 mL of the culture was used to prepare glycerol stocks, and the other 4 mL was used to extract plasmid. Recombinant plasmid was purified using a plasmid extraction Kit (Intron biotechnology) following the manufacturer’s instruction. The recombinant plasmids were screened by M13 forward and reverse primers with the following PCR protocol: 94°C for 5 minutes; 30 cycles of 94°C for 1 minutes, 55°C for 1 minute; 72°C for 1 minute; and 72°C for 7 minutes. PCR products were analyzed by agarose (Merck, Germany) gel electrophoresis. The plasmids containing unknown DNA were sequenced by M13 primers (Bioneer, Korea).

3.5. Nucleotide and Protein Sequence Accession Number

The partial nucleotide sequence of the C. krusei aconitase gene and its annotation was submitted to the GenBank under the assigned accession number KJ009334.

3.6. Molecular Phylogenetic Analyzes

The following aconitase amino acid sequences from different organisms were obtained from the GenBank database. Candida tropicalis, XP_002548754.1; C. albicans, XP_716225.1; C. parapsilosis, CCE42927.1; C. orthopsilosis, XP_003869555.1; C. dubliniensis, XP_002422320.1; C. glabrata, XP_445684.1; Pichia pastoris, XP_002489444.1; Scheffersomyces stipitis, XP_001386080.1; Lodderomyces elongisporus, XP_001525166.1; Clavispora lusitaniae, XP_002619843.1; Piromyces sp. CAA76360.1; Schizosaccharomyces pombe, NP_594031.2; Penicillium marneffei, XP_002147388.1; Aspergillus oryzae, XP_001819597.2; A. flavus, XP_002374915.1; A. terreus, AAC61778.1; A. fumigatus, XP_750430.1; A. clavatus, XP_001269477.1; A. niger, XP_001393703.1; Trichophyton rubrum, XP_003231764.1; Arthrodermaotae, XP_002847793.1; A. gypseum, XP_003177439.1; Coccidioides posadasii, ABH10644.1; Paracoccidioides brasiliensis, XP_002790154.1; Ajellomyces dermatitidis, XP_002627881.1; A. capsulatus, XP_001539799.1; Cryptococcus gattii, XP_003193698.1; C. neoformans, XP_570245.1; Homo sapiens, CU012932.1.

Nucleotide and amino acid sequences were aligned using ClustalW online software. Sequences were analyzed by the Unweighted-Pair-Group method with arithmetic mean (UPGMA) and neighbor joining (NJ) using Geneious R7 software. Bootstrap probability with 2000 replications was introduced to assess the statistical significance of the groups in the phylogenetic tree.

3.7. Active Site and Binding Site Characterization

3D structure, active site and binding site residues of C. krusei ACO was identified with COFACTOR online service (http://zhanglab.ccmb.med.umich.edu/COFACTOR/) to determine active site and binding site amino acid residues differences between C. krusei ACO and other Candida species and human aconitases to design new antifungal drugs.

4. Results

In the current study, 1419 nucleotides of C. kruseiACO gene clarified the first 300 nucleotides of 3’ region of C. krusei ACO gene obtained by 5’-RACE method. 1197 nucleotides of this gene were identified by PCR with 3ACO1F primer as both forward and reverse primers, whereas the degenerate PCR primer designed by alignment and comparison of previously published ACO sequences from five Candida species was unable to amplify this gene (Figure 1). Seventy-eight primary nucleotides overlapped with the sequence obtained from the 5’-RACE method. Finally, sequenced nucleotides were submitted to the GenBank as a partial sequence of C. krusei ACO gene.

Partial Nucleotide Sequence of the Aconitase Gene of Candida krusei
The nucleotide sequences were determined from the cloned PCR amplified Candida krusei aconitase gene.

Table 2 exhibits the pairwise nucleotide identities of Candida species calculated from the nucleotide sequences of the ACO genes. Candida albicans and C. dubliniensis (96%), C. albicans and C. tropicalis (90%), C. dubliniensis and C. tropicalis (90%) and C. orthopsilosis and C. parapsilosis (94%) had more than 90% identity. Candida krusei had maximum identity with C. tropicalis (84%).It also shows the pairwise amino acid identities among different Candida species. Based on the amino acid sequences, C. krusei had maximum identity with C. tropicalis (85%); however C. albicans and C. dubliniensis (99%), C. albicans and C. tropicalis (94%), C. dubliniensis and C. tropicalis (94%) and C. orthopsilosis and C. parapsilosis (98%) had more than 90% identity.

Table 2.

Levels of Aconitase Nucleotide and Amino Acid Sequence Similarities for Candida krusei and Five Other Candida Speciesa

SpeciesAconitase Sequence Similarity
C. albicansC. dubliniensisC. tropicalisC. parapsilosisC. orthopsilosisC. krusei
Amino acid sequence similarity
C. albicans
C. dubliniensis96
C. tropicalis9090
C. parapsilosis868986
C. orthopsilosis85898694
C. krusei8382848282
Nucleotide sequence similarity
C. albicans9994898984
C. dubliniensis94858583
C. tropicalis888885
C. parapsilosis9881
C. orthopsilosis82
C. krusei

The comparison of Candida species ACO genes showed that intra species variations between C. krusei and other Candida species were relatively higher than the variation between common Candida species. Multiple alignments of the nucleotide (Figure 2) and amino acid (Figure 3) sequences of ACO of C. krusei with other Candida species and Homo sapiens as an out-group showed that C. krusei possessed characteristic sequences.

Partial Multiple Alignment of ACO Nucleotide Sequences of Candida Species and Homo sapiens as an Out-Group
Partial Multiple Alignment of ACO Nucleotide Sequences of Candida Species and Homo sapiens as an Out-Group
Partial Multiple Alignment of ACO Amino Acid Sequences of Candida Species and Homo sapiens as an Out-Group
Partial Multiple Alignment of ACO Amino Acid Sequences of Candida Species and Homo sapiens as an Out-Group

The major pathogenic Candida species were distinctly positioned by phylogenetic trees constructed using UPGMA and NJ with H. sapiens as an out-group (Figures 4 - 6). Although there were some difference in the degree of resolution when nucleotide and amino acid sequences were applied, the topologies acquired by both sequence types and tree estimation algorithms were consistent.

Unweighted-Pair-Group Method With Arithmetic Mean Trees of the Fungal Species for the ACO Genes
Unweighted-pair-group method with arithmetic mean trees based on nucleotide (a) and amino acid (b) sequences. The Homo sapiens represents an out-group. Bars exhibit the numbers of nucleotide and amino acid substitutions per nucleotide and amino acid sites.
Unweighted-Pair-Group Method With Arithmetic Mean Trees of the Common Candida Species for the ACO Genes
Unweighted-pair-group method with arithmetic mean trees based on nucleotide (a) and amino acid (b) sequences. The Homo sapiens represents an out-group. Bars exhibit the numbers of nucleotide and amino acid substitutions per nucleotide and amino acid sites.
Neighbor Joining Trees of the Common Candida Species for the ACO Genes
Neighbor joining trees based on nucleotide (a) and amino acid (b) sequences. The Homo sapiens represents an out-group. Bars exhibit the numbers of nucleotide and amino acid substitutions per nucleotide and amino acid sites.

The active site and binding site of C. krusei aconitase contain six and ten residues, respectively (Table 3). The active site amino acid residues of ACO in six Candida species and human were similar; just the amino acids were located in different positions in protein structure. The binding site residues of ACO in six Candida species are identical but one isoleucine residue of C. kruseiACO was not present in aconitase binding site of human and there was one aspartate residue in human aconitase that was not present in C. krusei and other Candida species.

Table 3.

Residue Information for the Active and Binding Sites of Aconitase

SpeciesActive Site ResiduesBinding Site Residues
C. kruseiAsp123, His124, His170, Asp188, His190, Glu285His124, Ile168, Ile169, His170, His190, Ser381, Cys382, Cys445, Cys448, Ile449
C. albicansAsp113, His114, His160, Asp178, His180, Glu275His114, Ile158, Ile159, His160, His180, Ser371, Cys372, Cys435, Cys438, Ile439, Asn460
C. dubliniensisAsp111, His112, His158, Asp176, His178, Glu273, Ser654His112, Ile156, Ile157, His158, His178, Ser369, Cys370, Cys433, Cys436, Ile437, Asn458
C. tropicalisAsp113, His114, His160, Asp178, His180, Glu275His114, Ile158, Ile159, His160, His180, Ser371, Cys372, Cys435, Cys438, Ile439, Asn460
C. parapsilosisAsp112, His113, His159, Asp177, His179, Glu274, Ser655His113, Ile157, Ile158, His159, His179, Ser370, Cys371, Cys434, Cys437, Ile438, Asn459
C. orthopsilosisAsp113, His114, His160, Asp178, His180, Glu275His114, Ile158, Ile159, His160, His180, Ser371, Cys372, Cys435, Cys438, Ile439, Asn460
Homo sapiensAsp114, His115, His161, Asp179, His181, Glu276, Ser656His115, Ile159, His161, Asp179, His181, Ser371, Cys372, Cys435, Cys438, Ile439, Asn460

5. Discussion

The current study aimed to identify C. krusei aconitase gene sequence and evaluate the phylogenetic relationships of C. krusei with other organisms and Candida species based on nucleotide and amino acid sequences of this gene. Also, the differences of active and binding sites and amino acid residues between C. krusei and human aconitase protein were determined to design new antifungal drugs targeting C. krusei aconitase. Since the length of this gene in Candida species is approximately 2340 base pair, it is believed that the length of this gene in C. krusei is also 2340 bp. In the current study, 1419 bp of aconitase gene of C. krusei were identified and sequenced by different methods including 5’ RACE method and degenerate PCR. Partial sequence was submitted to the GenBank and the accession number KJ009334 was assigned.

The 5’ RACE method in the current study was a new method introduced by Xian Zong Shi and Jarvis to identify GC rich sequences in 2006 (25). This method was slightly modified and used to identify a part of C. kruseiACO gene. Although the ACO gene of C. krusei is not GC rich but the results of this method to identify 5’ end of mRNA of this gene was satisfactory. Since 3’ RACE method was unable to identify the other part of this gene, degenerate PCR was used to identify unknown parts of ACO gene. Sequence analysis indicated that 1197 bp of this gene were amplified using PCR with 3ACO1F as both forward and reverse primers and degenerate PCR was unable to amplify this gene.

The ACO gene encodes protein aconitase involved in energy production in the Krebs cycle (19). Since this gene encodes a critical enzyme in different organisms including Candida species, little difference is expected between the nucleotide and amino acid sequences among different species of a genus. The results showed that not only nucleotide and amino acid sequences of this gene in C. krusei were very different from other species of Candida but also there were considerable differences between the other species of Candida. However, differences of nucleotide and amino acid sequences of this gene in C. krusei were higher than those of the other Candida species, which supports the hypothesis that C. krusei may be taxonomically separated from Candida genus. Phylogenetic analysis of amino acid and nucleotide sequences of aconitase gene showed that although C. krusei takes place in the Hemiascomycetes clade (26), this species takes place beside P. pastoris in phylogenetic trees and separated from other Candida species in this clade because of more nucleotide differences with other species (Figure 4).

The sequences were analyzed using UPGMA and NJ. There was very little difference between UPGMA (Figure 5), NJ (Figure 6), and phylogenetic trees based on the nucleotide or amino acid sequences of the common Candida species ACO genes. It seems to be due to bootstrap values and differences in the analysis methods. Moreover, different nucleotide codons of an amino acid may cause differences in nucleotide- and amino acid-based phylogenetic trees. Despite some differences in the degree of resolution when nucleotide and amino acid sequences were applied, the topologies obtained with both sequence types and the tree estimation algorithms were mutually consistent.

Aconitase enzyme activates in the presence of a [4Fe-4S] cluster (27). There is a [3Fe-4S] cluster in inactive aconitase ligated to the active site via three cysteine residues, Cys382, Cys445, and Cys448 in C. krusei (Cys372, Cys435, and Cys438 in human) and this enzyme is activated upon addition of a fourth iron atom (termed Fea) (28). This fourth iron atom has a slightly different electronic character than the others, and is the active site of the enzymatic reaction (28). There is no adjacent cysteine residue to ligate the fourth Fe (Fea), which instead binds a hydroxyl ion as its fourth ligand. The hydroxyl ion is bound to histidine (His) and aspartate (Asp) and becomes water (H2O) upon addition of substrate where it is further hydrogen bonded to the carboxyl group of citrate or isocitrate. Isocitrate or citrate is bound to Fea via an oxygen on its carboxyl group as well as its hydroxyl group, changing the coordination state of Fea from 4- to 6-coordinate (29). Coordination of the hydroxyl group makes it a better leaving group and facilitates the dehydration reaction that forms cis-aconitate from either citrate or isocitrate (27, 30). Serine (Ser) 656 alkoxide (in human) attacks the hydrogen on hydroxyl group of citrate or isocitrate creating a carbanion intermediate. Depending on which carboxyl oxygen coordinates the iron atom, cis-aconitate can be converted to citrate or isocitrate upon addition of water (30). In the TCA cycle, cis-aconitate is converted to isocitrate.

This is report the first report on the identification, sequencing and phylogenetic analysis of C. krusei using ACO gene sequences. Differences of C. kruseiACO gene sequences from other Candida species are prominent and may facilitate the design of probes for rapid and specific identification of C. krusei from other common Candida species. Moreover, according to aconitase activity mechanism, there are not significant amino acid residue differences between C. krusei and human aconitases within the active site cleft, which may block the use of aconitase as a molecular target applied to design new anti-fungal drugs that selectively block this enzyme.

Acknowledgements

References

  • 1.

    Dixon DM, McNeil MM, Cohen ML, Gellin BG, La Montagne JR. Fungal infections: a growing threat. Public Health Rep. 1996;111(3):226-35. [PubMed ID: 8643813].

  • 2.

    Maccallum DM. Hosting infection: experimental models to assay Candida virulence. Int J Microbiol. 2012;2012:363764. [PubMed ID: 22235206]. https://doi.org/10.1155/2012/363764.

  • 3.

    Miceli MH, Diaz JA, Lee SA. Emerging opportunistic yeast infections. Lancet Infect Dis. 2011;11(2):142-51. [PubMed ID: 21272794]. https://doi.org/10.1016/s1473-3099(10)70218-8.

  • 4.

    Pfaller MA, Diekema DJ. Epidemiology of invasive candidiasis: a persistent public health problem. Clin Microbiol Rev. 2007;20(1):133-63. [PubMed ID: 17223626]. https://doi.org/10.1128/CMR.00029-06.

  • 5.

    Bodey G, Bueltmann B, Duguid W, Gibbs D, Hanak H, Hotchi M, et al. Fungal infections in cancer patients: an international autopsy survey. Eur J Clin Microbiol Infect Dis. 1992;11(2):99-109. [PubMed ID: 1396746].

  • 6.

    Boon AP, O'Brien D, Adams DH. 10 year review of invasive aspergillosis detected at necropsy. J Clin Pathol. 1991;44(6):452-4. [PubMed ID: 2066421].

  • 7.

    Denning DW, Stevens DA. Antifungal and surgical treatment of invasive aspergillosis: review of 2,121 published cases. Rev Infect Dis. 1990;12(6):1147-201. [PubMed ID: 2267490].

  • 8.

    Iwen PC, Reed EC, Armitage JO, Bierman PJ, Kessinger A, Vose JM, et al. Nosocomial invasive aspergillosis in lymphoma patients treated with bone marrow or peripheral stem cell transplants. Infect Control Hosp Epidemiol. 1993;14(3):131-9. [PubMed ID: 8478525].

  • 9.

    Morrison VA, Haake RJ, Weisdorf DJ. The spectrum of non-Candida fungal infections following bone marrow transplantation. Medicine (Baltimore). 1993;72(2):78-89. [PubMed ID: 8479326].

  • 10.

    Marichal P, Gorrens J, Coene MC, Le Jeune L, Vanden Bossche H. Origin of differences in susceptibility of Candida krusei to azole antifungal agents. Mycoses. 1995;38(3-4):111-7. [PubMed ID: 7477085].

  • 11.

    Marichal P, Vanden Bossche H. Mechanisms of resistance to azole antifungals. Acta Biochim Pol. 1995;42(4):509-16. [PubMed ID: 8852341].

  • 12.

    Samaranayake YH, Samaranayake LP. Candida krusei: biology, epidemiology, pathogenicity and clinical manifestations of an emerging pathogen. J Med Microbiol. 1994;41(5):295-310. [PubMed ID: 7966200]. https://doi.org/10.1099/00222615-41-5-295.

  • 13.

    Francois IE, Aerts AM, Cammue BP, Thevissen K. Currently used antimycotics: spectrum, mode of action and resistance occurrence. Curr Drug Targets. 2005;6(8):895-907. [PubMed ID: 16375673].

  • 14.

    Sheehan DJ, Hitchcock CA, Sibley CM. Current and emerging azole antifungal agents. Clin Microbiol Rev. 1999;12(1):40-79. [PubMed ID: 9880474].

  • 15.

    Wu JJ, Pang KR, Huang DB, Tyring SK. Therapy of systemic fungal infections. Dermatol Ther. 2004;17(6):532-8. [PubMed ID: 15571503]. https://doi.org/10.1111/j.1396-0296.2004.04057.x.

  • 16.

    Kontoyiannis DP, Mantadakis E, Samonis G. Systemic mycoses in the immunocompromised host: an update in antifungal therapy. J Hosp Infect. 2003;53(4):243-58. [PubMed ID: 12660121].

  • 17.

    Nelson DL, Cox MM. Lehninger Principles of Biochemistry. 4th ed. New York, USA: W.H. Freeman & Company; 2004.

  • 18.

    Regev-Rudzki N, Karniely S, Ben-Haim NN, Pines O. Yeast aconitase in two locations and two metabolic pathways: seeing small amounts is believing. Mol Biol Cell. 2005;16(9):4163-71. [PubMed ID: 15975908]. https://doi.org/10.1091/mbc.E04-11-1028.

  • 19.

    Tong WH, Rouault TA. Metabolic regulation of citrate and iron by aconitases: role of iron-sulfur cluster biogenesis. Biometals. 2007;20(3-4):549-64. [PubMed ID: 17205209]. https://doi.org/10.1007/s10534-006-9047-6.

  • 20.

    Doi M, Homma M, Chindamporn A, Tanaka K. Estimation of chromosome number and size by pulsed-field gel electrophoresis (PFGE) in medically important Candida species. J Gen Microbiol. 1992;138(10):2243-51. [PubMed ID: 1479351]. https://doi.org/10.1099/00221287-138-10-2243.

  • 21.

    Yamada YKK, editor. Taxonomic significance of the coenzyme Q system in yeasts and yeast-like fungi. In: Iwata K (ed) Yeasts and yeast like organisms in medical science. Tokyo International Special Symposium on Yeasts. 1972. p. 61-9.

  • 22.

    Schaefer BC. Revolutions in rapid amplification of cDNA ends: new strategies for polymerase chain reaction cloning of full-length cDNA ends. Anal Biochem. 1995;227(2):255-73. [PubMed ID: 7573945]. https://doi.org/10.1006/abio.1995.1279.

  • 23.

    Ohara O, Dorit RL, Gilbert W. One-sided polymerase chain reaction: the amplification of cDNA. Proc Natl Acad Sci U S A. 1989;86(15):5673-7. [PubMed ID: 2788276].

  • 24.

    Loh EY, Elliott JF, Cwirla S, Lanier LL, Davis MM. Polymerase chain reaction with single-sided specificity: analysis of T cell receptor delta chain. Science. 1989;243(4888):217-20. [PubMed ID: 2463672].

  • 25.

    Shi X, Jarvis DL. A new rapid amplification of cDNA ends method for extremely guanine plus cytosine-rich genes. Anal Biochem. 2006;356(2):222-8. [PubMed ID: 16875657]. https://doi.org/10.1016/j.ab.2006.06.028.

  • 26.

    Diezmann S, Cox CJ, Schonian G, Vilgalys RJ, Mitchell TG. Phylogeny and evolution of medical species of Candida and related taxa: a multigenic analysis. J Clin Microbiol. 2004;42(12):5624-35. [PubMed ID: 15583292]. https://doi.org/10.1128/JCM.42.12.5624-5635.2004.

  • 27.

    Garrett RHGC. Biochemistry. 4th ed. New York City: Saunders College Publishing; 1999.

  • 28.

    Emptage MH, Kent TA, Kennedy MC, Beinert H, Munck E. Mossbauer and EPR studies of activated aconitase: development of a localized valence state at a subsite of the [4Fe-4S] cluster on binding of citrate. Proc Natl Acad Sci U S A. 1983;80(15):4674-8. [PubMed ID: 6308639].

  • 29.

    Kennedy MC, Werst M, Telser J, Emptage MH, Beinert H, Hoffman BM. Mode of substrate carboxyl binding to the [4Fe-4S]+ cluster of reduced aconitase as studied by 17O and 13C electron-nuclear double resonance spectroscopy. Proc Natl Acad Sci U S A. 1987;84(24):8854-8. [PubMed ID: 3480514].

  • 30.

    Cleland WW, Kreevoy MM. Low-barrier hydrogen bonds and enzymic catalysis. Science. 1994;264(5167):1887-90. [PubMed ID: 8009219].