Genome organization and DNA accessibility control antigenic variation in trypanosomes

Protein phosphatase PP1 regulation of RNA polymerase II transcription termination and allelic exclusion of VSG genes in trypanosomes

The genomes of Leishmania and trypanosomes are organized into polycistronic transcription units flanked by a modified DNA base J involved in promoting RNA polymerase II (Pol II) termination. We recently characterized a Leishmania complex containing a J-binding protein, PP1 protein phosphatase 1, and PP1 regulatory protein (PNUTS) that controls transcription termination potentially via dephosphorylation of Pol II by PP1. While T. brucei contains eight PP1 isoforms, none purified with the PNUTS complex, complicating the analysis of PP1 function in termination. We now demonstrate that the PP1-binding motif of TbPNUTS is required for function in termination in vivo and that TbPP1-1 modulates Pol II termination in T. brucei and dephosphorylation of the large subunit of Pol II. PP1-1 knock-down results in increased cellular levels of phosphorylated RPB1 accompanied by readthrough transcription and aberrant transcription of the chromosome by Pol II, including Pol I transcribed loci that are typically silent, such as telomeric VSG expression sites involved in antigenic variation. These results provide important insights into the mechanism underlying Pol II transcription termination in primitive eukaryotes that rely on polycistronic transcription and maintain allelic exclusion of VSG genes.

Mono-allelic epigenetic regulation of bi-directional polycistronic transcription initiation by RNA Polymerase II in Trypanosoma brucei

Unique for a eukaryote, protein-coding genes in trypanosomes are arranged in polycistronic units (PTUs). This genome arrangement has led to a model where Pol II transcription of PTUs is unregulated. The initial step in trypanosome lytic factor (TLF) mediated lysis of Trypanosoma brucei requires high affinity haptoglobin/hemoglobin receptor (HpHbR) binding. Here we demonstrate that by in vitro selection with TLF, resistance is obtained in a stepwise process correlating with loss of HpHbR expression at an allelic level. RNA-seq, Pol II ChIP and run-on analysis indicate HpHbR silencing is at the transcriptional level, where loss of Pol II binding at the promoter region specifically shuts down transcription of the HpHbR containing gene cluster and the adjacent opposing gene cluster. Reversible transcriptional silencing of the divergent PTUs correlates with DNA base J modification of the shared promoter region. Therefore, epigenetic mechanisms exist to regulate gene expression via Pol II transcription initiation of gene clusters in a mono-allelic fashion. These findings suggest epigenetic chromatin-based regulation of gene expression is deeply conserved among eukaryotes, including primitive eukaryotes that rely on polycistronic transcription.

Trypanosoma brucei Invariant Surface Glycoprotein 75 Is an Immunoglobulin Fc Receptor Inhibiting Complement Activation and Antibody-Mediated Cellular Phagocytosis

Various subspecies of the unicellular parasite Trypanosoma brucei cause sleeping sickness, a neglected tropical disease affecting millions of individuals and domestic animals. Immune evasion mechanisms play a pivotal role in parasite survival within the host and enable the parasite to establish a chronic infection. In particular, the rapid switching of variant surface glycoproteins covering a large proportion of the parasite's surface enables the parasite to avoid clearance by the adaptive immune system of the host. In this article, we present the crystal structure and discover an immune-evasive function of the extracellular region of the T. brucei invariant surface gp75 (ISG75). Structural analysis determined that the ISG75 ectodomain is organized as a globular head domain and a long slender coiled-coil domain. Subsequent ligand screening and binding analysis determined that the head domain of ISG75 confers interaction with the Fc region of all subclasses of human IgG. Importantly, the ISG75-IgG interaction strongly inhibits both activation of the classical complement pathway and Ab-dependent cellular phagocytosis by competing with C1q and host cell FcγR CD32. Our data reveal a novel immune evasion mechanism of T. brucei, with ISG75 able to inactivate the activities of Abs recognizing the parasite surface proteins.

Smoother: on-the-fly processing of interactome data using prefix sums

Nucleic acid interactome data, such as chromosome conformation capture data and RNA-DNA interactome data, are currently analyzed via pipelines that must be rerun for each new parameter set. A more dynamic approach is desirable since the optimal parameter set is commonly unknown ahead of time and rerunning pipelines is a time-consuming process. We have developed an approach fast enough to process interactome data on-the-fly using a sparse prefix sum index. With this index, we created Smoother, a flexible, multifeatured visualization and analysis tool that allows interactive filtering, e.g. by mapping quality, almost instant comparisons between different normalization approaches, e.g. iterative correction, and ploidy correction. Further, Smoother can overlay other sequencing data or genomic annotations, compare different samples, and perform virtual 4C analysis. Smoother permits a novel way to interact with and explore interactome data, fostering comprehensive, high-quality data analysis. Smoother is available at https://github.com/Siegel-Lab/BioSmoother under the MIT license.

Two DOT1 enzymes cooperatively mediate efficient ubiquitin-independent histone H3 lysine 76 tri-methylation in kinetoplastids

In higher eukaryotes, a single DOT1 histone H3 lysine 79 (H3K79) methyltransferase processively produces H3K79me2/me3 through histone H2B mono-ubiquitin interaction, while the kinetoplastid Trypanosoma brucei di-methyltransferase DOT1A and tri-methyltransferase DOT1B efficiently methylate the homologous H3K76 without H2B mono-ubiquitination. Based on structural and biochemical analyses of DOT1A, we identify key residues in the methyltransferase motifs VI and X for efficient ubiquitin-independent H3K76 methylation in kinetoplastids. Substitution of a basic to an acidic residue within motif VI (Gx6K) is essential to stabilize the DOT1A enzyme-substrate complex, while substitution of the motif X sequence VYGE by CAKS renders a rigid active-site loop flexible, implying a distinct mechanism of substrate recognition. We further reveal distinct methylation kinetics and substrate preferences of DOT1A (H3K76me0) and DOT1B (DOT1A products H3K76me1/me2) in vitro, determined by a Ser and Ala residue within motif IV, respectively, enabling DOT1A and DOT1B to mediate efficient H3K76 tri-methylation non-processively but cooperatively, and suggesting why kinetoplastids have evolved two DOT1 enzymes.

Multiomics analysis reveals B. MO1 as a distinct Babesia species and provides insights into its evolution and virulence

Babesiosis, caused by protozoan parasites of the genus Babesia , is an emerging tick-borne disease of significance for both human and animal health. Babesia parasites infect erythrocytes of vertebrate hosts where they develop and multiply rapidly to cause the pathological symptoms associated with the disease. The identification of various Babesia species underscores the ongoing risk of new zoonotic pathogens capable of infecting humans, a concern amplified by anthropogenic activities and environmental shifts impacting the distribution and transmission dynamics of parasites, their vectors, and reservoir hosts. One such species, Babesia MO1, previously implicated in severe cases of human babesiosis in the midwestern United States, was initially considered closely related to B. divergens , the predominant agent of human babesiosis in Europe. Yet, uncertainties persist regarding whether these pathogens represent distinct variants of the same species or are entirely separate species. We show that although both B. MO1 and B. divergens share similar genome sizes, comprising three nuclear chromosomes, one linear mitochondrial chromosome, and one circular apicoplast chromosome, major differences exist in terms of genomic sequence divergence, gene functions, transcription profiles, replication rates and susceptibility to antiparasitic drugs. Furthermore, both pathogens have evolved distinct classes of multigene families, crucial for their pathogenicity and adaptation to specific mammalian hosts. Leveraging genomic information for B. MO1, B. divergens , and other members of the Babesiidae family within Apicomplexa provides valuable insights into the evolution, diversity, and virulence of these parasites. This knowledge serves as a critical tool in preemptively addressing the emergence and rapid transmission of more virulent strains.

Unwrap RAP1's Mystery at Kinetoplastid Telomeres

Although located at the chromosome end, telomeres are an essential chromosome component that helps maintain genome integrity and chromosome stability from protozoa to mammals. The role of telomere proteins in chromosome end protection is conserved, where they suppress various DNA damage response machineries and block nucleolytic degradation of the natural chromosome ends, although the detailed underlying mechanisms are not identical. In addition, the specialized telomere structure exerts a repressive epigenetic effect on expression of genes located at subtelomeres in a number of eukaryotic organisms. This so-called telomeric silencing also affects virulence of a number of microbial pathogens that undergo antigenic variation/phenotypic switching. Telomere proteins, particularly the RAP1 homologs, have been shown to be a key player for telomeric silencing. RAP1 homologs also suppress the expression of Telomere Repeat-containing RNA (TERRA), which is linked to their roles in telomere stability maintenance. The functions of RAP1s in suppressing telomere recombination are largely conserved from kinetoplastids to mammals. However, the underlying mechanisms of RAP1-mediated telomeric silencing have many species-specific features. In this review, I will focus on Trypanosoma brucei RAP1's functions in suppressing telomeric/subtelomeric DNA recombination and in the regulation of monoallelic expression of subtelomere-located major surface antigen genes. Common and unique mechanisms will be compared among RAP1 homologs, and their implications will be discussed.

Trypanosoma brucei proliferates normally even after losing all S-adenosylhomocysteine hydrolase genes

S-adenosylhomocysteine (SAH) hydrolase is the enzyme responsible for breaking down SAH into adenosine and homocysteine. It has long been believed that a deficiency of this enzyme leads to SAH accumulation, subsequently inhibiting methyltransferases responsible for nucleic acids and proteins, which severely affects cell proliferation. To investigate whether targeting this enzyme could be a viable strategy to combat Trypanosoma brucei, the causative agent of human African trypanosomiasis, we created a null mutant of the SAH hydrolase gene in T. brucei using the Cre/loxP system and conducted a phenotype analysis. Surprisingly, the null mutant, where all five SAH hydrolase gene loci were deleted, exhibited normal proliferation despite the observed SAH accumulation. These findings suggest that inhibiting SAH hydrolase may not be an effective approach to suppressing T. brucei proliferation, making the enzyme a less promising target for antitrypanosome drug development.

An allele-selective inter-chromosomal protein bridge supports monogenic antigen expression in the African trypanosome

UPF1-like helicases play roles in telomeric heterochromatin formation and X-chromosome inactivation, and also in monogenic variant surface glycoprotein (VSG) expression via VSG exclusion-factor-2 (VEX2), a UPF1-related protein in the African trypanosome. We show that VEX2 associates with chromatin specifically at the single active VSG expression site on chromosome 6, forming an allele-selective connection, via VEX1, to the trans-splicing locus on chromosome 9, physically bridging two chromosomes and the VSG transcription and splicing compartments. We further show that the VEX-complex is multimeric and self-regulates turnover to tightly control its abundance. Using single cell transcriptomics following VEX2-depletion, we observed simultaneous derepression of many other telomeric VSGs and multi-allelic VSG expression in individual cells. Thus, an allele-selective, inter-chromosomal, and self-limiting VEX1-2 bridge supports monogenic VSG expression and multi-allelic VSG exclusion.

Imaging of genomic loci in Trypanosoma brucei using an optimised LacO-LacI system

Visualisation of genomic loci by microscopy is essential for understanding nuclear organisation, particularly at the single cell level. One powerful technique for studying the positioning of genomic loci is through the Lac Operator-Lac Repressor (LacO-LacI) system, in which LacO repeats introduced into a specific genomic locus can be visualised through expression of a LacI-protein fused to a fluorescent tag. First utilised in Trypanosoma brucei over 20 years ago, we have now optimised this system with short, stabilised LacO repeats of less than 2 kb paired with a constitutively expressed mNeongreen::LacI fusion protein to facilitate visualisation of genomic loci. We demonstrate the compatibility of this system with super-resolution microscopy and propose its suitability for multiplexing with inducible RNAi or protein over expression which will allow analysis of nuclear organisation after perturbation of gene expression.

Pathogenicity and virulence of African trypanosomes: From laboratory models to clinically relevant hosts

African trypanosomes are vector-borne protozoa, which cause significant human and animal disease across sub-Saharan Africa, and animal disease across Asia and South America. In humans, infection is caused by variants of Trypanosoma brucei, and is characterized by varying rate of progression to neurological disease, caused by parasites exiting the vasculature and entering the brain. Animal disease is caused by multiple species of trypanosome, primarily T. congolense, T. vivax, and T. brucei. These trypanosomes also infect multiple species of mammalian host, and this complexity of trypanosome and host diversity is reflected in the spectrum of severity of disease in animal trypanosomiasis, ranging from hyperacute infections associated with mortality to long-term chronic infections, and is also a main reason why designing interventions for animal trypanosomiasis is so challenging. In this review, we will provide an overview of the current understanding of trypanosome determinants of infection progression and severity, covering laboratory models of disease, as well as human and livestock disease. We will also highlight gaps in knowledge and capabilities, which represent opportunities to both further our fundamental understanding of how trypanosomes cause disease, as well as facilitating the development of the novel interventions that are so badly needed to reduce the burden of disease caused by these important pathogens.

PI(3,4,5)P3 allosteric regulation of repressor activator protein 1 controls antigenic variation in trypanosomes

African trypanosomes evade host immune clearance by antigenic variation, causing persistent infections in humans and animals. These parasites express a homogeneous surface coat of variant surface glycoproteins (VSGs). They transcribe one out of hundreds of VSG genes at a time from telomeric expression sites (ESs) and periodically change the VSG expressed by transcriptional switching or recombination. The mechanisms underlying the control of VSG switching and its developmental silencing remain elusive. We report that telomeric ES activation and silencing entail an on/off genetic switch controlled by a nuclear phosphoinositide signaling system. This system includes a nuclear phosphatidylinositol 5-phosphatase (PIP5Pase), its substrate PI(3,4,5)P3, and the repressor-activator protein 1 (RAP1). RAP1 binds to ES sequences flanking VSG genes via its DNA binding domains and represses VSG transcription. In contrast, PI(3,4,5)P3 binds to the N-terminus of RAP1 and controls its DNA binding activity. Transient inactivation of PIP5Pase results in the accumulation of nuclear PI(3,4,5)P3, which binds RAP1 and displaces it from ESs, activating transcription of silent ESs and VSG switching. The system is also required for the developmental silencing of VSG genes. The data provides a mechanism controlling reversible telomere silencing essential for the periodic switching in VSG expression and its developmental regulation.

RAD51-mediated R-loop formation acts to repair transcription-associated DNA breaks driving antigenic variation in Trypanosoma brucei

RNA-DNA hybrids are epigenetic features of all genomes that intersect with many processes, including transcription, telomere homeostasis, and centromere function. Increasing evidence suggests that RNA-DNA hybrids can provide two conflicting roles in the maintenance and transmission of genomes: They can be the triggers of DNA damage, leading to genome change, or can aid the DNA repair processes needed to respond to DNA lesions. Evasion of host immunity by African trypanosomes, such as Trypanosoma brucei, relies on targeted recombination of silent Variant Surface Glycoprotein (VSG) genes into a specialized telomeric locus that directs transcription of just one VSG from thousands. How such VSG recombination is targeted and initiated is unclear. Here, we show that a key enzyme of T. brucei homologous recombination, RAD51, interacts with RNA-DNA hybrids. In addition, we show that RNA-DNA hybrids display a genome-wide colocalization with DNA breaks and that this relationship is impaired by mutation of RAD51. Finally, we show that RAD51 acts to repair highly abundant, localised DNA breaks at the single transcribed VSG and that mutation of RAD51 alters RNA-DNA hybrid abundance at 70 bp repeats both around the transcribed VSG and across the silent VSG archive. This work reveals a widespread, generalised role for RNA-DNA hybrids in directing RAD51 activity during recombination and uncovers a specialised application of this interplay during targeted DNA break repair needed for the critical T. brucei immune evasion reaction of antigenic variation.

Telomere maintenance in African trypanosomes

Telomere maintenance is essential for genome integrity and chromosome stability in eukaryotic cells harboring linear chromosomes, as telomere forms a specialized structure to mask the natural chromosome ends from DNA damage repair machineries and to prevent nucleolytic degradation of the telomeric DNA. In Trypanosoma brucei and several other microbial pathogens, virulence genes involved in antigenic variation, a key pathogenesis mechanism essential for host immune evasion and long-term infections, are located at subtelomeres, and expression and switching of these major surface antigens are regulated by telomere proteins and the telomere structure. Therefore, understanding telomere maintenance mechanisms and how these pathogens achieve a balance between stability and plasticity at telomere/subtelomere will help develop better means to eradicate human diseases caused by these pathogens. Telomere replication faces several challenges, and the "end replication problem" is a key obstacle that can cause progressive telomere shortening in proliferating cells. To overcome this challenge, most eukaryotes use telomerase to extend the G-rich telomere strand. In addition, a number of telomere proteins use sophisticated mechanisms to coordinate the telomerase-mediated de novo telomere G-strand synthesis and the telomere C-strand fill-in, which has been extensively studied in mammalian cells. However, we recently discovered that trypanosomes lack many telomere proteins identified in its mammalian host that are critical for telomere end processing. Rather, T. brucei uses a unique DNA polymerase, PolIE that belongs to the DNA polymerase A family (E. coli DNA PolI family), to coordinate the telomere G- and C-strand syntheses. In this review, I will first briefly summarize current understanding of telomere end processing in mammals. Subsequently, I will describe PolIE-mediated coordination of telomere G- and C-strand synthesis in T. brucei and implication of this recent discovery.

Vector-borne Trypanosoma brucei parasites develop in artificial human skin and persist as skin tissue forms

Transmission of Trypanosoma brucei by tsetse flies involves the deposition of the cell cycle-arrested metacyclic life cycle stage into mammalian skin at the site of the fly's bite. We introduce an advanced human skin equivalent and use tsetse flies to naturally infect the skin with trypanosomes. We detail the chronological order of the parasites' development in the skin by single-cell RNA sequencing and find a rapid activation of metacyclic trypanosomes and differentiation to proliferative parasites. Here we show that after the establishment of a proliferative population, the parasites enter a reversible quiescent state characterized by slow replication and a strongly reduced metabolism. We term these quiescent trypanosomes skin tissue forms, a parasite population that may play an important role in maintaining the infection over long time periods and in asymptomatic infected individuals.

Immunoprecipitation of RNA-DNA hybrid interacting proteins in Trypanosoma brucei reveals conserved and novel activities, including in the control of surface antigen expression needed for immune evasion by antigenic variation

RNA-DNA hybrids are epigenetic features of genomes that provide a diverse and growing range of activities. Understanding of these functions has been informed by characterising the proteins that interact with the hybrids, but all such analyses have so far focused on mammals, meaning it is unclear if a similar spectrum of RNA-DNA hybrid interactors is found in other eukaryotes. The African trypanosome is a single-cell eukaryotic parasite of the Discoba grouping and displays substantial divergence in several aspects of core biology from its mammalian host. Here, we show that DNA-RNA hybrid immunoprecipitation coupled with mass spectrometry recovers 602 putative interactors in T. brucei mammal- and insect-infective cells, some providing activities also found in mammals and some lineage-specific. We demonstrate that loss of three factors, two putative helicases and a RAD51 paralogue, alters T. brucei nuclear RNA-DNA hybrid and DNA damage levels. Moreover, loss of each factor affects the operation of the parasite immune survival mechanism of antigenic variation. Thus, our work reveals the broad range of activities contributed by RNA-DNA hybrids to T. brucei biology, including new functions in host immune evasion as well as activities likely fundamental to eukaryotic genome function.

Genome-wide chromatin interaction map for Trypanosoma cruzi

Trypanosomes are eukaryotic, unicellular parasites, such as Trypanosoma brucei, which causes sleeping sickness, and Trypanosoma cruzi, which causes Chagas disease. Genomes of these parasites comprise core regions and species-specific disruptive regions that encode multigene families of surface glycoproteins. Few transcriptional regulators have been identified in these parasites, and the role of spatial organization of the genome in gene expression is unclear. Here we mapped genome-wide chromatin interactions in T. cruzi using chromosome conformation capture (Hi-C), and we show that the core and disruptive regions form three-dimensional chromatin compartments named C and D. These chromatin compartments differ in levels of DNA methylation, nucleosome positioning and chromatin interactions, affecting genome expression dynamics. Our data reveal that the trypanosome genome is organized into chromatin-folding domains and transcription is affected by the local chromatin structure. We propose a model in which epigenetic mechanisms affect gene expression in trypanosomes.

Protein Phosphatase PP1 Regulation of Pol II Phosphorylation is Linked to Transcription Termination and Allelic Exclusion of VSG Genes and TERRA in Trypanosomes

The genomes of Leishmania and trypanosomes are organized into polycistronic transcription units flanked by a modified DNA base J involved in promoting RNA polymerase II (Pol II) termination. We recently characterized a Leishmania complex containing a J-binding protein, PP1 protein phosphatase 1, and PP1 regulatory protein (PNUTS) that controls transcription termination potentially via dephosphorylation of Pol II by PP1. While T. brucei contains eight PP1 isoforms, none purified with the PNUTS complex, suggesting a unique PP1-independent mechanism of termination. We now demonstrate that the PP1-binding motif of TbPNUTS is required for function in termination in vivo and that TbPP1-1 modulates Pol II termination in T. brucei involving dephosphorylation of the C-terminal domain of the large subunit of Pol II. PP1-1 knock-down results in increased cellular levels of phosphorylated large subunit of Pol II accompanied by readthrough transcription and pervasive transcription of the entire genome by Pol II, including Pol I transcribed loci that are typically silent, such as telomeric VSG expression sites involved in antigenic variation and production of TERRA RNA. These results provide important insights into the mechanism underlying Pol II transcription termination in primitive eukaryotes that rely on polycistronic transcription and maintain allelic exclusion of VSG genes.

Evolutionary, structural and functional insights in nuclear organisation and nucleocytoplasmic transport in trypanosomes

One of the remarkable features of eukaryotes is the nucleus, delimited by the nuclear envelope (NE), a complex structure and home to the nuclear lamina and nuclear pore complex (NPC). For decades, these structures were believed to be mainly architectural elements and, in the case of the NPC, simply facilitating nucleocytoplasmic trafficking. More recently, the critical roles of the lamina, NPC and other NE constituents in genome organisation, maintaining chromosomal domains and regulating gene expression have been recognised. Importantly, mutations in genes encoding lamina and NPC components lead to pathogenesis in humans, while pathogenic protozoa disrupt the progression of normal development and expression of pathogenesis-related genes. Here, we review features of the lamina and NPC across eukaryotes and discuss how these elements are structured in trypanosomes, protozoa of high medical and veterinary importance, highlighting lineage-specific and conserved aspects of nuclear organisation.

Histone divergence in trypanosomes results in unique alterations to nucleosome structure

Eukaryotes have a multitude of diverse mechanisms for organising and using their genomes, but the histones that make up chromatin are highly conserved. Unusually, histones from kinetoplastids are highly divergent. The structural and functional consequences of this variation are unknown. Here, we have biochemically and structurally characterised nucleosome core particles (NCPs) from the kinetoplastid parasite Trypanosoma brucei. A structure of the T. brucei NCP reveals that global histone architecture is conserved, but specific sequence alterations lead to distinct DNA and protein interaction interfaces. The T. brucei NCP is unstable and has weakened overall DNA binding. However, dramatic changes at the H2A-H2B interface introduce local reinforcement of DNA contacts. The T. brucei acidic patch has altered topology and is refractory to known binders, indicating that the nature of chromatin interactions in T. brucei may be unique. Overall, our results provide a detailed molecular basis for understanding evolutionary divergence in chromatin structure.

Decoding the impact of nuclear organization on antigenic variation in parasites

Antigenic variation as a strategy to evade the host adaptive immune response has evolved in divergent pathogens. Antigenic variation involves restricted, and often mutually exclusive, expression of dominant antigens and a periodic switch in antigen expression during infection. In eukaryotes, nuclear compartmentalization, including three-dimensional folding of the genome and physical separation of proteins in compartments or condensates, regulates mutually exclusive gene expression and chromosomal translocations. In this Review, we discuss the impact of nuclear organization on antigenic variation in the protozoan pathogens Trypanosoma brucei and Plasmodium falciparum. In particular, we highlight the relevance of nuclear organization in both mutually exclusive antigen expression and genome stability, which underlie antigenic variation.

Shifting landscapes: the role of 3D genomic organizations in gene regulatory strategies

3D genome folding enables the physical storage of chromosomes into the compact volume of a cell's nucleus, allows for the accurate segregation of chromatin to daughter cells, and has been shown to be tightly coupled to the way in which genetic information is converted into transcriptional programs [1-3]. Importantly, this link between chromatin architecture and gene regulation is a selectable feature in which modifications to chromatin organization accompany, or perhaps even drive the establishment of new regulatory strategies with enduring impacts on animal body plan complexity. Here, we discuss the nature of different 3D genome folding systems found across the tree of life, with particular emphasis on metazoans, and the relative influence of these systems on gene regulation. We suggest how the properties of these folding systems have influenced regulatory strategies employed by different lineages and may have catalyzed the partitioning and specialization of genetic programs that enabled multicellularity and organ-grade body plan complexity.

From contigs towards chromosomes: automatic improvement of long read assemblies (ILRA)

Recent advances in long read technologies not only enable large consortia to aim to sequence all eukaryotes on Earth, but they also allow individual laboratories to sequence their species of interest with relatively low investment. Long read technologies embody the promise of overcoming scaffolding problems associated with repeats and low complexity sequences, but the number of contigs often far exceeds the number of chromosomes and they may contain many insertion and deletion errors around homopolymer tracts. To overcome these issues, we have implemented the ILRA pipeline to correct long read-based assemblies. Contigs are first reordered, renamed, merged, circularized, or filtered if erroneous or contaminated. Illumina short reads are used subsequently to correct homopolymer errors. We successfully tested our approach by improving the genome sequences of Homo sapiens, Trypanosoma brucei, and Leptosphaeria spp., and by generating four novel Plasmodium falciparum assemblies from field samples. We found that correcting homopolymer tracts reduced the number of genes incorrectly annotated as pseudogenes, but an iterative approach seems to be required to correct more sequencing errors. In summary, we describe and benchmark the performance of our new tool, which improved the quality of novel long read assemblies up to 1 Gbp. The pipeline is available at GitHub: https://github.com/ThomasDOtto/ILRA.

How much (ATP) does it cost to build a trypanosome? A theoretical study on the quantity of ATP needed to maintain and duplicate a bloodstream-form Trypanosoma brucei cell

ATP hydrolysis is required for the synthesis, transport and polymerization of monomers for macromolecules as well as for the assembly of the latter into cellular structures. Other cellular processes not directly related to synthesis of biomass, such as maintenance of membrane potential and cellular shape, also require ATP. The unicellular flagellated parasite Trypanosoma brucei has a complex digenetic life cycle. The primary energy source for this parasite in its bloodstream form (BSF) is glucose, which is abundant in the host's bloodstream. Here, we made a detailed estimation of the energy budget during the BSF cell cycle. As glycolysis is the source of most produced ATP, we calculated that a single parasite produces 6.0 x 1011 molecules of ATP/cell cycle. Total biomass production (which involves biomass maintenance and duplication) accounts for ~63% of the total energy budget, while the total biomass duplication accounts for the remaining ~37% of the ATP consumption, with in both cases translation being the most expensive process. These values allowed us to estimate a theoretical YATP of 10.1 (g biomass)/mole ATP and a theoretical [Formula: see text] of 28.6 (g biomass)/mole ATP. Flagellar motility, variant surface glycoprotein recycling, transport and maintenance of transmembrane potential account for less than 30% of the consumed ATP. Finally, there is still ~5.5% available in the budget that is being used for other cellular processes of as yet unknown cost. These data put a new perspective on the assumptions about the relative energetic weight of the processes a BSF trypanosome undergoes during its cell cycle.

Profiling the bloodstream form and procyclic form Trypanosoma brucei cell cycle using single-cell transcriptomics

African trypanosomes proliferate as bloodstream forms (BSFs) and procyclic forms in the mammal and tsetse fly midgut, respectively. This allows them to colonise the host environment upon infection and ensure life cycle progression. Yet, understanding of the mechanisms that regulate and drive the cell replication cycle of these forms is limited. Using single-cell transcriptomics on unsynchronised cell populations, we have obtained high resolution cell cycle regulated (CCR) transcriptomes of both procyclic and slender BSF Trypanosoma brucei without prior cell sorting or synchronisation. Additionally, we describe an efficient freeze-thawing protocol that allows single-cell transcriptomic analysis of cryopreserved T. brucei. Computational reconstruction of the cell cycle using periodic pseudotime inference allowed the dynamic expression patterns of cycling genes to be profiled for both life cycle forms. Comparative analyses identify a core cycling transcriptome highly conserved between forms, as well as several genes where transcript levels dynamics are form specific. Comparing transcript expression patterns with protein abundance revealed that the majority of genes with periodic cycling transcript and protein levels exhibit a relative delay between peak transcript and protein expression. This work reveals novel detail of the CCR transcriptomes of both forms, which are available for further interrogation via an interactive webtool.

Drug resistance in Leishmania: does it really matter?

Treatment failure (TF) jeopardizes the management of parasitic diseases, including leishmaniasis. From the parasite's point of view, drug resistance (DR) is generally considered as central to TF. However, the link between TF and DR, as measured by in vitro drug susceptibility assays, is unclear, some studies revealing an association between treatment outcome and drug susceptibility, others not. Here we address three fundamental questions aiming to shed light on these ambiguities. First, are the right assays being used to measure DR? Second, are the parasites studied, which are generally those that adapt to in vitro culture, actually appropriate? Finally, are other parasite factors - such as the development of quiescent forms that are recalcitrant to drugs - responsible for TF without DR?

The RRM-mediated RNA binding activity in T. brucei RAP1 is essential for VSG monoallelic expression

Trypanosoma brucei is a protozoan parasite that causes human African trypanosomiasis. Its major surface antigen VSG is expressed from subtelomeric loci in a strictly monoallelic manner. We previously showed that the telomere protein TbRAP1 binds dsDNA through its 737RKRRR741 patch to silence VSGs globally. How TbRAP1 permits expression of the single active VSG is unknown. Through NMR structural analysis, we unexpectedly identify an RNA Recognition Motif (RRM) in TbRAP1, which is unprecedented for RAP1 homologs. Assisted by the 737RKRRR741 patch, TbRAP1 RRM recognizes consensus sequences of VSG 3'UTRs in vitro and binds the active VSG RNA in vivo. Mutating conserved RRM residues abolishes the RNA binding activity, significantly decreases the active VSG RNA level, and derepresses silent VSGs. The competition between TbRAP1's RNA and dsDNA binding activities suggests a VSG monoallelic expression mechanism in which the active VSG's abundant RNA antagonizes TbRAP1's silencing effect, thereby sustaining its full-level expression.

Targeting trypanosomes: how chemogenomics and artificial intelligence can guide drug discovery

Trypanosomatids are protozoan parasites that cause human and animal neglected diseases. Despite global efforts, effective treatments are still much needed. Phenotypic screens have provided several chemical leads for drug discovery, but the mechanism of action for many of these chemicals is currently unknown. Recently, chemogenomic screens assessing the susceptibility or resistance of parasites carrying genome-wide modifications started to define the mechanism of action of drugs at large scale. In this review, we discuss how genomics is being used for drug discovery in trypanosomatids, how integration of chemical and genomics data from these and other organisms has guided prioritisations of candidate therapeutic targets and additional chemical starting points, and how these data can fuel the expansion of drug discovery pipelines into the era of artificial intelligence.

Structural similarities between the metacyclic and bloodstream form variant surface glycoproteins of the African trypanosome

During infection of mammalian hosts, African trypanosomes thwart immunity using antigenic variation of the dense Variant Surface Glycoprotein (VSG) coat, accessing a large repertoire of several thousand genes and pseudogenes, and switching to antigenically distinct copies. The parasite is transferred to mammalian hosts by the tsetse fly. In the salivary glands of the fly, the pathogen adopts the metacyclic form and expresses a limited repertoire of VSG genes specific to that developmental stage. It has remained unknown whether the metacyclic VSGs possess distinct properties associated with this particular and discrete phase of the parasite life cycle. We present here three novel metacyclic form VSG N-terminal domain crystal structures (mVSG397, mVSG531, and mVSG1954) and show that they mirror closely in architecture, oligomerization, and surface diversity the known classes of bloodstream form VSGs. These data suggest that the mVSGs are unlikely to be a specialized subclass of VSG proteins, and thus could be poor candidates as the major components of prophylactic vaccines against trypanosomiasis.

The role of invariant surface glycoprotein 75 in xenobiotic acquisition by African trypanosomes

The surface proteins of parasitic protozoa mediate functions essential to survival within a host, including nutrient accumulation, environmental sensing and immune evasion. Several receptors involved in nutrient uptake and defence from the innate immune response have been described in African trypanosomes and, together with antigenic variation, contribute towards persistence within vertebrate hosts. Significantly, a superfamily of invariant surface glycoproteins (ISGs) populates the trypanosome surface, one of which, ISG75, is implicated in uptake of the century-old drug suramin. By CRISPR/Cas9 knockout and biophysical analysis, we show here that ISG75 directly binds suramin and mediates uptake of additional naphthol-related compounds, making ISG75 a conduit for entry of at least one structural class of trypanocidal compounds. However, ISG75 null cells present only modest attenuation of suramin sensitivity, have unaltered viability in vivo and in vitro and no alteration to suramin-invoked proteome responses. While ISG75 is demonstrated as a valid suramin cell entry pathway, we suggest the presence of additional mechanisms for suramin accumulation, further demonstrating the complexity of trypanosomatid drug interactions and potential for evolution of resistance.

VSGs Expressed during Natural T. b. gambiense Infection Exhibit Extensive Sequence Divergence and a Subspecies-Specific Bias towards Type B N-Terminal Domains

Trypanosoma brucei gambiense is the primary causative agent of human African trypanosomiasis (HAT), a vector-borne disease endemic to West and Central Africa. The extracellular parasite evades antibody recognition within the host bloodstream by altering its variant surface glycoprotein (VSG) coat through a process of antigenic variation. The serological tests that are widely used to screen for HAT use VSG as one of the target antigens. However, the VSGs expressed during human infection have not been characterized. Here, we use VSG sequencing (VSG-seq) to analyze the VSGs expressed in the blood of patients infected with T. b. gambiense and compared them to VSG expression in Trypanosoma brucei rhodesiense infections in humans as well as Trypanosoma brucei brucei infections in mice. The 44 VSGs expressed during T. b. gambiense infection revealed a striking bias toward expression of type B N termini (82% of detected VSGs). This bias is specific to T. b. gambiense, which is unique among T. brucei subspecies in its chronic clinical presentation and anthroponotic nature. The expressed T. b. gambiense VSGs also share very little similarity to sequences from 36 T. b. gambiense whole-genome sequencing data sets, particularly in areas of the VSG protein exposed to host antibodies, suggesting the antigen repertoire is under strong selective pressure to diversify. Overall, this work demonstrates new features of antigenic variation in T. brucei gambiense and highlights the importance of understanding VSG repertoires in nature. IMPORTANCE Human African trypanosomiasis is a neglected tropical disease primarily caused by the extracellular parasite Trypanosoma brucei gambiense. To avoid elimination by the host, these parasites repeatedly replace their variant surface glycoprotein (VSG) coat. Despite the important role of VSGs in prolonging infection, VSG expression during human infections is poorly understood. A better understanding of natural VSG gene expression dynamics can clarify the mechanisms that T. brucei uses to alter its VSG coat. We analyzed the expressed VSGs detected in the blood of patients with trypanosomiasis. Our findings indicate that there are features of antigenic variation unique to human-infective T. brucei subspecies and that natural VSG repertoires may vary more than previously expected.

The RNA-binding protein RBP33 dampens non-productive transcription in trypanosomes

In-depth analysis of the transcriptomes of several model organisms has revealed that genomes are pervasively transcribed, giving rise to an abundance of non-canonical and mainly antisense RNA polymerase II-derived transcripts that are produced from almost any genomic context. Pervasive RNAs are degraded by surveillance mechanisms, but the repertoire of proteins that control the fate of these non-productive transcripts is still incomplete. Trypanosomes are single-celled eukaryotes that show constitutive RNA polymerase II transcription and in which initiation and termination of transcription occur at a limited number of sites per chromosome. It is not known whether pervasive transcription exists in organisms with unregulated RNA polymerase II activity, and which factors could be involved in the process. We show here that depletion of RBP33 results in overexpression of ∼40% of all annotated genes in the genome, with a marked accumulation of sense and antisense transcripts derived from silenced regions. RBP33 loss does not result in a significant increase in chromatin accessibility. Finally, we have found that transcripts that increase in abundance upon RBP33 knockdown are significantly more stable in RBP33-depleted trypanosomes, and that the exosome complex is responsible for their degradation. Our results provide strong evidence that RBP33 dampens non-productive transcription in trypanosomes.

Uncovering Pseudogenes and Intergenic Protein-coding Sequences in TriTryps' Genomes

Trypanosomatids belong to a remarkable group of unicellular, parasitic organisms of the order Kinetoplastida, an early diverging branch of the phylogenetic tree of eukaryotes, exhibiting intriguing biological characteristics affecting gene expression (intronless polycistronic transcription, trans-splicing, and RNA editing), metabolism, surface molecules, and organelles (compartmentalization of glycolysis, variation of the surface molecules, and unique mitochondrial DNA), cell biology and life cycle (phagocytic vacuoles evasion and intricate patterns of cell morphogenesis). With numerous genomic-scale data of several trypanosomatids becoming available since 2005 (genomes, transcriptomes, and proteomes), the scientific community can further investigate the mechanisms underlying these unusual features and address other unexplored phenomena possibly revealing biological aspects of the early evolution of eukaryotes. One fundamental aspect comprises the processes and mechanisms involved in the acquisition and loss of genes throughout the evolutionary history of these primitive microorganisms. Here, we present a comprehensive in silico analysis of pseudogenes in three major representatives of this group: Leishmania major, Trypanosoma brucei, and Trypanosoma cruzi. Pseudogenes, DNA segments originating from altered genes that lost their original function, are genomic relics that can offer an essential record of the evolutionary history of functional genes, as well as clues about the dynamics and evolution of hosting genomes. Scanning these genomes with functional proteins as proxies to reveal intergenic regions with protein-coding features, relying on a customized threshold to distinguish statistically and biologically significant sequence similarities, and reassembling remnant sequences from their debris, we found thousands of pseudogenes and hundreds of open reading frames, with particular characteristics in each trypanosomatid: mutation profile, number, content, density, codon bias, average size, single- or multi-copy gene origin, number and type of mutations, putative primitive function, and transcriptional activity. These features suggest a common process of pseudogene formation, different patterns of pseudogene evolution and extant biological functions, and/or distinct genome organization undertaken by those parasites during evolution, as well as different evolutionary and/or selective pressures acting on distinct lineages.

Trypanosoma brucei histones are heavily modified with combinatorial post-translational modifications and mark Pol II transcription start regions with hyperacetylated H2A

Trypanosomes diverged from the main eukaryotic lineage about 600 million years ago, and display some unusual genomic and epigenetic properties that provide valuable insight into the early processes employed by eukaryotic ancestors to regulate chromatin-mediated functions. We analysed Trypanosoma brucei core histones by high mass accuracy middle-down mass spectrometry to map core histone post-translational modifications (PTMs) and elucidate cis-histone combinatorial PTMs (cPTMs). T. brucei histones are heavily modified and display intricate cPTMs patterns, with numerous hypermodified cPTMs that could contribute to the formation of non-repressive euchromatic states. The Trypanosoma brucei H2A C-terminal tail is hyperacetylated, containing up to five acetylated lysine residues. MNase-ChIP-seq revealed a striking enrichment of hyperacetylated H2A at Pol II transcription start regions, and showed that H2A histones that are hyperacetylated in different combinations localised to different genomic regions, suggesting distinct epigenetic functions. Our genomics and proteomics data provide insight into the complex epigenetic mechanisms used by this parasite to regulate a genome that lacks the transcriptional control mechanisms found in later-branched eukaryotes. The findings further demonstrate the complexity of epigenetic mechanisms that were probably shared with the last eukaryotic common ancestor.

Roles and interactions of the specialized initiation factors EIF4E2, EIF4E5 and EIF4E6 in Trypanosoma brucei: EIF4E2 maintains the abundances of S-phase mRNAs

Trypanosoma brucei has six versions of the cap-binding translation initiation factor EIF4E. We investigated the functions of EIF4E2, EIF4E3, EIF4E5 and EIF4E6 in bloodstream forms. We confirmed the protein associations previously found in procyclic forms, and detected specific co-purification of some RNA-binding proteins. Bloodstream forms lacking EIF4E5 grew normally and differentiated to replication-incompetent procyclic forms. Depletion of EIF4E6 inhibited bloodstream-form trypanosome growth and translation. EIF4E2 co-purified only the putative RNA binding protein SLBP2. Bloodstream forms lacking EIF4E2 multiplied slowly, had a low maximal cell density, and expressed the stumpy-form marker PAD1, but showed no evidence for enhanced stumpy-form signalling. EIF4E2 knock-out cells differentiated readily to replication-competent procyclic forms. EIF4E2 was strongly associated with a subset of mRNAs that are maximally abundant in S-phase, and these all had decreased abundances in EIF4E2 knock-out cells. Three EIF4E2 target mRNAs are also bound and stabilized by the Pumilio domain protein PUF9. Yeast 2-hybrid results suggested that PUF9 interacts directly with SLBP2, but PUF9 was not detected in EIF4E2 pull-downs. We speculate that the EIF4E2-SLBP2 complex might interact with its target mRNAs, perhaps via PUF9, only early during G1/S, stabilizing the mRNAs in preparation for translation later in S-phase or in early G2.

Invariant surface glycoprotein 65 of Trypanosoma brucei is a complement C3 receptor

African trypanosomes are extracellular pathogens of mammals and are exposed to the adaptive and innate immune systems. Trypanosomes evade the adaptive immune response through antigenic variation, but little is known about how they interact with components of the innate immune response, including complement. Here we demonstrate that an invariant surface glycoprotein, ISG65, is a receptor for complement component 3 (C3). We show how ISG65 binds to the thioester domain of C3b. We also show that C3 contributes to control of trypanosomes during early infection in a mouse model and provide evidence that ISG65 is involved in reducing trypanosome susceptibility to C3-mediated clearance. Deposition of C3b on pathogen surfaces, such as trypanosomes, is a central point in activation of the complement system. In ISG65, trypanosomes have evolved a C3 receptor which diminishes the downstream effects of C3 deposition on the control of infection.

Trypanosomes evade the immune response through antigenic variation of a surface coat containing variant surface glycoproteins (VSG). They also express invariant surface glycoproteins (ISGs), which are less well understood. Here, Macleod et al. show that ISG65 of T. brucei is a receptor for complement component 3. They provide the crystal structure of T. brucei ISG65 in complex with complement C3d and show evidence that ISG65 is involved in reducing trypanosome susceptibility to C3-mediated clearance in vivo.

Localization of Epigenetic Markers in Leishmania Chromatin

Eukaryotes use histone variants and post-translation modifications (PTMs), as well as DNA base modifications, to regulate DNA replication/repair, chromosome condensation, and gene expression. Despite the unusual organization of their protein-coding genes into large polycistronic transcription units (PTUs), trypanosomatid parasites also employ a “histone code” to control these processes, but the details of this epigenetic code are poorly understood. Here, we present the results of experiments designed to elucidate the distribution of histone variants and PTMs over the chromatin landscape of Leishmania tarentolae. These experiments show that two histone variants (H2A.Z and H2B.V) and three histone H3 PTMs (H3K4me3, H3K16ac, and H3K76me3) are enriched at transcription start sites (TSSs); while a histone variant (H3.V) and the trypanosomatid-specific hyper-modified DNA base J are located at transcription termination sites (TTSs). Reduced nucleosome density was observed at all TTSs and TSSs for RNA genes transcribed by RNA polymerases I (RNAPI) or RNAPIII; as well as (to a lesser extent) at TSSs for the PTUs transcribed by RNAPII. Several PTMs (H3K4me3, H3K16ac H3K20me2 and H3K36me3) and base J were enriched at centromeres, while H3K50ac was specifically associated with the periphery of these centromeric sequences. These findings significantly expand our knowledge of the epigenetic markers associated with transcription, DNA replication and/or chromosome segregation in these early diverging eukaryotes and will hopefully lay the groundwork for future studies to elucidate how they control these fundamental processes.

Stage-specific transcription activator ESB1 regulates monoallelic antigen expression in Trypanosoma brucei

Variant surface glycoprotein (VSG) coats bloodstream form Trypanosoma brucei parasites, and monoallelic VSG expression underpins the antigenic variation necessary for pathogenicity. One of thousands of VSG genes is transcribed by RNA polymerase I in a singular nuclear structure called the expression site body (ESB), but how monoallelic VSG transcription is achieved remains unclear. Using a localization screen of 153 proteins we found one, ESB-specific protein 1 (ESB1), that localized only to the ESB and is expressed only in VSG-expressing life cycle stages. ESB1 associates with DNA near the active VSG promoter and is necessary for VSG expression, with overexpression activating inactive VSG promoters. Mechanistically, ESB1 is necessary for recruitment of a subset of ESB components, including RNA polymerase I, revealing that the ESB has separately assembled subdomains. Because many trypanosomatid parasites have divergent ESB1 orthologues yet do not undergo antigenic variation, ESB1 probably represents an important class of transcription regulators.

A novel SNF2 ATPase complex in Trypanosoma brucei with a role in H2A.Z-mediated chromatin remodelling

A cascade of histone acetylation events with subsequent incorporation of a histone H2A variant plays an essential part in transcription regulation in various model organisms. A key player in this cascade is the chromatin remodelling complex SWR1, which replaces the canonical histone H2A with its variant H2A.Z. Transcriptional regulation of polycistronic transcription units in the unicellular parasite Trypanosoma brucei has been shown to be highly dependent on acetylation of H2A.Z, which is mediated by the histone-acetyltransferase HAT2. The chromatin remodelling complex which mediates H2A.Z incorporation is not known and an SWR1 orthologue in trypanosomes has not yet been reported. In this study, we identified and characterised an SWR1-like remodeller complex in T. brucei that is responsible for Pol II-dependent transcriptional regulation. Bioinformatic analysis of potential SNF2 DEAD/Box helicases, the key component of SWR1 complexes, identified a 1211 amino acids-long protein that exhibits key structural characteristics of the SWR1 subfamily. Systematic protein-protein interaction analysis revealed the existence of a novel complex exhibiting key features of an SWR1-like chromatin remodeller. RNAi-mediated depletion of the ATPase subunit of this complex resulted in a significant reduction of H2A.Z incorporation at transcription start sites and a subsequent decrease of steady-state mRNA levels. Furthermore, depletion of SWR1 and RNA-polymerase II (Pol II) caused massive chromatin condensation. The potential function of several proteins associated with the SWR1-like complex and with HAT2, the key factor of H2A.Z incorporation, is discussed.

Trypanosoma brucei is the causative agent of African trypanosomiasis (sleeping sickness) in humans and nagana in cattle. Its unusual genomic organisation featuring large polycistronic units requires a general mechanism of transcription initiation, because individual gene promoters are mostly absent. Despite the fact that the histone variant H2A.Z has previously been identified as a key player of transcription regulation, the complex responsible for correct H2A.Z incorporation at transcription start sites (TSS) remains elusive. In other eukaryotes, SWR1, a SNF2 ATPase-associated chromatin remodelling complex, is responsible for correct incorporation of this histone variant. This study identified a SWR1-like complex in T. brucei. Depletion of the SNF2 ATPase resulted in a reduction of H2A.Z incorporation at the TSS and decreased steady-state mRNA levels accompanied by chromatin condensation. In addition to the SWR1-like complex, we also identified a trypanosome-specific HAT2 complex that includes the histone acetyltransferases HAT2, a key player in the H2A.Z incorporation process. This complex has a trypanosome-specific composition that is different from the NuA4/TIP60 complex in Saccharomyces cerevisiae.

Transcriptomic profiling of Trypanosoma congolense mouthpart parasites from naturally infected flies

Background

Animal African trypanosomiasis, or nagana, is a veterinary disease caused by African trypanosomes transmitted by tsetse flies. In Africa, Trypanosoma congolense is one of the most pathogenic and prevalent causes of nagana in livestock, resulting in high animal morbidity and mortality and extensive production losses. In the tsetse fly, parasites colonise the midgut and eventually reach the mouthparts, from where they can be transmitted as the fly feeds on vertebrate hosts such as cattle. Despite the extreme importance of mouthpart-form parasites for disease transmission, very few global expression profile studies have been conducted in these parasite forms.

Methods

Here, we collected tsetse flies from the Shimba Hills National Reserve, a wildlife area in southeast Kenya, diagnosed T. congolense infections, and sequenced the transcriptomes of the T. congolense parasites colonising the mouthparts of the flies.

Results

We found little correlation between mouthpart parasites from natural and experimental fly infections. Furthermore, we performed differential gene expression analysis between mouthpart and bloodstream parasite forms and identified several surface-expressed genes and 152 novel hypothetical proteins differentially expressed in mouthpart parasites. Finally, we profiled variant antigen expression and observed that a variant surface glycoprotein (VSG) transcript belonging to T. congolense phylotype 8 (i.e. TcIL3000.A.H_000381200), previously observed to be enriched in metacyclic transcriptomes, was present in all wild-caught mouthpart samples as well as bloodstream-form parasites, suggestive of constitutive expression.

Conclusion

Our study provides transcriptomes of trypanosome parasites from naturally infected tsetse flies and suggests that a phylotype 8 VSG gene is constitutively expressed in metacyclic- and bloodstream-form parasites at the population level.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05258-y.

Pedal to the Metal: Nuclear Splicing Bodies Turbo-Charge VSG mRNA Production in African Trypanosomes

The African trypanosome Trypanosoma brucei is a parasite of the mammalian bloodstream and tissues, where an antigenically variable Variant Surface Glycoprotein (VSG) coat protects it from immune attack. This dense layer comprised of ∼107 VSG proteins, makes VSG by far the most abundant mRNA (7-10% total) and protein (∼10% total) in the bloodstream form trypanosome. How can such prodigious amounts of VSG be produced from a single VSG gene? Extremely high levels of RNA polymerase I (Pol I) transcription of the active VSG provide part of the explanation. However, recent discoveries highlight the role of pre-mRNA processing, both in maintaining high levels of VSG transcription, as well as its monoallelic expression. Trypanosome mRNAs are matured through trans-splicing a spliced leader (SL) RNA to the 5' end of precursor transcripts, meaning abundant SL RNA is required throughout the nucleus. However, requirement for SL RNA in the vicinity of the active VSG gene is so intense, that the cell reconfigures its chromatin architecture to facilitate interaction between the SL RNA genes and the active VSG. This presumably ensures that sufficient localised SL RNA is available, and not limiting for VSG mRNA expression. Recently, novel nuclear splicing bodies which appear to provide essential trans-splicing components, have been identified associating with the active VSG. These observations highlight the underappreciated role of pre-mRNA processing in modulating gene expression in trypanosomes. Dissecting the function of these nuclear RNA processing bodies should help us elucidate the mechanisms of both VSG expression and monoallelic exclusion in T. brucei.

VEX1 Influences mVSG Expression During the Transition to Mammalian Infectivity in Trypanosoma brucei

The Trypanosoma (T) brucei life cycle alternates between the tsetse fly vector and the mammalian host. In the insect, T. brucei undergoes several developmental stages until it reaches the salivary gland and differentiates into the metacyclic form, which is capable of infecting the next mammalian host. Mammalian infectivity is dependent on expression of the metacyclic variant surface glycoprotein genes as the cells develop into mature metacyclics. The VEX complex is essential for monoallelic variant surface glycoprotein expression in T. brucei bloodstream form, however, initiation of expression of the surface proteins genes during metacyclic differentiation is poorly understood. To better understand the transition to mature metacyclics and the control of metacyclic variant surface glycoprotein expression we examined the role of VEX1 in this process. We show that modulating VEX1 expression leads to a dysregulation of variant surface glycoprotein expression during metacyclogenesis, and that following both in vivo and in vitro metacyclic differentiation VEX1 relocalises from multiple nuclear foci in procyclic cells to one to two distinct nuclear foci in metacyclic cells - a pattern like the one seen in mammalian infective bloodstream forms. Our data suggest a role for VEX1 in the metacyclic differentiation process and their capacity to become infectious to the mammalian host.

Several different sequences are implicated in bloodstream-form-specific gene expression in Trypanosoma brucei

The parasite Trypanosoma brucei grows as bloodstream forms in mammalian hosts, and as procyclic forms in tsetse flies. In trypanosomes, gene expression regulation depends heavily on post-transcriptional mechanisms. Both the RNA-binding protein RBP10 and glycosomal phosphoglycerate kinase PGKC are expressed only in mammalian-infective forms. RBP10 targets procyclic-specific mRNAs for destruction, while PGKC is required for bloodstream-form glycolysis. Developmental regulation of both is essential: expression of either RBP10 or PGKC in procyclic forms inhibits their proliferation. We show that the 3’-untranslated region of the RBP10 mRNA is extraordinarily long—7.3kb—and were able to identify six different sequences, scattered across the untranslated region, which can independently cause bloodstream-form-specific expression. The 3’-untranslated region of the PGKC mRNA, although much shorter, still contains two different regions, of 125 and 153nt, that independently gave developmental regulation. No short consensus sequences were identified that were enriched either within these regulatory regions, or when compared with other mRNAs with similar regulation, suggesting that more than one regulatory RNA-binding protein is important for repression of mRNAs in procyclic forms. We also identified regions, including an AU repeat, that increased expression in bloodstream forms, or suppressed it in both forms. Trypanosome mRNAs that encode RNA-binding proteins often have extremely extended 3’-untranslated regions. We suggest that one function of this might be to act as a fail-safe mechanism to ensure correct regulation even if mRNA processing or expression of trans regulators is defective.

The parasite Trypanosoma brucei causes sleeping sickness in humans, and nagana in cattle, and is transmitted by Tsetse flies. It grows in the bloodstream and tissue fluids of mammalian hosts, as "bloodstream forms", and as "procyclic forms" in the midgut of tsetse flies. Several hundred proteins are expressed in a stage-specific fashion, and this is essential for parasite survival in the different environments. RBP10 is an RNA-binding protein that is expressed only in bloodstream forms. It binds to procyclic-specific mRNAs, and causes their destruction. PGKC is an enzyme that is also specifically expressed in bloodstream forms. Developmental regulation of both is essential: expression of either RBP10 or PGKC in procyclic forms prevents their growth. The mRNAs encoding both proteins are very unstable in procyclic forms, and the sequences responsible are in an "untranslated region" of the mRNA—sequences that follow the part that codes for protein. We here show that the mRNA encoding PGKC has two regions that independently cause developmental regulation, and that the very long untranslated region of the RBP10 mRNA has no fewer than six regulatory regions, but there were no obvious similarities between them. We suggest that the presence of several different regulatory sequences in trypanosome mRNAs might be a fail-safe mechanism to ensure correct regulation.

Shared Mechanisms for Mutually Exclusive Expression and Antigenic Variation by Protozoan Parasites

Cellular decision-making at the level of gene expression is a key process in the development and evolution of every organism. Variations in gene expression can lead to phenotypic diversity and the development of subpopulations with adaptive advantages. A prime example is the mutually exclusive activation of a single gene from within a multicopy gene family. In mammals, this ranges from the activation of one of the two immunoglobulin (Ig) alleles to the choice in olfactory sensory neurons of a single odorant receptor (OR) gene from a family of more than 1,000. Similarly, in parasites like Trypanosoma brucei, Giardia lamblia or Plasmodium falciparum, the process of antigenic variation required to escape recognition by the host immune system involves the monoallelic expression of vsg, vsp or var genes, respectively. Despite the importance of this process, understanding how this choice is made remains an enigma. The development of powerful techniques such as single cell RNA-seq and Hi-C has provided new insights into the mechanisms these different systems employ to achieve monoallelic gene expression. Studies utilizing these techniques have shown how the complex interplay between nuclear architecture, physical interactions between chromosomes and different chromatin states lead to single allele expression. Additionally, in several instances it has been observed that high-level expression of a single gene is preceded by a transient state where multiple genes are expressed at a low level. In this review, we will describe and compare the different strategies that organisms have evolved to choose one gene from within a large family and how parasites employ this strategy to ensure survival within their hosts.

Transcription Dependent Loss of an Ectopically Expressed Variant Surface Glycoprotein during Antigenic Variation in Trypanosoma brucei

In the mammalian host, Trypanosoma brucei is coated in a single-variant surface glycoprotein (VSG) species. Stochastic switching of the expressed VSG allows the parasite to escape detection by the host immune system. DNA double-strand breaks (DSB) trigger VSG switching, and repair via gene conversion results in an antigenically distinct VSG being expressed from the single active bloodstream-form expression site (BES). The single active BES is marked by VSG exclusion 2 (VEX2) protein. Here, we have disrupted monoallelic VSG expression by stably expressing a second telomeric VSG from a ribosomal locus. We found that cells expressing two VSGs contained one VEX2 focus that was significantly larger in size than the wild-type cells; this therefore suggests the ectopic VSG is expressed from the same nuclear position as the active BES. Unexpectedly, we report that in the double VSG-expressing cells, the DNA sequence of the ectopic copy is lost following a DSB in the active BES, despite it being spatially separated in the genome. The loss of the ectopic VSG is dependent on active transcription and does not disrupt the number or variety of templates used to repair a BES DSB and elicit a VSG switch. We propose that there are stringent mechanisms within the cell to reinforce monoallelic expression during antigenic variation. IMPORTANCE The single-cell parasite Trypanosoma brucei causes the fatal disease human African trypanosomiasis and is able to colonize the blood, fat, skin, and central nervous system. Trypanosomes survive in the mammalian host owing to a dense protective protein coat that consists of a single-variant surface glycoprotein species. Stochastic switching of one VSG for an immunologically distinct one enables the parasite to escape recognition by the host immune system. We have disrupted monoallelic antigen expression by expressing a second VSG and report that following DSB-triggered VSG switching, the DNA sequence of the ectopic VSG is lost in a transcription-dependent manner. We propose that there are strict requirements to ensure that only one variant antigen is expressed following a VSG switch, which has important implications for understanding how the parasite survives in the mammalian host.

Signatures of hybridization in Trypanosoma brucei

Genetic exchange among disease-causing micro-organisms can generate progeny that combine different pathogenic traits. Though sexual reproduction has been described in trypanosomes, its impact on the epidemiology of Human African Trypanosomiasis (HAT) remains controversial. However, human infective and non-human infective strains of Trypanosoma brucei circulate in the same transmission cycles in HAT endemic areas in subsaharan Africa, providing the opportunity for mating during the developmental cycle in the tsetse fly vector. Here we investigated inheritance among progeny from a laboratory cross of T. brucei and then applied these insights to genomic analysis of field-collected isolates to identify signatures of past genetic exchange. Genomes of two parental and four hybrid progeny clones with a range of DNA contents were assembled and analysed by k-mer and single nucleotide polymorphism (SNP) frequencies to determine heterozygosity and chromosomal inheritance. Variant surface glycoprotein (VSG) genes and kinetoplast (mitochondrial) DNA maxi- and minicircles were extracted from each genome to examine how each of these components was inherited in the hybrid progeny. The same bioinformatic approaches were applied to an additional 37 genomes representing the diversity of T. brucei in subsaharan Africa and T. evansi. SNP analysis provided evidence of crossover events affecting all 11 pairs of megabase chromosomes and demonstrated that polyploid hybrids were formed post-meiotically and not by fusion of the parental diploid cells. VSGs and kinetoplast DNA minicircles were inherited biparentally, with approximately equal numbers from each parent, whereas maxicircles were inherited uniparentally. Extrapolation of these findings to field isolates allowed us to distinguish clonal descent from hybridization by comparing maxicircle genotype to VSG and minicircle repertoires. Discordance between maxicircle genotype and VSG and minicircle repertoires indicated inter-lineage hybridization. Significantly, some of the hybridization events we identified involved human infective and non-human infective trypanosomes circulating in the same geographic areas.

Sexual reproduction allows genes from different individuals to be mixed up in the offspring. This is particularly important for disease-causing microbes, because new combinations of harmful traits can arise, potentially leading to more severe outbreaks of disease. Tsetse-transmitted trypanosomes are single-celled parasites that cause severe human and livestock diseases in tropical Africa. During their developmental cycle in the tsetse fly, trypanosomes can mate and produce hybrid trypanosomes, which have one set of chromosomes from each parent. But polyploid hybrids, with more than one set of chromosomes from one or both parents, are often observed too. Here we have investigated how these polyploid hybrids are formed by comparing the genomes of hybrid progeny with those of their parents. Analysis of the large, paired chromosomes of both diploid and polyploid hybrids showed frequent crossovers, which are the hallmark of meiosis, the special form of division that produces haploid gametes. This indicates that the polyploids were formed after meiosis rather than by fusion of the parental diploid cells. We also investigated the inheritance of two other features of trypanosomes: the large family of variant surface glycoprotein (VSG) genes, and the mitochondrial (kinetoplast) DNA. Hybrid clones had inherited about half the VSG genes from each parent, and also showed biparental inheritance of one component of the kinetoplast DNA, the minicircles. We assessed the relatedness of field-collected trypanosomes by comparing their VSG and minicircle repertoires, together with maxicircle genotype. While most isolates shared few VSGs or minicircles, a group of mostly human-infective strains from Uganda had a large proportion of their repertoires in common. Most of these trypanosomes were probably related by clonal descent, but we also identified that some were hybrids by the mismatch between their maxicircle genotype and their VSG and minicircle repertoires. These signals of hybridization were also detected in some of the other field-collected isolates, suggesting that genetic exchange is widespread in nature. Significantly, the hybridization events involved human infective and non-human infective trypanosomes circulating in the same geographic areas, providing a mechanism for the generation of new, potentially more pathogenic, trypanosome strains causing human disease.

H2B.V demarcates divergent strand-switch regions, some tDNA loci, and genome compartments in Trypanosoma cruzi and affects parasite differentiation and host cell invasion

Histone variants play a crucial role in chromatin structure organization and gene expression. Trypanosomatids have an unusual H2B variant (H2B.V) that is known to dimerize with the variant H2A.Z generating unstable nucleosomes. Previously, we found that H2B.V protein is enriched in tissue-derived trypomastigote (TCT) life forms, a nonreplicative stage of Trypanosoma cruzi, suggesting that this variant may contribute to the differences in chromatin structure and global transcription rates observed among parasite life forms. Here, we performed the first genome-wide profiling of histone localization in T. cruzi using epimastigotes and TCT life forms, and we found that H2B.V was preferentially located at the edges of divergent transcriptional strand switch regions, which encompass putative transcriptional start regions; at some tDNA loci; and between the conserved and disrupted genome compartments, mainly at trans-sialidase, mucin and MASP genes. Remarkably, the chromatin of TCT forms was depleted of H2B.V-enriched peaks in comparison to epimastigote forms. Interactome assays indicated that H2B.V associated specifically with H2A.Z, bromodomain factor 2, nucleolar proteins and a histone chaperone, among others. Parasites expressing reduced H2B.V levels were associated with higher rates of parasite differentiation and mammalian cell infectivity. Taken together, H2B.V demarcates critical genomic regions and associates with regulatory chromatin proteins, suggesting a scenario wherein local chromatin structures associated with parasite differentiation and invasion are regulated during the parasite life cycle.

Identification of sequence-specific promoters driving polycistronic transcription initiation by RNA polymerase II in trypanosomes

Protein-coding genes in trypanosomes occur in polycistronic transcription units (PTUs). How RNA polymerase II (Pol II) initiates transcription of PTUs has not been resolved; the current model favors chromatin modifications inducing transcription rather than sequence-specific promoters. Here, we uncover core promoters by functional characterization of Pol II peaks identified by chromatin immunoprecipitation sequencing (ChIP-seq). Two distinct promoters are located between divergent PTUs, each driving unidirectional transcription. Detailed analysis identifies a 75-bp promoter that is necessary and sufficient to drive full reporter expression and contains functional motifs. Analysis of further promoters suggests transcription initiation is regulated and promoters are either focused or dispersed. In contrast to the previous model of unregulated and promoter-independent transcription initiation, we find that sequence-specific promoters determine the initiation of Pol II transcription of protein-coding genes PTUs. These findings in Trypanosoma brucei suggest that in addition of chromatin modifications, promoter motifs-based regulation of gene expression is deeply conserved among eukaryotes.

May the Odds Be Ever in Your Favor: Non-deterministic Mechanisms Diversifying Cell Surface Molecule Expression

How does the information in the genome program the functions of the wide variety of cells in the body? While the development of biological organisms appears to follow an explicit set of genomic instructions to generate the same outcome each time, many biological mechanisms harness molecular noise to produce variable outcomes. Non-deterministic variation is frequently observed in the diversification of cell surface molecules that give cells their functional properties, and is observed across eukaryotic clades, from single-celled protozoans to mammals. This is particularly evident in immune systems, where random recombination produces millions of antibodies from only a few genes; in nervous systems, where stochastic mechanisms vary the sensory receptors and synaptic matching molecules produced by different neurons; and in microbial antigenic variation. These systems employ overlapping molecular strategies including allelic exclusion, gene silencing by constitutive heterochromatin, targeted double-strand breaks, and competition for limiting enhancers. Here, we describe and compare five stochastic molecular mechanisms that produce variety in pathogen coat proteins and in the cell surface receptors of animal immune and neuronal cells, with an emphasis on the utility of non-deterministic variation.