Horizontal gene transfer in eukaryotic parasites: a case study of Entamoeba histolytica and Trichomonas vaginalis

Phylogenetic Analysis of Pyruvate-Ferredoxin Oxidoreductase, a Redox Enzyme Involved in the Pharmacological Activation of Nitro-Based Prodrugs in Bacteria and Protozoa

The present frontrunners in the chemotherapy of infections caused by protozoa are nitro-based prodrugs that are selectively activated by PFOR-mediated redox reactions. This study seeks to analyze the distribution of PFOR in selected protozoa and bacteria by applying comparative genomics to test the hypothesis that PFOR in eukaryotes was acquired through horizontal gene transfer (HGT) from bacteria. Furthermore, to identify other putatively acquired genes, proteome-wide and gene enrichment analyses were used. A plausible explanation for the patchy occurrence of PFOR in protozoa is based on the hypothesis that bacteria are potential sources of genes that enhance the adaptation of protozoa in hostile environments. Comparative genomics of Entamoeba histolytica and the putative gene donor, Desulfovibrio vulgaris, identified eleven candidate genes for HGT involved in intermediary metabolism. If these results can be reproduced in other PFOR-possessing protozoa, it would provide more validated evidence to support the horizontal transfer of pfor from bacteria.

Pathogenicity and virulence of Entamoeba histolytica, the agent of amoebiasis

The amoeba parasite Entamoeba histolytica is the causative agent of human amebiasis, an enteropathic disease affecting millions of people worldwide. This ancient protozoan is an elementary example of how parasites evolve with humans, e.g. taking advantage of multiple mechanisms to evade immune responses, interacting with microbiota for nutritional and protective needs, utilizing host resources for growth, division, and encystation. These skills of E. histolytica perpetuate the species and incidence of infection. However, in 10% of infected cases, the parasite turns into a pathogen; the host-parasite equilibrium is then disorganized, and the simple lifecycle based on two cell forms, trophozoites and cysts, becomes unbalanced. Trophozoites acquire a virulent phenotype which, when non-controlled, leads to intestinal invasion with the onset of amoebiasis symptoms. Virulent E. histolytica must cross mucus, epithelium, connective tissue and possibly blood. This highly mobile parasite faces various stresses and a powerful host immune response, with oxidative stress being a challenge for its survival. New emerging research avenues and omics technologies target gene regulation to determine human or parasitic factors activated upon infection, their role in virulence activation, and in pathogenesis; this research bears in mind that E. histolytica is a resident of the complex intestinal ecosystem. The goal is to eradicate amoebiasis from the planet, but the parasitic life of E. histolytica is ancient and complex and will likely continue to evolve with humans. Advances in these topics are summarized here.

Presence of CRISPR CAS-Like Sequences as a Proposed Mechanism for Horizontal Genetic Exchanges between Trichomonas vaginalis and Its Associated Virus: A Comparative Genomic Analysis with the First Report of a Putative CRISPR CAS Structures in Eukaryotic Cells

Introduction

Trichomonas vaginalis genome is among the largest genome size and coding capacities. Combinations of gene duplications, transposon, repeated sequences, and lateral gene transfers (LGTs) have contributed to the unexpected large genomic size and diversity. This study is aimed at investigating genomic exchange and seeking for presence of the CRISPR CAS system as one of the possible mechanisms for some level of genetic exchange. Material and Methods. In this comparative analysis, 398 publicly available Trichomonas vaginalis complete genomes were investigated for the presence of CRISPR CAS. Spacer sequences were also analyzed for their origin using BLAST.

Results

We identified a CRISPR CAS (Cas3). CRISPR spacers are highly similar to transposable genetic elements such as viruses of protozoan parasites, especially megavirals, some transposons, and, interestingly, papillomavirus and HIV-1 in a few cases. Discussion. There is a striking similarity between the prokaryotes/Archaean CRISPR and what we find as eukaryotic CRISPR. About 5-10% of the 398 T. vaginalis possess a CRISPR structure.

Conclusion

According to sequences and their organization, we assume that these repeated sequences and spacer, along with their mentioned features, could be the eukaryotic homolog of prokaryotes and Archaean CRISPR systems and may involve in a process similar to the CRISPR function.

Old genes in new places: A taxon-rich analysis of interdomain lateral gene transfer events

Vertical inheritance is foundational to Darwinian evolution, but fails to explain major innovations such as the rapid spread of antibiotic resistance among bacteria and the origin of photosynthesis in eukaryotes. While lateral gene transfer (LGT) is recognized as an evolutionary force in prokaryotes, the role of LGT in eukaryotic evolution is less clear. With the exception of the transfer of genes from organelles to the nucleus, a process termed endosymbiotic gene transfer (EGT), the extent of interdomain transfer from prokaryotes to eukaryotes is highly debated. A common critique of studies of interdomain LGT is the reliance on the topology of single-gene trees that attempt to estimate more than one billion years of evolution. We take a more conservative approach by identifying cases in which a single clade of eukaryotes is found in an otherwise prokaryotic gene tree (i.e. exclusive presence). Starting with a taxon-rich dataset of over 13,600 gene families and passing data through several rounds of curation, we identify and categorize the function of 306 interdomain LGT events into diverse eukaryotes, including 189 putative EGTs, 52 LGTs into Opisthokonta (i.e. animals, fungi and their microbial relatives), and 42 LGTs nearly exclusive to anaerobic eukaryotes. To assess differential gene loss as an explanation for exclusive presence, we compare branch lengths within each LGT tree to a set of vertically-inherited genes subsampled to mimic gene loss (i.e. with the same taxonomic sampling) and consistently find shorter relative distance between eukaryotes and prokaryotes in LGT trees, a pattern inconsistent with gene loss. Our methods provide a framework for future studies of interdomain LGT and move the field closer to an understanding of how best to model the evolutionary history of eukaryotes.

Identification and Functional Characterization of Divergent 3'-Phosphate tRNA Ligase From Entamoeba histolytica

HSPC117/RtcB, 3'-phosphate tRNA ligase, is a critical enzyme involved in tRNA splicing and maturation. HSPC117/RtcB is also involved in mRNA splicing of some protein-coding genes including XBP-1. Entamoeba histolytica, a protozoan parasite responsible for human amebiasis, possesses two RtcB proteins (EhRtcB1 and 2), but their biological functions remain unknown. Both RtcBs show kinship with mammalian/archaeal type, and all amino acid residues present in the active sites are highly conserved, as suggested by protein alignment and phylogenetic analyses. EhRtcB1 was demonstrated to be localized to the nucleus, while EhRtcB2 was in the cytosol. EhRtcB1, but not EhRtcB2, was required for optimal growth of E. histolytica trophozoites. Both EhRtcB1 (in cooperation with EhArchease) and EhRtcB2 showed RNA ligation activity in vitro. The predominant role of EhRtcB1 in tRNA^Ile(UAU) processing in vivo was demonstrated in EhRtcB1- and 2-gene silenced strains. Taken together, we have demonstrated the conservation of tRNA splicing and functional diversification of RtcBs in this amoebozoan lineage.

The Mastigamoeba balamuthi Genome and the Nature of the Free-Living Ancestor of Entamoeba

The transition of free-living organisms to parasitic organisms is a mysterious process that occurs in all major eukaryotic lineages. Parasites display seemingly unique features associated with their pathogenicity; however, it is important to distinguish ancestral preconditions to parasitism from truly new parasite-specific functions. Here, we sequenced the genome and transcriptome of anaerobic free-living Mastigamoeba balamuthi and performed phylogenomic analysis of four related members of the Archamoebae, including Entamoeba histolytica, an important intestinal pathogen of humans. We aimed to trace gene histories throughout the adaptation of the aerobic ancestor of Archamoebae to anaerobiosis and throughout the transition from a free-living to a parasitic lifestyle. These events were associated with massive gene losses that, in parasitic lineages, resulted in a reduction in structural features, complete losses of some metabolic pathways, and a reduction in metabolic complexity. By reconstructing the features of the common ancestor of Archamoebae, we estimated preconditions for the evolution of parasitism in this lineage. The ancestor could apparently form chitinous cysts, possessed proteolytic enzyme machinery, compartmentalized the sulfate activation pathway in mitochondrion-related organelles, and possessed the components for anaerobic energy metabolism. After the split of Entamoebidae, this lineage gained genes encoding surface membrane proteins that are involved in host–parasite interactions. In contrast, gene gains identified in the M. balamuthi lineage were predominantly associated with polysaccharide catabolic processes. A phylogenetic analysis of acquired genes suggested an essential role of lateral gene transfer in parasite evolution (Entamoeba) and in adaptation to anaerobic aquatic sediments (Mastigamoeba).

Characterization of the BspA and Pmp protein family of trichomonads

Background

Trichomonas vaginalis is a human-infecting trichomonad and as such the best studied and the only for which the full genome sequence is available considering its parasitic lifestyle, T. vaginalis encodes an unusually high number of proteins. Many gene families are massively expanded and some genes are speculated to have been acquired from prokaryotic sources. Among the latter are two gene families that harbour domains which share similarity with proteins of Bacteroidales/Spirochaetales and Chlamydiales: the BspA and the Pmp proteins, respectively.

Results

We sequenced the transcriptomes of five trichomonad species and screened for the presence of BspA and Pmp domain-containing proteins and characterized individual candidate proteins from both families in T. vaginalis. Here, we demonstrate that (i) BspA and Pmp domain-containing proteins are universal to trichomonads, but specifically expanded in T. vaginalis; (ii) in line with a concurrent expansion of the endocytic machinery, there is a high number of BspA and Pmp proteins which carry C-terminal endocytic motifs; and (iii) both families traffic through the ER and have the ability to increase adhesion performance in a non-virulent T. vaginalis strain and Tetratrichomonas gallinarum by a so far unknown mechanism.

Conclusions

Our results initiate the functional characterization of these two broadly distributed protein families and help to better understand the origin and evolution of BspA and Pmp domains in trichomonads.

Electronic supplementary material

The online version of this article (10.1186/s13071-019-3660-z) contains supplementary material, which is available to authorized users.

Redox Pathways as Drug Targets in Microaerophilic Parasites

The microaerophilic parasites Entamoeba histolytica, Trichomonas vaginalis, and Giardia lamblia jointly cause hundreds of millions of infections in humans every year. Other microaerophilic parasites such as Tritrichomonas foetus and Spironucleus spp. pose a relevant health problem in veterinary medicine. Unfortunately, vaccines against these pathogens are unavailable, but their microaerophilic lifestyle opens opportunities for specifically developed chemotherapeutics. In particular, their high sensitivity towards oxygen can be exploited by targeting redox enzymes. This review focusses on the redox pathways of microaerophilic parasites and on drugs, either already in use or currently in the state of development, which target these pathways.

Genomic Data Quality Impacts Automated Detection of Lateral Gene Transfer in Fungi

Lateral gene transfer (LGT, also known as horizontal gene transfer), an atypical mechanism of transferring genes between species, has almost become the default explanation for genes that display an unexpected composition or phylogeny. Numerous methods of detecting LGT events all rely on two fundamental strategies: primary structure composition or gene tree/species tree comparisons. Discouragingly, the results of these different approaches rarely coincide. With the wealth of genome data now available, detection of laterally transferred genes is increasingly being attempted in large uncurated eukaryotic datasets. However, detection methods depend greatly on the quality of the underlying genomic data, which are typically complex for eukaryotes. Furthermore, given the automated nature of genomic data collection, it is typically impractical to manually verify all protein or gene models, orthology predictions, and multiple sequence alignments, requiring researchers to accept a substantial margin of error in their datasets. Using a test case comprising plant-associated genomes across the fungal kingdom, this study reveals that composition- and phylogeny-based methods have little statistical power to detect laterally transferred genes. In particular, phylogenetic methods reveal extreme levels of topological variation in fungal gene trees, the vast majority of which show departures from the canonical species tree. Therefore, it is inherently challenging to detect LGT events in typical eukaryotic genomes. This finding is in striking contrast to the large number of claims for laterally transferred genes in eukaryotic species that routinely appear in the literature, and questions how many of these proposed examples are statistically well supported.

Lateral gene transfer of p-cresol- and indole-producing enzymes from environmental bacteria to Mastigamoeba balamuthi

p-Cresol and indole are volatile biologically active products of the bacterial degradation of tyrosine and tryptophan respectively. They are typically produced by bacteria in animal intestines, soil and various sediments. Here, we demonstrate that the free-living eukaryote Mastigamoeba balamuthi and its pathogenic relative Entamoeba histolytica produce significant amounts of indole via tryptophanase activity. Unexpectedly, M. balamuthi also produces p-cresol in concentrations that are bacteriostatic to non-p-cresol-producing bacteria. The ability of M. balamuthi to produce p-cresol, which has not previously been observed in any eukaryotic microbe, was gained due to the lateral acquisition of a bacterial gene for 4-hydroxyphenylacetate decarboxylase (HPAD). In bacteria, the genes for HPAD and the S-adenosylmethionine-dependent activating enzyme (AE) are present in a common operon. In M. balamuthi, HPAD displays a unique fusion with the AE that suggests the operon-mediated transfer of genes from a bacterial donor. We also clarified that the tyrosine-to-4-hydroxyphenylacetate conversion proceeds via the Ehrlich pathway. The acquisition of the bacterial HPAD gene may provide M. balamuthi a competitive advantage over other microflora in its native habitat.

The E. histolytica Genome Structure and Virulence

The outcome of an Entamoeba histolytica infection is variable and the contribution of genetic diversity within E. histolytica to human disease is not fully understood. The information provided by the whole genome sequence of the E. histolytica reference laboratory strain (HM-1:IMSS) and thirteen additional laboratory strains have been made publically available. In this review theories on the source of the unexpected level of structural diversity found in E. histolytica will be discussed.

Bulky Trichomonad Genomes: Encoding a Swiss Army Knife

The trichomonads are a remarkably successful lineage of ancient, predominantly parasitic protozoa. Recent molecular analyses have revealed extensive duplication of certain genetic loci in trichomonads. Consequently, their genomes are exceptionally large compared to other parasitic protozoa. Retention of these large gene expansions across different trichomonad families raises the question: do these duplications afford an advantage? Many duplicated genes are linked to the parasitic lifestyle and some are regulated differently to their paralogues, suggesting they have acquired new functions. It is proposed that these large genomes encode a Swiss army knife of sorts, packed with a multitude of tools for use in many different circumstances. This may have bestowed trichomonads with the extraordinary versatility that has undoubtedly contributed to their success.

Nothing in Evolution Makes Sense Except in the Light of Genomics: Read-Write Genome Evolution as an Active Biological Process

The 21st century genomics-based analysis of evolutionary variation reveals a number of novel features impossible to predict when Dobzhansky and other evolutionary biologists formulated the neo-Darwinian Modern Synthesis in the middle of the last century. These include three distinct realms of cell evolution; symbiogenetic fusions forming eukaryotic cells with multiple genome compartments; horizontal organelle, virus and DNA transfers; functional organization of proteins as systems of interacting domains subject to rapid evolution by exon shuffling and exonization; distributed genome networks integrated by mobile repetitive regulatory signals; and regulation of multicellular development by non-coding lncRNAs containing repetitive sequence components. Rather than single gene traits, all phenotypes involve coordinated activity by multiple interacting cell molecules. Genomes contain abundant and functional repetitive components in addition to the unique coding sequences envisaged in the early days of molecular biology. Combinatorial coding, plus the biochemical abilities cells possess to rearrange DNA molecules, constitute a powerful toolbox for adaptive genome rewriting. That is, cells possess "Read-Write Genomes" they alter by numerous biochemical processes capable of rapidly restructuring cellular DNA molecules. Rather than viewing genome evolution as a series of accidental modifications, we can now study it as a complex biological process of active self-modification.

Recent events dominate interdomain lateral gene transfers between prokaryotes and eukaryotes and, with the exception of endosymbiotic gene transfers, few ancient transfer events persist

While there is compelling evidence for the impact of endosymbiotic gene transfer (EGT; transfer from either mitochondrion or chloroplast to the nucleus) on genome evolution in eukaryotes, the role of interdomain transfer from bacteria and/or archaea (i.e. prokaryotes) is less clear. Lateral gene transfers (LGTs) have been argued to be potential sources of phylogenetic information, particularly for reconstructing deep nodes that are difficult to recover with traditional phylogenetic methods. We sought to identify interdomain LGTs by using a phylogenomic pipeline that generated 13 465 single gene trees and included up to 487 eukaryotes, 303 bacteria and 118 archaea. Our goals include searching for LGTs that unite major eukaryotic clades, and describing the relative contributions of LGT and EGT across the eukaryotic tree of life. Given the difficulties in interpreting single gene trees that aim to capture the approximately 1.8 billion years of eukaryotic evolution, we focus on presence-absence data to identify interdomain transfer events. Specifically, we identify 1138 genes found only in prokaryotes and representatives of three or fewer major clades of eukaryotes (e.g. Amoebozoa, Archaeplastida, Excavata, Opisthokonta, SAR and orphan lineages). The majority of these genes have phylogenetic patterns that are consistent with recent interdomain LGTs and, with the notable exception of EGTs involving photosynthetic eukaryotes, we detect few ancient interdomain LGTs. These analyses suggest that LGTs have probably occurred throughout the history of eukaryotes, but that ancient events are not maintained unless they are associated with endosymbiotic gene transfer among photosynthetic lineages.

Horizontal Gene Transfers from Bacteria to Entamoeba Complex: A Strategy for Dating Events along Species Divergence

Horizontal gene transfer has proved to be relevant in eukaryotic evolution, as it has been found more often than expected and related to adaptation to certain niches. A relatively large list of laterally transferred genes has been proposed and evaluated for the parasite Entamoeba histolytica. The goals of this work were to elucidate the importance of lateral gene transfer along the evolutionary history of some members of the genus Entamoeba, through identifying donor groups and estimating the divergence time of some of these events. In order to estimate the divergence time of some of the horizontal gene transfer events, the dating of some Entamoeba species was necessary, following an indirect dating strategy based on the fossil record of plausible hosts. The divergence between E. histolytica and E. nuttallii probably occurred 5.93 million years ago (Mya); this lineage diverged from E. dispar 9.97 Mya, while the ancestor of the latter separated from E. invadens 68.18 Mya. We estimated times for 22 transferences; the most recent occurred 31.45 Mya and the oldest 253.59 Mya. Indeed, the acquisition of genes through lateral transfer may have triggered a period of adaptive radiation, thus playing a major role in the evolution of the Entamoeba genus.

Gene expression profiling gut microbiota in different races of humans

The gut microbiome is shaped and modified by the polymorphisms of microorganisms in the intestinal tract. Its composition shows strong individual specificity and may play a crucial role in the human digestive system and metabolism. Several factors can affect the composition of the gut microbiome, such as eating habits, living environment, and antibiotic usage. Thus, various races are characterized by different gut microbiome characteristics. In this present study, we studied the gut microbiomes of three different races, including individuals of Asian, European and American races. The gut microbiome and the expression levels of gut microbiome genes were analyzed in these individuals. Advanced feature selection methods (minimum redundancy maximum relevance and incremental feature selection) and four machine-learning algorithms (random forest, nearest neighbor algorithm, sequential minimal optimization, Dagging) were employed to capture key differentially expressed genes. As a result, sequential minimal optimization was found to yield the best performance using the 454 genes, which could effectively distinguish the gut microbiomes of different races. Our analyses of extracted genes support the widely accepted hypotheses that eating habits, living environments and metabolic levels in different races can influence the characteristics of the gut microbiome.

Analysis of the Bacterial Diversity in Liver Abscess: Differences Between Pyogenic and Amebic Abscesses

Several recent studies have demonstrated that virulence in Entamoeba histolytica is triggered in the presence of both pathogenic and nonpathogenic bacteria species using in vitro and in vivo experimental animal models. In this study, we examined samples aspirated from abscess material obtained from patients who were clinically diagnosed with amebic liver abscess (ALA) or pyogenic liver abscess (PLA). To determine the diversity of bacterial species in the abscesses, we performed partial 16S rRNA gene sequencing. In addition, the E. histolytica and Entamoeba dispar species were genotyped using tRNA-linked short tandem repeats as specific molecular markers. The association between clinical data and bacterial and parasite genotypes were examined through a correspondence analysis. The results showed the presence of numerous bacterial groups. These taxonomic groups constitute common members of the gut microbiota, although all of the detected bacterial species have a close phylogenetic relationship with bacterial pathogens. Furthermore, some patients clinically diagnosed with PLA and ALA were coinfected with E. dispar or E. histolytica, which suggests that the virulence of these parasites increased in the presence of bacteria. However, no specific bacterial groups were associated with this effect. Together, our results suggest a nonspecific mechanism of virulence modulation by bacteria in Entamoeba.

The Transcriptome Sequence of Dientamoeba fragilis Offers New Biological Insights on its Metabolism, Kinome, Degradome and Potential Mechanisms of Pathogenicity

Dientamoeba fragilis is a human bowel parasite with a worldwide distribution. Dientamoeba was once described as a rare and harmless commensal though recent reports suggest it is common and potentially pathogenic. Molecular data on Dientamoeba is scarce which limits our understanding of this parasite. To address this, sequencing of the Dientamoeba transcriptome was performed. Messenger RNA was extracted from cultured Dientamoeba trophozoites originating from clinical stool specimens, and sequenced using Roche GS FLX and Illumina HiSeq technologies. In total 6,595 Dientamoeba transcripts were identified. These sequences were analysed using the BLAST2GO software suite and via BLAST comparisons to sequences available from TrichDB, GenBank, MEROPS and kinase.com. Several novel KEGG pathway maps were generated and gene ontology analysis was also performed. These results are thoroughly discussed guided by knowledge available for other related protozoa. Attention is paid to the novel biological insights afforded by this data including peptidases and kinases of Dientamoeba, as well as its metabolism, novel chemotherapeutics and possible mechanisms of pathogenicity. Currently, this work represents the largest contribution to our understanding of Dientamoeba molecular biology and also represents a major contribution to our understanding of the trichomonads generally, many of which are important pathogens of humans and animals.

The FAD-dependent glycerol-3-phosphate dehydrogenase of Giardia duodenalis: an unconventional enzyme that interacts with the g14-3-3 and it is a target of the antitumoral compound NBDHEX

The flagellated protozoan Giardia duodenalis is a worldwide parasite causing giardiasis, an acute and chronic diarrheal disease. Metabolism in G. duodenalis has a limited complexity thus making metabolic enzymes ideal targets for drug development. However, only few metabolic pathways (i.e., carbohydrates) have been described so far. Recently, the parasite homolog of the mitochondrial-like glycerol-3-phosphate dehydrogenase (gG3PD) has been identified among the interactors of the g14-3-3 protein. G3PD is involved in glycolysis, electron transport, glycerophospholipids metabolism, and hyperosmotic stress response, and is emerging as promising target in tumor treatment. In this work, we demonstrate that gG3PD is a functional flavoenzyme able to convert glycerol-3-phosphate into dihydroxyacetone phosphate and that its activity and the intracellular glycerol level increase during encystation. Taking advantage of co-immunoprecipitation assays and deletion mutants, we provide evidence that gG3PD and g14-3-3 interact at the trophozoite stage, the intracellular localization of gG3PD is stage dependent and it partially co-localizes with mitosomes during cyst development. Finally, we demonstrate that the gG3PD activity is affected by the antitumoral compound 6-(7-nitro-2,1,3-benzoxadiazol-4-ylthio)hexanol, that results more effective in vitro at killing G. duodenalis trophozoites than the reference drug metronidazole. Overall, our results highlight the involvement of gG3PD in processes crucial for the parasite survival thus proposing this enzyme as target for novel antigiardial interventions.

The parasite Trichomonas vaginalis expresses thousands of pseudogenes and long non-coding RNAs independently from functional neighbouring genes

Background

The human pathogen Trichomonas vaginalis is a parabasalian flagellate that is estimated to infect 3% of the world’s population annually. With a 160 megabase genome and up to 60,000 genes residing in six chromosomes, the parasite has the largest genome among sequenced protists. Although it is thought that the genome size and unusual large coding capacity is owed to genome duplication events, the exact reason and its consequences are less well studied.

Results

Among transcriptome data we found thousands of instances, in which reads mapped onto genomic loci not annotated as genes, some reaching up to several kilobases in length. At first sight these appear to represent long non-coding RNAs (lncRNAs), however, about half of these lncRNAs have significant sequence similarities to genomic loci annotated as protein-coding genes. This provides evidence for the transcription of hundreds of pseudogenes in the parasite. Conventional lncRNAs and pseudogenes are expressed in Trichomonas through their own transcription start sites and independently from flanking genes in Trichomonas. Expression of several representative lncRNAs was verified through reverse-transcriptase PCR in different T. vaginalis strains and case studies exclude the use of alternative start codons or stop codon suppression for the genes analysed.

Conclusion

Our results demonstrate that T. vaginalis expresses thousands of intergenic loci, including numerous transcribed pseudogenes. In contrast to yeast these are expressed independently from neighbouring genes. Our results furthermore illustrate the effect genome duplication events can have on the transcriptome of a protist. The parasite’s genome is in a steady state of changing and we hypothesize that the numerous lncRNAs could offer a large pool for potential innovation from which novel proteins or regulatory RNA units could evolve.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-906) contains supplementary material, which is available to authorized users.

Phylogenomic study indicates widespread lateral gene transfer in Entamoeba and suggests a past intimate relationship with parabasalids

Lateral gene transfer (LGT) has impacted the evolutionary history of eukaryotes, though to a lesser extent than in bacteria and archaea. Detecting LGT and distinguishing it from single gene tree artifacts is difficult, particularly when considering very ancient events (i.e., over hundreds of millions of years). Here, we use two independent lines of evidence--a taxon-rich phylogenetic approach and an assessment of the patterns of gene presence/absence--to evaluate the extent of LGT in the parasitic amoebozoan genus Entamoeba. Previous work has suggested that a number of genes in the genome of Entamoeba spp. were acquired by LGT. Our approach, using an automated phylogenomic pipeline to build taxon-rich gene trees, suggests that LGT is more extensive than previously thought. Our analyses reveal that genes have frequently entered the Entamoeba genome via nonvertical events, including at least 116 genes acquired directly from bacteria or archaea, plus an additional 22 genes in which Entamoeba plus one other eukaryote are nested among bacteria and/or archaea. These genes may make good candidates for novel therapeutics, as drugs targeting these genes are less likely to impact the human host. Although we recognize the challenges of inferring intradomain transfers given systematic errors in gene trees, we find 109 genes supporting LGT from a eukaryote to Entamoeba spp., and 178 genes unique to Entamoeba spp. and one other eukaryotic taxon (i.e., presence/absence data). Inspection of these intradomain LGTs provide evidence of a common sister relationship between genes of Entamoeba (Amoebozoa) and parabasalids (Excavata). We speculate that this indicates a past close relationship (e.g., symbiosis) between ancestors of these extant lineages.

A recently transferred cluster of bacterial genes in Trichomonas vaginalis--lateral gene transfer and the fate of acquired genes

Background

Lateral Gene Transfer (LGT) has recently gained recognition as an important contributor to some eukaryote proteomes, but the mechanisms of acquisition and fixation in eukaryotic genomes are still uncertain. A previously defined norm for LGTs in microbial eukaryotes states that the majority are genes involved in metabolism, the LGTs are typically localized one by one, surrounded by vertically inherited genes on the chromosome, and phylogenetics shows that a broad collection of bacterial lineages have contributed to the transferome.

Results

A unique 34 kbp long fragment with 27 clustered genes (TvLF) of prokaryote origin was identified in the sequenced genome of the protozoan parasite Trichomonas vaginalis. Using a PCR based approach we confirmed the presence of the orthologous fragment in four additional T. vaginalis strains. Detailed sequence analyses unambiguously suggest that TvLF is the result of one single, recent LGT event. The proposed donor is a close relative to the firmicute bacterium Peptoniphilus harei. High nucleotide sequence similarity between T. vaginalis strains, as well as to P. harei, and the absence of homologs in other Trichomonas species, suggests that the transfer event took place after the radiation of the genus Trichomonas. Some genes have undergone pseudogenization and degradation, indicating that they may not be retained in the future. Functional annotations reveal that genes involved in informational processes are particularly prone to degradation.

Conclusions

We conclude that, although the majority of eukaryote LGTs are single gene occurrences, they may be acquired in clusters of several genes that are subsequently cleansed of evolutionarily less advantageous genes.

The scale and evolutionary significance of horizontal gene transfer in the choanoflagellate Monosiga brevicollis

Background

It is generally agreed that horizontal gene transfer (HGT) is common in phagotrophic protists. However, the overall scale of HGT and the cumulative impact of acquired genes on the evolution of these organisms remain largely unknown.

Results

Choanoflagellates are phagotrophs and the closest living relatives of animals. In this study, we performed phylogenomic analyses to investigate the scale of HGT and the evolutionary importance of horizontally acquired genes in the choanoflagellate Monosiga brevicollis. Our analyses identified 405 genes that are likely derived from algae and prokaryotes, accounting for approximately 4.4% of the Monosiga nuclear genome. Many of the horizontally acquired genes identified in Monosiga were probably acquired from food sources, rather than by endosymbiotic gene transfer (EGT) from obsolete endosymbionts or plastids. Of 193 genes identified in our analyses with functional information, 84 (43.5%) are involved in carbohydrate or amino acid metabolism, and 45 (23.3%) are transporters and/or involved in response to oxidative, osmotic, antibiotic, or heavy metal stresses. Some identified genes may also participate in biosynthesis of important metabolites such as vitamins C and K12, porphyrins and phospholipids.

Conclusions

Our results suggest that HGT is frequent in Monosiga brevicollis and might have contributed substantially to its adaptation and evolution. This finding also highlights the importance of HGT in the genome and organismal evolution of phagotrophic eukaryotes.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-729) contains supplementary material, which is available to authorized users.

Biochemical and functional characterization of novel NADH kinase in the enteric protozoan parasite Entamoeba histolytica

NAD(H) kinase catalyzes the phosphorylation of NAD(H) to form NADP(H) using ATP or inorganic polyphosphate as a phosphoryl donor. While the enzyme is conserved throughout prokaryotes and eukaryotes, remarkable differences in kinetic parameters including substrate preference, cation dependence, and physiological roles exist among the organisms. In the present study, we biochemically characterized NAD(H) kinase from the anaerobic/microaerophilic fermentative protozoan parasite Entamoeba histolytica, which lacks the conventional mitochondria capable of oxidative phosphorylation, leading to ATP. The kinetic properties of E. histolytica NAD(H) kinase recombinantly produced in Escherichia coli showed remarkable differences from those in bacteria and higher eukaryotes. Entamoeba NAD(H) kinase preferred NADH to NAD+ as the phosphoryl acceptor, utilized nucleoside triphosphates including ATP, GTP and deoxyATP, but not nucleoside di-, mono-phosphates, or inorganic polyphosphates, as the phosphoryl donor. To further understand the physiological roles in E. histolytica, we generated a stable transformant overexpressing NAD(H) kinase. Overexpression of NAD(H) kinase resulted in a 1.6-2 fold increase in the NADPH and NADP+ concentrations, a 40% reduction of the intracellular concentration of reactive oxygen species, and also led to increased tolerance toward hydrogen peroxide. These data, together with the essentially of NAD(H) kinase gene, underscore its significance as an NADP(H)-producing enzyme in this organism, and should help in designing of drugs targeting this enzyme.

Kinetic and phylogenetic analysis of plant polyamine uptake transporters

The rice gene Polyamine Uptake Transporter1 (PUT1) was originally identified based on its homology to the polyamine uptake transporters LmPOT1 and TcPAT12 in Leishmania major and Trypanosoma cruzi, respectively. Here we show that five additional transporters from rice and Arabidopsis that cluster in the same clade as PUT1 all function as high affinity spermidine uptake transporters. Yeast expression assays of these genes confirmed that uptake of spermidine was minimally affected by 166 fold or greater concentrations of amino acids. Characterized polyamine transporters from both Arabidopsis thaliana and Oryza sativa along with the two polyamine transporters from L. major and T. cruzi were aligned and used to generate a hidden Markov model. This model was used to identify significant matches to proteins in other angiosperms, bryophytes, chlorophyta, discicristates, excavates, stramenopiles and amoebozoa. No significant matches were identified in fungal or metazoan genomes. Phylogenic analysis showed that some sequences from the haptophyte, Emiliania huxleyi, as well as sequences from oomycetes and diatoms clustered closer to sequences from plant genomes than from a homologous sequence in the red algal genome Galdieria sulphuraria, consistent with the hypothesis that these polyamine transporters were acquired by horizontal transfer from green algae. Leishmania and Trypansosoma formed a separate cluster with genes from other Discicristates and two Entamoeba species. We surmise that the genes in Entamoeba species were acquired by phagotrophy of Discicristates. In summary, phylogenetic and functional analysis has identified two clades of genes that are predictive of polyamine transport activity.

Trichomonas vaginalis pathobiology new insights from the genome sequence

The draft genome of the common sexually transmitted pathogen Trichomonas vaginalis encodes one of the largest known proteome with 60,000 candidate proteins. This provides parasitologists and molecular cell biologists alike with exciting, yet challenging, opportunities to unravel the molecular features of the parasite's cellular systems and potentially the molecular basis of its pathobiology. Here, recent investigations addressing selected aspects of the parasite's molecular cell biology are discussed, including surface and secreted virulent factors, membrane trafficking, cell signalling, the degradome, and the potential role of RNA interference in the regulation of gene expression.

A novel extracellular metallopeptidase domain shared by animal host-associated mutualistic and pathogenic microbes

The mucosal microbiota is recognised as an important factor for our health, with many disease states linked to imbalances in the normal community structure. Hence, there is considerable interest in identifying the molecular basis of human-microbe interactions. In this work we investigated the capacity of microbes to thrive on mucosal surfaces, either as mutualists, commensals or pathogens, using comparative genomics to identify co-occurring molecular traits. We identified a novel domain we named M60-like/PF13402 (new Pfam entry PF13402), which was detected mainly among proteins from animal host mucosa-associated prokaryotic and eukaryotic microbes ranging from mutualists to pathogens. Lateral gene transfers between distantly related microbes explained their shared M60-like/PF13402 domain. The novel domain is characterised by a zinc-metallopeptidase-like motif and is distantly related to known viral enhancin zinc-metallopeptidases. Signal peptides and/or cell surface anchoring features were detected in most microbial M60-like/PF13402 domain-containing proteins, indicating that these proteins target an extracellular substrate. A significant subset of these putative peptidases was further characterised by the presence of associated domains belonging to carbohydrate-binding module family 5/12, 32 and 51 and other glycan-binding domains, suggesting that these novel proteases are targeted to complex glycoproteins such as mucins. An in vitro mucinase assay demonstrated degradation of mammalian mucins by a recombinant form of an M60-like/PF13402-containing protein from the gut mutualist Bacteroides thetaiotaomicron. This study reveals that M60-like domains are peptidases targeting host glycoproteins. These peptidases likely play an important role in successful colonisation of both vertebrate mucosal surfaces and the invertebrate digestive tract by both mutualistic and pathogenic microbes. Moreover, 141 entries across various peptidase families described in the MEROPS database were also identified with carbohydrate-binding modules defining a new functional context for these glycan-binding domains and providing opportunities to engineer proteases targeting specific glycoproteins for both biomedical and industrial applications.

A machine learning approach to identify hydrogenosomal proteins in Trichomonas vaginalis

The protozoan parasite Trichomonas vaginalis is the causative agent of trichomoniasis, the most widespread nonviral sexually transmitted disease in humans. It possesses hydrogenosomes-anaerobic mitochondria that generate H(2), CO(2), and acetate from pyruvate while converting ADP to ATP via substrate-level phosphorylation. T. vaginalis hydrogenosomes lack a genome and translation machinery; hence, they import all their proteins from the cytosol. To date, however, only 30 imported proteins have been shown to localize to the organelle. A total of 226 nuclear-encoded proteins inferred from the genome sequence harbor a characteristic short N-terminal presequence, reminiscent of mitochondrial targeting peptides, which is thought to mediate hydrogenosomal targeting. Recent studies suggest, however, that the presequences might be less important than previously thought. We sought to identify new hydrogenosomal proteins within the 59,672 annotated open reading frames (ORFs) of T. vaginalis, independent of the N-terminal targeting signal, using a machine learning approach. Our training set included 57 gene and protein features determined for all 30 known hydrogenosomal proteins and 576 nonhydrogenosomal proteins. Several classifiers were trained on this set to yield an import score for all proteins encoded by T. vaginalis ORFs, predicting the likelihood of hydrogenosomal localization. The machine learning results were tested through immunofluorescence assay and immunodetection in isolated cell fractions of 14 protein predictions using hemagglutinin constructs expressed under the homologous SCSα promoter in transiently transformed T. vaginalis cells. Localization of 6 of the 10 top predicted hydrogenosome-localized proteins was confirmed, and two of these were found to lack an obvious N-terminal targeting signal.

Evolutionary genomics of Entamoeba

Entamoeba histolytica is a human pathogen that causes amoebic dysentery and leads to significant morbidity and mortality worldwide. Understanding the genome and evolution of the parasite will help explain how, when and why it causes disease. Here we review current knowledge about the evolutionary genomics of Entamoeba: how differences between the genomes of different species may help explain different phenotypes, and how variation among E. histolytica parasites reveals patterns of population structure. The imminent expansion of the amount genome data will greatly improve our knowledge of the genus and of pathogenic species within it.

Horizontal gene transfers with or without cell fusions in all categories of the living matter

This article reviews the history of widespread exchanges of genetic segments initiated over 3 billion years ago, to be part of their life style, by sphero-protoplastic cells, the ancestors of archaea, prokaryota, and eukaryota. These primordial cells shared a hostile anaerobic and overheated environment and competed for survival. "Coexist with, or subdue and conquer, expropriate its most useful possessions, or symbiose with it, your competitor" remain cellular life's basic rules. This author emphasizes the role of viruses, both in mediating cell fusions, such as the formation of the first eukaryotic cell(s) from a united crenarchaeon and prokaryota, and the transfer of host cell genes integrated into viral (phages) genomes. After rising above the Darwinian threshold, rigid rules of speciation and vertical inheritance in the three domains of life were established, but horizontal gene transfers with or without cell fusions were never abolished. The author proves with extensive, yet highly selective documentation, that not only unicellular microorganisms, but the most complex multicellular entities of the highest ranks resort to, and practice, cell fusions, and donate and accept horizontally (laterally) transferred genes. Cell fusions and horizontally exchanged genetic materials remain the fundamental attributes and inherent characteristics of the living matter, whether occurring accidentally or sought after intentionally. These events occur to cells stagnating for some 3 milliard years at a lower yet amazingly sophisticated level of evolution, and to cells achieving the highest degree of differentiation, and thus functioning in dependence on the support of a most advanced multicellular host, like those of the human brain. No living cell is completely exempt from gene drains or gene insertions.

Trichomonas transmembrane cyclases result from massive gene duplication and concomitant development of pseudogenes

Background

Trichomonas vaginalis has an unusually large genome (∼160 Mb) encoding ∼60,000 proteins. With the goal of beginning to understand why some Trichomonas genes are present in so many copies, we characterized here a family of ∼123 Trichomonas genes that encode transmembrane adenylyl cyclases (TMACs).

Methodology/Principal Findings

The large family of TMACs genes is the result of recent duplications of a small set of ancestral genes that appear to be unique to trichomonads. Duplicated TMAC genes are not closely associated with repetitive elements, and duplications of flanking sequences are rare. However, there is evidence for TMAC gene replacements by homologous recombination. A high percentage of TMAC genes (∼46%) are pseudogenes, as they contain stop codons and/or frame shifts, or the genes are truncated. Numerous stop codons present in the genome project G3 strain are not present in orthologous genes of two other Trichomonas strains (S1 and B7RC2). Each TMAC is composed of a series of N-terminal transmembrane helices and a single C-terminal cyclase domain that has adenylyl cyclase activity. Multiple TMAC genes are transcribed by Trichomonas cloned by limiting dilution.

Conclusions/Significance

We conclude that one reason for the unusually large genome of Trichomonas is the presence of unstable families of genes such as those encoding TMACs that are undergoing massive gene duplication and concomitant development of pseudogenes.

Trichomonas vaginalis is the only medically important protist (single-cell eukaryote) that is sexually transmitted. The ∼160-Mb Trichomonas genome contains more predicted protein-encoding genes (∼60,000) than the human genome. To begin to understand why there are so many copies of some genes, we chose here to study a large family of genes encoding unique transmembrane cyclases. Our most important results include the following. More than 100 transmembrane cyclase genes do not result from chromosomal duplications, because for the most part only the coding regions of the genes, rather than flanking sequences, are duplicated. Almost half of the transmembrane cyclase genes are pseudogenes, and these pseudogenes are polymorphic among laboratory strains of Trichomonas. Messenger RNAs for numerous transmembrane cyclases are expressed simultaneously, and representative cyclase domains have adenylyl cyclase activity. In summary, the large family of Trichomonas genes encoding transmembrane adenylyl cyclases results from massive gene duplication and concomitant development of pseudogenes.

Trichomonas vaginalis vast BspA-like gene family: evidence for functional diversity from structural organisation and transcriptomics

Background

Trichomonas vaginalis is the most common non-viral human sexually transmitted pathogen and importantly, contributes to facilitating the spread of HIV. Yet very little is known about its surface and secreted proteins mediating interactions with, and permitting the invasion and colonisation of, the host mucosa. Initial annotations of T. vaginalis genome identified a plethora of candidate extracellular proteins.

Results

Data mining of the T. vaginalis genome identified 911 BspA-like entries (TvBspA) sharing TpLRR-like leucine-rich repeats, which represent the largest gene family encoding potential extracellular proteins for the pathogen. A broad range of microorganisms encoding BspA-like proteins was identified and these are mainly known to live on mucosal surfaces, among these T. vaginalis is endowed with the largest gene family. Over 190 TvBspA proteins with inferred transmembrane domains were characterised by a considerable structural diversity between their TpLRR and other types of repetitive sequences and two subfamilies possessed distinct classic sorting signal motifs for endocytosis. One TvBspA subfamily also shared a glycine-rich protein domain with proteins from Clostridium difficile pathogenic strains and C. difficile phages. Consistent with the hypothesis that TvBspA protein structural diversity implies diverse roles, we demonstrated for several TvBspA genes differential expression at the transcript level in different growth conditions. Identified variants of repetitive segments between several TvBspA paralogues and orthologues from two clinical isolates were also consistent with TpLRR and other repetitive sequences to be functionally important. For one TvBspA protein cell surface expression and antibody responses by both female and male T. vaginalis infected patients were also demonstrated.

Conclusions

The biased mucosal habitat for microbial species encoding BspA-like proteins, the characterisation of a vast structural diversity for the TvBspA proteins, differential expression of a subset of TvBspA genes and the cellular localisation and immunological data for one TvBspA; all point to the importance of the TvBspA proteins to various aspects of T. vaginalis pathobiology at the host-pathogen interface.

Mobile DNA and evolution in the 21st century

Scientific history has had a profound effect on the theories of evolution. At the beginning of the 21st century, molecular cell biology has revealed a dense structure of information-processing networks that use the genome as an interactive read-write (RW) memory system rather than an organism blueprint. Genome sequencing has documented the importance of mobile DNA activities and major genome restructuring events at key junctures in evolution: exon shuffling, changes in cis-regulatory sites, horizontal transfer, cell fusions and whole genome doublings (WGDs). The natural genetic engineering functions that mediate genome restructuring are activated by multiple stimuli, in particular by events similar to those found in the DNA record: microbial infection and interspecific hybridization leading to the formation of allotetraploids. These molecular genetic discoveries, plus a consideration of how mobile DNA rearrangements increase the efficiency of generating functional genomic novelties, make it possible to formulate a 21st century view of interactive evolutionary processes. This view integrates contemporary knowledge of the molecular basis of genetic change, major genome events in evolution, and stimuli that activate DNA restructuring with classical cytogenetic understanding about the role of hybridization in species diversification.

Integrating horizontal gene transfer and common descent to depict evolution and contrast it with "common design"

Horizontal gene transfer (HGT) and common descent interact in space and time. Because events of HGT co-occur with phylogenetic evolution, it is difficult to depict evolutionary patterns graphically. Tree-like representations of life's diversification are useful, but they ignore the significance of HGT in evolutionary history, particularly of unicellular organisms, ancestors of multicellular life. Here we integrate the reticulated-tree model, ring of life, symbiogenesis whole-organism model, and eliminative pattern pluralism to represent evolution. Using Entamoeba histolytica alcohol dehydrogenase 2 (EhADH2), a bifunctional enzyme in the glycolytic pathway of amoeba, we illustrate how EhADH2 could be the product of both horizontally acquired features from ancestral prokaryotes (i.e. aldehyde dehydrogenase [ALDH] and alcohol dehydrogenase [ADH]), and subsequent functional integration of these enzymes into EhADH2, which is now inherited by amoeba via common descent. Natural selection has driven the evolution of EhADH2 active sites, which require specific amino acids (cysteine 252 in the ALDH domain; histidine 754 in the ADH domain), iron- and NAD(+) as cofactors, and the substrates acetyl-CoA for ALDH and acetaldehyde for ADH. Alternative views invoking "common design" (i.e. the non-naturalistic emergence of major taxa independent from ancestry) to explain the interaction between horizontal and vertical evolution are unfounded.