Evolution of tubulin gene arrays in Trypanosomatid parasites: genomic restructuring in Leishmania

Genomics of Trypanosomatidae: Where We Stand and What Needs to Be Done?

Trypanosomatids are easy to cultivate and they are (in many cases) amenable to genetic manipulation. Genome sequencing has become a standard tool routinely used in the study of these flagellates. In this review, we summarize the current state of the field and our vision of what needs to be done in order to achieve a more comprehensive picture of trypanosomatid evolution. This will also help to illuminate the lineage-specific proteins and pathways, which can be used as potential targets in treating diseases caused by these parasites.

Targeted integration by homologous recombination enables in situ tagging and replacement of genes in the marine microeukaryote Diplonema papillatum

Diplonemids are a group of highly diverse and abundant marine microeukaryotes that belong to the phylum Euglenozoa and form a sister clade to the well-studied, mostly parasitic kinetoplastids. Very little is known about the biology of diplonemids, as few species have been formally described and just one, Diplonema papillatum, has been studied to a decent extent at the molecular level. Following up on our previous results showing stable but random integration of delivered extraneous DNA, we demonstrate here homologous recombination in D. papillatum. Targeting various constructs to the intended position in the nuclear genome was successful when 5' and 3' homologous regions longer than 1 kbp were used, achieving N-terminal tagging with mCherry and gene replacement of α- and β-tubulins. For more convenient genetic manipulation, we designed a modular plasmid, pDP002, which bears a protein-A tag and used it to generate and express a C-terminally tagged mitoribosomal protein. Lastly, we developed an improved transformation protocol for broader applicability across laboratories. Our robust methodology allows the replacement, integration as well as endogenous tagging of D. papillatum genes, thus opening the door to functional studies in this species and establishing a basic toolkit for reverse genetics of diplonemids in general.

Transcriptional and genomic parallels between the monoxenous parasite Herpetomonas muscarum and Leishmania

Trypanosomatid parasites are causative agents of important human and animal diseases such as sleeping sickness and leishmaniasis. Most trypanosomatids are transmitted to their mammalian hosts by insects, often belonging to Diptera (or true flies). These are called dixenous trypanosomatids since they infect two different hosts, in contrast to those that infect just insects (monoxenous). However, it is still unclear whether dixenous and monoxenous trypanosomatids interact similarly with their insect host, as fly-monoxenous trypanosomatid interaction systems are rarely reported and under-studied-despite being common in nature. Here we present the genome of monoxenous trypanosomatid Herpetomonas muscarum and discuss its transcriptome during in vitro culture and during infection of its natural insect host Drosophila melanogaster. The H. muscarum genome is broadly syntenic with that of human parasite Leishmania major. We also found strong similarities between the H. muscarum transcriptome during fruit fly infection, and those of Leishmania during sand fly infections. Overall this suggests Drosophila-Herpetomonas is a suitable model for less accessible insect-trypanosomatid host-parasite systems such as sand fly-Leishmania.

More than Microtubules: The Structure and Function of the Subpellicular Array in Trypanosomatids

The subpellicular microtubule array defines the wide range of cellular morphologies found in parasitic kinetoplastids (trypanosomatids). Morphological studies have characterized array organization, but little progress has been made towards identifying the molecular mechanisms that are responsible for array differentiation during the trypanosomatid life cycle, or the apparent stability and longevity of array microtubules. In this review, we outline what is known about the structure and biogenesis of the array, with emphasis on Trypanosoma brucei, Trypanosoma cruzi, and Leishmania, which cause life-threatening diseases in humans and livestock. We highlight unanswered questions about this remarkable cellular structure that merit new consideration in light of our recently improved understanding of how the 'tubulin code' influences microtubule dynamics to generate complex cellular structures.

Transcriptome, proteome and draft genome of Euglena gracilis

BACKGROUND

Photosynthetic euglenids are major contributors to fresh water ecosystems. Euglena gracilis in particular has noted metabolic flexibility, reflected by an ability to thrive in a range of harsh environments. E. gracilis has been a popular model organism and of considerable biotechnological interest, but the absence of a gene catalogue has hampered both basic research and translational efforts.

RESULTS

We report a detailed transcriptome and partial genome for E. gracilis Z1. The nuclear genome is estimated to be around 500 Mb in size, and the transcriptome encodes over 36,000 proteins and the genome possesses less than 1% coding sequence. Annotation of coding sequences indicates a highly sophisticated endomembrane system, RNA processing mechanisms and nuclear genome contributions from several photosynthetic lineages. Multiple gene families, including likely signal transduction components, have been massively expanded. Alterations in protein abundance are controlled post-transcriptionally between light and dark conditions, surprisingly similar to trypanosomatids.

CONCLUSIONS

Our data provide evidence that a range of photosynthetic eukaryotes contributed to the Euglena nuclear genome, evidence in support of the 'shopping bag' hypothesis for plastid acquisition. We also suggest that euglenids possess unique regulatory mechanisms for achieving extreme adaptability, through mechanisms of paralog expansion and gene acquisition.

Euglena gracilis Genome and Transcriptome: Organelles, Nuclear Genome Assembly Strategies and Initial Features

Euglena gracilis is a major component of the aquatic ecosystem and together with closely related species, is ubiquitous worldwide. Euglenoids are an important group of protists, possessing a secondarily acquired plastid and are relatives to the Kinetoplastidae, which themselves have global impact as disease agents. To understand the biology of E. gracilis, as well as to provide further insight into the evolution and origins of the Kinetoplastidae, we embarked on sequencing the nuclear genome; the plastid and mitochondrial genomes are already in the public domain. Earlier studies suggested an extensive nuclear DNA content, with likely a high degree of repetitive sequence, together with significant extrachromosomal elements. To produce a list of coding sequences we have combined transcriptome data from both published and new sources, as well as embarked on de novo sequencing using a combination of 454, Illumina paired end libraries and long PacBio reads. Preliminary analysis suggests a surprisingly large genome approaching 2 Gbp, with a highly fragmented architecture and extensive repeat composition. Over 80% of the RNAseq reads from E. gracilis maps to the assembled genome sequence, which is comparable with the well assembled genomes of T. brucei and T. cruzi. In order to achieve this level of assembly we employed multiple informatics pipelines, which are discussed here. Finally, as a preliminary view of the genome architecture, we discuss the tubulin and calmodulin genes, which highlight potential novel splicing mechanisms.

Six subgroups and extensive recent duplications characterize the evolution of the eukaryotic tubulin protein family

Tubulins belong to the most abundant proteins in eukaryotes providing the backbone for many cellular substructures like the mitotic and meiotic spindles, the intracellular cytoskeletal network, and the axonemes of cilia and flagella. Homologs have even been reported for archaea and bacteria. However, a taxonomically broad and whole-genome-based analysis of the tubulin protein family has never been performed, and thus, the number of subfamilies, their taxonomic distribution, and the exact grouping of the supposed archaeal and bacterial homologs are unknown. Here, we present the analysis of 3,524 tubulins from 504 species. The tubulins formed six major subfamilies, α to ζ. Species of all major kingdoms of the eukaryotes encode members of these subfamilies implying that they must have already been present in the last common eukaryotic ancestor. The proposed archaeal homologs grouped together with the bacterial TubZ proteins as sister clade to the FtsZ proteins indicating that tubulins are unique to eukaryotes. Most species contained α- and/or β-tubulin gene duplicates resulting from recent branch- and species-specific duplication events. This shows that tubulins cannot be used for constructing species phylogenies without resolving their ortholog–paralog relationships. The many gene duplicates and also the independent loss of the δ-, ε-, or ζ-tubulins, which have been shown to be part of the triplet microtubules in basal bodies, suggest that tubulins can functionally substitute each other.

A panel of Trypanosoma brucei strains tagged with blue and red-shifted luciferases for bioluminescent imaging in murine infection models

Background

Genetic engineering with luciferase reporter genes allows monitoring Trypanosoma brucei (T.b.) infections in mice by in vivo bioluminescence imaging (BLI). Until recently, luminescent T.b. models were based on Renilla luciferase (RLuc) activity. Our study aimed at evaluating red-shifted luciferases for in vivo BLI in a set of diverse T.b. strains of all three subspecies, including some recently isolated from human patients.

Methodology/Principal findings

We transfected T.b. brucei, T.b. rhodesiense and T.b. gambiense strains with either RLuc, click beetle red (CBR) or Photinus pyralis RE9 (PpyRE9) luciferase and characterised their in vitro luciferase activity, growth profile and drug sensitivity, and their potential for in vivo BLI. Compared to RLuc, the red-shifted luciferases, CBR and PpyRE9, allow tracking of T.b. brucei AnTaR 1 trypanosomes with higher details on tissue distribution, and PpyRE9 allows detection of the parasites with a sensitivity of at least one order of magnitude higher than CBR luciferase. With CBR-tagged T.b. gambiense LiTaR1, T.b. rhodesiense RUMPHI and T.b. gambiense 348 BT in an acute, subacute and chronic infection model respectively, we observed differences in parasite tropism for murine tissues during in vivo BLI. Ex vivo BLI on the brain confirmed central nervous system infection by all luminescent strains of T.b. brucei AnTaR 1, T.b. rhodesiense RUMPHI and T.b. gambiense 348 BT.

Conclusions/Significance

We established a genetically and phenotypically diverse collection of bioluminescent T.b. brucei, T.b. gambiense and T.b. rhodesiense strains, including drug resistant strains. For in vivo BLI monitoring of murine infections, we recommend trypanosome strains transfected with red-shifted luciferase reporter genes, such as CBR and PpyRE9. Red-shifted luciferases can be detected with a higher sensitivity in vivo and at the same time they improve the spatial resolution of the parasites in the entire body due to the better kinetics of their substrate D-luciferin.

Research on African trypanosomes heavily relies on rodent infection models. One way to reduce the number of laboratory rodents used in each experiment and effectively follow the progression of the infection in the same animals is to use genetically modified trypanosomes that allow monitoring of the infection over time with bioluminescence technology, without having to sacrifice the animals at multiple time points. In this study, we were able to establish a collection of bioluminescent strains of all three subspecies of Trypanosoma brucei (T.b.), including T.b. gambiense and T.b. rhodesiense that cause human African trypanosomiasis (HAT) or sleeping sickness. Making use of bioluminescence assays, we demonstrate the diversity of our collection in terms of in vitro and in vivo growth, drug sensitivity and in vivo parasite distribution, including central nervous system tropism. Growth characteristics and drug sensitivity are not affected by the genetic modification with luciferase reporter genes. Trypanosome strains transfected with red-shifted luciferase reporter genes have several advantages compared to the corresponding blue luciferase modified strains. Red light is less absorbed in the blood than blue light, which should lead to higher sensitivity of detection. Furthermore, the substrates that drive the light reaction are better distributed through the body for the red luciferase than for the blue luciferase, which greatly improves spatial resolution of the infection.

New insights about cross-reactive epitopes of six trypanosomatid genera revealed that Crithidia and Leptomonas have antigenic similarity to L. (L.) chagasi

We investigated whether ELISA using crude antigens from insect and plant trypanosomatids, which are non-pathogenic and easily cultivated in large scale, has the same positivity data as Leishmania (Leishmania) chagasi, the etiological agent of human visceral leishmaniasis (VL) or canine leishmaniasis (CanL), or as Trypanosoma cruzi, the etiological agent of Chagas disease (CD). The antigens from Crithidia fasciculata, Crithidia luciliae, and Leptomonas seymouri showed 100% cross-reactivity with VL and CanL samples, with no statistically titers differences from L. (L.) chagasi, however, 34% (17/50) of VL samples revealed higher titers using the insect trypanosomatids than the homologous antigen. On the other hand, antigens from Strigomonas culicis, Angomonas deanei, and Phytomonas serpens showed low cross-reactivity with VL and CanL samples. The sera from patients with American tegumentary leishmaniasis showed low levels of cross-reactivity with all trypanosomatids investigated, even with L. (L) chagasi, without titers dissimilarity among them. These parasites were also worthless as antigen source for detection of CD cases, which required homologous antigens to reach 100% positivity. This study showed, by ELISA, that crude extract of Crithidia and Leptomonas have epitopes similar to L. (L.) chagasi, which supports the idea of using them as antigens source for the serodiagnosis of visceral leishmaniasis.

Alpha tubulin genes from Leishmania braziliensis: genomic organization, gene structure and insights on their expression

Background

Alpha tubulin is a fundamental component of the cytoskeleton which is responsible for cell shape and is involved in cell division, ciliary and flagellar motility and intracellular transport. Alpha tubulin gene expression varies according to the morphological changes suffered by Leishmania in its life cycle. However, the objective of studying the mechanisms responsible for the differential expression has resulted to be a difficult task due to the complex genome organization of tubulin genes and to the non-conventional mechanisms of gene regulation operating in Leishmania.

Results

We started this work by analyzing the genomic organization of α-tubulin genes in the Leishmania braziliensis genome database. The genomic organization of L. braziliensis α-tubulin genes differs from that existing in the L. major and L. infantum genomes. Two loci containing α-tubulin genes were found in the chromosomes 13 and 29, even though the existence of sequence gaps does not allow knowing the exact number of genes at each locus. Southern blot assays showed that α-tubulin locus at chromosome 13 contains at least 8 gene copies, which are tandemly organized with a 2.08-kb repetition unit; the locus at chromosome 29 seems to contain a sole α-tubulin gene. In addition, it was found that L. braziliensis α-tubulin locus at chromosome 13 contains two types of α-tubulin genes differing in their 3′ UTR, each one presumably containing different regulatory motifs. It was also determined that the mRNA expression levels of these genes are controlled by post-transcriptional mechanisms tightly linked to the growth temperature. Moreover, the decrease in the α-tubulin mRNA abundance observed when promastigotes were cultured at 35°C was accompanied by parasite morphology alterations, similar to that occurring during the promastigote to amastigote differentiation.

Conclusions

Information found in the genome databases indicates that α-tubulin genes have been reorganized in a drastic manner along Leishmania speciation. In the L. braziliensis genome database, two loci containing α-tubulin sequences were found, but only the locus at chromosome 13 contains the prototypic α-tubulin genes, which are repeated in a head-to-tail manner. Also, we determined that the levels of α-tubulin mRNAs are down-regulated drastically in response to heat shock by a post-transcriptional mechanism which is dependent upon active protein synthesis.

Genome organization is a major component of gene expression control in response to stress and during the cell division cycle in trypanosomes

The trypanosome genome is characterized by RNA polymerase II-driven polycistronic transcription of protein-coding genes. Ten to hundreds of genes are co-transcribed from a single promoter; thus, selective regulation of individual genes via initiation is impossible. However, selective responses to external stimuli occur and post-transcriptional mechanisms are thought to account for all temporal gene expression patterns. We show that genes encoding mRNAs that are differentially regulated during the heat-shock response are selectively positioned in polycistronic transcription units; downregulated genes are close to transcription initiation sites and upregulated genes are distant. We demonstrate that the position of a reporter gene within a transcription unit is sufficient to reproduce this effect. Analysis of gene ontology annotations reveals that positional bias is not restricted to stress-response genes and that there is a genome-wide organization based on proximity to transcription initiation sites. Furthermore, we show that the relative abundance of mRNAs at different time points in the cell division cycle is dependent on the location of the corresponding genes to transcription initiation sites. This work provides evidence that the genome in trypanosomes is organized to facilitate co-coordinated temporal control of gene expression in the absence of selective promoters.

In Vitro Trypanocidal Activity of Antibodies to Bacterially Expressed Trypanosoma brucei Tubulin

BACKGROUND

There are only four drugs for treating African trypanosomiasis, a devastating disease in sub-Saharan Africa. With slow discovery of better drugs, vaccination is viewed as the best method of control. We previously showed that antibodies to native Trypanosoma brucei brucei tubulin inhibit the growth of trypanosomes in culture. Here, we aimed to determine the effect of antibodies to bacterially expressed trypanosome tubulin on T. brucei brucei growth.

METHODS

T. brucei brucei alpha and beta tubulin genes were individually expressed in Escherichia coli under the tryptophan promoter. Monoclonal tubulin antibodies reacted specifically with the expressed tubulins with no cross-reaction with the opposite tubulin. Rabbits were immunized with 450µg each of the concentrated recombinant tubulin, and production of antibodies assessed by ELISA and Western blotting. The effect of polyclonal antibodies on trypanosome growth was determined by culturing bloodstream T. brucei brucei in up to 25% of antisera.

RESULTS

Low antisera dilutions (25%) from the immunized rabbits inhibited trypanosome growth. The most cytotoxic antisera were from one rabbit immunized with a mixture of both alpha and beta tubulins. However, the result was not reproduced in other rabbits and there was no apparent effect on growth at higher antisera dilutions.

CONCLUSION

Antibodies to bacterially expressed trypanosome tubulin are not effective at killing cultured bloodstream trypanosomes.

Chromosome and gene copy number variation allow major structural change between species and strains of Leishmania

Leishmania parasites cause a spectrum of clinical pathology in humans ranging from disfiguring cutaneous lesions to fatal visceral leishmaniasis. We have generated a reference genome for Leishmania mexicana and refined the reference genomes for Leishmania major, Leishmania infantum, and Leishmania braziliensis. This has allowed the identification of a remarkably low number of genes or paralog groups (2, 14, 19, and 67, respectively) unique to one species. These were found to be conserved in additional isolates of the same species. We have predicted allelic variation and find that in these isolates, L. major and L. infantum have a surprisingly low number of predicted heterozygous SNPs compared with L. braziliensis and L. mexicana. We used short read coverage to infer ploidy and gene copy numbers, identifying large copy number variations between species, with 200 tandem gene arrays in L. major and 132 in L. mexicana. Chromosome copy number also varied significantly between species, with nine supernumerary chromosomes in L. infantum, four in L. mexicana, two in L. braziliensis, and one in L. major. A significant bias against gene arrays on supernumerary chromosomes was shown to exist, indicating that duplication events occur more frequently on disomic chromosomes. Taken together, our data demonstrate that there is little variation in unique gene content across Leishmania species, but large-scale genetic heterogeneity can result through gene amplification on disomic chromosomes and variation in chromosome number. Increased gene copy number due to chromosome amplification may contribute to alterations in gene expression in response to environmental conditions in the host, providing a genetic basis for disease tropism.

Proteome mapping of overexpressed membrane-enriched and cytosolic proteins in sodium antimony gluconate (SAG) resistant clinical isolate of Leishmania donovani

AIMS

This study aimed to identify differentially overexpressed membrane-enriched as well as cytosolic proteins in SAG sensitive and resistant clinical strains of L. donovani isolated from VL patients which are involved in the drug resistance mechanism.

METHODS

The proteins in the membrane-enriched as well as cytosolic fractions of drug-sensitive as well as drug-resistant clinical isolates were separated using two-dimensional gel electrophoresis and overexpressed identified protein spots of interest were excised and analysed using MALDI-TOF/TOF.

RESULTS

Six out of 12 overexpressed proteins were identified in the membrane-enriched fraction of the SAG resistant strain of L. donovani whereas 14 out of 18 spots were identified in the cytosolic fraction as compared with the SAG sensitive strain. The major proteins in the membrane-enriched fraction were ABC transporter, HSP-83, GPI protein transamidase, cysteine-leucine rich protein and 60S ribosomal protein L23a whereas in the cytosolic fraction proliferative cell nuclear antigen (PCNA), proteasome alpha 5 subunit, carboxypeptidase, HSP-70, enolase, fructose-1,6-bisphosphate aldolase, tubulin-beta chain have been identified. Most of these proteins have been reported as potential drug targets, except 60S ribosomal protein L23a and PCNA which have not been reported to date for their possible involvement in drug resistance against VL.

CONCLUSION

This study for the first time provided a cumulative proteomic analysis of proteins overexpressed in drug resistant clinical isolates of L. donovani indicating their possible role in antimony resistance of the parasite. Identified proteins provide a vast field to be exploited for novel treatment strategies against VL such as cloning and overexpression of these targets to produce recombinant therapeutic/prophylactic proteins.

BtubA-BtubB heterodimer is an essential intermediate in protofilament assembly

Background

BtubA and BtubB are two tubulin-like genes found in the bacterium Prosthecobacter. Our work and a previous crystal structure suggest that BtubB corresponds to α−tubulin and BtubA to β−tubulin. A 1∶1 mixture of the two proteins assembles into tubulin-like protofilaments, which further aggregate into pairs and bundles. The proteins also form a BtubA/B heterodimer, which appears to be a repeating subunit in the protofilament.

Methodology/Principal Findings

We have designed point mutations to disrupt the longitudinal interfaces bonding subunits into protofilaments. The mutants are in two classes, within dimers and between dimers. We have characterized one mutant of each class for BtubA and BtubB. When mixed 1∶1 with a wild type partner, none of the mutants were capable of assembly. An excess of between-dimer mutants could depolymerize preformed wild type polymers, while within-dimer mutants had no activity.

Conclusions

An essential first step in assembly of BtubA + BtubB is formation of a heterodimer. An excess of between-dimer mutants depolymerize wild type BtubA/B by sequestering the partner wild type subunit into inactive dimers. Within-dimer mutants cannot form dimers and have no activity.

Gene duplication, transfer, and evolution in the chloroplast genome

In addition to the nuclear genome, organisms have organelle genomes. Most of the DNA present in eukaryotic organisms is located in the cell nucleus. Chloroplasts have independent genomes which are inherited from the mother. Duplicated genes are common in the genomes of all organisms. It is believed that gene duplication is the most important step for the origin of genetic variation, leading to the creation of new genes and new gene functions. Despite the fact that extensive gene duplications are rare among the chloroplast genome, gene duplication in the chloroplast genome is an essential source of new genetic functions and a mechanism of neo-evolution. The events of gene transfer between the chloroplast genome and nuclear genome via duplication and subsequent recombination are important processes in evolution. The duplicated gene or genome in the nucleus has been the subject of several recent reviews. In this review, we will briefly summarize gene duplication and evolution in the chloroplast genome. Also, we will provide an overview of gene transfer events between chloroplast and nuclear genomes.

Insights into the genome sequence of a free-living Kinetoplastid: Bodo saltans (Kinetoplastida: Euglenozoa)

Background

Bodo saltans is a free-living kinetoplastid and among the closest relatives of the trypanosomatid parasites, which cause such human diseases as African sleeping sickness, leishmaniasis and Chagas disease. A B. saltans genome sequence will provide a free-living comparison with parasitic genomes necessary for comparative analyses of existing and future trypanosomatid genomic resources. Various coding regions were sequenced to provide a preliminary insight into the bodonid genome sequence, relative to trypanosomatid sequences.

Results

0.4 Mbp of B. saltans genome was sequenced from 12 distinct regions and contained 178 coding sequences. As in trypanosomatids, introns were absent and %GC was elevated in coding regions, greatly assisting in gene finding. In the regions studied, roughly 60% of all genes had homologs in trypanosomatids, while 28% were Bodo-specific. Intergenic sequences were typically short, resulting in higher gene density than in trypanosomatids. Although synteny was typically conserved for those genes with trypanosomatid homologs, strict colinearity was rarely observed because gene order was regularly disrupted by Bodo-specific genes.

Conclusion

The B. saltans genome contains both sequences homologous to trypanosomatids and sequences never seen before. Structural similarities suggest that its assembly should be solvable, and, although de novo assembly will be necessary, existing trypanosomatid projects will provide some guide to annotation. A complete genome sequence will provide an effective ancestral model for understanding the shared and derived features of known trypanosomatid genomes, but it will also identify those kinetoplastid genome features lost during the evolution of parasitism.

Evolutionary analysis of synteny and gene fusion for pyrimidine biosynthetic enzymes in Euglenozoa: an extraordinary gap between kinetoplastids and diplonemids

A unique feature of the genome architecture in the parasitic trypanosomatid protists is large-scale synteny. We addressed the evolutionary trait of synteny in the eukaryotic group, Euglenozoa, which consists of euglenoids (earliest branching), diplonemids, and kinetoplastids (trypanosomatids and bodonids). Synteny of the pyrimidine biosynthetic (pyr) gene cluster, which constitutes part of a large syntenic cluster in trypanosomatids and includes four separate genes (pyr1-pyr4) and one fused gene (pyr6/pyr5 fusion), was conserved in the bodonid, Parabodo caudatus. In the diplonemid, Diplonema papillatum, we identified pyr4 and pyr6 genes. Phylogenetic analyses of pyr4 and pyr6 showed the separate origin of each in kinetoplastids and euglenoids/diplonemids and suggested that kinetoplastids have acquired these genes via lateral gene transfer (LGT). Because replacement of genes by non-orthologs within the syntenic cluster is highly unlikely, we concluded that, after separation of the line leading to diplonemids, the syntenic pyr gene cluster was established in the common ancestor of kinetoplastids, preceded by their acquisition via LGT. Notably, we found that diplonemid pyr6 is a stand-alone gene, inconsistent with both euglenoid pyr5/pyr6 and kinetoplastid pyr6/pyr5 fusions. Our findings provide insights into the evolutionary gaps within Euglenozoa and the evolutionary trait of rearrangement of gene fusion in this lineage.