Genomic decoding of Theobroma grandiflorum (cupuassu) at chromosomal scale: evolutionary insights for horticultural innovation.

BACKGROUND: Theobroma grandiflorum (Malvaceae), known as cupuassu, is a tree indigenous to the Amazon basin, valued for its large fruits and seed pulp, contributing notably to the Amazonian bioeconomy. The seed pulp is utilized in desserts and beverages, and its seed butter is used in cosmetics. Here, we present the sequenced telomere-to-telomere genome of cupuassu, disclosing its genomic structure, evolutionary features, and phylogenetic relationships within the Malvaceae family. FINDINGS: The cupuassu genome spans 423 Mb, encodes 31,381 genes distributed in 10 chromosomes, and exhibits approximately 65% gene synteny with the Theobroma cacao genome, reflecting a conserved evolutionary history, albeit punctuated with unique genomic variations. The main changes are pronounced by bursts of long-terminal repeat retrotransposons at postspecies divergence, retrocopied and singleton genes, and gene families displaying distinctive patterns of expansion and contraction. Furthermore, positively selected genes are evident, particularly among retained and dispersed tandem and proximal duplicated genes associated with general fruit and seed traits and defense mechanisms, supporting the hypothesis of potential episodes of subfunctionalization and neofunctionalization following duplication, as well as impact from distinct domestication process. These genomic variations may underpin the differences observed in fruit and seed morphology, ripening, and disease resistance between cupuassu and the other Malvaceae species. CONCLUSIONS: The cupuassu genome offers a foundational resource for both breeding improvement and conservation biology, yielding insights into the evolution and diversity within the genus Theobroma.

Genomic analysis of two phlebotomine sand fly vectors of Leishmania from the New and Old World.

Phlebotomine sand flies are of global significance as important vectors of human disease, transmitting bacterial, viral, and protozoan pathogens, including the kinetoplastid parasites of the genus Leishmania, the causative agents of devastating diseases collectively termed leishmaniasis. More than 40 pathogenic Leishmania species are transmitted to humans by approximately 35 sand fly species in 98 countries with hundreds of millions of people at risk around the world. No approved efficacious vaccine exists for leishmaniasis and available therapeutic drugs are either toxic and/or expensive, or the parasites are becoming resistant to the more recently developed drugs. Therefore, sand fly and/or reservoir control are currently the most effective strategies to break transmission. To better understand the biology of sand flies, including the mechanisms involved in their vectorial capacity, insecticide resistance, and population structures we sequenced the genomes of two geographically widespread and important sand fly vector species: Phlebotomus papatasi, a vector of Leishmania parasites that cause cutaneous leishmaniasis, (distributed in Europe, the Middle East and North Africa) and Lutzomyia longipalpis, a vector of Leishmania parasites that cause visceral leishmaniasis (distributed across Central and South America). We categorized and curated genes involved in processes important to their roles as disease vectors, including chemosensation, blood feeding, circadian rhythm, immunity, and detoxification, as well as mobile genetic elements. We also defined gene orthology and observed micro-synteny among the genomes. Finally, we present the genetic diversity and population structure of these species in their respective geographical areas. These genomes will be a foundation on which to base future efforts to prevent vector-borne transmission of Leishmania parasites.

Comparative analyses of Theobroma cacao and T. grandiflorum mitogenomes reveal conserved gene content embedded within complex and plastic structures.

Unlike the chloroplast genomes (ptDNA), the plant mitochondrial genomes (mtDNA) are much more plastic in structure and size but maintain a conserved and essential gene set related to oxidative phosphorylation. Moreover, the plant mitochondrial genes and mtDNA are good markers for phylogenetic, evolutive, and comparative analyses. The two most known species in Theobroma L. (Malvaceae s.l.) genus are T. cacao, and T. grandiflorum. Besides the economic value, both species also show considerable biotechnology potential due to their other derived products, thus, aggregating additional economic value for the agroindustry. Here, we assembled and compared the mtDNA of Theobroma cacao and T. grandiflorum to generate a new genomics resource and unravel evolutionary trends. Graph-based analyses revealed that both mtDNA exhibit multiple alternative arrangements, confirming the dynamism commonly observed in plant mtDNA. The disentangled assembly graph revealed potential predominant circular molecules. The master circle molecules span 543,794 bp for T. cacao and 501,598 bp for T. grandiflorum, showing 98.9% of average sequence identity. Both mtDNA contains the same set of 39 plant mitochondrial genes, commonly found in other rosid mitogenomes. The main features are a duplicated copy of atp4, the absence of rpl6, rps2, rps8, and rps11, and the presence of two chimeric open-reading frames. Moreover, we detected few ptDNA integrations mainly represented by tRNAs, and no viral sequences were detected. Phylogenomics analyses indicate Theobroma spp. are nested in Malvaceae family. The main mtDNA differences are related to distinct structural rearrangements and exclusive regions associated with relics of Transposable Elements, supporting the hypothesis of dynamic mitochondrial genome maintenance and divergent evolutionary paths and pressures after species differentiation.

Assessment of Gene Flow to Wild Relatives and Nutritional Composition of Sugarcane in Brazil.

The commercial release of genetically modified organisms (GMO) requires a prior environmental and human/animal health risk assessment. In Brazil, the National Biotechnology Technical Commission (CTNBio) requires a survey of the area of natural occurrence of wild relatives of the GMO in the Brazilian ecosystems to evaluate the possibility of introgressive hybridization between sexually compatible species. Modern sugarcane cultivars, the focus of this study, derive from a series of hybridization and backcrossing events among Saccharum species. The so-called "Saccharum broad sense" group includes around 40 species from a few genera, including Erianthus, found in various tropical regions, particularly South-Eastern Asia. In Brazil, three native species, originally considered to belong to Erianthus, were reclassified as S. angustifolium (Nees) Trin., S. asperum (Nees) Steud., and S. villosum Steud., based on inflorescence morphology. Thus, we have investigated the potential occurrence of gene flow among the Brazilian Saccharum native species and commercial hybrids as a requisite for GMO commercial release. A comprehensive survey was carried out to map the occurrence of the three native Saccharum species in Brazil, concluding that they are sympatric with sugarcane cultivation only from around 14°S southwards, which precludes most Northeastern sugarcane-producing states from undergoing introgression. Based on phenology, we concluded that the Brazilian Saccharum species are unable to outcross naturally with commercial sugarcane since the overlap between the flowering periods of sugarcane and the native species is limited. A phylogenomic reconstruction based on the full plastid genome sequence showed that the three native Saccharum species are the taxa closest to sugarcane in Brazil, being closer than introduced Erianthus or Miscanthus. A 2-year study on eight nutritional composition traits of the 20 main sugarcane cultivars cultivated in Brazil was carried out in six environments. The minimum and maximum values obtained were, in percent: moisture (62.6-82.5); sucrose (9.65-21.76); crude fiber (8.06-21.03); FDN (7.20-20.68); FDA (4.55-16.90); lipids (0.06-1.59); ash (0.08-2.67); and crude protein (0.18-1.18). Besides a considerable amount of genetic variation and plastic responses, many instances of genotype-by-environment interaction were detected.

Genomic and Metabolomic Analyses of Natural Products in Nodularia spumigena Isolated from a Shrimp Culture Pond.

The bloom-forming cyanobacterium Nodularia spumigena CENA596 encodes the biosynthetic gene clusters (BGCs) of the known natural products nodularins, spumigins, anabaenopeptins/namalides, aeruginosins, mycosporin-like amino acids, and scytonemin, along with the terpenoid geosmin. Targeted metabolomics confirmed the production of these metabolic compounds, except for the alkaloid scytonemin. Genome mining of N. spumigena CENA596 and its three closely related Nodularia strains-two planktonic strains from the Baltic Sea and one benthic strain from Japanese marine sediment-revealed that the number of BGCs in planktonic strains was higher than in benthic one. Geosmin-a volatile compound with unpleasant taste and odor-was unique to the Brazilian strain CENA596. Automatic annotation of the genomes using subsystems technology revealed a related number of coding sequences and functional roles. Orthologs from the Nodularia genomes are involved in the primary and secondary metabolisms. Phylogenomic analysis of N. spumigena CENA596 based on 120 conserved protein sequences positioned this strain close to the Baltic Nodularia. Phylogeny of the 16S rRNA genes separated the Brazilian CENA596 strain from those of the Baltic Sea, despite their high sequence identities (99% identity, 100% coverage). The comparative analysis among planktic Nodularia strains showed that their genomes were considerably similar despite their geographically distant origin.

Metagenomic Approaches to Analyze Antimicrobial Resistance: An Overview.

Antimicrobial resistance is a major global public health problem, which develops when pathogens acquire antimicrobial resistance genes (ARGs), primarily through genetic recombination between commensal and pathogenic microbes. The resistome is a collection of all ARGs. In microorganisms, the primary method of ARG acquisition is horizontal gene transfer (HGT). Thus, understanding and identifying HGTs, can provide insight into the mechanisms of antimicrobial resistance transmission and dissemination. The use of high-throughput sequencing technologies has made the analysis of ARG sequences feasible and accessible. In particular, the metagenomic approach has facilitated the identification of community-based antimicrobial resistance. This approach is useful, as it allows access to the genomic data in an environmental sample without the need to isolate and culture microorganisms prior to analysis. Here, we aimed to reflect on the challenges of analyzing metagenomic data in the three main approaches for studying antimicrobial resistance: (i) analysis of microbial diversity, (ii) functional gene analysis, and (iii) searching the most complete and pertinent resistome databases.

Corrigendum: Genomic and Genotypic Characterization of Cylindrospermopsis raciborskii: Toward an Intraspecific Phylogenetic Evaluation by Comparative Genomics.

[This corrects the article on p. 306 in vol. 9, PMID: 29535689.].

Genomic and Genotypic Characterization of Cylindrospermopsis raciborskii: Toward an Intraspecific Phylogenetic Evaluation by Comparative Genomics.

Cylindrospermopsis raciborskii is a freshwater cyanobacterial species with increasing bloom reports worldwide that are likely due to factors related to climate change. In addition to the deleterious effects of blooms on aquatic ecosystems, the majority of ecotypes can synthesize toxic secondary metabolites causing public health issues. To overcome the harmful effects of C. raciborskii blooms, it is important to advance knowledge of diversity, genetic variation, and evolutionary processes within populations. An efficient approach to exploring this diversity and understanding the evolution of C. raciborskii is to use comparative genomics. Here, we report two new draft genomes of C. raciborskii (strains CENA302 and CENA303) from Brazilian isolates of different origins and explore their molecular diversity, phylogeny, and evolutionary diversification by comparing their genomes with sequences from other strains available in public databases. The results obtained by comparing seven C. raciborskii and the Raphidiopsis brookii D9 genomes revealed a set of conserved core genes and a variable set of accessory genes, such as those involved in the biosynthesis of natural products, heterocyte glycolipid formation, and nitrogen fixation. Gene cluster arrangements related to the biosynthesis of the antifungal cyclic glycosylated lipopeptide hassallidin were identified in four C. raciborskii genomes, including the non-nitrogen fixing strain CENA303. Shifts in gene clusters involved in toxin production according to geographic origins were observed, as well as a lack of nitrogen fixation (nif) and heterocyte glycolipid (hgl) gene clusters in some strains. Single gene phylogeny (16S rRNA sequences) was congruent with phylogeny based on 31 concatenated housekeeping protein sequences, and both analyses have shown, with high support values, that the species C. raciborskii is monophyletic. This comparative genomics study allowed a species-wide view of the biological diversity of C. raciborskii and in some cases linked genome differences to phenotype.

Quadruplex PCR assay for identification of Corynebacterium pseudotuberculosis differentiating biovar Ovis and Equi.

BACKGROUND: Corynebacterium pseudotuberculosis is classified into two biovars, nitrate-negative biovar Ovis which is the etiologic agent of caseous lymphadenitis in small ruminants and nitrate-positive biovar Equi, which causes abscesses and ulcerative lymphangitis in equines. The aim of this study was to develop a quadruplex PCR assay that would allow simultaneous detection and biovar-typing of C. pseudotuberculosis. METHODS: In the present study, genomes of C. pseudotuberculosis strains were used to identify the genes involved in the nitrate reduction pathway to improve a species identification three-primer multiplex PCR assay. The nitrate reductase gene (narG) was included in the PCR assay along with the 16S, rpoB and pld genes to enhance the diagnosis of the multiplex PCR at biovar level. RESULTS: A novel quadruplex PCR assay for C. pseudotuberculosis species and biovar identification was developed. The results of the quadruplex PCR of 348 strains, 346 previously well-characterized clinical isolates of C. pseudotuberculosis from different hosts (goats, sheep, horse, cattle, buffalo, llamas and humans), the vaccine strain 1002 and the type strain ATCC 19410T, were compared to the results of nitrate reductase identification by biochemical test. The McNemar's Chi-squared test used to compare the two methods used for C. pseudotuberculosis biovar identification showed no significant difference (P = 0.75) [95% CI for odds ratio (0.16-6.14)] between the quadruplex PCR and the nitrate biochemical test. Concordant results were observed for 97.13% (338 / 348) of the tested strains and the kappa value was 0.94 [95% CI (0.90-0.98)]. CONCLUSIONS: The ability of the quadruplex assay to discriminate between C. pseudotuberculosis biovar Ovis and Equi strains enhances its usefulness in the clinical microbiology laboratory.

Exploration of Nitrate Reductase Metabolic Pathway in Corynebacterium pseudotuberculosis.

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.

Draft Genome Assembly of the Bloom-Forming Cyanobacterium Nodularia spumigena Strain CENA596 in Shrimp Production Ponds.

We report here the draft genome assembly of the brackish cyanobacterium Nodularia spumigena strain CENA596 isolated from a shrimp production pond in Rio Grande do Sul, Brazil. The draft genome consists of 291 contigs with a total size of 5,189,679 bp. Secondary metabolite annotations resulted in several predicted gene clusters, including those responsible for encoding the hepatotoxin nodularin.

Draft Genome Sequences of Two Species of "Difficult-to-Identify" Human-Pathogenic Corynebacteria: Implications for Better Identification Tests.

Non-diphtheriae Corynebacterium species have been increasingly recognized as the causative agents of infections in humans. Differential identification of these bacteria in the clinical microbiology laboratory by the most commonly used biochemical tests is challenging, and normally requires additional molecular methods. Herein, we present the annotated draft genome sequences of two isolates of "difficult-to-identify" human-pathogenic corynebacterial species: C. xerosis and C. minutissimum. The genome sequences of ca. 2.7 Mbp, with a mean number of 2,580 protein encoding genes, were also compared with the publicly available genome sequences of strains of C. amycolatum and C. striatum. These results will aid the exploration of novel biochemical reactions to improve existing identification tests as well as the development of more accurate molecular identification methods through detection of species-specific target genes for isolate's identification or drug susceptibility profiling.

CMRegNet-An interspecies reference database for corynebacterial and mycobacterial regulatory networks.

BACKGROUND: Organisms utilize a multitude of mechanisms for responding to changing environmental conditions, maintaining their functional homeostasis and to overcome stress situations. One of the most important mechanisms is transcriptional gene regulation. In-depth study of the transcriptional gene regulatory network can lead to various practical applications, creating a greater understanding of how organisms control their cellular behavior. DESCRIPTION: In this work, we present a new database, CMRegNet for the gene regulatory networks of Corynebacterium glutamicum ATCC 13032 and Mycobacterium tuberculosis H37Rv. We furthermore transferred the known networks of these model organisms to 18 other non-model but phylogenetically close species (target organisms) of the CMNR group. In comparison to other network transfers, for the first time we utilized two model organisms resulting into a more diverse and complete network of the target organisms. CONCLUSION: CMRegNet provides easy access to a total of 3,103 known regulations in C. glutamicum ATCC 13032 and M. tuberculosis H37Rv and to 38,940 evolutionary conserved interactions for 18 non-model species of the CMNR group. This makes CMRegNet to date the most comprehensive database of regulatory interactions of CMNR bacteria. The content of CMRegNet is publicly available online via a web interface found at http://lgcm.icb.ufmg.br/cmregnet .

Genome Sequence of Lactococcus lactis subsp. lactis NCDO 2118, a GABA-Producing Strain.

Lactococcus lactis subsp. lactis NCDO 2118 is a nondairy lactic acid bacterium, a xylose fermenter, and a gamma-aminobutyric acid (GABA) producer isolated from frozen peas. Here, we report the complete genome sequence of L. lactis NCDO 2118, a strain with probiotic potential activity.

C. pseudotuberculosis Phop confers virulence and may be targeted by natural compounds.

The bacterial two-component system (TCS) regulates genes that are crucial for virulence in several pathogens. One of such TCS, the PhoPR system, consisting of a transmembrane sensory histidine kinase protein (PhoR) and an intracellular response regulator protein (PhoP), has been reported to have a major role in mycobacterial pathogenesis. We knocked out the phoP in C. pseudotuberculosis, the causal organism of caseous lymphadenitis (CLA), and using a combination of in vitro and in vivo mouse system, we showed for the first time, that the PhoP of C. pseudotuberculosis plays an important role in the virulence and pathogenicity of this bacterium. Furthermore, we modeled the PhoP of C. pseudotuberculosis and our docking results showed that several natural compounds including Rhein, an anthraquinone from Rheum undulatum, and some drug-like molecules may target PhoP to inhibit the TCS of C. pseudotuberculosis, and therefore may facilitate a remarkable attenuation of bacterial pathogenicity being the CLA. Experiments are currently underway to validate these in silico docking results.

Proteome scale comparative modeling for conserved drug and vaccine targets identification in Corynebacterium pseudotuberculosis.

Corynebacterium pseudotuberculosis (Cp) is a pathogenic bacterium that causes caseous lymphadenitis (CLA), ulcerative lymphangitis, mastitis, and edematous to a broad spectrum of hosts, including ruminants, thereby threatening economic and dairy industries worldwide. Currently there is no effective drug or vaccine available against Cp. To identify new targets, we adopted a novel integrative strategy, which began with the prediction of the modelome (tridimensional protein structures for the proteome of an organism, generated through comparative modeling) for 15 previously sequenced C. pseudotuberculosis strains. This pan-modelomics approach identified a set of 331 conserved proteins having 95-100% intra-species sequence similarity. Next, we combined subtractive proteomics and modelomics to reveal a set of 10 Cp proteins, which may be essential for the bacteria. Of these, 4 proteins (tcsR, mtrA, nrdI, and ispH) were essential and non-host homologs (considering man, horse, cow and sheep as hosts) and satisfied all criteria of being putative targets. Additionally, we subjected these 4 proteins to virtual screening of a drug-like compound library. In all cases, molecules predicted to form favorable interactions and which showed high complementarity to the target were found among the top ranking compounds. The remaining 6 essential proteins (adk, gapA, glyA, fumC, gnd, and aspA) have homologs in the host proteomes. Their active site cavities were compared to the respective cavities in host proteins. We propose that some of these proteins can be selectively targeted using structure-based drug design approaches (SBDD). Our results facilitate the selection of C. pseudotuberculosis putative proteins for developing broad-spectrum novel drugs and vaccines. A few of the targets identified here have been validated in other microorganisms, suggesting that our modelome strategy is effective and can also be applicable to other pathogens.

Exoproteome and secretome derived broad spectrum novel drug and vaccine candidates in Vibrio cholerae targeted by Piper betel derived compounds.

Vibrio cholerae is the causal organism of the cholera epidemic, which is mostly prevalent in developing and underdeveloped countries. However, incidences of cholera in developed countries are also alarming. Because of the emergence of new drug-resistant strains, even though several generic drugs and vaccines have been developed over time, Vibrio infections remain a global health problem that appeals for the development of novel drugs and vaccines against the pathogen. Here, applying comparative proteomic and reverse vaccinology approaches to the exoproteome and secretome of the pathogen, we have identified three candidate targets (ompU, uppP and yajC) for most of the pathogenic Vibrio strains. Two targets (uppP and yajC) are novel to Vibrio, and two targets (uppP and ompU) can be used to develop both drugs and vaccines (dual targets) against broad spectrum Vibrio serotypes. Using our novel computational approach, we have identified three peptide vaccine candidates that have high potential to induce both B- and T-cell-mediated immune responses from our identified two dual targets. These two targets were modeled and subjected to virtual screening against natural compounds derived from Piper betel. Seven compounds were identified first time from Piper betel to be highly effective to render the function of these targets to identify them as emerging potential drugs against Vibrio. Our preliminary validation suggests that these identified peptide vaccines and betel compounds are highly effective against Vibrio cholerae. Currently we are exhaustively validating these targets, candidate peptide vaccines, and betel derived lead compounds against a number of Vibrio species.

The pan-genome of the animal pathogen Corynebacterium pseudotuberculosis reveals differences in genome plasticity between the biovar ovis and equi strains.

Corynebacterium pseudotuberculosis is a facultative intracellular pathogen and the causative agent of several infectious and contagious chronic diseases, including caseous lymphadenitis, ulcerative lymphangitis, mastitis, and edematous skin disease, in a broad spectrum of hosts. In addition, Corynebacterium pseudotuberculosis infections pose a rising worldwide economic problem in ruminants. The complete genome sequences of 15 C. pseudotuberculosis strains isolated from different hosts and countries were comparatively analyzed using a pan-genomic strategy. Phylogenomic, pan-genomic, core genomic, and singleton analyses revealed close relationships among pathogenic corynebacteria, the clonal-like behavior of C. pseudotuberculosis and slow increases in the sizes of pan-genomes. According to extrapolations based on the pan-genomes, core genomes and singletons, the C. pseudotuberculosis biovar ovis shows a more clonal-like behavior than the C. pseudotuberculosis biovar equi. Most of the variable genes of the biovar ovis strains were acquired in a block through horizontal gene transfer and are highly conserved, whereas the biovar equi strains contain great variability, both intra- and inter-biovar, in the 16 detected pathogenicity islands (PAIs). With respect to the gene content of the PAIs, the most interesting finding is the high similarity of the pilus genes in the biovar ovis strains compared with the great variability of these genes in the biovar equi strains. Concluding, the polymerization of complete pilus structures in biovar ovis could be responsible for a remarkable ability of these strains to spread throughout host tissues and penetrate cells to live intracellularly, in contrast with the biovar equi, which rarely attacks visceral organs. Intracellularly, the biovar ovis strains are expected to have less contact with other organisms than the biovar equi strains, thereby explaining the significant clonal-like behavior of the biovar ovis strains.

Genome sequence of Corynebacterium pseudotuberculosis biovar equi strain 258 and prediction of antigenic targets to improve biotechnological vaccine production.

Corynebacterium pseudotuberculosis is the causative agent of several veterinary diseases in a broad range of economically important hosts, which can vary from caseous lymphadenitis in sheep and goats (biovar ovis) to ulcerative lymphangitis in cattle and horses (biovar equi). Existing vaccines against C. pseudotuberculosis are mainly intended for small ruminants and, even in these hosts, they still present remarkable limitations. In this study, we present the complete genome sequence of C. pseudotuberculosis biovar equi strain 258, isolated from a horse with ulcerative lymphangitis. The genome has a total size of 2,314,404 bp and contains 2088 predicted protein-coding regions. Using in silico analysis, eleven pathogenicity islands were detected in the genome sequence of C. pseudotuberculosis 258. The application of a reverse vaccinology strategy identified 49 putative antigenic proteins, which can be used as candidate vaccine targets in future works.

The Corynebacterium pseudotuberculosis in silico predicted pan-exoproteome.

BACKGROUND: Pan-genomic studies aim, for instance, at defining the core, dispensable and unique genes within a species. A pan-genomics study for vaccine design tries to assess the best candidates for a vaccine against a specific pathogen. In this context, rather than studying genes predicted to be exported in a single genome, with pan-genomics it is possible to study genes present in different strains within the same species, such as virulence factors. The target organism of this pan-genomic work here presented is Corynebacterium pseudotuberculosis, the etiologic agent of caseous lymphadenitis (CLA) in goat and sheep, which causes significant economic losses in those herds around the world. Currently, only a few antigens against CLA are known as being the basis of commercial and still ineffective vaccines. In this regard, the here presented work analyses, in silico, five C. pseudotuberculosis genomes and gathers data to predict common exported proteins in all five genomes. These candidates were also compared to two recent C. pseudotuberculosis in vitro exoproteome results. RESULTS: The complete genome of five C. pseudotuberculosis strains (1002, C231, I19, FRC41 and PAT10) were submitted to pan-genomics analysis, yielding 306, 59 and 12 gene sets, respectively, representing the core, dispensable and unique in silico predicted exported pan-genomes. These sets bear 150 genes classified as secreted (SEC) and 227 as potentially surface exposed (PSE). Our findings suggest that the main C. pseudotuberculosis in vitro exoproteome could be greater, appended by a fraction of the 35 proteins formerly predicted as making part of the variant in vitro exoproteome. These genomes were manually curated for correct methionine initiation and redeposited with a total of 1885 homogenized genes. CONCLUSIONS: The in silico prediction of exported proteins has allowed to define a list of putative vaccine candidate genes present in all five complete C. pseudotuberculosis genomes. Moreover, it has also been possible to define the in silico predicted dispensable and unique C. pseudotuberculosis exported proteins. These results provide in silico evidence to further guide experiments in the areas of vaccines, diagnosis and drugs. The work here presented is the first whole C. pseudotuberculosis in silico predicted pan-exoproteome completed till today.