Putative mobilized colistin resistance genes in the human gut microbiome.

BACKGROUND: The high incidence of bacterial genes that confer resistance to last-resort antibiotics, such as colistin, caused by mobilized colistin resistance (mcr) genes, poses an unprecedented threat to human health. Understanding the spread, evolution, and distribution of such genes among human populations will help in the development of strategies to diminish their occurrence. To tackle this problem, we investigated the distribution and prevalence of potential mcr genes in the human gut microbiome using a set of bioinformatics tools to screen the Unified Human Gastrointestinal Genome (UHGG) collection for the presence, synteny and phylogeny of putative mcr genes, and co-located antibiotic resistance genes. RESULTS: A total of 2079 antibiotic resistance genes (ARGs) were classified as mcr genes in 2046 metagenome assembled genomes (MAGs), distributed across 1596 individuals from 41 countries, of which 215 were identified in plasmidial contigs. The genera that presented the largest number of mcr-like genes were Suterella and Parasuterella. Other potential pathogens carrying mcr genes belonged to the genus Vibrio, Escherichia and Campylobacter. Finally, we identified a total of 22,746 ARGs belonging to 21 different classes in the same 2046 MAGs, suggesting multi-resistance potential in the corresponding bacterial strains, increasing the concern of ARGs impact in the clinical settings. CONCLUSION: This study uncovers the diversity of mcr-like genes in the human gut microbiome. We demonstrated the cosmopolitan distribution of these genes in individuals worldwide and the co-presence of other antibiotic resistance genes, including Extended-spectrum Beta-Lactamases (ESBL). Also, we described mcr-like genes fused to a PAP2-like domain in S. wadsworthensis. These novel sequences increase our knowledge about the diversity and evolution of mcr-like genes. Future research should focus on activity, genetic mobility and a potential colistin resistance in the corresponding strains to experimentally validate those findings.

Core of the saliva microbiome: an analysis of the MG-RAST data.

BACKGROUND: Oral microbiota is considered as the second most complex in the human body and its dysbiosis can be responsible for oral diseases. Interactions between the microorganism communities and the host allow establishing the microbiological proles. Identifying the core microbiome is essential to predicting diseases and changes in environmental behavior from microorganisms. METHODS: Projects containing the term "SALIVA", deposited between 2014 and 2019 were recovered on the MG-RAST portal. Quality (Failed), taxonomic prediction (Unknown and Predicted), species richness (Rarefaction), and species diversity (Alpha) were analyzed according to sequencing approaches (Amplicon sequencing and Shotgun metagenomics). All data were checked for normality and homoscedasticity. Metagenomic projects were compared using the Mann-Whitney U test and Spearman's correlation. Microbiome cores were inferred by Principal Component Analysis. For all statistical tests, p < 0.05 was used. RESULTS: The study was performed with 3 projects, involving 245 Amplicon and 164 Shotgun metagenome datasets. All comparisons of variables, according to the type of sequencing, showed significant differences, except for the Predicted. In Shotgun metagenomics datasets the highest correlation was between Rarefaction and Failed (r = - 0.78) and the lowest between Alpha and Unknown (r = - 0.12). In Amplicon sequencing datasets, the variables Rarefaction and Unknown (r = 0.63) had the highest correlation and the lowest was between Alpha and Predicted (r = - 0.03). Shotgun metagenomics datasets showed a greater number of genera than Amplicon. Propionibacterium, Lactobacillus, and Prevotella were the most representative genera in Amplicon sequencing. In Shotgun metagenomics, the most representative genera were Escherichia, Chitinophaga, and Acinetobacter. CONCLUSIONS: Core of the salivary microbiome and genera diversity are dependent on the sequencing approaches. Available data suggest that Shotgun metagenomics and Amplicon sequencing have similar sensitivities to detect the taxonomic level investigated, although Shotgun metagenomics allows a deeper analysis of the microorganism diversity. Microbiome studies must consider characteristics and limitations of the sequencing approaches. Were identified 20 genera in the core of saliva microbiome, regardless of the health condition of the host. Some bacteria of the core need further study to better understand their role in the oral cavity.

2020 taxonomic update for phylum Negarnaviricota (Riboviria: Orthornavirae), including the large orders Bunyavirales and Mononegavirales.

In March 2020, following the annual International Committee on Taxonomy of Viruses (ICTV) ratification vote on newly proposed taxa, the phylum Negarnaviricota was amended and emended. At the genus rank, 20 new genera were added, two were deleted, one was moved, and three were renamed. At the species rank, 160 species were added, four were deleted, ten were moved and renamed, and 30 species were renamed. This article presents the updated taxonomy of Negarnaviricota as now accepted by the ICTV.

SparkBLAST: scalable BLAST processing using in-memory operations.

BACKGROUND: The demand for processing ever increasing amounts of genomic data has raised new challenges for the implementation of highly scalable and efficient computational systems. In this paper we propose SparkBLAST, a parallelization of a sequence alignment application (BLAST) that employs cloud computing for the provisioning of computational resources and Apache Spark as the coordination framework. As a proof of concept, some radionuclide-resistant bacterial genomes were selected for similarity analysis. RESULTS: Experiments in Google and Microsoft Azure clouds demonstrated that SparkBLAST outperforms an equivalent system implemented on Hadoop in terms of speedup and execution times. CONCLUSIONS: The superior performance of SparkBLAST is mainly due to the in-memory operations available through the Spark framework, consequently reducing the number of local I/O operations required for distributed BLAST processing.

Bacterial and Archaeal Communities Variability Associated with Upwelling and Anthropogenic Pressures in the Protection Area of Arraial do Cabo (Cabo Frio region - RJ).

Upwelling systems contain a high diversity of pelagic microorganisms and their composition and activity are defined by factors like temperature and nutrient concentration. Denaturing gradient gel electrophoresis (DGGE) technique was used to verify the spatial and temporal genetic variability of Bacteria and Archaea in two stations of the Arraial do Cabo coastal region, one under upwelling pressure and another under anthropogenic pressure. In addition, biotic and abiotic variables were measured in surface and deep waters from three other stations between these stations. Six samplings were done during a year and adequately represented the degrees of upwelling and anthropogenic pressures to the system. Principal Component Analysis (PCA) showed negative correlations between the concentrations of ammonia and phosphorous with prokaryotic secondary production and the total heterotrophic bacteria. PCA also showed negative correlation between temperature and the abundance of prokaryotic cells. Bacterial and archaeal compositions were changeable as were the oceanographic conditions, and upwelling had a regional pressure while anthropogenic pressure was punctual. We suggest that the measurement of prokaryotic secondary production was associated with both Bacteria and Archaea activities, and that substrate availability and temperature determine nutrients cycling.

Metagenomic Analysis of Upwelling-Affected Brazilian Coastal Seawater Reveals Sequence Domains of Type I PKS and Modular NRPS.

Marine environments harbor a wide range of microorganisms from the three domains of life. These microorganisms have great potential to enable discovery of new enzymes and bioactive compounds for industrial use. However, only ~1% of microorganisms from the environment can currently be identified through cultured isolates, limiting the discovery of new compounds. To overcome this limitation, a metagenomics approach has been widely adopted for biodiversity studies on samples from marine environments. In this study, we screened metagenomes in order to estimate the potential for new natural compound synthesis mediated by diversity in the Polyketide Synthase (PKS) and Nonribosomal Peptide Synthetase (NRPS) genes. The samples were collected from the Praia dos Anjos (Angel's Beach) surface water-Arraial do Cabo (Rio de Janeiro state, Brazil), an environment affected by upwelling. In order to evaluate the potential for screening natural products in Arraial do Cabo samples, we used KS (keto-synthase) and C (condensation) domains (from PKS and NRPS, respectively) to build Hidden Markov Models (HMM) models. From both samples, a total of 84 KS and 46 C novel domain sequences were obtained, showing the potential of this environment for the discovery of new genes of biotechnological interest. These domains were classified by phylogenetic analysis and this was the first study conducted to screen PKS and NRPS genes in an upwelling affected sample.

A comparative genomics study of the microbiome and freshwater resistome in Southern Pantanal.

This study explores the resistome and bacterial diversity of two small lakes in the Southern Pantanal, one in Aquidauana sub-region, close to a farm, and one in Abobral sub-region, an environmentally preserved area. Shotgun metagenomic sequencing data from water column samples collected near and far from the floating macrophyte Eichhornia crassipes were used. The Abobral small lake exhibited the highest diversity and abundance of antibiotic resistance genes (ARGs), antibiotic resistance classes (ARGCs), phylum, and genus. RPOB2 and its resistance class, multidrug resistance, were the most abundant ARG and ARGC, respectively. Pseudomonadota was the dominant phylum across all sites, and Streptomyces was the most abundant genus considering all sites.

Genomic analysis of an outbreak of bovine tuberculosis in a man-made multi-host species system: A call for action on wildlife in Brazil.

We report on a 15-year-long outbreak of bovine tuberculosis (bTB) in wildlife from a Brazilian safari park. A timeline of diagnostic events and whole-genome sequencing (WGS) of 21 Mycobacterium bovis isolates from deer and llamas were analyzed. Accordingly, from 2003 to 2018, at least 16 animals, from eight species, died due to TB, which is likely an underestimated number. In three occasions since 2013, the deer presented positive tuberculin tests, leading to the park closure and culling of all deer. WGS indicated that multiple M. bovis strains were circulating, with at least three founding introductions since the park inauguration in 1977. Using a previously sequenced dataset of 71 M. bovis genomes from cattle, we found no recent transmission events between nearby farms and the park based on WGS. Lastly, by discussing socio-economic and environmental factors escaping current regulatory gaps that were determinant of this outbreak, we pledge for the development of a plan to report and control bTB in wildlife in Brazil.

Genetic Diversity and Potential Paths of Transmission of Mycobacterium bovis in the Amazon: The Discovery of M. bovis Lineage Lb1 Circulating in South America.

Bovine tuberculosis (bTB) has yet to be eradicated in Brazil. Herds of cattle and buffalo are important sources of revenue to people living in the banks of the Amazon River basin. A better understanding of Mycobacterium bovis (M. bovis) populational structure and transmission dynamics affecting these animals can significantly contribute in efforts to improve their sanitary status. Herein, we sequenced the whole genome of 22 M. bovis isolates (15 from buffalo and 7 from cattle) from 10 municipalities in the region of the Lower Amazon River Basin in Brazil and performed phylogenomic analysis and Single Nucleotide Polymorphism (SNP)-based transmission inference to evaluate population structure and transmission networks. Additionally, we compared these genomes to others obtained in unrelated studies in the Marajó Island (n = 15) and worldwide (n = 128) to understand strain diversity in the Amazon and to infer M. bovis lineages. Our results show a higher genomic diversity of M. bovis genomes obtained in the Lower Amazon River region when compared to the Marajó Island, while no significant difference was observed between M. bovis genomes obtained from cattle and buffalo (p ≥ 0.05). This high genetic diversity is reflected by the weak phylogenetic clustering of M. bovis from the Lower Amazon River region based on geographic proximity and in the detection of only two putative transmission clusters in the region. One of these clusters is the first description of inter-species transmission between cattle and buffalo in the Amazon, bringing implications to the bTB control program. Surprisingly, two M. bovis lineages were detected in our dataset, namely Lb1 and Lb3, constituting the first description of Lb1 in South America. Most of the strains of this study (13/22) and all 15 strains of the Marajó Island carried no clonal complex marker, suggesting that the recent lineage classification better describe the diversity of M. bovis in the Amazon.

ProtozoaDB: dynamic visualization and exploration of protozoan genomes.

ProtozoaDB (http://www.biowebdb.org/protozoadb) is being developed to initially host both genomics and post-genomics data from Plasmodium falciparum, Entamoeba histolytica, Trypanosoma brucei, T. cruzi and Leishmania major, but will hopefully host other protozoan species as more genomes are sequenced. It is based on the Genomics Unified Schema and offers a modern Web-based interface for user-friendly data visualization and exploration. This database is not intended to duplicate other similar efforts such as GeneDB, PlasmoDB, TcruziDB or even TDRtargets, but to be complementary by providing further analyses with emphasis on distant similarities (HMM-based) and phylogeny-based annotations including orthology analysis. ProtozoaDB will be progressively linked to the above-mentioned databases, focusing in performing a multi-source dynamic combination of information through advanced interoperable Web tools such as Web services. Also, to provide Web services will allow third-party software to retrieve and use data from ProtozoaDB in automated pipelines (workflows) or other interoperable Web technologies, promoting better information reuse and integration. We also expect ProtozoaDB to catalyze the development of local and regional bioinformatics capabilities (research and training), and therefore promote/enhance scientific advancement in developing countries.

Global ocean resistome revealed: Exploring antibiotic resistance gene abundance and distribution in TARA Oceans samples.

BACKGROUND: The rise of antibiotic resistance (AR) in clinical settings is of great concern. Therefore, the understanding of AR mechanisms, evolution, and global distribution is a priority for patient survival. Despite all efforts in the elucidation of AR mechanisms in clinical strains, little is known about its prevalence and evolution in environmental microorganisms. We used 293 metagenomic samples from the TARA Oceans project to detect and quantify environmental antibiotic resistance genes (ARGs) using machine learning tools. RESULTS: After manual curation of ARGs, their abundance and distribution in the global ocean are presented. Additionally, the potential of horizontal ARG transfer by plasmids and their correlation with environmental and geographical parameters is shown. A total of 99,205 environmental open reading frames (ORFs) were classified as 1 of 560 different ARGs conferring resistance to 26 antibiotic classes. We found 24,567 ORFs in putative plasmid sequences, suggesting the importance of mobile genetic elements in the dynamics of environmental ARG transmission. Moreover, 4,804 contigs with >=2 putative ARGs were found, including 2 plasmid-like contigs with 5 different ARGs, highlighting the potential presence of multi-resistant microorganisms in the natural ocean environment. Finally, we identified ARGs conferring resistance to some of the most relevant clinical antibiotics, revealing the presence of 15 ARGs similar to mobilized colistin resistance genes (mcr) with high abundance on polar biomes. Of these, 5 are assigned to Psychrobacter, a genus including opportunistic human pathogens. CONCLUSIONS: This study uncovers the diversity and abundance of ARGs in the global ocean metagenome. Our results are available on Zenodo in MySQL database dump format, and all the code used for the analyses, including a Jupyter notebook js avaliable on Github. We also developed a dashboard web application (http://www.resistomedb.com) for data visualization.

EST sequencing of blood-fed and Leishmania-infected midgut of Lutzomyia longipalpis, the principal visceral leishmaniasis vector in the Americas.

Leishmaniasis is an important worldwide public health problem. Visceral leishmaniasis caused by Leishmania infantum chagasi is mainly transmitted by Lutzomyia longipalpis in the Americas. Leishmania development within the sand fly vector is mostly restricted to the midgut. Thus, a comparative analysis of blood-fed versus infected midguts may provide an invaluable insight into various aspects of sand fly immunity, physiology of blood digestion, and, more importantly, of Leishmania development. To that end, we have engaged in a study to identify expressed sequenced tags (ESTs) from L. longipalpis cDNA libraries produced from midguts dissected at different times post blood meal and also after artificial infection with L. i. chagasi. A total of 2,520 ESTs were obtained and, according to the quality of the sequencing data obtained, assembled into 378 clusters and 1,526 individual sequences or singletons totalizing 1,904 sequences. Several sequences associated with defense, apoptosis, RNAi, and digestion processes were annotated. The data presented here increases current knowledge on the New World sand fly transcriptome, contributing to the understanding of various aspects of the molecular physiology of L. longipalpis, and mechanisms underlying the relationship of this sand fly species with L. i. chagasi.

Recovering Genomics Clusters of Secondary Metabolites from Lakes Using Genome-Resolved Metagenomics.

Metagenomic approaches became increasingly popular in the past decades due to decreasing costs of DNA sequencing and bioinformatics development. So far, however, the recovery of long genes coding for secondary metabolites still represents a big challenge. Often, the quality of metagenome assemblies is poor, especially in environments with a high microbial diversity where sequence coverage is low and complexity of natural communities high. Recently, new and improved algorithms for binning environmental reads and contigs have been developed to overcome such limitations. Some of these algorithms use a similarity detection approach to classify the obtained reads into taxonomical units and to assemble draft genomes. This approach, however, is quite limited since it can classify exclusively sequences similar to those available (and well classified) in the databases. In this work, we used draft genomes from Lake Stechlin, north-eastern Germany, recovered by MetaBat, an efficient binning tool that integrates empirical probabilistic distances of genome abundance, and tetranucleotide frequency for accurate metagenome binning. These genomes were screened for secondary metabolism genes, such as polyketide synthases (PKS) and non-ribosomal peptide synthases (NRPS), using the Anti-SMASH and NAPDOS workflows. With this approach we were able to identify 243 secondary metabolite clusters from 121 genomes recovered from our lake samples. A total of 18 NRPS, 19 PKS, and 3 hybrid PKS/NRPS clusters were found. In addition, it was possible to predict the partial structure of several secondary metabolite clusters allowing for taxonomical classifications and phylogenetic inferences. Our approach revealed a high potential to recover and study secondary metabolites genes from any aquatic ecosystem.

Improving model construction of profile HMMs for remote homology detection through structural alignment.

BACKGROUND: Remote homology detection is a challenging problem in Bioinformatics. Arguably, profile Hidden Markov Models (pHMMs) are one of the most successful approaches in addressing this important problem. pHMM packages present a relatively small computational cost, and perform particularly well at recognizing remote homologies. This raises the question of whether structural alignments could impact the performance of pHMMs trained from proteins in the Twilight Zone, as structural alignments are often more accurate than sequence alignments at identifying motifs and functional residues. Next, we assess the impact of using structural alignments in pHMM performance. RESULTS: We used the SCOP database to perform our experiments. Structural alignments were obtained using the 3DCOFFEE and MAMMOTH-mult tools; sequence alignments were obtained using CLUSTALW, TCOFFEE, MAFFT and PROBCONS. We performed leave-one-family-out cross-validation over super-families. Performance was evaluated through ROC curves and paired two tailed t-test. CONCLUSION: We observed that pHMMs derived from structural alignments performed significantly better than pHMMs derived from sequence alignment in low-identity regions, mainly below 20%. We believe this is because structural alignment tools are better at focusing on the important patterns that are more often conserved through evolution, resulting in higher quality pHMMs. On the other hand, sensitivity of these tools is still quite low for these low-identity regions. Our results suggest a number of possible directions for improvements in this area.

Molecular diversity of disintegrin-like domains within metalloproteinase precursors of Bothrops jararaca.

Disintegrins are small peptides isolated from the venom of several snake families which act as integrin-antagonists or agonists, interacting with a variety of biological processes mediated by integrins. In this work we describe five new disintegrin-like domains within metalloproteinase precursor sequences, obtained from a Bothrops jararaca venom gland cDNA library. Among the new disintegrin-like domains, four were contained in PIII metalloproteinase precursors, with three of them presenting ECD-motifs and one presenting a new KCD-motif. Moreover, we found three disintegrin-like domains within PII metalloproteinase precursors. Two of them are similar to the already described disintegrins jarastatin and jararacin. The third molecule is unusual, presenting some typical PIII metalloproteinase characteristics but lacking the cysteine-rich domain being, thus, classified as a PII metalloproteinase. Only few reports presented molecules with these characteristics. Sequence analysis suggests that these molecules are intermediate steps between the more ancient PIII and the more recent PII metalloproteinases. We also investigated disintegrin N-terminus diversity in B. jararaca crude venom by purifying jarastatin and jararacin and analyzing them by mass spectrometry.

Bothrops jararaca venom gland transcriptome: analysis of the gene expression pattern.

Bothrops jararaca is a pit viper responsible for the majority of snake envenoming accidents in Brazil. As an attempt to describe the transcriptional activity of the venom gland, ESTs of a cDNA library constructed from B. jararaca venom gland were generated and submitted to bioinformatics analysis. The results showed a clear predominance of transcripts coding for toxins instead of transcripts coding for proteins involved in cellular functions. Among toxins, the most frequent transcripts were from metalloproteinases (52.6%), followed by serine-proteinases (28.5%), C-type lectins (8.3%) and bradykinin-potentiating peptides (BPPs) (6.2%). Results were similar to that obtained from the transcriptome analysis of B. insularis, a phylogenetically close sister of B. jararaca, though some differences were observed and are pointed out, such as a higher amount of the hypotensive BPPs in B. insularis transcriptome (19.7%). Another striking difference observed is that PIII and PII-classes of metalloproteinases are similarly represented in B. jararaca in contrast to B. insularis, in which a predominance of PIII-class metalloproteinase, which present a more intense hemorrhagic action, is observed. These features may, in part, explain the higher potency of B. insularis venom. The results obtained can help in proteome studies, and the clones can be used to directly probe the genetic material from other snake species or to investigate differences in gene expression pattern in response to factors such as diet, aging and geographic localization.

Picoplankton Bloom in Global South? A High Fraction of Aerobic Anoxygenic Phototrophic Bacteria in Metagenomes from a Coastal Bay (Arraial do Cabo--Brazil).

Marine habitats harbor a great diversity of microorganism from the three domains of life, only a small fraction of which can be cultivated. Metagenomic approaches are increasingly popular for addressing microbial diversity without culture, serving as sensitive and relatively unbiased methods for identifying and cataloging the diversity of nucleic acid sequences derived from organisms in environmental samples. Aerobic anoxygenic phototrophic bacteria (AAP) play important roles in carbon and energy cycling in aquatic systems. In oceans, those bacteria are widely distributed; however, their abundance and importance are still poorly understood. The aim of this study was to estimate abundance and diversity of AAPs in metagenomes from an upwelling affected coastal bay in Arraial do Cabo, Brazil, using in silico screening for the anoxygenic photosynthesis core genes. Metagenomes from the Global Ocean Sample Expedition (GOS) were screened for comparative purposes. AAPs were highly abundant in the free-living bacterial fraction from Arraial do Cabo: 23.88% of total bacterial cells, compared with 15% in the GOS dataset. Of the ten most AAP abundant samples from GOS, eight were collected close to the Equator where solar irradiation is high year-round. We were able to assign most retrieved sequences to phylo-groups, with a particularly high abundance of Roseobacter in Arraial do Cabo samples. The high abundance of AAP in this tropical bay may be related to the upwelling phenomenon and subsequent picoplankton bloom. These results suggest a link between upwelling and light abundance and demonstrate AAP even in oligotrophic tropical and subtropical environments. Longitudinal studies in the Arraial do Cabo region are warranted to understand the dynamics of AAP at different locations and seasons, and the ecological role of these unique bacteria for biogeochemical and energy cycling in the ocean.

Distribution and Classification of Serine ß-Lactamases in Brazilian Hospital Sewage and Other Environmental Metagenomes Deposited in Public Databases.

ß-lactam is the most used antibiotic class in the clinical area and it acts on blocking the bacteria cell wall synthesis, causing cell death. However, some bacteria have evolved resistance to these antibiotics mainly due the production of enzymes known as ß-lactamases. Hospital sewage is an important source of dispersion of multidrug-resistant bacteria in rivers and oceans. In this work, we used next-generation DNA sequencing to explore the diversity and dissemination of serine ß-lactamases in two hospital sewage from Rio de Janeiro, Brazil (South Zone, SZ and North Zone, NZ), presenting different profiles, and to compare them with public environmental data available. Also, we propose a Hidden-Markov-Model approach to screen potential serine ß-lactamases genes (in public environments samples and generated hospital sewage data), exploring its evolutionary relationships. Due to the high variability in ß-lactamases, we used a position-specific scoring matrix search method (RPS-BLAST) against conserved domain database profiles (CDD, Pfam, and COG) followed by visual inspection to detect conserved motifs, to increase the reliability of the results and remove possible false positives. We were able to identify novel ß-lactamases from Brazilian hospital sewage and to estimate relative abundance of its types. The highest relative abundance found in SZ was the Class A (50%), while Class D is predominant in NZ (55%). CfxA (65%) and ACC (47%) types were the most abundant genes detected in SZ, while in NZ the most frequent were OXA-10 (32%), CfxA (28%), ACC (21%), CEPA (20%), and FOX (19%). Phylogenetic analysis revealed ß-lactamases from Brazilian hospital sewage grouped in the same clade and close to sequences belonging to Firmicutes and Bacteroidetes groups, but distant from potential ß-lactamases screened from public environmental data, that grouped closer to ß-lactamases of Proteobacteria. Our results demonstrated that HMM-based approach identified homologs of serine ß-lactamases, indicating the specificity and high sensitivity of this approach in large datasets, contributing for the identification and classification of a large number of homologous genes, comprising possible new ones. Phylogenetic analysis revealed the potential reservoir of ß-lactam resistance genes in the environment, contributing to understanding the evolution and dissemination of these genes.

Exploring the genome of Trypanosoma vivax through GSS and in silico comparative analysis.

A survey of the Trypanosoma vivax genome was carried out by the genome sequence survey (GSS) approach resulting in 1,086 genomic sequences. A total of 455 high-quality GSS sequences were generated, consisting of 331 non-redundant sequences distributed in 264 singlets and 67 clusters in a total of 135.5 Kb of the T. vivax genome. The estimation of the overall G+C content, and the prediction of the presence of ORFs and putative genes were carried out using the Glimmer and Jemboss packages. Analysis of the obtained sequences was carried out by BLAST programs against 12 different databases and also using the Conserved Domain Database, InterProScan, and tRNAscan-SE. Along with the existing 23 T. vivax entries in the GenBank, the 32 putative genes predicted and the 331 non-redundant GSS sequences reported herein represent new potential markers for the development of PCRbased assays for specific diagnosis and typing of Trypanosoma vivax.