Spathaspora marinasilvae sp. nov., a xylose-fermenting yeast isolated from galleries of passalid beetles and rotting wood in the Amazonian rainforest biome.

Four yeast isolates were obtained from rotting wood and galleries of passalid beetles collected in different sites of the Brazilian Amazonian Rainforest in Brazil. This yeast produces unconjugated allantoid asci each with a single elongated ascospore with curved ends. Sequence analysis of the internal transcribed spacer-5.8 S region and the D1/D2 domains of the large subunit ribosomal RNA (rRNA) gene showed that the isolates represent a novel species of the genus Spathaspora. The novel species is phylogenetically related to a subclade containing Spathaspora arborariae and Spathaspora suhii. Phylogenomic analysis based on 1884 single-copy orthologs for a set of Spathaspora species whose whole genome sequences are available confirmed that the novel species represented by strain UFMG-CM-Y285 is phylogenetically close to Sp. arborariae. The name Spathaspora marinasilvae sp. nov. is proposed to accommodate the novel species. The holotype of Sp. marinasilvae is CBS 13467 T (MycoBank 852799). The novel species was able to accumulate xylitol and produce ethanol from d-xylose, a trait of biotechnological interest common to several species of the genus Spathaspora.

Vitamin B12 attenuates leukocyte inflammatory signature in COVID-19 via methyl-dependent changes in epigenetic markings.

COVID-19 induces chromatin remodeling in host immune cells, and it had previously been shown that vitamin B12 downregulates some inflammatory genes via methyl-dependent epigenetic mechanisms. In this work, whole blood cultures from moderate or severe COVID-19 patients were used to assess the potential of B12 as adjuvant drug. The vitamin normalized the expression of a panel of inflammatory genes still dysregulated in the leukocytes despite glucocorticoid therapy during hospitalization. B12 also increased the flux of the sulfur amino acid pathway, that regulates the bioavailability of methyl. Accordingly, B12-induced downregulation of CCL3 strongly and negatively correlated with the hypermethylation of CpGs in its regulatory regions. Transcriptome analysis revealed that B12 attenuates the effects of COVID-19 on most inflammation-related pathways affected by the disease. As far as we are aware, this is the first study to demonstrate that pharmacological modulation of epigenetic markings in leukocytes favorably regulates central components of COVID-19 physiopathology.

Metformin Treatment Modulates Long Non-Coding RNA Isoforms Expression in Human Cells.

Long noncoding RNAs (lncRNAs) undergo splicing and have multiple transcribed isoforms. Nevertheless, for lncRNAs, as well as for mRNA, measurements of expression are routinely performed only at the gene level. Metformin is the first-line oral therapy for type 2 diabetes mellitus and other metabolic diseases. However, its mechanism of action remains not thoroughly explained. Transcriptomic analyses using metformin in different cell types reveal that only protein-coding genes are considered. We aimed to characterize lncRNA isoforms that were differentially affected by metformin treatment on multiple human cell types (three cancer, two non-cancer) and to provide insights into the lncRNA regulation by this drug. We selected six series to perform a differential expression (DE) isoform analysis. We also inferred the biological roles for lncRNA DE isoforms using in silico tools. We found the same isoform of an lncRNA (AC016831.6-205) highly expressed in all six metformin series, which has a second exon putatively coding for a peptide with relevance to the drug action. Moreover, the other two lncRNA isoforms (ZBED5-AS1-207 and AC125807.2-201) may also behave as cis-regulatory elements to the expression of transcripts in their vicinity. Our results strongly reinforce the importance of considering DE isoforms of lncRNA for understanding metformin mechanisms at the molecular level.

Yeast communities associated with cacti in Brazil and the description of Kluyveromyces starmeri sp. nov. based on phylogenomic analyses.

Yeast communities associated with cacti were studied in three ecosystems of Southeast, Central and North Brazil. A total of 473 yeast strains belonging to 72 species were isolated from 190 samples collected. Cactophilic yeast species were prevalent in necrotic tissues, flowers, fruits and insects of cacti collected in Southeast and North Brazil. Pichia cactophila, Candida sonorensis and species of the Sporopachydermia complex were the most prevalent cactophilic species in Southeast and Central regions. Kodamaea nitidulidarum, Candida restingae and Wickerhamiella cacticola were frequently associated with cactus flowers and fruits. The diversity of yeasts associated with the substrates studied was high. Twenty-one novel species were found. One is described here as Kluyveromyces starmeri sp. nov. based on 21 isolates obtained from necrotic tissues, flowers, fruits and associated insects of the columnar cacti Cereus saddianus, Micranthocereus dolichospermaticus and Pilosocereus arrabidae in two different ecosystems in Brazil. Phylogenetic analyses of sequences encoding the gene of the small subunit (SSU) rRNA gene, the internal transcribed spacer, the 5.8S rRNA gene and the D1/D2 domains of the large subunit (LSU) rRNA showed that the species is related to Kluyveromyces dobzhanskii, Kluyveromyces lactis and Kluyveromyces marxianus. Phylogenomic analyses based on 1264 conserved genes shared among the new species and 19 other members of the Saccharomycetaceae confirmed this phylogenetic relationship. The holotype is K. starmeri sp. nov. CBS 16103T (=UFMG-CM-Y3682T ). The Mycobank number is MB 836817.

Noncoding SNPs associated with increased GDF15 levels located in a metformin-activated enhancer region upstream of GDF15.

Aim: GDF15 levels are a biomarker for metformin use. We performed the functional annotation of noncoding genome-wide association study (GWAS) SNPs for GDF15 levels and the Genotype-Tissue Expression (GTEx)-expression quantitative trait loci (eQTLs) for GDF15 expression within metformin-activated enhancers around GDF15. Materials & methods: These enhancers were identified using chromatin immunoprecipitation followed by sequencing data for active (H3K27ac) and silenced (H3K27me3) histone marks on human hepatocytes treated with metformin, Encyclopedia of DNA Elements data and cis-regulatory elements assignment tools. Results: The GWAS lead SNP rs888663, the SNP rs62122429 associated with GDF15 levels in the Outcome Reduction with Initial Glargine Intervention trial, and the GTEx-expression quantitative trait locus rs4808791 for GDF15 expression in whole blood are located in a metformin-activated enhancer upstream of GDF15 and tightly linked in Europeans and East Asians. Conclusion: Noncoding variation within a metformin-activated enhancer may increase GDF15 expression and help to predict GDF15 levels.

Scheffersomyces stambukii f.a., sp. nov., a d-xylose-fermenting species isolated from rotting wood.

Two isolates representing a new species of Scheffersomyces were isolated from rotting wood samples collected in an Amazonian forest ecosystem in Brazil. Analysis of the sequences of the D1/D2 domains showed that this new species is phylogenetically related to Scheffersomyces NYMU 15730, a species without a formal description, and the two are in an early emerging position with respect to the xylose-fermenting subclade containing Scheffersomyces titanus and Scheffersomyces stipitis. Phylogenomic analyses using 474 orthologous genes placed the new species in an intermediary position between Scheffersomyces species and the larger genus Spathaspora and the Candida albicans/Lodderomyces clade. The novel species, Scheffersomyces stambukii f.a., sp. nov., is proposed to accommodate these isolates. The type strain of Scheffersomyces stambukii sp. nov. is UFMG-CM-Y427T (=CBS 14217T). The MycoBank number is MB 824093. In addition, we studied the xylose metabolism of this new species.

Spathaspora boniae sp. nov., a D-xylose-fermenting species in the Candida albicans/Lodderomyces clade.

Two yeast isolates producing asci-containing elongate ascospores with curved ends typical of the genus Spathaspora were isolated from rotting wood samples collected in an Atlantic rainforest ecosystem in Brazil. Phylogenetic analysis of the LSU rRNA gene D1/D2 domain sequences demonstrated that the strains represent a new species and placed it next to Candida blackwellae, in a clade that also contains Candida albicans and Candida dubliniensis. Other sequences of the ribosomal gene cluster supported same placementin the same clade, and a phylogenomic analysis placed this new species in an early emerging position relative to the larger C. albicans/Lodderomyces clade. One interpretation is that the genus Spathaspora is, in fact, paraphyletic. In conformity with this view, we propose the novel species Spathaspora boniae sp. nov. to accommodate the isolates. The type strain of Spathaspora boniae sp. nov. is UFMG-CM-Y306T (=CBS 13262T). The MycoBank number is MB 821297. A detailed analysis of xylose metabolism was conducted for the new species.

Draft Genome Sequence of Metschnikowia australis Strain UFMG-CM-Y6158, an Extremophile Marine Yeast Endemic to Antarctica.

Here we report the draft genome sequence of Metschnikowia australis strain UFMG-CM-Y6158, a yeast endemic to Antarctica. We isolated the strain from the marine seaweed Acrosiphonia arcta (Chlorophyta). The genome is 14.3 Mb long and contains 4,442 predicted protein-coding genes.

Draft genome sequence of Sugiyamaella xylanicola UFMG-CM-Y1884T, a xylan-degrading yeast species isolated from rotting wood samples in Brazil.

We present the draft genome sequence of the type strain of the yeast Sugiyamaella xylanicola UFMG-CM-Y1884T (= UFMG-CA-32.1T = CBS 12683T), a xylan-degrading species capable of fermenting d-xylose to ethanol. The assembled genome has a size of ~ 13.7 Mb and a GC content of 33.8% and contains 5971 protein-coding genes. We identified 15 genes with significant similarity to the d-xylose reductase gene from several other fungal species. The draft genome assembled from whole-genome shotgun sequencing of the yeast Sugiyamaella xylanicola UFMG-CM-Y1884T (= UFMG-CA-32.1T = CBS 12683T) has been deposited at DDBJ/ENA/GenBank under the accession number MQSX00000000 under version MQSX01000000.

Effect of ionizing radiation exposure on Trypanosoma cruzi ubiquitin-proteasome system.

In recent years, proteasome involvement in the damage response induced by ionizing radiation (IR) became evident. However, whether proteasome plays a direct or indirect role in IR-induced damage response still unclear. Trypanosoma cruzi is a human parasite capable of remarkable high tolerance to IR, suggesting a highly efficient damage response system. Here, we investigate the role of T. cruzi proteasome in the damage response induced by IR. We exposed epimastigotes to high doses of gamma ray and we analyzed the expression and subcellular localization of several components of the ubiquitin-proteasome system. We show that proteasome inhibition increases IR-induced cell growth arrest and proteasome-mediated proteolysis is altered after parasite exposure. We observed nuclear accumulation of 19S and 20S proteasome subunits in response to IR treatments. Intriguingly, the dynamic of 19S particle nuclear accumulation was more similar to the dynamic observed for Rad51 nuclear translocation than the observed for 20S. In the other hand, 20S increase and nuclear translocation could be related with an increase of its regulator PA26 and high levels of proteasome-mediated proteolysis in vitro. The intersection between the opposed peaks of 19S and 20S protein levels was marked by nuclear accumulation of both 20S and 19S together with Ubiquitin, suggesting a role of ubiquitin-proteasome system in the nuclear protein turnover at the time. Our results revealed the importance of proteasome-mediated proteolysis in T. cruzi IR-induced damage response suggesting that proteasome is also involved in T. cruzi IR tolerance. Moreover, our data support the possible direct/signaling role of 19S in DNA damage repair. Based on these results, we speculate that spatial and temporal differences between the 19S particle and 20S proteasome controls proteasome multiple roles in IR damage response.

A single FTO gene variant rs9939609 is associated with body weight evolution in a multiethnic extremely obese population that underwent bariatric surgery.

OBJECTIVES: The rs9939609 single nucleotide polymorphism (SNP) in the fat mass and obesity-associated (FTO) gene is involved in obesity. Few studies have been conducted on patients who underwent bariatric surgery. The aim of this study was to evaluate the influence of FTO SNPs on body weight, body composition, and weight regain during a 60-mo follow-up period after bariatric surgery. METHODS: The rs9939609 was genotyped in 146 individuals using a real-time polymerase chain reaction TaqMan assay. Data for lifestyle, comorbidities, body weight, body mass index (BMI), excess weight loss (EWL), and body composition were obtained before and 6, 12, 18, 24, 36, 48, and 60 mo after surgery. Data were analyzed by comparing two groups of patients according to rs9939609 FTO gene polymorphism. Mixed-regression models were constructed to evaluate the dynamics of body weight, BMI, and EWL over time in female patients. RESULTS: No differences were observed between the groups during the first 24 mo after surgery. After 36, 48, and 60 mo, body weight, fat mass, and BMI were higher, whereas fat-free mass and EWL were lower in the FTO-SNP patient group. Weight regain was more frequent and occurred sooner in the FTO-SNP group. CONCLUSIONS: There is a different evolution of weight loss in obese carriers of the FTO gene variant rs9939609 after bariatric surgery. However, this pattern was evident at only 2 y postbariatric surgery, inducing a lower proportion of surgery success and a greater and earlier weight regain.

Evaluation of the Schistosoma mansoni Y-box-binding protein (SMYB1) potential as a vaccine candidate against schistosomiasis.

Schistosomiasis is a neglected tropical disease, and after malaria, is the second most important tropical disease in public health. A vaccine that reduces parasitemia is desirable to achieve mass treatment with a low cost. Although potential antigens have been identified and tested in clinical trials, no effective vaccine against schistosomiasis is available. Y-box-binding proteins (YBPs) regulate gene expression and participate in a variety of cellular processes, including transcriptional and translational regulation, DNA repair, cellular proliferation, drug resistance, and stress responses. The Schistosoma mansoni ortholog of the human YB-1, SMYB1, is expressed in all stages of the parasite life cycle. Although SMYB1 binds to DNA or RNA oligonucleotides, immunohistochemistry assays demonstrated that it is primarily localized in the cytoplasm of parasite cells. In addition, SMYB1 interacts with a protein involved in mRNA processing, suggesting that SMYB1 functions in the turnover, transport, and/or stabilization of RNA molecules during post-transcriptional gene regulation. Here we report the potential of SMYB1 as a vaccine candidate. We demonstrate that recombinant SMYB1 stimulates the production of high levels of specific IgG1 antibodies in a mouse model. The observed levels of specific IgG1 and IgG2a antibodies indicate an actual protection against cercariae challenge. Animals immunized with rSMYB1 exhibited a 26% reduction in adult worm burden and a 28% reduction in eggs retained in the liver. Although proteins from the worm tegument are considered optimal targets for vaccine development, this study demonstrates that unexposed cytoplasmic proteins can reduce the load of intestinal worms and the number of eggs retained in the liver.

Draft Genome Sequence of the D-Xylose-Fermenting Yeast Spathaspora arborariae UFMG-HM19.1AT.

The draft genome sequence of the yeast Spathaspora arborariae UFMG-HM19.1A(T) (CBS 11463 = NRRL Y-48658) is presented here. The sequenced genome size is 12.7 Mb, consisting of 41 scaffolds containing a total of 5,625 predicted open reading frames, including many genes encoding enzymes and transporters involved in d-xylose fermentation.

The spliced leader trans-splicing mechanism in different organisms: molecular details and possible biological roles.

THE SPLICED LEADER (SL) IS A GENE THAT GENERATES A FUNCTIONAL NCRNA THAT IS COMPOSED OF TWO REGIONS: an intronic region of unknown function (SLi) and an exonic region (SLe), which is transferred to the 5' end of independent transcripts yielding mature mRNAs, in a process known as spliced leader trans-splicing (SLTS). The best described function for SLTS is to solve polycistronic transcripts into monocistronic units, specifically in Trypanosomatids. In other metazoans, it is speculated that the SLe addition could lead to increased mRNA stability, differential recruitment of the translational machinery, modification of the 5' region or a combination of these effects. Although important aspects of this mechanism have been revealed, several features remain to be elucidated. We have analyzed 157 SLe sequences from 148 species from seven phyla and found a high degree of conservation among the sequences of species from the same phylum, although no considerable similarity seems to exist between sequences of species from different phyla. When analyzing case studies, we found evidence that a given SLe will always be related to a given set of transcripts in different species from the same phylum, and therefore, different SLe sequences from the same species would regulate different sets of transcripts. In addition, we have observed distinct transcript categories to be preferential targets for the SLe addition in different phyla. This work sheds light into crucial and controversial aspects of the SLTS mechanism. It represents a comprehensive study concerning various species and different characteristics of this important post-transcriptional regulatory mechanism.

Oxidative stress and DNA lesions: the role of 8-oxoguanine lesions in Trypanosoma cruzi cell viability.

The main consequence of oxidative stress is the formation of DNA lesions, which can result in genomic instability and lead to cell death. Guanine is the base that is most susceptible to oxidation, due to its low redox potential, and 8-oxoguanine (8-oxoG) is the most common lesion. These characteristics make 8-oxoG a good cellular biomarker to indicate the extent of oxidative stress. If not repaired, 8-oxoG can pair with adenine and cause a G:C to T:A transversion. When 8-oxoG is inserted during DNA replication, it could generate double-strand breaks, which makes this lesion particularly deleterious. Trypanosoma cruzi needs to address various oxidative stress situations, such as the mammalian intracellular environment and the triatomine insect gut where it replicates. We focused on the MutT enzyme, which is responsible for removing 8-oxoG from the nucleotide pool. To investigate the importance of 8-oxoG during parasite infection of mammalian cells, we characterized the MutT gene in T. cruzi (TcMTH) and generated T. cruzi parasites heterologously expressing Escherichia coli MutT or overexpressing the TcMTH enzyme. In the epimastigote form, the recombinant and wild-type parasites displayed similar growth in normal conditions, but the MutT-expressing cells were more resistant to hydrogen peroxide treatment. The recombinant parasite also displayed significantly increased growth after 48 hours of infection in fibroblasts and macrophages when compared to wild-type cells, as well as increased parasitemia in Swiss mice. In addition, we demonstrated, using western blotting experiments, that MutT heterologous expression can influence the parasite antioxidant enzyme protein levels. These results indicate the importance of the 8-oxoG repair system for cell viability.

Identification of a new Schistosoma mansoni SMYB1 partner: putative roles in RNA metabolism.

SMYB1 is a Schistosoma mansoni protein highly similar to members of the Y-box binding protein family. Similar to other homologues, SMYB1 is able to bind double- and single-stranded DNA, as well as RNA molecules. The characterization of proteins involved in the regulation of gene expression in S. mansoni is of great importance for the understanding of molecular events that control morphological and physiological changes in this parasite. Here we demonstrate that SMYB1 is located in the cytoplasm of cells from different life-cycle stages of S. mansoni, suggesting that this protein is probably acting in mRNA metabolism in the cytoplasm and corroborating previous findings from our group that showed its ability to bind RNA. Protein-protein interactions are important events in all biological processes, since most proteins execute their functions through large supramolecular structures. Yeast two-hybrid screenings using SMYB1 as bait identified a partner in S. mansoni similar to the SmD3 protein of Drosophila melanogaster (SmRNP), which is important in the assembly of small nuclear ribonucleoprotein complexes. Also, pull-down assays were conducted using immobilized GST-SMYB1 proteins and confirmed the SMYB1-SmRNP interaction. The interaction of SMYB1 with a protein involved in mRNA processing suggests that it may act in processes such as turnover, transport and stabilization of RNA molecules.

Gene identification and comparative molecular modeling of a Trypanosoma rangeli major surface protease.

Trypanosoma rangeli is a hemoflagellate parasite which is able to infect humans. Distinct from Trypanosoma cruzi, the causative agent of Chagas disease, T. rangeli is non-pathogenic to the vertebrate host. The manner by which the T. rangeli interacts with the host is still unknown, but it certainly depends on the surface molecules. Major surface proteins (MSP) are GPI-anchored, zinc-dependent metalloproteases present in the surface of all trypanosomatids studied so far, which are implicated as virulence factors in pathogenic trypanosomatids, such as Leishmania spp and T. cruzi. The aims of this work were to generate the complete sequence of a T. rangeli MSP (TrMSP) gene and to determine the 3D-structure of the predicted protein by homology modeling. The plasmid bearing a complete copy of a TrMSP gene was completely sequenced and the predicted protein was modeled using Modeller software. Results indicate that TrMSP open reading frame (ORF) codes for a predicted 588 amino acid protein and shows all elements required for its posttranslational processing. Multiple sequence alignment of TrMSP with other trypanosomatids' MSPs showed an extensive conservation of the N-terminal and central regions and a more divergent C-terminal region. Leishmania major MSP (LmMSP), which had its crystal structure previously determined, has an overall 35% identity with TrMSP. This identity allowed the comparative molecular modeling of TrMSP, which demonstrated a high degree of structural conservation between MSPs from other trypanosomatids (TrypMSPs). All modeled MSPs have a conserved folding pattern, apart from structural divergences in the C-domain and discrete differences of charge and topology in the catalytic cleft, and present the same geometry of the canonical HEXXH zinc-binding motif. The determination of surface charges of the molecules revealed that TrMSP is a predominantly positive protein, whereas LmMSP and Trypanosoma cruzi MSP (TcMSP) are negative proteins, suggesting that substrates recognized by TcMSP and LmMSP could not interact with TrMSP. Moreover, the comparison between TrMSP and TcMSP protein sequences has revealed 45 non-neutral amino acid substitutions, which can be further assessed through protein engineering. The characteristics of TrMSP could explain, at least in part, the lack of pathogenicity of T. rangeli to humans and point to the necessity of identifying the biological targets of this enzyme.

KOMODO: a web tool for detecting and visualizing biased distribution of groups of homologous genes in monophyletic taxa.

The enrichment analysis is a standard procedure to interpret 'omics' experiments that generate large gene lists as outputs, such as transcriptomics and protemics. However, despite the huge success of enrichment analysis in these classes of experiments, there is a surprising lack of application of this methodology to survey other categories of large-scale biological data available. Here, we report Kegg Orthology enrichMent-Online DetectiOn (KOMODO), a web tool to systematically investigate groups of monophyletic genomes in order to detect significantly enriched groups of homologous genes in one taxon when compared with another. The results are displayed in their proper biochemical roles in a visual, explorative way, allowing users to easily formulate and investigate biological hypotheses regarding the taxonomical distribution of genomic elements. We validated KOMODO by analyzing portions of central carbon metabolism in two taxa extensively studied regarding their carbon metabolism profile (Enterobacteriaceae family and Lactobacillales order). Most enzymatic activities significantly biased were related to known key metabolic traits in these taxa, such as the distinct fates of pyruvate (the known tendency of lactate production in Lactobacillales and its complete oxidation in Enterobacteriaceae), demonstrating that KOMODO could detect biologically meaningful differences in the frequencies of shared genomic elements among taxa. KOMODO is freely available at http://komodotool.org.

Identification of the bacterial community responsible for traditional fermentation during sour cassava starch, cachaça and minas cheese production using culture-independent 16s rRNA gene sequence analysis

We used a cultivation-independent, clone library-based 16S rRNA gene sequence analysis to identify bacterial communities present during traditional fermentation in sour cassava starch, cachaça and cheese production in Brazil. Partial 16S rRNA gene clone sequences from sour cassava starch samples collected on day five of the fermentation process indicated that Leuconostoc citreum was the most prevalent species, representing 47.6 percent of the clones. After 27 days of fermentation, clones (GenBank accession numbers GQ999786 and GQ999788) related to unculturable bacteria were the most prevalent, representing 43.8 percent of the clones from the bacterial community analyzed. The clone represented by the sequence GQ999786 was the most prevalent at the end of the fermentation period. The majority of clones obtained from cachaça samples during the fermentation of sugar cane juice were from the genus Lactobacillus. Lactobacillus nagelli was the most prevalent at the beginning of the fermentation process, representing 76.9 percent of the clones analyzed. After 21 days, Lactobacillus harbinensis was the most prevalent species, representing 75 percent of the total clones. At the end of the fermentation period, Lactobacillus buchneri was the most prevalent species, representing 57.9 percent of the total clones. In the Minas cheese samples, Lactococcus lactis was the most prevalent species after seven days of ripening. After 60 days of ripening, Streptococcus salivarius was the most prevalent species. Our data show that these three fermentation processes are conducted by a succession of bacterial species, of which lactic acid bacteria are the most prevalent.

Evidence for reductive genome evolution and lateral acquisition of virulence functions in two Corynebacterium pseudotuberculosis strains.

BACKGROUND: Corynebacterium pseudotuberculosis, a gram-positive, facultative intracellular pathogen, is the etiologic agent of the disease known as caseous lymphadenitis (CL). CL mainly affects small ruminants, such as goats and sheep; it also causes infections in humans, though rarely. This species is distributed worldwide, but it has the most serious economic impact in Oceania, Africa and South America. Although C. pseudotuberculosis causes major health and productivity problems for livestock, little is known about the molecular basis of its pathogenicity. METHODOLOGY AND FINDINGS: We characterized two C. pseudotuberculosis genomes (Cp1002, isolated from goats; and CpC231, isolated from sheep). Analysis of the predicted genomes showed high similarity in genomic architecture, gene content and genetic order. When C. pseudotuberculosis was compared with other Corynebacterium species, it became evident that this pathogenic species has lost numerous genes, resulting in one of the smallest genomes in the genus. Other differences that could be part of the adaptation to pathogenicity include a lower GC content, of about 52%, and a reduced gene repertoire. The C. pseudotuberculosis genome also includes seven putative pathogenicity islands, which contain several classical virulence factors, including genes for fimbrial subunits, adhesion factors, iron uptake and secreted toxins. Additionally, all of the virulence factors in the islands have characteristics that indicate horizontal transfer. CONCLUSIONS: These particular genome characteristics of C. pseudotuberculosis, as well as its acquired virulence factors in pathogenicity islands, provide evidence of its lifestyle and of the pathogenicity pathways used by this pathogen in the infection process. All genomes cited in this study are available in the NCBI Genbank database (http://www.ncbi.nlm.nih.gov/genbank/) under accession numbers CP001809 and CP001829.