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

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.

Genome-wide localization of histone variants in Toxoplasma gondii implicates variant exchange in stage-specific gene expression

BACKGROUND

Toxoplasma gondii is a protozoan parasite that differentiates from acute tachyzoite stages to latent bradyzoite forms in response to environmental cues that modify the epigenome. We studied the distribution of the histone variants CenH3, H3.3, H2A.X, H2A.Z and H2B.Z, by genome-wide chromatin immunoprecipitation to understand the role of variant histones in developmental transitions of T. gondii parasites.

RESULTS

H3.3 and H2A.X were detected in telomere and telomere associated sequences, whereas H3.3, H2A.X and CenH3 were enriched in centromeres. Histones H2A.Z and H2B.Z colocalize with the transcriptional activation mark H3K4me3 in promoter regions surrounding the nucleosome-free region upstream of the transcription start site. The H2B.Z/H2A.Z histone pair also localizes to the gene bodies of genes that are silent but poised for activation, including bradyzoite stage-specific genes. The majority of H2A.X and H2A.Z/H2B.Z loci do not overlap, consistent with variant histones demarcating specific functional regions of chromatin. The extent of enrichment of H2A.Z/H2B.Z (and H3.3 and H2A.X) within the entire gene (5'UTR and gene body) reflects the timing of gene expression during the cell cycle, suggesting that dynamic turnover of H2B.Z/H2A.Z occurs during the tachyzoite cell cycle. Thus, the distribution of the variant histone H2A.Z/H2B.Z dimer defines active and developmentally silenced regions of the T. gondii epigenome including genes that are poised for expression.

CONCLUSIONS

Histone variants mark functional regions of parasite genomes with the dynamic placement of the H2A.Z/H2B.Z dimer implicated as an evolutionarily conserved regulator of parasite and eukaryotic differentiation.

Cdk7 mediates RPB1-driven mRNA synthesis in Toxoplasma gondii

Cyclin-dependent kinase 7 in conjunction with CyclinH and Mat1 activates cell cycle CDKs and is a part of the general transcription factor TFIIH. Role of Cdk7 is well characterized in model eukaryotes however its relevance in protozoan parasites has not been investigated. This important regulator of key processes warrants closer examination particularly in this parasite given its unique cell cycle progression and flexible mode of replication. We report functional characterization of TgCdk7 and its partners TgCyclinH and TgMat1. Recombinant Cdk7 displays kinase activity upon binding its cyclin partner and this activity is further enhanced in presence of Mat1. The activated kinase phosphorylates C-terminal domain of TgRPB1 suggesting its role in parasite transcription. Therefore, the function of Cdk7 in CTD phosphorylation and RPB1 mediated transcription was investigated using Cdk7 inhibitor. Unphosphorylated CTD binds promoter DNA while phosphorylation by Cdk7 triggers its dissociation from DNA with implications for transcription initiation. Inhibition of Cdk7 in the parasite led to strong reduction in Serine 5 phosphorylation of TgRPB1-CTD at the promoters of constitutively expressed actin1 and sag1 genes with concomitant reduction of both nascent RNA synthesis and 5′-capped transcripts. Therefore, we provide compelling evidence for crucial role of TgCdk7 kinase activity in mRNA synthesis.

Horizontal Gene Transfer Contributes to the Evolution of Arthropod Herbivory

Within animals, evolutionary transition toward herbivory is severely limited by the hostile characteristics of plants. Arthropods have nonetheless counteracted many nutritional and defensive barriers imposed by plants and are currently considered as the most successful animal herbivores in terrestrial ecosystems. We gather a body of evidence showing that genomes of various plant feeding insects and mites possess genes whose presence can only be explained by horizontal gene transfer (HGT). HGT is the asexual transmission of genetic information between reproductively isolated species. Although HGT is known to have great adaptive significance in prokaryotes, its impact on eukaryotic evolution remains obscure. Here, we show that laterally transferred genes into arthropods underpin many adaptations to phytophagy, including efficient assimilation and detoxification of plant produced metabolites. Horizontally acquired genes and the traits they encode often functionally diversify within arthropod recipients, enabling the colonization of more host plant species and organs. We demonstrate that HGT can drive metazoan evolution by uncovering its prominent role in the adaptations of arthropods to exploit plants.

Local admixture of amplified and diversified secreted pathogenesis determinants shapes mosaic Toxoplasma gondii genomes

Toxoplasma gondii is among the most prevalent parasites worldwide, infecting many wild and domestic animals and causing zoonotic infections in humans. T. gondii differs substantially in its broad distribution from closely related parasites that typically have narrow, specialized host ranges. To elucidate the genetic basis for these differences, we compared the genomes of 62 globally distributed T. gondii isolates to several closely related coccidian parasites. Our findings reveal that tandem amplification and diversification of secretory pathogenesis determinants is the primary feature that distinguishes the closely related genomes of these biologically diverse parasites. We further show that the unusual population structure of T. gondii is characterized by clade-specific inheritance of large conserved haploblocks that are significantly enriched in tandemly clustered secretory pathogenesis determinants. The shared inheritance of these conserved haploblocks, which show a different ancestry than the genome as a whole, may thus influence transmission, host range and pathogenicity.

Canonical histone H2Ba and H2A.X dimerize in an opposite genomic localization to H2A.Z/H2B.Z dimers in Toxoplasma gondii

Histone H2Ba of Toxoplasma gondii was expressed as recombinant protein (rH2Ba) and used to generate antibody in mouse that is highly specific. Antibody recognizing rH2Ba detects a single band in tachyzoite lysate of the expected molecular weight (12kDa). By indirect immunofluorescence (IFA) in in vitro grown tachyzoites and bradyzoites, the signal was detected only in the parasite nucleus. The nucleosome composition of H2Ba was determined through co-immunoprecipitation assays. H2Ba was detected in the same immunocomplex as H2A.X, but not with H2A.Z. Through chromatin immunoprecipitation (ChIP) assays and qPCR, it was observed that H2Ba is preferentially located at promoters of inactive genes and silent regions, accompanying H2A.X and opposed to H2A.Z/H2B.Z dimers.