G-Quadruplex (G4) Motifs in the Maize (Zea mays L.) Genome Are Enriched at Specific Locations in Thousands of Genes Coupled to Energy Status, Hypoxia, Low Sugar, and Nutrient Deprivation

The Characterization of G-Quadruplexes in Tobacco Genome and Their Function under Abiotic Stress

Tobacco is an ideal model plant in scientific research. G-quadruplex is a guanine-rich DNA structure, which regulates transcription and translation. In this study, the prevalence and potential function of G-quadruplexes in tobacco were systematically analyzed. In tobacco genomes, there were 2,924,271,002 G-quadruplexes in the nuclear genome, 430,597 in the mitochondrial genome, and 155,943 in the chloroplast genome. The density of the G-quadruplex in the organelle genome was higher than that in the nuclear genome. G-quadruplexes were abundant in the transcription regulatory region of the genome, and a difference in G-quadruplex density in two DNA strands was also observed. The promoter of 60.4% genes contained at least one G-quadruplex. Compared with up-regulated differentially expressed genes (DEGs), the G-quadruplex density in down-regulated DEGs was generally higher under drought stress and salt stress. The G-quadruplex formed by simple sequence repeat (SSR) and its flanking sequence in the promoter region of the NtBBX (Nitab4.5_0002943g0010) gene might enhance the drought tolerance of tobacco. This study lays a solid foundation for further research on G-quadruplex function in tobacco and other plants.

The Epigenomic Features and Potential Functions of PEG- and PDS-Favorable DNA G-Quadruplexes in Rice

A G-quadruplex (G4) is a typical non-B DNA structure and involved in various DNA-templated events in eukaryotic genomes. PEG and PDS chemicals have been widely applied for promoting the folding of in vivo or in vitro G4s. However, how PEG and PDS preferentially affect a subset of G4 formation genome-wide is still largely unknown. We here conducted a BG4-based IP-seq in vitro under K++PEG or K++PDS conditions in the rice genome. We found that PEG-favored IP-G4s+ have distinct sequence features, distinct genomic distributions and distinct associations with TEGs, non-TEGs and subtypes of TEs compared to PDS-favored ones. Strikingly, PEG-specific IP-G4s+ are associated with euchromatin with less enrichment levels of DNA methylation but with more enriched active histone marks, while PDS-specific IP-G4s+ are associated with heterochromatin with higher enrichment levels of DNA methylation and repressive marks. Moreover, we found that genes with PEG-specific IP-G4s+ are more expressed than those with PDS-specific IP-G4s+, suggesting that PEG/PDS-specific IP-G4s+ alone or coordinating with epigenetic marks are involved in the regulation of the differential expression of related genes, therefore functioning in distinct biological processes. Thus, our study provides new insights into differential impacts of PEG and PDS on G4 formation, thereby advancing our understanding of G4 biology.

A Molecular Orchestration of Plant Translation under Abiotic Stress

The complexities of translational strategies make this stage of implementing genetic information one of the most challenging to comprehend and, simultaneously, perhaps the most engaging. It is evident that this diverse range of strategies results not only from a long evolutionary history, but is also of paramount importance for refining gene expression and metabolic modulation. This notion is particularly accurate for organisms that predominantly exhibit biochemical and physiological reactions with a lack of behavioural ones. Plants are a group of organisms that exhibit such features. Addressing unfavourable environmental conditions plays a pivotal role in plant physiology. This is particularly evident with the changing conditions of global warming and the irrevocable loss or depletion of natural ecosystems. In conceptual terms, the plant response to abiotic stress comprises a set of elaborate and intricate strategies. This is influenced by a range of abiotic factors that cause stressful conditions, and molecular genetic mechanisms that fine-tune metabolic pathways allowing the plant organism to overcome non-standard and non-optimal conditions. This review aims to focus on the current state of the art in the field of translational regulation in plants under abiotic stress conditions. Different regulatory elements and patterns are being assessed chronologically. We deem it important to focus on significant high-performance techniques for studying the genetic information dynamics during the translation phase.

Weighted gene co-expression network analysis of nitrogen (N)-responsive genes and the putative role of G-quadruplexes in N use efficiency (NUE) in rice

Rice is an important target to improve crop nitrogen (N) use efficiency (NUE), and the identification and shortlisting of the candidate genes are still in progress. We analyzed data from 16 published N-responsive transcriptomes/microarrays to identify, eight datasets that contained the maximum number of 3020 common genes, referred to as N-responsive genes. These include different classes of transcription factors, transporters, miRNA targets, kinases and events of post-translational modifications. A Weighted gene co-expression network analysis (WGCNA) with all the 3020 N-responsive genes revealed 15 co-expression modules and their annotated biological roles. Protein-protein interaction network analysis of the main module revealed the hub genes and their functional annotation revealed their involvement in the ubiquitin process. Further, the occurrences of G-quadruplex sequences were examined, which are known to play important roles in epigenetic regulation but are hitherto unknown in N-response/NUE. Out of the 3020 N-responsive genes studied, 2298 contained G-quadruplex sequences. We compared these N-responsive genes containing G-quadruplex sequences with the 3601 genes we previously identified as NUE-related (for being both N-responsive and yield-associated). This analysis revealed 389 (17%) NUE-related genes containing G-quadruplex sequences. These genes may be involved in the epigenetic regulation of NUE, while the rest of the 83% (1811) genes may regulate NUE through genetic mechanisms and/or other epigenetic means besides G-quadruplexes. A few potentially important genes/processes identified as associated with NUE were experimentally validated in a pair of rice genotypes contrasting for NUE. The results from the WGCNA and G4 sequence analysis of N-responsive genes helped identify and shortlist six genes as candidates to improve NUE. Further, the hitherto unavailable segregation of genetic and epigenetic gene targets could aid in informed interventions through genetic and epigenetic means of crop improvement.

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.

Epigenomic Features and Potential Functions of K+ and Na+ Favorable DNA G-Quadruplexes in Rice

DNA G-quadruplexes (G4s) are non-canonical four-stranded DNA structures involved in various biological processes in eukaryotes. Molecularly crowded solutions and monovalent cations have been reported to stabilize in vitro and in vivo G4 formation. However, how K⁺ and Na⁺ affect G4 formation genome-wide is still unclear in plants. Here, we conducted BG4-DNA-IP-seq, DNA immunoprecipitation with anti-BG4 antibody coupled with sequencing, under K⁺ and Na⁺ + PEG conditions in vitro. We found that K⁺-specific IP-G4s had a longer peak size, more GC and PQS content, and distinct AT and GC skews compared to Na⁺-specific IP-G4s. Moreover, K⁺- and Na⁺-specific IP-G4s exhibited differential subgenomic enrichment and distinct putative functional motifs for the binding of certain trans-factors. More importantly, we found that K⁺-specific IP-G4s were more associated with active marks, such as active histone marks, and low DNA methylation levels, as compared to Na⁺-specific IP-G4s; thus, K⁺-specific IP-G4s in combination with active chromatin features facilitate the expression of overlapping genes. In addition, K⁺- and Na⁺-specific IP-G4 overlapping genes exhibited differential GO (gene ontology) terms, suggesting they may have distinct biological relevance in rice. Thus, our study, for the first time, explores the effects of K⁺ and Na⁺ on global G4 formation in vitro, thereby providing valuable resources for functional G4 studies in rice. It will provide certain G4 loci for the biotechnological engineering of rice in the future.

The potential interplay between G-quadruplex and p53: their roles in regulation of ferroptosis in cancer

Ferroptosis is a novel form of regulated cell death trigged by various biological processes, and p53 is involved in different ferroptosis regulations and functions as a crucial regulator. Both DNA and RNA can fold into G-quadruplex in GC-rich regions and increasing shreds of evidence demonstrate that G-quadruplexes have been associated with some important cellular events. Investigation of G-quadruplexes is thus vital to revealing their biological functions. Specific G-quadruplexes are investigated to discover new effective anticancer drugs. Multiple modulations have been discovered between the secondary structure G-quadruplex and p53, probably further influencing the ferroptosis in cancer. G-quadruplex binds to ferric iron-related structures directly and may affect the p53 pathways as well as ferroptosis in cancer. In addition, G-quadruplex also interacts with p53 indirectly, including iron-sulfur cluster metabolism, telomere homeostasis, lipid peroxidation, and glycolysis. In this review, we summarized the latent interplay between G-quadruplex and p53 which focused mainly on ferroptosis in cancer to provide the potential understanding and encourage future studies.

G4Boost: a machine learning-based tool for quadruplex identification and stability prediction

Background

G-quadruplexes (G4s), formed within guanine-rich nucleic acids, are secondary structures involved in important biological processes. Although every G4 motif has the potential to form a stable G4 structure, not every G4 motif would, and accurate energy-based methods are needed to assess their structural stability. Here, we present a decision tree-based prediction tool, G4Boost, to identify G4 motifs and predict their secondary structure folding probability and thermodynamic stability based on their sequences, nucleotide compositions, and estimated structural topologies.

Results

G4Boost predicted the quadruplex folding state with an accuracy greater then 93% and an F1-score of 0.96, and the folding energy with an RMSE of 4.28 and R² of 0.95 only by the means of sequence intrinsic feature. G4Boost was successfully applied and validated to predict the stability of experimentally-determined G4 structures, including for plants and humans.

Conclusion

G4Boost outperformed the three machine-learning based prediction tools, DeepG4, Quadron, and G4RNA Screener, in terms of both accuracy and F1-score, and can be highly useful for G4 prediction to understand gene regulation across species including plants and humans.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04782-z.

Epigenomic features of DNA G-quadruplexes and their roles in regulating rice gene transcription

A DNA G-quadruplex (G4) is a non-canonical four-stranded nucleic acid structure involved in many biological processes in mammals. The current knowledge on plant DNA G4s, however, is limited; whether and how DNA G4s impact gene expression in plants is still largely unknown. Here, we applied a protocol referred to as BG4-DNA-IP-seq followed by a comprehensive characterization of DNA G4s in rice (Oryza sativa L.); we next integrated dG4s (experimentally detectable G4s) with existing omics data and found that dG4s exhibited differential DNA methylation between transposable element (TE) and non-TE genes. dG4 regions displayed genic-dependent enrichment of epigenomic signatures; finally, we showed that these sites displayed a positive association with expression of DNA G4-containing genes when located at promoters, and a negative association when located in the gene body, suggesting localization-dependent promotional/repressive roles of DNA G4s in regulating gene transcription. This study reveals interrelations between DNA G4s and epigenomic signatures, as well as implicates DNA G4s in modulating gene transcription in rice. Our study provides valuable resources for the functional characterization or bioengineering of some of key DNA G4s in rice.

The Dynamics of NO3- and NH4+ Uptake in Duckweed Are Coordinated with the Expression of Major Nitrogen Assimilation Genes

Duckweed plants play important roles in aquatic ecosystems worldwide. They rapidly accumulate biomass and have potential uses in bioremediation of water polluted by fertilizer runoff or other chemicals. Here we studied the assimilation of two major sources of inorganic nitrogen, nitrate (NO3- ) and ammonium (NH4+), in six duckweed species: Spirodela polyrhiza, Landoltia punctata, Lemna aequinoctialis, Lemna turionifera, Lemna minor, and Wolffia globosa. All six duckweed species preferred NH4+ over NO3- and started using NO3- only when NH4+ was depleted. Using the available genome sequence, we analyzed the molecular structure and expression of eight key nitrogen assimilation genes in S. polyrhiza. The expression of genes encoding nitrate reductase and nitrite reductase increased about 10-fold when NO3- was supplied and decreased when NH4+ was supplied. NO3- and NH4+ induced the glutamine synthetase (GS) genes GS1;2 and the GS2 by 2- to 5-fold, respectively, but repressed GS1;1 and GS1;3. NH4+ and NO3- upregulated the genes encoding ferredoxin- and NADH-dependent glutamate synthases (Fd-GOGAT and NADH-GOGAT). A survey of nitrogen assimilation gene promoters suggested complex regulation, with major roles for NRE-like and GAATC/GATTC cis-elements, TATA-based enhancers, GA/CTn repeats, and G-quadruplex structures. These results will inform efforts to improve bioremediation and nitrogen use efficiency.

Searching for G-Quadruplex-Binding Proteins in Plants: New Insight into Possible G-Quadruplex Regulation

G-quadruplexes are four-stranded nucleic acid structures occurring in the genomes of all living organisms and viruses. It is increasingly evident that these structures play important molecular roles; generally, by modulating gene expression and overall genome integrity. For a long period, G-quadruplexes have been studied specifically in the context of human promoters, telomeres, and associated diseases (cancers, neurological disorders). Several of the proteins for binding G-quadruplexes are known, providing promising targets for influencing G-quadruplex-related processes in organisms. Nonetheless, in plants, only a small number of G-quadruplex binding proteins have been described to date. Thus, we aimed to bioinformatically inspect the available protein sequences to find the best protein candidates with the potential to bind G-quadruplexes. Two similar glycine and arginine-rich G-quadruplex-binding motifs were described in humans. The first is the so-called “RGG motif”-RRGDGRRRGGGGRGQGGRGRGGGFKG, and the second (which has been recently described) is known as the “NIQI motif”-RGRGRGRGGGSGGSGGRGRG. Using this general knowledge, we searched for plant proteins containing the above mentioned motifs, using two independent approaches (BLASTp and FIMO scanning), and revealed many proteins containing the G4-binding motif(s). Our research also revealed the core proteins involved in G4 folding and resolving in green plants, algae, and the key plant model organism, Arabidopsis thaliana. The discovered protein candidates were annotated using STRINGdb and sorted by their molecular and physiological roles in simple schemes. Our results point to the significant role of G4-binding proteins in the regulation of gene expression in plants.

G-Quadruplex Structures in Bacteria: Biological Relevance and Potential as an Antimicrobial Target

DNA strands consisting of multiple runs of guanines can adopt a noncanonical, four-stranded DNA secondary structure known as G-quadruplex or G4 DNA. G4 DNA is thought to play an important role in transcriptional and translational regulation of genes, DNA replication, genome stability, and oncogene expression in eukaryotic genomes. In other organisms, including several bacterial pathogens and some plant species, the biological roles of G4 DNA and G4 RNA are starting to be explored. Recent investigations showed that G4 DNA and G4 RNA are generally conserved across plant species. In silico analyses of several bacterial genomes identified putative guanine-rich, G4 DNA-forming sequences in promoter regions. The sequences were particularly abundant in certain gene classes, suggesting that these highly diverse structures can be employed to regulate the expression of genes involved in secondary metabolite synthesis and signal transduction. Furthermore, in the pathogen Mycobacterium tuberculosis, the distribution of G4 motifs and their potential role in the regulation of gene transcription advocate for the use of G4 ligands to develop novel antitubercular therapies. In this review, we discuss the various roles of G4 structures in bacterial DNA and the application of G4 DNA as inhibitors or therapeutic agents to address bacterial pathogens.

Quadruplex-Forming Motif Inserted into 3'UTR of Ty1his3-AI Retrotransposon Inhibits Retrotransposition in Yeast

Guanine quadruplexes (G4s) serve as regulators of replication, recombination and gene expression. G4 motifs have been recently identified in LTR retrotransposons, but their role in the retrotransposon life-cycle is yet to be understood. Therefore, we inserted G4s into the 3'UTR of Ty1his3-AI retrotransposon and measured the frequency of retrotransposition in yeast strains BY4741, Y00509 (without Pif1 helicase) and with G4-stabilization by N-methyl mesoporphyrin IX (NMM) treatment. We evaluated the impact of G4s on mRNA levels by RT-qPCR and products of reverse transcription by Southern blot analysis. We found that the presence of G4 inhibited Ty1his3-AI retrotransposition. The effect was stronger when G4s were on a transcription template strand which leads to reverse transcription interruption. Both NMM and Pif1p deficiency reduced the retrotransposition irrespective of the presence of a G4 motif in the Ty1his3-AI element. Quantity of mRNA and products of reverse transcription did not fully explain the impact of G4s on Ty1his3-AI retrotransposition indicating that G4s probably affect some other steps of the retrotransposon life-cycle (e.g., translation, VLP formation, integration). Our results suggest that G4 DNA conformation can tune the activity of mobile genetic elements that in turn contribute to shaping the eukaryotic genomes.

Genome-wide discovery of G-quadruplexes in barley

G-quadruplexes (G4s) are four-stranded nucleic acid structures with closely spaced guanine bases forming square planar G-quartets. Aberrant formation of G4 structures has been associated with genomic instability. However, most plant species are lacking comprehensive studies of G4 motifs. In this study, genome-wide identification of G4 motifs in barley was performed, followed by a comparison of genomic distribution and molecular functions to other monocot species, such as wheat, maize, and rice. Similar to the reports on human and some plants like wheat, G4 motifs peaked around the 5′ untranslated region (5′ UTR), the first coding domain sequence, and the first intron start sites on antisense strands. Our comparative analyses in human, Arabidopsis, maize, rice, and sorghum demonstrated that the peak points could be erroneously merged into a single peak when large window sizes are used. We also showed that the G4 distributions around genic regions are relatively similar in the species studied, except in the case of Arabidopsis. G4 containing genes in monocots showed conserved molecular functions for transcription initiation and hydrolase activity. Additionally, we provided examples of imperfect G4 motifs.

Sugar modulation of anaerobic-response networks in maize root tips

Sugar supply is a key component of hypoxia tolerance and acclimation in plants. However, a striking gap remains in our understanding of mechanisms governing sugar impacts on low-oxygen responses. Here, we used a maize (Zea mays) root-tip system for precise control of sugar and oxygen levels. We compared responses to oxygen (21 and 0.2%) in the presence of abundant versus limited glucose supplies (2.0 and 0.2%). Low-oxygen reconfigured the transcriptome with glucose deprivation enhancing the speed and magnitude of gene induction for core anaerobic proteins (ANPs). Sugar supply also altered profiles of hypoxia-responsive genes carrying G4 motifs (sources of regulatory quadruplex structures), revealing a fast, sugar-independent class followed more slowly by feast-or-famine-regulated G4 genes. Metabolite analysis showed that endogenous sugar levels were maintained by exogenous glucose under aerobic conditions and demonstrated a prominent capacity for sucrose re-synthesis that was undetectable under hypoxia. Glucose abundance had distinctive impacts on co-expression networks associated with ANPs, altering network partners and aiding persistence of interacting networks under prolonged hypoxia. Among the ANP networks, two highly interconnected clusters of genes formed around Pyruvate decarboxylase 3 and Glyceraldehyde-3-phosphate dehydrogenase 4. Genes in these clusters shared a small set of cis-regulatory elements, two of which typified glucose induction. Collective results demonstrate specific, previously unrecognized roles of sugars in low-oxygen responses, extending from accelerated onset of initial adaptive phases by starvation stress to maintenance and modulation of co-expression relationships by carbohydrate availability.

Genome-wide analysis of DNA G-quadruplex motifs across 37 species provides insights into G4 evolution

G-quadruplex (G4) structures have been predicted in the genomes of many organisms and proven to play regulatory roles in diverse cellular activities. However, there is little information on the evolutionary history and distribution characteristics of G4s. Here, whole-genome characteristics of potential G4s were studied in 37 evolutionarily representative species. During evolution, the number, length, and density of G4s generally increased. Immunofluorescence in seven species confirmed G4s' presence and evolutionary pattern. G4s tended to cluster in chromosomes and were enriched in genetic regions. Short-loop G4s were conserved in most species, while loop-length diversity also existed, especially in mammals. The proportion of G4-bearing genes and orthologue genes, which appeared to be increasingly enriched in transcription factors, gradually increased. The antagonistic relationship between G4s and DNA methylation sites was detected. These findings imply that organisms may have evolutionarily developed G4 into a novel reversible and elaborate transcriptional regulatory mechanism benefiting multiple physiological activities of higher organisms.

A G(enomic)P(ositioning)S(ystem) for Plant RNAPII Transcription

Post-translational modifications (PTMs) of histone residues shape the landscape of gene expression by modulating the dynamic process of RNA polymerase II (RNAPII) transcription. The contribution of particular histone modifications to the definition of distinct RNAPII transcription stages remains poorly characterized in plants. Chromatin immunoprecipitation combined with next-generation sequencing (ChIP-seq) resolves the genomic distribution of histone modifications. Here, we review histone PTM ChIP-seq data in Arabidopsis thaliana and find support for a Genomic Positioning System (GPS) that guides RNAPII transcription. We review the roles of histone PTM 'readers', 'writers', and 'erasers', with a focus on the regulation of gene expression and biological functions in plants. The distinct functions of RNAPII transcription during the plant transcription cycle may rely, in part, on the characteristic histone PTM profiles that distinguish transcription stages.

Genome-Wide Discovery of G-Quadruplexes in Wheat: Distribution and Putative Functional Roles

G-quadruplexes are nucleic acid secondary structures formed by a stack of square planar G-quartets. G-quadruplexes were implicated in many biological functions including telomere maintenance, replication, transcription, and translation, in many species including humans and plants. For wheat, however, though it is one of the world's most important staple food, no G-quadruplex studies have been reported to date. Here, we computationally identify putative G4 structures (G4s) in wheat genome for the first time and compare its distribution across the genome against five other genomes (human, maize, Arabidopsis, rice, and sorghum). We identified close to 1 million G4 motifs with a density of 76 G4s/Mb across the whole genome and 93 G4s/Mb over genic regions. Remarkably, G4s were enriched around three regions, two located on the antisense and one on the sense strand at the following positions: 1) the transcription start site (TSS) (antisense), 2) the first coding domain sequence (CDS) (antisense), and 3) the start codon (sense). Functional enrichment analysis revealed that the gene models containing G4 motifs within these peaks were associated with specific gene ontology (GO) terms, such as developmental process, localization, and cellular component organization or biogenesis. We investigated genes encoding MADS-box transcription factors and showed examples of G4 motifs within critical regulatory regions in the VRN-1 genes in wheat. Furthermore, comparison with other plants showed that monocots share a similar distribution of G4s, but Arabidopsis shows a unique G4 distribution. Our study shows for the first time the prevalence and possible functional roles of G4s in wheat.

The Importance of Phosphates for DNA G-Quadruplex Formation: Evaluation of Zwitterionic G-Rich Oligodeoxynucleotides

A quaternary ammonium butylsulfonyl phosphoramidate group (N+) was designed to replace all the phosphates in a G-rich oligodeoxynucleotide d(TG₄ T), resulting in a formally charge-neutral zwitterionic N+TG₄ T sequence. We evaluated the effects of N+phosphate modifications on the structural, thermodynamic and kinetic properties of the parallel G-quadruplexes (G4) formed by TG₄ T and compared them to the properties of the recently published phosphoryl guanidine d(TG₄ T) (PG-TG₄ T). Using size-exclusion chromatography, we established that, unlike PG-TG₄ T, which exists as a mixture of complexes of different molecularity in solution, N+TG₄ T forms an individual tetramolecular complex. In contrast to PG modifications that destabilized G4s, the presence of N+ modifications increased thermal stability relative to unmodified [d(TG₄ T)]₄ . The initial stage of assembly of N+TG₄ T proceeded faster in the presence of Na⁺ than K⁺ ions and, similarly to PG-TG₄ T, was independent of the salt concentration. However, after complex formation exceeded 75 %, N+TG₄ T in solution with Na⁺ showed slower association than with K⁺ . N+TG₄ T could also form G4s in solution with Li⁺ ions at a very low strand concentration (10 μM); something that has never been reported for the native d(TG₄ T). Charge-neutral PG-G4s can invade preformed native G4s, whereas no invasion was observed between N+and native G4s, possibly due to the increased thermal stability of [N+TG₄ T]₄ . The N+ modification makes d(TG₄ T) fully resistant to enzymatic digestion, which could be useful for intracellular application of N+-modified DNA or RNA.

Tight DNA-protein complexes isolated from barley seedlings are rich in potential guanine quadruplex sequences

Background

The concept of chromatin domains attached to the nuclear matrix is being revisited, with nucleus described as a set of topologically associating domains. The significance of the tightly bound to DNA proteins (TBP), a protein group that remains attached to DNA after its deproteinization should be also revisited, as the existence of these interactions is in good agreement with the concept of the topologically associating domain. The work aimed to characterize the DNA component of TBP isolated from barley seedlings.

Methods

The tight DNA-protein complexes from the first leaves, coleoptiles, and roots of barley seedlings were isolated by purification with chromatography on nitrocellulose or exhaustive digestion of DNA with DNase I. Cloning and transformation were performed using pMOSBBlue Blunt Ended Cloning Kit. Inserts were amplified by PCR, and sequencing was performed on the MegaBace 1000 Sequencing System. The BLAST search was performed using sequence databases at NCBI, CR-EST, and TREP and Ensembl Plants databases. Comparison to MAR/SAR sequences was performed using http://smartdb.bioinf.med.uni-goettingen.de/cgi-bin/SMARtDB/smar.cgi database. The prediction of G quadruplexes (GQ) was performed with the aid of R-studio library pqsfinder. CD spectra were recorded on a Chirascan CS/3D spectrometer.

Results

Although the barley genome is AT-rich (43% of GC pairs), most DNA fragments associated with TBP were GC-rich (up to 70% in some fractions). Both fractionation procedures yielded a high proportion of CT-motif sequences presented predominantly by the 16-bp CC(TCTCCC)₂ TC fragment present in clones derived from the TBP-bound DNA and absent in free DNA. BLAST analysis revealed alignment with different barley repeats. Some clones, however, aligned with both nuclear and chloroplast structural genes. Alignments with MAR/SAR motifs were very few. The analysis produced by the pqsfinder program revealed numerous potential quadruplex-forming sites in the TBP-bound sequences. A set of oligonucleotides containing sites of possible GQs were designed and ordered. Three of them represented the minus strand of the CT-repeat. Two were derived from sequences of two clones of nitrocellulose retained fraction from leaves and contained GC-rich motifs different from the CT motif. Circular dichroism spectroscopy revealed profound changes in spectra when oligonucleotides were incubated with 100 mM KCl. There was either an increase of positive band in the area of 260 nm or the formation of a positive band at 290 nm. In the former case, changes are typical for parallel G-quadruplexes and, in the latter, 3 + 1 structures.

Discussion

The G-quadruplexes anchor proteins are probably involved in the maintenance of the topologically associated domain structure.

Fertility restoration of maize CMS-C altered by a single amino acid substitution within the Rf4 bHLH transcription factor

Type C cytoplasmic male sterility (CMS-C) is the most commonly used form of CMS in maize hybrid seed production. Restorer of fertility 4 (Rf4), the major fertility restorer gene of CMS-C, is located on chromosome 8S. To positionally clone Rf4, a large F3 population derived from a cross between a non-restorer and restorer (n = 5104) was screened for recombinants and then phenotyped for tassel fertility, resulting in a final map-based cloning interval of 12 kb. Within this 12-kb interval, the only likely candidate for Rf4 was GRMZM2G021276, a basic helix-loop-helix (bHLH) transcription factor with tassel-specific expression. The Rf4 gene product contains a nuclear localization signal and is likely to not interact directly with the mitochondria. Sequence analysis of Rf4 revealed four encoded amino acid substitutions between restoring and non-restoring inbreds, however only one substitution, F187Y, was within the highly conserved bHLH domain. The hypothesis that Rf4 restoration is altered by a single amino acid was tested by using clustered regularly interspaced short palindromic repeat (CRISPR)-CRISPR associated protein 9 (Cas9) homology directed repair (HDR) to create isogenic lines that varied for the F187Y substitution. In a population of these CRISPR-Cas9 edited plants (n = 780) that was phenotyped for tassel fertility, plants containing F187 were completely fertile, indicating fertility restoration, and plants containing Y187 were sterile, indicating lack of fertility restoration. Structural modeling shows that this amino acid residue 187 is located within the four helix bundle core, a critical region for stabilizing dimer conformation and affecting interaction partner selection.

The Interplay between G-quadruplex and Transcription

G4 DNA is a non-canonical DNA structure consisting of a stacked array of Gquartets held together by base pairing between guanine bases. The formation of G4 DNA requires a cluster of guanine-runs within a strand of DNA. Even though the chemistry of this remarkable DNA structure has been under investigation for decades, evidence supporting the biological relevance of G4 DNA has only begun to emerge and point to very important and conserved biological functions. This review will specifically focus on the interplay between transcription and G4 DNA and discuss two alternative but interconnected perspectives. The first part of the review will describe the evidence substantiating the intriguing idea that a shift in DNA structural conformation could be another layer of non-genetic or epigenetic regulator of gene expression and thereby an important determinant of cell fate. The second part will describe the recent genetic studies showing that those genomic loci containing G4 DNA-forming guanine-rich sequences are potential hotspots of genome instability and that the level and orientation of transcription is critical in the materialization of genome instability associated with these sequences.

DNA G-quadruplexes: functional significance in plant and farm animal science

G-quadruplexes (G4s) are non-canonical structures that can be formed in DNA and RNA sequences which carry four short runs of guanines. They are distributed in the whole genome but are enriched in gene promoter regions, gene UTRs and chromosome telomeres. The whole array of their functional roles is not fully explored yet but there is solid evidence supporting their implication in a number of processes like regulation of transcription, replication and telomere organization, among others. During the last decade, there is an increased research interest for G4s that has resulted in a better understanding of their role in several physiological and pathological conditions. On the other hand, these structures are poorly studied in plant species and animals of agricultural interest. Here, we summarize the current methods that are used for studying G4s, we review the studies concerning plants and farm animals and we discuss the advantages of a more thorough inclusion of G4s research in the agricultural sciences.

Sequence Dynamics of Pre-mRNA G-Quadruplexes in Plants

Intramolecular G-quadruplexes (G4s) are secondary structures that may form within G-rich stretches of nucleic acids. Although their presence has been associated with genomic instability and mutagenicity, recent reports suggest their involvement in regulation of diverse cellular events, including transcription and translation. The majority of data regarding G4s stems from mammalian and yeast studies, leaving the plant G4s almost unexplored. Using the publicly available Arabidopsis thaliana and Oryza sativa WGS data, we examined the single nucleotide variability of sequences predicted to form G4s (pG4s) structures. We focused our analysis on protein coding transcripts and compared the results to well-characterized Homo sapiens data. We demonstrate that the overall high variability of pG4s is not uniform and differs between gene structural elements. Specifically, plant AUG-containing pG4s, located within 5'UTR/CDS junctions, are abundant and appear not to be affected by a higher frequency of sequence change, indicating their functional relevance. Furthermore, we show that substitutions lowering the probability of G4s' formation are preferred over neutral or stabilizing modifications.

Bulged and Canonical G-Quadruplex Conformations Determine NDPK Binding Specificity

Guanine-rich DNA strands can adopt tertiary structures known as G-quadruplexes (G4s) that form when Hoogsteen base-paired guanines assemble as planar stacks, stabilized by a central cation like K⁺. In this study, we investigated the conformational heterogeneity of a G-rich sequence from the 5′ untranslated region of the Zea mayshexokinase4 gene. This sequence adopted an extensively polymorphic G-quadruplex, including non-canonical bulged G-quadruplex folds that co-existed in solution. The nature of this polymorphism depended, in part, on the incorporation of different sets of adjacent guanines into a quadruplex core, which permitted the formation of the different conformations. Additionally, we showed that the maize homolog of the human nucleoside diphosphate kinase (NDPK) NM23-H2 protein—ZmNDPK1—specifically recognizes and promotes formation of a subset of these conformations. Heteromorphic G-quadruplexes play a role in microorganisms’ ability to evade the host immune system, so we also discuss how the underlying properties that determine heterogeneity of this sequence could apply to microorganism G4s.

Intensive Distribution of G₂-Quaduplexes in the Pseudorabies Virus Genome and Their Sensitivity to Cations and G-Quadruplex Ligands

Guanine-rich sequences in the genomes of herpesviruses can fold into G-quadruplexes. Compared with the widely-studied G₃-quadruplexes, the dynamic G₂-quadruplexes are more sensitive to the cell microenvironment, but they attract less attention. Pseudorabies virus (PRV) is the model species for the study of the latency and reactivation of herpesvirus in the nervous system. A total of 1722 G₂-PQSs and 205 G₃-PQSs without overlap were identified in the PRV genome. Twelve G₂-PQSs from the CDS region exhibited high conservation in the genomes of the Varicellovirus genus. Eleven G₂-PQSs were 100% conserved in the repeated region of the annotated PRV genomes. There were 212 non-redundant G₂-PQSs in the 3′ UTR and 19 non-redundant G₂-PQSs in the 5′ UTR, which would mediate gene expression in the post-transcription and translation processes. The majority of examined G₂-PQSs formed parallel structures and exhibited different sensitivities to cations and small molecules in vitro. Two G₂-PQSs, respectively, from 3′ UTR of UL5 (encoding helicase motif) and UL9 (encoding sequence-specific ori-binding protein) exhibited diverse regulatory activities with/without specific ligands in vivo. The G-quadruplex ligand, NMM, exhibited a potential for reducing the virulence of the PRV Ea strain. The systematic analysis of the distribution of G₂-PQSs in the PRV genomes could guide further studies of the G-quadruplexes’ functions in the life cycle of herpesviruses.

Above and Beyond Watson and Crick: Guanine Quadruplex Structures and Microbes

Advances in understanding mechanisms of nucleic acids have revolutionized molecular biology and medicine, but understanding of nontraditional nucleic acid conformations is less developed. The guanine quadruplex (G4) alternative DNA structure was first described in the 1960s, but the existence of G4 structures (G4-S) and their participation in myriads of biological functions are still underappreciated. Despite many tools to study G4s and many examples of roles for G4s in eukaryotic molecular processes and issues with uncontrolled G4-S formation, there is relatively little knowledge about the roles of G4-S in viral or prokaryotic systems. This review summarizes the state of the art with regard to G4-S in eukaryotes and their potential roles in human disease before discussing the evidence that G4-S have equivalent importance in affecting viral and bacterial life.

Review: Plant G-quadruplex (G4) motifs in DNA and RNA; abundant, intriguing sequences of unknown function

DNA sequences capable of forming G-quadruplex (G4) structures can be predicted and mapped in plant genomes using computerized pattern search programs. Non-telomeric G4 motifs have recently been found to number in the thousands across many plant species and enriched around gene promoters, prompting speculation that they may represent a newly uncovered and ubiquitous family of cis-acting elements. Comparative analysis shows that monocots exhibit five to ten times higher G4 motif density than eudicots, but the significance of this difference has not been determined. The vast scale and complexity of G4 functions, actual or theoretical, are reviewed in relation to the multiple modes of action and myriad genetic functions for which G4s have been implicated in DNA and RNA. Future experimental strategies and opportunities include identifying plant G4-interactomes, resolving the structures of G4s with and without their binding partners, and defining molecular mechanisms through reporter gene, genetic, or genome editing approaches. Given the global importance of plants for food, clothing, medicine, and energy, together with the potential role of G4 motifs as a widely conserved set of DNA sequences that could coordinate gene regulation, future plant G4 research holds great potential for use in plant improvement strategies.

Quadruplex DNA in long terminal repeats in maize LTR retrotransposons inhibits the expression of a reporter gene in yeast

BACKGROUND

Many studies have shown that guanine-rich DNA sequences form quadruplex structures (G4) in vitro but there is scarce evidence of guanine quadruplexes in vivo. The majority of potential quadruplex-forming sequences (PQS) are located in transposable elements (TEs), especially close to promoters within long terminal repeats of plant LTR retrotransposons.

RESULTS

In order to test the potential effect of G4s on retrotransposon expression, we cloned the long terminal repeats of selected maize LTR retrotransposons upstream of the lacZ reporter gene and measured its transcription and translation in yeast. We found that G4s had an inhibitory effect on translation in vivo since "mutants" (where guanines were replaced by adenines in PQS) showed higher expression levels than wild-types. In parallel, we confirmed by circular dichroism measurements that the selected sequences can indeed adopt G4 conformation in vitro. Analysis of RNA-Seq of polyA RNA in maize seedlings grown in the presence of a G4-stabilizing ligand (NMM) showed both inhibitory as well as stimulatory effects on the transcription of LTR retrotransposons.

CONCLUSIONS

Our results demonstrate that quadruplex DNA located within long terminal repeats of LTR retrotransposons can be formed in vivo and that it plays a regulatory role in the LTR retrotransposon life-cycle, thus also affecting genome dynamics.

G Quadruplex in Plants: A Ubiquitous Regulatory Element and Its Biological Relevance

G quadruplexes (G4) are higher-order DNA and RNA secondary structures formed by G-rich sequences that are built around tetrads of hydrogen-bonded guanine bases. Potential G4 quadruplex sequences have been identified in G-rich eukaryotic non-telomeric and telomeric genomic regions. Upon function, G4 formation is known to involve in chromatin remodeling, gene regulation and has been associated with genomic instability, genetic diseases and cancer progression. The natural role and biological validation of G4 structures is starting to be explored, and is of particular interest for the therapeutic interventions for human diseases. However, the existence and physiological role of G4 DNA and G4 RNA in plants species have not been much investigated yet and therefore, is of great interest for the development of improved crop varieties for sustainable agriculture. In this context, several recent studies suggests that these highly diverse G4 structures in plants can be employed to regulate expression of genes involved in several pathophysiological conditions including stress response to biotic and abiotic stresses as well as DNA damage. In the current review, we summarize the recent findings regarding the emerging functional significance of G4 structures in plants and discuss their potential value in the development of improved crop varieties.

Engineering the expression level of cytosolic nucleoside diphosphate kinase in transgenic Solanum tuberosum roots alters growth, respiration and carbon metabolism

Nucleoside diphosphate kinase (NDPK) is a ubiquitous enzyme that catalyzes the transfer of the γ-phosphate from a donor nucleoside triphosphate to an acceptor nucleoside diphosphate. In this study we used a targeted metabolomic approach and measurement of physiological parameters to report the effects of the genetic manipulation of cytosolic NDPK (NDPK1) expression on physiology and carbon metabolism in potato (Solanum tuberosum) roots. Sense and antisense NDPK1 constructs were introduced in potato using Agrobacterium rhizogenes to generate a population of root clones displaying a 40-fold difference in NDPK activity. Root growth, O₂ uptake, flux of carbon between sucrose and CO₂ , levels of reactive oxygen species and some tricarboxylic acid cycle intermediates were positively correlated with levels of NDPK1 expression. In addition, NDPK1 levels positively affected UDP-glucose and cellulose contents. The activation state of ADP-glucose pyrophosphorylase, a key enzyme in starch synthesis, was higher in antisense roots than in roots overexpressing NDPK1. Further analyses demonstrated that ADP-glucose pyrophosphorylase was more oxidized, and therefore less active, in sense clones than antisense clones. Consequently, antisense NDPK1 roots accumulated more starch and the starch to cellulose ratio was negatively affected by the level of NDPK1. These data support the idea that modulation of NDPK1 affects the distribution of carbon between starch and cellulose biosynthetic pathways.

RNA G-Quadruplexes as Key Motifs of the Transcriptome

G-Quadruplexes are non-canonical secondary structures that can be adopted under physiological conditions by guanine-rich DNA and RNA molecules. They have been reported to occur, and to perform multiple biological functions, in the genomes and transcriptomes of many species, including humans. This chapter focuses specifically on RNA G-quadruplexes and reviews the most recent discoveries in the field, as well as addresses the upcoming challenges researchers studying these structures face.

Discovery of a Family of Genomic Sequences Which Interact Specifically with the c-MYC Promoter to Regulate c-MYC Expression

G-quadruplex forming sequences are particularly enriched in the promoter regions of eukaryotic genes, especially of oncogenes. One of the most well studied G-quadruplex forming sequences is located in the nuclease hypersensitive element (NHE) III₁ of the c-MYC promoter region. The oncoprotein c-MYC regulates a large array of genes which play important roles in growth regulation and metabolism. It is dysregulated in >70% of human cancers. The silencer NHEIII₁ located upstream of the P1 promoter regulates up-to 80% of c-MYC transcription and includes a G-quadruplex structure (Pu27) that is required for promoter inhibition. We have identified, for the first time, a family of seventeen G-quadruplex-forming motifs with >90% identity with Pu27, located on different chromosomes throughout the human genome, some found near or within genes involved in stem cell maintenance or neural cell development. Notably, all members of the Pu27 family interact specifically with NHEIII₁ sequence, in vitro. Crosslinking studies demonstrate that Pu27 oligonucleotide binds specifically to the C-rich strand of the NHEIII₁ resulting in the G-quadruplex structure stabilization. Pu27 homologous sequences (Pu27-HS) significantly inhibit leukemic cell lines proliferation in culture. Exposure of U937 cells to the Pu27-HS induces cell growth inhibition associated with cell cycle arrest that is most likely due to downregulation of c-MYC expression at the RNA and/or protein levels. Expression of SOX2, another gene containing a Pu27-HS, was affected by Pu27-HS treatment as well. Our data suggest that the oligonucleotides encoding the Pu27 family target complementary DNA sequences in the genome, including those of the c-MYC and SOX2 promoters. This effect is most likely cell type and cell growth condition dependent. The presence of genomic G-quadruplex-forming sequences homologous to Pu27 of c-MYC silencer and the fact that they interact specifically with the parent sequence suggest a common regulatory mechanism for genes whose promoters contain these sequences.

G-Quadruplexes Involving Both Strands of Genomic DNA Are Highly Abundant and Colocalize with Functional Sites in the Human Genome

The G-quadruplex is a non-canonical DNA structure biologically significant in DNA replication, transcription and telomere stability. To date, only G4s with all guanines originating from the same strand of DNA have been considered in the context of the human nuclear genome. Here, I discuss interstrand topological configurations of G-quadruplex DNA, consisting of guanines from both strands of genomic DNA; an algorithm is presented for predicting such structures. I have identified over 550,000 non-overlapping interstrand G-quadruplex forming sequences in the human genome—significantly more than intrastrand configurations. Functional analysis of interstrand G-quadruplex sites shows strong association with transcription initiation, the results are consistent with the XPB and XPD transcriptional helicases binding only to G-quadruplex DNA with interstrand topology. Interstrand quadruplexes are also enriched in origin of replication sites. Several topology classes of interstrand quadruplex-forming sequences are possible, and different topologies are enriched in different types of structural elements. The list of interstrand quadruplex forming sequences, and the computer program used for their prediction are available at the web address http://moment.utmb.edu/allquads.

MaizeGDB update: new tools, data and interface for the maize model organism database

MaizeGDB is a highly curated, community-oriented database and informatics service to researchers focused on the crop plant and model organism Zea mays ssp. mays. Although some form of the maize community database has existed over the last 25 years, there have only been two major releases. In 1991, the original maize genetics database MaizeDB was created. In 2003, the combined contents of MaizeDB and the sequence data from ZmDB were made accessible as a single resource named MaizeGDB. Over the next decade, MaizeGDB became more sequence driven while still maintaining traditional maize genetics datasets. This enabled the project to meet the continued growing and evolving needs of the maize research community, yet the interface and underlying infrastructure remained unchanged. In 2015, the MaizeGDB team completed a multi-year effort to update the MaizeGDB resource by reorganizing existing data, upgrading hardware and infrastructure, creating new tools, incorporating new data types (including diversity data, expression data, gene models, and metabolic pathways), and developing and deploying a modern interface. In addition to coordinating a data resource, the MaizeGDB team coordinates activities and provides technical support to the maize research community. MaizeGDB is accessible online at http://www.maizegdb.org.

Guanine quadruplex structures localize to heterochromatin

Increasing amounts of data support a role for guanine quadruplex (G4) DNA and RNA structures in various cellular processes. We stained different organisms with monoclonal antibody 1H6 specific for G4 DNA. Strikingly, immuno-electron microscopy showed exquisite specificity for heterochromatin. Polytene chromosomes from Drosophila salivary glands showed bands that co-localized with heterochromatin proteins HP1 and the SNF2 domain-containing protein SUUR. Staining was retained in SUUR knock-out mutants but lost upon overexpression of SUUR. Somatic cells in Macrostomum lignano were strongly labeled, but pluripotent stem cells labeled weakly. Similarly, germline stem cells in Drosophila ovaries were weakly labeled compared to most other cells. The unexpected presence of G4 structures in heterochromatin and the difference in G4 staining between somatic cells and stem cells with germline DNA in ciliates, flatworms, flies and mammals point to a conserved role for G4 structures in nuclear organization and cellular differentiation.