Phylogenomics-based reconstruction of protozoan species tree

Commensal Intestinal Protozoa-Underestimated Members of the Gut Microbial Community

The human gastrointestinal microbiota contains a diverse consortium of microbes, including bacteria, protozoa, viruses, and fungi. Through millennia of co-evolution, the host-microbiota interactions have shaped the immune system to both tolerate and maintain the symbiotic relationship with commensal microbiota, while exerting protective responses against invading pathogens. Microbiome research is dominated by studies describing the impact of prokaryotic bacteria on gut immunity with a limited understanding of their relationship with other integral microbiota constituents. However, converging evidence shows that eukaryotic organisms, such as commensal protozoa, can play an important role in modulating intestinal immune responses as well as influencing the overall health of the host. The presence of several protozoa species has recently been shown to be a common occurrence in healthy populations worldwide, suggesting that many of these are commensals rather than invading pathogens. This review aims to discuss the most recent, conflicting findings regarding the role of intestinal protozoa in gut homeostasis, interactions between intestinal protozoa and the bacterial microbiota, as well as potential immunological consequences of protozoa colonization.

Genomic analyses reveal moderate levels of ploidy, high heterozygosity and structural variations in a Colombian isolate of Leishmania (Leishmania) amazonensis

Leishmania amazonensis is one of the causative agents of the different forms of cutaneous leishmaniasis present in Latin America. This species has been isolated from humans and animals (canine/feline) in some endemic regions of Colombia. Therefore, L. amazonensis is of great relevance at the clinical and epidemiological levels in medicine and veterinary science. Until now, very few genomes from this species are available. Here, we report the complete genome sequence of a laboratory-adapted L. amazonensis strain isolated from a human patient with clinical manifestation of cutaneous leishmaniasis in Colombia. The genome sequence not only allowed inter and intra-species comparative analyses to be performed with the sequenced genomes of L. amazonensis strains from different geographical regions, but also increased our knowledge about the genomic behavior of this L. amazonensis Colombian strain. This isolate was also characterized in terms of single nucleotide polymorphisms, chromosome and gene copy number variations (CNVs). The results revealed moderate aneuploidy, CNVs in genes involved in the virulence, growth, and survival of the parasite, and in the distributions of the multicopy genes on some chromosomes, as well as a high level of heterozygosity. The data confirmed previous reports that identified unique variations in L. amazonensis, suggesting aneuploidy may not have a high fitness cost and may allow the rapid generation of diversity in Leishmania parasites growing normally.

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.

An orthology-based analysis of pathogenic protozoa impacting global health: an improved comparative genomics approach with prokaryotes and model eukaryote orthologs

A key focus in 21(st) century integrative biology and drug discovery for neglected tropical and other diseases has been the use of BLAST-based computational methods for identification of orthologous groups in pathogenic organisms to discern orthologs, with a view to evaluate similarities and differences among species, and thus allow the transfer of annotation from known/curated proteins to new/non-annotated ones. We used here a profile-based sensitive methodology to identify distant homologs, coupled to the NCBI's COG (Unicellular orthologs) and KOG (Eukaryote orthologs), permitting us to perform comparative genomics analyses on five protozoan genomes. OrthoSearch was used in five protozoan proteomes showing that 3901 and 7473 orthologs can be identified by comparison with COG and KOG proteomes, respectively. The core protozoa proteome inferred was 418 Protozoa-COG orthologous groups and 704 Protozoa-KOG orthologous groups: (i) 31.58% (132/418) belongs to the category J (translation, ribosomal structure, and biogenesis), and 9.81% (41/418) to the category O (post-translational modification, protein turnover, chaperones) using COG; (ii) 21.45% (151/704) belongs to the categories J, and 13.92% (98/704) to the O using KOG. The phylogenomic analysis showed four well-supported clades for Eukarya, discriminating Multicellular [(i) human, fly, plant and worm] and Unicellular [(ii) yeast, (iii) fungi, and (iv) protozoa] species. These encouraging results attest to the usefulness of the profile-based methodology for comparative genomics to accelerate semi-automatic re-annotation, especially of the protozoan proteomes. This approach may also lend itself for applications in global health, for example, in the case of novel drug target discovery against pathogenic organisms previously considered difficult to research with traditional drug discovery tools.

Phylogeny of bacterial and archaeal genomes using conserved genes: supertrees and supermatrices

Over 3000 microbial (bacterial and archaeal) genomes have been made publically available to date, providing an unprecedented opportunity to examine evolutionary genomic trends and offering valuable reference data for a variety of other studies such as metagenomics. The utility of these genome sequences is greatly enhanced when we have an understanding of how they are phylogenetically related to each other. Therefore, we here describe our efforts to reconstruct the phylogeny of all available bacterial and archaeal genomes. We identified 24, single-copy, ubiquitous genes suitable for this phylogenetic analysis. We used two approaches to combine the data for the 24 genes. First, we concatenated alignments of all genes into a single alignment from which a Maximum Likelihood (ML) tree was inferred using RAxML. Second, we used a relatively new approach to combining gene data, Bayesian Concordance Analysis (BCA), as implemented in the BUCKy software, in which the results of 24 single-gene phylogenetic analyses are used to generate a "primary concordance" tree. A comparison of the concatenated ML tree and the primary concordance (BUCKy) tree reveals that the two approaches give similar results, relative to a phylogenetic tree inferred from the 16S rRNA gene. After comparing the results and the methods used, we conclude that the current best approach for generating a single phylogenetic tree, suitable for use as a reference phylogeny for comparative analyses, is to perform a maximum likelihood analysis of a concatenated alignment of conserved, single-copy genes.