Abstract

Based on the ability of nitrate reductase synthesis, Corynebacterium pseudotuberculosis is classified into two biovars: Ovis and Equi. Due to the presence of nitrate reductase, the Equi biovar can survive in absence of oxygen. On the other hand, Ovis biovar that does not have nitrate reductase is able to adapt to various ecological niches and can grow on certain carbon sources. Apart from these two biovars, some other strains are also able to carry out the reduction of nitrate. The enzymes that are involved in electron transport chain are also identified by in silico methods. Findings about pathogen metabolism can contribute to the identification of relationship between nitrate reductase and the C. pseudotuberculosis pathogenicity, virulence factors, and discovery of drug targets.

1. Introduction

Corynebacterium pseudotuberculosis is a Gram-positive facultative intracellular pathogen [111]. C. pseudotuberculosis can be classified into two biovars, based on their ability to convert nitrate to nitrite. The nitrate-positive biovar is Equi, which causes ulcerative lymphangitis in equines, while the nitrate-negative biovar is known as Ovis, which is the etiologic agent of caseous lymphadenitis in small ruminants [9]. Both diseases are globally distributed and cause large economic losses to goat, sheep, horse, and cattle farmers.

The nitrate reduction is associated with the bacterium’s ability to breathe in the absence of oxygen and having two different metabolic pathways, (1) respiratory nitrate reductase and (2) dissimilatory nitrate reduction. In the first pathway, the denitrification process takes place where the nitrate is sequentially reduced to nitrite, nitric oxide, nitrous oxide, and finally to dinitrogen [12, 13]. In the second pathway nitrate is directly converted into ammonia, which is secreted from the cell; this process can be performed by organisms with the nrf gene. This is a less common method of nitrate reduction than denitrification in most ecosystems [1214]. There are four classes of nitrate reductases, one in eukaryotes and three are in prokaryotes. Prokaryotic nitrate reductases include a class of assimilatory enzymes and two classes of respiratory enzymes; all contain a guanylate molybdenum cofactor but differ in their substructures, cellular location, and requirement for cofactor. Variability among enzyme is also found into the classes [12, 13].

Among the Actinobacteria, respiratory nitrate reductase (Nar) is observed [1518]. This class of enzymes is mostly membrane-bound and is having three different subunits: respiratory nitrate reductase from Pseudomonas stutzeri [19, 20], nitrate reductase A, and nitrate reductase Z from Escherichia coli K-12 [2123]. In Corynebacterium glutamicum and Mycobacterium tuberculosis the nitrate reduction reaction is catalyzed by the operons narKGHJI and narGHJI, respectively [15, 16]. In these organisms, the first step in nitrate assimilation is nitrate reduction () to nitrite () and the second step is nitrite reduced () to ammonia () [15, 16].

Nitrate respiration is associated with virulence and adoption of pathogenic organisms. In Brucella suis it helps in intramacrophagic multiplication, in M. bovis it allows survival inside the host, and in M. tuberculosis it confers resistance to stress [17, 24, 25]. However, the role of nitrate respiration is not known in C. pseudotuberculosis. We presume that Nar may be associated with pathogenicity in C. pseudotuberculosis biovar Equi, providing a generate energy pathway when it grows under oxygen deprived condition allowing the survival of C. pseudotuberculosis inside the host.

With the increase in genome sequencing projects, understanding of the biological capabilities of the organisms has also increased. The availability of the annotated genomes allows the information of metabolic pathways to be interpreted and used for computational reconstruction at genomic scale. Thus the reconstruction of metabolic pathways through comparative analysis can help identify potential targets and find the mechanisms that cause a disease.

In the present study, the genomes of C. pseudotuberculosis strains deposited in GenBank were used to identify the genes involved in the nitrate reduction pathway; a comparative genomic study was performed on 15 C. pseudotuberculosis strains of the biovars Ovis and Equi to identify the molecular bases of nitrate reduction in each C. pseudotuberculosis biovar Equi. Thus, we hope that understanding of Nar’s role in C. pseudotuberculosis can help in identifying potent targets both for the development of more effective diagnostics and therapeutics and for the treatment and prevention of disease.

2. Material and Methods

2.1. Identification of Nitrate Reductase Genes, Metabolic Pathways, and Nitrate Reduction Enzymes

C. pseudotuberculosis genomes Table 1 were obtained from the NCBI GenBank. The sequence and arrangement analysis of nitrate reductase genes in narKGHJI operon was performed using the Artemis software (http://www.sanger.ac.uk/Software/Artemis) and manual curation based annotations are performed. Pathway Tools software developed by SRI International [45] and MetaCyc database [46] were used for computational prediction of metabolic pathways in C. pseudotuberculosis.

2.2. narKGHJI Operon Analysis

The families of proteins as well as the conserved protein domains present in the operon narKGHJI and adjacent genes were analyzed using InterProScan database (https://www.ebi.ac.uk/interpro). Prediction of transmembrane helices in protein was performed using TMHMM Server v. 2.0 (http://www.cbs.dtu.dk/services/TMHMM/). The similarity analysis of genes involved in the nitrate reduction pathway was performed using the UniProtKB database in category BLASTp. Clustal Omega was used to multiple sequence alignment between operon narKGHJI of C. pseudotuberculosis biovar Equi and homologous sequences of multiple organisms (http://www.ebi.ac.uk/Tools/msa/clustalo/), and Jalview was used to analyze similarity (http://www.jalview.org/) among these sequences.

3. Results and Discussion

3.1. Correlation between Nitrate Reduction Phenotype and Genotype

Nineteen C. pseudotuberculosis strains that had genome information in GenBank are presented Table 1. To perform the phenotyping analysis, 18 strains that were sequenced in our lab were used. All these strains were biochemically tested in order to verify their ability of nitrate reduction; see Table 1. Among the C. pseudotuberculosis strains, 31, 258, 262, and MB20 strains showed concordance between phenotypic and genotypic results. Strains 1/06-A, 316, 162, and CIP52.97 are found positive phenotype for reduction of nitrate; however they do not have genes involved in nitrate reduction in their genome.

3.2. Genes Involved in Respiratory Nitrate Reductase Metabolic Pathway

We first compared the autoassigned pathways by the Pathologic software with the pathways in databases such as KEGG (http://www.genome.jp/kegg/) and resolved possible discrepancies using manual curation. Comparative genome analysis revealed that C. pseudotuberculosis biovar Equi possess narKGHJI gene clusters that participate via the respiratory anaerobic process of the nitrate reduction similarly to E. coli (Figure S8) [47]. The C. pseudotuberculosis narKGHJI gene cluster showed great similarity with the protein sequences found in other Actinomycetes, such as C. diphtheriae, C. glutamicum, and M. tuberculosis (Table 2). All biovar Ovis strains do not present any gene of the narKGHJI operon in their genomes (Figure 1).

The nitrate locus that contains narKGHJI operon in C. pseudotuberculosis is composed of gene cluster encoding the molybdopterins moeB, moaE, molB, molA, moeY, moaC, moeA, and moa (Figure 2). Molybdopterin is a cofactor that is indispensable for activity of nitrate reductase [48]. The organization and conservation of the narKGHJI gene cluster was assessed among C. pseudotuberculosis and C. diphtheriae, C. glutamicum, M. tuberculosis and, E. coli. The analysis showed that narKGHJI gene cluster is conserved among these species; however gene orientation of the five genes is only conserved among C. pseudotuberculosis, C. diphtheriae, and C. glutamicum. In addition, in C. pseudotuberculosis we found narT gene is located at the upstream of the narKGHJI cluster. The narT gene encodes a nitrate transporter and is not present in E. coli and C. glutamicum.

3.3. Enzymes Involved in Electron Transport Chain

We found conserved functional domains of the genes present in the narGHI complex (involved in electron transport chain) in C. pseudotuberculosis biovar Equi (Table 3). To infer the functionality of narGHI complex, we used tridimensional structure from E. coli (PDB accession 1Q16) as a model for nar genes, as the structure of these enzymes is resolved in these bacteria.

In the narGHI complex in Figure 2, the narG gene is a member of a superfamily of enzymes that use a Mo-bisMGD cofactor (bisMGD) for their catalytic activity. Residues present in the active site of the molybdenum atom (Mo) are highly conserved in NarG subunits of Gram-negative and Gram-positive bacterial species [21]. The NarH belongs to the superfamily of electron transfer subunit (ferredoxins) of bacterial oxidoreductases. NarH has three clusters (FS1, FS2, and FS3) and a cluster (FS4) [21, 42, 44]. The general structure of NarH in E. coli is composed of a central region of a linker region that is inserted between the two subdomains in coordination of the iron-sulfur clusters and forms an extended connection between the NarG and NarH subunits and an extended C-terminus. In C. pseudotuberculosis it was observed that the nine cysteine residues responsible for coordinating three clusters and one cluster are conserved in all analyzed sequences. NarI is the transmembrane subunit that anchors NarGH to the cytoplasmic side of the membrane and provides the connection between the quinone and the site of oxidation, referred to as “Site Q.” This subunit coordinates two hemes, one of which is located towards the side of cytoplasmic end of NarI and is referred to as the proximal heme (heme ), while the other is located towards the periplasmic side of the same and is referred to as the distal heme (heme ). Such hemes act as mediators in electron transfer from Site Q to the clusters [Fe-S] in NarH. NarI of C. pseudotuberculosis presented five transmembrane helices S1 Figure, which is in agreement with the literature when compared to the E. coli crystal structure [49, 50]. This C-terminal segment presents highly conserved residues in the family NarI proteins and also is present in C. pseudotuberculosis. It is interesting to note that NarI is considerably preserved in relation to its primary sequence between orthologous proteins found in species of the Corynebacterium genus, particularly in relation to C. diphtheriae, which has 96% identity and 98% similarity. These results suggest that the enzymatic complex narGHI is probably very similar with respect to its three dimensional structure among the species of C. pseudotuberculosis biovar Equi and C. diphtheriae in view of the high identity found among NarG, NaH, and NarI subunits of these species. The C. pseudotuberculosis biovar Equi presents 48% and 43% identity in relation to the orthologous proteins from M. tuberculosis and B. subtilis, respectively.

3.4. NarJ (Molybdenum-Cofactor-Assembly Chaperone)

The NarJ is a chaperone that is involved in folding, maturation, and molybdenum cofactor insertion of nitrate reductase in E. coli [51, 52]. The similarity values found between C. pseudotuberculosis biovar Equi and other sequences of orthologous proteins by BLAST, including the genus Corynebacterium, were slightly lower than those found in other proteins (NarG, NarH, and NarI) that form the three subunits of the enzyme nitrate reductase (Table 2). However, narJ remained conserved between C. pseudotuberculosis biovar Equi strains, and results of biochemical tests (Table 1) indicate that its function has remained preserved, since NarJ is important for nitrate reductase expression in others organisms [53].

3.5. Nitrate and Nitrite Transporters

The nitrate and nitrite export is mediated by nitrate and nitrite transporters that are members of the major facilitator superfamily (MFS) [54]. In bacteria such carriers are the narK family represented by NarK and NarU proteins in E. coli that transport both nitrate and nitrite [47]. However, NarK is expressed at higher levels than NarU and it remains unclear till date if NarK is a nitrate/nitrite antiporter or symporter [55]. In C. pseudotuberculosis biovar Equi, we found that the two genes narK and narT are responsible for coding two probable nitrate and nitrite carriers and in C. pseudotuberculosis biovar Equi is present at the upstream of the narGHJI operon (Figure 2). Both proteins are found to be composed of 12 transmembrane helices similar to other members of the superfamily of bacterial transporters; see S2 and S3 Figures. In C. glutamicum the narK is in close proximity to narG, narH, narJ, and narI, which together form the narKGHJI operon [56], but in E. coli, narGHJI operon is transcribed individually and narT is upstream of the operon narKGHJI [57]. Since we did only the in silico analysis, we cannot state the precise transcription and regulation of these genes in vivo. We presume that this narKGHJI cluster in C. pseudotuberculosis biovar Equi may be regulated together as found in C. glutamicum. However, experimental confirmation is required.

3.6. Nitrite Reductase in C. pseudotuberculosis

Genes in C. pseudotuberculosis are found to encode two different types of nitrite reductase enzymes. One is the nitric oxide-forming enzyme that is encoded by nirK gene and shown to be a copper dependent enzyme. The second groups of enzymes are ammonium nitrite reductase that are encoded by nrfA gene that encodes catalytic subunit and nrfH that encodes for a membrane associated electron transfer subunit. In all C. pseudotuberculosis strains we found NrfA and NrfH subunits of nitrite reductase NrfAH complex that catalyze the reduction of nitrite to ammonium. Beside C. pseudotuberculosis, such genes are found only in C. ulcerans and are absent in all other species from this genus. In C. pseudotuberculosis, the two genes nrfA and nrfH are located in a “cluster” of genes, S4 Figure, responsible for coding biogenesis system I proteins cytochrome c that appears to be a common feature among bacterial genomes. C. pseudotuberculosis also has a nitrite reductase enzyme likely dependent on copper (CuNIR). The gene for CuNIR is nirK that encodes a protein composed of 882 residues and is conserved in all strains of C. pseudotuberculosis; see Figure S5 of the Supplementary Material available online at https://doi.org/10.1155/2017/9481756. The sequence similarity analysis of nirK was performed using the BLAST and it was observed that such protein has an N-terminal region of about 500 amino acids, S6 Figure, compared to the most homologous proteins that have been experimentally well characterized [5861]. Within the genus Corynebacterium, the CuNIR proteins are found to have similar sizes to that of C. pseudotuberculosis, such as Corynebacterium durum (WP_006063127.1) and Corynebacterium vitaeruminis (WP_051483598.1), with 40% identity between the sequences; see Table S1. Furthermore, it was found that the N-terminal region of about 500 residues is probably composed of 13 transmembrane helices in C. pseudotuberculosis, S7 Figure, which indicates that it can be a membrane associated protein. To verify the residues conservation, we used G. kaustophilus CuNIR proteins as a reference sequence because its crystal structure has recently been determined (pdb∣3WIA∣A). All residues responsible for binding to the copper atoms and the catalytic site were analyzed. It is observed that the C-terminal region is 100% conserved in C. pseudotuberculosis; see S6 Figure.

3.7. Reconstruction of C. pseudotuberculosis Biovar Equi Metabolic Network

A metabolic network reconstruction is assembled piece-by-piece by compiling data on known enzymes, and annotated genome together with data from the literature is used to assemble a network reconstruction [62]. In C. pseudotuberculosis after the similarity analysis we observed that genomic sequences participating in the nitrate reduction pathway are conserved, and the residues were analyzed, indicating a possible conserved function. Thus we can infer that C. pseudotuberculosis biovar Equi is different from that occurring in C. glutamicum, where nitrate is excreted out of the cells as the major end product under the conditions of anaerobic respiration, and nitrate and nitrite are not used as sources of nitrogen for aerobic or anaerobic growth [18] so C. pseudotuberculosis biovar Equi has a likely route of denitrification composed of the nitrate reductase Nar, nitrite reductase dependent copper, and nitric oxide reductase quinol-dependent, and this is a way of incomplete denitrification, in which nitrous oxide is the final product, and, therefore, there is no presence of nitrous oxide reductase enzyme and subsequent reduction to nitrogen gas; see Figure 3. Furthermore C. pseudotuberculosis biovar Equi also shows the pathway of ammonification in dissimilatory nitrate reduction to ammonium, what has been proposed as catalyzed by the enzyme nitrate reductase Nar and nitrite reductase complex NrfAH, a membrane-associated cytochrome; see Figure 3. E. coli not only reduce nitrite to ammonia but also catalyze the reduction of nitrite to NO, which can be used as an electron acceptor for anaerobic respiration [63]. The presence of enzymes and pathways of nitrate reduction proposals was also compared between biovar Ovis and biovar Equi strains and with respect to enzymes, the main difference was the presence of the nitrate locus that contains the genes encoding the molybdopterin and the genes encoding the nitrate reductase and a region situated upstream of the ansA gene and this locus has a 19,606 pb length.

4. Conclusions

In this study, we found that C. pseudotuberculosis strains 1/06-A, 316, 162, and CIP52.97 are able to do nitrate reduction but in their genomes they do not have genes associated with nitrate reduction, warning us about resequencing of these genomes. The C. pseudotuberculosis nitrate locus is composed of narKGHJI gene cluster and by genes encoding the molybdopterins moeB, moaE, molB, molA, moeY, moaC, moeA, and moa genes are similar to the other bacteria under Actinomycetes and all the biovar Ovis strains that lack the nitrate locus in their genomes. The narGHI complex in C. pseudotuberculosis biovar Equi shows conserved functional domains. narK and narT genes in C. pseudotuberculosis biovar Equi encode two putative nitrate and nitrite carriers, respectively. Finally, the nirK, nrfA, and nrfH encoding nitrite reductase enzymes are conserved in C. pseudotuberculosis. Therefore, the study of pathogen metabolism can contribute to the identification of pathogen virulence, discovery of drug targets, host response, and generating hypotheses for experimental investigation.

Competing Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors would like to acknowledge the help of all the team members and the financing agencies. CNPq (Conselho Nacional de Desenvolvimento Científico e Tecnológico, Brazil), CAPES (Coordenação de Aperfeiçoamento de Pessoal de Nível Superior, Brazil), and FAPEMIG (Fundação de Amparo à Pesquisa do Estado de Minas Gerais, Brazil). Andrey P. Lage was also supported by the Programa Pesquisador Mineiro (PPM) (00923-15) from FAPEMIG.

Supplementary Materials

Supplementary Figures 1 to 3 and Figure 7: Transmembrane helices prediction to NarI, NarK, NarT and NirK proteins, respectively, from C. pseudotuberculosis biovar Equi. Supplementary Figure 4: Transcription unit for the nrfHA genes into C. pseudotuberculosis biovar Equi. Supplementary Figure 5: Multiple sequence alignment of C. pseudotuberculosis biovar Equi nirK gene. Supplementary Figure 6: Amino acid sequence alignment of bacterial NirKs. Supplementary Figure 8: Comparative nitrogen metabolism between E.coli and C. pseudotuberculosis biovar Equi. In red pathway to E. coli, in green to C. pseudotuberculosis.  Supplementary Table 1:  Similarity sequence analysis of nirK gene of C. pseudotuberculosis and bacterial CuNiRs

  1. Supplementary Material