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.

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.

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.

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.

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.

Chromatin immunoprecipitation indirect peaks highlight long-range interactions of insulator proteins and Pol II pausing.

Eukaryotic chromosomes are partitioned into topologically associating domains (TADs) that are demarcated by distinct insulator-binding proteins (IBPs) in Drosophila. Whether IBPs regulate specific long-range contacts and how this may impact gene expression remains unclear. Here we identify "indirect peaks" of multiple IBPs that represent their distant sites of interactions through long-range contacts. Indirect peaks depend on protein-protein interactions among multiple IBPs and their common cofactors, including CP190, as confirmed by high-resolution analyses of long-range contacts. Mutant IBPs unable to interact with CP190 impair long-range contacts as well as the expression of hundreds of distant genes that are specifically flanked by indirect peaks. Regulation of distant genes strongly correlates with RNAPII pausing, highlighting how this key transcriptional stage may trap insulator-based long-range interactions. Our data illustrate how indirect peaks may decipher gene regulatory networks through specific long-range interactions.

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.

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.

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.

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.

Re-evaluation of G-quadruplex propensity with G4Hunter.

Critical evidence for the biological relevance of G-quadruplexes (G4) has recently been obtained in seminal studies performed in a variety of organisms. Four-stranded G-quadruplex DNA structures are promising drug targets as these non-canonical structures appear to be involved in a number of key biological processes. Given the growing interest for G4, accurate tools to predict G-quadruplex propensity of a given DNA or RNA sequence are needed. Several algorithms such as Quadparser predict quadruplex forming propensity. However, a number of studies have established that sequences that are not detected by these tools do form G4 structures (false negatives) and that other sequences predicted to form G4 structures do not (false positives). Here we report development and testing of a radically different algorithm, G4Hunter that takes into account G-richness and G-skewness of a given sequence and gives a quadruplex propensity score as output. To validate this model, we tested it on a large dataset of 392 published sequences and experimentally evaluated quadruplex forming potential of 209 sequences using a combination of biophysical methods to assess quadruplex formation in vitro. We experimentally validated the G4Hunter algorithm on a short complete genome, that of the human mitochondria (16.6 kb), because of its relatively high GC content and GC skewness as well as the biological relevance of these quadruplexes near instability hotspots. We then applied the algorithm to genomes of a number of species, including humans, allowing us to conclude that the number of sequences capable of forming stable quadruplexes (at least in vitro) in the human genome is significantly higher, by a factor of 2-10, than previously thought.

New insights on the mechanism of quinoline-based DNA Methyltransferase inhibitors.

Among the epigenetic marks, DNA methylation is one of the most studied. It is highly deregulated in numerous diseases, including cancer. Indeed, it has been shown that hypermethylation of tumor suppressor genes promoters is a common feature of cancer cells. Because DNA methylation is reversible, the DNA methyltransferases (DNMTs), responsible for this epigenetic mark, are considered promising therapeutic targets. Several molecules have been identified as DNMT inhibitors and, among the non-nucleoside inhibitors, 4-aminoquinoline-based inhibitors, such as SGI-1027 and its analogs, showed potent inhibitory activity. Here we characterized the in vitro mechanism of action of SGI-1027 and two analogs. Enzymatic competition studies with the DNA substrate and the methyl donor cofactor, S-adenosyl-l-methionine (AdoMet), displayed AdoMet non-competitive and DNA competitive behavior. In addition, deviations from the Michaelis-Menten model in DNA competition experiments suggested an interaction with DNA. Thus their ability to interact with DNA was established; although SGI-1027 was a weak DNA ligand, analog 5, the most potent inhibitor, strongly interacted with DNA. Finally, as 5 interacted with DNMT only when the DNA duplex was present, we hypothesize that this class of chemical compounds inhibit DNMTs by interacting with the DNA substrate.

Fluorescence-based duplex-quadruplex competition test to screen for telomerase RNA quadruplex ligands.

RNA and DNA guanine-rich sequences can adopt unusual structures called Guanine quadruplexes (G4). A quadruplex-prone RNA sequence is present at the 5'-end of the 451-nt-long RNA component of telomerase, hTERC. As this quadruplex may interfere with P1 helix formation, a key structural element for this RNA, we are seeking molecules that would alter this RNA duplex-quadruplex equilibrium. In this work, we present a fluorescence-based test designed to identify G4 ligands specific for the hTERC G-rich motif and that can prevent P1 helix formation. From an initial panel of 169 different molecules, 11 were found to be excellent P1 duplex inhibitors. Interestingly, some of the compounds not only exhibit a strong selectivity for quadruplexes over duplexes, but also demonstrated a preference for G4-RNA over all other quadruplexes. This test may easily be adapted to almost any quadruplex-forming sequence and converted into HTS format.

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.

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.

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.

Sequence variant (CTAGGG)n in the human telomere favors a G-quadruplex structure containing a G.C.G.C tetrad.

Short contiguous arrays of variant CTAGGG repeats in the human telomere are unstable in the male germline and somatic cells, suggesting formation of unusual structures by this repeat type. Here, we report on the structure of an intramolecular G-quadruplex formed by DNA sequences containing four human telomeric variant CTAGGG repeats in potassium solution. Our results reveal a new robust antiparallel G-quadruplex fold involving two G-tetrads sandwiched between a G.C base pair and a G.C.G.C tetrad, which could represent a new platform for drug design targeted to human telomeric DNA.

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.