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.

G4access identifies G-quadruplexes and their associations with open chromatin and imprinting control regions.

Metazoan promoters are enriched in secondary DNA structure-forming motifs, such as G-quadruplexes (G4s). Here we describe 'G4access', an approach to isolate and sequence G4s associated with open chromatin via nuclease digestion. G4access is antibody- and crosslinking-independent and enriches for computationally predicted G4s (pG4s), most of which are confirmed in vitro. Using G4access in human and mouse cells, we identify cell-type-specific G4 enrichment correlated with nucleosome exclusion and promoter transcription. G4access allows measurement of variations in G4 repertoire usage following G4 ligand treatment, HDAC and G4 helicases inhibitors. Applying G4access to cells from reciprocal hybrid mouse crosses suggests a role for G4s in the control of active imprinting regions. Consistently, we also observed that G4access peaks are unmethylated, while methylation at pG4s correlates with nucleosome repositioning on DNA. Overall, our study provides a new tool for studying G4s in cellular dynamics and highlights their association with open chromatin, transcription and their antagonism to DNA methylation.

Guidelines for G-quadruplexes: I. In vitro characterization.

Besides the well-known DNA double-helix, non-canonical nucleic acid structures regulate crucial biological activities. Among these oddities, guanine-rich DNA sequences can form unusual four-stranded secondary structures called G-quadruplexes (G4s). G4-prone sequences have been found in the genomes of most species, and G4s play important roles in essential processes such as transcription, replication, genome integrity and epigenetic regulation. Here, we present a short overview of G-quadruplexes followed by a detailed description of the biophysical and biochemical methods used to characterize G4s in vitro. The principles, experimental details and possible shortcomings of each method are discussed to provide a comprehensive view of the techniques used to study these structures. We aim to provide a set of guidelines for standardizing research on G-quadruplexes; these guidelines are not meant to be a dogmatic set of rules, but should rather provide useful information on the methods currently used to study these fascinating motifs.

Status does not predict stress among Hadza hunter-gatherer men.

In recent years there has been much research regarding the extent to which social status is related to long-term indices of health. The majority of studies looking at the interplay between social status and health have been conducted in industrialized societies. However, it has been argued that most of human evolution took place in small, mobile and egalitarian hunter-gatherer groups where individuals exhibited very little variation in terms of material wealth or possessions. In this study, we looked at the extent to which two domains of social status, hunting reputation (being perceived as a good hunter) and popularity (being perceived as a friend), are related to physiological stress levels among Hadza men, hunter-gatherers living in Northern Tanzania. The results of our study show that neither hunting reputation nor popularity is associated with stress levels. Overall, our data suggest that, in at least some traditional small-scale societies exhibiting an egalitarian social model, such as the Hadza, the variation in social status measures based on both popularity and hunting reputation does not translate into one of the commonly used indices of wellbeing.

Relationship between proximity and physiological stress levels in hunter-gatherers: The Hadza.

In recent years there has been a great deal of documentation on how social relationships are related to various aspects of human wellbeing. However, until recently most studies investigating the effects of social relationships on wellbeing have applied social network measures to reported social contacts. Recent advances in the application of bio-loggers in biological studies have now made it possible to quantify social relationships based on in-person, rather than self-reported, social interactions. We used GPS-derived in-camp and out-of-camp proximity data to analyse how in-person proximity is related to Hair Cortisol Concentration (HCC) among Hadza hunter-gatherers. Time spent in close proximity to other camp members was associated with higher HCC, especially in women. In contrast, individuals who spent more time in close out-of-camp proximity to their best friend experienced lower HCC. Our study suggests that physiological costs related to group living might be mitigated by in-person interactions with close friends. We also find that the location (i.e., in-camp vs out-of-camp) of proximity to others and self-perceived friends is associated with HCC among the Hadza.

The DEAD box RNA helicase DDX42 is an intrinsic inhibitor of positive-strand RNA viruses.

Genome-wide screens are powerful approaches to unravel regulators of viral infections. Here, a CRISPR screen identifies the RNA helicase DDX42 as an intrinsic antiviral inhibitor of HIV-1. Depletion of endogenous DDX42 increases HIV-1 DNA accumulation and infection in cell lines and primary cells. DDX42 overexpression inhibits HIV-1 infection, whereas expression of a dominant-negative mutant increases infection. Importantly, DDX42 also restricts LINE-1 retrotransposition and infection with other retroviruses and positive-strand RNA viruses, including CHIKV and SARS-CoV-2. However, DDX42 does not impact the replication of several negative-strand RNA viruses, arguing against an unspecific effect on target cells, which is confirmed by RNA-seq analysis. Proximity ligation assays show DDX42 in the vicinity of viral elements, and cross-linking RNA immunoprecipitation confirms a specific interaction of DDX42 with RNAs from sensitive viruses. Moreover, recombinant DDX42 inhibits HIV-1 reverse transcription in vitro. Together, our data strongly suggest a direct mode of action of DDX42 on viral ribonucleoprotein complexes. Our results identify DDX42 as an intrinsic viral inhibitor, opening new perspectives to target the life cycle of numerous RNA viruses.

Genome-wide mapping of individual replication fork velocities using nanopore sequencing.

Little is known about replication fork velocity variations along eukaryotic genomes, since reference techniques to determine fork speed either provide no sequence information or suffer from low throughput. Here we present NanoForkSpeed, a nanopore sequencing-based method to map and extract the velocity of individual forks detected as tracks of the thymidine analogue bromodeoxyuridine incorporated during a brief pulse-labelling of asynchronously growing cells. NanoForkSpeed retrieves previous Saccharomyces cerevisiae mean fork speed estimates (≈2 kb/min) in the BT1 strain exhibiting highly efficient bromodeoxyuridine incorporation and wild-type growth, and precisely quantifies speed changes in cells with altered replisome progression or exposed to hydroxyurea. The positioning of >125,000 fork velocities provides a genome-wide map of fork progression based on individual fork rates, showing a uniform fork speed across yeast chromosomes except for a marked slowdown at known pausing sites.

Effects of social tolerance on stress: hair cortisol concentrations in the tolerant Tonkean macaques (Macaca tonkeana) and the despotic long-tailed macaques (Macaca fascicularis).

Group living is a source of stress and an individuals' social environment has been shown to have a significant effect on its health and well-being. However, little is known about how different social organizations affect the stress levels of their members. Is living in a hierarchical society more or less stressful than living in a more tolerant structure? Here, we assess cortisol concentrations in the hair of two macaque species with radically different dominance styles: the egalitarian Tonkean macaque (Macaca Tonkeana) and the despotic long-tailed macaques (Macaca fascicularis). Hair was sampled in winter and again in late spring in two mixed-sex groups of 22 Tonkean macaques and 9 long-tailed macaques; Hair cortisol concentrations were significantly higher in the egalitarian Tonkean macaques than in the despotic long-tailed macaques, ranging from 161.13 to 938.8 pg/mg (mean ± SD 349.67 ± 126.22) and from 134.46 to 339.86 pg/mg (mean ± SD 231.2 ± 44.24), respectively. There was no difference between male and female cortisol concentrations, but hair cortisol increased with age in males. Dominance rank certainty was lower among female Tonkean macaques compared to long-tailed macaques. Our results suggest that species differences in dominance styles may translate into differences in long-term cortisol concentrations. We suggest that the higher cortisol concentrations in Tonkean macaques could be linked to the instability and lack of predictability and control around social relationships.

LAY SUMMARYBeing part of a social group can be very stressful, depending on the social structure of this group. We found that the more socially tolerant Tonkean macaques exhibited on average higher hair cortisol concentrations than more despotic long-tailed macaques. Males and females exhibited similar concentrations of hair cortisol in both species, but in male Tonkean macaques concentrations increased with age. The finding that overall cortisol levels were higher in the more tolerant species suggests that psychological arousal due to unpredictable social relations and mating competition may be an important driver of cortisol release in this species despite its overall tolerant social organization.

DNA molecular combing-based replication fork directionality profiling.

The replication strategy of metazoan genomes is still unclear, mainly because definitive maps of replication origins are missing. High-throughput methods are based on population average and thus may exclusively identify efficient initiation sites, whereas inefficient origins go undetected. Single-molecule analyses of specific loci can detect both common and rare initiation events along the targeted regions. However, these usually concentrate on positioning individual events, which only gives an overview of the replication dynamics. Here, we computed the replication fork directionality (RFD) profiles of two large genes in different transcriptional states in chicken DT40 cells, namely untranscribed and transcribed DMD and CCSER1 expressed at WT levels or overexpressed, by aggregating hundreds of oriented replication tracks detected on individual DNA fibres stretched by molecular combing. These profiles reconstituted RFD domains composed of zones of initiation flanking a zone of termination originally observed in mammalian genomes and were highly consistent with independent population-averaging profiles generated by Okazaki fragment sequencing. Importantly, we demonstrate that inefficient origins do not appear as detectable RFD shifts, explaining why dispersed initiation has remained invisible to population-based assays. Our method can both generate quantitative profiles and identify discrete events, thereby constituting a comprehensive approach to study metazoan genome replication.

Human ORC/MCM density is low in active genes and correlates with replication time but does not delimit initiation zones.

Eukaryotic DNA replication initiates during S phase from origins that have been licensed in the preceding G1 phase. Here, we compare ChIP-seq profiles of the licensing factors Orc2, Orc3, Mcm3, and Mcm7 with gene expression, replication timing, and fork directionality profiles obtained by RNA-seq, Repli-seq, and OK-seq. Both, the origin recognition complex (ORC) and the minichromosome maintenance complex (MCM) are significantly and homogeneously depleted from transcribed genes, enriched at gene promoters, and more abundant in early- than in late-replicating domains. Surprisingly, after controlling these variables, no difference in ORC/MCM density is detected between initiation zones, termination zones, unidirectionally replicating regions, and randomly replicating regions. Therefore, ORC/MCM density correlates with replication timing but does not solely regulate the probability of replication initiation. Interestingly, H4K20me3, a histone modification proposed to facilitate late origin licensing, was enriched in late-replicating initiation zones and gene deserts of stochastic replication fork direction. We discuss potential mechanisms specifying when and where replication initiates in human cells.

Thermal and pH Stabilities of i-DNA: Confronting in vitro Experiments with Models and In-Cell NMR Data.

Recent studies indicate that i-DNA, a four-stranded cytosine-rich DNA also known as the i-motif, is actually formed in vivo; however, a systematic study on sequence effects on stability has been missing. Herein, an unprecedented number of different sequences (271) bearing four runs of 3-6 cytosines with different spacer lengths has been tested. While i-DNA stability is nearly independent on total spacer length, the central spacer plays a special role on stability. Stability also depends on the length of the C-tracts at both acidic and neutral pHs. This study provides a global picture on i-DNA stability thanks to the large size of the introduced data set; it reveals unexpected features and allows to conclude that determinants of i-DNA stability do not mirror those of G-quadruplexes. Our results illustrate the structural roles of loops and C-tracts on i-DNA stability, confirm its formation in cells, and allow establishing rules to predict its stability.

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

Mechanisms of transcriptional control in malaria parasites are still not fully understood. The positioning patterns of G-quadruplex (G4) DNA motifs in the parasite's AT-rich genome, especially within the var gene family which encodes virulence factors, and in the vicinity of recombination hotspots, points towards a possible regulatory role of G4 in gene expression and genome stability. Here, we carried out the most comprehensive genome-wide survey, to date, of G4s in the Plasmodium falciparum genome using G4Hunter, which identifies G4 forming sequences (G4FS) considering their G-richness and G-skewness. We show an enrichment of G4FS in nucleosome-depleted regions and in the first exon of var genes, a pattern that is conserved within the closely related Laverania Plasmodium parasites. Under G4-stabilizing conditions, i.e., following treatment with pyridostatin (a high affinity G4 ligand), we show that a bona fide G4 found in the non-coding strand of var promoters modulates reporter gene expression. Furthermore, transcriptional profiling of pyridostatin-treated parasites, shows large scale perturbations, with deregulation affecting for instance the ApiAP2 family of transcription factors and genes involved in ribosome biogenesis. Overall, our study highlights G4s as important DNA secondary structures with a role in Plasmodium gene expression regulation, sub-telomeric recombination and var gene biology.

Author Correction: Involvement of G-quadruplex regions in mammalian replication origin activity.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

FORK-seq: replication landscape of the Saccharomyces cerevisiae genome by nanopore sequencing.

Genome replication mapping methods profile cell populations, masking cell-to-cell heterogeneity. Here, we describe FORK-seq, a nanopore sequencing method to map replication of single DNA molecules at 200-nucleotide resolution. By quantifying BrdU incorporation along pulse-chased replication intermediates from Saccharomyces cerevisiae, we orient 58,651 replication tracks reproducing population-based replication directionality profiles and map 4964 and 4485 individual initiation and termination events, respectively. Although most events cluster at known origins and fork merging zones, 9% and 18% of initiation and termination events, respectively, occur at many locations previously missed. Thus, FORK-seq reveals the full extent of cell-to-cell heterogeneity in DNA replication.

Status does not predict stress: Women in an egalitarian hunter-gatherer society.

It is widely believed that there is strong association between physiological stress and an individual's social status in their social hierarchy. This has been claimed for all humans cross-culturally, as well as in non-human animals living in social groups. However, the relationship between stress and social status has not been explored in any egalitarian hunter-gatherer society; it is also under investigated in exclusively female social groups. Most of human evolutionary history was spent in small, mobile foraging bands of hunter-gatherers with little economic differentiation - egalitarian societies. We analysed women's hair cortisol concentration along with two domains of women's social status (foraging reputation and popularity) in an egalitarian hunter-gatherer society, the Hadza. We hypothesized that higher social status would be associated with lower physiological indicators of stress in these women. Surprisingly, we did not find any association between either foraging reputation or popularity and hair cortisol concentration. The results of our study suggest that social status is not a consistent or powerful predictor of physiological stress levels in women in an egalitarian social structure. This challenges the notion that social status has the same basic physiological implications across all demographics and in all human societies.

Pharmacological characterisation of MDI-222, a novel AMPA receptor positive allosteric modulator with an improved safety profile.

PURPOSE: There is considerable interest in positive allosteric modulators (PAMs) of the α-amino-3-hydroxyl-5-methyl-4-isoxazole-propionate (AMPA) subtype of ionotropic glutamate receptors as therapeutic agents for a range of cognitive and mood disorders. However, the challenge is to increase AMPA receptor (AMPAR) function sufficient to enhance cognitive function but not to the extent that there are mechanism-related pro-convulsant or convulsant side effects. In this present study, we report the preclinical pharmacology data for MDI-222, an AMPAR PAM which enhances cognition but has a much reduced side-effect (i.e. convulsant) liability relative to other molecules of this mechanism. METHODS: The pharmacological effects of MDI-222 were characterised in in vitro and in vivo preclinical electrophysiology, efficacy (cognition), side-effect (pro-convulsant/convulsant), tolerability and toxicity assays. RESULTS: We demonstrate that MDI-222 is an AMPAR PAM, since it enhanced AMPAR function in vitro at human (hGluA1-4) and rat (rGluA2) homomeric receptors, and potentiated hetero-oligomeric AMPARs in rat neurons. MDI-222 enhanced electrically evoked AMPAR-mediated synaptic transmission in the anaesthetised rat at 10 mg/kg (administered intravenously) and did not significantly lower the seizure threshold in the pro-convulsant maximal electroshock threshold test (MEST) at any dose tested up to a maximum of 30 mg/kg (administered by oral gavage (p.o.)). MDI-222 reversed a delay-induced deficit in novel object recognition (NOR) in rats with a minimum effective dose (MED) of 0.3 mg/kg (p.o.) following acute administration, which was reduced to 0.1 mg/kg following sub-chronic administration, and improved passive avoidance performance in scopolamine-impaired rats with a MED of 10 mg/kg p.o. On the other hand, MDI-222 was not pro-convulsant in the MEST, resulting in a therapeutic window between plasma concentrations that enhanced cognitive performance and those associated with mechanism-related side effects of ⩾1000-fold. Unfortunately, despite the excellent preclinical profile of this compound, further development had to be halted due to non-mechanism-related issues. CONCLUSIONS: We conclude that MDI-222 is an AMPAR PAM which enhances cognitive performance in rats and has a significantly improved safety profile in preclinical species.

Involvement of G-quadruplex regions in mammalian replication origin activity.

Genome-wide studies of DNA replication origins revealed that origins preferentially associate with an Origin G-rich Repeated Element (OGRE), potentially forming G-quadruplexes (G4). Here, we functionally address their requirements for DNA replication initiation in a series of independent approaches. Deletion of the OGRE/G4 sequence strongly decreased the corresponding origin activity. Conversely, the insertion of an OGRE/G4 element created a new replication origin. This element also promoted replication of episomal EBV vectors lacking the viral origin, but not if the OGRE/G4 sequence was deleted. A potent G4 ligand, PhenDC3, stabilized G4s but did not alter the global origin activity. However, a set of new, G4-associated origins was created, whereas suppressed origins were largely G4-free. In vitro Xenopus laevis replication systems showed that OGRE/G4 sequences are involved in the activation of DNA replication, but not in the pre-replication complex formation. Altogether, these results converge to the functional importance of OGRE/G4 elements in DNA replication initiation.

Mapping and characterization of G-quadruplexes in the genome of the social amoeba Dictyostelium discoideum.

G-quadruplexes (G4) are non-canonical DNA and/or RNA secondary structures formed in guanine-rich regions. Given their over-representation in specific regions in the genome such as promoters and telomeres, they are likely to play important roles in key processes such as transcription, replication or RNA maturation. Putative G4-forming sequences (G4FS) have been reported in humans, yeast, bacteria, viruses and many organisms. Here we present the first mapping of G-quadruplex sequences in Dictyostelium discoideum, the social amoeba. 'Dicty' is an ameboid protozoan with a small (34 Mb) and extremely AT rich genome (78%). As a consequence, very few G4-prone motifs are expected. An in silico analysis of the Dictyostelium genome with the G4Hunter software detected 249-1055 G4-prone motifs, depending on G4Hunter chosen threshold. Interestingly, despite an even lower GC content (as compared to the whole Dicty genome), the density of G4 motifs in Dictyostelium promoters and introns is significantly higher than in the rest of the genome. Fourteen selected sequences located in important genes were characterized by a combination of biophysical and biochemical techniques. Our data show that these sequences form highly stable G4 structures under physiological conditions. Five Dictyostelium genes containing G4-prone motifs in their promoters were studied for the effect of a new G4-binding porphyrin derivative on their expression. Our results demonstrated that the new ligand significantly decreased their expression. Overall, our results constitute the first step to adopt Dictyostelium discoideum as a 'G4-poor' model for studies on G-quadruplexes.

G4HunterApps.

MOTIVATION: In order to help G4Hunter users and make it more accessible, I have developed a set of small applications within the Shiny/R framework. RESULTS: Each application fulfils simple tasks ranging from computing the G4Hunter score for a sequence or a list of sequence to extracting sequences with a G4Hunter score above a threshold for a sequence up to 5 Mb or a list of short sequences. The application can be installed either on the user computer within Rstudio or on a Rstudio server. AVAILABILITY AND IMPLEMENTATION: The source code for the ShinyApps is available on GitHub (https://github.com/LacroixLaurent).

Quadruplexes in 'Dicty': crystal structure of a four-quartet G-quadruplex formed by G-rich motif found in the Dictyostelium discoideum genome.

Guanine-rich DNA has the potential to fold into non-canonical G-quadruplex (G4) structures. Analysis of the genome of the social amoeba Dictyostelium discoideum indicates a low number of sequences with G4-forming potential (249-1055). Therefore, D. discoideum is a perfect model organism to investigate the relationship between the presence of G4s and their biological functions. As a first step in this investigation, we crystallized the dGGGGGAGGGGTACAGGGGTACAGGGG sequence from the putative promoter region of two divergent genes in D. discoideum. According to the crystal structure, this sequence folds into a four-quartet intramolecular antiparallel G4 with two lateral and one diagonal loops. The G-quadruplex core is further stabilized by a G-C Watson-Crick base pair and a A-T-A triad and displays high thermal stability (Tm > 90°C at 100 mM KCl). Biophysical characterization of the native sequence and loop mutants suggests that the DNA adopts the same structure in solution and in crystalline form, and that loop interactions are important for the G4 stability but not for its folding. Four-tetrad G4 structures are sparse. Thus, our work advances understanding of the structural diversity of G-quadruplexes and yields coordinates for in silico drug screening programs and G4 predictive tools.