1471-21482013131190http://www.biomedcentral.com/1471-2148/13/1902401577810.1186/1471-2148-13-19012620136920139920132013Alves et al.; licensee BioMed Central Ltd.This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.EndosymbiosisTrypanosomatidsAmino acid biosynthesisPhylogenyGenomic analysesMetabolic pathway evolutionProteobacteria

Abstract

Background

Trypanosomatids of the genera Angomonas and Strigomonas live in a mutualistic association characterized by extensive metabolic cooperation with obligate endosymbiotic Betaproteobacteria. However, the role played by the symbiont has been more guessed by indirect means than evidenced. Symbiont-harboring trypanosomatids, in contrast to their counterparts lacking symbionts, exhibit lower nutritional requirements and are autotrophic for essential amino acids. To evidence the symbiont’s contributions to this autotrophy, entire genomes of symbionts and trypanosomatids with and without symbionts were sequenced here.

Results

Analyses of the essential amino acid pathways revealed that most biosynthetic routes are in the symbiont genome. By contrast, the host trypanosomatid genome contains fewer genes, about half of which originated from different bacterial groups, perhaps only one of which (ornithine cyclodeaminase, EC:4.3.1.12) derived from the symbiont. Nutritional, enzymatic, and genomic data were jointly analyzed to construct an integrated view of essential amino acid metabolism in symbiont-harboring trypanosomatids. This comprehensive analysis showed perfect concordance among all these data, and revealed that the symbiont contains genes for enzymes that complete essential biosynthetic routes for the host amino acid production, thus explaining the low requirement for these elements in symbiont-harboring trypanosomatids. Phylogenetic analyses show that the cooperation between symbionts and their hosts is complemented by multiple horizontal gene transfers, from bacterial lineages to trypanosomatids, that occurred several times in the course of their evolution. Transfers occur preferentially in parts of the pathways that are missing from other eukaryotes.

Conclusion

We have herein uncovered the genetic and evolutionary bases of essential amino acid biosynthesis in several trypanosomatids with and without endosymbionts, explaining and complementing decades of experimental results. We uncovered the remarkable plasticity in essential amino acid biosynthesis pathway evolution in these protozoans, demonstrating heavy influence of horizontal gene transfer events, from Bacteria to trypanosomatid nuclei, in the evolution of these pathways.

Background

Many protozoan and metazoan cells harbor vertically inherited endosymbionts in their cytoplasm. Prominent among them are the associations between Alphaproteobacteria and leguminous root cells, as well as Gammaproteobacteria and cells lining the digestive tube of insects. Comprehensive reviews have covered most aspects of these ancient mutualistic relationships, including metabolism, genetics, and evolutionary history of the endosymbiont/host cell associations 1234567. Much less is known about the relationship between protists and their bacterial endosymbionts 8910, including the symbiosis between trypanosomatids and Betaproteobacteria, herein examined 11121314.

The Trypanosomatidae (Euglenozoa, Kinetoplastea) are well studied mainly because species of the genera Trypanosoma and Leishmania are pathogenic in humans and domestic animals 15. However, despite their importance, these pathogens are a minority within the family, and most species are non-pathogenic commensals in the digestive tube of insects 161718. Usually, trypanosomatids are nutritionally fastidious and require very rich and complex culture media, however a very small group of these protozoa can be cultivated in very simple and defined media 1920212223. This reduced group of insect trypanosomatids carries cytoplasmic endosymbionts and is known as symbiont-harboring trypanosomatids, to distinguish them from regular insect trypanosomatids naturally lacking symbionts. Symbiont-harboring trypanosomatids belong to the genera Strigomonas and Angomonas24, and their lesser nutritional requirements indicate that they have enhanced biosynthetic capabilities. In a few cases, it has been shown that the symbiotic bacterium contains enzymes involved in host biosynthetic pathways, but in most cases the metabolic contribution of the endosymbiont has been inferred from nutritional data rather than genetically demonstrated 1214.

Each symbiont-harboring trypanosomatid carries just one symbiont in its cytoplasm, which divides synchronously with other host cell structures and is vertically transmitted 14. The endosymbionts’ original association with an ancestral trypanosomatid is thought to have occurred sometime in the Cretaceous 13. This long partnership has led to considerable changes on the endosymbiont genomes including gene loss, with clear preferential retention of genes involved in metabolic collaboration with the host, and consequent genomic size reduction 2526, as seen in other obligatory symbiotic associations 127272829.

Extensive comparative studies between symbiont-harboring trypanosomatids (wild and cured strains, obtained after antibiotic treatment) and regular trypanosomatids have permitted inferences about the symbiont dependence and contribution in the overall metabolism, in particular phospholipid 303132 and amino acid 333435363738394041 production of the host cell. Previous comparative studies on these organisms, often involving trace experiments using radioactive compounds, reported the requirement, substitution, and sparing of amino acids in culture media 22334243444546474849505152535455. Nutritional data revealed that, as for most animals, including humans, the amino acids lysine, histidine, threonine, isoleucine, leucine, methionine, cysteine, tryptophan, valine, phenylalanine, tyrosine, and arginine/citrulline are essential for regular trypanosomatids. However, similar analyses showed that symbiont-harboring trypanosomatids require only methionine or tyrosine in culture media, suggesting that they possess the necessary enzymatic equipment to synthesize most amino acids 20222334. Unfortunately, besides the symbiont-harboring trypanosomatids, most of these studies were performed only on Crithidia fasciculata, largely ignoring other trypanosomatids. Of the hundreds of enzymes known to be involved in the synthesis of essential amino acids in other organisms, only a few, i.e., diaminopimelic decarboxylase, threonine deaminase, ornithine carbamoyl transferase, argininosuccinate lyase, citrulline hydrolase, ornithine acetyl transferase, acetyl ornithinase, and arginase have been identified and characterized in trypanosomatids 2133353637383940414656. Thus, in contrast to the advanced state of knowledge of genes involved in amino acid biosynthesis in many microorganisms 57, the potential for amino acid synthesis in trypanosomatids remains largely unknown. In symbiont-harboring trypanosomatids, nutritional inferences provided little information about the effective participation of the symbiotic bacterium in the various metabolic pathways of the host protozoan. This contrasts with the advancement of knowledge about the presence/absence of genes for complete pathways for amino acid synthesis in many microorganisms.

Herein, we have identified the genes involved in the biosynthetic pathways of the essential amino acids in the genomes of symbiont-harboring and regular trypanosomatids of different genera (see Methods), through the characterization of each gene by similarity searches and protein domain analyses. We apply extensive phylogenetic inferences to determine the most likely origins of these genes, as it has been previously shown that other important metabolic enzymes in trypanosomatids have been transferred from bacteria, other than the present symbiont 58. Although detection of a gene with a presumed function does not definitely prove its activity, the association of its presence with complementary nutritional and biochemical data supports the conclusion that it functions as predicted. In the present work, we establish the clear and defined contribution of endosymbionts to the amino acid metabolism of their trypanosomatid hosts, which is related to high amounts of lateral transfer of genes from diverse bacterial lineages to trypanosomatid genomes.

Results and discussion

In this work, the presence or absence of a given gene for a particular enzyme was verified in the genomes of endosymbionts, symbiont-harboring and regular trypanosomatids and then compared with the available nutritional and enzymatic data on essential amino acid biosynthesis in insect trypanosomatids. Extensive phylogenetic analyses were also performed on most of the identified trypanosomatid genes, in addition to some symbiont genes of interest. Data are mostly limited to the regular and symbiont-harboring trypanosomatid and endosymbiont genomes that have been sequenced here. Although the genomes of all available symbiont-harboring trypanosomatids and endosymbionts have been examined, only a very limited sample of regular trypanosomatid genomes (H. muscarum and C. acanthocephali) was included in these analyses, precluding generalizations about trypanosomatids as a whole. Data on the genomes of leishmaniae and trypanosomes available in KEGG were also used for comparison, but a wider sampling of genomes from more diverse groups of Trypanosomatidae and other, more distant Kinetoplastida will be necessary to enable more generalizing conclusions on the evolution of essential amino acid synthesis pathways in these organisms.

Given the incomplete nature of the trypanosomatid genomes sequenced here and the possibility of contaminant sequences, we have taken extensive precautions before including each gene in our analyses (see Methods). Our genomic context analyses of the genes identified as horizontally transferred show (Additional file 1) that genes used in this work occurred, with one exception, in long contigs presenting the typical trypanosomatid architecture of long stretches of genes in the same orientation. Moreover, all these genes overwhelmingly matched those from previously sequenced trypanosomatids. The one exception is a gene (2.7.1.100, see below) that occurs only in the two regular trypanosomatids sequenced here, and whose sequences are isolated in short contigs. As described below, they form a monophyletic group in the phylogeny. GC percent (Additional file 1) and sequencing coverage (Additional file 2) analyses also show that all genes identified in this work present statistics typical of other genes from these organisms. In short, these data show that the trypanosomatid genes employed here are highly unlikely to be contaminants.

Additional file 1

Genomic context and GC content for candidate HGT genes in the Trypanosomatidae analyzed in this work. Arrows show TBLASTN alignments of the genome against UniRef100 and KEGG proteins. Alignment orientation is displayed in blue or red, except for the alignment for the gene currently in focus, which is colored black. Coordinates are in kilobases.

Click here for file

Additional file 2

Summary of phylogenetic and genome coverage analyses of the candidate HGT genes in the Trypanosomatidae analyzed in this work and a few other genes of interest.

Click here for file

Pathways of amino acid synthesis

Lysine

Lysine, as well as methionine and threonine, are essential amino acids generated from aspartate, a non-essential amino acid, which is synthesized from oxaloacetate that is produced in the Krebs cycle. There are two main routes for the biosynthesis of lysine: the diaminopimelate (DAP) and the aminoadipate (AA) pathways. The former is largely confined to bacteria, algae, some fungi, and plants, whereas the latter is described in fungi and euglenids 5960616263.

Early nutritional studies 46 showed that lysine is essential for the growth of regular trypanosomatids, but could be efficiently replaced by DAP. In accordance, radioactive tracer and enzymatic experiments revealed that DAP is readily incorporated as lysine into proteins. Moreover, DAP-decarboxylase (EC:4.1.1.20), the enzyme that converts DAP into lysine, was detected in cell homogenates of C. fasciculata46. Nevertheless, either lysine or DAP were always necessary for growth of these flagellates in defined medium, indicating that the lysine pathway was somehow incomplete. In contrast, symbiont-harboring trypanosomatids required neither lysine nor DAP to grow in defined media 1920212223. Interestingly, the genes encoding the nine enzymes of the bacterial-type DAP pathway, leading from aspartate to lysine, were identified in the genomes of all endosymbionts (Figure 1). In contrast, only the final gene of the DAP pathway was found in the genomes of the symbiont-harboring trypanosomatids, and the final two found in one regular trypanosomatid examined (H. muscarum), which explains why DAP could substitute for lysine in growth media of some regular trypanosomatids. There are no genes for lysine biosynthesis annotated in the leishmaniae and trypanosomes present in KEGG. It is worth mentioning that, with respect to the alternative AA pathway, we were unable to find any genes for the synthesis of lysine in any of the endosymbiont, symbiont-harboring or regular trypanosomatid genomes analyzed.

In summary, our findings using comparative genomics are in concordance with the data from previous nutritional and enzymatic studies, showing that only symbiont-harboring trypanosomatids, and not regular ones, are autotrophic for lysine and that this autonomy is provided by the DAP pathway present in their symbionts. The presence of DAP-decarboxylase in symbiont-harboring trypanosomatids may suggest that although the symbiont contains the great majority of genes for the lysine production, the host protozoan somehow controls the production of this essential amino acid.

Methionine and cysteine

Methionine is included in all defined media designed for the growth of trypanosomatids with or without symbionts 202243, suggesting that these protozoans are incapable of methionine synthesis. However, experimental evidence has shown that homocysteine and/or cystathionine could substitute for methionine in culture media for trypanosomatids 434564.

Our analyses suggest that regular trypanosomatids and symbiont-harboring trypanosomatids have the necessary genes to produce cystathionine, homocysteine, and methionine from homoserine (Figure 2), whereas the endosymbiont genomes have no gene for the enzymes involved in the synthesis of methionine from homoserine. However, homoserine is produced from aspartate semialdehyde through the mediation of homocysteine methyltransferase (EC:1.1.1.3), which is universally present in the genomes of all the endosymbionts, symbiont-harboring and regular trypanosomatids examined.

With respect to cysteine synthesis, it has been shown that the incubation of cell homogenates of C. fasciculata with 35S-methionine produced radioactive adenosyl-methionine (SAM), adenosyl-homocysteine (SAH), homocysteine, cystathionine, and cysteine 45. Thus, this trypanosomatid is fully equipped to methylate methionine to produce homocysteine and, thereon, to convert homocysteine into cysteine through the trans-sulfuration pathway. However, with respect to the cystathionine/cysteine interconversion, there is some ambiguity concerning the presence or absence of cystathionine gamma-lyase (EC:4.4.1.1) in regular trypanosomatids. Many sulfhydrolases have a domain composition very similar to that of EC:4.4.1.1, which makes a definitive in silico function assignment to any of them difficult. Specifically, the enzymes cystathionine gamma-synthase (EC:2.5.1.48) and O-acetylhomoserine aminocarboxypropyltransferase (EC:2.5.1.49), and the two versions of cystathionine beta-lyase (EC:4.4.1.8) are possible candidates to mediate the trans-sulfuration step attributed to EC:4.4.1.1, but further research is required to establish which of these enzymes, if any, performs that reaction. We also found that, in addition to the standard pathway for methionine/cysteine synthesis (Figure 2, compounds III-X), all symbiont-harboring and regular trypanosomatids examined had the genes to produce cysteine from serine in a simple two-step reaction, with acetylserine as an intermediate (Figure 2, I-III).

In summary, if regular and symbiont-harboring trypanosomatids are capable of interconverting methionine and cysteine, as shown for C. fasciculata43, none of these two amino acids can be considered essential for trypanosomatids as the presence of one renders the other unnecessary. In that case, both can be synthesized by trypanosomatids, without any participation of their symbionts, except in the optional production of aspartate semialdehyde and homoserine. However, the expression of these genes remains to be confirmed.

Threonine

In trypanosomatids, initial investigations about the nutritional requirements for threonine were controversial. Most results suggested that this amino acid is essential 434548646566, but other studies considered the addition of threonine to the growth media of regular trypanosomatids unnecessary 33. Our genomic analysis favors the latter observations.

Threonine, one of the precursors of isoleucine, can be produced by different biosynthetic pathways. We have examined two of these possible routes, one starting from glycine and the other from aspartate, as presented in Figure 3. The conversion of glycine plus acetoaldehyde into threonine is mediated by threonine aldolase (EC:4.1.2.5). The gene for this enzyme is absent from endosymbionts but present in the genomes of symbiont-harboring trypanosomatids and C. acanthocephali, but not Herpetomonas. It is also absent from the genomes of trypanosomes but present in the genome of Leishmania major (KEGG data).

The pathway from aspartate utilizes the first two enzymes (EC:2.7.2.4 and EC:1.2.1.11) of the DAP pathway from lysine synthesis for the production of aspartate semialdehyde. These genes are present exclusively in the symbiont genomes. Aspartate semialdehyde is then sequentially converted into homoserine, phosphohomoserine, and threonine. The gene encoding homoserine dehydrogenase (EC:1.1.1.3) is universally present in the genomes of the endosymbionts, symbiont-harboring and regular trypanosomatids. It is also present in the genomes of T. cruzi and Leishmania spp. In contrast, the genes for the enzymes leading from homoserine to threonine via phosphohomoserine (EC:2.7.1.39 and EC:4.2.3.1) are present in the genomes of all insect trypanosomatids (including symbiont-harboring ones), of Trypanosoma spp., and Leishmania spp., but totally absent from the endosymbiont genomes.

Thus, the genetic constitution of regular trypanosomatids is consistent with earlier nutritional data showing the insect trypanosomatids, with or without symbionts, to be autotrophic for threonine. This observation suggests that endosymbionts are able to enhance the host cell threonine synthesis by producing the metabolic precursor aspartate semialdehyde that is also involved in other metabolic pathways.

The overall genomic and enzymatic picture is in apparent contradiction with early nutritional findings showing that threonine promoted the growth of trypanosomatids in culture 67. This contradiction might find its basis in the fact that endogenously produced threonine is required by many metabolic processes, such that supplementation of the culture media could enhance the growth of the trypanosomatids.

Isoleucine, valine, and leucine

Isoleucine, valine, and leucine are considered essential nutrients for the growth of all trypanosomatids, except symbiont-harboring ones. The canonic pathway for the synthesis of isoleucine is depicted in Figure 4. Oxobutanoate (alpha-ketoglutaric acid) is the starting point of the pathway, and can be produced in two ways: from threonine (Figure 4, compounds II-III) or from pyruvate (Figure 4, compounds I, IX). The conversion of threonine into oxobutanoate is mediated by threonine deaminase (EC:4.3.1.19). The specific activity of this enzyme was higher in symbiont-enriched subcellular fractions of symbiont-harboring trypanosomatid homogenates than in any other cell fraction or in the cytosol, suggesting that this enzyme was located in the symbiont 33. However, genes for EC:4.3.1.19 are present in the genomes of endosymbionts, as well as those of symbiont-harboring and regular trypanosomatids (except Leishmania and Trypanosoma), contrasting with enzymatic determinations showing the absence of enzyme activity in regular trypanosomatids 33. Since the presence of the gene does not guarantee the functionality of the enzyme for that specific reaction, the issue remains to be experimentally verified. The next enzymatic step, the transference of the acetaldehyde from pyruvate to oxobutanoate, is mediated by the enzyme acetolactate synthase (EC:2.2.1.6), which is present exclusively in the genomes of endosymbionts. Also present only in symbionts are the genes for the next four enzymes of the pathway, which are common for valine and isoleucine synthesis. However, the gene for a branched-chain amino acid transaminase (EC:2.6.1.42), mediating the last step in the synthesis of isoleucine, valine, and leucine, is present in the genomes of symbiont-harboring and regular trypanosomatids, but not endosymbionts.

The first step of the valine pathway is the conversion of pyruvate into hydroxymethyl ThPP, mediated by an enzyme of the pyruvate dehydrogenase complex (EC:1.2.4.1) whose gene is present in the genomes of endosymbionts and symbiont-harboring and regular trypanosomatids. The next reaction, leading to acetolactate, is mediated by acetolactate synthase (EC:2.2.1.6), whose gene is present exclusively in the genomes of the endosymbionts. The reactions that follow from acetoacetate into valine involve the same endosymbiont genes from isoleucine synthesis.

Synthesis of leucine uses oxoisovalerate, an intermediate metabolite of the valine pathway that is converted into isopropylmalate by 2-isopropylmalate synthase (EC:2.3.3.13), encoded by a gene present only in the endosymbionts – as are the genes for the enzymes catalyzing the next three steps for leucine biosynthesis. The presence of the gene for this branched-chain amino acid transaminase (EC:2.6.1.42) in the genomes of regular trypanosomatids explains the earlier finding that oxopentanoate and oxoisovalerate, the immediate precursors of isoleucine, valine, and leucine could substitute for these amino acids when added to regular trypanosomatid synthetic culture media 43. Interestingly, this gene is present in all symbiont-harboring and regular trypanosomatid genomes examined, but absent from endosymbiont genomes (Figure 4). It is also present in the genomes of T. brucei and the leishmaniae available from KEGG. In addition to isoleucine, valine, and leucine biosynthesis, this enzyme also participates in the degradation of these amino acids for their use in other metabolic processes in the cell, which might explain the presence of this enzyme as the only representative of the pathway in all regular trypanosomatids examined.

A coupled biosynthetic pathway of the branched-chain amino acids was also described for the symbiotic bacterium Buchnera and its aphid host, where the symbiont has the capability to synthesize the carbon skeleton of these amino acids but lacks the genes for the terminal transaminase reactions 6869. The aphid possesses genes hypothesized to accomplish these missing steps, even if orthologs of those are found in other insects and carry out different functions 70. The branched-chain amino acid transaminase (EC:2.6.1.42) encoded by an aphid gene was shown to be up-regulated in the bacteriocytes, supporting the cooperation of Buchnera and its host in the synthesis of essential amino acids 71. Since this transamination involves the incorporation of amino-N and the aphid diet is low in nitrogen, the host mediation of this step would be a way of maintaining a balanced profile of amino acids through transamination between those that are over abundant and those that are rare 7172.

In summary, the presence in endosymbionts of most genes involved in isoleucine, valine and leucine synthesis explains why symbiont-harboring trypanosomatids, but not regular ones, are autotroph for these essential amino acids. However, it is worth noting that the presence of the branched-chain amino acid transaminase in trypanosomatids indicates that the host might control amino acid production according to their necessity and the nutrient availability in the medium.

Phenylalanine, tyrosine, and tryptophan

There are no enzymatic data concerning the synthesis of phenylalanine, tryptophan, and tyrosine in trypanosomatids. However, it is well known that these amino acids are essential in defined culture media designed for regular trypanosomatids, but not for symbiont-harboring ones 20224344. The biosynthetic routes for these three amino acids use chorismate, which is produced from phosphoenolpyruvate (PEP) via the shikimate pathway, as a common substrate. The genomes of all endosymbionts contain the genes for this route, while the genomes of symbiont-harboring and regular trypanosomatids do not (Figure 5).

The genes for the enzymes converting chorismate into prephenate and for transforming this compound into phenylalanine and tyrosine are present in all endosymbiont genomes. Symbiont-harboring and regular trypanosomatid genomes also have the genes for the last step in the synthesis of phenylalanine and tyrosine, but it is not known whether all of these enzymes are functional. The gene for phenylalanine-4-hydroxylase (EC:1.14.16.1), which converts phenylalanine into tyrosine, is present in symbiont-harboring and regular trypanosomatids, including the leishmaniae, but not in endosymbionts. Similarly, this enzyme is present only in the aphid concerning the metabolic partnership between Buchnera and its insect host. Furthermore, the gene encoding this enzyme is up-regulated in bacteriocytes, thus enhancing the production and interconversion of such amino acids 71. On the other hand, endosymbionts have an additional route for the synthesis of phenylalanine from prephenate, involving enzymes aromatic-amino-acid aminotransferase (EC:2.6.1.57) and prephenate dehydratase (EC:4.2.1.51), whose genes are absent in symbiont-harboring and regular trypanosomatid genomes.

The case of the last enzyme of the tryptophan pathway is rather interesting. Tryptophan synthase (EC:4.2.1.20) possesses two subunits. This bi-enzyme complex (a tetramer of two alpha and two beta subunits) channels the product of the alpha subunit (indole) to the beta subunit, which condenses indole and serine into tryptophan 73. Both subunits are present in the endosymbionts, whereas the genomes of symbiont-harboring trypanosomatids and H. muscarum have only the beta subunit. None of the other trypanosomatid genomes examined presented either subunit of tryptophan synthase.

In summary, the endosymbionts have all the genes for the different routes leading from chorismate to tryptophan, tyrosine, and phenylalanine, which are absent from symbiont-harboring and regular trypanosomatid genomes. This obviously prevents regular trypanosomatids from synthesizing any of these three amino acids and growing without supplementation. It is worth observing that the presence of phenylalanine hydroxylase, which converts phenylalanine into tyrosine, in trypanosomatids but not in endosymbionts indicates that the host might control tyrosine production.

Histidine

Histidine is derived from three precursors: the ATP purine ring furnishes a nitrogen and a carbon, the glutamine contributes with the second ring nitrogen, while PRPP donates five carbons. Histidine is a truly essential amino acid for most trypanosomatids, as corroborated by its obligatory presence in every synthetic media so far devised for regular trypanosomatid growth 224344. Accordingly, symbiont-harboring and regular trypanosomatid genomes do not seem to carry a single gene for histidine synthesis (Figure 6). All genes for the enzymes that participate in its biosynthesis, except the gene for histidinol-phosphate phosphatase (HPP, EC:3.1.3.15), which converts histidinol phosphate into histidinol, are present in the endosymbiont genomes. Since symbiont-harboring trypanosomatids do not require histidine, it is presumed that the absent EC:3.1.3.15 is replaced by an equivalent enzyme yet to be characterized (see Other observation on amino acid pathway peculiarities).

Organisms autotrophic for ornithine use the glutamate pathway 74 for its synthesis via acetylated compounds as represented in Figure 7 (I-VI). All genes for this pathway are present in the genomes of endosymbionts. The last step in the synthesis of ornithine can also be performed by the enzymes aminoacylase (EC:3.5.1.14) or acetylornithine deacetylase (EC:3.5.1.16), which convert acetylornithine into ornithine and are present in the genomes of symbiont-harboring and regular trypanosomatids, but not endosymbionts.

As represented in Figure 7, organisms lacking the glutamate pathway for the synthesis of ornithine can nevertheless produce it by different routes utilizing either citrulline or arginine 373954. Ornithine can be produced from the hydrolysis of citrulline mediated by citrulline hydrolase (EC:3.5.1.20). This activity is present in cell homogenates of all trypanosomatids, except the leishmaniae and trypanosomes, but the corresponding gene has not yet been identified to date in any organism, making it impossible to perform similarity searches. Ornithine can also be produced from arginine by means of arginase (EC:3.5.3.1), which splits arginine into ornithine and urea. The gene for arginase is present in the genomes of symbiont-harboring trypanosomatids and some regular trypanosomatids (Leishmania and C. acanthocephali), but not in the genomes of endosymbionts or H. muscarum – although a fragment was found in the later (see HGT and arginine and ornithine biosynthesis).

Arginine can be synthesized from ornithine through a recognized universal enzymatic pathway 74, the first step of which is the conversion of ornithine and carbamoyl phosphate into citrulline mediated by OCT (ornithine carbamoyl transferase, EC:2.1.3.3). The gene for OCT was found in the genomes of all endosymbionts and also in Herpetomonas, but was absent from other regular, as well as symbiont-harboring, trypanosomatid genomes examined. These findings confirm earlier immunocytochemical ultrastructural experiments showing the presence of OCT in the symbiont of Angomonas deanei36. The absence of the OCT gene renders most trypanosomatids unable to make citrulline from ornithine 75. However, the genes for the remaining enzymes leading from citrulline into arginine are all present in the genomes of all regular and symbiont-harboring trypanosomatids, but absent from the endosymbiont genomes. These data are in full accordance with earlier enzymatic determinations for argininosuccinate synthase (EC:6.3.4.5), argininosuccinate lyase (EC:4.3.2.1), and arginase (EC:3.5.3.1) in cell homogenates of trypanosomatids 383956.

Taking all these data together, we can conclude that regular trypanosomatids require exogenous sources of arginine or citrulline in their culture medium to produce ornithine. This is related to the fact that regular trypanosomatids lack the glutamate pathway for ornithine synthesis. Furthermore, ornithine cannot substitute for arginine or citrulline because most regular trypanosomatids lack OCT. Conversely, symbiont-harboring trypanosomatids are autotrophic for ornithine. This is due to the fact that, although the symbiont lacks most genes for ornithine production, it contains sequences for key enzymes such as those for the glutamate route and OCT, which converts ornithine into citrulline thus completing the urea cycle.

Polyamines

As shown in Figure 7, putrescine, a polyamine associated with cell proliferation, can be produced from ornithine in a one-step reaction mediated by ODC (ornithine decarboxylase, EC:4.1.1.17), whose gene is present in the genomes from the genus Angomonas and regular trypanosomatids, but not in endosymbionts or Strigomonas. Interestingly, it was proposed that the symbiont can enhance the ODC activity of A. deanei by producing protein factors that increase the production of polyamines in the host trypanosomatid 76. Such high ODC activity may be directly connected to the lowest generation time described for trypanosomatids that is equivalent to 6 hours 13. Putrescine could also be produced from agmatine since the genomes of regular and symbiont-harboring trypanosomatids have the gene for agmatinase (EC:3.5.3.11), converting agmatine into putrescine. However, the gene for the enzyme arginine decarboxylase (EC:4.1.1.19), which synthesizes agmatine, is present solely in the genomes of endosymbionts, thus completing the biosynthetic route for this polyamine, via agmatinase, in symbiont-harboring trypanosomatids. Putrescine is then converted to spermidine and spermine by enzymes S-adenosylmethionine decarboxylase (EC:4.1.1.50) and spermidine synthase (EC:2.5.1.16). The genes for these enzymes are present in the regular and symbiont-harboring trypanosomatids, but not in endosymbionts (Figure 7). Enzyme EC:2.5.1.16, converting S-adenosylmethioninamine and putrescine into S-methyl-5’-thioadenosine and spermidine, also participates in a reaction from the methionine salvage pathway. This pathway is present, complete in all symbiont-harboring and regular trypanosomatids examined (Additional file 3), although there are questions regarding the step catalyzed by acireductone synthase (EC:3.1.3.77, see HGT and methionine and cysteine biosynthesis).

Our data on the phylogeny of the genes for essential amino acid biosynthesis have clearly shown that the genes present in the symbionts are of betaproteobacterial origin (for an illustrative example, see Figure 8), as shown before for the genes of heme synthesis 58 and many others across the endosymbiont genomes 25. The symbiont-harboring and regular trypanosomatid genomes, on the other hand, present a rather different situation. Thus, 18 of the 39 genes required for the biosynthesis of essential amino acids exhibited at least some phylogenetic evidence of having been horizontally transferred from a bacterial group to a trypanosomatid group, with three other genes presenting undetermined affiliation (see Additional file 2 for a summary of the phylogenetic analyses results). As detailed below, horizontal gene transfer (HGT) events seem to have originated from a few different bacterial taxa, although in some cases the exact relationship was not completely clear. Also, while some transfers are common to all trypanosomatid groups examined, others were found to be specific to certain subgroups. This could be due to multiple HGT events from associated bacteria at different points of the family’s evolutionary history or, alternatively, to HGT events that occurred in the common ancestor of all trypanosomatids, whose corresponding genes were later differentially lost in certain taxa. Given the low number of genomes currently known in the family, it is difficult to assign greater probability to either scenario.

Maximum likelihood phylogenetic tree of homoserine dehydrogenase (EC:1.1.1.3). A – overall tree, colored according to taxonomic affiliation of each taxon, as per the legend on the right; distance bar only applies to panel A. B – details of the region of the tree where the Ca. Kinetoplastibacterium spp. are placed. C – details of the region of the tree where the Trypanosomatidae are placed. Values on nodes represent bootstrap support (only 50 or greater shown). Panels B and C are meant to only represent the branching patterns and do not portray estimated distances between sequences. Abbreviations as for Figure 1.

Regarding the taxonomic affiliation of the putative origin of these HGT events, it is possible to notice a preponderance of bacteria from a few phyla with three or more genes transferred, i.e. Firmicutes, Bacteroidetes, and Gammaproteobacteria, plus a few other phyla with two or less genes represented, like Actinobacteria, Betaproteobacteria, Acidobacteria, and Alphaproteobacteria. In a few other cases, the trypanosomatid genes grouped inside diverse bacterial phyla, in which case the assignment of a definite originating phylum was not possible. However, given the sometimes high rate of HGT in prokaryotic groups, it is difficult to assess with confidence the correct number of putative HGT events from Bacteria to Trypanosomatidae. It is possible that some of the genes that seem to have originated from different phyla could actually have come from one bacterial line that was itself the recipient of one or more previous HGT events from other bacteria.

Analysis of all generated phylogenetic inferences has uncovered a clear pattern for the HGT events, which were shown to be concentrated preferentially in pathways or enzymatic steps that are usually reported to be absent in eukaryotes, particularly animals and fungi. Thus, the HGT events identified in this study involve pathways for the synthesis of lysine, cysteine, methionine, threonine, tryptophan, ornithine, and arginine (Figures 1, 2, 3, 5, and 7) and also the synthesis of a few non-essential amino acids such as glycine, serine, and proline. The detailed analysis of these events in different genes and pathways follows.

HGT of homoserine dehydrogenase

Some enzymes are common to a number of pathways involving key precursors to many compounds. Homoserine dehydrogenase (EC:1.1.1.3), for example, participates in the aspartate semialdehyde pathway for the synthesis of lysine, cysteine, methionine, and threonine (Figures 1, 2, and 3). The gene for EC:1.1.1.3 present in symbiont-harboring and regular trypanosomatid genomes seems to have been transferred from a member of the Firmicutes, clustering most closely with Solibacillus silvestris, Lysinibacillus fusiformis, and L. sphaericus with bootstrap support value (BSV) of 100 (Figure 8). On the other hand, the endosymbiont ortholog groups deep within the Betaproteobacteria, more specifically in the Alcaligenaceae family, as expected in the case of no HGT of this gene into the endosymbiont genomes.

HGT and lysine biosynthesis

The two genes of the lysine pathway (Figure 1) that were found in trypanosomatid genomes presented evidence of HGT. H. muscarum was the only trypanosomatid analyzed containing the next to last gene, for diaminopimelate epimerase (EC:5.1.1.7), which clusters strongly with the phylum Bacteroidetes, with BSV of 99 (Additional file 4). The last gene, for diaminopimelate decarboxylase (EC:4.1.1.20), was present in the symbiont-harboring and regular trypanosomatids. In the phylogeny, this particular gene has Actinobacteria as sister group (BSV of 79), although also grouping with a few other eukaryotic genera, most closely Dictyostelium, Polysphondylium, and Capsaspora, with BSV of 65 (Additional file 5). There are, overall, very few Eukaryota in the tree for 4.1.1.20, making it hard to reach a definite conclusion on the direction of transfer for this gene, since other eukaryotes are also present basally to this substantially large group of Actinobacteria plus Trypanosomatidae, with the high BV of 98.

Additional file 4

Maximum likelihood phylogeny of diaminopimelate epimerase (EC:5.1.1.7). Overall tree colored according to taxonomic affiliation of sequences. Values on nodes represent bootstrap support (only 50 or greater shown) and distance bar only applies to the overall tree and not to the detailed regions.

Click here for file

Additional file 5

Maximum likelihood phylogeny of diaminopimelate decarboxylase (EC:4.1.1.20). Overall tree colored according to taxonomic affiliation of sequences. Values on nodes represent bootstrap support (only 50 or greater shown) and distance bar only applies to the overall tree and not to the detailed regions.

Click here for file

Using the C. acanthocephali gene for EC:4.1.1.20 in a manual search against the L. major genome has shown a small fragment with significant similarity (57% identity and 67% similarity, from amino acid 177 to 227), but containing stop codons. Search against predicted L. major proteins yielded no results. These sequence remains suggest that Leishmania could have lost DAP-decarboxylase in a relatively recent past.

HGT and methionine and cysteine biosynthesis

The pathways for cysteine and methionine synthesis (Figure 2) present the highest number of HGT events identified among the pathways studied here. The gene for the enzyme EC:2.3.1.30, necessary for the conversion of serine to cysteine, seems to have been transferred from Bacteria to the genomes of host trypanosomatids. EC:2.3.1.30 of symbiont-harboring and regular trypanosomatids grouped inside a large cluster of diverse Bacteria (predominantly Bacteroidetes and Betaproteobacteria), with BSV of 80 (Additional file 6). An even deeper branch, which separates the subtree containing the trypanosomatids from the rest of the tree, has BSV of 97. The evolutionary history of the other enzyme with the same functionality, EC:2.5.1.47, is unclear and can not be considered a case of HGT given current results. Its gene is present in symbiont-harboring and regular trypanosomatids (including one sequence from T. cruzi CL Brener) and clusters as a sister group of Actinobacteria, although with low BSV (Additional file 7). Although there are many other eukaryotes in the tree, they are not particularly close to the subtree containing the Trypanosomatidae. Interestingly, one Entamoeba dispar sequence is a sister group to the Trypanosomatidae, although with low BSV, raising the possibility of eukaryote-to-eukaryote HGT, as previously reported (reviewed in 77).

Additional file 6

Maximum likelihood phylogeny of serine O-acetyltransferase (EC:2.3.1.30). Overall tree colored according to taxonomic affiliation of sequences. Values on nodes represent bootstrap support (only 50 or greater shown) and distance bar only applies to the overall tree and not to the detailed regions.

Click here for file

Additional file 7

Maximum likelihood phylogeny of cysteine synthase (EC:2.5.1.47). Overall tree colored according to taxonomic affiliation of sequences. Values on nodes represent bootstrap support (only 50 or greater shown) and distance bar only applies to the overall tree and not to the detailed regions.

Click here for file

The gene for EC:2.3.1.46, the first in the pathway converting homoserine to cystathionine, is present in all symbiont-harboring trypanosomatids and Herpetomonas, but in no other regular trypanosomatid examined. This trypanosomatid gene groups within Bacteroidetes, with BSV of 53 and, in a deeper branch, BSV of 89, still clustering with Bacteroidetes only (Additional file 8).

Additional file 8

Maximum likelihood phylogeny of homoserine O-succinyltransferase (EC:2.3.1.46). Overall tree colored according to taxonomic affiliation of sequences. Values on nodes represent bootstrap support (only 50 or greater shown) and distance bar only applies to the overall tree and not to the detailed regions.

Click here for file

The gene for EC:2.1.1.37, responsible for the first step in the conversion of S-adenosylmethionine into homocysteine, is present in all symbiont-harboring and regular trypanosomatids, although the sequence is still partial in the genome sequences of the Angomonas species. Almost all organisms in the tree are Bacteria of several different phyla (Additional file 9), with the few Eukaryota present forming a weakly supported clade. KEGG shows that many Eukaryota do possess a gene for enzyme EC:2.1.1.37, but their sequences are very different from that present in the trypanosomatids (and other eukaryotes) studied here. This therefore suggests a bacterial origin for the EC:2.1.1.37 from the Eukaryota in our tree, although the specific donor group cannot be currently determined with confidence. It is interesting to note that, besides the Trypanosomatidae, the clade of eukaryotes is composed of Stramenopiles and green algae (both groups that have, or once had, plastids), with a Cyanobacteria close to the base of the group. Although the BSV of 54 does not allow strong conclusions regarding this group, it is interesting to speculate about the possibility of eukaryote-to-eukaryote gene transfer, as previously seen (reviewed in 77), after the acquisition of this gene from a so-far unidentified bacterium.

Additional file 9

Maximum likelihood phylogeny of DNA (cytosine-5-)-methyltransferase (EC:2.1.1.37). Overall tree colored according to taxonomic affiliation of sequences. Values on nodes represent bootstrap support (only 50 or greater shown) and distance bar only applies to the overall tree and not to the detailed regions.

Click here for file

The genes for EC:2.5.1.48, EC:2.5.1.49, and EC:4.4.1.8 (two versions) are quite similar in sequence and domain composition. Therefore, similarity searches with any one of these genes also retrieves the other three. In spite of the similarities, these genes are found in rather different phyletic and phylogenetic patterns on the trypanosomatids (Additional file 10). EC:2.5.1.48 is present in all symbiont-harboring and regular trypanosomatids examined, plus Trypanosoma sp. and a few other Eukaryota (mostly Apicomplexa and Stramenopiles), all within a group of Acidobacteria (BSV of 94). The gene for EC:2.5.1.49 is present in the symbiont-harboring trypanosomatids and Herpetomonas, but in none of the other regular trypanosomatids examined. This trypanosomatid gene also clusters with diverse groups of Bacteria, although low BSV makes it hard to confidently identify its most likely nearest neighbor, and it is not possible to conclude with reasonable certainty that this gene is derived from HGT. The gene for EC:4.4.1.8 occurs, in symbiont-harboring and regular trypanosomatids, as two orthologs presenting very different evolutionary histories. One of the orthologs clusters with eukaryotes, with BSV of 95, while the other seems to be of bacterial descent, grouping mostly with Alphaproteobacteria of the Rhizobiales order, with BSV of 99.

Additional file 10

Maximum likelihood phylogeny of cystathionine gamma-synthase, O-acetylhomoserine aminocarboxypropyltransferase, and cystathionine beta-lyase (EC:2.5.1.48, EC:2.5.1.49, and EC:4.4.1.8). Overall tree colored according to taxonomic affiliation of sequences. Values on nodes represent bootstrap support (only 50 or greater shown) and distance bar only applies to the overall tree and not to the detailed regions.

Click here for file

The presence of two genes identified as EC:4.4.1.8 raises the possibility of them performing different enzymatic reactions. Given the overall domain composition similarities of several of the genes of the methionine and cysteine synthesis pathways, it is possible that one of the enzymes identified as EC:4.4.1.8 is actually the enzyme EC:4.4.1.1, for which no gene has been found in our searches of the Trypanosomatidae genomes, as detailed above (Methionine and cysteine).

Genes for two of the enzymes for the last step in the methionine synthesis, EC:2.1.1.10 and EC:2.1.1.14 (Additional files 11 and 12), are present in all regular and symbiont-harboring trypanosomatids (except for Herpetomonas, which lacks the latter). EC:2.1.1.14 appears to be of bacterial origin, grouping within the Gammaproteobacteria with moderate (74) bootstrap support. While EC:2.1.1.10 also groups near Gammaproteobacteria, BSV is low and this gene cannot be considered a case of HGT given current data.

Additional file 11

Maximum likelihood phylogeny of homocysteine S-methyltransferase (EC:2.1.1.10). Overall tree colored according to taxonomic affiliation of sequences. Values on nodes represent bootstrap support (only 50 or greater shown) and distance bar only applies to the overall tree and not to the detailed regions.

Click here for file

Additional file 12

Maximum likelihood phylogeny of 5-methyltetrahydropteroyltriglutamate-homocysteine S-methyltransferase (EC:2.1.1.14). Overall tree colored according to taxonomic affiliation of sequences. Values on nodes represent bootstrap support (only 50 or greater shown) and distance bar only applies to the overall tree and not to the detailed regions.

Click here for file

As seen above, most genes in the de novo methionine synthesis pathway seem to have originated in one or more HGT events. Enzymes from the methionine salvage pathway (Additional file 3), on the other hand, are notably different. Of these, only S-methyl-5-thioribose kinase (EC:2.7.1.100), found in C. acanthocephali and Herpetomonas but not in the endosymbionts and symbiont-harboring trypanosomatids, seems to have originated in a bacterial group (Additional file 13). These two organisms’ enzymes group deep within the Gammaproteobacteria, with BSV of 97.

Additional file 13

Maximum likelihood phylogeny of S-methyl-5-thioribose kinase (EC:2.7.1.100). Overall tree colored according to taxonomic affiliation of sequences. Values on nodes represent bootstrap support (only 50 or greater shown) and distance bar only applies to the overall tree and not to the detailed regions.

Click here for file

Enzyme acireductone synthase (EC:3.1.3.77) presents an intriguing case, being the only methionine salvage pathway enzyme absent from the symbiont-harboring trypanosomatid genomes. This enzyme is of eukaryotic origin (not shown), and present in both H. muscarum and C. acanthocephali, but was not found in any other of the regular trypanosomatids available from KEGG. Interestingly, KEGG data for Trypanosoma brucei also shows the two enzymes preceding EC:3.1.3.77 as missing, which raises the question of whether this important pathway is in the process of being lost in trypanosomatids. If that is not the case, and given that all other enzymes from the pathway are present, the Trypanosomatidae must have a different enzyme (or enzymes) to perform the required reactions.

HGT and threonine biosynthesis

The gene for the enzyme that interconverts glycine and threonine (Figure 3), EC:4.1.2.5, was identified in all symbiont-harboring and regular trypanosomatids (except Herpetomonas), but the evolutionary histories of symbiont-harboring and regular trypanosomatid genes are very different (Additional file 14). The gene found in the regular trypanosomatids Leishmania sp. and C. acanthocephali groups deep within the Firmicutes, most closely Clostridium, with BSV of 63. The symbiont-harboring trypanosomatid genes, on the other hand, cluster as the most basal clade of one of the two large assemblages of eukaryotes present in this phylogeny, although all BSV are low and there is a large group of Bacteria from diverse phyla between the symbiont-harboring trypanosomatids plus a few other eukaryotic groups and the other eukaryotes in this part of the tree. It is therefore difficult to conclude whether the symbiont-harboring trypanosomatid gene is of bacterial or eukaryotic origin.

Additional file 14

Maximum likelihood phylogeny of L-threonine aldolase (EC:4.1.2.5). Overall tree colored according to taxonomic affiliation of sequences. Values on nodes represent bootstrap support (only 50 or greater shown) and distance bar only applies to the overall tree and not to the detailed regions.

Click here for file

HGT and tryptophan biosynthesis

Tryptophan synthase beta subunit (EC:4.2.1.20), present in the symbiont-harboring trypanosomatids and Herpetomonas, is the last enzyme of the tryptophan biosynthesis pathway, and the only one present in trypanosomatids for this pathway. Its gene groups robustly (BSV of 97) with the Bacteroidetes phylum (Additional file 15). It is also highly similar (around 80% identity and 90% similarity) to the corresponding genes of this phylum, suggesting either a very recent transfer or high sequence conservation. Given that the protein alignment of the orthologs (not shown) presents a maximum patristic distance value of 84.04% and a median of 47.22%, it is therefore likely that the transfer of EC:4.2.1.20 to the Trypanosomatidae is relatively recent.

Additional file 15

Maximum likelihood phylogeny of (EC:4.2.1.20). Overall tree colored according to taxonomic affiliation of sequences. Values on nodes represent bootstrap support (only 50 or greater shown) and distance bar only applies to the overall tree and not to the detailed regions.

Click here for file

HGT and arginine and ornithine biosynthesis

The arginine and ornithine synthesis pathway has been influenced by HGT events in a few key steps. As discussed above, one of the entry points for the urea cycle is through ornithine synthesized from glutamate. The last step, converting N-acetylornithine to ornithine, can be performed by either EC:3.5.1.14 or EC:3.5.1.16 (Figure 7). We have found that the genes for both enzymes, present in all symbiont-harboring and regular trypanosomatid genomes, originated from HGT events. All gene copies for EC:3.5.1.14 group as one clade with a gammaproteobacterium (BSV of 98), and with Bacteria of different phyla (predominantly Firmicutes) as nearest sister group, although with low BSV (Additional file 16). The few other eukaryotic groups present in the tree are very distant from the trypanosomatid group. The multiple copies of the gene for EC:3.5.1.16 in symbiont-harboring and regular trypanosomatids group together in a monophyletic clade (Additional file 17), which clusters within a large group of mostly Betaproteobacteria with BSV of 80, including the Alcaligenaceae, the family to which the endosymbionts belong. However, it seems highly unlikely that this sequence has been transferred from the endosymbiont genomes to their hosts genomes because the nuclear sequences are firmly removed from the Alcaligenaceae, and many regular trypanosomatids (including Trypanosoma spp.) also present this gene in the same part of the tree.

Additional file 16

Maximum likelihood phylogeny of aminoacylase (EC:3.5.1.14). Overall tree colored according to taxonomic affiliation of sequences. Values on nodes represent bootstrap support (only 50 or greater shown) and distance bar only applies to the overall tree and not to the detailed regions.

Click here for file

Additional file 17

Maximum likelihood phylogeny of acetylornithine deacetylase (EC:3.5.1.16). Overall tree colored according to taxonomic affiliation of sequences. Values on nodes represent bootstrap support (only 50 or greater shown) and distance bar only applies to the overall tree and not to the detailed regions.

Click here for file

The only trypanosomatid analyzed which presented ornithine carbamoyl transferase (OCT, EC:2.1.3.3) was Herpetomonas muscarum. Our phylogenetic analysis of this gene indicates that it is of eukaryotic origin (not shown). The symbiont-harboring trypanosomatids utilize the OCT provided by their endosymbionts, and their OCT genes group firmly inside the Alcaligenaceae family, next to Taylorella and Advenella, as expected.

The genes for EC:6.3.4.5 and EC:4.3.2.1 present similar evolutionary patterns: both are absent from endosymbiont genomes and present in all symbiont-harboring and regular trypanosomatid genomes – the only exception being the lack of the latter in Leishmania spp. The Trypanosomatid genes form monophyletic groups in their respective trees, grouping within Firmicutes in both cases (Additional files 18 and 19). BSV is higher (82) in the tree of EC:4.3.2.1 than in that of EC:6.3.4.5 (69). In both cases, support falls for deeper branches in the trees. Although the host genomic sequences are still incomplete and in varying degrees of contiguity, it is interesting to note that the genes for EC:6.3.4.5 and EC:4.3.2.1 are present in tandem in one contig in all symbiont-harboring trypanosomatids (Additional file 1). The flanking genes are eukaryotic: terbinafine resistance locus protein and a multidrug resistance ABC transporter. As seen in the genome browser at TriTrypDB (http://tritrypdb.org), Leishmania spp. have most of these same genes, although in a slightly different order (EC:6.3.4.5 occurring after the two eukaryotic genes instead of between them) and lacking EC:4.3.2.1. L. braziliensis seems to be in the process of additionally losing EC:6.3.4.5, which is annotated as a pseudogene. These phylogenetic and genomic data strongly suggest that EC:4.3.2.1 and EC:6.3.4.5 have been transferred together from a Firmicutes bacterium to the common ancestor of the symbiont-harboring and regular trypanosomatids studied, and that these transferred genes have been or are being lost from Leishmania at least.

Additional file 18

Maximum likelihood phylogeny of argininosuccinate synthase (EC:6.3.4.5). Overall tree colored according to taxonomic affiliation of sequences. Values on nodes represent bootstrap support (only 50 or greater shown) and distance bar only applies to the overall tree and not to the detailed regions.

Click here for file

Additional file 19

Maximum likelihood phylogeny of argininosuccinate lyase (EC:4.3.2.1). Overall tree colored according to taxonomic affiliation of sequences. Values on nodes represent bootstrap support (only 50 or greater shown) and distance bar only applies to the overall tree and not to the detailed regions.

Click here for file

The final enzyme in the urea cycle, arginase (EC:3.5.3.1), is present in all symbiont-harboring and regular trypanosomatids examined here. However, the sequence from Herpetomonas presents a partial arginase domain; while the protein sequence length is as expected, the domain match starts only after 70 amino acids. We speculate that this divergence could be responsible for the lack of arginase activity previously seen in Herpetomonas. Differently from most other enzymes in this work, there are different evolutionary histories for the arginase genes: all trypanosomatid genes but that from Herpetomonas cluster together with very high bootstrap support of 98, within Eukaryota (Additional file 20). The sequence from Herpetomonas on the other hand is the sister group (BSV of 79) of a large assemblage of Bacteria from several different phyla, but predominantly Deltaproteobacteria, Firmicutes, Actinobacteria, and Cyanobacteria. It is therefore clear that Herpetomonas must have acquired a different arginase than that present in the other trypanosomatids studied, which possess eukaryotic genes. Furthermore, this gene seems to be undergoing a process of decay, given its lack of significant similarity to the known arginase domain in a significant portion of the protein.

Additional file 20

Maximum likelihood phylogeny of arginase (EC:3.5.3.1). Overall tree colored according to taxonomic affiliation of sequences. Values on nodes represent bootstrap support (only 50 or greater shown) and distance bar only applies to the overall tree and not to the detailed regions.

Click here for file

HGT in other pathways: possible symbiont to host transfer

Ornithine cyclodeaminase (EC:4.3.1.12) converts ornithine directly into proline, a non-essential amino acid. In our analyses, we have found that the gene for EC:4.3.1.12 of symbiont-harboring trypanosomatid genomes is very similar to those from Betaproteobacteria of the Alcaligenaceae family, to which the endosymbionts belong. Regular trypanosomatid and endosymbiont genomes do not contain the gene for this enzyme. Accordingly, the phylogeny shows the symbiont-harboring trypanosomatid gene grouping close to several Alcaligenaceae, although the clade is not monophyletic and presents BSV of 63 (Additional file 21). This grouping, together with the gene presence in symbiont-harboring trypanosomatid genomes only, poses the possibility that EC:4.3.1.12 has been transferred from the ancestral endosymbiont to the corresponding host, before the radiation of symbiont-harboring trypanosomatids into the two genera and five species analyzed here.

Additional file 21

Maximum likelihood phylogeny of ornithine cyclodeaminase (EC:4.3.1.12). Overall tree colored according to taxonomic affiliation of sequences. Values on nodes represent bootstrap support (only 50 or greater shown) and distance bar only applies to the overall tree and not to the detailed regions.

Click here for file

Other observation on amino acid pathway peculiarities

Some interesting peculiarities of specific genes from a few pathways deserve to be discussed. Interestingly, the gene for branched-chain-amino-acid transaminase (EC:2.6.1.42), the last step in the synthesis of isoleucine, valine, and leucine (Figure 4), was identified in all bacteria of the Alcaligenaceae family present in KEGG, except for the endosymbionts’ closest relatives, Taylorella spp. (parasitic) and Advenella kashmirensis (free-living), which also lack the gene. The question is raised then of whether the common ancestor of Taylorella and the endosymbionts, which are sister groups 25, had already lost the gene. Another possibility is that independent loses occurred in endosymbionts, Taylorella, and Advenella. Considering that the rest of the pathway is present in these organisms and that the free-living Advenella would need the last gene to complete synthesis of these amino acids, it is reasonable to speculate that their EC:2.6.1.42 is novel or at least very different and thus could not be identified by similarity searches.

As mentioned above, the histidine pathway biosynthesis is performed by the endosymbionts and all enzymes, with the exception of histidinol-phosphate phosphatase (HPP, EC:3.1.3.15), have been identified. This is also the only enzyme of this pathway missing in other Betaproteobacteria available in KEGG. Recently, it was reported that such a gap in the histidine biosynthesis pathway in other organisms was completed by novel HPP families 7879. Our searches for the novel C. glutamicum HPP (cg0910, an inositol monophosphatase-like gene) have identified two possible candidate genes in the endosymbionts (BCUE_0333 and BCUE_0385, in C. K. blastocrithidii). As in Corynebacterium, neither of these genes is in the same operon as the known histidine synthesis genes. Given the absence of any other inositol phosphate metabolism genes in the endosymbiont genomes, except for these two IMPases, it is reasonable to hypothesize that at least one of the two aforementioned candidates could be the HPP.

Conclusions

In the present paper, we have put together nutritional, biochemical, and genomic data in order to describe how the metabolic co-evolution between the symbiont and the host trypanosomatid is reflected in amino acid production (Figure 9). In fact, amino acid biosynthetic pathways in symbiont-harboring trypanosomatids are frequently chimeras of host and endosymbiont encoded enzymes, with predominance of the latter in the synthesis of essential amino acids. After a careful analysis of different routes, it becomes clear that the symbiotic bacterium completes and/or potentiates most pathways of the host protozoa that are involved in amino acid production, as previously seen in other systems 7.

Figure 9

Overview of the biosynthetic pathways of essential amino acids in trypanosomatids

Overview of the biosynthetic pathways of essential amino acids in trypanosomatids. Dashed arrows: metabolite import; dotted arrows: reaction present in only some of the organisms analyzed; solid arrows: other reactions (a single arrow can summarize multiple steps); arrows surrounded by a gray box: enzymes possibly acquired through horizontal transfer from Bacteria to trypanosomatids (see main text). A. Contribution of symbiont-harboring trypanosomatids (SHT) and endosymbionts (TPE) based on the analysis of gene content in the genomes of A. deanei, A. desouzai, S. culicis, S. oncopelti, S. galati, and respective endosymbionts. B. Biochemical capability of trypanosomatids without symbionts (RT), based on the analysis of genomic data from H. muscarum, C. acanthocephali, and L. major.

Sometimes, as in the lysine and histidine synthesis, the symbionts contain all genes for enzymes that compose the metabolic route. By contrast, in the cysteine and methionine pathways the bacterium lacks most genes involved in amino acid interconversion, which are present in host trypanosomatids. Interestingly, the last step of some metabolic routes such as those for lysine and tryptophan, contains two genes; one in the host genome, the other in the endosymbiont genome. This phenomenon has also been observed in the synthesis of heme 5880, but the reasons for this peculiarity remain obscure. However, we have to consider the possibility that HGT events preceded the colonization of symbiont-harboring trypanosomatids by their endosymbionts, and that the genes present in the host genomes are just relics of previous HGT event(s). Alternatively, these genes could have been recruited to perform functions, as the control of amino acid production by the host trypanosomatid. This same strategy can be considered in isoleucine, valine, and leucine production, but in this case endosymbionts lack the enzyme for the last step, the branched-chain amino acid transaminase (EC:2.6.1.42).

A clear example of the integration of earlier nutritional and enzymatic data with the present gene screening is the synthesis of arginine and ornithine in trypanosomatids. Differently from other members of the family, the urea cycle is complete in symbiont-harboring trypanosomatids by the presence of the OCT gene (EC:2.1.3.3) in symbionts, making these protozoa entirely autotrophic for ornithine, citrulline, and arginine, as previously known from nutritional data 19224452. Symbiont-bearing trypanosomatids contain genes for all enzymes leading from glutamate to arginine. The corresponding genes are located partly in the genomes of their endosymbionts and partly in the protozoan nucleus; in this last case, genes are of bacterial origin, resulting from HGT and including at least one transfer of two genes at once (EC:4.3.2.1 and EC:6.3.4.5), as demonstrated in our phylogenies. Furthermore, endosymbionts also contain most genes for the glutamate pathway, thus enhancing synthesis of ornithine, that once decarboxylated generates polyamine, which is related to cell proliferation and to the low generation time displayed by symbiont-harboring trypanosomatids. Results in this study confirm previous findings 2558 showing the betaproteobacterial origin of the genes of endosymbionts. The nuclear genes, on the other hand, present a much more convoluted evolutionary picture, with probably numerous ancient HGT events shaping the amino acid metabolism in trypanosomatids. A few pathways in particular have been heavily affected, i.e. methionine/cysteine and arginine/ornithine synthesis. Transferred genes originated preferentially from three bacterial phyla, namely Firmicutes, Bacteroidetes, and Gammaproteobacteria, although possible transfers from other phyla of Bacteria have also been uncovered. Especially interesting was the finding of a gene, coding for ornithine cyclodeaminase (EC:4.3.1.12), which closely groups with the Alcaligenaceae family of the Betaprotebacteria and that is likely to have been transferred from the endosymbiont to the host genome. Accordingly, it is present only in symbiont-harboring trypanosomatid nuclear genomes and not in any of the currently sequenced regular trypanosomatid genomes. During review of this work, a very recent report of a similar situation of multiple lineages contributing to the metabolism in the symbiosis of mealybugs, involving the three interacting partners and genes acquired through HGT from other bacterial sources (mainly Alphaproteobacteria, but also Gammaproteobacteria and Bacteroidetes) to the insect host, has been published 81. This suggests that this phenomenon could be widespread and of great importance in genomic and metabolic evolution.

Having been detected in about half of the genes analyzed in this work, HGT events seem to have been fundamental in the genomic evolution of the Trypanosomatidae analyzed, and further phylogenetic studies of the whole host genomes should show the complete extent of this process and which additional pathways could be affected. Synthesis of vitamins (Klein et al., personal communication), heme, and amino acids have already been shown to benefit from bacterial-to-trypanosomatid HGT; many other processes in Trypanosomatidae metabolism might also be subjected to this evolutionary process.

Total genomic DNA was extracted by the phenol-chloroform method 83. We applied kDNA depletion methods to minimize the presence of this type of molecule, as previously described 58, which result in less than about 5% of remaining kDNA in the sample. After kDNA depletion, about 5 μg of DNA were submitted to each Roche 454 shotgun sequencing run, according to the manufacturer’s protocols. Different genomes have so far been sequenced to different levels of draft quality, with estimated coverages of 15X to 23X (considering a genome of ~30 Mbp). Sequences were assembled using the Newbler assembler version 2.3, provided by Roche. Resulting assemblies are available from GenBank, as detailed in “Availability of Supporting Data” below. Endosymbiont genomes were finished to a closed circle as previously described 25.

Gene discovery and annotation

Endosymbiont genes were used as previously published 25. In an initial scan of the genome, trypanosomatid genes were discovered and mapped to metabolic pathways using ASGARD 84, employing as reference the UniRef100 85 and the Kyoto Encyclopedia of Genes and Genomes, KEGG 86 databases. The identified segments of DNA were then extracted from the genome and manually curated for completion and proper location of start and stop codons by using the GBrowse genome browser 87. Putative sequence functions were confirmed by domain searches against NCBI’s Conserved Domain Database 88. Genes and annotations from other trypanosomatids were used when needed and as available at KEGG. All trypanosomatid genes characterized in this study have been submitted to NCBI’s GenBank and accession numbers are available from Additional file 22. All endosymbiont genes analyzed here have been previously sequenced 25; gene identifiers are available from Additional file 23.

Additional file 22

Genbank accession numbers for Trypanosomatidae genes characterized in this study.

Click here for file

Additional file 23

GenBank locus tags for the Ca.Kinetoplastibacterium genes analyzed in this study.

Click here for file

Due to the incomplete nature of our trypanosomatid assemblies, a set of criteria were used to avoid including contaminant sequences in our analyses. A gene was accepted as legitimate only when satisfying at least two of the following: genomic context compatible with a trypanosomatid gene (i.e. long stretches of genes in the same orientation in the contig, most neighboring genes similar to genes from other, previously sequenced trypanosomatids); sequencing coverage in the gene similar to, or higher than, that of the gene and genome averages (since contaminants that are difficult to detect will almost always be in small contigs of low coverage); GC percent content consistent with that of the neighboring genes, and of the overall genome; and phylogenetic congruence (i.e. whether genes from more than one trypanosomatid formed monophyletic assemblages). Genomic context and GC content graphs were drawn by GBrowse 87 and graphically edited for better use of space.

Phylogenetic analyses

For phylogenetic analysis of each enzyme characterized in this work, corresponding putative orthologous genes from all domains of life were collected from the public databases by BLAST search (E-value cutoff of 1e-10, maximum of 10,000 matches accepted) against the full NCBI NR protein database, collecting sequences from as widespread taxonomic groups as possible and keeping one from each species (except for alignments with more than ~1,500 sequences, in which case one organism per genus was kept). Only sequences that were complete and aligned along at least 75% of the length of the query were selected. All analyses were performed at the protein sequence level. Sequences were aligned by Muscle v. 3.8.31 89. Phylogenetic inferences were performed by the maximum likelihood method, using RAxML v. 7.2.8 90 and employing the WAG amino acid substitution model 91, with four gamma-distributed substitution rate heterogeneity categories and empirically determined residue frequencies (model PROTGAMMAWAGF). Each alignment was submitted to bootstrap analysis with 100 pseudo-replicates. Trees were initially drawn and formatted using TreeGraph2 92 and Dendroscope 93, with subsequent cosmetic adjustments performed with the Inkscape vector image editor (http://inkscape.org). Phylogenetic conclusions have been displayed as strong in the summary table for phylogenetic results (Additional file 2) if BSV was 80 or greater, and moderate if BSV was between 50 and 80 – with one exception, EC:2.1.1.37, described in the results. Resulting phylogenetic trees are available from TreeBase (accession number S14564), as detailed in “Availability of Supporting Data” below.

JMPA generated and analyzed metabolic pathway data, performed phylogenetic analyses, and wrote the manuscript. CCK generated and analyzed metabolic pathway data, and edited the manuscript. FMS and AGCM analyzed metabolic pathway data. MGS and GAB generated and analyzed genomic data. ATRV, MFS, and MMGT, participated in study design and data interpretation. MCM and EPC conceived the study and wrote the manuscript. All authors read and approved the manuscript.

Acknowledgments

We would like to thank Marta Campaner and Carmen C. Takata (USP), Vladimir Lee, Andrey Matveyev, and Yingping Wang (VCU) for technical support, and Carlisle G. Childress Jr. and J. Michael Davis (VCU Center for High-Performance Computing). Sequencing was performed in the Nucleic Acids Research Facilities, and analyses were performed in the Bioinformatics Computational Core Lab and the Center for High Performance Computing at VCU. The research leading to these results was funded by: the European Research Council under the European Community’s Seventh Framework Programme (FP7/2007-2013)/ERC grant agreement n° [247073]10; the French project ANR MIRI BLAN08-1335497; FAPERJ grant coordinated by Dr. Cristina Motta and the FAPERJ-INRIA project RAMPA; by the Laboratoire International Associé (LIA) LIRIO co-coordinated by Ana Tereza R. de Vasconcelos (Labinfo, LNCC, Brazil) and Marie-France Sagot (LBBE, UCBL-CNRS-INRIA, France); the National Science Foundation [USA, grant number NSF DEB-0830056 to Gregory Buck]; and Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq, Brazil) and Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES, Brazil) to Cristina Motta, Marta M.G. Teixeira and Erney P. Camargo.

Nutritional requirements of Blastocrithidia culicis, a trypanosomatid with an endosymbiont

MenezesMCNDRoitmanIJ Eukaryotic Microbiol199138122123

Phylogenetic validation of the genera Angomonas and Strigomonas of trypanosomatids harboring bacterial endosymbionts with the description of New species of trypanosomatids and of proteobacterial symbionts