Gene synteny and evolution of genome architecture in trypanosomatids

Complete assembly, annotation of virulence genes and CRISPR editing of the genome of Leishmania amazonensis PH8 strain

We report the sequencing and assembly of the PH8 strain of Leishmania amazonensis one of the etiological agents of leishmaniasis. After combining data from long Pacbio reads, short Illumina reads and synteny with the Leishmania mexicana genome, the sequence of 34 chromosomes with 8317 annotated genes was generated. Multigene families encoding three virulence factors, A2, amastins and the GP63 metalloproteases, were identified and compared to their annotation in other Leishmania species. As they have been recently recognized as virulence factors essential for disease establishment and progression of the infection, we also identified 14 genes encoding proteins involved in parasite iron and heme metabolism and compared to genes from other Trypanosomatids. To follow these studies with a genetic approach to address the role of virulence factors, we tested two CRISPR-Cas9 protocols to generate L. amazonensis knockout cell lines, using the Miltefosine transporter gene as a proof of concept.

A resource for improved predictions of Trypanosoma and Leishmania protein three-dimensional structure

AlphaFold2 and RoseTTAfold represent a transformative advance for predicting protein structure. They are able to make very high-quality predictions given a high-quality alignment of the protein sequence with related proteins. These predictions are now readily available via the AlphaFold database of predicted structures and AlphaFold or RoseTTAfold Colaboratory notebooks for custom predictions. However, predictions for some species tend to be lower confidence than model organisms. Problematic species include Trypanosoma cruzi and Leishmania infantum: important unicellular eukaryotic human parasites in an early-branching eukaryotic lineage. The cause appears to be due to poor sampling of this branch of life (Discoba) in the protein sequences databases used for the AlphaFold database and ColabFold. Here, by comprehensively gathering openly available protein sequence data for Discoba species, significant improvements to AlphaFold2 protein structure prediction over the AlphaFold database and ColabFold are demonstrated. This is made available as an easy-to-use tool for the parasitology community in the form of Colaboratory notebooks for generating multiple sequence alignments and AlphaFold2 predictions of protein structure for Trypanosoma, Leishmania and related species.

Rapid protein evolution, organellar reductions, and invasive intronic elements in the marine aerobic parasite dinoflagellate Amoebophrya spp

BACKGROUND

Dinoflagellates are aquatic protists particularly widespread in the oceans worldwide. Some are responsible for toxic blooms while others live in symbiotic relationships, either as mutualistic symbionts in corals or as parasites infecting other protists and animals. Dinoflagellates harbor atypically large genomes (~ 3 to 250 Gb), with gene organization and gene expression patterns very different from closely related apicomplexan parasites. Here we sequenced and analyzed the genomes of two early-diverging and co-occurring parasitic dinoflagellate Amoebophrya strains, to shed light on the emergence of such atypical genomic features, dinoflagellate evolution, and host specialization.

RESULTS

We sequenced, assembled, and annotated high-quality genomes for two Amoebophrya strains (A25 and A120), using a combination of Illumina paired-end short-read and Oxford Nanopore Technology (ONT) MinION long-read sequencing approaches. We found a small number of transposable elements, along with short introns and intergenic regions, and a limited number of gene families, together contribute to the compactness of the Amoebophrya genomes, a feature potentially linked with parasitism. While the majority of Amoebophrya proteins (63.7% of A25 and 59.3% of A120) had no functional assignment, we found many orthologs shared with Dinophyceae. Our analyses revealed a strong tendency for genes encoded by unidirectional clusters and high levels of synteny conservation between the two genomes despite low interspecific protein sequence similarity, suggesting rapid protein evolution. Most strikingly, we identified a large portion of non-canonical introns, including repeated introns, displaying a broad variability of associated splicing motifs never observed among eukaryotes. Those introner elements appear to have the capacity to spread over their respective genomes in a manner similar to transposable elements. Finally, we confirmed the reduction of organelles observed in Amoebophrya spp., i.e., loss of the plastid, potential loss of a mitochondrial genome and functions.

CONCLUSION

These results expand the range of atypical genome features found in basal dinoflagellates and raise questions regarding speciation and the evolutionary mechanisms at play while parastitism was selected for in this particular unicellular lineage.

Identification of Novel Interspersed DNA Repetitive Elements in the Trypanosoma cruzi Genome Associated with the 3'UTRs of Surface Multigenic Families

Trypanosoma cruzi is the etiological agent of Chagas disease, which affects millions of people in Latin America. No transcriptional control of gene expression has been demonstrated in this organism, and 50% of its genome consists of repetitive elements and members of multigenic families. In this study, we applied a novel bioinformatics approach to predict new repetitive elements in the genome sequence of T. cruzi. A new repetitive sequence measuring 241 nt was identified and found to be interspersed along the genome sequence from strains of different DTUs. This new repeat was mostly on intergenic regions, and upstream and downstream regions of the 241 nt repeat were enriched in surface protein genes. RNAseq analysis revealed that the repeat was part of processed mRNAs and was predominantly found in the 3′ untranslated regions (UTRs) of genes of multigenic families encoding surface proteins. Moreover, we detected a correlation between the presence of the repeat in the 3′UTR of multigenic family genes and the level of differential expression of these genes when comparing epimastigote and trypomastigote transcriptomes. These data suggest that this sequence plays a role in the posttranscriptional regulation of the expression of multigenic families.

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.

A Trypanosoma cruzi Genome Tandem Repetitive Satellite DNA Sequence as a Molecular Marker for a LAMP Assay for Diagnosing Chagas' Disease

Chagas' disease is a neglected tropical disease caused by Trypanosoma cruzi which is endemic throughout Latin America and is spread by worldwide migration. Diagnosis is currently limited to serological and molecular techniques having variations regarding their sensitivity and specificity. This work was aimed at developing a new sensitive, applicable, and cost-effective molecular diagnosis technique for loop-mediated isothermal amplification-based detection of T. cruzi (Tc-LAMP). The results led to determining a highly homologous satellite repeat region (231 bp) among parasite strains as a molecular marker for diagnosing the disease. Tc-LAMP was performed correctly for detecting parasite DNA (5 fg for the CL Brener strain and 50 fg for the DM28, TcVI, and TcI strains). Assay results proved negative for DNA from 16 helminth species and 7 protozoa, including Leishmania spp. Tc-LAMP based on the highly repeated T. cruzi satellite region is thus proposed as an important alternative for diagnosing T. cruzi infection, overcoming other methods' limitations such as their analytic capability, speed, and requiring specialized equipment or highly trained personnel. Tc-LAMP could be easily adapted for point-of-care testing in areas having limited resources.

Evolutionary Insight into the Trypanosomatidae Using Alignment-Free Phylogenomics of the Kinetoplast

Advancements in next-generation sequencing techniques have led to a substantial increase in the genomic information available for analyses in evolutionary biology. As such, this data requires the exponential growth in bioinformatic methods and expertise required to understand such vast quantities of genomic data. Alignment-free phylogenomics offer an alternative approach for large-scale analyses that may have the potential to address these challenges. The evolutionary relationships between various species within the trypanosomatid family, specifically members belonging to the genera Leishmania and Trypanosoma have been extensively studies over the last 30 years. However, there is a need for a more exhaustive analysis of the Trypanosomatidae, summarising the evolutionary patterns amongst the entire family of these important protists. The mitochondrial DNA of the trypanosomatids, better known as the kinetoplast, represents a valuable taxonomic marker given its unique presence across all kinetoplastid protozoans. The aim of this study was to validate the reliability and robustness of alignment-free approaches for phylogenomic analyses and its applicability to reconstruct the evolutionary relationships between the trypanosomatid family. In the present study, alignment-free analyses demonstrated the strength of these methods, particularly when dealing with large datasets compared to the traditional phylogenetic approaches. We present a maxicircle genome phylogeny of 46 species spanning the trypanosomatid family, demonstrating the superiority of the maxicircle for the analysis and taxonomic resolution of the Trypanosomatidae.

The Tritryps Comparative Repeatome: Insights on Repetitive Element Evolution in Trypanosomatid Pathogens

The major human pathogens Trypanosoma cruzi, Trypanosoma brucei, and Leishmania major are collectively known as the Tritryps. The initial comparative analysis of their genomes has uncovered that Tritryps share a great number of genes, but repetitive DNA seems to be extremely variable between them. However, the in-depth characterization of repetitive DNA in these pathogens has been in part neglected, mainly due to the well-known technical challenges of studying repetitive sequences from de novo assemblies using short reads. Here, we compared the repetitive DNA repertories between the Tritryps genomes using genome-wide, low-coverage Illumina sequencing coupled to RepeatExplorer analysis. Our work demonstrates that this extensively implemented approach for studying higher eukaryote repeatomes is also useful for protozoan parasites like trypanosomatids, as we recovered previously observed differences in the presence and amount of repetitive DNA families. Additionally, our estimations of repetitive DNA abundance were comparable to those obtained from enhanced-quality assemblies using longer reads. Importantly, our methodology allowed us to describe a previously undescribed transposable element in Leishmania major (TATE element), highlighting its potential to accurately recover distinctive features from poorly characterized repeatomes. Together, our results support the application of this low-cost, low-coverage sequencing approach for the extensive characterization of repetitive DNA evolutionary dynamics in trypanosomatid and other protozoan genomes.

Basalin is an evolutionarily unconstrained protein revealed via a conserved role in flagellum basal plate function

Most motile flagella have an axoneme that contains nine outer microtubule doublets and a central pair (CP) of microtubules. The CP coordinates the flagellar beat and defects in CP projections are associated with motility defects and human disease. The CP nucleate near a ‘basal plate’ at the distal end of the transition zone (TZ). Here, we show that the trypanosome TZ protein ‘basalin’ is essential for building the basal plate, and its loss is associated with CP nucleation defects, inefficient recruitment of CP assembly factors to the TZ, and flagellum paralysis. Guided by synteny, we identified a highly divergent basalin ortholog in the related Leishmania species. Basalins are predicted to be highly unstructured, suggesting they may act as ‘hubs’ facilitating many protein-protein interactions. This raises the general concept that proteins involved in cytoskeletal functions and appearing organism-specific, may have highly divergent and cryptic orthologs in other 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.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5112-0) contains supplementary material, which is available to authorized users.

Recent advances in trypanosomatid research: genome organization, expression, metabolism, taxonomy and evolution

Unicellular flagellates of the family Trypanosomatidae are obligatory parasites of invertebrates, vertebrates and plants. Dixenous species are aetiological agents of a number of diseases in humans, domestic animals and plants. Their monoxenous relatives are restricted to insects. Because of the high biological diversity, adaptability to dramatically different environmental conditions, and omnipresence, these protists have major impact on all biotic communities that still needs to be fully elucidated. In addition, as these organisms represent a highly divergent evolutionary lineage, they are strikingly different from the common 'model system' eukaryotes, such as some mammals, plants or fungi. A number of excellent reviews, published over the past decade, were dedicated to specialized topics from the areas of trypanosomatid molecular and cell biology, biochemistry, host-parasite relationships or other aspects of these fascinating organisms. However, there is a need for a more comprehensive review that summarizing recent advances in the studies of trypanosomatids in the last 30 years, a task, which we tried to accomplish with the current paper.

Biology of Trypanosoma cruzi Retrotransposons: From an Enzymatic to a Structural Point of View

INTRODUCTION

An important portion of the Trypanosoma cruzi genome is composed of mobile genetic elements, which are interspersed with genes on all chromosomes. The L1Tc non-LTR retrotransposon and its truncated version NARTc are the most highly represented and best studied of these elements. L1Tc is actively transcribed in all three forms of the Trypanosoma parasite and encodes the proteins that enable it to autonomously mobilize. This mini review discusses the enzymatic properties of L1Tc that enable its mobilization and possibly the mobilization of other non-autonomous retrotransposons in Trypanosoma. We also briefly review the Hepatitis Delta Virus-like autocatalytic and 2A self-cleaving viral-like sequences contained in L1Tc that regulate post-transcriptional properties such as relative protein abundance and mRNA stability. Special emphasis is placed on the Pr77 dual system, which is based on the RNA pol II-dependent internal promoter of L1Tc and NARTc and the HDV-like ribozyme activity encoded by the first 77 nucleotides of the element's DNA and RNA. The high degree of conservation of the Pr77 sequence, referred to as the "Pr77-hallmark", among different trypanosomatid retroelements suggests that these mobile elements are responsible for the distribution of regulatory sequences within the genome they inhabit.

CONCLUSION

We also discuss how the involvement of L1Tc and NARTc in the gene regulatory processes of these parasites could justify their domestication and long-term coexistence in these ancient organisms.

PhyloChromoMap, a Tool for Mapping Phylogenomic History along Chromosomes, Reveals the Dynamic Nature of Karyotype Evolution in Plasmodium falciparum

The genome of Plasmodium falciparum, the causative agent of malaria in Africa, has been extensively studied since it was first fully sequenced in 2002. However, many open questions remain, including understanding the chromosomal context of molecular evolutionary changes (e.g., relationship between chromosome map and phylogenetic conservation, patterns of gene duplication, and patterns of selection). Here, we present PhyloChromoMap, a method that generates a phylogenomic map of chromosomes from a custom-built bioinformatics pipeline. Using P. falciparum 3D7 as a model, we analyze 2,116 genes with homologs in up to 941 diverse eukaryotic, bacterial and archaeal lineages. We estimate the level of conservation along chromosomes based on conservation across clades, and identify “young” regions (i.e., those with recent or fast evolving genes) that are enriched in subtelomeric regions as compared with internal regions. We also demonstrate that patterns of molecular evolution for paralogous genes differ significantly depending on their location as younger paralogs tend to be found in subtelomeric regions whereas older paralogs are enriched in internal regions. Combining these observations with analyses of synteny, we demonstrate that subtelomeric regions are actively shuffled among chromosome ends, which is consistent with the hypothesis that these regions are prone to ectopic recombination. We also assess patterns of selection by comparing dN/dS ratios of gene family members in subtelomeric versus internal regions, and we include the important antigenic gene family var. These analyses illustrate the highly dynamic nature of the karyotype of P. falciparum, and provide a method for exploring genome dynamics in other lineages.

Identification of the centromeres of Leishmania major: revealing the hidden pieces

Leishmania affects millions of people worldwide. Its genome undergoes constitutive mosaic aneuploidy, a type of genomic plasticity that may serve as an adaptive strategy to survive distinct host environments. We previously found high rates of asymmetric chromosome allotments during mitosis that lead to the generation of such ploidy. However, the underlying molecular events remain elusive. Centromeres and kinetochores most likely play a key role in this process, yet their identification has failed using classical methods. Our analysis of the unconventional kinetochore complex recently discovered in Trypanosoma brucei (KKTs) leads to the identification of a Leishmania KKT gene candidate (LmKKT1). The GFP-tagged LmKKT1 displays "kinetochore-like" dynamics of intranuclear localization throughout the cell cycle. By ChIP-Seq assay, one major peak per chromosome is revealed, covering a region of 4 ±2 kb. We find two largely conserved motifs mapping to 14 of 36 chromosomes while a higher density of retroposons are observed in 27 of 36 centromeres. The identification of centromeres and of a kinetochore component of Leishmania chromosomes opens avenues to explore their role in mosaic aneuploidy.

Long-read sequencing improves assembly of Trichinella genomes 10-fold, revealing substantial synteny between lineages diverged over 7 million years

Genome assemblies can form the basis of comparative analyses fostering insight into the evolutionary genetics of a parasite's pathogenicity, host-pathogen interactions, environmental constraints and invasion biology; however, the length and complexity of many parasite genomes has hampered the development of well-resolved assemblies. In order to improve Trichinella genome assemblies, the genome of the sylvatic encapsulated species Trichinella murrelli was sequenced using third-generation, long-read technology and, using syntenic comparisons, scaffolded to a reference genome assembly of Trichinella spiralis, markedly improving both. A high-quality draft assembly for T. murrelli was achieved that totalled 63·2 Mbp, half of which was condensed into 26 contigs each longer than 571 000 bp. When compared with previous assemblies for parasites in the genus, ours required 10-fold fewer contigs, which were five times longer, on average. Better assembly across repetitive regions also enabled resolution of 8 Mbp of previously indeterminate sequence. Furthermore, syntenic comparisons identified widespread scaffold misassemblies in the T. spiralis reference genome. The two new assemblies, organized for the first time into three chromosomal scaffolds, will be valuable resources for future studies linking phenotypic traits within each species to their underlying genetic bases.

Genomic cartography and proposal of nomenclature for the repeated, interspersed elements of the Leishmania major SIDER2 family and identification of SIDER2-containing transcripts

The genomes of most eukaryotic organisms contain a large number of transposable elements that are able to move from one genomic site to another either by transferring of DNA mobile elements (transposons) or transpose via reverse transcription of an RNA intermediate (retroposons). An exception to this rule is found in protists of the subgenus Leishmania, in which active retroposons degenerated after a flourishing era, leaving only retroposon remains; these have been classified into two families: SIDER1 and SIDER2. In this work, we have re-examined the elements belonging to the family SIDER2 present in the genome of Leishmania major with the aim of providing a nomenclature that will facilitate a future reference to particular elements. According to sequence conservation, the 1100 SIDER2 elements have been grouped into subfamilies, and the inferred taxonomic relationships have also been incorporated into the nomenclature. Additionally, we are providing detailed data regarding the genomic distribution of these elements and their association with specific transcripts, based on the recently established transcriptome for L. major. Thus, the presented data can help to study and better understand the roles played by these degenerated retroposons in both regulation of gene expression and genome plasticity.

Leishmania DNA Replication Timing: A Stochastic Event?

For eukaryotic genomes, DNA synthesis initiates at multiple discrete regions known as replication origins in a dynamic yet regulated manner to ensure genomic stability. Two recent studies using different approaches reveal few Leishmania origins and that origin firing may proceed in a mainly stochastic manner.

The Trypanosomatid Pr77-hallmark contains a downstream core promoter element essential for transcription activity of the Trypanosoma cruzi L1Tc retrotransposon

Background

Trypanosomatid genomes are highly colonized by non-LTR retroelements that make up to 5 % of the nuclear genome. These elements are mainly accumulated in the strand switch regions (SSRs) where polycistronic transcription is initiated and have a 77 nt-long sequence - Pr77 - at their 5′ ends. L1Tc is the best represented retrotransposon in the Trypanosoma cruzi genome and is a potentially functional autonomous element that encodes its own retrotransposition machinery. The Pr77 of the T. cruzi L1Tc element activates gene transcription via RNA polymerase II, generating abundant, unspliced transcripts which are translated.

Results

The present manuscript describes the identification of a downstream core promoter element (DPE) in the L1Tc Pr77 sequence. Just four nucleotides long (CGTG), it covers in Pr77 positions +25 to +28 of the described L1Tc transcription start site. The Pr77-DPE motif is conserved in terms of sequence composition and position in the Pr77 of most trypanosomatid non-LTR retrotransposons, independent of the coding or non-coding capacity of these retroelements. Transcription assays in T. cruzi stable transfectants with vector containing point mutations at 17 locations of the Pr77 nucleotide sequence evidence that the DPE motif is essential for the promoter function of Pr77. Furthermore, the obtained data show that other nucleotides also contributed to the promoter function of Pr77. In addition, the presented results indicate that parasite nuclear proteins specifically bind to different regions of the Pr77 sequence although the strongest binding is to the DPE motif. Moreover, it is shown that the DPE sense single-stranded sequence is being required in DNA-protein recognition of nuclear factors.

Conclusions

The Pr77 sequence present in most of non-LTR retrotransposons of trypanosomatids contains a downstream core promoter element (DPE) which is conserved in terms of nucleotide composition and location. The Pr77-DPE motif is essential for the transcriptional activity of Pr77 although other nucleotides are also involved. DPE has a high affinity binding for nuclear proteins in T. cruzi. The wide retroelement-mediated distribution of Pr77 suggests that it may represent an important tool for regulating gene expression in trypanosomatids.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2427-6) contains supplementary material, which is available to authorized users.

Molecular Characterization of a Novel Family of Trypanosoma cruzi Surface Membrane Proteins (TcSMP) Involved in Mammalian Host Cell Invasion

Background

The surface coat of Trypanosoma cruzi is predominantly composed of glycosylphosphatidylinositol-anchored proteins, which have been extensively characterized. However, very little is known about less abundant surface proteins and their role in host-parasite interactions.

Methodology/ Principal Findings

Here, we described a novel family of T. cruzi surface membrane proteins (TcSMP), which are conserved among different T. cruzi lineages and have orthologs in other Trypanosoma species. TcSMP genes are densely clustered within the genome, suggesting that they could have originated by tandem gene duplication. Several lines of evidence indicate that TcSMP is a membrane-spanning protein located at the cellular surface and is released into the extracellular milieu. TcSMP exhibited the key elements typical of surface proteins (N-terminal signal peptide or signal anchor) and a C-terminal hydrophobic sequence predicted to be a trans-membrane domain. Immunofluorescence of live parasites showed that anti-TcSMP antibodies clearly labeled the surface of all T. cruzi developmental forms. TcSMP peptides previously found in a membrane-enriched fraction were identified by proteomic analysis in membrane vesicles as well as in soluble forms in the T. cruzi secretome. TcSMP proteins were also located intracellularly likely associated with membrane-bound structures. We demonstrated that TcSMP proteins were capable of inhibiting metacyclic trypomastigote entry into host cells. TcSMP bound to mammalian cells and triggered Ca²⁺ signaling and lysosome exocytosis, events that are required for parasitophorous vacuole biogenesis. The effects of TcSMP were of lower magnitude compared to gp82, the major adhesion protein of metacyclic trypomastigotes, suggesting that TcSMP may play an auxiliary role in host cell invasion.

Conclusion/Significance

We hypothesized that the productive interaction of T. cruzi with host cells that effectively results in internalization may depend on diverse adhesion molecules. In the metacyclic forms, the signaling induced by TcSMP may be additive to that triggered by the major surface molecule gp82, further increasing the host cell responses required for infection.

Trypanosoma cruzi is the etiologic agent of Chagas’ disease, which infects 6–7 million people worldwide, mostly in Latin America. Currently, there are no vaccines available, and the drugs used for treatment are toxic and are not fully effective. To infect mammalian hosts, T. cruzi relies on the ability to invade host cells, replicate intracellularly and spread the infection in different organs of the mammalian host. Knowledge of the structure and function of T. cruzi surface molecules is fundamental to understanding the mechanisms by which the parasite interacts with its host. T. cruzi infective forms engage a repertoire of surface and secreted molecules, some of which are involved in triggering signaling pathways both in the parasite and the host cell, leading to intracellular Ca²⁺ mobilization, a process essential for parasite internalization. Here, we described a novel family of T. cruzi surface membrane proteins (TcSMP), including their genomic distribution, expression and cellular localization. We studied the mechanism of action of TcSMP in host-cell invasion and proposed a triggering role for TcSMP in host-cell lysosome exocytosis during metacyclic internalization. TcSMP genes are conserved among different T. cruzi lineages and share orthologs in other Trypanosoma species. These results suggest that the diversification of TcSMP genes in mammalian trypanosomes occurred after continental drift. In T. cruzi this gene family expanded by gene duplication.

Potential Role of Acetyl-CoA Synthetase (acs) and Malate Dehydrogenase (mae) in the Evolution of the Acetate Switch in Bacteria and Archaea

Although many Archaea have AMP-Acs (acetyl-coenzyme A synthetase) and ADP-Acs, the extant methanogenic genus Methanosarcina is the only identified Archaeal genus that can utilize acetate via acetate kinase (Ack) and phosphotransacetylase (Pta). Despite the importance of ack as the potential urkinase in the ASKHA phosphotransferase superfamily, an origin hypothesis does not exist for the acetate kinase in Bacteria, Archaea, or Eukarya. Here we demonstrate that Archaeal AMP-Acs and ADP-Acs contain paralogous ATPase motifs previously identified in Ack, which demonstrate a novel relation between these proteins in Archaea. The identification of ATPase motif conservation and resulting structural features in AMP- and ADP-acetyl-CoA synthetase proteins in this study expand the ASKHA superfamily to include acetyl-CoA synthetase. Additional phylogenetic analysis showed that Pta and MaeB sequences had a common ancestor, and that the Pta lineage within the halophilc archaea was an ancestral lineage. These results suggested that divergence of a duplicated maeB within an ancient halophilic, archaeal lineage formed a putative pta ancestor. These results provide a potential scenario for the establishment of the Ack/Pta pathway and provide novel insight into the evolution of acetate metabolism for all three domains of life.

Chromosomal copy number variation reveals differential levels of genomic plasticity in distinct Trypanosoma cruzi strains

BACKGROUND

Trypanosoma cruzi, the etiologic agent of Chagas disease, is currently divided into six discrete typing units (DTUs), named TcI-TcVI. CL Brener, the reference strain of the T. cruzi genome project, is a hybrid with a genome assembled into 41 putative chromosomes. Gene copy number variation (CNV) is well documented as an important mechanism to enhance gene expression and variability in T. cruzi. Chromosomal CNV (CCNV) is another level of gene CNV in which whole blocks of genes are expanded simultaneously. Although the T. cruzi karyotype is not well defined, several studies have demonstrated a significant variation in the size and content of chromosomes between different T. cruzi strains. Despite these studies, the extent of diversity in CCNV among T. cruzi strains based on a read depth coverage analysis has not been determined.

RESULTS

We identify the CCNV in T. cruzi strains from the TcI, TcII and TcIII DTUs, by analyzing the depth coverage of short reads from these strains using the 41 CL Brener chromosomes as reference. This study led to the identification of a broader extent of CCNV in T. cruzi than was previously speculated. The TcI DTU strains have very few aneuploidies, while the strains from TcII and TcIII DTUs present a high degree of chromosomal expansions. Chromosome 31, which is the only chromosome that is supernumerary in all six T. cruzi samples evaluated in this study, is enriched with genes related to glycosylation pathways, highlighting the importance of glycosylation to parasite survival.

CONCLUSIONS

Increased gene copy number due to chromosome amplification may contribute to alterations in gene expression, which represents a strategy that may be crucial for parasites that mainly depend on post-transcriptional mechanisms to control gene expression.

Putative interchromosomal rearrangements in the hexaploid wheat (Triticum aestivum L.) genotype 'Chinese Spring' revealed by gene locations on homoeologous chromosomes

Background

Chromosomal rearrangements are a major driving force in shaping genome during evolution. Previous studies show that translocated genes could undergo elevated rates of evolution and recombination frequencies around these genes can be altered. Based on the recently released genome sequences of Triticum urartu, Aegilops tauschii, Brachypodium distachyon and bread wheat, an analysis of interchromosomal translocations in the hexaploid wheat genotype ‘Chinese Spring’ (‘CS’) was conducted based on chromosome shotgun sequences from individual chromosome arms of this genotype.

Results

A total of 720 genes representing putative interchromosomal rearrangements was identified. They were distributed across the 42 chromosome arms. About 59% of these translocated genes were those involved in the well-characterized translocations involving chromosomes 4A, 5A and 7B. The other 41% of the genes represent a large numbers of putative interchromosomal rearrangements which have not yet been described. The number of the putative translocation events in the D subgenome was about half of those presented in either the A or B subgenomes, which agreed well with that the times of interaction between the A and B subgenomes almost doubled that between either of them and the D subgenome.

Conclusions

The possible existence of a large number of interchromosomal rearrangements detected in this study provide further evidence that caution should be taken when using synteny in ordering sequence contigs or in cloning genes in hexaploid wheat. The identification of these putative translocations in ‘CS’ also provide a base for a systematic evaluation of their presence or absence in the full spectrum of bread wheat and its close relatives, which could have significant implications in a wide array of fields ranging from studies of systematics and evolution to practical breeding.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0313-5) contains supplementary material, which is available to authorized users.

Identification and analysis of ingi-related retroposons in the trypanosomatid genomes

Transposable elements (TE), defined as discrete pieces of DNA that can move from one site to another site in genomes, represent significant components of eukaryotic genomes, including trypanosomatids. Up to 5% of the trypanosomatid genome content is composed of retroposons of the ingi clade, further divided into subclades and subfamilies ranging from short extinct truncated elements (SIDER) to long active elements (ingi). Important differences in ingi-related retroposon content have been reported between trypanosomatid species. For instance, Leishmania spp. have expanded and recycled a whole SIDER family to fulfill an important biological pathway, i.e., regulation of gene expression, while trypanosome genomes are primarily composed of active elements. Here, we present an overview of the computational methods used to identify, annotate, and analyze ingi-related retroposons for providing a comprehensive picture of all these TE families in newly available trypanosomatid genome sequences.

In silico work flow for scaffold hopping in Leishmania

BACKGROUND

Leishmaniasis,a broad spectrum of diseases caused by several sister species of protozoa belonging to family trypanosomatidae and genus leishmania , generally affects poorer sections of the populace in third world countries. With the emergence of strains resistant to traditional therapies and the high cost of second line drugs which generally have severe side effects, it becomes imperative to continue the search for alternative drugs to combat the disease. In this work, the leishmanial genomes and the human genome have been compared to identify proteins unique to the parasite and whose structures (or those of close homologues) are available in the Protein Data Bank. Subsequent to the prioritization of these proteins (based on their essentiality, virulence factor etc.), inhibitors have been identified for a subset of these prospective drug targets by means of an exhaustive literature survey. A set of three dimensional protein-ligand complexes have been assembled from the list of leishmanial drug targets by culling structures from the Protein Data Bank or by means of template based homology modeling followed by ligand docking with the GOLD software. Based on these complexes several structure based pharmacophores have been designed and used to search for alternative inhibitors in the ZINC database.

RESULT

This process led to a list of prospective compounds which could serve as potential antileishmanials. These small molecules were also used to search the Drug Bank to identify prospective lead compounds already in use as approved drugs. Interestingly, paromomycin which is currently being used as an antileishmanial drug spontaneously appeared in the list, probably giving added confidence to the 'scaffold hopping' computational procedures adopted in this work.

CONCLUSIONS

The report thus provides the basis to experimentally verify several lead compounds for their predicted antileishmanial activity and includes several useful data bases of prospective drug targets in leishmania, their inhibitors and protein--inhibitor three dimensional complexes.

Pr77 and L1TcRz: A dual system within the 5'-end of L1Tc retrotransposon, internal promoter and HDV-like ribozyme

The sequence corresponding to the first 77 nucleotides of the L1Tc and NARTc non-LTR retrotransposons from Trypanosoma cruzi is an internal promoter (Pr77) that generates abundant, although poorly translatable, un-spliced transcripts. It has been recently described that L1TcRz, an HDV-like ribozyme, resides within the 5'-end of the RNA from the L1Tc and NARTc retrotransposons. Remarkably, the same first 77 nucleotides of L1Tc/NARTc elements comprise both the Pr77 internal promoter and the HDV-like L1TcRz. The L1TcRz cleaves on the 5'-side of the +1 nucleotide of the L1Tc element insuring that the promoter and the ribozyme functions travel with the transposon during retrotransposition. The ribozyme activity would prevent the mobilization of upstream sequences and insure the individuality of the L1Tc/NARTc copies transcribed from associated tandems. The Pr77/L1TcRz sequence is also found in other trypanosomatid's non-LTR retrotransposons and degenerated retroposons. The possible conservation of the ribozyme activity in a widely degenerated retrotransposon, as the Leishmania SIDERs, could indicate that the presence of this element and the catalytic activity could play some favorable genetic regulation. The functional implications of the Pr77/L1TcRz dual system in the regulation of the L1Tc/NARTc retrotransposons and in the gene expression of trypanosomatids are also discussed in this paper.

Genome of the avirulent human-infective trypanosome--Trypanosoma rangeli

Background

Trypanosoma rangeli is a hemoflagellate protozoan parasite infecting humans and other wild and domestic mammals across Central and South America. It does not cause human disease, but it can be mistaken for the etiologic agent of Chagas disease, Trypanosoma cruzi. We have sequenced the T. rangeli genome to provide new tools for elucidating the distinct and intriguing biology of this species and the key pathways related to interaction with its arthropod and mammalian hosts.

Methodology/Principal Findings

The T. rangeli haploid genome is ∼24 Mb in length, and is the smallest and least repetitive trypanosomatid genome sequenced thus far. This parasite genome has shorter subtelomeric sequences compared to those of T. cruzi and T. brucei; displays intraspecific karyotype variability and lacks minichromosomes. Of the predicted 7,613 protein coding sequences, functional annotations could be determined for 2,415, while 5,043 are hypothetical proteins, some with evidence of protein expression. 7,101 genes (93%) are shared with other trypanosomatids that infect humans. An ortholog of the dcl2 gene involved in the T. brucei RNAi pathway was found in T. rangeli, but the RNAi machinery is non-functional since the other genes in this pathway are pseudogenized. T. rangeli is highly susceptible to oxidative stress, a phenotype that may be explained by a smaller number of anti-oxidant defense enzymes and heat-shock proteins.

Conclusions/Significance

Phylogenetic comparison of nuclear and mitochondrial genes indicates that T. rangeli and T. cruzi are equidistant from T. brucei. In addition to revealing new aspects of trypanosome co-evolution within the vertebrate and invertebrate hosts, comparative genomic analysis with pathogenic trypanosomatids provides valuable new information that can be further explored with the aim of developing better diagnostic tools and/or therapeutic targets.

Comparative genomics is a powerful tool that affords detailed study of the genetic and evolutionary basis for aspects of lifecycles and pathologies caused by phylogenetically related pathogens. The reference genome sequences of three trypanosomatids, T. brucei, T. cruzi and L. major, and subsequent addition of multiple Leishmania and Trypanosoma genomes has provided data upon which large-scale investigations delineating the complex systems biology of these human parasites has been built. Here, we compare the annotated genome sequence of T. rangeli strain SC-58 to available genomic sequence and annotation data from related species. We provide analysis of gene content, genome architecture and key characteristics associated with the biology of this non-pathogenic trypanosome. Moreover, we report striking new genomic features of T. rangeli compared with its closest relative, T. cruzi, such as (1) considerably less amplification on the gene copy number within multigene virulence factor families such as MASPs, trans-sialidases and mucins; (2) a reduced repertoire of genes encoding anti-oxidant defense enzymes; and (3) the presence of vestigial orthologs of the RNAi machinery, which are insufficient to constitute a functional pathway. Overall, the genome of T. rangeli provides for a much better understanding of the identity, evolution, regulation and function of trypanosome virulence determinants for both mammalian host and insect vector.

Genome evolution in trypanosomatid parasites

A decade of genome sequencing has transformed our understanding of how trypanosomatid parasites have evolved and provided fresh impetus to explaining the origins of parasitism in the Kinetoplastida. In this review, I will consider the many ways in which genome sequences have influenced our view of genomic reduction in trypanosomatids; how species-specific genes, and the genomic domains they occupy, have illuminated the innovations in trypanosomatid genomes; and how comparative genomics has exposed the molecular mechanisms responsible for innovation and adaptation to a parasitic lifestyle.

Trypanosoma cruzi: Genome characterization of phosphatidylinositol kinase gene family (PIK and PIK-related) and identification of a novel PIK gene

Chagas disease is caused by the protozoan Trypanosoma cruzi which affects 10 million people worldwide. Very few kinases have been characterized in this parasite, including the phosphatidylinositol kinases (PIKs) that are at the heart of one of the major pathways of intracellular signal transduction. Recently, we have classified the PIK family in T. cruzi using five different models based on the presence of PIK conserved domains. In this study, we have mapped PIK genes to the chromosomes of two different T. cruzi lineages (G and CL Brener) and determined the cellular localization of two PIK members. The kinases have crucial roles in metabolism and are assumed to be conserved throughout evolution. For this reason, they should display a conserved localization within the same eukaryotic species. In spite of this, there is an extensive polymorphism regarding PIK localization at both genomic and cellular levels, among different T. cruzi isolates and between T. cruzi and Trypanosomabrucei, respectively. We showed in this study that the cellular localization of two PIK-related proteins (TOR1 and 2) in the T. cruzi lineage is distinct from that previously observed in T. brucei. In addition, we identified a new PIK gene with peculiar feature, that is, it codes for a FYVE domain at N-terminal position. FYVE-PIK genes are phylogenetically distant from the groups containing exclusively the FYVE or PIK domain. The FYVE-PIK architecture is only present in trypanosomatids and in virus such as Acanthamoeba mimivirus, suggesting a horizontal acquisition. Our Bayesian phylogenetic inference supports this hypothesis. The exact functions of this FYVE-PIK gene are unknown, but the presence of FYVE domain suggests a role in membranous compartments, such as endosome. Taken together, the data presented here strengthen the possibility that trypanosomatids are characterized by extensive genomic plasticity that may be considered in designing drugs and vaccines for prevention of Chagas disease.

Genomic and proteomic approaches for Chagas’ disease: critical analysis of diagnostic methods

Trypanosoma cruzi is the etiologic agent of Chagas' disease, a chronic inflammatory condition that results in heart and digestive complications. The first draft of the parasite genome is now complete and it is expected that, along with the published genomic and proteomic analyses discussed herein, it will lead to the identification of crucial genes and proteins directly associated with disease. This article reviews the current research trends addressing host-parasite interaction, parasite genetic variability and diagnosis. These advances will certainly bring about major developments not only in our understanding of Trypanosoma cruzi biology, but also in the application of new technologies to disease prevention and control.

A comparative analysis of trypanosomatid SNARE proteins

The Kinetoplastida are flagellated protozoa evolutionary distant and divergent from yeast and humans. Kinetoplastida include trypanosomatids, and a number of important pathogens. Trypanosoma brucei, Trypanosoma cruzi and Leishmania spp. inflict significant morbidity and mortality on humans and livestock as the etiological agents of human African trypanosomiasis, Chagas' disease and leishmaniasis respectively. For all of these organisms, intracellular trafficking is vital for maintenance of the host–pathogen interface, modulation/evasion of host immune system responses and nutrient uptake. Soluble N-ethylmaleimide-sensitive factor attachment protein receptors (SNAREs) are critical components of the intracellular trafficking machinery in eukaryotes, mediating membrane fusion and contributing to organelle specificity. We asked how the SNARE complement evolved across the trypanosomatids. An in silico search of the predicted proteomes of T. b. brucei and T. cruzi was used to identify candidate SNARE sequences. Phylogenetic analysis, including comparisons with yeast and human SNAREs, allowed assignment of trypanosomatid SNAREs to the Q or R subclass, as well as identification of several SNAREs orthologous with those of opisthokonts. Only limited variation in number and identity of SNAREs was found, with Leishmania major having 27 and T. brucei 26, suggesting a stable SNARE complement post-speciation. Expression analysis of T. brucei SNAREs revealed significant differential expression between mammalian and insect infective forms, especially within R and Qb-SNARE subclasses, suggesting possible roles in adaptation to different environments. For trypanosome SNAREs with clear orthologs in opisthokonts, the subcellular localization of TbVAMP7C is endosomal while both TbSyn5 and TbSyn16B are at the Golgi complex, which suggests conservation of localization and possibly also function. Despite highly distinct life styles, the complement of trypanosomatid SNAREs is quite stable between the three pathogenic lineages, suggesting establishment in the last common ancestor of trypanosomes and Leishmania. Developmental changes to SNARE mRNA levels between blood steam and procyclic life stages suggest that trypanosomes modulate SNARE functions via expression. Finally, the locations of some conserved SNAREs have been retained across the eukaryotic lineage.

•

SNARE proteins are essential components of intracellular transport.

•

These proteins exhibit considerable conservation across pathogenic trypanosomes.

•

Some trypanosome SNARE families are expanded or lost.

•

Developmental changes in trypanosome SNARE expression are apparent.

•

Orthologous SNAREs demonstrate conserved locations and hence function.

Genomic analysis of sequence-dependent DNA curvature in Leishmania

Leishmania major is a flagellated protozoan parasite of medical importance. Like other members of the Trypanosomatidae family, it possesses unique mechanisms of gene expression such as constitutive polycistronic transcription of directional gene clusters, gene amplification, mRNA trans-splicing, and extensive editing of mitochondrial transcripts. The molecular signals underlying most of these processes remain under investigation. In order to investigate the role of DNA secondary structure signals in gene expression, we carried out a genome-wide in silico analysis of the intrinsic DNA curvature. The L. major genome revealed a lower frequency of high intrinsic curvature regions as well as inter- and intra- chromosomal distribution heterogeneity, when compared to prokaryotic and eukaryotic organisms. Using a novel method aimed at detecting region-integrated intrinsic curvature (RIIC), high DNA curvature was found to be associated with regions implicated in transcription initiation. Those include divergent strand-switch regions between directional gene clusters and regions linked to markers of active transcription initiation such as acetylated H3 histone, TRF4 and SNAP50. These findings suggest a role for DNA curvature in transcription initiation in Leishmania supporting the relevance of DNA secondary structures signals.

The mitochondrial genome and a 60-kb nuclear DNA segment from Naegleria fowleri, the causative agent of primary amoebic meningoencephalitis

Naegleria fowleri is a unicellular eukaryote causing primary amoebic meningoencephalitis, a neuropathic disease killing 99% of those infected, usually within 7-14 days. Naegleria fowleri is found globally in regions including the US and Australia. The genome of the related nonpathogenic species Naegleria gruberi has been sequenced, but the genetic basis for N. fowleri pathogenicity is unclear. To generate such insight, we sequenced and assembled the mitochondrial genome and a 60-kb segment of nuclear genome from N. fowleri. The mitochondrial genome is highly similar to its counterpart in N. gruberi in gene complement and organization, while distinct lack of synteny is observed for the nuclear segments. Even in this short (60-kb) segment, we identified examples of potential factors for pathogenesis, including ten novel N. fowleri-specific genes. We also identified a homolog of cathepsin B; proteases proposed to be involved in the pathogenesis of diverse eukaryotic pathogens, including N. fowleri. Finally, we demonstrate a likely case of horizontal gene transfer between N. fowleri and two unrelated amoebae, one of which causes granulomatous amoebic encephalitis. This initial look into the N. fowleri nuclear genome has revealed several examples of potential pathogenesis factors, improving our understanding of a neglected pathogen of increasing global importance.

Anatomy and evolution of telomeric and subtelomeric regions in the human protozoan parasite Trypanosoma cruzi

BACKGROUND

The subtelomeres of many protozoa are highly enriched in genes with roles in niche adaptation. T. cruzi trypomastigotes express surface proteins from Trans-Sialidase (TS) and Dispersed Gene Family-1 (DGF-1) superfamilies which are implicated in host cell invasion. Single populations of T. cruzi may express different antigenic forms of TSs. Analysis of TS genes located at the telomeres suggests that chromosome ends could have been the sites where new TS variants were generated. The aim of this study is to characterize telomeric and subtelomeric regions of T. cruzi available in TriTrypDB and connect the sequences of telomeres to T. cruzi working draft sequence.

RESULTS

We first identified contigs carrying the telomeric repeat (TTAGGG). Of 49 contigs identified, 45 have telomeric repeats at one end, whereas in four contigs the repeats are located internally. All contigs display a conserved telomeric junction sequence adjacent to the hexamer repeats which represents a signature of T. cruzi chromosome ends. We found that 40 telomeric contigs are located on T. cruzi chromosome-sized scaffolds. In addition, we were able to map several telomeric ends to the chromosomal bands separated by pulsed-field gel electrophoresis.The subtelomeric sequence structure varies widely, mainly as a result of large differences in the relative abundance and organization of genes encoding surface proteins (TS and DGF-1), retrotransposon hot spot genes (RHS), retrotransposon elements, RNA-helicase and N-acetyltransferase genes. While the subtelomeric regions are enriched in pseudogenes, they also contain complete gene sequences matching both known and unknown expressed genes, indicating that these regions do not consist of nonfunctional DNA but are instead functional parts of the expressed genome. The size of the subtelomeric regions varies from 5 to 182 kb; the smaller of these regions could have been generated by a recent chromosome breakage and telomere healing event.

CONCLUSIONS

The lack of synteny in the subtelomeric regions suggests that genes located in these regions are subject to recombination, which increases their variability, even among homologous chromosomes. The presence of typical subtelomeric genes can increase the chance of homologous recombination mechanisms or microhomology-mediated end joining, which may use these regions for the pairing and recombination of free ends.

Genome size, karyotype polymorphism and chromosomal evolution in Trypanosoma cruzi

Background

The Trypanosoma cruzi genome was sequenced from a hybrid strain (CL Brener). However, high allelic variation and the repetitive nature of the genome have prevented the complete linear sequence of chromosomes being determined. Determining the full complement of chromosomes and establishing syntenic groups will be important in defining the structure of T. cruzi chromosomes. A large amount of information is now available for T. cruzi and Trypanosoma brucei, providing the opportunity to compare and describe the overall patterns of chromosomal evolution in these parasites.

Methodology/Principal Findings

The genome sizes, repetitive DNA contents, and the numbers and sizes of chromosomes of nine strains of T. cruzi from four lineages (TcI, TcII, TcV and TcVI) were determined. The genome of the TcI group was statistically smaller than other lineages, with the exception of the TcI isolate Tc1161 (José-IMT). Satellite DNA content was correlated with genome size for all isolates, but this was not accompanied by simultaneous amplification of retrotransposons. Regardless of chromosomal polymorphism, large syntenic groups are conserved among T. cruzi lineages. Duplicated chromosome-sized regions were identified and could be retained as paralogous loci, increasing the dosage of several genes. By comparing T. cruzi and T. brucei chromosomes, homologous chromosomal regions in T. brucei were identified. Chromosomes Tb9 and Tb11 of T. brucei share regions of syntenic homology with three and six T. cruzi chromosomal bands, respectively.

Conclusions

Despite genome size variation and karyotype polymorphism, T. cruzi lineages exhibit conservation of chromosome structure. Several syntenic groups are conserved among all isolates analyzed in this study. The syntenic regions are larger than expected if rearrangements occur randomly, suggesting that they are conserved owing to positive selection. Mapping of the syntenic regions on T. cruzi chromosomal bands provides evidence for the occurrence of fusion and split events involving T. brucei and T. cruzi chromosomes.

Pathogenesis of chagas' disease: parasite persistence and autoimmunity

Acute Trypanosoma cruzi infections can be asymptomatic, but chronically infected individuals can die of Chagas' disease. The transfer of the parasite mitochondrial kinetoplast DNA (kDNA) minicircle to the genome of chagasic patients can explain the pathogenesis of the disease; in cases of Chagas' disease with evident cardiomyopathy, the kDNA minicircles integrate mainly into retrotransposons at several chromosomes, but the minicircles are also detected in coding regions of genes that regulate cell growth, differentiation, and immune responses. An accurate evaluation of the role played by the genotype alterations in the autoimmune rejection of self-tissues in Chagas' disease is achieved with the cross-kingdom chicken model system, which is refractory to T. cruzi infections. The inoculation of T. cruzi into embryonated eggs prior to incubation generates parasite-free chicks, which retain the kDNA minicircle sequence mainly in the macrochromosome coding genes. Crossbreeding transfers the kDNA mutations to the chicken progeny. The kDNA-mutated chickens develop severe cardiomyopathy in adult life and die of heart failure. The phenotyping of the lesions revealed that cytotoxic CD45, CD8(+) γδ, and CD8α(+) T lymphocytes carry out the rejection of the chicken heart. These results suggest that the inflammatory cardiomyopathy of Chagas' disease is a genetically driven autoimmune disease.

Cell biology of the trypanosome genome

Trypanosomes are a group of protozoan eukaryotes, many of which are major parasites of humans and livestock. The genomes of trypanosomes and their modes of gene expression differ in several important aspects from those of other eukaryotic model organisms. Protein-coding genes are organized in large directional gene clusters on a genome-wide scale, and their polycistronic transcription is not generally regulated at initiation. Transcripts from these polycistrons are processed by global trans-splicing of pre-mRNA. Furthermore, in African trypanosomes, some protein-coding genes are transcribed by a multifunctional RNA polymerase I from a specialized extranucleolar compartment. The primary DNA sequence of the trypanosome genomes and their cellular organization have usually been treated as separate entities. However, it is becoming increasingly clear that in order to understand how a genome functions in a living cell, we will need to unravel how the one-dimensional genomic sequence and its trans-acting factors are arranged in the three-dimensional space of the eukaryotic nucleus. Understanding this cell biology of the genome will be crucial if we are to elucidate the genetic control mechanisms of parasitism. Here, we integrate the concepts of nuclear architecture, deduced largely from studies of yeast and mammalian nuclei, with recent developments in our knowledge of the trypanosome genome, gene expression, and nuclear organization. We also compare this nuclear organization to those in other systems in order to shed light on the evolution of nuclear architecture in eukaryotes.

The biology and evolution of transposable elements in parasites

Transposable elements (TEs) are dynamic elements that can reshape host genomes by generating rearrangements with the potential to create or disrupt genes, to shuffle existing genes, and to modulate their patterns of expression. In the genomes of parasites that infect mammals several TEs have been identified that probably have been maintained throughout evolution due to their contribution to gene function and regulation of gene expression. This review addresses how TEs are organized, how they colonize the genomes of mammalian parasites, the functional role these elements play in parasite biology, and the interactions between these elements and the parasite genome.

Gene expression in trypanosomatid parasites

The parasites Leishmania spp., Trypanosoma brucei, and Trypanosoma cruzi are the trypanosomatid protozoa that cause the deadly human diseases leishmaniasis, African sleeping sickness, and Chagas disease, respectively. These organisms possess unique mechanisms for gene expression such as constitutive polycistronic transcription of protein-coding genes and trans-splicing. Little is known about either the DNA sequences or the proteins that are involved in the initiation and termination of transcription in trypanosomatids. In silico analyses of the genome databases of these parasites led to the identification of a small number of proteins involved in gene expression. However, functional studies have revealed that trypanosomatids have more general transcription factors than originally estimated. Many posttranslational histone modifications, histone variants, and chromatin modifying enzymes have been identified in trypanosomatids, and recent genome-wide studies showed that epigenetic regulation might play a very important role in gene expression in this group of parasites. Here, we review and comment on the most recent findings related to transcription initiation and termination in trypanosomatid protozoa.

Cascades of convergent evolution: the corresponding evolutionary histories of euglenozoans and dinoflagellates

The majority of eukaryotic diversity is hidden in protists, yet our current knowledge of processes and structures in the eukaryotic cell is almost exclusively derived from multicellular organisms. The increasing sensitivity of molecular methods and growing interest in microeukaryotes has only recently demonstrated that many features so far considered to be universal for eukaryotes actually exist in strikingly different versions. In other words, during their long evolutionary histories, protists have solved general biological problems in many more ways than previously appreciated. Interestingly, some groups have broken more rules than others, and the Euglenozoa and the Alveolata stand out in this respect. A review of the numerous odd features in these 2 groups allows us to draw attention to the high level of convergent evolution in protists, which perhaps reflects the limits that certain features can be altered. Moreover, the appearance of one deviation in an ancestor can constrain the set of possible downstream deviations in its descendents, so features that might be independent functionally, can still be evolutionarily linked. What functional advantage may be conferred by the excessive complexity of euglenozoan and alveolate gene expression, organellar genome structure, and RNA editing and processing has been thoroughly debated, but we suggest these are more likely the products of constructive neutral evolution, and as such do not necessarily confer any selective advantage at all.

Organization and evolution of two SIDER retroposon subfamilies and their impact on the Leishmania genome

Background

We have recently identified two large families of extinct transposable elements termed Short Interspersed DEgenerated Retroposons (SIDERs) in the parasitic protozoan Leishmania major. The characterization of SIDER elements was limited to the SIDER2 subfamily, although members of both subfamilies have been shown to play a role in the regulation of gene expression at the post-transcriptional level. Apparent functional domestication of SIDERs prompted further investigation of their characterization, dissemination and evolution throughout the Leishmania genus, with particular attention to the disregarded SIDER1 subfamily.

Results

Using optimized statistical profiles of both SIDER1 and SIDER2 subgroups, we report the first automated and highly sensitive annotation of SIDERs in the genomes of L. infantum, L. braziliensis and L. major. SIDER annotations were combined to in-silico mRNA extremity predictions to generate a detailed distribution map of the repeat family, hence uncovering an enrichment of antisense-oriented SIDER repeats between the polyadenylation and trans-splicing sites of intergenic regions, in contrast to the exclusive sense orientation of SIDER elements within 3'UTRs. Our data indicate that SIDER elements are quite uniformly dispersed throughout all three genomes and that their distribution is generally syntenic. However, only 47.4% of orthologous genes harbor a SIDER element in all three species. There is evidence for species-specific enrichment of SIDERs and for their preferential association, especially for SIDER2s, with different metabolic functions. Investigation of the sequence attributes and evolutionary relationship of SIDERs to other trypanosomatid retroposons reveals that SIDER1 is a truncated version of extinct autonomous ingi-like retroposons (DIREs), which were functional in the ancestral Leishmania genome.

Conclusion

A detailed characterization of the sequence traits for both SIDER subfamilies unveils major differences. The SIDER1 subfamily is more heterogeneous and shows an evolutionary link with vestigial DIRE retroposons as previously observed for the ingi/RIME and L1Tc/NARTc couples identified in the T. brucei and T. cruzi genomes, whereas no identified DIREs are related to SIDER2 sequences. Although SIDER1s and SIDER2s display equivalent genomic distribution globally, the varying degrees of sequence conservation, preferential genomic disposition, and differential association to orthologous genes allude to an intricate web of SIDER assimilation in these parasitic organisms.

An analysis of synteny of Arachis with Lotus and Medicago sheds new light on the structure, stability and evolution of legume genomes

Background

Most agriculturally important legumes fall within two sub-clades of the Papilionoid legumes: the Phaseoloids and Galegoids, which diverged about 50 Mya. The Phaseoloids are mostly tropical and include crops such as common bean and soybean. The Galegoids are mostly temperate and include clover, fava bean and the model legumes Lotus and Medicago (both with substantially sequenced genomes). In contrast, peanut (Arachis hypogaea) falls in the Dalbergioid clade which is more basal in its divergence within the Papilionoids. The aim of this work was to integrate the genetic map of Arachis with Lotus and Medicago and improve our understanding of the Arachis genome and legume genomes in general. To do this we placed on the Arachis map, comparative anchor markers defined using a previously described bioinformatics pipeline. Also we investigated the possible role of transposons in the patterns of synteny that were observed.

Results

The Arachis genetic map was substantially aligned with Lotus and Medicago with most synteny blocks presenting a single main affinity to each genome. This indicates that the last common whole genome duplication within the Papilionoid legumes predated the divergence of Arachis from the Galegoids and Phaseoloids sufficiently that the common ancestral genome was substantially diploidized. The Arachis and model legume genomes comparison made here, together with a previously published comparison of Lotus and Medicago allowed all possible Arachis-Lotus-Medicago species by species comparisons to be made and genome syntenies observed. Distinct conserved synteny blocks and non-conserved regions were present in all genome comparisons, implying that certain legume genomic regions are consistently more stable during evolution than others. We found that in Medicago and possibly also in Lotus, retrotransposons tend to be more frequent in the variable regions. Furthermore, while these variable regions generally have lower densities of single copy genes than the more conserved regions, some harbor high densities of the fast evolving disease resistance genes.

Conclusion

We suggest that gene space in Papilionoids may be divided into two broadly defined components: more conserved regions which tend to have low retrotransposon densities and are relatively stable during evolution; and variable regions that tend to have high retrotransposon densities, and whose frequent restructuring may fuel the evolution of some gene families.

Role of transposable elements in trypanosomatids

Transposable elements constitute 2-5% of the genome content in trypanosomatid parasites. Some of them are involved in critical cellular functions, such as the regulation of gene expression in Leishmania spp. In this review, we highlight the remarkable role extinct transposable elements can play as the source of potential new functions.

In silico, biologically-inspired modelling of genomic variation generation in surface proteins of Trypanosoma cruzi

BACKGROUND

Protozoan parasites improve the likelihood of invading or adapting to the host through their capacity to present a large repertoire of surface molecules. The understanding of the mechanisms underlying the generation of antigenic diversity is crucial to aid in the development of therapies and the study of evolution. Despite advances driven by molecular biology and genomics, there is a need to gain a deeper understanding of key properties that may facilitate variation generation, models for explaining the role of genomic re-arrangements and the characterisation of surface protein families on the basis of their capacity to generate variation. Computer models may be implemented to explore, visualise and estimate the variation generation capacity of gene families in a dynamic fashion. In this paper we report the dynamic simulation of genomic variation using real T. cruzi coding sequences as inputs to a computational simulation system. The effects of random, multiple-point mutations and gene conversions on genomic variation generation were quantitatively estimated and visualised. Simulations were also implemented to investigate the potential role of pseudogenes as a source of antigenic variation in T. cruzi.

RESULTS

Computational models of variation generation were applied to real coding sequences from surface proteins in T. cruzi: trans-sialidase-like proteins and putative surface protein dispersed gene family-1. In the simulations the sequences self-replicated, mutated and re-arranged during thousands of generations. Simulations were implemented for different mutation rates to estimate the relative robustness of the protein families in the face of DNA multiple-point mutations and sequence re-arrangements. The gene super-families and families showed distinguishing evolutionary responses, which may be used to characterise them on the basis of their capacity to generate variability. The simulations showed that sequences from T. cruzi nuclear genes tend to be relatively more robust against random, multiple-point mutations than those obtained from surface protein genes. Simulations also showed that a gene conversion model may act as an effective variation generation mechanism. Differential variation responses can be used to characterise the sequence groups under study. For example, unlike other families, sequences from the DGF1 family have the capacity to maximise variation at the amino acid level under relatively low mutation rates and through gene conversion. However, in relation to the other protein families, they exhibit more robust behaviour in response to more severe modifications through intra-family genomic sequence exchange. Independent simulations indicate that DGF1 pseudogenes might play a role in the generation of greater genomic variation in the DFG1 gene family through gene conversion under different experimental conditions.

CONCLUSION

Digital, dynamic simulations may be implemented to characterise gene families on the basis of their capacity to generate variation in the face of genomic perturbations. Such simulations may be useful to explore antigenic variation mechanisms and hypotheses about robustness at the genomic level. This investigation illustrated how sequences derived from surface protein genes and computer simulations can be used to investigate variation generation mechanisms. Such in silico experiments of self-replicating sequences undergoing random mutations and genomic re-arrangements can offer insights into the diversity generation potential of the genes under study. Biologically-inspired simulations may support the study of genomic variation mechanisms in pathogens whose genomes have been recently sequenced.

Comparative genomics: from genotype to disease phenotype in the leishmaniases

Recent progress in sequencing the genomes of several Leishmania species, causative agents of cutaneous, mucocutaneous and visceral leishmaniasis, is revealing unusual features of potential relevance to parasite virulence and pathogenesis in the host. While the genomes of Leishmania major, Leishmania braziliensis and Leishmania infantum are highly similar in content and organisation, species-specific genes and mechanisms distinguish one from another. In particular, the presence of retrotransposons and the components of a putative RNA interference machinery in L. braziliensis suggest the potential for both greater diversity and more tractable experimentation in this Leishmania Viannia species.

A postgenomic view of the heat shock proteins in kinetoplastids

The kinetoplastids Leishmania major, Trypanosoma brucei and Trypanosoma cruzi are causative agents of a diverse spectrum of human diseases: leishmaniasis, sleeping sickness and Chagas' disease, respectively. These protozoa possess digenetic life cycles that involve development in mammalian and insect hosts. It is generally accepted that temperature is a triggering factor of the developmental programme allowing the adaptation of the parasite to the mammalian conditions. The heat shock response is a general homeostatic mechanism that protects cells from the deleterious effects of environmental stresses, such as heat. This response is universal and includes the synthesis of the heat-shock proteins (HSPs). In this review, we summarize the salient features of the different HSP families and describe their main cellular functions. In parallel, we analyse the composition of these families in kinetoplastids according to literature data and our understanding of genome sequence data. The genome sequences of these parasites have been recently completed. The HSP families described here are: HSP110, HSP104, group I chaperonins, HSP90, HSP70, HSP40 and small HSPs. All these families are widely represented in these parasites. In particular, kinetoplastids possess an unprecedented number of members of the HSP70, HSP60 and HSP40 families, suggesting key roles for these HSPs in their biology.