A central role for Plasmodium falciparum subtelomeric regions in spatial positioning and telomere length regulation

Chromatin structure and var2csa - a tango in regulation of var gene expression in the human malaria parasite, Plasmodium falciparum?

Over the last few decades, novel methods have been developed to study how chromosome positioning within the nucleus may play a role in gene regulation. Adaptation of these methods in the human malaria parasite, Plasmodium falciparum, has recently led to the discovery that the three-dimensional structure of chromatin within the nucleus may be critical in controlling expression of virulence genes (var genes). Recent work has implicated an unusual, highly conserved var gene called var2csa in contributing to coordinated transcriptional switching, however how this gene functions in this capacity is unknown. To further understand how var2csa influences var gene switching, targeted DNA double-strand breaks (DSBs) within the sub-telomeric region of chromosome 12 were used to delete the gene and the surrounding chromosomal region. To characterize the changes in chromatin architecture stemming from this deletion and how these changes could affect var gene expression, we used a combination of RNA-seq, Chip-seq and Hi-C to pinpoint epigenetic and chromatin structural modifications in regions of differential gene expression. We observed a net gain of interactions in sub-telomeric regions and internal var gene regions following var2csa knockout, indicating an increase of tightly controlled heterochromatin structures. Our results suggest that disruption of var2csa results not only in changes in var gene transcriptional regulation but also a significant tightening of heterochromatin clusters thereby disrupting coordinated activation of var genes throughout the genome. Altogether our result confirms a strong link between the var2csa locus, chromatin structure and var gene expression.

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.

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.

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.

Chromosome splitting of Plasmodium berghei using the CRISPR/Cas9 system

Spatial arrangement of chromosomes is responsible for gene expression in Plasmodium parasites. However, methods for rearranging chromosomes have not been established, which makes it difficult to investigate its role in detail. Here, we report a method for splitting chromosome in rodent malaria parasite by CRISPR/Cas9 system using fragments in which a telomere and a centromere were incorporated. The resultant split chromosomes segregated accurately into daughter parasites by the centromere. In addition, elongation of de novo telomeres were observed, indicating its proper function. Furthermore, chromosome splitting had no effect on development of parasites. Splitting of the chromosome is expected to alter its spatial arrangement, and our method will thus be useful for investigating its biological role related with gene expression.

Plasmodium falciparum GBP2 Is a Telomere-Associated Protein That Binds to G-Quadruplex DNA and RNA

In the early-diverging protozoan parasite Plasmodium, few telomere-binding proteins have been identified and several are unique. Plasmodium telomeres, like those of most eukaryotes, contain guanine-rich repeats that can form G-quadruplex structures. In model systems, quadruplex-binding drugs can disrupt telomere maintenance and some quadruplex-binding drugs are potent anti-plasmodial agents. Therefore, telomere-interacting and quadruplex-interacting proteins may offer new targets for anti-malarial therapy. Here, we report that P. falciparum GBP2 is such a protein. It was identified via 'Proteomics of Isolated Chromatin fragments', applied here for the first time in Plasmodium. In vitro, PfGBP2 binds specifically to G-rich telomere repeats in quadruplex form and it can also bind to G-rich RNA. In vivo, PfGBP2 partially colocalises with the known telomeric protein HP1 but is also found in the cytoplasm, probably due to its affinity for RNA. Consistently, its interactome includes numerous RNA-associated proteins. PfGBP2 is evidently a multifunctional DNA/RNA-binding factor in Plasmodium.

Telomere Roles in Fungal Genome Evolution and Adaptation

Telomeres form the ends of linear chromosomes and usually comprise protein complexes that bind to simple repeated sequence motifs that are added to the 3′ ends of DNA by the telomerase reverse transcriptase (TERT). One of the primary functions attributed to telomeres is to solve the “end-replication problem” which, if left unaddressed, would cause gradual, inexorable attrition of sequences from the chromosome ends and, eventually, loss of viability. Telomere-binding proteins also protect the chromosome from 5′ to 3′ exonuclease action, and disguise the chromosome ends from the double-strand break repair machinery whose illegitimate action potentially generates catastrophic chromosome aberrations. Telomeres are of special interest in the blast fungus, Pyricularia, because the adjacent regions are enriched in genes controlling interactions with host plants, and the chromosome ends show enhanced polymorphism and genetic instability. Previously, we showed that telomere instability in some P. oryzae strains is caused by novel retrotransposons (MoTeRs) that insert in telomere repeats, generating interstitial telomere sequences that drive frequent, break-induced rearrangements. Here, we sought to gain further insight on telomeric involvement in shaping Pyricularia genome architecture by characterizing sequence polymorphisms at chromosome ends, and surrounding internalized MoTeR loci (relics) and interstitial telomere repeats. This provided evidence that telomere dynamics have played historical, and likely ongoing, roles in shaping the Pyricularia genome. We further demonstrate that even telomeres lacking MoTeR insertions are poorly preserved, such that the telomere-adjacent sequences exhibit frequent presence/absence polymorphism, as well as exchanges with the genome interior. Using TERT knockout experiments, we characterized chromosomal responses to failed telomere maintenance which suggested that much of the MoTeR relic-/interstitial telomere-associated polymorphism could be driven by compromised telomere function. Finally, we describe three possible examples of a phenomenon known as “Adaptive Telomere Failure,” where spontaneous losses of telomere maintenance drive rapid accumulation of sequence polymorphism with possible adaptive advantages. Together, our data suggest that telomere maintenance is frequently compromised in Pyricularia but the chromosome alterations resulting from telomere failure are not as catastrophic as prior research would predict, and may, in fact, be potent drivers of adaptive polymorphism.

Telomere length dynamics in response to DNA damage in malaria parasites

Malaria remains a major cause of morbidity and mortality in the developing world. Recent work has implicated chromosome end stability and the repair of DNA breaks through telomere healing as potent drivers of variant antigen diversification, thus associating basic mechanisms for maintaining genome integrity with aspects of host-parasite interactions. Here we applied long-read sequencing technology to precisely examine the dynamics of telomere addition and chromosome end stabilization in response to double-strand breaks within subtelomeric regions. We observed that the process of telomere healing induces the initial synthesis of telomere repeats well in excess of the minimal number required for end stability. However, once stabilized, these newly created telomeres appear to function normally, eventually returning to a length nearing that of intact chromosome ends. These results parallel recent observations in humans, suggesting an evolutionarily conserved mechanism for chromosome end repair.

PfGBP2 is a novel G-quadruplex binding protein in Plasmodium falciparum

Guanine-quadruplexes (G4s) are non-canonical DNA structures that can regulate key biological processes such as transcription, replication and telomere maintenance in several organisms including eukaryotes, prokaryotes and viruses. Recent reports have identified the presence of G4s within the AT-rich genome of Plasmodium falciparum, the protozoan parasite causing malaria. In Plasmodium, potential G4-forming sequences (G4FS) are enriched in the telomeric and sub-telomeric regions of the genome where they are associated with telomere maintenance and recombination events within virulence genes. However, there is a little understanding about the biological role of G4s and G4-binding proteins. Here, we provide the first snapshot of G4-interactome in P. falciparum using DNA pull-down assay followed by LC-MS/MS. Interestingly, we identified ~24 potential G4-binding proteins (G4-BP) that bind to a stable G4FS (AP2_G4). Furthermore, we characterised the role of G-strand binding protein 2 (PfGBP2), a putative telomere-binding protein in P. falciparum. We validated the interaction of PfGBP2 with G4 in vitro as well as in vivo. PfGBP2 is expressed throughout the intra-erythrocytic developmental cycle and is essential for the parasites in the presence of G4-stabilising ligand, pyridostatin. Gene knockout studies showed the role of PfGBP2 in the expression of var genes. Taken together, this study suggests that PfGBP2 is a bona fide G4-binding protein, which is likely to be involved in the regulation of G4-related functions in these malarial parasites. In addition, this study sheds light on this understudied G4 biology in P. falciparum.

Subtelomeric Transcription and its Regulation

The subtelomeres, highly heterogeneous repeated sequences neighboring telomeres, are transcribed into coding and noncoding RNAs in a variety of organisms. Telomereproximal subtelomeric regions produce non-coding transcripts i.e., ARRET, αARRET, subTERRA, and TERRA, which function in telomere maintenance. The role and molecular mechanisms of the majority of subtelomeric transcripts remain unknown. This review depicts the current knowledge and puts into perspective the results obtained in different models from yeasts to humans.

Genetic polymorphism of histidine rich protein 2 in Plasmodium falciparum isolates from different infection sources in Yunnan Province, China

Background

Failed diagnoses of some falciparum malaria cases by RDTs are constantly reported in recent years. Plasmodium falciparum histidine-rich protein 2 (pfhpr2) gene deficiency has been found to be the major reason of RDTs failure in many countries. This article analysed the deletion of pfhpr2 gene of falciparum malaria cases isolated in Yunnan Province, China.

Methods

Blood samples from falciparum malaria cases diagnosed in Yunnan Province were collected. Plasmodium genomic DNA was extracted and the pfhrp2 gene exon2 region was amplified via nested PCR. The haplotype of the DNA sequence, the nucleic acid diversity index (PI) and expected heterozygosity (He) were analyzed. Count PfHRP2 amino acid peptide sequence repeat and its times, and predict the properties of PfHRP2 peptide chain reaction to RDTs testing.

Results

A total of 306 blood samples were collected, 84.9% (259/306) from which pfhrp2 PCR amplification products (gene exon2) were obtained, while the remaining 47 samples were false amplification. The length of the 250 DNA sequences ranged from 345 - 927 bp, with 151 haplotypes, with PI and He values of 0.169 and 0.983, respectively. The length of the PfHRP2 peptide chain translated from 250 DNA sequences ranged from 115 to 309 aa. All peptide chains had more than an amino acid codon deletion. All 250 PfHRP2 strands ended with a type 12 amino acid repeat, 98.0% (245/250) started with a type 1 repetition and 2.0% (5/250) with a type 2 repetition. The detection rate for type 2 duplicates was 100% (250/250). Prediction of RDT sensitivity of PfHRP2 peptide chains based on type 2 and type 7 repeats showed that 9.60% (24/250), 50.0% (125/250), 13.20% (33/250) and 27.20.5% (68/250) of the 250 peptide chains were very sensitive, sensitive, borderline and non-sensitive, respectively.

Conclusion

The diversified polymorphism of the pfhrp2 gene deletion from different infection sources in the Yunnan province are extremely complex. The cause of the failure of pfhrp2 exon2 amplification is still to be investigated. The results of this study appeal to Yunnan Province for a timely evaluation of the effectiveness and applicability of RDTs in the diagnosis of malaria.

Identification of a TRBD zinc finger-interacting protein in Giardia duodenalis and its regulation of telomerase

Background

Giardia duodenalis causes giardiasis, with diarrhea as the primary symptom. The trophozoite proliferation of this zoonotic parasite is mainly affected by telomerase, although the mechanism of telomerase regulation has not been thoroughly analyzed.

Methods

This study was performed to identify the telomerase RNA-binding domain (TRBD)-interacting protein in G. duodenalis and its regulation of telomerase. Interaction between TRBD and interacting proteins was verified via pulldown assays and co-immunoprecipitation (co-IP) techniques, and the subcellular localization of the protein interactions was determined in vivo via split SNAP-tag labeling. The hammerhead ribozyme was designed to deplete the mRNA of TRBD-interacting proteins.

Results

Using TRBD as bait, we identified zinc-finger domain (ZFD)-containing proteins and verified it via pulldown and co-IP experiments. Protein-protein interaction occurred in the nuclei of 293T cells and both nuclei of G. duodenalis. The hammerhead ribozyme depleted ZFD mRNA levels, which reduced the reproduction rate of G. duodenalis, telomerase activity and telomere length.

Conclusions

Our findings suggest that ZFD may regulate telomere function in G. duodenalis nuclei.

Chromatin architecture and virulence-related gene expression in eukaryotic microbial pathogens

A fundamental question in biology is to understand how appropriate transcriptional regulation and dense packaging of the genetic material within the eukaryotic nucleus are achieved. The exquisite gene expression control and other metabolic processes of DNA require a highly complex, multilayered, three-dimensional architecture of the chromatin and its specific compartmentalization within the nucleus. Some of these architectural and sub-nuclear positioning mechanisms have been extensively co-opted by eukaryotic pathogens to keep fine expression control and expansion of virulence-related gene families in Plasmodium falciparum, Trypanosoma brucei and Candida glabrata. For example non-linear interactions between distant cis-acting regions and the formation of chromatin loops are required for appropriate regulation of the expression of virulence-related multi-gene families encoding cell surface proteins. These gene families are located near the chromosome ends and tethered to the nuclear periphery. Consequently, only one or very few genes of the family are expressed at a time. These genes are involved in antigenic variation in parasites and the generation of subpopulations of cells with diverse antigenic proteins at the surface in some pathogenic fungi, making them highly efficient pathogens.

G-Quadruplex DNA Motifs in the Malaria Parasite Plasmodium falciparum and Their Potential as Novel Antimalarial Drug Targets

G-quadruplexes are DNA or RNA secondary structures that can be formed from guanine-rich nucleic acids. These four-stranded structures, composed of stacked quartets of guanine bases, can be highly stable and have been demonstrated to occur in vivo in the DNA of human cells and other systems, where they play important biological roles, influencing processes such as telomere maintenance, DNA replication and transcription, or, in the case of RNA G-quadruplexes, RNA translation and processing. We report for the first time that DNA G-quadruplexes can be detected in the nuclei of the malaria parasite Plasmodium falciparum, which has one of the most A/T-biased genomes sequenced and therefore possesses few guanine-rich sequences with the potential to form G-quadruplexes. We show that despite this paucity of putative G-quadruplex-forming sequences, P. falciparum parasites are sensitive to several G-quadruplex-stabilizing drugs, including quarfloxin, which previously reached phase 2 clinical trials as an anticancer drug. Quarfloxin has a rapid initial rate of kill and is active against ring stages as well as replicative stages of intraerythrocytic development. We show that several G-quadruplex-stabilizing drugs, including quarfloxin, can suppress the transcription of a G-quadruplex-containing reporter gene in P. falciparum but that quarfloxin does not appear to disrupt the transcription of rRNAs, which was proposed as its mode of action in both human cells and trypanosomes. These data suggest that quarfloxin has potential for repositioning as an antimalarial with a novel mode of action. Furthermore, G-quadruplex biology in P. falciparum may present a target for development of other new antimalarial drugs.

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

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

Transcriptional Regulation of Telomeric Expression Sites and Antigenic Variation in Trypanosomes

INTRODUCTION

Trypanosoma brucei uses antigenic variation to evade the host antibody clearance by periodically changing its Variant Surface Glycoprotein (VSGs) coat. T. brucei encode over 2,500 VSG genes and pseudogenes, however they transcribe only one VSG gene at time from one of the 20 telomeric Expression Sites (ESs). VSGs are transcribed in a monoallelic fashion by RNA polymerase I from an extranucleolar site named ES body. VSG antigenic switching occurs by transcriptional switching between telomeric ESs or by recombination of the VSG gene expressed. VSG expression is developmentally regulated and its transcription is controlled by epigenetic mechanisms and influenced by a telomere position effect.

CONCLUSION

Here, we discuss 1) the molecular basis underlying transcription of telomeric ESs and VSG antigenic switching; 2) the current knowledge of VSG monoallelic expression; 3) the role of inositol phosphate pathway in the regulation of VSG expression and switching; and 4) the developmental regulation of Pol I transcription of procyclin and VSG genes.

PfAP2Tel, harbouring a non-canonical DNA-binding AP2 domain, binds to Plasmodium falciparum telomeres

The telomeres of the malaria parasite Plasmodium falciparum are essential not only for chromosome end maintenance during blood stage development in humans but also to generate genetic diversity by facilitating homologous recombination of subtelomeric, multigene virulence families such as var and rifin. However, other than the telomerase PfTERT, proteins that act at P. falciparum telomeres are poorly characterised. To isolate components that bind to telomeres, we performed oligonucleotide pulldowns and electromobility shift assays with a telomeric DNA probe and identified a non-canonical member of the ApiAP2 family of transcription factors, PfAP2Tel (encoded by PF3D7_0622900), as a component of the P. falciparum telomere-binding protein complex. PfAP2Tel is expressed throughout the intra-erythrocytic life cycle and localises to the nuclear periphery, co-localising with telomeric clusters. Furthermore, EMSAs using the recombinant protein demonstrated direct binding of PfAP2Tel to telomeric repeats in vitro, while genome-wide chromatin immunoprecipitation followed by next generation sequencing corroborated the high specificity of this protein to telomeric ends of all 14 chromosomes in vivo. Taken together, our data describe a novel function for ApiAP2 proteins at chromosome ends and open new avenues to study the molecular machinery that regulates telomere function in P. falciparum.

Malaria parasites possess a telomere repeat-binding protein that shares ancestry with transcription factor IIIA

Telomere repeat-binding factors (TRFs) are essential components of the molecular machinery that regulates telomere function. TRFs are widely conserved across eukaryotes and bind duplex telomere repeats via a characteristic MYB-type domain. Here, we identified the telomere repeat-binding protein PfTRZ in the malaria parasite Plasmodium falciparum, a member of the Alveolate phylum for which TRFs have not been described so far. PfTRZ lacks an MYB domain and binds telomere repeats via a C₂H₂-type zinc finger domain instead. In vivo, PfTRZ binds with high specificity to the telomeric tract and to interstitial telomere repeats upstream of subtelomeric virulence genes. Conditional depletion experiments revealed that PfTRZ regulates telomere length homeostasis and is required for efficient cell cycle progression. Intriguingly, we found that PfTRZ also binds to and regulates the expression of 5S rDNA genes. Combined with detailed phylogenetic analyses, our findings identified PfTRZ as a remote functional homologue of the basic transcription factor TFIIIA, which acquired a new function in telomere maintenance early in the apicomplexan lineage. Our work sheds unexpected new light on the evolution of telomere repeat-binding proteins and paves the way for dissecting the presumably divergent mechanisms regulating telomere functionality in one of the most deadly human pathogens.

Activation and clustering of a Plasmodium falciparum var gene are affected by subtelomeric sequences

The Plasmodium falciparum var multigene family encodes the cytoadhesive, variant antigen PfEMP1. P. falciparum antigenic variation and cytoadhesion specificity are controlled by epigenetic switching between the single, or few, simultaneously expressed var genes. Most var genes are maintained in perinuclear clusters of heterochromatic telomeres. The active var gene(s) occupy a single, perinuclear var expression site. It is unresolved whether the var expression site forms in situ at a telomeric cluster or whether it is an extant compartment to which single chromosomes travel, thus controlling var switching. Here we show that transcription of a var gene did not require decreased colocalisation with clusters of telomeres, supporting var expression site formation in situ. However following recombination within adjacent subtelomeric sequences, the same var gene was persistently activated and did colocalise less with telomeric clusters. Thus, participation in stable, heterochromatic, telomere clusters and var switching are independent but are both affected by subtelomeric sequences. The var expression site colocalised with the euchromatic mark H3K27ac to a greater extent than it did with heterochromatic H3K9me3. H3K27ac was enriched within the active var gene promoter even when the var gene was transiently repressed in mature parasites and thus H3K27ac may contribute to var gene epigenetic memory.

Genome-wide regulatory dynamics of G-quadruplexes in human malaria parasite Plasmodium falciparum

The AT-rich genome of P. falciparum has uniquely localized G-rich stretches that have propensity to form G-quadruplexes. However, their global occurrence and potential biological roles in the parasite are poorly explored. Our genome-wide analysis revealed unique enrichment of quadruplexes in P. falciparum genome which was remarkably different from other Plasmodium species. A distinct predominance of quadruplexes was observed in nuclear and organellar genes that participate in antigenic variation, pathogenesis, DNA/RNA regulation, metabolic and protein quality control processes. Data also suggested association of quadruplexes with SNPs and DNA methylation. Furthermore, analysis of steady state mRNA (RNA-seq) and polysome-associated mRNA (Ribosome profiling) data revealed stage-specific differences in translational efficiency of quadruplex harboring genes. Taken together, our findings hint towards existence of regulatory dynamics associated with quadruplexes that may modulate translational efficiency of quadruplex harboring genes to provide survival advantage to the parasite against host immune response and antimalarial drug pressure.

Trans-acting GC-rich non-coding RNA at var expression site modulates gene counting in malaria parasite

Monoallelic expression of the var multigene family enables immune evasion of the malaria parasite Plasmodium falciparum in its human host. At a given time only a single member of the 60-member var gene family is expressed at a discrete perinuclear region called the 'var expression site'. However, the mechanism of var gene counting remains ill-defined. We hypothesize that activation factors associating specifically with the expression site play a key role in this process. Here, we investigate the role of a GC-rich non-coding RNA (ncRNA) gene family composed of 15 highly homologous members. GC-rich genes are positioned adjacent to var genes in chromosome-central gene clusters but are absent near subtelomeric var genes. Fluorescence in situ hybridization demonstrates that GC-rich ncRNA localizes to the perinuclear expression site of central and subtelomeric var genes in trans. Importantly, overexpression of distinct GC-rich ncRNA members disrupts the gene counting process at the single cell level and results in activation of a specific subset of var genes in distinct clones. We identify the first trans-acting factor targeted to the elusive perinuclear var expression site and open up new avenues to investigate ncRNA function in antigenic variation of malaria and other protozoan pathogens.

The nuclear envelope and gene organization in parasitic protozoa: Specializations associated with disease

The parasitic protozoa Trypanosoma brucei and Plasmodium falciparum are lethal human parasites that have developed elegant strategies of immune evasion by antigenic variation. Despite the vast evolutionary distance between the two taxa, both parasites employ strict monoallelic expression of their membrane proteins, variant surface glycoproteins in Trypanosomes and the var, rif and stevor genes in Plasmodium, in order to evade their host's immune system. Additionally, both telomeric location and epigenetic controls are prominent features of these membrane proteins. As such, telomeres, chromatin structure and nuclear organization all contribute to control of gene expression and immune evasion. Here, we discuss the importance of epigenetics and sub-nuclear context for the survival of these disease-causing parasites.

Complete telomere-to-telomere de novo assembly of the Plasmodium falciparum genome through long-read (>11 kb), single molecule, real-time sequencing

The application of next-generation sequencing to estimate genetic diversity of Plasmodium falciparum, the most lethal malaria parasite, has proved challenging due to the skewed AT-richness [∼80.6% (A + T)] of its genome and the lack of technology to assemble highly polymorphic subtelomeric regions that contain clonally variant, multigene virulence families (Ex: var and rifin). To address this, we performed amplification-free, single molecule, real-time sequencing of P. falciparum genomic DNA and generated reads of average length 12 kb, with 50% of the reads between 15.5 and 50 kb in length. Next, using the Hierarchical Genome Assembly Process, we assembled the P. falciparum genome de novo and successfully compiled all 14 nuclear chromosomes telomere-to-telomere. We also accurately resolved centromeres [∼90–99% (A + T)] and subtelomeric regions and identified large insertions and duplications that add extra var and rifin genes to the genome, along with smaller structural variants such as homopolymer tract expansions. Overall, we show that amplification-free, long-read sequencing combined with de novo assembly overcomes major challenges inherent to studying the P. falciparum genome. Indeed, this technology may not only identify the polymorphic and repetitive subtelomeric sequences of parasite populations from endemic areas but may also evaluate structural variation linked to virulence, drug resistance and disease transmission.

Recombination and Diversification of the Variant Antigen Encoding Genes in the Malaria Parasite Plasmodium falciparum

The most severe form of human malaria is caused by the protozoan parasite Plasmodium falciparum. These parasites invade and replicate within the circulating red blood cells of infected individuals leading to numerous disease manifestations, including severe anemia, altered circulation, and tissue inflammation. Malaria parasites are also known for their ability to maintain a chronic infection through antigenic variation, the ability to systematically alter the antigens displayed on the surface of infected cells and thereby avoid clearance by the host's antibody response. The genome of P. falciparum includes several large, multicopy gene families that encode highly variable forms of the surface proteins that are the targets of host immunity. Alterations in expression of genes within these families are responsible for antigenic variation. This process requires the continuous generation of new antigenic variants within these gene families, and studies have shown that new variants arise through extensive recombination and gene conversion events between family members. Malaria parasites possess an unusual complement of DNA repair pathways, thus the study of recombination between variant antigen encoding genes provides a unique view into the evolution of mobile DNA in an organism distantly related to the more closely studied model eukaryotes.

Subtelomere organization in the genome of the microsporidian Encephalitozoon cuniculi: patterns of repeated sequences and physicochemical signatures

Background

The microsporidian Encephalitozoon cuniculi is an obligate intracellular eukaryotic pathogen with a small nuclear genome (2.9 Mbp) consisting of 11 chromosomes. Although each chromosome end is known to contain a single rDNA unit, the incomplete assembly of subtelomeric regions following sequencing of the genome identified only 3 of the 22 expected rDNA units. While chromosome end assembly remains a difficult process in most eukaryotic genomes, it is of significant importance for pathogens because these regions encode factors important for virulence and host evasion.

Results

Here we report the first complete assembly of E. cuniculi chromosome ends, and describe a novel mosaic structure of segmental duplications (EXT repeats) in these regions. EXT repeats range in size between 3.5 and 23.8 kbp and contain four multigene families encoding membrane associated proteins. Twenty-one recombination sites were identified in the sub-terminal region of E. cuniculi chromosomes. Our analysis suggests that these sites contribute to the diversity of chromosome ends organization through Double Strand Break repair mechanisms. The region containing EXT repeats at chromosome extremities can be differentiated based on gene composition, GC content, recombination sites density and chromosome landscape.

Conclusion

Together this study provides the complete structure of the chromosome ends of E. cuniculi GB-M1, and identifies important factors, which could play a major role in parasite diversity and host-parasite interactions. Comparison with other eukaryotic genomes suggests that terminal regions could be distinguished precisely based on gene content, genetic instability and base composition biais. The diversity of processes assciated with chromosome extremities and their biological consequences, as they are presented in the present study, emphasize the fact that great effort will be necessary in the future to characterize more carefully these regions during whole genome sequencing efforts.

Electronic supplementary material

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

Malaria Parasite Proteins and Their Role in Alteration of the Structure and Function of Red Blood Cells

Malaria, caused by Plasmodium spp., continues to be a major threat to human health and a significant cause of socioeconomic hardship in many countries. Almost half of the world's population live in malaria-endemic regions and many of them suffer one or more, often life-threatening episodes of malaria every year, the symptoms of which are attributable to replication of the parasite within red blood cells (RBCs). In the case of Plasmodium falciparum, the species responsible for most malaria-related deaths, parasite replication within RBCs is accompanied by striking alterations to the morphological, biochemical and biophysical properties of the host cell that are essential for the parasites' survival. To achieve this, the parasite establishes a unique and extensive protein export network in the infected RBC, dedicating at least 6% of its genome to the process. Understanding the full gamut of proteins involved in this process and the mechanisms by which P. falciparum alters the structure and function of RBCs is important both for a more complete understanding of the pathogenesis of malaria and for development of new therapeutic strategies to prevent or treat this devastating disease. This review focuses on what is currently known about exported parasite proteins, their interactions with the RBC and their likely pathophysiological consequences.

Antigenic Variation in Plasmodium falciparum

Plasmodium falciparum is the protozoan parasite that causes most malaria-associated morbidity and mortality in humans with over 500,000 deaths annually. The disease symptoms are associated with repeated cycles of invasion and asexual multiplication inside red blood cells of the parasite. Partial, non-sterile immunity to P. falciparum malaria develops only after repeated infections and continuous exposure. The successful evasion of the human immune system relies on the large repertoire of antigenically diverse parasite proteins displayed on the red blood cell surface and on the merozoite membrane where they are exposed to the human immune system. Expression switching of these polymorphic proteins between asexual parasite generations provides an efficient mechanism to adapt to the changing environment in the host and to maintain chronic infection. This chapter discusses antigenic diversity and variation in the malaria parasite and our current understanding of the molecular mechanisms that direct the expression of these proteins.

Functional analysis of sirtuin genes in multiple Plasmodium falciparum strains

Plasmodium falciparum, the causative agent of severe human malaria, employs antigenic variation to avoid host immunity. Antigenic variation is achieved by transcriptional switching amongst polymorphic var genes, enforced by epigenetic modification of chromatin. The histone-modifying 'sirtuin' enzymes PfSir2a and PfSir2b have been implicated in this process. Disparate patterns of var expression have been reported in patient isolates as well as in cultured strains. We examined var expression in three commonly used laboratory strains (3D7, NF54 and FCR-3) in parallel. NF54 parasites express significantly lower levels of var genes compared to 3D7, despite the fact that 3D7 was originally a clone of the NF54 strain. To investigate whether this was linked to the expression of sirtuins, genetic disruption of both sirtuins was attempted in all three strains. No dramatic changes in var gene expression occurred in NF54 or FCR-3 following PfSir2b disruption, contrasting with previous observations in 3D7. In 3D7, complementation of the PfSir2a genetic disruption resulted in a significant decrease in previously-elevated var gene expression levels, but with the continued expression of multiple var genes. Finally, rearranged chromosomes were observed in the 3D7 PfSir2a knockout line. Our results focus on the potential for parasite genetic background to contribute to sirtuin function in regulating virulence gene expression and suggest a potential role for sirtuins in maintaining genome integrity.

Genome-wide patterns of genetic polymorphism and signatures of selection in Plasmodium vivax

Plasmodium vivax is the most prevalent human malaria parasite outside of Africa. Yet, studies aimed to identify genes with signatures consistent with natural selection are rare. Here, we present a comparative analysis of the pattern of genetic variation of five sequenced isolates of P. vivax and its divergence with two closely related species, Plasmodium cynomolgi and Plasmodium knowlesi, using a set of orthologous genes. In contrast to Plasmodium falciparum, the parasite that causes the most lethal form of human malaria, we did not find significant constraints on the evolution of synonymous sites genome wide in P. vivax. The comparative analysis of polymorphism and divergence across loci allowed us to identify 87 genes with patterns consistent with positive selection, including genes involved in the “exportome” of P. vivax, which are potentially involved in evasion of the host immune system. Nevertheless, we have found a pattern of polymorphism genome wide that is consistent with a significant amount of constraint on the replacement changes and prevalent negative selection. Our analyses also show that silent polymorphism tends to be larger toward the ends of the chromosomes, where many genes involved in antigenicity are located, suggesting that natural selection acts not only by shaping the patterns of variation within the genes but it also affects genome organization.

P. berghei telomerase subunit TERT is essential for parasite survival

Telomeres define the ends of chromosomes protecting eukaryotic cells from chromosome instability and eventual cell death. The complex regulation of telomeres involves various proteins including telomerase, which is a specialized ribonucleoprotein responsible for telomere maintenance. Telomeres of chromosomes of malaria parasites are kept at a constant length during blood stage proliferation. The 7-bp telomere repeat sequence is universal across different Plasmodium species (GGGTTT/CA), though the average telomere length varies. The catalytic subunit of telomerase, telomerase reverse transcriptase (TERT), is present in all sequenced Plasmodium species and is approximately three times larger than other eukaryotic TERTs. The Plasmodium RNA component of TERT has recently been identified in silico. A strategy to delete the gene encoding TERT via double cross-over (DXO) homologous recombination was undertaken to study the telomerase function in P. berghei. Expression of both TERT and the RNA component (TR) in P. berghei blood stages was analysed by Western blotting and Northern analysis. Average telomere length was measured in several Plasmodium species using Telomere Restriction Fragment (TRF) analysis. TERT and TR were detected in blood stages and an average telomere length of ∼950 bp established. Deletion of the tert gene was performed using standard transfection methodologies and we show the presence of tert⁻ mutants in the transfected parasite populations. Cloning of tert- mutants has been attempted multiple times without success. Thorough analysis of the transfected parasite populations and the parasite obtained from extensive parasite cloning from these populations provide evidence for a so called delayed death phenotype as observed in different organisms lacking TERT. The findings indicate that TERT is essential for P. berghei cell survival. The study extends our current knowledge on telomere biology in malaria parasites and validates further investigations to identify telomerase inhibitors to induce parasite cell death.

Noncoding RNAs as emerging regulators of Plasmodium falciparum virulence gene expression

The eukaryotic unicellular pathogen Plasmodium falciparum tightly regulates gene expression, both during development and in adaptation to dynamic host environments. This regulation is evident in the mutually exclusive expression of members of clonally variant virulence multigene families. While epigenetic regulators have been selectively identified at active or repressed virulence genes, their specific recruitment remains a mystery. In recent years, noncoding RNAs (ncRNAs) have emerged as lynchpins of eukaryotic gene regulation; by binding to epigenetic regulators, they provide target specificity to otherwise non-specific enzyme complexes. Not surprisingly, there is great interest in understanding the role of ncRNA in P. falciparum, in particular, their contribution to the mutually exclusive expression of virulence genes. The current repertoire of P. falciparum ncRNAs includes, but is not limited to, subtelomeric ncRNAs, virulence gene-associated ncRNAs and natural antisense RNA transcripts. Continued improvement in high-throughput sequencing methods is sure to expand this repertoire. Here, we summarize recent advances in P. falciparum ncRNA biology, with an emphasis on ncRNA-mediated epigenetic modes of gene regulation.

Analysis of pfhrp2 genetic diversity in Senegal and implications for use of rapid diagnostic tests

Background

The Senegalese National Malaria Control Programme has recommended use of rapid diagnostic tests (RDTs) that target the histidine-rich protein 2 (HRP2), specific to Plasmodium falciparum, to diagnose malaria cases. The target antigen has been shown to be polymorphic, which may explain the variability in HRP2-based RDT results reported in field studies. The genetic diversity of the pfhrp2 gene has not been investigated in depth in many African countries. The goal of this study is to determine the extent of polymorphism in pfhrp2 among Senegal, Mali and Uganda parasite populations, and discuss the implications of these findings on the utility of RDTs that are based on HRP2 detection.

Methods

Sequencing data from the pfhrp2 locus were used to analyze the genetic diversity of this gene among three populations, with different transmission dynamics and malaria parasite ecologies. Nucleotide diversity (π) and non-synonymous nucleotide diversity (π_NS) were studied in the pfhrp2 gene from isolates obtained in Senegal. Amino acid repeat length polymorphisms in the PfHRP2 antigen were characterized and parameters of genetic diversity, such as frequency and correlation between repeats in these populations, were assessed.

Results

The diversity survey of the pfhrp2 gene identified 29 SNPs as well as insertion and deletion polymorphisms within a 918 bp region. The Senegal pfhrp2 exhibited a substantial level of diversity [π = 0.00559 and π_NS = 0.014111 (π_S = 0.0291627)], similar to several polymorphic genes, such as msp1, involved in immune responses, and the gene encoding the SURFIN polymorphic antigen, which are surface exposed parasite proteins. Extensive repeat length polymorphisms in PfHRP2, as well as similar patterns in the number, organization and the type of predicted amino acid repeats were observed among the three populations, characterized by an occurrence of Type 2, Type 4 and Type 7 repeats.

Conclusions

These results warrant deeper monitoring of the RDT target antigen diversity and emphasize that development of other essential genes as a target for diagnostic tools is critical.

Epigenetic regulation of the Plasmodium falciparum genome

Recent research has highlighted some unique aspects of chromatin biology in the malaria parasite Plasmodium falciparum. During its erythrocytic lifecycle P. falciparum maintains its genome primarily as unstructured euchromatin. Indeed there is no clear role for chromatin-mediated silencing of the majority of the developmentally expressed genes in P. falciparum. However discontinuous stretches of heterochromatin are critical for variegated expression of contingency genes that mediate key pathogenic processes in malaria. These range from invasion of erythrocytes and antigenic variation to solute transport and growth adaptation in response to environmental changes. Despite lack of structure within euchromatin the nucleus maintains functional compartments that regulate expression of many genes at the nuclear periphery, particularly genes with clonally variant expression. The typical components of the chromatin regulatory machinery are present in P. falciparum; however, some of these appear to have evolved novel species-specific functions, e.g. the dynamic regulation of histone variants at virulence gene promoters. The parasite also appears to have repeatedly acquired chromatin regulatory proteins through lateral transfer from endosymbionts and from the host. P. falciparum chromatin regulators have been successfully targeted with multiple drugs in laboratory studies; hopefully their functional divergence from human counterparts will allow the development of parasite-specific inhibitors.

The Toxoplasma gondii kinetochore is required for centrosome association with the centrocone (spindle pole)

The kinetochore is a multi-protein structure assembled on eukaryotic centromeres mediating chromosome attachment to spindle microtubules. Here we identified the kinetochore proteins Nuf2 and Ndc80 in the apicomplexan parasite Toxoplasma gondii. Localization revealed that kinetochores remain clustered throughout the cell cycle and colocalize with clustered centromeres at the centrocone, a structure containing the spindle pole embedded in the nuclear envelope. Pharmacological disruption of microtubules resulted in partial loss of some kinetochore and centromere clustering, indicating microtubules are necessary but not strictly required for kinetochore clustering. Generation of a TgNuf2 conditional knock-down strain revealed it is essential for chromosome segregation, but dispensable for centromere clustering. The centromeres actually remained associated with the centrocone suggesting microtubule binding is not required for their interaction with the spindle pole. The most striking observation upon TgNuf2 depletion was that the centrosome behaved normally, but that it lost its association with the centrocone. This suggests that microtubules are essential to maintain contact between the centrosome and chromosomes, and this interaction is critical for the partitioning of the nuclei into the two daughter parasites. Finally, genetic complementation experiments with mutated TgNuf2 constructs highlighted an apicomplexan-specific motif with a putative role in nuclear localization.

Whole genome mapping and re-organization of the nuclear and mitochondrial genomes of Babesia microti isolates

Babesia microti is the primary causative agent of human babesiosis, an emerging pathogen that causes a malaria-like illness with possible fatal outcome in immunocompromised patients. The genome sequence of the B. microti R1 strain was reported in 2012 and revealed a distinct evolutionary path for this pathogen relative to that of other apicomplexa. Lacking from the first genome assembly and initial molecular analyses was information about the terminal ends of each chromosome, and both the exact number of chromosomes in the nuclear genome and the organization of the mitochondrial genome remained ambiguous. We have now performed various molecular analyses to characterize the nuclear and mitochondrial genomes of the B. microti R1 and Gray strains and generated high-resolution Whole Genome maps. These analyses show that the genome of B. microti consists of four nuclear chromosomes and a linear mitochondrial genome present in four different structural types. Furthermore, Whole Genome mapping allowed resolution of the chromosomal ends, identification of areas of misassembly in the R1 genome, and genomic differences between the R1 and Gray strains, which occur primarily in the telomeric regions. These studies set the stage for a better understanding of the evolution and diversity of this important human pathogen.

Nuclear pores and perinuclear expression sites of var and ribosomal DNA genes correspond to physically distinct regions in Plasmodium falciparum

The human malaria parasite Plasmodium falciparum modifies the erythrocyte it infects by exporting variant proteins to the host cell surface. The var gene family that codes for a large, variant adhesive surface protein called P. falciparum erythrocyte membrane protein 1 (PfEMP1) plays a particular role in this process, which is linked to pathogenesis and immune evasion. A single member of this gene family is highly transcribed while the other 59 members remain silenced. Importantly, var gene transcription occurs at a spatially restricted, but yet undefined, perinuclear site that is distinct from repressed var gene clusters. To advance our understanding of monoallelic expression, we investigated whether nuclear pores associate with the var gene expression site. To this end, we studied the nuclear pore organization during the asexual blood stage using a specific antibody directed against a subunit of the nuclear pore, P. falciparum Nup116 (PfNup116). Ring and schizont stage parasites showed highly polarized nuclear pore foci, whereas in trophozoite stage nuclear pores redistributed over the entire nuclear surface. Colocalization studies of var transcripts and anti-PfNup116 antibodies showed clear dissociation between nuclear pores and the var gene expression site in ring stage. Similar results were obtained for another differentially transcribed perinuclear gene family, the ribosomal DNA units. Furthermore, we show that in the poised state, the var gene locus is not physically linked to nuclear pores. Our results indicate that P. falciparum does form compartments of high transcriptional activity at the nuclear periphery which are, unlike the case in yeast, devoid of nuclear pores.

Mitotic evolution of Plasmodium falciparum shows a stable core genome but recombination in antigen families

Malaria parasites elude eradication attempts both within the human host and across nations. At the individual level, parasites evade the host immune responses through antigenic variation. At the global level, parasites escape drug pressure through single nucleotide variants and gene copy amplification events conferring drug resistance. Despite their importance to global health, the rates at which these genomic alterations emerge have not been determined. We studied the complete genomes of different Plasmodium falciparum clones that had been propagated asexually over one year in the presence and absence of drug pressure. A combination of whole-genome microarray analysis and next-generation deep resequencing (totaling 14 terabases) revealed a stable core genome with only 38 novel single nucleotide variants appearing in seventeen evolved clones (avg. 5.4 per clone). In clones exposed to atovaquone, we found cytochrome b mutations as well as an amplification event encompassing the P. falciparum multidrug resistance associated protein (mrp1) on chromosome 1. We observed 18 large-scale (>1 kb on average) deletions of telomere-proximal regions encoding multigene families, involved in immune evasion (9.5×10⁻⁶ structural variants per base pair per generation). Six of these deletions were associated with chromosomal crossovers generated during mitosis. We found only minor differences in rates between genetically distinct strains and between parasites cultured in the presence or absence of drug. Using these derived mutation rates for P. falciparum (1.0–9.7×10⁻⁹ mutations per base pair per generation), we can now model the frequency at which drug or immune resistance alleles will emerge under a well-defined set of assumptions. Further, the detection of mitotic recombination events in var gene families illustrates how multigene families can arise and change over time in P. falciparum. These results will help improve our understanding of how P. falciparum evolves to evade control efforts within both the individual hosts and large populations.

Malaria is one of the six diseases that together are responsible for 90% of all infectious disease deaths throughout the world. The five species of Plasmodium that cause human malaria take over 655,000 lives each year. Parasites evade the immune response through antigenic variation and develop resistance to anti-malarial drugs through genetic changes in either the drug target or genes conferring resistance. We used whole-genome sequencing and microarray techniques to study evolution in P. falciparum parasites propagated in vitro for up to 180 generations. We determined the mutation rate and found that the core genome of a single clone is stable, while the subtelomeric regions are prone to acquire structural variants. These changes occur mainly in multigene families involved in immune evasion. Our findings indicate that the parasite specifically increases the sequence variability in multigene families through mitotic recombination. This high plasticity of the parasite genome suggests that multiple haplotypes will be present in a natural infection initiated by a single parasite.

Impact of chromosome ends on the biology and virulence of Plasmodium falciparum

In recent years, many studies have focused on heterochromatin located at chromosome ends, which plays an important role in regulating gene expression in many organisms ranging from yeast to humans. Similarly, in the protozoan Plasmodium falciparum, which is the most virulent human malaria parasite, the heterochromatin present in telomeres and subtelomeric regions exerts a silencing effect on the virulence gene families located therein. Studies addressing P. falciparum chromosome ends have demonstrated that these regions participate in other functions, such as the formation of the T-loop structure, the replication of telomeric regions, the regulation of telomere length and the formation of telomeric heterochromatin. In addition, telomeres are involved in anchoring chromosome ends to the nuclear periphery, thereby playing an important role in nuclear architecture and gene expression regulation. Here, we review the current understanding of chromosome ends, the proteins that bind to these regions and their impact on the biology and virulence of P. falciparum.

Two long non-coding RNAs generated from subtelomeric regions accumulate in a novel perinuclear compartment in Plasmodium falciparum

Chromosome ends have been implicated in the default silencing of clonally variant gene families in the human malaria parasite Plasmodium falciparum. These chromosome regions are organized into heterochromatin, as defined by the presence of a repressive histone H3 lysine 9 trimethylated marker and heterochromatin protein 1. Here, we show that the non-coding subtelomeric region adjacent to virulence genes forms facultative heterochromatin in a cell cycle-dependent manner. We demonstrate that telomere-associated repeat elements (TAREs) and telomeres are transcribed as long non-coding RNAs (lncRNAs) during schizogony. Northern blot assays revealed two classes of lncRNAs: a ~4-kb transcript composed of telomere sequences and a TARE-3 element, and a >6-kb transcript composed of 21-bp repeats from TARE-6. These lncRNAs are transcribed by RNA polymerase II as single-stranded molecules. RNA-FISH analysis showed that these lncRNAs form several nuclear foci during the schizont stage, whereas in the ring stage, they are located in a single perinuclear compartment that does not co-localize with any known nuclear subcompartment. Furthermore, the TARE-6 lncRNA is predicted to form a stable and repetitive hairpin structure that is able to bind histones. Consequently, the characterization of the molecular interactions of these lncRNAs with nuclear proteins may reveal novel modes of gene regulation and nuclear function in P. falciparum.

Epigenetic dysregulation of virulence gene expression in severe Plasmodium falciparum malaria

Chronic infections with the human malaria parasite Plasmodium falciparum depend on antigenic variation. P. falciparum erythrocyte membrane protein 1 (PfEMP1), the major erythrocyte surface antigen mediating parasite sequestration in the microvasculature, is encoded in parasites by a highly diverse family of var genes. Antigenic switching is mediated by clonal variation in var expression, and recent in vitro studies have demonstrated a role for epigenetic processes in var regulation. Expression of particular PfEMP1 variants may result in parasite enrichment in different tissues, a factor in the development of severe disease. Here, we study in vivo human infections and provide evidence that infection-induced stress responses in the host can modify PfEMP1 expression via the perturbation of epigenetic mechanisms. Our work suggests that severe disease may not be the direct result of an adaptive virulence strategy to maximize parasite survival but that it may indicate a loss of control of the carefully regulated process of antigenic switching that maintains chronic infections.

The role of N-terminus of Plasmodium falciparum ORC1 in telomeric localization and var gene silencing

Plasmodium falciparum origin recognition complex 1 (ORC1) protein has been implicated in DNA replication and silencing var gene family. However, the mechanism and the domain structure of ORC1 related to the regulation of var gene family are unknown. Here we show that the unique N-terminus of PfORC1 (PfORC1N_1–238) is targeted to the nuclear periphery in vivo and this region binds to the telomeric DNA in vitro due to the presence of a leucine heptad repeats. Like PfORC1N_1–238, endogenous full length ORC1, was found to be associated with sub telomeric repeat regions and promoters of various var genes. Additionally, binding and propagation of ORC1 to telomeric and subtelomeric regions was severely compromised in PfSir2 deficient parasites suggesting the dependence of endogenous ORC1 on Sir2 for var gene regulation. This feature is not previously described for Plasmodium ORC1 and contrary to yeast Saccharomyces cerevisiae where ORC function as a landing pad for Sir proteins. Interestingly, the overexpression of ORC1N_1–238 compromises the binding of Sir2 at the subtelomeric loci and var gene promoters consistent with de-repression of some var genes. These results establish role of the N-terminus of PfORC1 in heterochromatin formation and regulation of var gene expression in co-ordination with Sir2 in P. falciparum.

PfAlbas constitute a new eukaryotic DNA/RNA-binding protein family in malaria parasites

In Plasmodium falciparum, perinuclear subtelomeric chromatin conveys monoallelic expression of virulence genes. However, proteins that directly bind to chromosome ends are poorly described. Here we identify a novel DNA/RNA-binding protein family that bears homology to the archaeal protein Alba (Acetylation lowers binding affinity). We isolated three of the four PfAlba paralogs as part of a molecular complex that is associated with the P. falciparum-specific TARE6 (Telomere-Associated Repetitive Elements 6) subtelomeric region and showed in electromobility shift assays (EMSAs) that the PfAlbas bind to TARE6 repeats. In early blood stages, the PfAlba proteins were enriched at the nuclear periphery and partially co-localized with PfSir2, a TARE6-associated histone deacetylase linked to the process of antigenic variation. The nuclear location changed at the onset of parasite proliferation (trophozoite-schizont), where the PfAlba proteins were also detectable in the cytoplasm in a punctate pattern. Using single-stranded RNA (ssRNA) probes in EMSAs, we found that PfAlbas bind to ssRNA, albeit with different binding preferences. We demonstrate for the first time in eukaryotes that Alba-like proteins bind to both DNA and RNA and that their intracellular location is developmentally regulated. Discovery of the PfAlbas may provide a link between the previously described subtelomeric non-coding RNA and the regulation of antigenic variation.

Clustering of dispersed ribosomal DNA and its role in gene regulation and chromosome-end associations in malaria parasites

Dynamic changes in gene positioning contribute to differential expression of virulence-related gene families in protozoan pathogens; however, the role of nuclear architecture in gene expression in the human malaria parasite Plasmodium falciparum remains poorly understood. Here we investigated the developmentally regulated ribosomal RNA (rRNA) gene family in P. falciparum, which, unlike that in most eukaryotes, contains only a few unlinked copies of rRNA genes scattered over the subtelomeric regions of several chromosomes. We show that active and silent members of this gene family cluster in a single perinuclear nucleolus. This rDNA nuclear confinement is DNA sequence dependent, as plasmids carrying rDNA fragments are targeted to the nucleolus. Likewise, insertion of an rDNA sequence into a subtelomere from a chromosome lacking rRNA genes leads to repositioning in the nucleolus. Furthermore, we observed that rDNA spatial organization restricted interchromosomal interactions, as chromosome end-bearing rRNA genes were found to be preferentially juxtaposed, demonstrating nonrandom association of telomeres. Using Br-UTP incorporation, we observed two alpha-amanitin-resistant nucleolar transcription sites that disappeared when the rDNA cluster broke up in the replicative blood stages. Taken together, our results provide conceptual insights into functionally differentiated nuclear territories and their role in gene expression in malaria parasites.

Global sequence variation in the histidine-rich proteins 2 and 3 of Plasmodium falciparum: implications for the performance of malaria rapid diagnostic tests

Background

Accurate diagnosis is essential for prompt and appropriate treatment of malaria. While rapid diagnostic tests (RDTs) offer great potential to improve malaria diagnosis, the sensitivity of RDTs has been reported to be highly variable. One possible factor contributing to variable test performance is the diversity of parasite antigens. This is of particular concern for Plasmodium falciparum histidine-rich protein 2 (PfHRP2)-detecting RDTs since PfHRP2 has been reported to be highly variable in isolates of the Asia-Pacific region.

Methods

The pfhrp2 exon 2 fragment from 458 isolates of P. falciparum collected from 38 countries was amplified and sequenced. For a subset of 80 isolates, the exon 2 fragment of histidine-rich protein 3 (pfhrp3) was also amplified and sequenced. DNA sequence and statistical analysis of the variation observed in these genes was conducted. The potential impact of the pfhrp2 variation on RDT detection rates was examined by analysing the relationship between sequence characteristics of this gene and the results of the WHO product testing of malaria RDTs: Round 1 (2008), for 34 PfHRP2-detecting RDTs.

Results

Sequence analysis revealed extensive variations in the number and arrangement of various repeats encoded by the genes in parasite populations world-wide. However, no statistically robust correlation between gene structure and RDT detection rate for P. falciparum parasites at 200 parasites per microlitre was identified.

Conclusions

The results suggest that despite extreme sequence variation, diversity of PfHRP2 does not appear to be a major cause of RDT sensitivity variation.

A modified fluorescence in situ hybridization protocol for Plasmodium falciparum greatly improves nuclear architecture conservation

Fluorescence in situ hybridization (FISH) has been used extensively in the study of nuclear organization and gene positioning in Plasmodium falciparum. While performing FISH with published protocols, we observed large variations in parasite nuclear morphology. We hypothesized that these inconsistencies might be due to the type of parasite preparation prior to FISH, which commonly involves air-drying, prompting us to develop a new fixation protocol. Here we show both qualitatively and quantitatively that compared to air-dried and briefly fixed parasites, longer fixation in suspension leads to improved conservation of nuclear structure and lower intra-population variation of nuclear shape as well as area after FISH development. While the fixation protocol per se does not cause detectable disruptions in nuclear morphology, it greatly influences the conservation of nuclear shape and size during the most stringent steps of FISH. The type of fixation used also influences the detection of telomeric clusters, and we show that the new fixation protocol permits improved conservation of the chromosome end cluster perinuclear distribution and higher colocalization indexes for two adjacent chromosome end probes, Rep20 and telomere. Overall, the results indicate that our alternative protocol dramatically improves conservation of the nuclear architecture compared to previously reported Plasmodium DNA-FISH protocols and highlights the necessity of carefully choosing the fixation protocol for FISH.

A major role for the Plasmodium falciparum ApiAP2 protein PfSIP2 in chromosome end biology

The heterochromatic environment and physical clustering of chromosome ends at the nuclear periphery provide a functional and structural framework for antigenic variation and evolution of subtelomeric virulence gene families in the malaria parasite Plasmodium falciparum. While recent studies assigned important roles for reversible histone modifications, silent information regulator 2 and heterochromatin protein 1 (PfHP1) in epigenetic control of variegated expression, factors involved in the recruitment and organization of subtelomeric heterochromatin remain unknown. Here, we describe the purification and characterization of PfSIP2, a member of the ApiAP2 family of putative transcription factors, as the unknown nuclear factor interacting specifically with cis-acting SPE2 motif arrays in subtelomeric domains. Interestingly, SPE2 is not bound by the full-length protein but rather by a 60kDa N-terminal domain, PfSIP2-N, which is released during schizogony. Our experimental re-definition of the SPE2/PfSIP2-N interaction highlights the strict requirement of both adjacent AP2 domains and a conserved bipartite SPE2 consensus motif for high-affinity binding. Genome-wide in silico mapping identified 777 putative binding sites, 94% of which cluster in heterochromatic domains upstream of subtelomeric var genes and in telomere-associated repeat elements. Immunofluorescence and chromatin immunoprecipitation (ChIP) assays revealed co-localization of PfSIP2-N with PfHP1 at chromosome ends. Genome-wide ChIP demonstrated the exclusive binding of PfSIP2-N to subtelomeric SPE2 landmarks in vivo but not to single chromosome-internal sites. Consistent with this specialized distribution pattern, PfSIP2-N over-expression has no effect on global gene transcription. Hence, contrary to the previously proposed role for this factor in gene activation, our results provide strong evidence for the first time for the involvement of an ApiAP2 factor in heterochromatin formation and genome integrity. These findings are highly relevant for our understanding of chromosome end biology and variegated expression in P. falciparum and other eukaryotes, and for the future analysis of the role of ApiAP2-DNA interactions in parasite biology.

Telomeric heterochromatin in Plasmodium falciparum

Until very recently, little was known about the chromatin structure of the telomeres and subtelomeric regions in Plasmodium falciparum. In yeast and Drosophila melanogaster, chromatin structure has long been known to be an important aspect in the regulation and functioning of these regions. Telomeres and subtelomeric regions are enriched in epigenetic marks that are specific to heterochromatin, such as methylation of lysine 9 of histone H3 and lysine 20 of histone H4. In P. falciparum, histone modifications and the presence of both the heterochromatin "writing" (PfSir2, PKMT) and "reading" (PfHP1) machinery at telomeric and subtelomeric regions indicate that these regions are likely to have heterochromatic structure that is epigenetically regulated. This structure may be important for telomere functions such as the silencing of the var gene family implicated in the cytoadherence and antigenic variation of these parasites.

Comparative transcriptional and genomic analysis of Plasmodium falciparum field isolates

Mechanisms for differential regulation of gene expression may underlie much of the phenotypic variation and adaptability of malaria parasites. Here we describe transcriptional variation among culture-adapted field isolates of Plasmodium falciparum, the species responsible for most malarial disease. It was found that genes coding for parasite protein export into the red cell cytosol and onto its surface, and genes coding for sexual stage proteins involved in parasite transmission are up-regulated in field isolates compared with long-term laboratory isolates. Much of this variability was associated with the loss of small or large chromosomal segments, or other forms of gene copy number variation that are prevalent in the P. falciparum genome (copy number variants, CNVs). Expression levels of genes inside these segments were correlated to that of genes outside and adjacent to the segment boundaries, and this association declined with distance from the CNV boundary. This observation could not be explained by copy number variation in these adjacent genes. This suggests a local-acting regulatory role for CNVs in transcription of neighboring genes and helps explain the chromosomal clustering that we observed here. Transcriptional co-regulation of physical clusters of adaptive genes may provide a way for the parasite to readily adapt to its highly heterogeneous and strongly selective environment.