G4Hunter web application: a web server for G-quadruplex prediction

Computational design and experimental verification of pseudoknotted ribozymes

The design of new RNA sequences that retain the function of a model RNA structure is a challenge in bioinformatics because of the structural complexity of these molecules. RNA can fold into its secondary and tertiary structures by forming stem-loops and pseudoknots. A pseudoknot is a set of base pairs between a region within a stem-loop and nucleotides outside of this stem-loop; this motif is very important for numerous functional structures. It is important for any computational design algorithm to take into account these interactions to give a reliable result for any structures that include pseudoknots. In our study, we experimentally validated synthetic ribozymes designed by Enzymer, which implements algorithms allowing for the design of pseudoknots. Enzymer is a program that uses an inverse folding approach to design pseudoknotted RNAs; we used it in this study to design two types of ribozymes. The ribozymes tested were the hammerhead and the glmS, which have a self-cleaving activity that allows them to liberate the new RNA genome copy during rolling-circle replication or to control the expression of the downstream genes, respectively. We demonstrated the efficiency of Enzymer by showing that the pseudoknotted hammerhead and glmS ribozymes sequences it designed were extensively modified compared to wild-type sequences and were still active.

Presence, Location and Conservation of Putative G-Quadruplex Forming Sequences in Arboviruses Infecting Humans

Guanine quadruplexes (G4s) are non-canonical nucleic acid structures formed by guanine (G)-rich tracts that assemble into a core of stacked planar tetrads. G4s are found in the human genome and in the genomes of human pathogens, where they are involved in the regulation of gene expression and genome replication. G4s have been proposed as novel pharmacological targets in humans and their exploitation for antiviral therapy is an emerging research topic. Here, we report on the presence, conservation and localization of putative G4-forming sequences (PQSs) in human arboviruses. The prediction of PQSs was performed on more than twelve thousand viral genomes, belonging to forty different arboviruses that infect humans, and revealed that the abundance of PQSs in arboviruses is not related to the genomic GC content, but depends on the type of nucleic acid that constitutes the viral genome. Positive-strand ssRNA arboviruses, especially Flaviviruses, are significantly enriched in highly conserved PQSs, located in coding sequences (CDSs) or untranslated regions (UTRs). In contrast, negative-strand ssRNA and dsRNA arboviruses contain few conserved PQSs. Our analyses also revealed the presence of bulged PQSs, accounting for 17-26% of the total predicted PQSs. The data presented highlight the presence of highly conserved PQS in human arboviruses and present non-canonical nucleic acid-structures as promising therapeutic targets in arbovirus infections.

Dual mode of IFI16 binding to supercoiled and linear DNA: A closer insight

IFI16 (Interferon inducible protein 16) is a DNA sensor responsible for innate immune response stimulation and a direct viral restriction by modulating gene expression and replication. Many IFI16-DNA binding properties were described - length-dependent and sequence-independent binding, oligomerization of IFI16 upon recognition, sliding on the DNA, and preference for supercoiled DNA. However, the question of the role of IFI16-DNA binding in distinct IFI16 functions remains unclear. Here we demonstrate two modes of IFI16 binding to DNA using atomic force microscopy and electrophoretic mobility shift assays. In our study, we show that IFI16 can bind to DNA in the form of globular complexes or oligomers depending on DNA topology and molar ratios. The stability of the complexes is different in higher salt concentrations. In addition, we observed no preferential binding with the HIN-A or HIN-B domains to supercoiled DNA, revealing the importance of the whole protein for this specificity. These results provide more profound insight into IFI16-DNA interactions and may be important in answering the question of self- and non-self-DNA binding by the IFI16 protein and potentially could shed light on the role of DNA binding in distinct IFI16 functions.

Co-evolution of large inverted repeats and G-quadruplex DNA in fungal mitochondria may facilitate mitogenome stability: the case of Malassezia

Mitogenomes are essential due to their contribution to cell respiration. Recently they have also been implicated in fungal pathogenicity mechanisms. Members of the basidiomycetous yeast genus Malassezia are an important fungal component of the human skin microbiome, linked to various skin diseases, bloodstream infections, and they are increasingly implicated in gut diseases and certain cancers. In this study, the comparative analysis of Malassezia mitogenomes contributed to phylogenetic tree construction for all species. The mitogenomes presented significant size and gene order diversity which correlates to their phylogeny. Most importantly, they showed the inclusion of large inverted repeats (LIRs) and G-quadruplex (G4) DNA elements, rendering Malassezia mitogenomes a valuable test case for elucidating the evolutionary mechanisms responsible for this genome diversity. Both LIRs and G4s coexist and convergently evolved to provide genome stability through recombination. This mechanism is common in chloroplasts but, hitherto, rarely found in mitogenomes.

Dynamic alternative DNA structures in biology and disease

Repetitive elements in the human genome, once considered 'junk DNA', are now known to adopt more than a dozen alternative (that is, non-B) DNA structures, such as self-annealed hairpins, left-handed Z-DNA, three-stranded triplexes (H-DNA) or four-stranded guanine quadruplex structures (G4 DNA). These dynamic conformations can act as functional genomic elements involved in DNA replication and transcription, chromatin organization and genome stability. In addition, recent studies have revealed a role for these alternative structures in triggering error-generating DNA repair processes, thereby actively enabling genome plasticity. As a driving force for genetic variation, non-B DNA structures thus contribute to both disease aetiology and evolution.

A novel hypervariable variable number tandem repeat in the dopamine transporter gene (SLC6A3)

The dopamine transporter gene, SLC6A3, has received substantial attention in genetic association studies of various phenotypes. Although some variable number tandem repeats (VNTRs) present in SLC6A3 have been tested in genetic association studies, results have not been consistent. VNTRs in SLC6A3 that have not been examined genetically were characterized. The Tandem Repeat Annotation Library was used to characterize the VNTRs of 64 unrelated long-read haplotype-phased SLC6A3 sequences. Sequence similarity of each repeat unit of the five VNTRs is reported, along with the correlations of SNP-SNP, SNP-VNTR, and VNTR-VNTR alleles across the gene. One of these VNTRs is a novel hyper-VNTR (hyVNTR) in intron 8 of SLC6A3, which contains a range of 3.4-133.4 repeat copies and has a consensus sequence length of 38 bp, with 82% G+C content. The 38-base repeat was predicted to form G-quadruplexes in silico and was confirmed by circular dichroism spectroscopy. In addition, this hyVNTR contains multiple putative binding sites for PRDM9, which, in combination with low levels of linkage disequilibrium around the hyVNTR, suggests it might be a recombination hotspot.

Crosstalk between G-quadruplex and ROS

The excessive production of reactive oxygen species (ROS) can lead to single nucleic acid base damage, DNA strand breakage, inter- and intra-strand cross-linking of nucleic acids, and protein-DNA cross-linking involved in the pathogenesis of cancer, neurodegenerative diseases, and aging. G-quadruplex (G4) is a stacked nucleic acid structure that is ubiquitous across regulatory regions of multiple genes. Abnormal formation and destruction of G4s due to multiple factors, including cations, helicases, transcription factors (TFs), G4-binding proteins, and epigenetic modifications, affect gene replication, transcription, translation, and epigenetic regulation. Due to the lower redox potential of G-rich sequences and unique structural characteristics, G4s are highly susceptible to oxidative damage. Additionally, the formation, stability, and biological regulatory role of G4s are affected by ROS. G4s are involved in regulating gene transcription, translation, and telomere length maintenance, and are therefore key players in age-related degeneration. Furthermore, G4s also mediate the antioxidant process by forming stress granules and activating Nrf2, which is suggestive of their involvement in developing ROS-related diseases. In this review, we have summarized the crosstalk between ROS and G4s, and the possible regulatory mechanisms through which G4s play roles in aging and age-related diseases.

Analysis of G-Quadruplex-Forming Sequences in Drought Stress-Responsive Genes, and Synthesis Genes of Phenolic Compounds in Arabidopsis thaliana

Sequences of nucleic acids with the potential to form four-stranded G-quadruplex structures are intensively studied mainly in the context of human diseases, pathogens, or extremophile organisms; nonetheless, the knowledge about their occurrence and putative role in plants is still limited. This work is focused on G-quadruplex-forming sites in two gene sets of interest: drought stress-responsive genes, and genes related to the production/biosynthesis of phenolic compounds in the model plant organism Arabidopsis thaliana. In addition, 20 housekeeping genes were analyzed as well, where the constitutive gene expression was expected (with no need for precise regulation depending on internal or external factors). The results have shown that none of the tested gene sets differed significantly in the content of G-quadruplex-forming sites, however, the highest frequency of G-quadruplex-forming sites was found in the 5'-UTR regions of phenolic compounds' biosynthesis genes, which indicates the possibility of their regulation at the mRNA level. In addition, mainly within the introns and 1000 bp flanks downstream gene regions, G-quadruplex-forming sites were highly underrepresented. Finally, cluster analysis allowed us to observe similarities between particular genes in terms of their PQS characteristics. We believe that the original approach used in this study may become useful for further and more comprehensive bioinformatic studies in the field of G-quadruplex genomics.

Beyond the Primary Structure of Nucleic Acids: Potential Roles of Epigenetics and Noncanonical Structures in the Regulations of Plant Growth and Stress Responses

Epigenetics deals with changes in gene expression that are not caused by modifications in the primary sequence of nucleic acids. These changes beyond primary structures of nucleic acids not only include DNA/RNA methylation, but also other reversible conversions, together with histone modifications or RNA interference. In addition, under particular conditions (such as specific ion concentrations or protein-induced stabilization), the right-handed double-stranded DNA helix (B-DNA) can form noncanonical structures commonly described as "non-B DNA" structures. These structures comprise, for example, cruciforms, i-motifs, triplexes, and G-quadruplexes. Their formation often leads to significant differences in replication and transcription rates. Noncanonical RNA structures have also been documented to play important roles in translation regulation and the biology of noncoding RNAs. In human and animal studies, the frequency and dynamics of noncanonical DNA and RNA structures are intensively investigated, especially in the field of cancer research and neurodegenerative diseases. In contrast, noncanonical DNA and RNA structures in plants have been on the fringes of interest for a long time and only a few studies deal with their formation, regulation, and physiological importance for plant stress responses. Herein, we present a review focused on the main fields of epigenetics in plants and their possible roles in stress responses and signaling, with special attention dedicated to noncanonical DNA and RNA structures.

Inverted repeats in the monkeypox virus genome are hot spots for mutation

The current monkeypox virus (MPXV) strain differs from the strain arising in 2018 by 50+ single nucleotide polymorphisms (SNPs) and is mutating much faster than expected. The cytidine deaminase apolipoprotein B messenger RNA editing enzyme, catalytic subunit B (APOBEC3) was hypothesized to be driving this increased mutation. APOBEC has recently been identified to preferentially mutate cruciform DNA secondary structures formed by inverted repeats (IRs). IRs were recently identified as hot spots for mutation in severe acute respiratory syndrome coronavirus 2, and we aimed to identify whether IRs were also hot spots for mutation within MPXV genomes. We found that MPXV genomes were replete with IR sequences. Of the 50+ SNPs identified in the 2022 outbreak strain, 63.9% of these were found to have arisen within IR regions in the 2018 reference strain (MT903344.1). Notably, IR sequences found in the 2018 reference strain were significantly lost over time, with an average of 32.5% of these sequences being conserved in the 2022 MPXV genomes. This evidence was highly indicative that mutations were arising within IRs. This data provides further support to the hypothesis that APOBEC may be driving MPXV mutation and highlights the necessity for greater surveillance of IRs of MPXV genomes to detect new mutations.

Targeting G-quadruplex for rescuing impaired chondrogenesis in WRN-deficient stem cells

BACKGROUND

Pathogenic mutations in WRN are a cause of premature aging disease Werner syndrome (WS). Besides accelerated aging phenotypes and cancer predisposition, patients with WS also display underdevelopment in the skeletal system, characterized by short stature, light body weight and unusually thin extremities. The reasons for these developmental defects are not completely understood and the underlying molecular mechanism remains to be elucidated.

RESULTS

In this study, WRN was found to modulate transcription of short stature homeobox gene SHOX. Loss of WRN resulted in insufficient expression of SHOX, the gene dose of which is critical for driving chondrocyte differentiation. WRN could bind the G-quadruplex (G4) structures in the SHOX promoter and stimulate transcription. Aberrant formation of G4 structures in WRN-deficient cells impeded normal transcription of SHOX, thus resulting in impaired chondrogenesis. Chondrogenesis could be rescued by overexpression of WRN helicase or SHOX, suggesting that SHOX is a downstream target of WRN. Gene editing of the G4 structures in the SHOX promoter could increase SHOX expression, therefore rescuing the impaired chondrogenesis in WRN-deficient cells.

CONCLUSIONS

Our data suggest that dysgenesis of the developing bone in WS might be caused by SHOX insufficiency. Aberrant formation of G4 structures in SHOX promoter suppresses SHOX expression and impairs chondrogenesis. Targeted mutagenesis in the G4 structures enhances SHOX expression and thus providing an opportunity to rescue the chondrogenic defect.

Unlocking G-Quadruplexes as Targets and Tools against COVID-19

The applicability of G-quadruplexes (G4s) as antiviral targets, therapeutic agents and diagnostic tools for coronavirus disease 2019 (COVID-19) is currently being evaluated, which has drawn the extensive attention of the scientific community. During the COVID-19 pandemic, research in this field is rapidly accumulating. In this review, we summarize the latest achievements and breakthroughs in the use of G4s as antiviral targets, therapeutic agents and diagnostic tools for COVID-19, particularly using G4 ligands. Finally, strength and weakness regarding G4s in anti-SARS-CoV-2 field are highlighted for prospective future projects.

Dysregulation of Plasma miR-146a and miR-155 Expression Profile in Mycosis Fungoides Is Associated with rs2910164 and rs767649 Polymorphisms

Diagnosis of Mycosis Fungoides (MF) may be challenging, due to its polymorphic nature. The use of miRNAs as biomarkers to assist in diagnosis has been investigated, mainly in skin lesion biopsies. The purpose of this study is to evaluate the plasma levels of miR-146a and miR-155 in MF patients and to investigate their association with SNPs of their genes. Plasma miRNAs were quantified by RT-qPCR. Genomic DNA was used for SNPs’ genotyping by Sanger sequencing. Plasma levels of miR-146a and miR-155 were significantly higher in patients vs. controls, in early MF patients vs. controls, and in advanced vs. early MF patients. Both miRNAs’ levels were significantly higher in stage IIB vs. early-stage patients. miR-155 plasma levels were significantly higher in patients with skin tumors or erythroderma. CC genotype (rs2910164 C>G) was significantly more frequent in healthy controls and associated with lower MF risk and lower miR-146a levels. The AA genotype (rs767649 T>A) was significantly more frequent in patients and correlated with increased MF risk and increased miR-155 levels. The combination of GG+AA was only detected in patients and was correlated with higher MF susceptibility. Increased mir-146a and mir-155 plasma levels in MF is an important finding to establish putative noninvasive biomarkers. The presence of SNPs is closely associated with miRs’ expression, and possibly with disease susceptibility.

Interaction of C-terminal p53 isoforms depends strongly upon DNA sequence and topology

The p53 protein is a key tumor suppressor and the most commonly mutated and down-regulated protein in human tumors. It functions mainly through interaction with DNA, and p53 acts as a transcription factor that recognizes the so-called p53 target sites on the promoters of various genes. P53 has been shown to exist as many isoforms, including three C-terminal isoforms that are produced by alternative splicing. Because the C-terminal domain is responsible for sequence-nonspecific binding and regulation of p53 binding, we have analyzed DNA recognition by these C-terminal isoforms. Using atomic force microscopy, we show for the first time that all C-terminal isoforms recognize superhelical DNA. It is particularly noteworthy that a sequence-specific p53 consensus binding site is bound by p53α and β isoforms with similar affinities, whilst p53α shows higher binding to a quadruplex sequence than both p53β and p53γ, and p53γ loses preferential binding to both the consensus binding sequence and the quadruplex-forming sequence. These results show the important role of the variable p53 C-terminal amino acid sequences for DNA recognition.

Deciphering RNA G-quadruplex function during the early steps of HIV-1 infection

G-quadruplexes (G4s) are four-stranded nucleic acid structures formed by the stacking of G-tetrads. Here we investigated their formation and function during HIV-1 infection. Using bioinformatics and biophysics analyses we first searched for evolutionary conserved G4-forming sequences in HIV-1 genome. We identified 10 G4s with conservation rates higher than those of HIV-1 regulatory sequences such as RRE and TAR. We then used porphyrin-based G4-binders to probe the formation of the G4s during infection of human cells by native HIV-1. The G4-binders efficiently inhibited HIV-1 infectivity, which is attributed to the formation of G4 structures during HIV-1 replication. Using a qRT-PCR approach, we showed that the formation of viral G4s occurs during the first 2 h post-infection and their stabilization by the G4-binders prevents initiation of reverse transcription. We also used a G4-RNA pull-down approach, based on a G4-specific biotinylated probe, to allow the direct detection and identification of viral G4-RNA in infected cells. Most of the detected G4-RNAs contain crucial regulatory elements such as the PPT and cPPT sequences as well as the U3 region. Hence, these G4s would function in the early stages of infection when the viral RNA genome is being processed for the reverse transcription step.

QUADRatlas: the RNA G-quadruplex and RG4-binding proteins database

RNA G-quadruplexes (RG4s) are non-canonical, disease-associated post-transcriptional regulators of gene expression whose functions are driven by RNA-binding proteins (RBPs). Being able to explore transcriptome-wide RG4 formation and interaction with RBPs is thus paramount to understanding how they are regulated and exploiting them as potential therapeutic targets. Towards this goal, we present QUADRatlas (https://rg4db.cibio.unitn.it), a database of experimentally-derived and computationally predicted RG4s in the human transcriptome, enriched with biological function and disease associations. As RBPs are key to their function, we mined known interactions of RG4s with such proteins, complemented with an extensive RBP binding sites dataset. Users can thus intersect RG4s with their potential regulators and effectors, enabling the formulation of novel hypotheses on RG4 regulation, function and pathogenicity. To support this capability, we provide analysis tools for predicting whether an RBP can bind RG4s, RG4 enrichment in a gene set, and de novo RG4 prediction. Genome-browser and table views allow exploring, filtering, and downloading the data quickly for individual genes and in batch. QUADRatlas is a significant step forward in our ability to understand the biology of RG4s, offering unmatched data content and enabling the integrated analysis of RG4s and their interactions with RBPs.

Telomeres and Their Neighbors

Telomeres are essential structures formed from satellite DNA repeats at the ends of chromosomes in most eukaryotes. Satellite DNA repeat sequences are useful markers for karyotyping, but have a more enigmatic role in the eukaryotic cell. Much work has been done to investigate the structure and arrangement of repetitive DNA elements in classical models with implications for species evolution. Still more is needed until there is a complete picture of the biological function of DNA satellite sequences, particularly when considering non-model organisms. Celebrating Gregor Mendel’s anniversary by going to the roots, this review is designed to inspire and aid new research into telomeres and satellites with a particular focus on non-model organisms and accessible experimental and in silico methods that do not require specialized equipment or expensive materials. We describe how to identify telomere (and satellite) repeats giving many examples of published (and some unpublished) data from these techniques to illustrate the principles behind the experiments. We also present advice on how to perform and analyse such experiments, including details of common pitfalls. Our examples are a selection of recent developments and underexplored areas of research from the past. As a nod to Mendel’s early work, we use many examples from plants and insects, especially as much recent work has expanded beyond the human and yeast models traditional in telomere research. We give a general introduction to the accepted knowledge of telomere and satellite systems and include references to specialized reviews for the interested reader.

Identification and structural analysis of novel malathion-specific DNA aptameric sensors designed for food testing

Malathion is an organophosphate chemical (OPC) and a toxic contaminant that adversely impacts food quality, human health, biodiversity, and the environment. Due to its small size and unavailability of sensitive sensors, detection of malathion remains a challenging task. Often chromatographic methods employed to analyze OPCs suffer from several shortcomings, including cost, immobility, laboriousness, and unsuitability for point-of-care settings. Hence, developing a specific and sensitive diagnostic sensor for quick and inexpensive food testing is essential. We discovered four unique malathion-specific ssDNA aptamers; designed two independent sensing strategies using fluorescence labeling and Thioflavin T (ThT) displacement. Selected aptamers formed the G4-quadruplex-like (G4Q) structure, which helped develop a label-free detection approach with a 2.01 ppb limit of detection. Additionally, 3D structures of aptamers were generated and validated using a series of computational modeling programs. Furthermore, we explored structural features using CD spectroscopy and molecular docking, probing ligands' binding mode, and revealed vital intermolecular interactions with aptamers. Subsequently, the novel sensors were optimized to detect malathion from food samples. The novel sensors could be further developed to meet the demands of sensing and quantifying toxic contaminants from real food samples in field conditions.

The Newly Sequenced Genome of Pisum sativum Is Replete with Potential G-Quadruplex-Forming Sequences-Implications for Evolution and Biological Regulation

G-quadruplexes (G4s) have been long considered rare and physiologically unimportant in vitro curiosities, but recent methodological advances have proved their presence and functions in vivo. Moreover, in addition to their functional relevance in bacteria and animals, including humans, their importance has been recently demonstrated in evolutionarily distinct plant species. In this study, we analyzed the genome of Pisum sativum (garden pea, or the so-called green pea), a unique member of the Fabaceae family. Our results showed that this genome contained putative G4 sequences (PQSs). Interestingly, these PQSs were located nonrandomly in the nuclear genome. We also found PQSs in mitochondrial (mt) and chloroplast (cp) DNA, and we experimentally confirmed G4 formation for sequences found in these two organelles. The frequency of PQSs for nuclear DNA was 0.42 PQSs per thousand base pairs (kbp), in the same range as for cpDNA (0.53/kbp), but significantly lower than what was found for mitochondrial DNA (1.58/kbp). In the nuclear genome, PQSs were mainly associated with regulatory regions, including 5'UTRs, and upstream of the rRNA region. In contrast to genomic DNA, PQSs were located around RNA genes in cpDNA and mtDNA. Interestingly, PQSs were also associated with specific transposable elements such as TIR and LTR and around them, pointing to their role in their spreading in nuclear DNA. The nonrandom localization of PQSs uncovered their evolutionary and functional significance in the Pisum sativum genome.

Short LNA-modified oligonucleotide probes as efficient disruptors of DNA G-quadruplexes

G-quadruplexes (G4s) are well known non-canonical DNA secondary structures that can form in human cells. Most of the tools available to investigate G4-biology rely on small molecule ligands that stabilise these structures. However, the development of probes that disrupt G4s is equally important to study their biology. In this study, we investigated the disruption of G4s using Locked Nucleic Acids (LNA) as invader probes. We demonstrated that strategic positioning of LNA-modifications within short oligonucleotides (10 nts.) can significantly accelerate the rate of G4-disruption. Single-molecule experiments revealed that short LNA-probes can promote disruption of G4s with mechanical stability sufficient to stall polymerases. We corroborated this using a single-step extension assay, revealing that short LNA-probes can relieve replication dependent polymerase-stalling at G4 sites. We further demonstrated the potential of such LNA-based probes to study G4-biology in cells. By using a dual-luciferase assay, we found that short LNA probes can enhance the expression of c-KIT to levels similar to those observed when the c-KIT promoter is mutated to prevent the formation of the c-KIT1 G4. Collectively, our data suggest a potential use of rationally designed LNA-modified oligonucleotides as an accessible chemical-biology tool for disrupting individual G4s and interrogating their biological functions in cells.

Small molecule-based detection of non-canonical RNA G-quadruplex structures that modulate protein translation

Tandem repeats of guanine-rich sequences in RNA often form thermodynamically stable four-stranded RNA structures. Such RNA G-quadruplexes have long been considered to be linked to essential biological processes, yet their physiological significance in cells remains unclear. Here, we report a approach that permits the detection of RNA G-quadruplex structures that modulate protein translation in mammalian cells. The approach combines antibody arrays and RGB-1, a small molecule that selectively stabilizes RNA G-quadruplex structures. Analysis of the protein and mRNA products of 84 cancer-related human genes identified Nectin-4 and CapG as G-quadruplex-controlled genes whose mRNAs harbor non-canonical G-quadruplex structures on their 5′UTR region. Further investigations revealed that the RNA G-quadruplex of CapG exhibits a structural polymorphism, suggesting a possible mechanism that ensures the translation repression in a KCl concentration range of 25–100 mM. The approach described in the present study sets the stage for further discoveries of RNA G-quadruplexes.

G4detector: Convolutional Neural Network to Predict DNA G-Quadruplexes

G-quadruplexes (G4s) are nucleic acid secondary structures that form within guanine-rich DNA or RNA sequences. G4 formation can affect chromatin architecture and gene regulation, and has been associated with genomic instability, genetic diseases, and cancer progression. The experimental data produced by the G4-seq experiment provides unprecedented details on G4 formation in the genome. Still, running the experimental protocol on a whole genome is an expensive and time-consuming process. Thus, it is highly desirable to have a computational method to predict G4 formation in new DNA sequences or whole genomes. Here, we present G4detector, a new method based on a convolutional neural network to predict G4s from DNA sequences. On top of the sequence information, we improved prediction accuracy by the addition of RNA secondary structure information. To train and test G4detector, we compiled novel high-throughput benchmarks over multiple species genomes measured by the G4-seq protocol. We show that G4detector outperforms extant methods for the same task on all benchmark datasets, can detect G4s genome-wide with high accuracy, and is able to extrapolate human-trained measurements to various non-human species. The code and benchmarks are publicly available on github.com/OrensteinLab/G4detector.

Ni+2 permease system of Helicobacter pylori contains highly conserved G-quadruplex motifs

The genome of a micro-organism contains all the information required for its survival inside its host cells. The guanine rich regions of the genome can form stable G-quadruplex structures that act as the regulators of gene expression. Herein, the completely sequenced genomes of Helicobacter pylori were explored for the identification and characterization of the conserved G-quadruplex motifs in this gastrointestinal pathogen. Initial in silico analysis revealed the presence of ~8241 GQ motifs in the H. pylori genome. Metal binding proteins of H. pylori are significantly enriched in the GQ motifs. Our study emphasizes the identification and characterization of four highly conserved G-quadruplex forming motifs (HPGQs) in the nickel transporter genes (nixA, niuB1, niuB2, and niuD) of the H. pylori. Nickel is a virulence determinant in H. pylori and is required as a co-factor for the urease and [NiFe] hydrogenase enzymes that are crucial for its survival in the stomach lining of humans. The presence of GQ motifs in these nickel transporter genes can affect their expression and may alter the functioning of Urease and [NiFe] hydrogenase. Similar to human and virus G-quadruplexes, targeting these conserved PGQs with bioactive molecules may represent a novel therapeutic avenue for combating infection of H. pylori. The identified HPGQs were characterized in-vitro by using CD spectroscopy, electrophoresis technique, and NMR spectroscopy at both acidic (4.5) and neutral pH (7.0). ITC revealed the specific interaction of these HPGQs with high affinity to the known G-quadruplex binding ligand, TMPyP4. The mTFP based reporter assay showed decrease in the gene expression of mTFP in the TMPyP4 treated cells as compared to the untreated and further affirmed the formation of stable G-quadruplex structures in the HPGQ motifs in vivo. This is the first report for characterizing G-quadruplex motifs in nickel transport-associated genes in the H. pylori bacterium.

Selection and characterization of DNA aptamers inhibiting a druggable target of osteoarthritis, ADAMTS-5

There is a critical need for the development of more potent inhibitors for osteoarthritis (OA) therapy given the poor life quality of arthritis patients. Aggrecanase ADAMTS-5 (a disintegrin and metalloproteinase with thrombospondin motifs 5) is an established drug target identified for osteoarthritis. In this study, we evolved and characterized two new DNA aptamer inhibitors of ADAMTS-5, namely apt21 and apt25. The aptamers exhibited nanomolar binding affinity and high specificity against ADAMTS-5. K_D values of apt21 and apt25 were determined by the Enzyme-linked Oligonucleotide Assay (ELONA) at 1.54 ± 0.16 nM and 1.79 ± 0.08 nM, respectively. Circular Dichroism (CD) analysis demonstrated that both aptamers formed monovalent cation dependent G-quadruplex structures. Calcium ions did not affect the binding of the aptamers to ADAMTS-5. The inhibitory effects of apt21 and apt25 on ADAMTS-5 were evaluated by the Förster Resonance Energy Transfer (FRET) assay, in which IC₅₀ values of apt21 and apt25 were estimated at 52.76 ± 6.70 μM and 61.14 ± 9.67 μM, respectively. These two aptamers are the first DNA G-quadruplex aptamers demonstrated to inhibit ADAMTS-5 and could have value for OA therapy.

Topoisomerase 1 inhibits MYC promoter activity by inducing G-quadruplex formation

We have investigated the function of human topoisomerase 1 (TOP1) in regulation of G-quadruplex (G4) formation in the Pu27 region of the MYC P1 promoter. Pu27 is among the best characterized G4 forming sequences in the human genome and it is well known that promoter activity is inhibited upon G4 formation in this region. We found that TOP1 downregulation stimulated transcription from a promoter with wildtype Pu27 but not if the G4 motif in Pu27 was interrupted by mutation(s). The effect was not specific to the MYC promoter and similar results were obtained for the G4 forming promoter element WT21. The other major DNA topoisomerases with relaxation activity, topoisomerases 2α and β, on the other hand, did not affect G4 dependent promoter activity. The cellular studies were supported by in vitro investigations demonstrating a high affinity of TOP1 for wildtype Pu27 but not for mutant sequences unable to form G4. Moreover, TOP1 was able to induce G4 formation in Pu27 inserted in double stranded plasmid DNA in vitro. This is the first time TOP1 has been demonstrated capable of inducing G4 formation in double stranded DNA and of influencing G4 formation in cells.

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.

Chromosome organization affects genome evolution in Sulfolobus archaea

In all organisms, the DNA sequence and the structural organization of chromosomes affect gene expression. The extremely thermophilic crenarchaeon Sulfolobus has one circular chromosome with three origins of replication. We previously revealed that this chromosome has defined A and B compartments that have high and low gene expression, respectively. As well as higher levels of gene expression, the A compartment contains the origins of replication. To evaluate the impact of three-dimensional organization on genome evolution, we characterized the effect of replication origins and compartmentalization on primary sequence evolution in eleven Sulfolobus species. Using single-nucleotide polymorphism analyses, we found that distance from an origin of replication was associated with increased mutation rates in the B but not in the A compartment. The enhanced polymorphisms distal to replication origins suggest that replication termination may have a causal role in their generation. Further mutational analyses revealed that the sequences in the A compartment are less likely to be mutated, and that there is stronger purifying selection than in the B compartment. Finally, we applied the Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq) to show that the B compartment is less accessible than the A compartment. Taken together, our data suggest that compartmentalization of chromosomal DNA can influence chromosome evolution in Sulfolobus. We propose that the A compartment serves as a haven for stable maintenance of gene sequences, while sequences in the B compartment can be diversified.

Development of RNA G-quadruplex (rG4)-targeting L-RNA aptamers by rG4-SELEX

RNA G-quadruplex (rG4)-SELEX is a method that generates L-RNA aptamers to target an rG4 structure of interest, which can be applied to inhibit G-quadruplex-mediated interactions that have important roles in gene regulation and function. Here we present a Protocol Extension substantially modifying an existing SELEX protocol to describe in detail the procedures involved in performing rG4-SELEX to identify rG4-specific binders that can effectively suppress rG4-peptide and rG4-protein associations. This Protocol Extension improves the speed of aptamer discovery and identification, offering a suite of techniques to characterize the aptamer secondary structure and monitor binding affinity and specificity, and demonstrating the utility of the L-RNA aptamer. The previous protocol mainly describes the identification of RNA aptamers against proteins of interest, whereas in this Protocol Extension we present the development of an unnatural RNA aptamer against an RNA structure of interest, with the potential to be applicable to other nucleic acid motifs or biomolecules. rG4-SELEX starts with a random D-RNA library incubated with the L-rG4 target of interest, followed by binding, washing and elution of the library. Enriched D-aptamer candidates are sequenced and structurally characterized. Then, the L-aptamer is synthesized and used for different applications. rG4-SELEX can be carried out by an experienced molecular biologist with a basic understanding of nucleic acids. The development of rG4-targeting L-RNA aptamers expands the current rG4 toolkit to explore innovative rG4-related applications, and opens new doors to discovering novel rG4 biology in the near future. The duration of each selection cycle as outlined in the protocol is ~2 d.

Development of novel fluorescence-based and label-free noncanonical G4-quadruplex-like DNA biosensor for facile, specific, and ultrasensitive detection of fipronil

Fipronil is a broad-spectrum insecticide widely used in agriculture and residential areas; its indiscriminate use leads to environmental pollution and poses health hazards. Early detection of fipronil is critical to prevent the deleterious effects. However, current insecticide analysis methods such as HPLC, LC/MS, and GC/MS are incompetent; they are costly, immobile, time-consuming, laborious, and need skilled technicians. Hence, a sensitive, specific, and cheap biosensor are essential to containing the contamination. Here, we designed two novel biosensors-the first design relied on fluorescent labeling/quenching, while the second sensor focused on label-free detection using Thioflavin T displacement. Altogether, we identified four candidate aptamers, predicted secondary structures, and performed 3D molecular modeling to predict the binding pocket of fipronil in FiPA6B aptamer. Furthermore, the aptameric sensors showed high sensitivity to fipronil of sub-ppb level LOD, attributed to stringent experimental design. The biosensors displayed high specificity against other phenylpyrazole insecticides and demonstrated robust sensitivity for fipronil in real samples like cabbage and cucumber. Notably, to the best of our knowledge, this is the first demonstration of noncanonical G4-quadruplex-like aptamer binding to fipronil, verified using CD spectroscopy. Such aptasensors possess considerable potential for real-time measurements of hazardous insecticides as point-of-care technology.

From telomere to telomere: The transcriptional and epigenetic state of human repeat elements

Mobile elements and repetitive genomic regions are sources of lineage-specific genomic innovation and uniquely fingerprint individual genomes. Comprehensive analyses of such repeat elements, including those found in more complex regions of the genome, require a complete, linear genome assembly. We present a de novo repeat discovery and annotation of the T2T-CHM13 human reference genome. We identified previously unknown satellite arrays, expanded the catalog of variants and families for repeats and mobile elements, characterized classes of complex composite repeats, and located retroelement transduction events. We detected nascent transcription and delineated CpG methylation profiles to define the structure of transcriptionally active retroelements in humans, including those in centromeres. These data expand our insight into the diversity, distribution, and evolution of repetitive regions that have shaped the human genome.

A new lineage of non-photosynthetic green algae with extreme organellar genomes

Background

The plastid genomes of the green algal order Chlamydomonadales tend to expand their non-coding regions, but this phenomenon is poorly understood. Here we shed new light on organellar genome evolution in Chlamydomonadales by studying a previously unknown non-photosynthetic lineage. We established cultures of two new Polytoma-like flagellates, defined their basic characteristics and phylogenetic position, and obtained complete organellar genome sequences and a transcriptome assembly for one of them.

Results

We discovered a novel deeply diverged chlamydomonadalean lineage that has no close photosynthetic relatives and represents an independent case of photosynthesis loss. To accommodate these organisms, we establish the new genus Leontynka, with two species (L. pallida and L. elongata) distinguishable through both their morphological and molecular characteristics. Notable features of the colourless plastid of L. pallida deduced from the plastid genome (plastome) sequence and transcriptome assembly include the retention of ATP synthase, thylakoid-associated proteins, the carotenoid biosynthesis pathway, and a plastoquinone-based electron transport chain, the latter two modules having an obvious functional link to the eyespot present in Leontynka. Most strikingly, the ~362 kbp plastome of L. pallida is by far the largest among the non-photosynthetic eukaryotes investigated to date due to an extreme proliferation of sequence repeats. These repeats are also present in coding sequences, with one repeat type found in the exons of 11 out of 34 protein-coding genes, with up to 36 copies per gene, thus affecting the encoded proteins. The mitochondrial genome of L. pallida is likewise exceptionally large, with its >104 kbp surpassed only by the mitogenome of Haematococcus lacustris among all members of Chlamydomonadales hitherto studied. It is also bloated with repeats, though entirely different from those in the L. pallida plastome, which contrasts with the situation in H. lacustris where both the organellar genomes have accumulated related repeats. Furthermore, the L. pallida mitogenome exhibits an extremely high GC content in both coding and non-coding regions and, strikingly, a high number of predicted G-quadruplexes.

Conclusions

With its unprecedented combination of plastid and mitochondrial genome characteristics, Leontynka pushes the frontiers of organellar genome diversity and is an interesting model for studying organellar genome evolution.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01263-w.

The Ribosomal DNA Loci of the Ancient Monocot Pistia stratiotes L. (Araceae) Contain Different Variants of the 35S and 5S Ribosomal RNA Gene Units

The freshwater plant water lettuce (Pistia stratiotes L.) grows in warm climatic zones and is used for phytoremediation and biomass production. P. stratiotes belongs to the Araceae, an ecologically and structurally diverse early monocot family, but the phylogenetic relationships among Araceae members are poorly understood. Ribosomal DNAs (rDNAs), including the 35S and 5S rDNA, encode the RNA components of ribosomes and are widely used in phylogenetic and evolutionary studies of various plant taxa. Here, we comprehensively characterized the chromosomal locations and molecular organization of 35S and 5S rDNA genes in water lettuce using karyological and molecular methods. Fluorescence in situ hybridization revealed a single location for the 35S and 5S rDNA loci, each on a different pair of the species' 28 chromosomes. Molecular cloning and nucleotide sequencing of 35S rDNA of P. stratiotes, the first representative Araceae sensu stricto in which such a study was performed, displayed typical structural characteristics. The full-length repeat showed high sequence conservation of the regions producing the 18S, 5.8S, and 25S rRNAs and divergence of the internal transcribed spacers ITS1 and ITS2 as well as the large intergenic spacer (IGS). Alignments of the deduced sequence of 18S rDNA with the sequences available for other Araceae and representatives of other clades were used for phylogenetic analysis. Examination of 11 IGS sequences revealed significant intra-genomic length variability due to variation in subrepeat number, with four types of units detected within the 35S rDNA locus of the P. stratiotes genome (estimated size 407 Mb/1C). Similarly, the 5S rDNA locus harbors gene units comprising a conserved 119-bp sequence encoding 5S rRNA and two types of non-transcribed spacer (NTS) sequences. Type I was classified into four subtypes, which apparently originated via progressive loss of subrepeats within the duplicated NTS region containing the 3' part of the 5S rRNA gene. The minor Type II NTS is shorter than Type I and differs in nucleotide composition. Some DNA clones containing two or three consecutive 5S rDNA repeats harbored 5S rDNA genes with different types of NTSs, confirming the mosaic composition of the 5S rDNA locus.

G-quadruplexes in helminth parasites

Parasitic helminths infecting humans are highly prevalent infecting ∼2 billion people worldwide, causing inflammatory responses, malnutrition and anemia that are the primary cause of morbidity. In addition, helminth infections of cattle have a significant economic impact on livestock production, milk yield and fertility. The etiological agents of helminth infections are mainly Nematodes (roundworms) and Platyhelminths (flatworms). G-quadruplexes (G4) are unusual nucleic acid structures formed by G-rich sequences that can be recognized by specific G4 ligands. Here we used the G4Hunter Web Tool to identify and compare potential G4 sequences (PQS) in the nuclear and mitochondrial genomes of various helminths to identify G4 ligand targets. PQS are nonrandomly distributed in these genomes and often located in the proximity of genes. Unexpectedly, a Nematode, Ascaris lumbricoides, was found to be highly enriched in stable PQS. This species can tolerate high-stability G4 structures, which are not counter selected at all, in stark contrast to most other species. We experimentally confirmed G4 formation for sequences found in four different parasitic helminths. Small molecules able to selectively recognize G4 were found to bind to Schistosoma mansoni G4 motifs. Two of these ligands demonstrated potent activity both against larval and adult stages of this parasite.

Reversal of G-Quadruplexes' Role in Translation Control When Present in the Context of an IRES

G-quadruplexes (GQs) are secondary nucleic acid structures that play regulatory roles in various cellular processes. G-quadruplex-forming sequences present within the 5' UTR of mRNAs can function not only as repressors of translation but also as elements required for optimum function. Based upon previous reports, the majority of the 5' UTR GQ structures inhibit translation, presumably by blocking the ribosome scanning process that is essential for detection of the initiation codon. However, there are certain mRNAs containing GQs that have been identified as positive regulators of translation, as they are needed for translation initiation. While most cellular mRNAs utilize the 5' cap structure to undergo cap-dependent translation initiation, many rely on cap-independent translation under certain conditions in which the cap-dependent initiation mechanism is not viable or slowed down, for example, during development, under stress and in many diseases. Cap-independent translation mainly occurs via Internal Ribosomal Entry Sites (IRESs) that are located in the 5' UTR of mRNAs and are equipped with structural features that can recruit the ribosome or other factors to initiate translation without the need for a 5' cap. In this review, we will focus only on the role of RNA GQs present in the 5' UTR of mRNAs, where they play a critical role in translation initiation, and discuss the potential mechanism of this phenomenon, which is yet to be fully delineated.

Hidden Information Revealed Using the Orthogonal System of Nucleic Acids

In this study, the organization of genetic information in nucleic acids is defined using a novel orthogonal representation. Clearly defined base pairing in DNA allows the linear base chain and sequence to be mathematically transformed into an orthogonal representation where the G–C and A–T pairs are displayed in different planes that are perpendicular to each other. This form of base allocation enables the evaluation of any nucleic acid and predicts the likelihood of a particular region to form non-canonical motifs. The G4Hunter algorithm is currently a popular method of identifying G-quadruplex forming sequences in nucleic acids, and offers promising scores despite its lack of a substantial rational basis. The orthogonal representation described here is an effort to address this incongruity. In addition, the orthogonal display facilitates the search for other sequences that are capable of adopting non-canonical motifs, such as direct and palindromic repeats. The technique can also be used for various RNAs, including any aptamers. This powerful tool based on an orthogonal system offers considerable potential for a wide range of applications.

Pre-miRNA-149 G-quadruplex as a molecular agent to capture nucleolin

One of the most significant challenges in capturing and detecting biomarkers is the choice of an appropriate biomolecular receptor. Recently, RNA G-quadruplexes emerged as plausible receptors due to their ability to recognize with high-affinity proteins. Herein, we have unveiled and characterized the capability of the precursor microRNA 149 to form a G-quadruplex structure and determined the role that some ligands may have in its folding and binding capacity to nucleolin. The G-quadruplex formation was induced by K+ ions and stabilized by ligands, as demonstrated by nuclear magnetic resonance and circular dichroism experiments. Surface plasmon resonance measurements showed a binding affinity of precursor microRNA 149 towards ligands in the micromolar range (10-5-10-6 M) and a strong binding affinity to nucleolin RNA-binding domains 1 and 2 (8.38 × 10-10 M). Even in the presence of the ligand PhenDC3, the binding remains almost identical and in the same order of magnitude (4.46 × 10-10 M). The molecular interactions of the RNA G-quadruplex motif found in precursor miRNA 149 (5'-GGGAGGGAGGGACGGG- 3') and nucleolin RNA-binding domains 1 and 2 were explored by means of molecular docking and molecular dynamics studies. The results showed that RNA G-quadruplex binds to a cavity between domains 1 and 2 of the protein. Then, complex formation was also evaluated through polyacrylamide gel electrophoresis. The results suggest that precursor microRNA 149/ligands and precursor microRNA 149/nucleolin RNA-binding domains 1 and 2 form stable molecular complexes. The in vitro co-localization of precursor microRNA 149 and nucleolin in PC3 cells was demonstrated using confocal microscopy. Finally, a rapid and straightforward microfluidic strategy was employed to check the ability of precursor microRNA 149 to capture nucleolin RNA-binding domains 1 and 2. The results revealed that precursor microRNA 149 can capture nucleolin RNA-binding domains 1 and 2 labeled with Fluorescein 5-isothiocyanate in a concentration-dependent manner, but PhenDC3 complexation seems to decrease the ability of precursor microRNA 149 to capture the protein. Overall, our results proved the formation of the G-quadruplex structure in the precursor microRNA 149 and the ability to recognize and detect nucleolin. This proof-of-concept study could open up a new framework for developing new strategies to design improved molecular receptors for capture and detection of nucleolin in complex biological samples.

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.

Conservation and over-representation of G-quadruplex sequences in regulatory regions of mitochondrial DNA across distinct taxonomic sub-groups

G-quadruplexes have important regulatory roles in the nuclear genome but their distribution and potential roles in mitochondrial DNA (mtDNA) are poorly understood. We analysed 11883 mtDNA sequences from 18 taxonomic sub-groups and identified their frequency and location within mtDNA. Large differences in both the frequency and number of putative quadruplex-forming sequences (PQS) were observed amongst all the organisms and PQS frequency was negatively correlated with an increase in evolutionary age. PQS were over-represented in the 3'UTRs, D-loops, replication origins, and stem loops, indicating regulatory roles for quadruplexes in mtDNA. Variations of the G-quadruplex-forming sequence in the conserved sequence block II (CSBII) region of the human D-loop were conserved amongst other mammals, amphibians, birds, reptiles, and fishes. This D-loop PQS was conserved in the duplicated control regions of some birds and reptiles, indicating its importance to mitochondrial function. The guanine tracts in these PQS also displayed significant length heterogeneity and the length of these guanine tracts were generally longest in bird mtDNA. This information provides further insights into how G4s may contribute to the regulation and function of mtDNA and acts as a database of information for future studies investigating mitochondrial G4s in organisms other than humans.

R-Loop Tracker: Web Access-Based Tool for R-Loop Detection and Analysis in Genomic DNA Sequences

R-loops are common non-B nucleic acid structures formed by a three-stranded nucleic acid composed of an RNA–DNA hybrid and a displaced single-stranded DNA (ssDNA) loop. Because the aberrant R-loop formation leads to increased mutagenesis, hyper-recombination, rearrangements, and transcription-replication collisions, it is regarded as important in human diseases. Therefore, its prevalence and distribution in genomes are studied intensively. However, in silico tools for R-loop prediction are limited, and therefore, we have developed the R-loop tracker tool, which was implemented as a part of the DNA Analyser web server. This new tool is focused upon (1) prediction of R-loops in genomic DNA without length and sequence limitations; (2) integration of R-loop tracker results with other tools for nucleic acids analyses, including Genome Browser; (3) internal cross-evaluation of in silico results with experimental data, where available; (4) easy export and correlation analyses with other genome features and markers; and (5) enhanced visualization outputs. Our new R-loop tracker tool is freely accessible on the web pages of DNA Analyser tools, and its implementation on the web-based server allows effective analyses not only for DNA segments but also for full chromosomes and genomes.

G3BP1 binds to guanine quadruplexes in mRNAs to modulate their stabilities

RNA guanine quadruplexes (rG4) assume important roles in post-transcriptional regulations of gene expression, which are often modulated by rG4-binding proteins. Hence, understanding the biological functions of rG4s requires the identification and functional characterizations of rG4-recognition proteins. By employing a bioinformatic approach based on the analysis of overlap between peaks obtained from rG4-seq analysis and those detected in >230 eCLIP-seq datasets for RNA-binding proteins generated from the ENCODE project, we identified a large number of candidate rG4-binding proteins. We showed that one of these proteins, G3BP1, is able to bind directly to rG4 structures with high affinity and selectivity, where the binding entails its C-terminal RGG domain and is further enhanced by its RRM domain. Additionally, our seCLIP-Seq data revealed that pyridostatin, a small-molecule rG4 ligand, could displace G3BP1 from mRNA in cells, with the most pronounced effects being observed for the 3′-untranslated regions (3′-UTR) of mRNAs. Moreover, luciferase reporter assay results showed that G3BP1 positively regulates mRNA stability through its binding with rG4 structures. Together, we identified a number of candidate rG4-binding proteins and validated that G3BP1 can bind directly with rG4 structures and regulate the stabilities of mRNAs.

New telomere to telomere assembly of human chromosome 8 reveals a previous underestimation of G-quadruplex forming sequences and inverted repeats

Taking advantage of evolving and improving sequencing methods, human chromosome 8 is now available as a gapless, end-to-end assembly. Thanks to advances in long-read sequencing technologies, its centromere, telomeres, duplicated gene families and repeat-rich regions are now fully sequenced. We were interested to assess if the new assembly altered our understanding of the potential impact of non-B DNA structures within this completed chromosome sequence. It has been shown that non-B secondary structures, such as G-quadruplexes, hairpins and cruciforms, have important regulatory functions and potential as targeted therapeutics. Therefore, we analysed the presence of putative G-quadruplex forming sequences and inverted repeats in the current human reference genome (GRCh38) and in the new end-to-end assembly of chromosome 8. The comparison revealed that the new assembly contains significantly more inverted repeats and G-quadruplex forming sequences compared to the current reference sequence. This observation can be explained by improved accuracy of the new sequencing methods, particularly in regions that contain extensive repeats of bases, as is preferred by many non-B DNA structures. These results show a significant underestimation of the prevalence of non-B DNA secondary structure in previous assembly versions of the human genome and point to their importance being not fully appreciated. We anticipate that similar observations will occur as the improved sequencing technologies fill in gaps across the genomes of humans and other organisms.

Analysis of nucleoid-associated protein-binding regions reveals DNA structural features influencing genome organization in Mycobacterium tuberculosis

Nucleoid-associated proteins (NAPs) maintain bacterial nucleoid configuration through their architectural properties of DNA bending, wrapping, and bridging. However, the contribution of DNA structural alterations to DNA-NAP recognition at the genomic scale remains unresolved. Present work dissects the DNA sequence, shape and altered structural preferences at a genomic scale for six NAPs in Mycobacterium tuberculosis. Results suggest narrower minor groove width (MGW) and higher DNA rigidity are marked for the binding sites of EspR and Lsr2, while mIHF, MtHU and NapM have heterogeneous DNA structural predilections. In contrast, WhiB4-DNA-binding sites were characterized by wider MGW, highly deformable and less curved DNA. This work provides systematic insight into NAP-mediated genome organization as a function of DNA structural features.

G-quadruplex DNA: a novel target for drug design

G-quadruplex (G4) DNA is a type of quadruple helix structure formed by a continuous guanine-rich DNA sequence. Emerging evidence in recent years authenticated that G4 DNA structures exist both in cell-free and cellular systems, and function in different diseases, especially in various cancers, aging, neurological diseases, and have been considered novel promising targets for drug design. In this review, we summarize the detection method and the structure of G4, highlighting some non-canonical G4 DNA structures, such as G4 with a bulge, a vacancy, or a hairpin. Subsequently, the functions of G4 DNA in physiological processes are discussed, especially their regulation of DNA replication, transcription of disease-related genes (c-MYC, BCL-2, KRAS, c-KIT et al.), telomere maintenance, and epigenetic regulation. Typical G4 ligands that target promoters and telomeres for drug design are also reviewed, including ellipticine derivatives, quinoxaline analogs, telomestatin analogs, berberine derivatives, and CX-5461, which is currently in advanced phase I/II clinical trials for patients with hematologic cancer and BRCA1/2-deficient tumors. Furthermore, since the long-term stable existence of G4 DNA structures could result in genomic instability, we summarized the G4 unfolding mechanisms emerged recently by multiple G4-specific DNA helicases, such as Pif1, RecQ family helicases, FANCJ, and DHX36. This review aims to present a general overview of the field of G-quadruplex DNA that has progressed in recent years and provides potential strategies for drug design and disease treatment.

G-quadruplexes in genomes of viruses infecting eukaryotes or prokaryotes are under different selection pressures from hosts

G-quadruplexes in viral genomes can be applied as the targets of antiviral therapies, which has attracted wide interest. However, it is still not clear whether the pervasive number of such elements in the viral world is the result of natural selection for functionality. In this study, we identified putative quadruplex-forming sequences (PQSs) across the known viral genomes and analyzed the abundance, structural stability, and conservation of viral PQSs. A Viral Putative G-quadruplex Database (ViPGD,http://jsjds.hzau.edu.cn/MBPC/ViPGD/index.php/home/index) was constructed to collect the details of each viral PQS, which provides guidance for selecting the desirable PQS. The PQS with two putative G-tetrads (G₂-PQS) was significantly enriched in both eukaryotic viruses and prokaryotic viruses, while the PQSs with three putative G-tetrads (G₃-PQS) were only enriched in eukaryotic viruses and depleted in prokaryotic viruses. The structural stability of PQSs in prokaryotic viruses was significantly lower than that in eukaryotic viruses. Conservation analysis showed that the G₂-PQS, instead of G₃-PQS, was highly conserved within the genus. This suggested that the G₂-quadruplex might play an important role in viral biology, and the difference in the occurrence of G-quadruplex between eukaryotic viruses and prokaryotic viruses may result from the different selection pressures from hosts.

Novel G-quadruplex prone sequences emerge in the complete assembly of the human X chromosome

G-quadruplexes are non-B secondary structures with regulatory functions and therapeutic potential. Improvements in sequencing methods recently allowed the completion of the first human chromosome which is now available as a gapless, end-to-end assembly, with the previously remaining spaces filled and newly identified regions added. We compared the presence of G-quadruplex forming sequences in the current human reference genome (GRCh38) and in the new end-to-end assembly of the X chromosome constructed by high-coverage ultra-long-read nanopore sequencing. This comparison revealed that, even though the corrected length of the chromosome X assembly is surprisingly 1.14% shorter than expected, the number of G-quadruplex forming sequences found in this gapless chromosome is significantly higher, with 493 new motifs having G4Hunter scores above 1.4 and 23 new sequences with G4Hunter scores above 3.5. This observation reflects an improved precision of the new sequencing approaches and points to an underestimation of G-quadruplex propensity in the previous, widely used version of the human genome assembly, especially for motifs with a high G4Hunter score, expected to be very stable. These G-quadruplex forming sequences probably remained undiscovered in earlier genome datasets due to previously unsolved G-rich and repetitive genomic regions. These observations allow a precise targeting of these important regulatory regions.

G-Quadruplexes and Their Ligands: Biophysical Methods to Unravel G-Quadruplex/Ligand Interactions

Progress in the design of G-quadruplex (G4) binding ligands relies on the availability of approaches that assess the binding mode and nature of the interactions between G4 forming sequences and their putative ligands. The experimental approaches used to characterize G4/ligand interactions can be categorized into structure-based methods (circular dichroism (CD), nuclear magnetic resonance (NMR) spectroscopy and X-ray crystallography), affinity and apparent affinity-based methods (surface plasmon resonance (SPR), isothermal titration calorimetry (ITC) and mass spectrometry (MS)), and high-throughput methods (fluorescence resonance energy transfer (FRET)-melting, G4-fluorescent intercalator displacement assay (G4-FID), affinity chromatography and microarrays. Each method has unique advantages and drawbacks, which makes it essential to select the ideal strategies for the biological question being addressed. The structural- and affinity and apparent affinity-based methods are in several cases complex and/or time-consuming and can be combined with fast and cheap high-throughput approaches to improve the design and development of new potential G4 ligands. In recent years, the joint use of these techniques permitted the discovery of a huge number of G4 ligands investigated for diagnostic and therapeutic purposes. Overall, this review article highlights in detail the most commonly used approaches to characterize the G4/ligand interactions, as well as the applications and types of information that can be obtained from the use of each technique.

A CpG island promoter drives the CXXC5 gene expression

CXXC5 is a member of the zinc-finger CXXC family that binds to unmethylated CpG dinucleotides. CXXC5 modulates gene expressions resulting in diverse cellular events mediated by distinct signaling pathways. However, the mechanism responsible for CXXC5 expression remains largely unknown. We found here that of the 14 annotated CXXC5 transcripts with distinct 5' untranslated regions encoding the same protein, transcript variant 2 with the highest expression level among variants represents the main transcript in cell models. The DNA segment in and at the immediate 5'-sequences of the first exon of variant 2 contains a core promoter within which multiple transcription start sites are present. Residing in a region with high G-C nucleotide content and CpG repeats, the core promoter is unmethylated, deficient in nucleosomes, and associated with active RNA polymerase-II. These findings suggest that a CpG island promoter drives CXXC5 expression. Promoter pull-down revealed the association of various transcription factors (TFs) and transcription co-regulatory proteins, as well as proteins involved in histone/chromatin, DNA, and RNA processing with the core promoter. Of the TFs, we verified that ELF1 and MAZ contribute to CXXC5 expression. Moreover, the first exon of variant 2 may contain a G-quadruplex forming region that could modulate CXXC5 expression.