Genome wide distribution of G-quadruplexes and their impact on gene expression in malaria parasites

Deciphering the Plasmodium falciparum perinuclear var gene expression site

The protozoan parasite Plasmodium falciparum, responsible for the deadliest form of human malaria, employs antigenic variation via monoallelic expression as a key survival strategy. The selective activation of one out of the 60-member var gene family is key to understanding the parasite's ability to cause severe disease and evade the host immune response. var gene activation is initiated by its relocation to a specialized expression site. While the perinuclear expression site (PES) plays a crucial role in enabling the expression of a single allele, the characteristics of this PES remain largely obscure. Recent breakthroughs in genome editing tools and the discovery of regulatory noncoding RNAs have shed light on this intriguing biological feature, offering significant insights into the mechanisms of pathogen virulence.

DNA minicircles capable of forming a variety of non-canonical structural motifs

Although more than 10% of the human genome has the potential to fold into non-B DNA, the formation of non-canonical structural motifs as part of long dsDNA chains are usually considered as unfavorable from a thermodynamic point of view. However, recent experiments have confirmed that non-canonical motifs do exist and are non-randomly distributed in genomic DNA. This distribution is highly dependent not only on the DNA sequence but also on various other factors such as environmental conditions, DNA topology and the expression of specific cellular factors in different cell types. In this study, we describe a new strategy used in the preparation of DNA minicircles containing different non-canonical motifs which arise as a result of imperfect base pairing between complementary strands. The approach exploits the fact that imperfections in the pairing of complementary strands thermodynamically weaken the dsDNA structure at the expense of enhancing the formation of non-canonical motifs. In this study, a completely different concept of stable integration of a non-canonical motif into dsDNA is presented. Our approach allows the integration of various types of non-canonical motifs into the dsDNA structure such as hairpin, cruciform, G-quadruplex and i-motif forms but also combinations of these forms. Small DNA minicircles have recently become the subject of considerable interest in both fundamental research and in terms of their potential therapeutic applications.

The genetic landscape of origins of replication in P. falciparum

Various origin mapping approaches have enabled genome-wide identification of origins of replication (ORI) in model organisms, but only a few studies have focused on divergent organisms. By employing three complementary approaches we provide a high-resolution map of ORIs in Plasmodium falciparum, the deadliest human malaria parasite. We profiled the distribution of origin of recognition complex (ORC) binding sites by ChIP-seq of two PfORC subunits and mapped active ORIs using NFS and SNS-seq. We show that ORIs lack sequence specificity but are not randomly distributed, and group in clusters. Licensing is biased towards regions of higher GC content and associated with G-quadruplex forming sequences (G4FS). While strong transcription likely enhances firing, active origins are depleted from transcription start sites. Instead, most accumulate in transcriptionally active gene bodies. Single molecule analysis of nanopore reads containing multiple initiation events, which could have only come from individual nuclei, showed a relationship between the replication fork pace and the distance to the nearest origin. While some similarities were drawn with the canonic eukaryote model, the distribution of ORIs in P. falciparum is likely shaped by unique genomic features such as extreme AT-richness-a product of evolutionary pressure imposed by the parasitic lifestyle.

Effects of the G-quadruplex-binding drugs quarfloxin and CX-5461 on the malaria parasite Plasmodium falciparum

Plasmodium falciparum is the deadliest causative agent of human malaria. This parasite has historically developed resistance to most drugs, including the current frontline treatments, so new therapeutic targets are needed. Our previous work on guanine quadruplexes (G4s) in the parasite's DNA and RNA has highlighted their influence on parasite biology, and revealed G4 stabilising compounds as promising candidates for repositioning. In particular, quarfloxin, a former anticancer agent, kills blood-stage parasites at all developmental stages, with fast rates of kill and nanomolar potency. Here we explored the molecular mechanism of quarfloxin and its related derivative CX-5461. In vitro, both compounds bound to P. falciparum-encoded G4 sequences. In cellulo, quarfloxin was more potent than CX-5461, and could prevent establishment of blood-stage malaria in vivo in a murine model. CX-5461 showed clear DNA damaging activity, as reported in human cells, while quarfloxin caused weaker signatures of DNA damage. Both compounds caused transcriptional dysregulation in the parasite, but the affected genes were largely different, again suggesting different modes of action. Therefore, CX-5461 may act primarily as a DNA damaging agent in both Plasmodium parasites and mammalian cells, whereas the complete antimalarial mode of action of quarfloxin may be parasite-specific and remains somewhat elusive.

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.

The role of long noncoding RNAs in malaria parasites

The human malaria parasites, including Plasmodium falciparum, persist as a major cause of global morbidity and mortality. The recent stalling of progress toward malaria elimination substantiates a need for novel interventions. Controlled gene expression is central to the parasite's numerous life cycle transformations and adaptation. With few specific transcription factors (TFs) identified, crucial roles for chromatin states and epigenetics in parasite transcription have become evident. Although many chromatin-modifying enzymes are known, less is known about which factors mediate their impacts on transcriptional variation. Like those of higher eukaryotes, long noncoding RNAs (lncRNAs) have recently been shown to have integral roles in parasite gene regulation. This review aims to summarize recent developments and key findings on the role of lncRNAs in P. falciparum.

Noncanonical DNA structures are drivers of genome evolution

In addition to the canonical right-handed double helix, other DNA structures, termed 'non-B DNA', can form in the genomes across the tree of life. Non-B DNA regulates multiple cellular processes, including replication and transcription, yet its presence is associated with elevated mutagenicity and genome instability. These discordant cellular roles fuel the enormous potential of non-B DNA to drive genomic and phenotypic evolution. Here we discuss recent studies establishing non-B DNA structures as novel functional elements subject to natural selection, affecting evolution of transposable elements (TEs), and specifying centromeres. By highlighting the contributions of non-B DNA to repeated evolution and adaptation to changing environments, we conclude that evolutionary analyses should include a perspective of not only DNA sequence, but also its structure.

Discovery of RUF6 ncRNA-interacting proteins involved in P. falciparum immune evasion

Non-coding RNAs (ncRNAs) are emerging regulators of immune evasion and transmission of Plasmodium falciparum RUF6 is an ncRNA gene family that is transcribed by RNA polymerase III but actively regulates the Pol II-transcribed var virulence gene family. Understanding how RUF6 ncRNA connects to downstream effectors is lacking. We developed an RNA-directed proteomic discovery (ChIRP-MS) protocol to identify in vivo RUF6 ncRNA-protein interactions. The RUF6 ncRNA interactome was purified with biotinylated antisense oligonucleotides. Quantitative label-free mass spectrometry identified several unique proteins linked to gene transcription including RNA Pol II subunits, nucleosome assembly proteins, and a homologue of DEAD box helicase 5 (DDX5). Affinity purification of Pf-DDX5 identified proteins originally found by our RUF6-ChIRP protocol, validating the technique's robustness for identifying ncRNA interactomes in P. falciparum Inducible displacement of nuclear Pf-DDX5 resulted in significant down-regulation of the active var gene. Our work identifies a RUF6 ncRNA-protein complex that interacts with RNA Pol II to sustain the var gene expression, including a helicase that may resolve G-quadruplex secondary structures in var genes to facilitate transcriptional activation and progression.

Evolutionary Invariant of the Structure of DNA Double Helix in RNAP II Core Promoters

Eukaryotic and archaeal RNA polymerase II (POL II) machinery is highly conserved, regardless of the extreme changes in promoter sequences in different organisms. The goal of our work is to find the cause of this conservatism. The representative sets of aligned promoter sequences of fifteen organisms belonging to different evolutional stages were studied. Their textual profiles, as well as profiles of the indexes that characterize the secondary structure and the mechanical and physicochemical properties, were analyzed. The evolutionarily stable, extremely heterogeneous special secondary structure of POL II core promoters was revealed, which includes two singular regions—hexanucleotide “INR” around TSS and octanucleotide “TATA element” of about −28 bp upstream. Such structures may have developed at some stage of evolution. It turned out to be so well matched for the pre-initiation complex formation and the subsequent initiation of transcription for POL II machinery that in the course of evolution there were selected only those nucleotide sequences that were able to reproduce these structural properties. The individual features of specific sequences representing the singular region of the promoter of each gene can affect the kinetics of DNA-protein complex formation and facilitate strand separation in double-stranded DNA at the TSS position.

Iso-FRET: an isothermal competition assay to analyze quadruplex formation in vitro

Algorithms have been widely used to predict G-quadruplexes (G4s)-prone sequences. However, an experimental validation of these predictions is generally required. We previously reported a high-throughput technique to evidence G4 formation in vitro called FRET-MC. This method, while convenient and reproducible, has one known weakness: its inability to pin point G4 motifs of low thermal stability. As such quadruplexes may still be biologically relevant if formed at physiological temperature, we wanted to develop an independent assay to overcome this limitation. To this aim, we introduced an isothermal version of the competition assay, called iso-FRET, based on a duplex-quadruplex competition and a well-characterized bis-quinolinium G4 ligand, PhenDC3. G4-forming competitors act as decoys for PhenDC3, lowering its ability to stabilize the G4-forming motif reporter oligonucleotide conjugated to a fluorescence quencher (37Q). The decrease in available G4 ligand concentration restores the ability of 37Q to hybridize to its FAM-labeled short complementary C-rich strand (F22), leading to a decrease in fluorescence signal. In contrast, when no G4-forming competitor is present, PhenDC3 remains available to stabilize the 37Q quadruplex, preventing the formation of the F22 + 37Q complex. Iso-FRET was first applied to a reference panel of 70 sequences, and then used to investigate 23 different viral sequences.

G-quadruplex occurrence and conservation: more than just a question of guanine–cytosine content

G-quadruplexes are motifs found in DNA and RNA that can fold into tertiary structures. Until now, they have been studied experimentally mainly in humans and a few other species. Recently, predictions have been made with bacterial and archaeal genomes. Nevertheless, a global comparison of predicted G4s (pG4s) across and within the three living kingdoms has not been addressed. In this study, we aimed to predict G4s in genes and transcripts of all kingdoms of living organisms and investigated the differences in their distributions. The relation of the predictions with GC content was studied. It appears that GC content is not the only parameter impacting G4 predictions and abundance. The distribution of pG4 densities varies depending on the class of transcripts and the group of species. Indeed, we have observed that, in coding transcripts, there are more predicted G4s than expected for eukaryotes but not for archaea and bacteria, while in noncoding transcripts, there are as many or fewer predicted G4s in all species groups. We even noticed that some species with the same GC content presented different pG4 profiles. For instance, Leishmania major and Chlamydomonas reinhardtii both have 60% of GC content, but the former has a pG4 density of 0.07 and the latter 1.16.

Major Achievements in the Design of Quadruplex-Interactive Small Molecules

Organic small molecules that can recognize and bind to G-quadruplex and i-Motif nucleic acids have great potential as selective drugs or as tools in drug target discovery programs, or even in the development of nanodevices for medical diagnosis. Hundreds of quadruplex-interactive small molecules have been reported, and the challenges in their design vary with the intended application. Herein, we survey the major achievements on the therapeutic potential of such quadruplex ligands, their mode of binding, effects upon interaction with quadruplexes, and consider the opportunities and challenges for their exploitation in drug discovery.

Identification of an Aptamer With Binding Specificity to Tumor-Homing Myeloid-Derived Suppressor Cells

Myeloid-derived suppressor cells (MDSCs) play a critical role in tumor growth and metastasis. Since they constantly infiltrate into the tumor tissue, these cells are considered as an ideal carrier for tumor-targeted drug delivery. We recently identified a DNA-based thioaptamer (T1) with tumor accumulating activity, demonstrated its potential on tumor targeting and drug delivery. In the current study, we have carried out structure-activity relationship analysis to further optimize the aptamer. In the process, we have identified a sequence-modified aptamer (M1) that shows an enhanced binding affinity to MDSCs over the parental T1 aptamer. In addition, M1 can penetrate into the tumor tissue more effectively by hitchhiking on MDSCs. Taken together, we have identified a new reagent for enhanced tumor-targeted drug delivery.

Structured Waters Mediate Small Molecule Binding to G-Quadruplex Nucleic Acids

The role of G-quadruplexes in human cancers is increasingly well-defined. Accordingly, G-quadruplexes can be suitable drug targets and many small molecules have been identified to date as G-quadruplex binders, some using computer-based design methods and co-crystal structures. The role of bound water molecules in the crystal structures of G-quadruplex-small molecule complexes has been analyzed in this study, focusing on the water arrangements in several G-quadruplex ligand complexes. One is the complex between the tetrasubstituted naphthalene diimide compound MM41 and a human intramolecular telomeric DNA G-quadruplex, and the others are in substituted acridine bimolecular G-quadruplex complexes. Bridging water molecules form most of the hydrogen-bond contacts between ligands and DNA in the parallel G-quadruplex structures examined here. Clusters of structured water molecules play essential roles in mediating between ligand side chain groups/chromophore core and G-quadruplex. These clusters tend to be conserved between complex and native G-quadruplex structures, suggesting that they more generally serve as platforms for ligand binding, and should be taken into account in docking and in silico studies.

G-quadruplex RNA motifs influence gene expression in the malaria parasite Plasmodium falciparum

G-quadruplexes are non-helical secondary structures that can fold in vivo in both DNA and RNA. In human cells, they can influence replication, transcription and telomere maintenance in DNA, or translation, transcript processing and stability of RNA. We have previously showed that G-quadruplexes are detectable in the DNA of the malaria parasite Plasmodium falciparum, despite a very highly A/T-biased genome with unusually few guanine-rich sequences. Here, we show that RNA G-quadruplexes can also form in P. falciparum RNA, using rG4-seq for transcriptome-wide structure-specific RNA probing. Many of the motifs, detected here via the rG4seeker pipeline, have non-canonical forms and would not be predicted by standard in silico algorithms. However, in vitro biophysical assays verified formation of non-canonical motifs. The G-quadruplexes in the P. falciparum transcriptome are frequently clustered in certain genes and associated with regions encoding low-complexity peptide repeats. They are overrepresented in particular classes of genes, notably those that encode PfEMP1 virulence factors, stress response genes and DNA binding proteins. In vitro translation experiments and in vivo measures of translation efficiency showed that G-quadruplexes can influence the translation of P. falciparum mRNAs. Thus, the G-quadruplex is a novel player in post-transcriptional regulation of gene expression in this major human pathogen.

Guanine quadruplexes and their roles in molecular processes

The role of guanine quadruplexes (G4) in fundamental biological processes like DNA replication, transcription, translation and telomere maintenance is recognized. G4 structure dynamics is regulated by G4 structure binding proteins and is thought to be crucial for the maintenance of genome integrity in both prokaryotic and eukaryotic cells. Growing research over the last decade has expanded the existing knowledge of the functional diversity of G4 (DNA and RNA) structures across the working models. The control of G4 structure dynamics using G4 binding drugs has been suggested as the putative targets in the control of cancer and bacterial pathogenesis. This review has brought forth the collections of recent information that indicate G4 (mostly G4 DNA) roles in microbial pathogenesis, DNA damaging stress response in bacteria and mammalian cells. Studies in mitochondrial gene function regulation by G4s have also been underscored. Finally, the interdependence of G4s and epigenetic modifications and their speculated medical implications through G4 interacting proteins has been discussed.

Unlocking G-Quadruplexes as Antiviral Targets

Guanine-rich DNA and RNA sequences can fold into noncanonical nucleic acid structures called G-quadruplexes (G4s). Since the discovery that these structures may act as scaffolds for the binding of specific ligands, G4s aroused the attention of a growing number of scientists. The versatile roles of G4 structures in viral replication, transcription, and translation suggest direct applications in therapy or diagnostics. G4-interacting molecules (proteins or small molecules) may also affect the balance between latent and lytic phases, and increasing evidence reveals that G4s are implicated in generally suppressing viral processes, such as replication, transcription, translation, or reverse transcription. In this review, we focus on the discovery of G4s in viruses and the role of G4 ligands in the antiviral drug discovery process. After assessing the role of viral G4s, we argue that host G4s participate in immune modulation, viral tumorigenesis, cellular pathways involved in virus maturation, and DNA integration of viral genomes, which can be potentially employed for antiviral therapeutics. Furthermore, we scrutinize the impediments and shortcomings in the process of studying G4 ligands and drug discovery. Finally, some unanswered questions regarding viral G4s are highlighted for prospective future projects. SIGNIFICANCE STATEMENT: G-quadruplexes (G4s) are noncanonical nucleic acid structures that have gained increasing recognition during the last few decades. First identified as relevant targets in oncology, their importance in virology is now increasingly clear. A number of G-quadruplex ligands are known: viral transcription and replication are the main targets of these ligands. Both viral and cellular G4s may be targeted; this review embraces the different aspects of G-quadruplexes in both host and viral contexts.

Identifying and validating the presence of Guanine-Quadruplexes (G4) within the blood fluke parasite Schistosoma mansoni

Schistosomiasis is a neglected tropical disease that currently affects over 250 million individuals worldwide. In the absence of an immunoprophylactic vaccine and the recognition that mono-chemotherapeutic control of schistosomiasis by praziquantel has limitations, new strategies for managing disease burden are urgently needed. A better understanding of schistosome biology could identify previously undocumented areas suitable for the development of novel interventions. Here, for the first time, we detail the presence of G-quadruplexes (G4) and putative quadruplex forming sequences (PQS) within the Schistosoma mansoni genome. We find that G4 are present in both intragenic and intergenic regions of the seven autosomes as well as the sex-defining allosome pair. Amongst intragenic regions, G4 are particularly enriched in 3´ UTR regions. Gene Ontology (GO) term analysis evidenced significant G4 enrichment in the wnt signalling pathway (p<0.05) and PQS oligonucleotides synthetically derived from wnt-related genes resolve into parallel and anti-parallel G4 motifs as elucidated by circular dichroism (CD) spectroscopy. Finally, utilising a single chain anti-G4 antibody called BG4, we confirm the in situ presence of G4 within both adult female and male worm nuclei. These results collectively suggest that G4-targeted compounds could be tested as novel anthelmintic agents and highlights the possibility that G4-stabilizing molecules could be progressed as candidates for the treatment of schistosomiasis.

Evaluating the Influence of a G-Quadruplex Prone Sequence on the Transactivation Potential by Wild-Type and/or Mutant P53 Family Proteins through a Yeast-Based Functional Assay

P53, P63, and P73 proteins belong to the P53 family of transcription factors, sharing a common gene organization that, from the P1 and P2 promoters, produces two groups of mRNAs encoding proteins with different N-terminal regions; moreover, alternative splicing events at C-terminus further contribute to the generation of multiple isoforms. P53 family proteins can influence a plethora of cellular pathways mainly through the direct binding to specific DNA sequences known as response elements (REs), and the transactivation of the corresponding target genes. However, the transcriptional activation by P53 family members can be regulated at multiple levels, including the DNA topology at responsive promoters. Here, by using a yeast-based functional assay, we evaluated the influence that a G-quadruplex (G4) prone sequence adjacent to the p53 RE derived from the apoptotic PUMA target gene can exert on the transactivation potential of full-length and N-terminal truncated P53 family α isoforms (wild-type and mutant). Our results show that the presence of a G4 prone sequence upstream or downstream of the P53 RE leads to significant changes in the relative activity of P53 family proteins, emphasizing the potential role of structural DNA features as modifiers of P53 family functions at target promoter sites.

G-quadruplexes in H1N1 influenza genomes

BACKGROUND

Influenza viruses are dangerous pathogens. Seventy-Seven genomes of recently emerged genotype 4 reassortant Eurasian avian-like H1N1 virus (G4-EA-H1N1) are currently available. We investigated the presence and variation of potential G-quadruplex forming sequences (PQS), which can serve as targets for antiviral treatment.

RESULTS

PQS were identified in all 77 genomes. The total number of PQS in G4-EA-H1N1 genomes was 571. Interestingly, the number of PQS per genome in individual close relative viruses varied from 4 to 12. PQS were not randomly distributed in the 8 segments of the G4-EA-H1N1 genome, the highest frequency of PQS being found in the NP segment (1.39 per 1000 nt), which is considered a potential target for antiviral therapy. In contrast, no PQS was found in the NS segment. Analyses of variability pointed the importance of some PQS; even if genome variation of influenza virus is extreme, the PQS with the highest G4Hunter score is the most conserved in all tested genomes. G-quadruplex formation in vitro was experimentally confirmed using spectroscopic methods.

CONCLUSIONS

The results presented here hint several G-quadruplex-forming sequences in G4-EA-H1N1 genomes, that could provide good therapeutic targets.

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.

G-Quadruplexes in the Archaea Domain

The importance of unusual DNA structures in the regulation of basic cellular processes is an emerging field of research. Amongst local non-B DNA structures, G-quadruplexes (G4s) have gained in popularity during the last decade, and their presence and functional relevance at the DNA and RNA level has been demonstrated in a number of viral, bacterial, and eukaryotic genomes, including humans. Here, we performed the first systematic search of G4-forming sequences in all archaeal genomes available in the NCBI database. In this article, we investigate the presence and locations of G-quadruplex forming sequences using the G4Hunter algorithm. G-quadruplex-prone sequences were identified in all archaeal species, with highly significant differences in frequency, from 0.037 to 15.31 potential quadruplex sequences per kb. While G4 forming sequences were extremely abundant in Hadesarchaea archeon (strikingly, more than 50% of the Hadesarchaea archaeon isolate WYZ-LMO6 genome is a potential part of a G4-motif), they were very rare in the Parvarchaeota phylum. The presence of G-quadruplex forming sequences does not follow a random distribution with an over-representation in non-coding RNA, suggesting possible roles for ncRNA regulation. These data illustrate the unique and non-random localization of G-quadruplexes in Archaea.