Frequent ectopic recombination of virulence factor genes in telomeric chromosome clusters of P. falciparum

Role for gene conversion in the evolution of cell-surface antigens of the malaria parasite Plasmodium falciparum

While the malaria parasite Plasmodium falciparum has low average genome-wide diversity levels, likely due to its recent introduction from a gorilla-infecting ancestor (approximately 10,000 to 50,000 years ago), some genes display extremely high diversity levels. In particular, certain proteins expressed on the surface of human red blood cell-infecting merozoites (merozoite surface proteins (MSPs)) possess exactly 2 deeply diverged lineages that have seemingly not recombined. While of considerable interest, the evolutionary origin of this phenomenon remains unknown. In this study, we analysed the genetic diversity of 2 of the most variable MSPs, DBLMSP and DBLMSP2, which are paralogs (descended from an ancestral duplication). Despite thousands of available Illumina WGS datasets from malaria-endemic countries, diversity in these genes has been hard to characterise as reads containing highly diverged alleles completely fail to align to the reference genome. To solve this, we developed a pipeline leveraging genome graphs, enabling us to genotype them at high accuracy and completeness. Using our newly- resolved sequences, we found that both genes exhibit 2 deeply diverged lineages in a specific protein domain (DBL) and that one of the 2 lineages is shared across the genes. We identified clear evidence of nonallelic gene conversion between the 2 genes as the likely mechanism behind sharing, leading us to propose that gene conversion between diverged paralogs, and not recombination suppression, can generate this surprising genealogy; a model that is furthermore consistent with high diversity levels in these 2 genes despite the strong historical P. falciparum transmission bottleneck.

Subtelomeric plasticity contributes to gene family expansion in the human parasitic flatworm Schistosoma mansoni

BACKGROUND

The genomic region that lies between the telomere and chromosome body, termed the subtelomere, is heterochromatic, repeat-rich, and frequently undergoes rearrangement. Within this region, large-scale structural changes enable gene diversification, and, as such, large multicopy gene families are often found at the subtelomere. In some parasites, genes associated with proliferation, invasion, and survival are often found in these regions, where they benefit from the subtelomere's highly plastic, rapidly changing nature. The increasing availability of complete (or near complete) parasite genomes provides an opportunity to investigate these typically poorly defined and overlooked genomic regions and potentially reveal relevant gene families necessary for the parasite's lifestyle.

RESULTS

Using the latest chromosome-scale genome assembly and hallmark repeat richness observed at chromosome termini, we have identified and characterised the subtelomeres of Schistosoma mansoni, a metazoan parasitic flatworm that infects over 250 million people worldwide. Approximately 12% of the S. mansoni genome is classified as subtelomeric, and, in line with other organisms, we find these regions to be gene-poor but rich in transposable elements. We find that S. mansoni subtelomeres have undergone extensive interchromosomal recombination and that these sites disproportionately contribute to the 2.3% of the genome derived from segmental duplications. This recombination has led to the expansion of subtelomeric gene clusters containing 103 genes, including the immunomodulatory annexins and other gene families with unknown roles. The largest of these is a 49-copy plexin domain-containing protein cluster, exclusively expressed in the tegument-the tissue located at the host-parasite physical interface-of intramolluscan life stages.

CONCLUSIONS

We propose that subtelomeric regions act as a genomic playground for trial-and-error of gene duplication and subsequent divergence. Owing to the importance of subtelomeric genes in other parasites, gene families implicated in this subtelomeric expansion within S. mansoni warrant further characterisation for a potential role in parasitism.

Hyper-diverse antigenic variation and resilience to transmission-reducing intervention in falciparum malaria

Intervention against falciparum malaria in high transmission regions remains challenging, with relaxation of control efforts typically followed by rapid resurgence. Resilience to intervention co-occurs with incomplete immunity, whereby children eventually become protected from severe disease but not infection and a large transmission reservoir results from high asymptomatic prevalence across all ages. Incomplete immunity relates to the vast antigenic variation of the parasite, with the major surface antigen of the blood stage of infection encoded by the multigene family known as var. Recent deep sampling of var sequences from individual isolates in northern Ghana showed that parasite population structure exhibited persistent features of high-transmission regions despite the considerable decrease in prevalence during transient intervention with indoor residual spraying (IRS). We ask whether despite such apparent limited impact, the transmission system had been brought close to a transition in both prevalence and resurgence ability. With a stochastic agent-based model, we investigate the existence of such a transition to pre-elimination with intervention intensity, and of molecular indicators informative of its approach. We show that resurgence ability decreases sharply and nonlinearly across a narrow region of intervention intensities in model simulations, and identify informative molecular indicators based on var gene sequences. Their application to the survey data indicates that the transmission system in northern Ghana was brought close to transition by IRS. These results suggest that sustaining and intensifying intervention would have pushed malaria dynamics to a slow-rebound regime with an increased probability of local parasite extinction.

A paradoxical population structure of var DBLα types in Africa

The var multigene family encodes the P. falciparum erythrocyte membrane protein 1 (PfEMP1), which is important in host-parasite interaction as a virulence factor and major surface antigen of the blood stages of the parasite, responsible for maintaining chronic infection. Whilst important in the biology of P. falciparum, these genes (50 to 60 genes per parasite genome) are routinely excluded from whole genome analyses due to their hyper-diversity, achieved primarily through recombination. The PfEMP1 head structure almost always consists of a DBLα-CIDR tandem. Categorised into different groups (upsA, upsB, upsC), different head structures have been associated with different ligand-binding affinities and disease severities. We study how conserved individual DBLα types are at the country, regional, and local scales in Sub-Saharan Africa. Using publicly-available sequence datasets and a novel ups classification algorithm, cUps, we performed an in silico exploration of DBLα conservation through time and space in Africa. In all three ups groups, the population structure of DBLα types in Africa consists of variants occurring at rare, low, moderate, and high frequencies. Non-rare variants were found to be temporally stable in a local area in endemic Ghana. When inspected across different geographical scales, we report different levels of conservation; while some DBLα types were consistently found in high frequencies in multiple African countries, others were conserved only locally, signifying local preservation of specific types. Underlying this population pattern is the composition of DBLα types within each isolate DBLα repertoire, revealed to also consist of a mix of types found at rare, low, moderate, and high frequencies in the population. We further discuss the adaptive forces and balancing selection, including host genetic factors, potentially shaping the evolution and diversity of DBLα types in Africa.

Fungal antigenic variation using mosaicism and reassortment of subtelomeric genes' repertoires

Surface antigenic variation is crucial for major pathogens that infect humans. To escape the immune system, they exploit various mechanisms. Understanding these mechanisms is important to better prevent and fight the deadly diseases caused. Those used by the fungus Pneumocystis jirovecii that causes life-threatening pneumonia in immunocompromised individuals remain poorly understood. Here, though this fungus is currently not cultivable, our detailed analysis of the subtelomeric sequence motifs and genes encoding surface proteins suggests that the system involves the reassortment of the repertoire of ca. 80 non-expressed genes present in each strain, from which single genes are retrieved for mutually exclusive expression. Dispersion of the new repertoires, supposedly by healthy carrier individuals, appears very efficient because identical alleles are observed in patients from different countries. Our observations reveal a unique strategy of antigenic variation. They also highlight the possible role in genome rearrangements of small imperfect mirror sequences forming DNA triplexes.

Epigenetic Regulation and Chromatin Remodeling in Malaria Parasites

Plasmodium falciparum, the human malaria parasite, infects two hosts and various cell types, inducing distinct morphological and physiological changes in the parasite in response to different environmental conditions. These variations required the parasite to adapt and develop elaborate molecular mechanisms to ensure its spread and transmission. Recent findings have significantly improved our understanding of the regulation of gene expression in P. falciparum. Here, we provide an up-to-date overview of technologies used to highlight the transcriptomic adjustments occurring in the parasite throughout its life cycle. We also emphasize the complementary and complex epigenetic mechanisms regulating gene expression in malaria parasites. This review concludes with an outlook on the chromatin architecture, the remodeling systems, and how this 3D genome organization is critical in various biological processes.

Uncovering the role of Rad51 in homologous recombination-mediated antigenic diversification in the human malaria parasite Plasmodium falciparum

The human malaria parasite Plasmodium falciparum maintains the chronicity of infections through antigenic variation, a well-coordinated immune evasion mechanism. The most prominent molecular determinant of antigenic variation in this parasite includes the members of the var multigene family. Homologous recombination (HR)-mediated genomic rearrangements have been implicated to play a major role in var gene diversification. However, the key molecular factors involved in the generation of diversity at var loci are less known. Here, we tested the hypothesis that PfRad51 could carry out recombination between var genes that are not homologous but homeologous in nature. We employed the whole-genome sequencing (WGS) approach to investigate recombination events among var sequences over 100 generations and compared the rate of sequence rearrangement at the var loci in both PfRad51-proficient and -deficient parasite lines. This brief report provides evidence that the loss of the key recombinase function renders the parasite with inefficient HR and results in fewer recombination events among the var sequences, thereby impacting the diversification of the var gene repertoire.

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.

Plasmodium falciparum gametocytes display global chromatin remodelling during sexual differentiation

BACKGROUND

The protozoan malaria parasite Plasmodium falciparum has a complex life cycle during which it needs to differentiate into multiple morphologically distinct life forms. A key process for transmission of the disease is the development of male and female gametocytes in the human blood, yet the mechanisms determining sexual dimorphism in these haploid, genetically identical sexual precursor cells remain largely unknown. To understand the epigenetic program underlying the differentiation of male and female gametocytes, we separated the two sexual forms by flow cytometry and performed RNAseq as well as comprehensive ChIPseq profiling of several histone variants and modifications.

RESULTS

We show that in female gametocytes the chromatin landscape is globally remodelled with respect to genome-wide patterns and combinatorial usage of histone variants and histone modifications. We identified sex specific differences in heterochromatin distribution, implicating exported proteins and ncRNAs in sex determination. Specifically in female gametocytes, the histone variants H2A.Z/H2B.Z were highly enriched in H3K9me3-associated heterochromatin. H3K27ac occupancy correlated with stage-specific gene expression, but in contrast to asexual parasites this was unlinked to H3K4me3 co-occupancy at promoters in female gametocytes.

CONCLUSIONS

Collectively, we defined novel combinatorial chromatin states differentially organising the genome in gametocytes and asexual parasites and unravelled fundamental, sex-specific differences in the epigenetic code. Our chromatin maps represent an important resource for future understanding of the mechanisms driving sexual differentiation in P. falciparum.

Analysis of genome instability and implications for the consequent phenotype in Plasmodium falciparum containing mutated MSH2-1 (P513T)

Malarial parasites exhibit extensive genomic plasticity, which induces the antigen diversification and the development of antimalarial drug resistance. Only a few studies have examined the genome maintenance mechanisms of parasites. The study aimed at elucidating the impact of a mutation in a DNA mismatch repair gene on genome stability by maintaining the mutant and wild-type parasites through serial in vitro cultures for approximately 400 days and analysing the subsequent spontaneous mutations. A P513T mutant of the DNA mismatch repair protein PfMSH2-1 from Plasmodium falciparum 3D7 was created. The mutation did not influence the base substitution rate but significantly increased the insertion/deletion (indel) mutation rate in short tandem repeats (STRs) and minisatellite loci. STR mutability was affected by allele size, genomic category and certain repeat motifs. In the mutants, significant telomere healing and homologous recombination at chromosomal ends caused extensive gene loss and generation of chimeric genes, resulting in large-scale chromosomal alteration. Additionally, the mutant showed increased tolerance to N-methyl-N'-nitro-N-nitrosoguanidine, suggesting that PfMSH2-1 was involved in recognizing DNA methylation damage. This work provides valuable insights into the role of PfMSH2-1 in genome stability and demonstrates that the genomic destabilization caused by its dysfunction may lead to antigen diversification.

PfSET2 Is Involved in Genome Organization of Var Gene Family in Plasmodium falciparum

The three-dimensional (3D) genome structure of human malaria parasite Plasmodium falciparum is highly organized and plays important roles in regulating coordinated expression patterns of specific genes such as virulence genes which are involved in antigenic variation and immune escape. However, the molecular mechanisms that control 3D genome of the parasite remain elusive. Here, by analyzing genome organization of P. falciparum, we identify high-interacting regions (HIRs) with strong chromatin interactions at telomeres and virulence genes loci. Specifically, HIRs are highly enriched with repressive histone marks (H3K36me3 and H3K9me3) and form the transcriptional repressive center. Deletion of PfSET2, which controls H3K36me3 level, results in marked reduction of both intrachromosomal and interchromosomal interactions for HIRs. Importantly, such chromatin reorganization coordinates with dynamic changes in epigenetic feature in HIRs and transcriptional activation of var genes. Additionally, different cluster of var genes based on the pattern of chromatin interactions show distinct transcriptional activation potential after deletion of PfSET2. Our results uncover a fundamental mechanism that the epigenetic factor PfSET2 controls the 3D organization of heterochromatin to regulate the transcription activities of var genes family in P. falciparum. IMPORTANCE PfSET2 has been reported to play key role in silencing var genes in Plasmodium falciparum, while the underlying molecular mechanisms remain unclear. Here, we provide evidence that PfSET2 is essential to maintain 3D genome organization of heterochromatin region to keep var genes in transcription repressive state. These findings can contribute better understanding of the regulation of high-order chromatin structure in P. falciparum.

Life in plastic, it's fantastic! How Leishmania exploit genome instability to shape gene expression

Leishmania are kinetoplastid pathogens that cause leishmaniasis, a debilitating and potentially life-threatening infection if untreated. Unusually, Leishmania regulate their gene expression largely post-transcriptionally due to the arrangement of their coding genes into polycistronic transcription units that may contain 100s of functionally unrelated genes. Yet, Leishmania are capable of rapid and responsive changes in gene expression to challenging environments, often instead correlating with dynamic changes in their genome composition, ranging from chromosome and gene copy number variations to the generation of extrachromosomal DNA and the accumulation of point mutations. Typically, such events indicate genome instability in other eukaryotes, coinciding with genetic abnormalities, but for Leishmania, exploiting these products of genome instability can provide selectable substrates to catalyse necessary gene expression changes by modifying gene copy number. Unorthodox DNA replication, DNA repair, replication stress factors and DNA repeats are recognised in Leishmania as contributors to this intrinsic instability, but how Leishmania regulate genome plasticity to enhance fitness whilst limiting toxic under- or over-expression of co-amplified and co-transcribed genes is unclear. Herein, we focus on fresh, and detailed insights that improve our understanding of genome plasticity in Leishmania. Furthermore, we discuss emerging models and factors that potentially circumvent regulatory issues arising from polycistronic transcription. Lastly, we highlight key gaps in our understanding of Leishmania genome plasticity and discuss future studies to define, in higher resolution, these complex regulatory interactions.

Unravelling var complexity: Relationship between DBLα types and var genes in Plasmodium falciparum

The enormous diversity and complexity of var genes that diversify rapidly by recombination has led to the exclusion of assembly of these genes from major genome initiatives (e.g., Pf6). A scalable solution in epidemiological surveillance of var genes is to use a small ‘tag’ region encoding the immunogenic DBLα domain as a marker to estimate var diversity. As var genes diversify by recombination, it is not clear the extent to which the same tag can appear in multiple var genes. This relationship between marker and gene has not been investigated in natural populations. Analyses of in vitro recombination within and between var genes have suggested that this relationship would not be exclusive. Using a dataset of publicly-available assembled var sequences, we test this hypothesis by studying DBLα-var relationships for four study sites in four countries: Pursat (Cambodia) and Mae Sot (Thailand), representing low malaria transmission, and Navrongo (Ghana) and Chikwawa (Malawi), representing high malaria transmission. In all study sites, DBLα-var relationships were shown to be predominantly 1-to-1, followed by a second largest proportion of 1-to-2 DBLα-var relationships. This finding indicates that DBLα tags can be used to estimate not just DBLα diversity but var gene diversity when applied in a local endemic area. Epidemiological applications of this result are discussed.

Genome-wide signatures of adaptation to extreme environments in red algae

The high temperature, acidity, and heavy metal-rich environments associated with hot springs have a major impact on biological processes in resident cells. One group of photosynthetic eukaryotes, the Cyanidiophyceae (Rhodophyta), has successfully thrived in hot springs and associated sites worldwide for more than 1 billion years. Here, we analyze chromosome-level assemblies from three representative Cyanidiophyceae species to study environmental adaptation at the genomic level. We find that subtelomeric gene duplication of functional genes and loss of canonical eukaryotic traits played a major role in environmental adaptation, in addition to horizontal gene transfer events. Shared responses to environmental stress exist in Cyanidiales and Galdieriales, however, most of the adaptive genes (e.g., for arsenic detoxification) evolved independently in these lineages. Our results underline the power of local selection to shape eukaryotic genomes that may face vastly different stresses in adjacent, extreme microhabitats.

Repurposing of Plasmodium falciparum var genes beyond the blood stage

A commonly observed survival strategy in protozoan parasites is the sequential expression of clonally variant-surface antigens to avoid elimination by the host's immune response. In malaria-causing P. falciparum, the immunovariant erythrocyte-membrane protein-1 (PfEMP1) adhesin family, encoded by var genes, is responsible for both antigenic variation and cytoadherence of infected erythrocytes to the microvasculature. Until recently, the biological function of these variant genes was believed to be restricted to intraerythrocytic developmental stages. With the advent of new technologies, var gene expression has been confirmed in transmission and pre-erythrocytic stages. Here, we discuss how repurposing of var gene expression beyond chronic blood-stage infection may be critical for successful transmission.

Chromatin structure can introduce systematic biases in genome-wide analyses of Plasmodium falciparum

Background: The maintenance, regulation, and dynamics of heterochromatin in the human malaria parasite, Plasmodium falciparum, has drawn increasing attention due to its regulatory role in mutually exclusive virulence gene expression and the silencing of key developmental regulators. The advent of genome-wide analyses such as chromatin-immunoprecipitation followed by sequencing (ChIP-seq) has been instrumental in understanding chromatin composition; however, even in model organisms, ChIP-seq experiments are susceptible to intrinsic experimental biases arising from underlying chromatin structure. Methods: We performed a control ChIP-seq experiment, re-analyzed previously published ChIP-seq datasets and compared different analysis approaches to characterize biases of genome-wide analyses in P. falciparum. Results: We found that heterochromatic regions in input control samples used for ChIP-seq normalization are systematically underrepresented in regard to sequencing coverage across the P. falciparum genome. This underrepresentation, in combination with a non-specific or inefficient immunoprecipitation, can lead to the identification of false enrichment and peaks across these regions. We observed that such biases can also be seen at background levels in specific and efficient ChIP-seq experiments. We further report on how different read mapping approaches can also skew sequencing coverage within highly similar subtelomeric regions and virulence gene families. To ameliorate these issues, we discuss orthogonal methods that can be used to characterize bona fide chromatin-associated proteins. Conclusions: Our results highlight the impact of chromatin structure on genome-wide analyses in the parasite and the need for caution when characterizing chromatin-associated proteins and features.

Differential Trafficking and Expression of PIR Proteins in Acute and Chronic Plasmodium Infections

Plasmodium multigene families are thought to play important roles in the pathogenesis of malaria. Plasmodium interspersed repeat (pir) genes comprise the largest multigene family in many Plasmodium species. However, their expression pattern and localisation remain to be elucidated. Understanding protein subcellular localisation is fundamental to reveal the functional importance and cell-cell interactions of the PIR proteins. Here, we use the rodent malaria parasite, Plasmodium chabaudi chabaudi, as a model to investigate the localisation pattern of this gene family. We found that most PIR proteins are co-expressed in clusters during acute and chronic infection; members of the S7 clade are predominantly expressed during the acute-phase, whereas members of the L1 clade dominate the chronic-phase of infection. Using peptide antisera specific for S7 or L1 PIRS, we show that these PIRs have different localisations within the infected red blood cells. S7 PIRs are exported into the infected red blood cell cytoplasm where they are co-localised with parasite-induced host cell modifications termed Maurer’s clefts, whereas L1 PIRs are localised on or close to the parasitophorous vacuolar membrane. This localisation pattern changes following mosquito transmission and during progression from acute- to chronic-phase of infection. The presence of PIRs in Maurer’s clefts, as seen for Plasmodium falciparum RIFIN and STEVOR proteins, might suggest trafficking of the PIRs on the surface of the infected erythrocytes. However, neither S7 nor L1 PIR proteins detected by the peptide antisera are localised on the surface of infected red blood cells, suggesting that they are unlikely to be targets of surface variant-specific antibodies or to be directly involved in adhesion of infected red blood cells to host cells, as described for Plasmodium falciparum VAR proteins. The differences in subcellular localisation of the two major clades of Plasmodium chabaudi PIRs across the blood cycle, and the apparent lack of expression on the red cell surface strongly suggest that the function(s) of this gene family may differ from those of other multigene families of Plasmodium, such as the var genes of Plasmodium falciparum.

DNA replication dynamics during erythrocytic schizogony in the malaria parasites Plasmodium falciparum and Plasmodium knowlesi

Malaria parasites are unusual, early-diverging protozoans with non-canonical cell cycles. They do not undergo binary fission, but divide primarily by schizogony. This involves the asynchronous production of multiple nuclei within the same cytoplasm, culminating in a single mass cytokinesis event. The rate and efficiency of parasite reproduction is fundamentally important to malarial disease, which tends to be severe in hosts with high parasite loads. Here, we have studied for the first time the dynamics of schizogony in two human malaria parasite species, Plasmodium falciparum and Plasmodium knowlesi. These differ in their cell-cycle length, the number of progeny produced and the genome composition, among other factors. Comparing them could therefore yield new information about the parameters and limitations of schizogony. We report that the dynamics of schizogony differ significantly between these two species, most strikingly in the gap phases between successive nuclear multiplications, which are longer in P. falciparum and shorter, but more heterogenous, in P. knowlesi. In both species, gaps become longer as schizogony progresses, whereas each period of active DNA replication grows shorter. In both species there is also extreme variability between individual cells, with some schizonts producing many more nuclei than others, and some individual nuclei arresting their DNA replication for many hours while adjacent nuclei continue to replicate. The efficiency of schizogony is probably influenced by a complex set of factors in both the parasite and its host cell.

Malaria parasites are unusual, early-diverging single-celled organisms. One of their atypical features is their mode of producing new cells. Most cells replicate their genome, segregate the copies into two nuclei and then split the whole cell in two: a process called binary fission. Malaria parasites, by contrast, multiply primarily by schizogony. Each cell replicates its genome several times over, asynchronously, generating many nuclei within the same cytoplasm, before splitting into many new daughter cells in a single mass event. All malaria species do this, but there are stark differences between species in how long schizogony takes and how many progeny each cell can produce. Understanding this is important because the rate of parasite reproduction is fundamentally important to malarial disease: humans who have a high burden of parasites in their blood are most likely to suffer severe malaria. Here we compare the process of schizogony in two different species, by developing cell lines in which we can follow de novo DNA replication at high spatial and temporal resolution, at both whole-cell and single-molecule levels. We establish the dynamics of schizogony, highlighting similarities and differences between two species and setting parameters for future studies of interventions that could interfere with parasite reproduction.

Integration of DNA Repair, Antigenic Variation, Cytoadhesion, and Chance in Babesia Survival: A Perspective

Apicomplexan parasites live in hostile environments in which they are challenged chemically and their hosts attempt in many ways to kill them. In response, the parasites have evolved multiple mechanisms that take advantage of these challenges to enhance their survival. Perhaps the most impressive example is the evolutionary co-option of DNA repair mechanisms by the parasites as a means to rapidly manipulate the structure, antigenicity, and expression of the products of specific multigene families. The purpose of variant proteins that mediate cytoadhesion has long been thought to be primarily the avoidance of splenic clearance. Based upon known biology, I present an alternative perspective in which it is survival of the oxidative environment within which Babesia spp. parasites live that has driven integration of DNA repair, antigenic variation, and cytoadhesion, and speculate on how genome organization affects that integration. This perspective has ramifications for the development of parasite control strategies.

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.

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.

An accurate method for identifying recent recombinants from unaligned sequences

Motivation

Recombination is a fundamental process in molecular evolution, and the identification of recombinant sequences is thus of major interest. However, current methods for detecting recombinants are primarily designed for aligned sequences. Thus they struggle with analyses of highly diverse genes, such as the var genes of the malaria parasite Plasmodium falciparum, which are known to diversify primarily through recombination.

Results

We introduce an algorithm to detect recent recombinant sequences from a dataset without a full multiple alignment. Our algorithm can handle thousands of gene-length sequences without the need for a reference panel. We demonstrate the accuracy of our algorithm through extensive numerical simulations; in particular, it maintains its effectiveness in the presence of insertions and deletions. We apply our algorithm to a dataset of 17,335 DBLα types in var genes from Ghana, observing that sequences belonging to the same ups group or domain subclass recombine amongst themselves more frequently, and that non-recombinant DBLα types are more conserved than recombinant ones.

Availability

Source code is freely available at https://github.com/qianfeng2/detREC_program.

Supplementary information

Supplementary data are available at Bioinformatics online.

OUP accepted manuscript

Heterochromatin-associated gene silencing controls multiple physiological processes in malaria parasites, however, little is known concerning the regulatory network and cis-acting sequences involved in the organization of heterochromatin and how they modulate heterochromatic gene expression. Based on systematic profiling of genome-wide occupancy of eighteen Apicomplexan AP2 transcription factors by ChIP-seq analysis, we identify and characterize eight heterochromatin-associated factors (PfAP2-HFs), which exhibit preferential enrichment within heterochromatic regions but with differential coverage profiles. Although these ApiAP2s target euchromatic gene loci via specific DNA motifs, they are likely integral components of heterochromatin independent of DNA motif recognition. Systematic knockout screenings of ApiAP2 factors coupled with RNA-seq transcriptomic profiling revealed three activators and three repressors of heterochromatic gene expression including four PfAP2-HFs. Notably, expression of virulence genes is either completely silenced or significantly reduced upon the depletion of PfAP2-HC. Integrated multi-omics analyses reveal autoregulation and feed-forward loops to be common features of the ApiAP2 regulatory network, in addition to the occurrence of dynamic interplay between local chromatin structure and ApiAP2s in transcriptional control. Collectively, this study provides a valuable resource describing the genome-wide landscape of the ApiAP2 family and insights into functional divergence and cooperation within this family during the blood-stage development of malaria parasites.

Resetting var Gene Transcription in Plasmodium falciparum

One of the key mechanisms contributing to the virulence of Plasmodium falciparum is its ability to undergo antigenic switching among antigenically distinct variants of the PfEMP1 adhesive proteins, encoded by the var gene family. To avoid premature exposure of its antigenic repertoire, the parasite transcribes its var genes in a mutually exclusive manner, and switch expression at a very slow rate. This process is epigenetically regulated and it relies on "epigenetic memory," which imprints the single active var gene to remain active for multiple replication cycles. Erasing this epigenetic memory in parasites grown in culture resembles parasites, which egress from the liver. It could therefore be of interest for investigating var switching patterns at the onset of malaria infections. In addition, this procedure could be used for creating heterogeneity of var expression among parasite populations. The methodology described here for resetting of var gene expression is based on promoter titration, also known as molecular sponging.

Analysis of pir gene expression across the Plasmodium life cycle

Background

Plasmodium interspersed repeat (pir) is the largest multigene family in the genomes of most Plasmodium species. A variety of functions for the PIR proteins which they encode have been proposed, including antigenic variation, immune evasion, sequestration and rosetting. However, direct evidence for these is lacking. The repetitive nature of the family has made it difficult to determine function experimentally. However, there has been some success in using gene expression studies to suggest roles for some members in virulence and chronic infection.

Methods

Here pir gene expression was examined across the life cycle of Plasmodium berghei using publicly available RNAseq data-sets, and at high resolution in the intraerythrocytic development cycle using new data from Plasmodium chabaudi.

Results

Expression of pir genes is greatest in stages of the parasite which invade and reside in red blood cells. The marked exception is that liver merozoites and male gametocytes produce a very large number of pir gene transcripts, notably compared to female gametocytes, which produce relatively few. Within the asexual blood stages different subfamilies peak at different times, suggesting further functional distinctions. Representing a subfamily of its own, the highly conserved ancestral pir gene warrants further investigation due to its potential tractability for functional investigation. It is highly transcribed in multiple life cycle stages and across most studied Plasmodium species and thus is likely to play an important role in parasite biology.

Conclusions

The identification of distinct expression patterns for different pir genes and subfamilies is likely to provide a basis for the design of future experiments to uncover their function.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-021-03979-6.

Subtelomeric Chromatin in the Fission Yeast S. pombe

Telomeres play important roles in safeguarding the genome. The specialized repressive chromatin that assembles at telomeres and subtelomeric domains is key to this protective role. However, in many organisms, the repetitive nature of telomeric and subtelomeric sequences has hindered research efforts. The fission yeast S. pombe has provided an important model system for dissection of chromatin biology due to the relative ease of genetic manipulation and strong conservation of important regulatory proteins with higher eukaryotes. Telomeres and the telomere-binding shelterin complex are highly conserved with mammals, as is the assembly of constitutive heterochromatin at subtelomeres. In this review, we seek to summarize recent work detailing the assembly of distinct chromatin structures within subtelomeric domains in fission yeast. These include the heterochromatic SH subtelomeric domains, the telomere-associated sequences (TAS), and ST chromatin domains that assemble highly condensed chromatin clusters called knobs. Specifically, we review new insights into the sequence of subtelomeric domains, the distinct types of chromatin that assemble on these sequences and how histone H3 K36 modifications influence these chromatin structures. We address the interplay between the subdomains of chromatin structure and how subtelomeric chromatin is influenced by both the telomere-bound shelterin complexes and by euchromatic chromatin regulators internal to the subtelomeric domain. Finally, we demonstrate that telomere clustering, which is mediated via the condensed ST chromatin knob domains, does not depend on knob assembly within these domains but on Set2, which mediates H3K36 methylation.

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.

A unique GCN5 histone acetyltransferase complex controls erythrocyte invasion and virulence in the malaria parasite Plasmodium falciparum

The histone acetyltransferase GCN5-associated SAGA complex is evolutionarily conserved from yeast to human and functions as a general transcription co-activator in global gene regulation. In this study, we identified a divergent GCN5 complex in Plasmodium falciparum, which contains two plant homeodomain (PHD) proteins (PfPHD1 and PfPHD2) and a plant apetela2 (AP2)-domain transcription factor (PfAP2-LT). To dissect the functions of the PfGCN5 complex, we generated parasite lines with either the bromodomain in PfGCN5 or the PHD domain in PfPHD1 deleted. The two deletion mutants closely phenocopied each other, exhibiting significantly reduced merozoite invasion of erythrocytes and elevated sexual conversion. These domain deletions caused dramatic decreases not only in histone H3K9 acetylation but also in H3K4 trimethylation, indicating synergistic crosstalk between the two euchromatin marks. Domain deletion in either PfGCN5 or PfPHD1 profoundly disturbed the global transcription pattern, causing altered expression of more than 60% of the genes. At the schizont stage, these domain deletions were linked to specific down-regulation of merozoite genes involved in erythrocyte invasion, many of which contain the AP2-LT binding motif and are also regulated by AP2-I and BDP1, suggesting targeted recruitment of the PfGCN5 complex to the invasion genes by these specific factors. Conversely, at the ring stage, PfGCN5 or PfPHD1 domain deletions disrupted the mutually exclusive expression pattern of the entire var gene family, which encodes the virulent factor PfEMP1. Correlation analysis between the chromatin state and alteration of gene expression demonstrated that up- and down-regulated genes in these mutants are highly correlated with the silent and active chromatin states in the wild-type parasite, respectively. Collectively, the PfGCN5 complex represents a novel HAT complex with a unique subunit composition including an AP2 transcription factor, which signifies a new paradigm for targeting the co-activator complex to regulate general and parasite-specific cellular processes in this low-branching parasitic protist.

Epigenetic regulation of gene expression plays essential roles in orchestrating the general and parasite-specific cellular pathways in the malaria parasite Plasmodium falciparum. To better understand the epigenetic mechanisms in this parasite, we characterized the histone acetyltransferase GCN5-mediated transcription regulation during intraerythrocytic development of the parasite. Using tandem affinity purification and proteomic characterization, we identified that the PfGCN5-associated complex contains nine core components, including two PHD domain proteins (PfPHD1 and PfPHD2) and an AP2-domain transcription factor, which is divergent from the canonical GCN5 complexes evolutionarily conserved from yeast to human. To understand the functions of the PfGCN5 complex, we performed domain deletions in two subunits of this complex, PfGCN5 and PfPHD1. We found that the two deletion mutants displayed very similar growth phenotypes, including significantly reduced merozoite invasion rates and elevated sexual conversion. These two mutants were associated with dramatic decreases in histone H3K9 acetylation and H3K4 trimethylation, which led to global changes in chromatin states and gene expression. Consistent with the phenotypes, genes significantly affected by the PfGCN5 and PfPHD1 gene disruption include those participating in parasite-specific pathways such as invasion, virulence, and sexual development. In conclusion, this study presents a new model of the PfGCN5 complex for targeting the co-activator complex to regulate general and parasite-specific cellular processes in this low-branching parasitic protist.

Cell wall protein variation, break-induced replication, and subtelomere dynamics in Candida glabrata

Candida glabrata is an opportunistic pathogen of humans, responsible for up to 30% of disseminated candidiasis. Adherence of C. glabrata to host cells is mediated by adhesin-like proteins (ALPs), about half of which are encoded in the subtelomeres. We performed a de novo assembly of two C. glabrata strains, BG2 and BG3993, using long single-molecule real-time (SMRT) reads, and constructed high-quality telomere-to-telomere assemblies of all 13 chromosomes to assess differences between C. glabrata strains. We documented variation between strains, and in agreement with earlier studies, found high (~0.5%-1%) frequencies of SNVs across the genome, including within subtelomeric regions. We documented changes in ALP gene structure and complement: there are large length differences in ALP genes in different strains, resulting from copy number variation in tandem repeats. We compared strains to characterize chromosome rearrangement events including within the poorly characterized subtelomeric regions. We show that rearrangements within the subtelomere regions all affect ALP-encoding genes, and 14/16 involve just the most terminal ALP gene. We present evidence that these rearrangements are mediated by break-induced replication. This study highlights the constrained nature of subtelomeric changes impacting ALP gene complement and subtelomere structure.

Evolutionary analyses of the major variant surface antigen-encoding genes reveal population structure of Plasmodium falciparum within and between continents

Malaria remains a major public health problem in many countries. Unlike influenza and HIV, where diversity in immunodominant surface antigens is understood geographically to inform disease surveillance, relatively little is known about the global population structure of PfEMP1, the major variant surface antigen of the malaria parasite Plasmodium falciparum. The complexity of the var multigene family that encodes PfEMP1 and that diversifies by recombination, has so far precluded its use in malaria surveillance. Recent studies have demonstrated that cost-effective deep sequencing of the region of var genes encoding the PfEMP1 DBLα domain and subsequent classification of within host sequences at 96% identity to define unique DBLα types, can reveal structure and strain dynamics within countries. However, to date there has not been a comprehensive comparison of these DBLα types between countries. By leveraging a bioinformatic approach (jumping hidden Markov model) designed specifically for the analysis of recombination within var genes and applying it to a dataset of DBLα types from 10 countries, we are able to describe population structure of DBLα types at the global scale. The sensitivity of the approach allows for the comparison of the global dataset to ape samples of Plasmodium Laverania species. Our analyses show that the evolution of the parasite population emerging out of Africa underlies current patterns of DBLα type diversity. Most importantly, we can distinguish geographic population structure within Africa between Gabon and Ghana in West Africa and Uganda in East Africa. Our evolutionary findings have translational implications in the context of globalization. Firstly, DBLα type diversity can provide a simple diagnostic framework for geographic surveillance of the rapidly evolving transmission dynamics of P. falciparum. It can also inform efforts to understand the presence or absence of global, regional and local population immunity to major surface antigen variants. Additionally, we identify a number of highly conserved DBLα types that are present globally that may be of biological significance and warrant further characterization.

The Role of the Histone Methyltransferase PfSET10 in Antigenic Variation by Malaria Parasites: a Cautionary Tale

The identification of specific epigenetic regulatory proteins in infectious organisms has become a high-profile research topic and a focus for several drug development initiatives. However, studies that define specific roles for different epigenetic modifiers occasionally report differing results, and we similarly provide evidence regarding the histone methyltransferase PfSET10 that is in stark contrast with previously published results.

The virulence of the malaria parasite Plasmodium falciparum is due in large part to its ability to avoid immune destruction through antigenic variation. This results from changes in expression within the multicopy var gene family that encodes the surface antigen P. falciparum erythrocyte protein one (PfEMP1). Understanding the mechanisms underlying this process has been a high-profile research focus for many years. The histone methyltransferase PfSET10 was previously identified as a key enzyme required both for parasite viability and for regulating var gene expression, thus making it a prominent target for developing antimalarial intervention strategies and the subject of considerable research focus. Here, however, we show that disruption of the gene encoding PfSET10 is not lethal and has no effect on var gene expression, in sharp contrast with previously published reports. The contradictory findings highlight the importance of reevaluating previous conclusions when new technologies become available and suggest the possibility of a previously unappreciated plasticity in epigenetic gene regulation in P. falciparum.

IMPORTANCE The identification of specific epigenetic regulatory proteins in infectious organisms has become a high-profile research topic and a focus for several drug development initiatives. However, studies that define specific roles for different epigenetic modifiers occasionally report differing results, and we similarly provide evidence regarding the histone methyltransferase PfSET10 that is in stark contrast with previously published results. We believe that the conflicting results, rather than suggesting erroneous conclusions, instead reflect the importance of revisiting previous conclusions using newly developed methodologies, as well as caution in interpreting seemingly contrary results in fields that are known to display considerable plasticity, for example metabolism and epigenetics.

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.

Trypanosoma Cruzi Genome: Organization, Multi-Gene Families, Transcription, and Biological Implications

Chagas disease caused by the parasite Trypanosoma cruzi affects millions of people. Although its first genome dates from 2005, its complexity hindered a complete assembly and annotation. However, the new sequencing methods have improved genome annotation of some strains elucidating the broad genetic diversity and complexity of this parasite. Here, we reviewed the genomic structure and regulation, the genetic diversity, and the analysis of the principal multi-gene families of the recent genomes for several strains. The telomeric and sub-telomeric regions are sites with high recombination events, the genome displays two different compartments, the core and the disruptive, and the genome plasticity seems to play a key role in the survival and the infection process. Trypanosoma cruzi (T. cruzi) genome is composed mainly of multi-gene families as the trans-sialidases, mucins, and mucin-associated surface proteins. Trans-sialidases are the most abundant genes in the genome and show an important role in the effectiveness of the infection and the parasite survival. Mucins and MASPs are also important glycosylated proteins of the surface of the parasite that play a major biological role in both insect and mammal-dwelling stages. Altogether, these studies confirm the complexity of T. cruzi genome revealing relevant concepts to better understand Chagas disease.

The relationship between protein domains and homopeptides in the Plasmodium falciparum proteome

The proteome of the malaria parasite Plasmodium falciparum is notable for the pervasive occurrence of homopeptides or low-complexity regions (i.e., regions that are made from a small subset of amino-acid residue types). The most prevalent of these are made from residues encoded by adenine/thymidine (AT)-rich codons, in particular asparagine. We examined homopeptide occurrences within protein domains in P. falciparum. Homopeptide enrichments occur for hydrophobic (e.g., valine), or small residues (alanine or glycine) in short spans (<5 residues), but these enrichments disappear for longer lengths. We observe that short asparagine homopeptides (<10 residues long) have a dramatic relative depletion inside protein domains, indicating some selective constraint to keep them from forming. We surmise that this is possibly linked to co-translational protein folding, although there are specific protein domains that are enriched in longer asparagine homopeptides (≥10 residues) indicating a functional linkage for specific poly-asparagine tracts. Top gene ontology functional category enrichments for homopeptides associated with diverse protein domains include “vesicle-mediated transport”, and “DNA-directed 5′-3′ RNA polymerase activity”, with various categories linked to “binding” evidencing significant homopeptide depletions. Also, in general homopeptides are substantially enriched in the parts of protein domains that are near/in IDRs. The implications of these findings are discussed.

Transposon-mediated telomere destabilization: a driver of genome evolution in the blast fungus

The fungus Magnaporthe oryzae causes devastating diseases of crops, including rice and wheat, and in various grasses. Strains from ryegrasses have highly unstable chromosome ends that undergo frequent rearrangements, and this has been associated with the presence of retrotransposons (Magnaporthe oryzae Telomeric Retrotransposons-MoTeRs) inserted in the telomeres. The objective of the present study was to determine the mechanisms by which MoTeRs promote telomere instability. Targeted cloning, mapping, and sequencing of parental and novel telomeric restriction fragments (TRFs), along with MinION sequencing of genomic DNA allowed us to document the precise molecular alterations underlying 109 newly-formed TRFs. These included truncations of subterminal rDNA sequences; acquisition of MoTeR insertions by 'plain' telomeres; insertion of the MAGGY retrotransposons into MoTeR arrays; MoTeR-independent expansion and contraction of subtelomeric tandem repeats; and a variety of rearrangements initiated through breaks in interstitial telomere tracts that are generated during MoTeR integration. Overall, we estimate that alterations occurred in approximately sixty percent of chromosomes (one in three telomeres) analyzed. Most importantly, we describe an entirely new mechanism by which transposons can promote genomic alterations at exceptionally high frequencies, and in a manner that can promote genome evolution while minimizing collateral damage to overall chromosome architecture and function.

Humanized Mice and the Rebirth of Malaria Genetic Crosses

The first experimental crosses carried out with the human malaria parasite Plasmodium falciparum played a key role in determining the genetic loci responsible for drug resistance, virulence, invasion, growth rate, and transmission. These crosses relied on splenectomized chimpanzees to complete the liver stage of the parasite's life cycle and the subsequent transition to asexual blood stage culture followed by cloning of recombinant progeny in vitro. Crosses can now be routinely carried out using human-liver-chimeric mice infused with human erythrocytes to generate hundreds of unique recombinant progeny for genetic linkage mapping, bulk segregant analysis, and high-throughput 'omics readouts. The high number of recombinant progeny should allow for unprecedented power and efficiency in the execution of a systems genetics approach to study P. falciparum biology.

Detection of simple and complex de novo mutations with multiple reference sequences

The characterization of de novo mutations in regions of high sequence and structural diversity from whole-genome sequencing data remains highly challenging. Complex structural variants tend to arise in regions of high repetitiveness and low complexity, challenging both de novo assembly, in which short reads do not capture the long-range context required for resolution, and mapping approaches, in which improper alignment of reads to a reference genome that is highly diverged from that of the sample can lead to false or partial calls. Long-read technologies can potentially solve such problems but are currently unfeasible to use at scale. Here we present Corticall, a graph-based method that combines the advantages of multiple technologies and prior data sources to detect arbitrary classes of genetic variant. We construct multisample, colored de Bruijn graphs from short-read data for all samples, align long-read–derived haplotypes and multiple reference data sources to restore graph connectivity information, and call variants using graph path-finding algorithms and a model for simultaneous alignment and recombination. We validate and evaluate the approach using extensive simulations and use it to characterize the rate and spectrum of de novo mutation events in 119 progeny from four Plasmodium falciparum experimental crosses, using long-read data on the parents to inform reconstructions of the progeny and to detect several known and novel nonallelic homologous recombination events.

Rapid activation of distinct members of multigene families in Plasmodium spp

The genomes of Plasmodium spp. encode a number of different multigene families that are thought to play a critical role for survival. However, with the exception of the P. falciparum var genes, very little is known about the biological roles of any of the other multigene families. Using the recently developed Selection Linked Integration method, we have been able to activate the expression of a single member of a multigene family of our choice in Plasmodium spp. from its endogenous promoter. We demonstrate the usefulness of this approach by activating the expression of a unique var, rifin and stevor in P. falciparum as well as yir in P. yoelii. Characterization of the selected parasites reveals differences between the different families in terms of mutual exclusive control, co-regulation, and host adaptation. Our results further support the application of the approach for the study of multigene families in Plasmodium and other organisms.

Omelianczyk, Loh et al. activate the expression of a single member of a multigene family in Plasmodium spp. from its endogenous promoter, identifying differences between the different families. This study supports the application of the Selection Linked Integration method for studying multigene families in Plasmodium.

A Targeted "Next-Generation" Sequencing-Informatic Approach to Define Genetic Diversity in Theileria orientalis Populations within Individual Cattle: Proof-of-Principle

Oriental theileriosis is an economically important tickborne disease of bovines, caused by some members of the Theileria orientalis complex. Currently, 11 distinct operational taxonomic units (OTUs), or genotypes, are recognized based on their major piroplasm surface protein (MPSP) gene sequences. Two of these genotypes (i.e., chitose and ikeda) are recognized as pathogenic in cattle, causing significant disease in countries of the Asia-Pacific region. However, the true extent of genetic variation and associated virulence/pathogenicity within this complex is unknown. Here, we undertook a proof-of-principle study of a small panel of genomic DNAs (n = 13) from blood samples originating from individual cattle known to harbor T. orientalis, in order to assess the performance of a targeted “next-generation” sequencing-informatic approach to identify genotypes. Five genotypes (chitose, ikeda, buffeli, type 4, and type 5) were defined; multiple genotypes were found within individual samples, with dominant and minor sequence types representing most genotypes. This study indicates that this sequencing-informatic workflow could be useful to assess the nature and extent of genetic variation within and among populations of T. orientalis on a large scale, and to potentially employ panels of distinct gene markers for expanded molecular epidemiological investigations of socioeconomically important protistan pathogens more generally.

The role of epigenetics and chromatin structure in transcriptional regulation in malaria parasites

Due to the unique selective pressures and extreme changes faced by the human malaria parasite Plasmodium falciparum throughout its life cycle, the parasite has evolved distinct features to alter its gene expression patterns. Along with classical gene regulation by transcription factors (TFs), of which only one family, the AP2 TFs, has been described in the parasite genome, a large body of evidence points toward chromatin structure and epigenetic factors mediating the changes in gene expression associated with parasite life cycle stages. These attributes may be critically important for immune evasion, host cell invasion and development of the parasite in its two hosts, the human and the Anopheles vector. Thus, the factors involved in the maintenance and regulation of chromatin and epigenetic features represent potential targets for antimalarial drugs. In this review, we discuss the mechanisms in P. falciparum that regulate chromatin structure, nucleosome landscape, the 3-dimensional structure of the genome and additional distinctive features created by parasite-specific genes and gene families. We review conserved traits of chromatin in eukaryotes in order to highlight what is unique in the parasite.

Transcriptional variation in malaria parasites: why and how

Transcriptional differences enable the generation of alternative phenotypes from the same genome. In malaria parasites, transcriptional plasticity plays a major role in the process of adaptation to fluctuations in the environment. Multiple studies with culture-adapted parasites and field isolates are starting to unravel the different transcriptional alternatives available to Plasmodium falciparum and the underlying molecular mechanisms. Here we discuss how epigenetic variation, directed transcriptional responses and also genetic changes that affect transcript levels can all contribute to transcriptional variation and, ultimately, parasite survival. Some transcriptional changes are driven by stochastic events. These changes can occur spontaneously, resulting in heterogeneity within parasite populations that provides the grounds for adaptation by dynamic natural selection. However, transcriptional changes can also occur in response to external cues. A better understanding of the mechanisms that the parasite has evolved to alter its transcriptome may ultimately contribute to the design of strategies to combat malaria to which the parasite cannot adapt.

Conserved associations between G-quadruplex-forming DNA motifs and virulence gene families in malaria parasites

Background

The Plasmodium genus of malaria parasites encodes several families of antigen-encoding genes. These genes tend to be hyper-variable, highly recombinogenic and variantly expressed. The best-characterized family is the var genes, exclusively found in the Laveranian subgenus of malaria parasites infecting humans and great apes. Var genes encode major virulence factors involved in immune evasion and the maintenance of chronic infections. In the human parasite P. falciparum, var gene recombination and diversification appear to be promoted by G-quadruplex (G4) DNA motifs, which are strongly associated with var genes in P. falciparum. Here, we investigated how this association might have evolved across Plasmodium species – both Laverania and also more distantly related species which lack vars but encode other, more ancient variant gene families.

Results

The association between var genes and G4-forming motifs was conserved across Laverania, spanning ~ 1 million years of evolutionary time, with suggestive evidence for evolution of the association occurring within this subgenus. In rodent malaria species, G4-forming motifs were somewhat associated with pir genes, but this was not conserved in the Laverania, nor did we find a strong association of these motifs with any gene family in a second outgroup of avian malaria parasites. Secondly, we compared two different G4 prediction algorithms in their performance on extremely A/T-rich Plasmodium genomes, and also compared these predictions with experimental data from G4-seq, a DNA sequencing method for identifying G4-forming motifs. We found a surprising lack of concordance between the two algorithms and also between the algorithms and G4-seq data.

Conclusions

G4-forming motifs are uniquely strongly associated with Plasmodium var genes, suggesting a particular role for G4s in recombination and diversification of these genes. Secondly, in the A/T-rich genomes of Plasmodium species, the choice of prediction algorithm may be particularly influential when studying G4s in these important protozoan pathogens.

Diversity and Complexity of the Large Surface Protein Family in the Compacted Genomes of Multiple Pneumocystis Species

Pneumocystis, a major opportunistic pathogen in patients with a broad range of immunodeficiencies, contains abundant surface proteins encoded by a multicopy gene family, termed the major surface glycoprotein (Msg) gene superfamily. This superfamily has been identified in all Pneumocystis species characterized to date, highlighting its important role in Pneumocystis biology. In this report, through a comprehensive and in-depth characterization of 459 msg genes from 7 Pneumocystis species, we demonstrate, for the first time, the phylogeny and evolution of conserved domains in Msg proteins and provide a detailed description of the classification, unique characteristics, and phylogenetic relatedness of five Msg families. We further describe, for the first time, the relative expression levels of individual msg families in two rodent Pneumocystis species, the substantial variability of the msg repertoires in P. carinii from laboratory and wild rats, and the distinct features of the expression site for the classic msg genes in Pneumocystis from 8 mammalian host species. Our analysis suggests multiple functions for this superfamily rather than just conferring antigenic variation to allow immune evasion as previously believed. This study provides a rich source of information that lays the foundation for the continued experimental exploration of the functions of the Msg superfamily in Pneumocystis biology.IMPORTANCEPneumocystis continues to be a major cause of disease in humans with immunodeficiency, especially those with HIV/AIDS and organ transplants, and is being seen with increasing frequency worldwide in patients treated with immunodepleting monoclonal antibodies. Annual health care associated with Pneumocystis pneumonia costs ∼$475 million dollars in the United States alone. In addition to causing overt disease in immunodeficient individuals, Pneumocystis can cause subclinical infection or colonization in healthy individuals, which may play an important role in species preservation and disease transmission. Our work sheds new light on the diversity and complexity of the msg superfamily and strongly suggests that the versatility of this superfamily reflects multiple functions, including antigenic variation to allow immune evasion and optimal adaptation to host environmental conditions to promote efficient infection and transmission. These findings are essential to consider in developing new diagnostic and therapeutic strategies.

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.

PhyloToL: A Taxon/Gene-Rich Phylogenomic Pipeline to Explore Genome Evolution of Diverse Eukaryotes

Estimating multiple sequence alignments (MSAs) and inferring phylogenies are essential for many aspects of comparative biology. Yet, many bioinformatics tools for such analyses have focused on specific clades, with greatest attention paid to plants, animals, and fungi. The rapid increase in high-throughput sequencing (HTS) data from diverse lineages now provides opportunities to estimate evolutionary relationships and gene family evolution across the eukaryotic tree of life. At the same time, these types of data are known to be error-prone (e.g., substitutions, contamination). To address these opportunities and challenges, we have refined a phylogenomic pipeline, now named PhyloToL, to allow easy incorporation of data from HTS studies, to automate production of both MSAs and gene trees, and to identify and remove contaminants. PhyloToL is designed for phylogenomic analyses of diverse lineages across the tree of life (i.e., at scales of >100 My). We demonstrate the power of PhyloToL by assessing stop codon usage in Ciliophora, identifying contamination in a taxon- and gene-rich database and exploring the evolutionary history of chromosomes in the kinetoplastid parasite Trypanosoma brucei, the causative agent of African sleeping sickness. Benchmarking PhyloToL’s homology assessment against that of OrthoMCL and a published paper on superfamilies of bacterial and eukaryotic organellar outer membrane pore-forming proteins demonstrates the power of our approach for determining gene family membership and inferring gene trees. PhyloToL is highly flexible and allows users to easily explore HTS data, test hypotheses about phylogeny and gene family evolution and combine outputs with third-party tools (e.g., PhyloChromoMap, iGTP).

Evolutionary analysis of the most polymorphic gene family in falciparum malaria

The var gene family of the human malaria parasite Plasmodium falciparum encode proteins that are crucial determinants of both pathogenesis and immune evasion and are highly polymorphic. Here we have assembled nearly complete var gene repertoires from 2398 field isolates and analysed a normalised set of 714 from across 12 countries. This therefore represents the first large scale attempt to catalogue the worldwide distribution of var gene sequences We confirm the extreme polymorphism of this gene family but also demonstrate an unexpected level of sequence sharing both within and between continents. We show that this is likely due to both the remnants of selective sweeps as well as a worrying degree of recent gene flow across continents with implications for the spread of drug resistance. We also address the evolution of the var repertoire with respect to the ancestral genes within the Laverania and show that diversity generated by recombination is concentrated in a number of hotspots. An analysis of the subdomain structure indicates that some existing definitions may need to be revised From the analysis of this data, we can now understand the way in which the family has evolved and how the diversity is continuously being generated. Finally, we demonstrate that because the genes are distributed across the genome, sequence sharing between genotypes acts as a useful population genetic marker.

Three-Dimensional Genome Organization and Virulence in Apicomplexan Parasites

Mounting evidence supports the idea that epigenetic, and the overall 3-dimensional (3D) architecture of the genome, plays an important role in gene expression for eukaryotic organisms. We recently used Hi-C methodologies to generate and compare the 3D genome of 7 different apicomplexan parasites, including several pathogenic and less pathogenic malaria parasites as well as related human parasites Babesia microti and Toxoplasma gondii. Our goal was to understand the possible relationship between genome organization, gene expression, and pathogenicity of these infectious agents. Collectively, our results demonstrate that spatial genome organization in most Plasmodium species is constrained by the colocalization of virulence genes that are unique in their effect on chromosome folding, indicating a link between genome organization and gene expression in more virulent pathogens.

Homologous recombination shapes the genetic diversity of African swine fever viruses

The African swine fever virus (ASFV) has severely influenced the swine industry of the world. Currently, there is no effective vaccine or drugs against the ASFV. How to effectively control the virus is challenging. In this study, we have analyzed all the publicly available ASFV genomes and demonstrated that there was a large genetic diversity of ASFV genomes. Interestingly, the genetic diversity was mainly caused by extensive genomic insertions and/or deletions (indels) instead of the point mutations. Further analyses showed that the indels may be attributed much to the homologous recombination, as supported by significant associations between the occurrence of extensive recombination events and the indels in the ASFV genomes. Besides, the homologous recombination also led to changes of gene content of ASFVs. Finally, repeated elements of dozens of nucleotides in length were observed to widely distribute and cluster in the adjacent positions of ASFV genomes, which may facilitate the occurrence of homologous recombination. This work highlighted the importance of homologous recombination in shaping the genetic diversity of the ASFVs, and could help understand the evolution of the virus.

Plasmodium Genomics and Genetics: New Insights into Malaria Pathogenesis, Drug Resistance, Epidemiology, and Evolution

Protozoan Plasmodium parasites are the causative agents of malaria, a deadly disease that continues to afflict hundreds of millions of people every year. Infections with malaria parasites can be asymptomatic, with mild or severe symptoms, or fatal, depending on many factors such as parasite virulence and host immune status. Malaria can be treated with various drugs, with artemisinin-based combination therapies (ACTs) being the first-line choice. Recent advances in genetics and genomics of malaria parasites have contributed greatly to our understanding of parasite population dynamics, transmission, drug responses, and pathogenesis. However, knowledge gaps in parasite biology and host-parasite interactions still remain. Parasites resistant to multiple antimalarial drugs have emerged, while advanced clinical trials have shown partial efficacy for one available vaccine. Here we discuss genetic and genomic studies of Plasmodium biology, host-parasite interactions, population structures, mosquito infectivity, antigenic variation, and targets for treatment and immunization. Knowledge from these studies will advance our understanding of malaria pathogenesis, epidemiology, and evolution and will support work to discover and develop new medicines and vaccines.