- Extended genome report
- Open Access
First genome sequencing and comparative analyses of Corynebacterium pseudotuberculosis strains from Mexico
Standards in Genomic Sciences volume 13, Article number: 21 (2018)
Corynebacterium pseudotuberculosis is a pathogenic bacterium which has been rapidly spreading all over the world, causing economic losses in the agricultural sector and sporadically infecting humans. Six C. pseudotuberculosis strains were isolated from goats, sheep, and horses with distinct abscess locations. For the first time, Mexican genomes of this bacterium were sequenced and studied in silico. All strains were sequenced using Ion Personal Genome Machine sequencer, assembled using Newbler and SPAdes software. The automatic genome annotation was done using the software RAST and in-house scripts for transference, followed by manual curation using Artemis software and BLAST against NCBI and UniProt databases. The six genomes are publicly available in NCBI database. The analysis of nucleotide sequence similarity and the generated phylogenetic tree led to the observation that the Mexican strains are more similar between strains from the same host, but the genetic structure is probably more influenced by transportation of animals between farms than host preference. Also, a putative drug target was predicted and in silico analysis of 46 strains showed two gene clusters capable of differentiating the biovars equi and ovis: Restriction Modification system and CRISPR-Cas cluster.
Corynebacterium pseudotuberculosis is a Gram-positive bacterium that infects several different species of mammals. Strains of the biovar ovis infect sheep and goats, and strains of the biovar equi infect larger mammals such as horses, camels, and buffaloes. The manifestation of the infection depends on the host [1,2,3,4]. This bacterium causes significant economic loss to animal production all over the world due to reduced production of wool, milk and meat, carcass condemnation, as well as the death of infected animals [4,5,6]. C. pseudotuberculosis can also affect humans, causing distinct kinds of lymphadenitis. Contamination occurs through contact with infected animals and consumption of infected food [4, 5, 7].
This organism affects several countries such as Australia, Brazil, Canada, Egypt, Israel, New Zealand, South Africa, United Kingdom and United States [4, 8,9,10,11,12,13,14,15,16,17]. Cases in other countries such as Portugal , Mexico  and Equatorial Guinea  have been reported in the recent years. In the United States, C. pseudotuberculosis infections are reemerging and considered endemic , and the state with the highest number of cases of this bacterium was Texas, which borders Mexico . The spread of C. pseudotuberculosis to other countries brings out the importance of improving the understanding of this bacterium. In the present study, six Mexican C. pseudotuberculosis strains were investigated, two from the biovar equi and four from the biovar ovis. This is the first time that strains of this bacterium, isolated in Mexico, have been completely sequenced. Among those strains, these are the first isolates of the biovar equi coming from this country . The characterization of these strains is important for achieving a better understanding of this species, considering they can present relevant features not yet identified in other strains.
C. pseudotuberculosis is a pathogenic bacterium that belongs to the CMNR ( Corynebacterium , Mycobacterium , Nocardia , and Rhodococcus ) group. This group is characterized by high GC content (46–74%) and by the structure of the cell wall which is mainly composed of peptidoglycan, arabinogalactan and mycolic acids [4, 22]. C. pseudotuberculosis is placed in the phylum Actinobacteria , class Actinobacteria , order Actinomycetales , suborder Corynebacterineae and genus Corynebacterium [23,24,25,26,27,28,29,30]. The species is considered a facultative intracellular pathogen [4, 31] which is Gram-positive, pleomorphic, non-motile, non-sporulating, mesophilic and can survive both in the host and in the soil [25, 31,32,33,34,35]. Its strains are classified into two biovars, ovis and equi, according to its host preference and nitrate reduction capacity, which is identified through the presence or absence of the narG gene in a PCR Multiplex test . The biovar equi can reduce nitrate and affects mostly large ruminants. The biovar ovis is not able to reduce nitrate and affects mostly small ruminants . More information about classification, general features of this species and some details about the target strains are shown in Table 1 (Additional file 1).
Six C. pseudotuberculosis strains were isolated in Mexico from different hosts and biovars. The strain MEX1 was isolated from a retropharyngeal abscess in a goat. The strain MEX9 was isolated from a prescapular abscess in a goat. The strain MEX25 was isolated from a parotidean abscess in a sheep. The strain MEX29 was isolated from a retropharyngeal abscess in a sheep. These four strains presented negative result for the presence of the narG gene in the PCR multiplex test and were classified as belonging to the biovar ovis. All ovis strains were obtained from outbreaks occurred relatively close to Mexico City. MEX30 and MEX31 were isolated from abscesses in the pectoral muscles of two horses . These two strains were positive for the presence of the narG gene in PCR Multiplex. Consequently, they were classified as belonging to the biovar equi. Although both equi strains were obtained in the same city, they could be considered as isolated cases.
To verify the phylogenetic relationship of these strains to other strains of C. pseudotuberculosis , we generated a phylogenetic tree (Fig. 1) based on the core proteome and progressive refinement, using a bootstrap value of 100. The tree was generated using the PEPR software (https://github.com/enordber/pepr.git) with the Maximum-Likelihood method. The Mexican strains were clustered according to the respective biovars and host preferences, as shown in previous works) [1, 37].
MEX30 and MEX31 were isolated in Valparaiso, in the first reported case of infection of horses in Mexico . They clustered together probably because they came from the same source, that could be transported infected animals. Affected horses were identified in all regions of the US and the state of Texas, which borders Mexico, has the highest number of cases) [9, 21].
Ovis strains were isolated in Tlaxcala (MEX1) and Rio Frio de Juárez (MEX29), with a 50 Km distance from each other, and Guanajuato (MEX9 and MEX25), within a 400–450 Km distance from the two other isolation localities. However, the strains cluster by host rather than locality of isolation. MEX1 and MEX9 were isolated from goat and MEX25 and MEX29 were isolated from sheep. However, MEX25 and MEX29 (goat) clustered with isolates from lhama (USA) and cow (Israel), while MEX1 and MEX9 (sheep) clustered with isolates from goat and sheep (Brazil), all with a 100% bootstrap. Strains of Ovis biovar are more clonal but does not show the same degree of clustering by the host as Equi [1, 37]. Considering a maximum distance of 450 Km between localities of isolation, this genetic structure could better be explained by farming history than host preference. The goat and sheep farms could have different sources of Ovis strains. Transportation of infected animals and further contact and transmission of the disease probably occurred between farms of the same host species [38,39,40].
Genome sequencing information
Genome project history
The present project is a collaboration between the National Autonomous University of Mexico (UNAM), Mexico City, Mexico, and the Federal University of Minas Gerais (UFMG), Belo Horizonte, Minas Gerais, Brazil. The six C. pseudotuberculosis strains were isolated by UNAM researchers. Sequencing was performed at the National Reference Laboratory for Aquatic Animal Diseases (AQUACEN), and the two processes of assembly and annotation were performed at the Laboratory of Cellular and Molecular Genetics (LGCM), both laboratories located at UFMG. All genomes are complete and available at the National Center for Biotechnology Information (NCBI). This information is shown in Table 2 and conforms with MIGS recommendations . As mentioned above, the present study presents the first sequencing of C. pseudotuberculosis , and the first isolation of the biovar equi, from Mexico. This data can provide new insights into the diagnosis and treatment of diseases caused by this organism.
Growth conditions and genomic DNA preparation
The samples used in the present study are in the sample collection of LGCM. All six strains were grown in a brain-heart-infusion media (BHI-HiMedia Laboratories Pvt. Ltd., India) with 1.5% of bacteriological agar and supplemented with 0.5% of Tween 80, at 37 °C for 72 h under rotation. Genomic DNA was extracted following the protocol of Pacheco et al. .
Genome sequencing and assembly
The first step in sequencing each genome was the library construction, following manufacturer’s recommendations (IonXpress™ Plus gDNA Fragment Library Preparation). This was performed in three steps: (i) DNA fragmentation using the Ion Shear™ Plus Reagents Kit, (ii) addition of adapters using Ion Xpress™ Barcode Adapters and (iii) library amplification using the Ion PGM™ Template OT2 200 kit (all kits from Thermo Fisher Scientific, USA). The resulting library was put on the semiconductor chip Ion 318 Chip Kit v2 (Thermo Fisher Scientific) and then into the sequencer Ion Personal Genome Machine™ (Thermo Fisher Scientific). The number of reads and the mean read length of MEX1, MEX9, MEX25, MEX29, MEX30 and MEX31 strains are respectively: 1,100,551 and 244; 1,496,261 and 201; 1,117,243 and 206; 1,371,907 and 230; 1,127,325 and 186; and, 1,262,316 and 230.
The assembly process was managed using SIMBA software . The quality assessment of the reads was performed using FastQC software . The assemblies were performed using SPAdes version 3.6  on MEX1 and MEX31; and, Newbler version 2.9 (Roche, USA) on MEX9, MEX25, MEX29, and MEX30. This produced the following contigs: 6 on MEX1, 7 on MEX9, 7 on MEX25, 9 on MEX29, 33 on MEX30 and 13 on MEX31. The N50 s were: 543,202 on MEX1, 372,309 on MEX9, 543,326 on MEX25, 367,275 on MEX29, 103,276 on MEX30 and 535,978 on MEX31. The QUAST software  was used to evaluate the quality of the assemblies for all strains. The scaffolds were constructed using CONTIGuator software version 2.0  with C. pseudotuberculosis strain 29,156 (CP010795.1) as a reference to MEX9, MEX25 and MEX29, C. pseudotuberculosis strain MEX9 as a reference to MEX1, C. pseudotuberculosis strain 316 (CP003077.1) as a reference to MEX30 and C. pseudotuberculosis strain E19 (CP012136.1) as a reference to MEX31. Gap closure was performed using CLC Genomics Workbench 7 (Qiagen, USA). This process resulted in six complete genome sequences.
Genome annotation was performed in two steps: automatic annotation and manual curation. The RAST  and tRNAscan-SE  software were used in the automated annotation. An in-house script was also employed to transfer the annotation from a reference genome. The Artemis software version 16.0.0 , the UniProt  and the National Center for Biotechnology Information (NCBI) databases  were used in the manual curation. Putative frameshifts were analyzed using CLC Genomics Workbench 7 (Qiagen, USA) and fixed whenever possible.
Genome sizes of the respective strains are: 2,337,090 bp (base pairs) on MEX1, 2,337,578 bp on MEX9, 2,337,529 bp on MEX25, 2,337,866 bp on MEX29, 2,368,140 bp on MEX30 and 2,367,880 bp on MEX31. The respective percentages of the predicted coding regions are: 86.16% on MEX1, 86.33% on MEX9, 85.94% on MEX25, 86.66% on MEX29, 83.06% on MEX30 and 86.64% on MEX31. These genome sizes and the G + C content (~ 52%) are consistent with other C. pseudotuberculosis studies [2, 6, 52]. There are 64 predicted RNA genes in strains of the biovar ovis (MEX1, MEX9, MEX25 and MEX29) and 66 from the biovar equi (MEX30 and MEX31). The numbers (and percentages) of predicted protein coding genes and pseudogenes of MEX1, MEX9, MEX25, MEX29, MEX30 and MEX31 strains are, respectively: 2021 (94.22%) and 60 (2.80%); 2025 (94.36%) and 57 (2.66%); 2016 (94.07%) and 63 (2.94%); 2032 (94.73%) and 49 (2.28%); 2008 (91.77%) and 114 (5.21%); and 2058 (94.32%) and 61 (2.80%). Table 3 shows detailed information about properties and statistics of these genomes. The number of genes associated with general COG functional categories [53, 54] was generated with the in-house script Blast Cog (https://github.com/aquacen/blast_cog) and are summarized in Table 4. The circular maps of C. pseudotuberculosis MEX1 and MEX30 strains in comparison with the other strains of the present study are shown in Figs. 2 and 3, respectively.
Insights from the genome sequence
The nucleotide sequences, analyzed using the Gegenees software version 2.1 , show high similarity (> 92%) between the strains. Higher similarity (> = 99.7%) within strains belonging to the same biovar was found (Fig. 4). This is consistent with a previous study , using 15 strains of C. pseudotuberculosis , that shows similarity greater than 99% within the biovar ovis strains and at least 95% of sequencing similarity within the biovar equi strains. Moreover, the sequencing similarity among strains isolated from the same host is higher than the similarity among strains isolated from different hosts (Figs. 1 and 4).
Traditionally, the two biovars are differentiated using a nitrate reduction test, in which equi is positive, and ovis is negative . Figure 3 highlights the cluster of genes related to nitrate reduction in Mexican equi strains with the black rectangle. The Protein Family Sorter tool  was used to search for genes or clusters of genes that may be used to differentiate the biovars. Within the six genomes of the present study, we found the cluster of genes that is related to proteins of type III restriction-modification (RM) systems [58, 59] exclusively in the biovar ovis (highlighted in blue in Fig. 2). A cluster of genes related to the proteins of Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR-Cas) systems, probably belonging to type I-E , was found exclusively in the biovar equi (highlighted in blue in Fig. 3). Both systems work as protection systems, defending the bacteria against exogenous DNA. We analyzed 40 other sequenced strains of C. pseudotuberculosis to confirm these results in other strains. The same pattern was observed.
RM systems have two main components, a DNA methyltransferase, and a restriction endonuclease. The first one methylates the DNA in possible cleavage sites; the second one is responsible for the cleavage of DNA from external sources . A good review of RM systems can be found in . CRISPR-Cas systems are adaptive immune systems in bacteria and archaea. They use a complex of proteins known as Cas that are responsible for acquiring new, short sequences of external sources (exogenous genetic elements). These short sequences are incorporated into the bacterial chromosome and are called CRISPRs. The CRISPRs are transcribed into small RNAs that guide the Cas proteins to recognize and cleave foreign DNA, protecting the bacterial genome . Reviews of CRISPR-Cas systems can be found in [63,64,65].
Possible new drug targets were predicted using the Specialty Genes Search from the Pathosystems Resource Integration Center (PATRIC) bioinformatics resource center . The result shows a new putative target, the gene nrdF2, for five of the six strains used in the present study. In the C. pseudotuberculosis MEX30 strain, this gene is annotated as a pseudogene, which can explain why it was not considered a putative target. The product of this gene is the small subunit of ribonucleotide reductase (RNR) which is involved in dNTP (deoxynucleotide triphosphate) synthesis that reduces ribonucleotides to nucleotides. The RNRs can be classified into three classes (I, II and III). Class I is oxygen dependent and has two subclasses (Ia and Ib). Class Ia is coded by nrdA and nrdB genes; class Ib is coded by nrdE and nrdF. Therefore, the RNR found in the biovar ovis strains belongs to class Ib . Previous studies [68,69,70] show the importance of this gene for growth under normal conditions (in vitro) in Mycobacterium tuberculosis , Corynebacterium ammoniagenes and Corynebacterium glutamicum . Additionally, other studies have pointed to this gene as a potential target of M. tuberculosis vaccine [70,71,72].
In the present study, we investigated six strains of C. pseudotuberculosis from different hosts and their sequenced genomes, the first whole-genome investigation of this organism from Mexico. The phylogenomic analysis suggested that the genetic structure of Ovis is more influenced by animal transportation than host preference. An in silico analysis of protein families showed two important clusters that may differentiate the biovars equi and ovis. Also, the present work identified a new putative drug target against C. pseudotuberculosis , the gene nrdF2, which has been previously described as a potential vaccine target [70,71,72]. Further in silico and in vitro analyses are required to validate these findings. Those results could provide a better understanding of this organism and its mechanisms of virulence and pathogenesis, as well as develop new diagnoses, vaccines, and treatments.
National Reference Laboratory for Aquatic Animal Diseases
Brain heart infusion
Clustered regularly interspaced short palindromic repeats
Laboratory of Cellular and Molecular Genetics
Pathosystems Resource Integration Center
Phylogenomic Estimation with Progressive Refinement
- RM systems:
Federal University of Minas Gerais
National Autonomous University of Mexico
Soares SC, Silva A, Trost E, Blom J, Ramos R, Carneiro A, Ali A, et al. The pan-genome of the animal pathogen Corynebacterium pseudotuberculosis reveals differences in genome plasticity between the biovar ovis and equi strains. PLoS One. 2013;8(1):e53818.
Baraúna RA, Ramos RT, Veras AA, de Sá PH, Guimarães LC, das Graças DA, et al. Genomic analysis of four strains of Corynebacterium pseudotuberculosis bv. Equi isolated from horses showing distinct signs of infection. Stand Genomic Sci. 2017;12(1):16.
Barba M, Stewart AJ, Passler T, Wooldridge AA, Santen E, Chamorro MF, et al. Experimental transmission of Corynebacterium pseudotuberculosis biovar equi in horses by house flies. J Vet Intern Med. 2015;29(2):636–43.
Dorella FA, Pacheco LG, Oliveira SC, Miyoshi A, Azevedo V. Corynebacterium pseudotuberculosis: microbiology, biochemical properties, pathogenesis and molecular studies of virulence. Vet Res. 2006;37(2):201–18.
Trost E, Ott L, Schneider J, Schröder J, Jaenicke S, Goesmann A, et al. The complete genome sequence of Corynebacterium pseudotuberculosis FRC41 isolated from a 12-year-old girl with necrotizing lymphadenitis reveals insights into gene-regulatory networks contributing to virulence. BMC Genomics. 2010;11(1):728.
Almeida S, Tiwari S, Mariano D, Souza F, Jamal SB, Coimbra N, et al. The genome anatomy of Corynebacterium pseudotuberculosis VD57 a highly virulent strain causing Caseous lymphadenitis. Stand Genomic Sci. 2016;11(1):29.
Peel MM, Palmer GG, Stacpoole AM, Kerr TG. Human lymphadenitis due to Corynebacterium pseudotuberculosis: report of ten cases from Australia and review. Clin Infect Dis. 1997;24(2):185–91.
House RW, Schousboe M, Allen JP, Grant CC. Corynebacterium ovis (pseudo-tuberculosis) lymphadenitis in a sheep farmer: a new occupational disease in New Zealand. The New Zealand medical journal. 1986;99(809):659–62.
Kilcoyne I, Spier SJ, Carter CN, Smith JL, Swinford AK, Cohen ND. Frequency of Corynebacterium pseudotuberculosis infection in horses across the United States during a 10-year period. J Am Vet Med Assoc. 2014;245(3):309–14.
Seyffert N, Guimarães AS, Pacheco LG, Portela RW, Bastos BL, Dorella FA, et al. High seroprevalence of caseous lymphadenitis in Brazilian goat herds revealed by Corynebacterium pseudotuberculosis secreted proteins-based ELISA. Res Vet Sci. 2010;88(1):50–5.
Müller B, de Klerk-Lorist LM, Henton MM, Lane E, Parsons S, van Pittius NC, et al. Mixed infections of Corynebacterium pseudotuberculosis and non-tuberculous mycobacteria in south African antelopes presenting with tuberculosis-like lesions. Vet Microbiol. 2011;147(3):340–5.
Connor KM, Quirie MM, Baird G, Donachie W. Characterization of United Kingdom isolates of Corynebacterium pseudotuberculosis using pulsed-field gel electrophoresis. J Clin Microbiol. 2000;38(7):2633–7.
Ben SM, Ben MH, Benzarti M, Messadi L, Rejeb A, Amara A. Epidemiological and clinical studies of ovine caseous lymphadenitis. Arch Inst Pasteur Tunis. 2002;79(1–4):51–7.
Arsenault J, Girard C, Dubreuil P, Daignault D, Galarneau JR, Boisclair J, et al. Prevalence of and carcass condemnation from maedi–visna, paratuberculosis and caseous lymphadenitis in culled sheep from Quebec, Canada. Preventive veterinary medicine. 2003;59(1):67–81.
Paton MW, Walker SB, Rose IR, Watt GF. Prevalence of caseous lymphadenitis and usage of caseous lymphadenitis vaccines in sheep flocks. Aust Vet J. 2003;81(1–2):91–5.
Selim SA. Oedematous skin disease of buffalo in Egypt. Zoonoses Public Health. 2001;48(4):241–58.
Silva A, Schneider MP, Cerdeira L, Barbosa MS, Ramos RT, Carneiro AR, et al. Complete genome sequence of Corynebacterium pseudotuberculosis I19, a strain isolated from a cow in Israel with bovine mastitis. J Bacteriol. 2011;193(1):323–4.
Oliveira M, Barroco C, Mottola C, Santos R, Lemsaddek A, Tavares L, et al. First report of Corynebacterium pseudotuberculosis from caseous lymphadenitis lesions in black Alentejano pig (Sus scrofa domesticus). BMC Vet Res. 2014;10(1):218.
Muñoz-Bucio AV, Cortés-Pérez YA, Arellano-Reynoso B, Hernández-Gil M, Hernández-Castro R, Díaz-Aparicio E. Identification of Corynebacterium pseudotuberculosis isolated from muscular abscesses in two horses: first report in Mexico. Equine Veterinary Education. 2016;29:431-5.
Loureiro D, Portela RW, Sousa TJ, Rocha F, Pereira FL, Dorella FA, et al. Complete genome sequence of Corynebacterium pseudotuberculosis viscerotropic strain N1. Genome announcements. 2016;4(1):e01673–15.
Spier SJ, Azevedo V. Corynebacterium pseudotuberculosis infection in horses: increasing frequency and spread to new regions of North America. Equine Veterinary Education. 2016;29:436-9.
Tauch A, Sandbote J. The family Corynebacteriaceae. In: The prokaryotes: Springer Berlin Heidelberg; 2014. p. 239:277.
Woese CR, Kandler O, Wheelis ML. Towards a natural system of organisms: proposal for the domains archaea, Bacteria, and Eucarya. Proc Natl Acad Sci. 1990;87(12):4576–9.
Garrity GM, Holt JG. The road map to the manual. Bergey's Manual of Systematics of Archaea and Bacteria. 2001;1:119–66.
Stackebrandt E, Rainey FA, Ward-Rainey NL. Proposal for a new hierarchic classification system, Actinobacteria classis nov. Int J Syst Evol Microbiol. 1997;47(2):479–91.
Skerman VB, McGowan V, Sneath PH. Approved lists of bacterial names. Int J Syst Evol Microbiol. 1980;30(1):225–420.
Buchanan RE. Studies in the nomenclature and classification of the bacteria: II. The primary subdivisions of the schizomycetes. J Bacteriol. 1917;2(2):155.
Zhi XY, Li WJ, Stackebrandt E. An update of the structure and 16S rRNA gene sequence-based definition of higher ranks of the class Actinobacteria, with the proposal of two new suborders and four new families and emended descriptions of the existing higher taxa. Int J Syst Evol Microbiol. 2009;59(3):589–608.
Lehmann KB, Neumann R. Lehmann's Medizin, Handatlanten. X Atlas und Grundriss der Bakteriologie und Lehrbuch der speziellen bakteriologischen Diagnostik. 1907;2:270.
Bernard KA, Wiebe D, Burdz T, Reimer A, Ng B, Singh C, et al. Assignment of Brevibacterium stationis (ZoBell and Upham 1944) breed 1953 to the genus Corynebacterium, as Corynebacterium stationis comb. nov., and emended description of the genus Corynebacterium to include isolates that can alkalinize citrate. Int J Syst Evol Microbiol. 2010;60(4):874–9.
Moura-Costa LF, Bahia RC, Carminati R, Vale VL, Paule BJ, Portela RW, et al. Evaluation of the humoral and cellular immune response to different antigens of Corynebacterium pseudotuberculosis in Canindé goats and their potential protection against caseous lymphadenitis. Vet Immunol Immunopathol. 2008;126(1):131–41.
Pinto AC, de Sá PH, Ramos RT, Barbosa S, Barbosa HP, Ribeiro AC, et al. Differential transcriptional profile of Corynebacterium pseudotuberculosis in response to abiotic stresses. BMC Genomics. 2014;15(1):14.
Aleman M, Spier SJ, Wilson WD, Doherr M. Corynebacterium pseudotuberculosis infection in horses: 538 cases (1982-1993). J Am Vet Med Assoc. 1996;209(4):804–9.
Spier SJ, Toth B, Edman J, Quave A, Habasha F, Garrick M, et al. Survival of Corynebacterium pseudotuberculosis biovar equi in soil. Vet Rec. 2012;170:180-1. vetrec-2011.
Merchant IA, Packer RA. Veterinary bacteriology and virology. Ames: Iowa State University Press; 1967. p. 752.
Pacheco LG, Pena RR, Castro TL, Dorella FA, Bahia RC, Carminati R, et al. Multiplex PCR assay for identification of Corynebacterium pseudotuberculosis from pure cultures and for rapid detection of this pathogen in clinical samples. J Med Microbiol. 2007;56(4):480–6.
Viana MV, Figueiredo H, Ramos R, Guimarães LC, Pereira FL, Dorella FA, et al. Comparative genomic analysis between Corynebacterium pseudotuberculosis strains isolated from buffalo. PLoS One. 2017;12(4):e0176347.
Gebrekidan H, Nelson L, Smith G, Gasser RB, Jabbar A. An outbreak of oriental theileriosis in dairy cattle imported to Vietnam from Australia. Parasitology. 2017;144(6):738–46.
Keim PS, Wagner DM. Humans and evolutionary and ecological forces shaped the phylogeography of recently emerged diseases. Nat Rev Microbiol. 2009;7(11):813–21.
Skuce RA, Mallon TR, McCormick CM, McBride SH, Clarke G, Thompson A, Couzens C, Gordon AW, McDowell SW. Mycobacterium bovis genotypes in Northern Ireland: herd-level surveillance (2003 to 2008). The Veterinary record. 2010;167(18):684–9.
Field D, Garrity G, Gray T, Morrison N, Selengut J, Sterk P, et al. The minimum information about a genome sequence (MIGS) specification. Nat Biotechnol. 2008;26(5):541–7.
Mariano DC, Pereira FL, Aguiar EL, Oliveira LC, Benevides L, Guimarães LC, et al. SIMBA: a web tool for managing bacterial genome assembly generated by ion PGM sequencing technology. BMC bioinformatics. 2016;17(18):456.
Babraham Bioinformatics: FastQC. https://github.com/s-andrews/FastQC. Accessed 14 June 2017.
Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 2012;19(5):455–77.
Gurevich A, Saveliev V, Vyahhi N, Tesler G. QUAST: quality assessment tool for genome assemblies. Bioinformatics. 2013;29(8):1072–5.
Galardini M, Biondi EG, Bazzicalupo M, Mengoni A. CONTIGuator: a bacterial genomes finishing tool for structural insights on draft genomes. Source code for biology and medicine. 2011;6(1):11.
Aziz RK, Bartels D, Best AA, DeJongh M, Disz T, Edwards RA, et al. The RAST server: rapid annotations using subsystems technology. BMC Genomics. 2008;9(1):75.
Lowe TM, Eddy SR. tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997;25(5):955–64.
Rutherford K, Parkhill J, Crook J, Horsnell T, Rice P, Rajandream MA, et al. Artemis: sequence visualization and annotation. Bioinformatics. 2000;16(10):944–5.
NCBI. https://www.ncbi.nlm.nih.gov/. Accessed in 15 Jun 2017.
Hassan SS, Guimarães LC, Pereira UP, Islam A, Ali A, Bakhtiar SM, et al. Complete genome sequence of Corynebacterium pseudotuberculosis biovar ovis strain P54B96 isolated from antelope in South Africa obtained by rapid next generation sequencing technology. Stand Genomic Sci. 2012;7(2):189.
Galperin MY, Makarova KS, Wolf YI, Koonin EV. Expanded microbial genome coverage and improved protein family annotation in the COG database. Nucleic acids res. 2015;43:D261–9.
Tatusov RL. A genomic perspective on protein families. Science (80-. ). 1997;278:631–7.
Ågren J, Sundström A, Håfström T, Segerman B. Gegenees: fragmented alignment of multiple genomes for determining phylogenomic distances and genetic signatures unique for specified target groups. PLoS One. 2012;7(6):e39107.
Biberstein EL, Knight HD, Jang S. Two biotypes of Corynebacterium pseudotuberculosis. Vet Rec. 1971;89:691-2.
Davis JJ, Gerdes S, Olsen GJ, Olson R, Pusch GD, Shukla M, et al. PATtyFams: protein families for the microbial genomes in the PATRIC database. Front Microbiol. 2016;7:118.
Kwiatek A, Mrozek A, Bacal P, Piekarowicz A, Adamczyk-Popławska M. Type III methyltransferase M. NgoAX from Neisseria gonorrhoeae FA1090 regulates biofilm formation and interactions with human cells. Front Microbiol. 2015;6:1426.
Tan A, Hill DM, Harrison OB, Srikhanta YN, Jennings MP, Maiden MC, et al. Distribution of the type III DNA methyltransferases modA, modB and modD among Neisseria meningitidis genotypes: implications for gene regulation and virulence. Sci Rep. 2016;6:21015.
Sangal V, Fineran PC, Hoskisson PA. Novel configurations of type I and II CRISPR–Cas systems in Corynebacterium diphtheriae. Microbiology. 2013;159(10):2118–26.
Blow MJ, Clark TA, Daum CG, Deutschbauer AM, Fomenkov A, Fries R, et al. The epigenomic landscape of prokaryotes. PLoS Genet. 2016;12(2):e1005854.
Suzuki H. Host-mimicking strategies in DNA methylation for improved bacterial transformation. In: Methylation-from DNA, RNA and histones to diseases and treatment. London: InTech Open Access Publisher; 2012.
Bhaya D, Davison M, Barrangou R. CRISPR-Cas systems in bacteria and archaea: versatile small RNAs for adaptive defense and regulation. Annu Rev Genet. 2011;45:273–97.
Makarova KS, Haft DH, Barrangou R, Brouns SJ, Charpentier E, Horvath P, et al. Evolution and classification of the CRISPR–Cas systems. Nat Rev Microbiol. 2011;9(6):467–77.
Makarova KS, Wolf YI, Alkhnbashi OS, Costa F, Shah SA, Saunders SJ, et al. An updated evolutionary classification of CRISPR-Cas systems. Nat Rev Microbiol. 2015;13(11):722–36.
Mao C, Abraham D, Wattam AR, Wilson MJ, Shukla M, Yoo HS, et al. Curation, integration and visualization of bacterial virulence factors in PATRIC. Bioinformatics. 2014;31(2):252–8.
Kolberg M, Strand KR, Graff P, Andersson KK. Structure, function, and mechanism of ribonucleotide reductases. Biochimica et Biophysica Acta (BBA)-proteins and. Proteomics. 2004;1699(1):1–34.
Oehlmann W, Auling G. Ribonucleotide reductase (RNR) of Corynebacterium glutamicum ATCC 13032–genetic characterization of a second class IV enzyme. Microbiology. 1999;145(7):1595–604.
Torrents E, Roca I, Gibert I. Corynebacterium ammoniagenes class Ib ribonucleotide reductase: transcriptional regulation of an atypical genomic organization in the nrd cluster. Microbiology. 2003;149(4):1011–20.
Mowa MB. Function and expression of class I ribonucleotide reductase small subunit-encoding genes in Mycobacterium tuberculosis and Mycobacterium smegmatis (doctoral dissertation). 2009.
Yang F, Curran SC, Li LS, Avarbock D, Graf JD, Chua MM, et al. Characterization of two genes encoding the Mycobacterium tuberculosis ribonucleotide reductase small subunit. J Bacteriol. 1997;179(20):6408–15.
Nurbo J, Roos AK, Muthas D, Wahlström E, Ericsson DJ, Lundstedt T, et al. Design, synthesis and evaluation of peptide inhibitors of Mycobacterium tuberculosis ribonucleotide reductase. J Pept Sci. 2007;13(12):822–32.
Buxton A, Corynebacterium FG. Animal Microbiology. Edinburgh: Blackwell Scientific Publications; 1977. p. 177.
Brown CC, Olander HJ, Alves SF. Synergistic hemolysis-inhibition titers associated with caseous lymphadenitis in a slaughterhouse survey of goats and sheep in northeastern Brazil. Can J Vet Res. 1987;51(1):46.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. Nat Genet. 2000;25(1):25–9.
The authors thank the Brazilian funding agencies, CNPq, CAPES, and FAPEMIG, for providing financial support to the National Institute of Science and Technology in Theranostics and Nanobiotechnology – INCT-TeraNano (CNPq/CAPES/FAPEMIG, Grant numbers CNPq-465669/2014-0 and FAPEMIG-CBB-APQ-03613-17). Debmalya Barh acknowledges the “TWAS-CNPq Postdoctoral Fellowship Programme” for granting a fellowship for postdoctoral studies.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.