Structure analysis of two Toxoplasma gondii and Neospora caninum satellite DNA families and evolution of their common monomeric sequence

Mechanisms of adaptation and evolution in Toxoplasma gondii

Toxoplasma has high host flexibility, infecting all nucleated cells of mammals and birds. This implies that during its infective process the parasite must constantly adapt to different environmental situations, which in turn leads to modifications in its metabolism, regulation of gene transcription, translation of mRNAs and stage specific factors. There are conserved pathways that support these adaptations, which we aim to elucidate in this review. We begin by exploring the widespread epigenetic mechanisms and transcription regulators, continue with the supportive role of Heat Shock Proteins (Hsp), the translation regulation, stress granules, and finish with the emergence of contingency genes in highly variable genomic domains, such as subtelomeres. Within epigenetics, the discovery of a new histone variant of the H2B family (H2B.Z), contributing to T. gondii virulence and differentiation, but also gene expression regulation and its association with the metabolic state of the parasite, is highlighted. Associated with the regulation of gene expression are transcription factors (TFs). An overview of the main findings on TF and development is presented. We also emphasize the role of Hsp90 and Tgj1 in T. gondii metabolic fitness and the regulation of protein translation. Translation regulation is also highlighted as a mechanism for adaptation to conditions encountered by the parasite as well as stress granules containing mRNA and proteins generated in the extracellular tachyzoite. Another important aspect in evolution and adaptability are the subtelomeres because of their high variability and gene duplication rate. Toxoplasma possess multigene families of membrane proteins and contingency genes that are associated with different metabolic stresses. Among them parasite differentiation and environmental stresses stand out, including those that lead tachyzoite to bradyzoite conversion. Finally, we are interested in positioning protozoa as valuable evolution models, focusing on research related to the Extended Evolutionary Synthesis, based on models recently generated, such as extracellular adaptation and ex vivo cyst recrudescence.

Architecture, Chromatin and Gene Organization of Toxoplasma gondii Subtelomeres

Subtelomeres (ST) are chromosome regions that separate telomeres from euchromatin and play relevant roles in various biological processes of the cell. While their functions are conserved, ST structure and genetic compositions are unique to each species. This study aims to identify and characterize the subtelomeric regions of the 13 Toxoplasma gondii chromosomes of the Me49 strain. Here, STs were defined at chromosome ends based on poor gene density. The length of STs ranges from 8.1 to 232.4 kbp, with a gene density of 0.049 genes/kbp, lower than the Me49 genome (0.15 kbp). Chromatin organization showed that H3K9me3, H2A.X, and H3.3 are highly enriched near telomeres and the 5' end of silenced genes, decaying in intensity towards euchromatin. H3K4me3 and H2A.Z/H2B.Z are shown to be enriched in the 5' end of the ST genes. Satellite DNA was detected in almost all STs, mainly the sat350 family and a novel satellite named sat240. Beyond the STs, only short dispersed fragments of sat240 and sat350 were found. Within STs, there were 12 functional annotated genes, 59 with unknown functions (Hypothetical proteins), 15 from multigene FamB, and 13 from multigene family FamC. Some genes presented low interstrain synteny associated with the presence of satellite DNA. Orthologues of FamB and FamC were also detected in Neospora caninum and Hammondia hammondi. A re-analysis of previous transcriptomic data indicated that ST gene expression is strongly linked to the adaptation to different situations such as extracellular passage (evolve and resequencing study) and changes in metabolism (lack of acetyl-CoA cofactor). In conclusion, the ST region of the T. gondii chromosomes was defined, the STs genes were determined, and it was possible to associate them with high interstrain plasticity and a role in the adaptability of T. gondii to environmental changes.

Third-generation sequencing revises the molecular karyotype for Toxoplasma gondii and identifies emerging copy number variants in sexual recombinants

Toxoplasma gondii is a useful model for intracellular parasitism given its ease of culture in the laboratory and genomic resources. However, as for many other eukaryotes, the T. gondii genome contains hundreds of sequence gaps owing to repetitive and/or unclonable sequences that disrupt the assembly process. Here, we use the Oxford Nanopore Minion platform to generate near-complete de novo genome assemblies for multiple strains of T. gondii and its near relative, N. caninum. We significantly improved T. gondii genome contiguity (average N50 of ∼6.6 Mb) and added ∼2 Mb of newly assembled sequence. For all of the T. gondii strains that we sequenced (RH, ME49, CTG, II×III progeny clones CL13, S27, S21, S26, and D3X1), the largest contig ranged in size between 11.9 and 12.1 Mb in size, which is larger than any previously reported T. gondii chromosome, and found to be due to a consistent fusion of Chromosomes VIIb and VIII. These data were validated by mapping existing T. gondii ME49 Hi-C data to our assembly, providing parallel lines of evidence that the T. gondii karyotype consists of 13, rather than 14, chromosomes. By using this technology, we also resolved hundreds of tandem repeats of varying lengths, including in well-known host-targeting effector loci like rhoptry protein 5 (ROP5) and ROP38. Finally, when we compared T. gondii with N. caninum, we found that although the 13-chromosome karyotype was conserved, extensive, previously unappreciated chromosome-scale rearrangements had occurred in T. gondii and N. caninum since their most recent common ancestry.

Characterization of Toxoplasma gondii subtelomeric-like regions: identification of a long-range compositional bias that is also associated with gene-poor regions

BACKGROUND

Chromosome ends are composed of telomeric repeats and subtelomeric regions, which are patchworks of genes interspersed with repeated elements. Although chromosome ends display similar arrangements in different species, their sequences are highly divergent. In addition, these regions display a particular nucleosomal composition and bind specific factors, therefore producing a special kind of heterochromatin. Using data from currently available draft genomes we have characterized these putative Telomeric Associated Sequences in Toxoplasma gondii.

RESULTS

An all-vs-all pairwise comparison of T. gondii assembled chromosomes revealed the presence of conserved regions of ∼ 30 Kb located near the ends of 9 of the 14 chromosomes of the genome of the ME49 strain. Sequence similarity among these regions is ∼ 70%, and they are also highly conserved in the GT1 and VEG strains. However, they are unique to Toxoplasma with no detectable similarity in other Apicomplexan parasites. The internal structure of these sequences consists of 3 repetitive regions separated by high-complexity sequences without annotated genes, except for a gene from the Toxoplasma Specific Family. ChIP-qPCR experiments showed that nucleosomes associated to these sequences are enriched in histone H4 monomethylated at K20 (H4K20me1), and the histone variant H2A.X, suggesting that they are silenced sequences (heterochromatin). A detailed characterization of the base composition of these sequences, led us to identify a strong long-range compositional bias, which was similar to that observed in other genomic silenced fragments such as those containing centromeric sequences, and was negatively correlated to gene density.

CONCLUSIONS

We identified and characterized a region present in most Toxoplasma assembled chromosomes. Based on their location, sequence features, and nucleosomal markers we propose that these might be part of subtelomeric regions of T. gondii. The identified regions display a unique trinucleotide compositional bias, which is shared (despite the lack of any detectable sequence similarity) with other silenced sequences, such as those making up the chromosome centromeres. We also identified other genomic regions with this compositional bias (but no detectable sequence similarity) that might be functionally similar.

A complex small RNA repertoire is generated by a plant/fungal-like machinery and effected by a metazoan-like Argonaute in the single-cell human parasite Toxoplasma gondii

In RNA silencing, small RNAs produced by the RNase-III Dicer guide Argonaute-like proteins as part of RNA-induced silencing complexes (RISC) to regulate gene expression transcriptionally or post-transcriptionally. Here, we have characterized the RNA silencing machinery and exhaustive small RNAome of Toxoplasma gondii, member of the Apicomplexa, a phylum of animal- and human-infecting parasites that cause extensive health and economic damages to human populations worldwide. Remarkably, the small RNA-generating machinery of Toxoplasma is phylogenetically and functionally related to that of plants and fungi, and accounts for an exceptionally diverse array of small RNAs. This array includes conspicuous populations of repeat-associated small interfering RNA (siRNA), which, as in plants, likely generate and maintain heterochromatin at DNA repeats and satellites. Toxoplasma small RNAs also include many microRNAs with clear metazoan-like features whose accumulation is sometimes extremely high and dynamic, an unexpected finding given that Toxoplasma is a unicellular protist. Both plant-like heterochromatic small RNAs and metazoan-like microRNAs bind to a single Argonaute protein, Tg-AGO. Toxoplasma miRNAs co-sediment with polyribosomes, and thus, are likely to act as translational regulators, consistent with the lack of catalytic residues in Tg-AGO. Mass spectrometric analyses of the Tg-AGO protein complex revealed a common set of virtually all known RISC components so far characterized in human and Drosophila, as well as novel proteins involved in RNA metabolism. In agreement with its loading with heterochromatic small RNAs, Tg-AGO also associates substoichiometrically with components of known chromatin-repressing complexes. Thus, a puzzling patchwork of silencing processor and effector proteins from plant, fungal and metazoan origin accounts for the production and action of an unsuspected variety of small RNAs in the single-cell parasite Toxoplasma and possibly in other apicomplexans. This study establishes Toxoplasma as a unique model system for studying the evolution and molecular mechanisms of RNA silencing among eukaryotes.

Toxoplasma gondii is an important human parasite that causes life-threatening diseases in developing fetuses and in immunocompromised individuals, especially AIDS and transplant patients. Curiously, the Toxoplasma genome is deprived of most of the basic transcription factors that regulate gene expression in other eukaryotic cells. Therefore, alternative strategies must exist to modulate the many phases of the Toxoplasma complex life cycle that includes invasion of several hosts. Here, we investigate one of these strategies, by studying the repertoire of Toxoplasma silencing small RNAs (sRNAs). In eukaryotes, most of these regulatory molecules, 20–30nt-long, are produced by members of the Dicer RNase-III family, and exert their various functions through ubiquitous proteins called Argonaute (Ago). The surprising diversity of the Toxoplasma sRNAome uncovered in our study is consistent with those molecules exerting key functions during the parasite's life cycle, including, possibly, during virulent infection. The study also unravels an unsuspected level of complexity in the origin and mechanisms of action of the factors that generate and affect Toxoplasma sRNA, prompting a re-evaluation of our current views on RNA silencing in eukaryotes.

Toxoplasma H2A variants reveal novel insights into nucleosome composition and functions for this histone family

Toxoplasma gondii is an obligate intracellular parasite. Toxoplasmosis is incurable because of its ability to differentiate from the rapidly replicating tachyzoite stage into a latent cyst form (bradyzoite stage). Gene regulation pertinent to Toxoplasma differentiation involves histone modification, but very little is known about the histone proteins in this early branching eukaryote. Here, we report the characterization of three H2A histones, variants H2AX and H2AZ, and a canonical H2A1. H2AZ is the minor parasite H2A member. H2A1 and H2AX both have an SQ motif, but only H2AX has a complete SQ(E/D)varphi (where varphi denotes a hydrophobic residue) known to be phosphorylated in response to DNA damage. We show that a novel H2B variant interacts with H2AZ and H2A1 but not with H2AX. Chromatin immunoprecipitation (ChIP) revealed that H2AZ and H2Bv are enriched at active genes while H2AX is enriched at repressed genes as well as the silent TgIRE repeat element. During DNA damage, we detected an increase in H2AX phosphorylation as well as increases in h2a1 and h2ax transcription. We found that expression of h2ax, but not h2a1 or h2az, increases in bradyzoites generated in vitro. Similar analysis performed on mature bradyzoites generated in vivo, which are arrested in G0, showed that h2az and h2ax are expressed but h2a1 is not, consistent with the idea that h2a1 is the canonical histone orthologue in the parasite. The increase of H2AX, which localizes to silenced areas during bradyzoite differentiation, is consistent with the quiescent nature of this stage of the life cycle. Our results indicate that the early-branching eukaryotic parasite Toxoplasma contains nucleosomes of novel composition, which is likely to impact multiple facets of parasite biology, including the clinically important process of bradyzoite differentiation.

SET8-mediated methylations of histone H4 lysine 20 mark silent heterochromatic domains in apicomplexan genomes

Posttranslational histone modifications modulate chromatin-templated processes in various biological systems. H4K20 methylation is considered to have an evolutionarily ancient role in DNA repair and genome integrity, while its function in heterochromatin function and gene expression is thought to have arisen later during evolution. Here, we identify and characterize H4K20 methylases of the Set8 family in Plasmodium and Toxoplasma, two medically important members of the protozoan phylum Apicomplexa. Remarkably, parasite Set8-related proteins display H4K20 mono-, di-, and trimethylase activities, in striking contrast to the monomethylase-restricted human Set8. Structurally, few residues forming the substrate-specific channel dictate enzyme methylation multiplicity. These enzymes are cell cycle regulated and focally enriched at pericentric and telomeric heterochromatin in both parasites. Collectively, our findings provide new insights into the evolution of Set8-mediated biochemical pathways, suggesting that the heterochromatic function of the marker is not restricted to metazoans. Thus, these lower eukaryotes have developed a diverse panel of biological stages through their high capacity to differentiate, and epigenetics only begins to emerge as a strong determinant of their biology.