Baseline gene signatures of reactogenicity to Ebola vaccination: a machine learning approach across multiple cohorts.

Introduction: The rVSVDG-ZEBOV-GP (Ervebo®) vaccine is both immunogenic and protective against Ebola. However, the vaccine can cause a broad range of transient adverse reactions, from headache to arthritis. Identifying baseline reactogenicity signatures can advance personalized vaccinology and increase our understanding of the molecular factors associated with such adverse events. Methods: In this study, we developed a machine learning approach to integrate prevaccination gene expression data with adverse events that occurred within 14 days post-vaccination. Results and Discussion: We analyzed the expression of 144 genes across 343 blood samples collected from participants of 4 phase I clinical trial cohorts: Switzerland, USA, Gabon, and Kenya. Our machine learning approach revealed 22 key genes associated with adverse events such as local reactions, fatigue, headache, myalgia, fever, chills, arthralgia, nausea, and arthritis, providing insights into potential biological mechanisms linked to vaccine reactogenicity.

Erratum: Prior upregulation of interferon pathways in the nasopharynx impacts viral shedding following live attenuated influenza vaccine challenge in children.

[This corrects the article DOI: 10.1016/j.xcrm.2021.100465.].

Comparison of carbohydrate ABC importers from Mycobacterium tuberculosis.

BACKGROUND: Mycobacterium tuberculosis, the etiological agent of tuberculosis, has at least four ATP-Binding Cassette (ABC) transporters dedicated to carbohydrate uptake: LpqY/SugABC, UspABC, Rv2038c-41c, and UgpAEBC. LpqY/SugABC transporter is essential for M. tuberculosis survival in vivo and potentially involved in the recycling of cell wall components. The three-dimensional structures of substrate-binding proteins (SBPs) LpqY, UspC, and UgpB were described, however, questions about how these proteins interact with the cognate transporter are still being explored. Components of these transporters, such as SBPs, show high immunogenicity and could be used for the development of diagnostic and therapeutic tools. In this work, we used a phylogenetic and structural bioinformatics approach to compare the four systems, in an attempt to predict functionally important regions. RESULTS: Through the analysis of the putative orthologs of the carbohydrate ABC importers in species of Mycobacterium genus it was shown that Rv2038c-41c and UgpAEBC systems are restricted to pathogenic species. We showed that the components of the four ABC importers are phylogenetically separated into four groups defined by structural differences in regions that modulate the functional activity or the interaction with domain partners. The regulatory region in nucleotide-binding domains, the periplasmic interface in transmembrane domains and the ligand-binding pocket of the substrate-binding proteins define their substrates and segregation in different branches. The interface between transmembrane domains and nucleotide-binding domains show conservation of residues and charge. CONCLUSIONS: The presence of four ABC transporters in M. tuberculosis dedicated to uptake and transport of different carbohydrate sources, and the exclusivity of at least two of them being present only in pathogenic species of Mycobacterium genus, highlights their relevance in virulence and pathogenesis. The significant differences in the SBPs, not present in eukaryotes, and in the regulatory region of NBDs can be explored for the development of inhibitory drugs targeting the bacillus. The possible promiscuity of NBDs also contributes to a less specific and more comprehensive control approach.

In-depth analysis of laboratory parameters reveals the interplay between sex, age, and systemic inflammation in individuals with COVID-19.

BACKGROUND: The progression and severity of COVID-19 vary significantly in the population. While the hallmarks of SARS-CoV-2 and severe COVID-19 within routine laboratory parameters are emerging, the impact of sex and age on these profiles is still unknown. METHODS: A multidimensional analysis was performed involving millions of records of laboratory parameters and diagnostic tests for 178 887 individuals from Brazil, of whom 33 266 tested positive for SARS-CoV-2. Analyzed data included those relating to complete blood cell count, electrolytes, metabolites, arterial blood gases, enzymes, hormones, cancer biomarkers, and others. FINDINGS: COVID-19 induced similar alterations in laboratory parameters in males and females. CRP and ferritin were increased, especially in older men with COVID-19, whereas abnormal liver function tests were common across several age groups, except for young women. Low peripheral blood basophils and eosinophils were more common in the elderly with COVID-19. Both male and female COVID-19 patients admitted to intensive care units displayed alterations in the coagulation system, and higher values for neutrophils, CRP, and lactate dehydrogenase. CONCLUSIONS: Our study uncovered the laboratory profiles of a large cohort of COVID-19 patients, which formed the basis of discrepancies influenced by aging and biological sex. These profiles directly linked COVID-19 disease presentation to an intricate interplay between sex, age, and immune activation.

Prior upregulation of interferon pathways in the nasopharynx impacts viral shedding following live attenuated influenza vaccine challenge in children.

In children lacking influenza-specific adaptive immunity, upper respiratory tract innate immune responses may influence viral replication and disease outcome. We use trivalent live attenuated influenza vaccine (LAIV) as a surrogate challenge model in children aged 24-59 months to identify pre-infection mucosal transcriptomic signatures associated with subsequent viral shedding. Upregulation of interferon signaling pathways prior to LAIV is significantly associated with lower strain-specific viral loads (VLs) at days 2 and 7. Several interferon-stimulated genes are differentially expressed in children with pre-LAIV asymptomatic respiratory viral infections and negatively correlated with LAIV VLs. Upregulation of genes enriched in macrophages, neutrophils, and eosinophils is associated with lower VLs and found more commonly in children with asymptomatic viral infections. Variability in pre-infection mucosal interferon gene expression in children may impact the course of subsequent influenza infections. This variability may be due to frequent respiratory viral infections, demonstrating the potential importance of mucosal virus-virus interactions in children.

Network vaccinology.

The structure and function of the immune system is governed by complex networks of interactions between cells and molecular components. Vaccination perturbs these networks, triggering specific pathways to induce cellular and humoral immunity. Systems vaccinology studies have generated vast data sets describing the genes related to vaccination, motivating the use of new approaches to identify patterns within the data. Here, we describe a framework called Network Vaccinology to explore the structure and function of biological networks responsible for vaccine-induced immunity. We demonstrate how the principles of graph theory can be used to identify modules of genes, proteins, and metabolites that are associated with innate and adaptive immune responses. Network vaccinology can be used to assess specific and shared molecular mechanisms of different types of vaccines, adjuvants, and routes of administration to direct rational design of the next generation of vaccines.

Genomic Organization and Generation of Genetic Variability in the RHS (Retrotransposon Hot Spot) Protein Multigene Family in Trypanosoma cruzi.

Retrotransposon Hot Spot (RHS) is the most abundant gene family in Trypanosoma cruzi, with unknown function in this parasite. The aim of this work was to shed light on the organization and expression of RHS in T. cruzi. The diversity of the RHS protein family in T. cruzi was demonstrated by phylogenetic and recombination analyses. Transcribed sequences carrying the RHS domain were classified into ten distinct groups of monophyletic origin. We identified numerous recombination events among the RHS and traced the origins of the donors and target sequences. The transcribed RHS genes have a mosaic structure that may contain fragments of different RHS inserted in the target sequence. About 30% of RHS sequences are located in the subtelomere, a region very susceptible to recombination. The evolution of the RHS family has been marked by many events, including gene duplication by unequal mitotic crossing-over, homologous, as well as ectopic recombination, and gene conversion. The expression of RHS was analyzed by immunofluorescence and immunoblotting using anti-RHS antibodies. RHS proteins are evenly distributed in the nuclear region of T. cruzi replicative forms (amastigote and epimastigote), suggesting that they could be involved in the control of the chromatin structure and gene expression, as has been proposed for T. brucei.

Mapping benznidazole resistance in trypanosomatids and exploring evolutionary histories of nitroreductases and ABCG transporter protein sequences.

The nitro-heterocyclic compound benznidazole (BZ) is the first-line drug for the treatment of Chagas disease, caused by the protozoan Trypanosoma cruzi. However, therapeutic failures are common for reasons that include the influences of parasite and host genetics, the effects of toxicity on adherence to treatment, and difficulties in demonstrating parasitological cure. To obtain information on the origin of the resistance to BZ and eliminate from the scenery the participation of the host, initially we mapped the susceptibility to the drug in thirteen species of seven genera of the family Trypanosomatidae. We verified that all Trypanosoma species are sensitive to low concentrations of the drug (IC50 2.7 to 25 µM) while Non-Trypanosoma species are highly resistant to these concentrations. The two groups of parasites correspond to the major phylogenetic lineages of trypanosomatids. Next, we searched in the trypanosomatid genome databases homologs of two type-I nitroreductases (NTR-1 and OYE) and an ABC transporter (ABCG1) that have been associated with BZ resistance in T. cruzi. The predicted proteins were characterized regarding domains and used for phylogenetic analyses. Homologous NTR-1 genes were found in all trypanosomatids investigated and the structural characteristics of the enzyme suggest that it may be functional. OYE genes were absent in BZ-sensitive African trypanosomes, which excludes the participation of this enzyme in BZ bio-activation. Two copies of ABCG1 genes were observed in most BZ resistant species, while Trypanosoma species exhibit only one copy per haploid genome. Functional studies are required to verify the involvement of these genes in BZ resistance. In addition, since multiple mechanisms can contribute to BZ susceptibility, our study poses a range of organisms highly resistant to BZ in which these aspects can be investigated. Preliminary studies on BZ uptake indicate marked differences between BZ-sensitive and BZ-resistant species.

Genomic comparison of Trypanosoma conorhini and Trypanosoma rangeli to Trypanosoma cruzi strains of high and low virulence.

BACKGROUND: Trypanosoma conorhini and Trypanosoma rangeli, like Trypanosoma cruzi, are kinetoplastid protist parasites of mammals displaying divergent hosts, geographic ranges and lifestyles. Largely nonpathogenic T. rangeli and T. conorhini represent clades that are phylogenetically closely related to the T. cruzi and T. cruzi-like taxa and provide insights into the evolution of pathogenicity in those parasites. T. rangeli, like T. cruzi is endemic in many Latin American countries, whereas T. conorhini is tropicopolitan. T. rangeli and T. conorhini are exclusively extracellular, while T. cruzi has an intracellular stage in the mammalian host. RESULTS: Here we provide the first comprehensive sequence analysis of T. rangeli AM80 and T. conorhini 025E, and provide a comparison of their genomes to those of T. cruzi G and T. cruzi CL, respectively members of T. cruzi lineages TcI and TcVI. We report de novo assembled genome sequences of the low-virulent T. cruzi G, T. rangeli AM80, and T. conorhini 025E ranging from ~ 21-25 Mbp, with ~ 10,000 to 13,000 genes, and for the highly virulent and hybrid T. cruzi CL we present a ~ 65 Mbp in-house assembled haplotyped genome with ~ 12,500 genes per haplotype. Single copy orthologs of the two T. cruzi strains exhibited ~ 97% amino acid identity, and ~ 78% identity to proteins of T. rangeli or T. conorhini. Proteins of the latter two organisms exhibited ~ 84% identity. T. cruzi CL exhibited the highest heterozygosity. T. rangeli and T. conorhini displayed greater metabolic capabilities for utilization of complex carbohydrates, and contained fewer retrotransposons and multigene family copies, i.e. trans-sialidases, mucins, DGF-1, and MASP, compared to T. cruzi. CONCLUSIONS: Our analyses of the T. rangeli and T. conorhini genomes closely reflected their phylogenetic proximity to the T. cruzi clade, and were largely consistent with their divergent life cycles. Our results provide a greater context for understanding the life cycles, host range expansion, immunity evasion, and pathogenesis of these trypanosomatids.

Genome-wide identification of evolutionarily conserved Small Heat-Shock and eight other proteins bearing α-crystallin domain-like in kinetoplastid protists.

Small Heat-Shock Proteins (sHSPs) and other proteins bearing alpha-crystallin domains (ACD) participate in defense against heat and oxidative stress and play important roles in cell cycle, cytoskeleton dynamics, and immunological and pathological mechanisms in eukaryotes. However, little is known about these proteins in early-diverging lineages of protists such as the kinetoplastids. Here, ACD-like proteins (ACDp) were investigated in genomes of 61 species of 12 kinetoplastid genera, including Trypanosoma spp. (23 species of mammals, reptiles and frogs), Leishmania spp. (mammals and lizards), trypanosomatids of insects, Phytomonas spp. of plants, and bodonids. Comparison of ACDps based on domain architecture, predicted tertiary structure, phylogeny and genome organization reveals a kinetoplastid evolutionarily conserved repertoire, which diversified prior to trypanosomatid adaptation to parasitic life. We identified 9 ACDp orthologs classified in 8 families of TryACD: four previously recognized (HSP20, Tryp23A, Tryp23B and ATOM69), and four characterized for the first time in kinetoplastids (TryACDP, TrySGT1, TryDYX1C1 and TryNudC). A single copy of each ortholog was identified in each genome alongside TryNudC1/TrypNudC2 homologs and, overall, ACDPs were under strong selection pressures at main phylogenetic lineages. Transcripts of all ACDPs were identified across the life stages of T. cruzi, T. brucei and Leishmania spp., but proteomic profiles suggested that most ACDPs may be species- and stage-regulated. Our findings establish the basis for functional studies, and provided evolutionary and structural support for an underestimated repertoire of ACDps in the kinetoplastids.

Trypanosoma rangeli is phylogenetically closer to Old World trypanosomes than to Trypanosoma cruzi.

Trypanosoma rangeli and Trypanosoma cruzi are generalist trypanosomes sharing a wide range of mammalian hosts; they are transmitted by triatomine bugs, and are the only trypanosomes infecting humans in the Neotropics. Their origins, phylogenetic relationships, and emergence as human parasites have long been subjects of interest. In the present study, taxon-rich analyses (20 trypanosome species from bats and terrestrial mammals) using ssrRNA, glycosomal glyceraldehyde-3-phosphate dehydrogenase (gGAPDH), heat shock protein-70 (HSP70) and Spliced Leader RNA sequences, and multilocus phylogenetic analyses using 11 single copy genes from 15 selected trypanosomes, provide increased resolution of relationships between species and clades, strongly supporting two main sister lineages: lineage Schizotrypanum, comprising T. cruzi and bat-restricted trypanosomes, and Tra[Tve-Tco] formed by T. rangeli, Trypanosoma vespertilionis and Trypanosoma conorhini clades. Tve comprises European T. vespertilionis and African T. vespertilionis-like of bats and bat cimicids characterised in the present study and Trypanosoma sp. Hoch reported in monkeys and herein detected in bats. Tco included the triatomine-transmitted tropicopolitan T. conorhini from rats and the African NanDoum1 trypanosome of civet (carnivore). Consistent with their very close relationships, Tra[Tve-Tco] species shared highly similar Spliced Leader RNA structures that were highly divergent from those of Schizotrypanum. In a plausible evolutionary scenario, a bat trypanosome transmitted by cimicids gave origin to the deeply rooted Tra[Tve-Tco] and Schizotrypanum lineages, and bat trypanosomes of diverse genetic backgrounds jumped to new hosts. A long and independent evolutionary history of T. rangeli more related to Old World trypanosomes from bats, rats, monkeys and civets than to Schizotrypanum spp., and the adaptation of these distantly related trypanosomes to different niches of shared mammals and vectors, is consistent with the marked differences in transmission routes, life-cycles and host-parasite interactions, resulting in T. cruzi (but not T. rangeli) being pathogenic to humans.

New insights from Gorongosa National Park and Niassa National Reserve of Mozambique increasing the genetic diversity of Trypanosoma vivax and Trypanosoma vivax-like in tsetse flies, wild ungulates and livestock from East Africa.

BACKGROUND: Trypanosoma (Duttonella) vivax is a major pathogen of livestock in Africa and South America (SA), and genetic studies limited to small sampling suggest greater diversity in East Africa (EA) compared to both West Africa (WA) and SA. METHODS: Multidimensional scaling and phylogenetic analyses of 112 sequences of the glycosomal glyceraldehyde phosphate dehydrogenase (gGAPDH) gene and 263 sequences of the internal transcribed spacer of rDNA (ITS rDNA) were performed to compare trypanosomes from tsetse flies from Gorongosa National Park and Niassa National Reserve of Mozambique (MZ), wild ungulates and livestock from EA, and livestock isolates from WA and SA. RESULTS: Multidimensional scaling (MDS) supported Tvv (T. vivax) and TvL (T. vivax-like) evolutionary lineages: 1) Tvv comprises two main groups, TvvA/B (all SA and WA isolates plus some isolates from EA) and TvvC/D (exclusively from EA). The network revealed five ITS-genotypes within Tvv: Tvv1 (WA/EA isolates), Tvv2 (SA) and Tvv3-5 (EA). EA genotypes of Tvv ranged from highly related to largely different from WA/SA genotypes. 2) TvL comprises two gGAPDH-groups formed exclusively by EA sequences, TvLA (Tanzania/Kenya) and TvLB-D (MZ). This lineage contains more than 11 ITS-genotypes, seven forming the lineage TvL-Gorongosa that diverged from T. vivax Y486 enough to be identified as another species of the subgenus Duttonella. While gGAPDH sequences were fundamental for classification at the subgenus, major evolutionary lineages and species levels, ITS rDNA sequences permitted identification of known and novel genotypes. CONCLUSIONS: Our results corroborate a remarkable diversity of Duttonella trypanosomes in EA, especially in wildlife conservation areas, compared to the moderate diversity in WA. Surveys in wilderness areas in WA may reveal greater diversity. Biogeographical and phylogenetic data point to EA as the place of origin, diversification and spread of Duttonella trypanosomes across Africa, providing relevant insights towards the understanding of T. vivax evolutionary history.

The Evolutionary Loss of RNAi Key Determinants in Kinetoplastids as a Multiple Sporadic Phenomenon.

We screened the genomes of a broad panel of kinetoplastid protists for genes encoding proteins associated with the RNA interference (RNAi) system using probes from the Argonaute (AGO1), Dicer1 (DCL1), and Dicer2 (DCL2) genes of Leishmania brasiliensis and Crithidia fasciculata. We identified homologs for all the three of these genes in the genomes of a subset of these organisms. However, several of these organisms lacked evidence for any of these genes, while others lacked only DCL2. The open reading frames encoding these putative proteins were structurally analyzed in silico. The alignments indicated that the genes are homologous with a high degree of confidence, and three-dimensional structural models strongly supported a functional relationship to previously characterized AGO1, DCL1, and DCL2 proteins. Phylogenetic analysis of these putative proteins showed that these genes, when present, evolved in parallel with other nuclear genes, arguing that the RNAi system genes share a common progenitor, likely across all Kinetoplastea. In addition, the genome segments bearing these genes are highly conserved and syntenic, even among those taxa in which they are absent. However, taxa in which these genes are apparently absent represent several widely divergent branches of kinetoplastids, arguing that these genes were independently lost at least six times in the evolutionary history of these organisms. The mechanisms responsible for the apparent coordinate loss of these RNAi system genes independently in several lineages of kinetoplastids, while being maintained in other related lineages, are currently unknown.

Phylogenetic and syntenic data support a single horizontal transference to a Trypanosoma ancestor of a prokaryotic proline racemase implicated in parasite evasion from host defences.

BACKGROUND: Proline racemase (PRAC) enzymes of Trypanosoma cruzi (TcPRAC), the agent of Chagas disease, and Trypanosoma vivax (TvPRAC), the agent of livestock trypanosomosis, have been implicated in the B-cells polyclonal activation contributing to immunosuppression and the evasion of host defences. The similarity to prokaryotic PRAC and the absence in Trypanosoma brucei and Trypanosoma congolense have raised many questions about the origin, evolution, and functions of trypanosome PRAC (TryPRAC) enzymes. FINDINGS: We identified TryPRAC homologs as single copy genes per haploid genome in 12 of 15 Trypanosoma species, including T. cruzi and T. cruzi marinkellei, T. dionisii, T. erneyi, T. rangeli, T. conorhini and T. lewisi, all parasites of mammals. Polymorphisms in TcPRAC genes matched T. cruzi genotypes: TcI-TcIV and Tcbat have unique genes, while the hybrids TcV and TcVI contain TcPRACA and TcPRACB from parental TcII and TcIII, respectively. PRAC homologs were identified in trypanosomes from anurans, snakes, crocodiles, lizards, and birds. Most trypanosomes have intact PRAC genes. T. rangeli possesses only pseudogenes, maybe in the process of being lost. T. brucei, T. congolense and their allied species, except the more distantly related T. vivax, have completely lost PRAC genes. CONCLUSIONS: The genealogy of TryPRAC homologs supports an evolutionary history congruent with the Trypanosoma phylogeny. This finding, together with the synteny of PRAC loci, the relationships with prokaryotic PRAC inferred by taxon-rich phylogenetic analysis, and the absence in trypanosomatids of any other genera or in bodonids or euglenids suggest that a common ancestor of Trypanosoma gained PRAC gene by a single and ancient horizontal gene transfer (HGT) from a Firmicutes bacterium more closely related to Gemella and other species of Bacilli than to Clostridium as previously suggested. Our broad phylogenetic study allowed investigation of TryPRAC evolution over long and short timescales. TryPRAC genes diverged to become species-specific and genotype-specific for T. cruzi and T. rangeli, with resulting genealogies congruent with those obtained using vertically inherited genes. The inventory of TryPRAC genes described here is the first step toward the understanding of the roles of PRAC enzymes in trypanosomes differing in life cycles, virulence, and infection and immune evasion strategies.

Congopain genes diverged to become specific to Savannah, Forest and Kilifi subgroups of Trypanosoma congolense, and are valuable for diagnosis, genotyping and phylogenetic inferences.

Trypanosoma congolense is the most important agent of nagana, a wasting livestock trypanosomosis in sub-Saharan Africa. This species is a complex of three subgroups (Savannah, Forest and Kilifi) that differ in virulence, pathogenicity, drug resistance, vectors, and geographical distribution. Congopain, the major Cathepsin L-like cysteine protease (CP2) of T. congolense, has been extensively investigated as a pathogenic factor and target for drugs and vaccines, but knowledge about this enzyme is mostly restricted to the reference strain IL3000, which belongs to the Savannah subgroup. In this work we compared sequences of congopain genes from IL3000 genome database and isolates of the three subgroups of T. congolense. Results demonstrated that the congopain genes diverged into three subclades consistent with the three subgroups within T. congolense. Laboratory and field isolates of Savannah exhibited a highly polymorphic repertoire both inter- and intra-isolates: sequences sharing the archetypical catalytic triad clustered into SAV1-SAV3 groups, whereas polymorphic sequences that, in general, exhibited unusual catalytic triad (variants) assigned to SAV4 or not assigned to any group. Congopain homologous genes from Forest and Kilifi isolates showed, respectively, moderate and limited diversity. In the phylogenetic tree based on congopain and homologues, Savannah was closer to Forest than to Kilifi. All T. congolense subgroup nested into a single clade, which together with the sister clade formed by homologues from Trypanosoma simiae and Trypanosoma godfreyi formed a clade supporting the subgenus Nannomonas. A single PCR targeting congopain sequences was developed for the diagnosis of T. congolense isolates of the three subgroups. Our findings demonstrated that congopain genes are valuable targets for the diagnosis, genotyping, and phylogenetic and taxonomic inferences among T. congolense isolates and other members of the subgenus Nannomonas.

Endosymbiosis in trypanosomatids: the genomic cooperation between bacterium and host in the synthesis of essential amino acids is heavily influenced by multiple horizontal gene transfers.

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.