Outbreak report of polymyxin-carbapenem-resistant Klebsiella pneumoniae causing untreatable infections evidenced by synergy tests and bacterial genomes.

Polymyxin-carbapenem-resistant Klebsiella pneumoniae (PCR-Kp) with pan (PDR)- or extensively drug-resistant phenotypes has been increasingly described worldwide. Here, we report a PCR-Kp outbreak causing untreatable infections descriptively correlated with bacterial genomes. Hospital-wide surveillance of PCR-Kp was initiated in December-2014, after the first detection of a K. pneumoniae phenotype initially classified as PDR, recovered from close spatiotemporal cases of a sentinel hospital in Rio de Janeiro. Whole-genome sequencing of clinical PCR-Kp was performed to investigate similarities and dissimilarities in phylogeny, resistance and virulence genes, plasmid structures and genetic polymorphisms. A target phenotypic profile was detected in 10% (12/117) of the tested K. pneumoniae complex bacteria recovered from patients (8.5%, 8/94) who had epidemiological links and were involved in intractable infections and death, with combined therapeutic drugs failing to meet synergy. Two resistant bacterial clades belong to the same transmission cluster (ST437) or might have different sources (ST11). The severity of infection was likely related to patients' comorbidities, lack of antimicrobial therapy and predicted bacterial genes related to high resistance, survival, and proliferation. This report contributes to the actual knowledge about the natural history of PCR-Kp infection, while reporting from a time when there were no licensed drugs in the world to treat some of these infections. More studies comparing clinical findings with bacterial genetic markers during clonal spread are needed.

Uncovering Pseudogenes and Intergenic Protein-coding Sequences in TriTryps' Genomes.

Trypanosomatids belong to a remarkable group of unicellular, parasitic organisms of the order Kinetoplastida, an early diverging branch of the phylogenetic tree of eukaryotes, exhibiting intriguing biological characteristics affecting gene expression (intronless polycistronic transcription, trans-splicing, and RNA editing), metabolism, surface molecules, and organelles (compartmentalization of glycolysis, variation of the surface molecules, and unique mitochondrial DNA), cell biology and life cycle (phagocytic vacuoles evasion and intricate patterns of cell morphogenesis). With numerous genomic-scale data of several trypanosomatids becoming available since 2005 (genomes, transcriptomes, and proteomes), the scientific community can further investigate the mechanisms underlying these unusual features and address other unexplored phenomena possibly revealing biological aspects of the early evolution of eukaryotes. One fundamental aspect comprises the processes and mechanisms involved in the acquisition and loss of genes throughout the evolutionary history of these primitive microorganisms. Here, we present a comprehensive in silico analysis of pseudogenes in three major representatives of this group: Leishmania major, Trypanosoma brucei, and Trypanosoma cruzi. Pseudogenes, DNA segments originating from altered genes that lost their original function, are genomic relics that can offer an essential record of the evolutionary history of functional genes, as well as clues about the dynamics and evolution of hosting genomes. Scanning these genomes with functional proteins as proxies to reveal intergenic regions with protein-coding features, relying on a customized threshold to distinguish statistically and biologically significant sequence similarities, and reassembling remnant sequences from their debris, we found thousands of pseudogenes and hundreds of open reading frames, with particular characteristics in each trypanosomatid: mutation profile, number, content, density, codon bias, average size, single- or multi-copy gene origin, number and type of mutations, putative primitive function, and transcriptional activity. These features suggest a common process of pseudogene formation, different patterns of pseudogene evolution and extant biological functions, and/or distinct genome organization undertaken by those parasites during evolution, as well as different evolutionary and/or selective pressures acting on distinct lineages.

Complexation of a Polypeptide-Polyelectrolytes Bioparticle as a Biomaterial of Antibacterial Activity.

The development of biomaterials to enable application of antimicrobial peptides represents a strategy of high and current interest. In this study, a bioparticle was produced by the complexation between an antimicrobial polypeptide and the biocompatible and biodegradable polysaccharides chitosan-N-arginine and alginate, giving rise to a colloidal polyelectrolytic complex of pH-responsive properties. The inclusion of the polypeptide in the bioparticle structure largely increases the binding sites of complexation during the bioparticles production, leading to its effective incorporation. After lyophilization, detailed evaluation of colloidal structure of redispersed bioparticles evidenced nano or microparticles with size, polydispersity and zeta potential dependent on pH and ionic strength, and the dependence was not withdrawn with the polypeptide inclusion. Significant increase of pore edge tension in giant vesicles evidenced effective interaction of the polypeptide-bioparticle with lipid model membrane. Antibacterial activity against Aeromonas dhakensis was effective at 0.1% and equal for the isolated polypeptide and the same complexed in bioparticle, which opens perspectives to the composite material as an applicable antibacterial system.

The synthetic antimicrobial peptide IKR18 displays anti-infectious properties in Galleria mellonella in vivo model.

Antimicrobial peptides (AMPs) are promising tools for developing new antibiotics. We described the design of IKR18, an AMP designed with the aid of computational tools. IKR18 showed antimicrobial activity against Gram-negative and Gram-positive bacteria, including methicillin-resistant Staphylococcus aureus (MRSA). CD studies revealed that IKR18 assumes an alpha-helical structure in the membrane-mimetic environment. The action mechanism IKR18 involves damage to the bacteria membrane, as demonstrated by Sytox green uptake. Furthermore, IKR18 displayed synergic and additive effects in combination with antibiotics ciprofloxacin and vancomycin. The peptide showed anti-biofilm activity in concentration and efficiency compared with commercial antibiotics, involving the direct death of bacteria, as confirmed by scanning electron microscopy. The anti-infective activity of IKR18 was demonstrated in the Galleria mellonella model infected with S. aureus, MRSA, and Acinetobacter baumannii. The novel bioinspired peptide, IKR18, proved to be effective in the control of bacterial infection, opening opportunities for the development of further assays, including preclinical models.

ProteinWorldDB: querying radical pairwise alignments among protein sets from complete genomes.

MOTIVATION: Many analyses in modern biological research are based on comparisons between biological sequences, resulting in functional, evolutionary and structural inferences. When large numbers of sequences are compared, heuristics are often used resulting in a certain lack of accuracy. In order to improve and validate results of such comparisons, we have performed radical all-against-all comparisons of 4 million protein sequences belonging to the RefSeq database, using an implementation of the Smith-Waterman algorithm. This extremely intensive computational approach was made possible with the help of World Community Grid, through the Genome Comparison Project. The resulting database, ProteinWorldDB, which contains coordinates of pairwise protein alignments and their respective scores, is now made available. Users can download, compare and analyze the results, filtered by genomes, protein functions or clusters. ProteinWorldDB is integrated with annotations derived from Swiss-Prot, Pfam, KEGG, NCBI Taxonomy database and gene ontology. The database is a unique and valuable asset, representing a major effort to create a reliable and consistent dataset of cross-comparisons of the whole protein content encoded in hundreds of completely sequenced genomes using a rigorous dynamic programming approach. AVAILABILITY: The database can be accessed through http://proteinworlddb.org

Sequence analysis of NAT2 gene in Brazilians: identification of undescribed single nucleotide polymorphisms and molecular modeling of the N-acetyltransferase 2 protein structure.

N-Acetyltransferase 2 (NAT2) metabolizes a variety of xenobiotics that includes many drugs, chemicals and carcinogens. This enzyme is genetically variable in human populations and polymorphisms in the NAT2 gene have been associated with drug toxicity and efficacy as well as cancer susceptibility. Here, we have focused on the identification of NAT2 variants in Brazilian individuals from two different regions, Rio de Janeiro and Goiás, by direct sequencing, and on the characterization of new haplotypes after cloning and re-sequencing. Upon analysis of DNA samples from 404 individuals, six new SNPs (c.29T>C, c.152G>T, c.203G>A, c.228C>T, c.458C>T and c.600A>G) and seven new NAT2 alleles were identified with different frequencies in Rio de Janeiro and Goiás. All new SNPs were found as singletons (observed only once in 808 genes) and were confirmed by three independent technical replicates. Molecular modeling and structural analysis suggested that p.Gly51Val variant may have an important effect on substrate recognition by NAT2. We also observed that amino acid change p.Cys68Tyr would affect acetylating activity due to the resulting geometric restrictions and incompatibility of the functional group in the Tyr side chain with the admitted chemical mechanism for catalysis by NATs. Moreover, other variants, such like p.Thr153Ile, p.Thr193Met, p.Pro228Leu and p.Val280Met, may lead to the presence of hydrophobic residues on NAT2 surface involved in protein aggregation and/or targeted degradation. Finally, the new alleles NAT2*6H and NAT2*5N, which showed the highest frequency in the Brazilian populations considered in this study, may code for a slow activity. Functional studies are needed to clarify the mechanisms by which new SNPs interfere with acetylation.

Antimicrobial compounds produced by Lactobacillus sakei subsp. sakei 2a, a bacteriocinogenic strain isolated from a Brazilian meat product.

Bacteriocins produced by lactic acid bacteria are gaining increased importance due to their activity against undesirable microorganisms in foods. In this study, a concentrated acid extract of a culture of Lactobacillus sakei subsp. sakei 2a, a bacteriocinogenic strain isolated from a Brazilian pork product, was purified by cation exchange and reversed-phase chromatographic methods. The amino acid sequences of the active antimicrobial compounds determined by Edman degradation were compared to known protein sequences using the BLAST-P software. Three different antimicrobial compounds were obtained, P1, P2 and P3, and mass spectrometry indicated molecular masses of 4.4, 6.8 and 9.5 kDa, respectively. P1 corresponds to classical sakacin P, P2 is identical to the 30S ribosomal protein S21 of L. sakei subsp. sakei 23 K, and P3 is identical to a histone-like DNA-binding protein HV produced by L. sakei subsp. sakei 23 K. Total genomic DNA was extracted and used as target DNA for PCR amplification of the genes sak, lis and his involved in the synthesis of P1, P2 and P3. The fragments were cloned in pET28b expression vector and the resulting plasmids transformed in E. coli KRX competent cells. The transformants were active against Listeria monocytogenes, indicating that the activity of the classical sakacin P produced by L. sakei 2a can be complemented by other antimicrobial proteins.

AnEnPi: identification and annotation of analogous enzymes.

BACKGROUND: Enzymes are responsible for the catalysis of the biochemical reactions in metabolic pathways. Analogous enzymes are able to catalyze the same reactions, but they present no significant sequence similarity at the primary level, and possibly different tertiary structures as well. They are thought to have arisen as the result of independent evolutionary events. A detailed study of analogous enzymes may reveal new catalytic mechanisms, add information about the origin and evolution of biochemical pathways and disclose potential targets for drug development. RESULTS: In this work, we have constructed and implemented a new approach, AnEnPi (the Analogous Enzyme Pipeline), using a combination of bioinformatics tools like BLAST, HMMer, and in-house scripts, to assist in the identification, annotation, comparison and study of analogous and homologous enzymes. The algorithm for the detection of analogy is based i) on the construction of groups of homologous enzymes and ii) on the identification of cases where a given enzymatic activity is performed by two or more proteins without significant similarity between their primary structures. We applied this approach to a dataset obtained from KEGG Comprising all annotated enzymes, which resulted in the identification of 986 EC classes where putative analogy was detected (40.5% of all EC classes). AnEnPi is of considerable value in the construction of initial datasets that can be further curated, particularly in gene and genome annotation, in studies involving molecular evolution and metabolism and in the identification of new potential drug targets. CONCLUSION: AnEnPi is an efficient tool for detection and annotation of analogous enzymes and other enzymes in whole genomes. It is available for academic use at: http://bioinfo.pdtis.fiocruz.br/AnEnPi/

ReRep: computational detection of repetitive sequences in genome survey sequences (GSS).

BACKGROUND: Genome survey sequences (GSS) offer a preliminary global view of a genome since, unlike ESTs, they cover coding as well as non-coding DNA and include repetitive regions of the genome. A more precise estimation of the nature, quantity and variability of repetitive sequences very early in a genome sequencing project is of considerable importance, as such data strongly influence the estimation of genome coverage, library quality and progress in scaffold construction. Also, the elimination of repetitive sequences from the initial assembly process is important to avoid errors and unnecessary complexity. Repetitive sequences are also of interest in a variety of other studies, for instance as molecular markers. RESULTS: We designed and implemented a straightforward pipeline called ReRep, which combines bioinformatics tools for identifying repetitive structures in a GSS dataset. In a case study, we first applied the pipeline to a set of 970 GSSs, sequenced in our laboratory from the human pathogen Leishmania braziliensis, the causative agent of leishmaniosis, an important public health problem in Brazil. We also verified the applicability of ReRep to new sequencing technologies using a set of 454-reads of an Escheria coli. The behaviour of several parameters in the algorithm is evaluated and suggestions are made for tuning of the analysis. CONCLUSION: The ReRep approach for identification of repetitive elements in GSS datasets proved to be straightforward and efficient. Several potential repetitive sequences were found in a L. braziliensis GSS dataset generated in our laboratory, and further validated by the analysis of a more complete genomic dataset from the EMBL and Sanger Centre databases. ReRep also identified most of the E. coli K12 repeats prior to assembly in an example dataset obtained by automated sequencing using 454 technology. The parameters controlling the algorithm behaved consistently and may be tuned to the properties of the dataset, in particular to the length of sequencing reads and the genome coverage. ReRep is freely available for academic use at http://bioinfo.pdtis.fiocruz.br/ReRep/.

Mining of potential drug targets through the identification of essential and analogous enzymes in the genomes of pathogens of Glycine max, Zea mays and Solanum lycopersicum.

Pesticides are one of the most widely used pest and disease control measures in plant crops and their indiscriminate use poses a direct risk to the health of populations and environment around the world. As a result, there is a great need for the development of new, less toxic molecules to be employed against plant pathogens. In this work, we employed an in silico approach to study the genes coding for enzymes of the genomes of three commercially important plants, soybean (Glycine max), tomato (Solanum lycopersicum) and corn (Zea mays), as well as 15 plant pathogens (4 bacteria and 11 fungi), focusing on revealing a set of essential and non-homologous isofunctional enzymes (NISEs) that could be prioritized as drug targets. By combining sequence and structural data, we obtained an initial set of 568 cases of analogy, of which 97 were validated and further refined, revealing a subset of 29 essential enzymatic activities with a total of 119 different structural forms, most belonging to central metabolic routes, including the carbohydrate metabolism, the metabolism of amino acids, among others. Further, another subset of 26 enzymatic activities possess a tertiary structure specific for the pathogen, not present in plants, men and Apis mellifera, which may be of importance for the development of specific enzymatic inhibitors against plant diseases that are less harmful to humans and the environment.

Systematic Identification and Classification of ß-Lactamases Based on Sequence Similarity Criteria: ß-Lactamase Annotation.

ß-lactamases, the enzymes responsible for resistance to ß-lactam antibiotics, are widespread among prokaryotic genera. However, current ß-lactamase classification schemes do not represent their present diversity. Here, we propose a workflow to identify and classify ß-lactamases. Initially, a set of curated sequences was used as a model for the construction of profiles Hidden Markov Models (HMM), specific for each ß-lactamase class. An extensive, nonredundant set of ß-lactamase sequences was constructed from 7 different resistance proteins databases to test the methodology. The profiles HMM were improved for their specificity and sensitivity and then applied to fully assembled genomes. Five hierarchical classification levels are described, and a new class of ß-lactamases with fused domains is proposed. Our profiles HMM provide a better annotation of ß-lactamases, with classes and subclasses defined by objective criteria such as sequence similarity. This classification offers a solid base to the elaboration of studies on the diversity, dispersion, prevalence, and evolution of the different classes and subclasses of this critical enzymatic activity.

Functional Analogy in Human Metabolism: Enzymes with Different Biological Roles or Functional Redundancy?

Since enzymes catalyze almost all chemical reactions that occur in living organisms, it is crucial that genes encoding such activities are correctly identified and functionally characterized. Several studies suggest that the fraction of enzymatic activities in which multiple events of independent origin have taken place during evolution is substantial. However, this topic is still poorly explored, and a comprehensive investigation of the occurrence, distribution, and implications of these events has not been done so far. Fundamental questions, such as how analogous enzymes originate, why so many events of independent origin have apparently occurred during evolution, and what are the reasons for the coexistence in the same organism of distinct enzymatic forms catalyzing the same reaction, remain unanswered. Also, several isofunctional enzymes are still not recognized as nonhomologous, even with substantial evidence indicating different evolutionary histories. In this work, we begin to investigate the biological significance of the cooccurrence of nonhomologous isofunctional enzymes in human metabolism, characterizing functional analogous enzymes identified in metabolic pathways annotated in the human genome. Our hypothesis is that the coexistence of multiple enzymatic forms might not be interpreted as functional redundancy. Instead, these enzymatic forms may be implicated in distinct (and probably relevant) biological roles.

Differentiation of Leishmania major is impaired by over-expression of pyroglutamyl peptidase I.

Pyroglutamyl peptidases I (PPI) are cysteine peptidases of the clan CF, family C15, which hydrolyse N-terminal l-pyroglutamyl residues (l-pGlu). The l-pGlu modification is a post-transcriptional modification that confers relative aminopeptidase resistance and, in some cases, is essential to the modified peptides' biological activity. PPIs have been identified in a variety of organisms, although definitive biological functions have yet to be attributed to them. The L. major PPI was expressed in Escherichia coli as active recombinant enzyme, and shown to have biochemical properties more similar to mammalian than bacterial PPIs. The LmPPI active site catalytic triad of E101, C210, and H234 was confirmed by mutagenesis. PPI activity was detected in L. major promastigotes, and the enzyme localised to the parasite cytosol. No detectable phenotype could be observed for L. major PPI-deficient mutants, which retained infectivity to macrophages in vitro and mice. However, over-expression of the active PPI, but not inactive PPI(C210A), in L. major impaired differentiation from the procyclic promastigote to the infective metacyclic promastigote. Susceptibility to a natural l-pGlu-modified antimicrobial peptide, gomesin, was tested using the different cell lines, which were all equally susceptible. Whilst PPI is widespread through the eukaryotic kingdom, this study now suggests that the enzyme is not essential for normal eukaryotic cell function. However, PPI could be involved in regulating the action of l-pGlu-modified peptides required for differentiation of L. major.

BioParser: a tool for processing of sequence similarity analysis reports.

UNLABELLED: The widely used programs BLAST (in this article, 'BLAST' includes both the National Center for Biotechnology Information [NCBI] BLAST and the Washington University version WU BLAST) and FASTA for similarity searches in nucleotide and protein databases usually result in copious output. However, when large query sets are used, human inspection rapidly becomes impractical. BioParser is a Perl program for parsing BLAST and FASTA reports. Making extensive use of the BioPerl toolkit, the program filters, stores and returns components of these reports in either ASCII or HTML format. BioParser is also capable of automatically feeding a local MySQL database with the parsed information, allowing subsequent filtering of hits and/or alignments with specific attributes. For this reason, BioParser is a valuable tool for large-scale similarity analyses by improving the access to the information present in BLAST or FASTA reports, facilitating extraction of useful information of large sets of sequence alignments, and allowing for easy handling and processing of the data. AVAILABILITY: BioParser is licensed under the Creative Commons Attribution-NonCommercial-NoDerivs 2.0 license terms (http://creativecommons.org/licenses/by-nc-nd/2.0/) and is available upon request. Additional information can be found at the BioParser website (http://www.dbbm.fiocruz.br/BioParser.html).

Squid - a simple bioinformatics grid.

BACKGROUND: BLAST is a widely used genetic research tool for analysis of similarity between nucleotide and protein sequences. This paper presents a software application entitled "Squid" that makes use of grid technology. The current version, as an example, is configured for BLAST applications, but adaptation for other computing intensive repetitive tasks can be easily accomplished in the open source version. This enables the allocation of remote resources to perform distributed computing, making large BLAST queries viable without the need of high-end computers. RESULTS: Most distributed computing / grid solutions have complex installation procedures requiring a computer specialist, or have limitations regarding operating systems. Squid is a multi-platform, open-source program designed to "keep things simple" while offering high-end computing power for large scale applications. Squid also has an efficient fault tolerance and crash recovery system against data loss, being able to re-route jobs upon node failure and recover even if the master machine fails. Our results show that a Squid application, working with N nodes and proper network resources, can process BLAST queries almost N times faster than if working with only one computer. CONCLUSION: Squid offers high-end computing, even for the non-specialist, and is freely available at the project web site. Its open-source and binary Windows distributions contain detailed instructions and a "plug-n-play" instalation containing a pre-configured example.

Surgical management of cystic dilatation bile ducts in adults.

BACKGROUND: The cystic dilatation of the biliary tract is a rare disease and uncertain origin. It is recognized more frequently in children; however, its incidence comes increasing in adults, representing 20% of the cases. AIM: To evaluate morbimortality rates, evolution and handing of patients with cystic dilatation bile ducts in adults. METHODS: Were evaluated, retrospectively, five adults who had the diagnosis of choledochal cyst and that had been submitted to some surgical procedure. RESULTS: Abdominal pain was the commonest complain to all patients. Jaundice was present in 80%. Ultrasound scanning was done in all the cases as initial examination. CT scan, magnetic resonance imaging and endoscopic retrograde cholangiopancreatography were also done in some patients; however, the diagnosis was established intra-operatively in all cases. The cyst resection with reconstruction of the biliary tract was done in 60%; the cystojejunostomy in 20%; and in 20% biliary tract drainage. CONCLUSIONS: Biliary tract cystic dilatation is a rare disease. However, its incidence is increasing in the adult population, so, it must be thought as differential diagnosis when facing obstructive jaundice.

Genome Comparison of an Ancestral Isolate and a Modern Isolate of Mycobacterium tuberculosis of the Beijing Lineage from São Paulo, Brazil.

Mycobacterium tuberculosis of the Bejing subtype (MtbB) is transmitted efficiently in high burden countries for this genotype. A higher virulence was associated with isolates of the "modern" Beijing genotype sub-lineages when compared to "ancient" ones. Here, we report the full genomes of the strain representing these two genotypes from Brazil, a country with a low incidence of MtbB.

Identification of proteins in Streptococcus pneumoniae by reverse vaccinology and genetic diversity of these proteins in clinical isolates.

Streptococcus pneumoniae is a major cause of morbidity and mortality worldwide. Virulence-associated proteins common and conserved among all capsular types now represent the best strategy to combat pneumococcal infections. Our aim was to identify conserved targets in pneumococci that showed positive prediction for lipoprotein and extracellular subcellular location using bioinformatics programs and verify the distribution and the degree of conservation of these targets in pneumococci. These targets can be considered potential vaccine candidate to be evaluated in the future. A set of 13 targets were analyzed and confirmed the presence in all pneumococci tested. These 13 genes were highly conserved showing around >96 % of amino acid and nucleotide identity, but they were also present and show high identity in the closely related species Streptococcus mitis, Streptococcus oralis, and Streptococcus pseudopneumoniae. S. oralis clusters away from S. pneumoniae, while S. pseudopneumoniae and S. mitis cluster closer. The divergence between the selected targets was too small to be observed consistently in phylogenetic groups between the analyzed genomes of S. pneumoniae. The proteins analyzed fulfill two of the initial criteria of a vaccine candidate: targets are present in a variety of different pneumococci strains including different serotypes and are conserved among the samples evaluated.

The Purine Bias of Coding Sequences is Determined by Physicochemical Constraints on Proteins.

For this report, we analyzed protein secondary structures in relation to the statistics of three nucleotide codon positions. The purpose of this investigation was to find which properties of the ribosome, tRNA or protein level, could explain the purine bias (Rrr) as it is observed in coding DNA. We found that the Rrr pattern is the consequence of a regularity (the codon structure) resulting from physicochemical constraints on proteins and thermodynamic constraints on ribosomal machinery. The physicochemical constraints on proteins mainly come from the hydropathy and molecular weight (MW) of secondary structures as well as the energy cost of amino acid synthesis. These constraints appear through a network of statistical correlations, such as (i) the cost of amino acid synthesis, which is in favor of a higher level of guanine in the first codon position, (ii) the constructive contribution of hydropathy alternation in proteins, (iii) the spatial organization of secondary structure in proteins according to solvent accessibility, (iv) the spatial organization of secondary structure according to amino acid hydropathy, (v) the statistical correlation of MW with protein secondary structures and their overall hydropathy, (vi) the statistical correlation of thymine in the second codon position with hydropathy and the energy cost of amino acid synthesis, and (vii) the statistical correlation of adenine in the second codon position with amino acid complexity and the MW of secondary protein structures. Amino acid physicochemical properties and functional constraints on proteins constitute a code that is translated into a purine bias within the coding DNA via tRNAs. In that sense, the Rrr pattern within coding DNA is the effect of information transfer on nucleotide composition from protein to DNA by selection according to the codon positions. Thus, coding DNA structure and ribosomal machinery co-evolved to minimize the energy cost of protein coding given the functional constraints on proteins.

The Comparative Genomics and Phylogenomics of Leishmania amazonensis Parasite.

Leishmaniasis is an infectious disease caused by Leishmania species. Leishmania amazonensis is a New World Leishmania species belonging to the Mexicana complex, which is able to cause all types of leishmaniasis infections. The L. amazonensis reference strain MHOM/BR/1973/M2269 was sequenced identifying 8,802 codifying sequences (CDS), most of them of hypothetical function. Comparative analysis using six Leishmania species showed a core set of 7,016 orthologs. L. amazonensis and Leishmania mexicana share the largest number of distinct orthologs, while Leishmania braziliensis presented the largest number of inparalogs. Additionally, phylogenomic analysis confirmed the taxonomic position for L. amazonensis within the "Mexicana complex", reinforcing understanding of the split of New and Old World Leishmania. Potential non-homologous isofunctional enzymes (NISE) were identified between L. amazonensis and Homo sapiens that could provide new drug targets for development.