QGRS Mapper: a web-based server for predicting G-quadruplexes in nucleotide sequences

An RNA-binding Protein, Lin28, Recognizes and Remodels G-quartets in the MicroRNAs (miRNAs) and mRNAs It Regulates

Lin28 is an evolutionarily conserved RNA-binding protein that inhibits processing of pre-let-7 microRNAs (miRNAs) and regulates translation of mRNAs that control developmental timing, pluripotency, metabolism, and tumorigenesis. The RNA features that mediate Lin28 binding to the terminal loops of let-7 pre-miRNAs and to Lin28-responsive elements (LREs) in mRNAs are not well defined. Here we show that Lin28 target datasets are enriched for RNA sequences predicted to contain stable planar structures of 4 guanines known as G-quartets (G4s). The imino NMR spectra of pre-let-7 loops and LREs contain resonances characteristic of G4 hydrogen bonds. These sequences bind to a G4-binding fluorescent dye, N-methyl-mesoporphyrin IX (NMM). Mutations and truncations in the RNA sequence that prevent G4 formation also prevent Lin28 binding. The addition of Lin28 to a pre-let-7 loop or an LRE reduces G4 resonance intensity and NMM binding, suggesting that Lin28 may function to remodel G4s. Further, we show that NMM inhibits Lin28 binding. Incubation of a human embryonal carcinoma cell line with NMM reduces its stem cell traits. In particular it increases mature let-7 levels, decreases OCT4, HMGA1, CCNB1, CDK4, and Lin28A protein, decreases sphere formation, and inhibits colony formation. Our results suggest a previously unknown structural feature of Lin28 targets and a new strategy for manipulating Lin28 function.

Mapping of Replication Origins in the X Inactivation Center of Vole Microtus levis Reveals Extended Replication Initiation Zone

DNA replication initiates at specific positions termed replication origins. Genome-wide studies of human replication origins have shown that origins are organized into replication initiation zones. However, only few replication initiation zones have been described so far. Moreover, few origins were mapped in other mammalian species besides human and mouse. Here we analyzed pattern of short nascent strands in the X inactivation center (XIC) of vole Microtus levis in fibroblasts, trophoblast stem cells, and extraembryonic endoderm stem cells and confirmed origins locations by ChIP approach. We found that replication could be initiated in a significant part of XIC. We also analyzed state of XIC chromatin in these cell types. We compared origin localization in the mouse and vole XIC. Interestingly, origins associated with gene promoters are conserved in these species. The data obtained allow us to suggest that the X inactivation center of M. levis is one extended replication initiation zone.

G-quadruplexes mediate local translation in neurons

There has recently been a huge increase in interest in the formation of stable G-quadruplex structures in mRNAs and their functional significance. In neurons, local translation of mRNA is essential for normal neuronal behaviour. It has been discovered that local translation of specific mRNAs encoding some of the best known synaptic proteins is dependent on the presence of a G-quadruplex. The recognition of G-quadruplexes in mRNAs, their transport as repressed complexes and the control of their translation at their subcellular destinations involves a diversity of proteins, including those associated with disease pathologies. This is an exciting field, with rapid improvements to our knowledge and understanding. Here, we discuss some of the recent work on how G-quadruplexes mediate local translation in neurons.

Sites of instability in the human TCF3 (E2A) gene adopt G-quadruplex DNA structures in vitro

The formation of highly stable four-stranded DNA, called G-quadruplex (G4), promotes site-specific genome instability. G4 DNA structures fold from repetitive guanine sequences, and increasing experimental evidence connects G4 sequence motifs with specific gene rearrangements. The human transcription factor 3 (TCF3) gene (also termed E2A) is subject to genetic instability associated with severe disease, most notably a common translocation event t(1;19) associated with acute lymphoblastic leukemia. The sites of instability in TCF3 are not randomly distributed, but focused to certain sequences. We asked if G4 DNA formation could explain why TCF3 is prone to recombination and mutagenesis. Here we demonstrate that sequences surrounding the major t(1;19) break site and a region associated with copy number variations both contain G4 sequence motifs. The motifs identified readily adopt G4 DNA structures that are stable enough to interfere with DNA synthesis in physiological salt conditions in vitro. When introduced into the yeast genome, TCF3 G4 motifs promoted gross chromosomal rearrangements in a transcription-dependent manner. Our results provide a molecular rationale for the site-specific instability of human TCF3, suggesting that G4 DNA structures contribute to oncogenic DNA breaks and recombination.

Duplex stem-loop-containing quadruplex motifs in the human genome: a combined genomic and structural study

Duplex stem-loops and four-stranded G-quadruplexes have been implicated in (patho)biological processes. Overlap of stem-loop- and quadruplex-forming sequences could give rise to quadruplex-duplex hybrids (QDH), which combine features of both structural forms and could exhibit unique properties. Here, we present a combined genomic and structural study of stem-loop-containing quadruplex sequences (SLQS) in the human genome. Based on a maximum loop length of 20 nt, our survey identified 80 307 SLQS, embedded within 60 172 unique clusters. Our analysis suggested that these should cover close to half of total SLQS in the entire genome. Among these, 48 508 SLQS were strand-specifically located in genic/promoter regions, with the majority of genes displaying a low number of SLQS. Notably, genes containing abundant SLQS clusters were strongly associated with brain tissues. Enrichment analysis of SLQS-positive genes and mapping of SLQS onto transcriptional/mutagenesis hotspots and cancer-associated genes, provided a statistical framework supporting the biological involvements of SLQS. In vitro formation of diverse QDH by selective SLQS hits were successfully verified by nuclear magnetic resonance spectroscopy. Folding topologies of two SLQS were elucidated in detail. We also demonstrated that sequence changes at mutation/single-nucleotide polymorphism loci could affect the structural conformations adopted by SLQS. Thus, our predicted SLQS offer novel insights into the potential involvement of QDH in diverse (patho)biological processes and could represent novel regulatory signals.

Short loop length and high thermal stability determine genomic instability induced by G-quadruplex-forming minisatellites

G-quadruplexes (G4) are polymorphic four-stranded structures formed by certain G-rich nucleic acids, with various biological roles. However, structural features dictating their formation and/or function in vivo are unknown. In S. cerevisiae, the pathological persistency of G4 within the CEB1 minisatellite induces its rearrangement during leading-strand replication. We now show that several other G4-forming sequences remain stable. Extensive mutagenesis of the CEB25 minisatellite motif reveals that only variants with very short (≤ 4 nt) G4 loops preferentially containing pyrimidine bases trigger genomic instability. Parallel biophysical analyses demonstrate that shortening loop length does not change the monomorphic G4 structure of CEB25 variants but drastically increases its thermal stability, in correlation with the in vivo instability. Finally, bioinformatics analyses reveal that the threat for genomic stability posed by G4 bearing short pyrimidine loops is conserved in C. elegans and humans. This work provides a framework explanation for the heterogeneous instability behavior of G4-forming sequences in vivo, highlights the importance of structure thermal stability, and questions the prevailing assumption that G4 structures with short or longer loops are as likely to form in vivo.

Non-coding RNA Generated following Lariat Debranching Mediates Targeting of AID to DNA

Transcription through immunoglobulin switch (S) regions is essential for class switch recombination (CSR), but no molecular function of the transcripts has been described. Likewise, recruitment of activation-induced cytidine deaminase (AID) to S regions is critical for CSR; however, the underlying mechanism has not been fully elucidated. Here, we demonstrate that intronic switch RNA acts in trans to target AID to S region DNA. AID binds directly to switch RNA through G-quadruplexes formed by the RNA molecules. Disruption of this interaction by mutation of a key residue in the putative RNA-binding domain of AID impairs recruitment of AID to S region DNA, thereby abolishing CSR. Additionally, inhibition of RNA lariat processing leads to loss of AID localization to S regions and compromises CSR; both defects can be rescued by exogenous expression of switch transcripts in a sequence-specific manner. These studies uncover an RNA-mediated mechanism of targeting AID to DNA.

Retargeting pre-existing human antibodies to a bacterial pathogen with an alpha-Gal conjugated aptamer

The ever-increasing threat of multi-drug resistant bacterial infections has spurred renewed interest in alternative approaches to classical antibiotic therapy. In contrast to other mammals, humans do not express the galactose-α-1,3-galactosyl-β-1,4-N-acetyl-glucosamine (α-Gal) epitope. As a result of exposure of humans to α-Gal in the environment, a large proportion of circulating antibodies are specific for the trisaccharide. In this study, we examine whether these anti-Gal antibodies can be recruited and redirected to exert anti-bacterial activity. We show that a specific DNA aptamer conjugated to an α-Gal epitope at its 5' end, herein termed an alphamer, can bind to group A Streptococcus (GAS) bacteria by recognition of a conserved region of the surface-anchored M protein. The anti-GAS alphamer was shown to recruit anti-Gal antibodies to the streptococcal surface in an α-Gal-specific manner, elicit uptake and killing of the bacteria by human phagocytes, and slow growth of invasive GAS in human whole blood. These studies provide a first in vitro proof of concept that alphamers have the potential to redirect pre-existing antibodies to bacteria in a specific manner and trigger an immediate antibacterial immune response. Further validation of this novel therapeutic approach of applying α-Gal technology in in vivo models of bacterial infection is warranted.

KEY MESSAGES

. α-Gal-tagged aptamers lead to GAS opsonization with anti-Gal antibodies. . α-Gal-tagged aptamers confer phagocytosis and killing of GAS cells by human phagocytes. . α-Gal-tagged aptamers reduces replication of GAS in human blood. . α-Gal-tagged aptamers may have the potential to be used as novel passive immunization drugs.

Targeted Detection of G-Quadruplexes in Cellular RNAs

The G-quadruplex (G4) is a non-canonical nucleic acid structure which regulates important cellular processes. RNA G4s have recently been shown to exist in human cells and be biologically significant. Described herein is a new approach to detect and map RNA G4s in cellular transcripts. This method exploits the specific control of RNA G4-cation and RNA G4-ligand interactions during reverse transcription, by using a selective reverse transcriptase to monitor RNA G4-mediated reverse transcriptase stalling (RTS) events. Importantly, a ligation-amplification strategy is coupled with RTS, and enables detection and mapping of G4s in important, low-abundance cellular RNAs. Strong evidence is provided for G4 formation in full-length cellular human telomerase RNA, offering important insights into its cellular function.

Isolation of single-stranded DNA aptamers that distinguish influenza virus hemagglutinin subtype H1 from H5

Surface protein hemagglutinin (HA) mediates the binding of influenza virus to host cell receptors containing sialic acid, facilitating the entry of the virus into host cells. Therefore, the HA protein is regarded as a suitable target for the development of influenza virus detection devices. In this study, we isolated single-stranded DNA (ssDNA) aptamers binding to the HA1 subunit of subtype H1 (H1-HA1), but not to the HA1 subunit of subtype H5 (H5-HA1), using a counter-systematic evolution of ligands by exponential enrichment (counter-SELEX) procedure. Enzyme-linked immunosorbent assay and surface plasmon resonance studies showed that the selected aptamers bind tightly to H1-HA1 with dissociation constants in the nanomolar range. Western blot analysis demonstrated that the aptamers were binding to H1-HA1 in a concentration-dependent manner, yet were not binding to H5-HA1. Interestingly, the selected aptamers contained G-rich sequences in the central random nucleotides region. Further biophysical analysis showed that the G-rich sequences formed a G-quadruplex structure, which is a distinctive structure compared to the starting ssDNA library. Using flow cytometry analysis, we found that the aptamers did not bind to the receptor-binding site of H1-HA1. These results indicate that the selected aptamers that distinguish H1-HA1 from H5-HA1 can be developed as unique probes for the detection of the H1 subtype of influenza virus.

A stable RNA G-quadruplex within the 5'-UTR of Arabidopsis thaliana ATR mRNA inhibits translation

Guanine quadruplex structures (GQSs) play important roles in the regulation of gene expression and cellular processes. Recent studies provide strong evidence for the formation and function of DNA and RNA GQSs in human cells. However, whether GQSs form and are functional in plants remains essentially unexplored. On the basis of circular dichroism (CD)-detected titration, UV-detected melting, in-line probing (ILP) and reporter gene assay studies, we report the first example of a plant RNA GQS that inhibits translation. This GQS is located within the 5'-UTR of the ATAXIA TELANGIECTASIA-MUTATED AND RAD3-RELATED (ATR) mRNA of Arabidopsis thaliana (mouse-ear cress). We show that this GQS is highly stable and is thermodynamically favoured over a competing hairpin structure in the 5'-UTR at physiological K⁺ and Mg²⁺ concentrations. Results from ILP reveal the secondary structure of the RNA and support formation of the GQS in vitro in the context of the complete 5'-UTR. Transient reporter gene assays performed in living plants reveal that the GQS inhibits translation but not transcription, implicating this GQS as a translational repressor in vivo. Our results provide the first complete demonstration of the formation and function of a regulatory RNA GQS in plants and open new avenues to explore potential functional roles of GQS in the plant kingdom.

Aptamers targeting different functional groups of 17β-estradiol

Aptamers, short synthetic ssDNA or RNA molecules with a specific three-dimensional structure, are promising recognition elements in biosensor technology. In vitro generation of aptamers with high sensitivity and specificity toward a broad range of analytes has been achieved using the systematic evolution of ligands by exponential enrichment (SELEX) process. This iterative pathway of aptamer generation consists of sequential positive and counterselection steps. The present research aimed to select two sets of ssDNA aptamers which both are able to bind to different functional groups on the cyclopentanoperhydrophenanthrene ring of 17β-estradiol (E2). By repetitively switching between positive selection steps using E2 as target molecule and counterselection steps with nortestosterone as countermolecule, aptamers were successfully selected against the hydroxylated aromatic A ring of E2. Additionally, an aptamer which binds the upper segments of the B, C and D ring of the cyclopentanoperhydrophenanthrene ring of E2 was generated after repetitively swapping between positive selection steps with E2 as target molecule and counterselection steps with dexamethasone as countermolecule. Epitope specificity of the aptamers was demonstrated by evaluating their binding responses toward a number of steroid hormones structurally related to E2. The selected aptamers with affinities for different functional groups of E2 can potentially be applied to develop a cross-reactive aptasensor. This aptasensor introduces a promising tool for the future of in-field real-time monitoring of a wide range of steroid hormones.

Two-quartet G-quadruplexes formed by DNA sequences containing four contiguous GG runs

The DNA sequence containing four contiguous GG runs (G2NxG2NyG2NzG2, G2 sequence) has the potential to form a two-quartet G-quadruplex. However, the prevalence, structure, and function of G2 sequences have not been well-studied. Here, bioinformatics analysis reveals the abundance of G2 sequences in the human genome and their enrichment in promoter regions. The density of G2 sequences in the genome and promoters is much higher than that of the G3 sequence (G3NxG3NyG3NzG3). Experiments show that the conformations and thermal stabilities of the two-quartet G-quadruplexes of G2 sequences are highly sensitive to the length and composition of the loops. Among the two-quartet G-quadruplexes, the parallel G-quadruplex with a loop length of 1 and the antiparallel G-quadruplex with a loop length of 3 show high thermal stabilities. Additionally, the stable parallel G-quadruplexes are stacked into intermolecular higher-order structures. This work determines the prevalence of G2 sequences in the human genome and demonstrates that the G-quadruplex structures for certain loop lengths and compositions may be stable in vivo. Thus, more attention should be paid to the structure and function of the two-quartet G-quadruplex.

Characterizing and controlling intrinsic biases of lambda exonuclease in nascent strand sequencing reveals phasing between nucleosomes and G-quadruplex motifs around a subset of human replication origins

Nascent strand sequencing (NS-seq) is used to discover DNA replication origins genome-wide, allowing identification of features for their specification. NS-seq depends on the ability of lambda exonuclease (λ-exo) to efficiently digest parental DNA while leaving RNA-primer protected nascent strands intact. We used genomics and biochemical approaches to determine if λ-exo digests all parental DNA sequences equally. We report that λ-exo does not efficiently digest G-quadruplex (G4) structures in a plasmid. Moreover, λ-exo digestion of nonreplicating genomic DNA (LexoG0) enriches GC-rich DNA and G4 motifs genome-wide. We used LexoG0 data to control for nascent strand–independent λ-exo biases in NS-seq and validated this approach at the rDNA locus. The λ-exo–controlled NS-seq peaks are not GC-rich, and only 35.5% overlap with 6.8% of all G4s, suggesting that G4s are not general determinants for origin specification but may play a role for a subset. Interestingly, we observed a periodic spacing of G4 motifs and nucleosomes around the peak summits, suggesting that G4s may position nucleosomes at this subset of origins. Finally, we demonstrate that use of Na⁺ instead of K⁺ in the λ-exo digestion buffer reduced the effect of G4s on λ-exo digestion and discuss ways to increase both the sensitivity and specificity of NS-seq.

G-quadruplexes in pathogens: a common route to virulence control?

DNA can form several secondary structures besides the classic double helix: one that has received much attention in recent years is the G-quadruplex (G4). This is a stable four-stranded structure formed by the stacking of quartets of guanine bases. Recent work has convincingly shown that G4s can form in vivo as well as in vitro and can affect both replication and transcription of DNA. They also play important roles at G-rich telomeres. Now, a spate of exciting reports has begun to reveal roles for G4 structures in virulence processes in several important microbial pathogens of humans. Interestingly, these come from a range of kingdoms—bacteria and protozoa as well as viruses—and all facilitate immune evasion in different ways. In particular, roles for G4s have been posited in the antigenic variation systems of bacteria and protozoa, as well as in the silencing of at least two major human viruses, human immunodeficiency virus (HIV) and Epstein-Barr virus (EBV). Although antigenic variation and the silencing of latent viruses are quite distinct from one another, both are routes to immune evasion and the maintenance of chronic infections. Thus, highly disparate pathogens can use G4 motifs to control DNA/RNA dynamics in ways that are relevant to common virulence phenotypes. This review explores the evidence for G4 biology in such processes across a range of important human pathogens.

A potassium ion-dependent RNA structural switch regulates human pre-miRNA 92b maturation

MicroRNAs (miRNAs) are an important set of oligonucleotide sequences with a biogenesis that involves Dicer-mediated cleavage as a critical step. Dicer cleaves the precursor miRNA (pre-miRNA) stem-loop structure to produce the mature miRNA. Using bioinformatics analysis, we discovered the presence of putative G-quadruplex (GQ)-forming sequences in about 16% of pre-miRNAs. We report the existence of a GQ as an alternative to the canonical stem-loop structure in the clinically important human pre-miRNA 92b. GQ formation led to unwinding of the stem-loop structure imparting resistance to Dicer-mediated cleavage both in vitro and in vivo. A potential K(+) ion-dependent equilibrium between GQ and the stem-loop structure has the ability to regulate the Dicer-mediated maturation of pre-miRNA 92b, which consequently affects target gene silencing. These findings unravel a new mechanism of regulation in pre-miRNA maturation, albeit at the RNA structure level.

Aggregating distyrylpyridinium dye as a bimodal structural probe for G-quadruplex DNA

A coumarin-substituted distyrylpyridinium dye BCVP is proposed as a universal colorimetric and fluorescent probe for the detection of G-quadruplex DNA structures.

Using hidden Markov models to investigate G-quadruplex motifs in genomic sequences

BACKGROUND

G-quadruplexes are four-stranded structures formed in guanine-rich nucleotide sequences. Several functional roles of DNA G-quadruplexes have so far been investigated, where their putative functional roles during DNA replication and transcription have been suggested. A necessary condition for G-quadruplex formation is the presence of four regions of tandem guanines called G-runs and three nucleotide subsequences called loops that connect G-runs. A simple computational way to detect potential G-quadruplex regions in a given genomic sequence is pattern matching with regular expression. Although many putative G-quadruplex motifs can be found in most genomes by the regular expression-based approach, the majority of these sequences are unlikely to form G-quadruplexes because they are unstable as compared with canonical double helix structures.

RESULTS

Here we present elaborate computational models for representing DNA G-quadruplex motifs using hidden Markov models (HMMs). Use of HMMs enables us to evaluate G-quadruplex motifs quantitatively by a probabilistic measure. In addition, the parameters of HMMs can be trained by using experimentally verified data. Computational experiments in discriminating between positive and negative G-quadruplex sequences as well as reducing putative G-quadruplexes in the human genome were carried out, indicating that HMM-based models can discern bona fide G-quadruplex structures well and one of them has the possibility of reducing false positive G-quadruplexes predicted by existing regular expression-based methods. Furthermore, our results show that one of our models can be specialized to detect G-quadruplex sequences whose functional roles are expected to be involved in DNA transcription.

CONCLUSIONS

The HMM-based method along with the conventional pattern matching approach can contribute to reducing costly and laborious wet-lab experiments to perform functional analysis on a given set of potential G-quadruplexes of interest. The C++ and Perl programs are available at http://tcs.cira.kyoto-u.ac.jp/~ykato/program/g4hmm/.

G-quadruplex structures and CpG methylation cause drop-out of the maternal allele in polymerase chain reaction amplification of the imprinted MEST gene promoter

We observed apparent non-Mendelian behaviour of alleles when genotyping a region in a CpG island at the 5′ end of the maternally imprinted human MEST isoform. This region contains three single nucleotide polymorphisms (SNPs) in total linkage disequilibrium, such that only two haplotypes occur in the human population. Only one haplotype was detectable in each subject, never both, despite the use of multiple primers and several genotyping methods. We observed that this region contains motifs capable of forming several G-quadruplex structures. Circular dichroism spectroscopy and native polyacrylamide gel electrophoresis confirmed that at least three G-quadruplexes form in vitro in the presence of potassium ions, and one of these structures has a T_m of greater than 99°C in polymerase chain reaction (PCR) buffer. We demonstrate that it is the methylated maternal allele that is always lost during PCR amplification, and that formation of G-quadruplexes and presence of methylated cytosines both contributed to this phenomenon. This observed parent-of-origin specific allelic drop-out has important implications for analysis of imprinted genes in research and diagnostic settings.

MOV10 and FMRP regulate AGO2 association with microRNA recognition elements

The fragile X mental retardation protein FMRP regulates translation of its bound mRNAs through incompletely defined mechanisms. FMRP has been linked to the microRNA pathway, and we show here that it associates with the RNA helicase MOV10, also associated with the microRNA pathway. FMRP associates with MOV10 directly and in an RNA-dependent manner and facilitates MOV10's association with RNAs in brain and cells, suggesting a cooperative interaction. We identified the RNAs recognized by MOV10 using RNA immunoprecipitation and iCLIP. Examination of the fate of MOV10 on RNAs revealed a dual function for MOV10 in regulating translation: it facilitates microRNA-mediated translation of some RNAs, but it also increases expression of other RNAs by preventing AGO2 function. The latter subset was also bound by FMRP in close proximity to the MOV10 binding site, suggesting that FMRP prevents MOV10-mediated microRNA suppression. We have identified a mechanism for FMRP-mediated translational regulation through its association with MOV10.

Targeting G-quadruplex nucleic acids with heterocyclic alkaloids and their derivatives

G-Quadruplex nucleic acids or G-quadruplexes (G4s) are four-stranded DNA or RNA secondary structures that are formed in guanine-rich sequences. They are widely distributed in functional regions of the human genome, such as telomeres, ribosomal DNA (rDNA), transcription start sites, promoter regions and untranslated regions of mRNA, suggesting that G-quadruplex structures may play a pivotal role in the control of a variety of cellular processes. G-Quadruplexes are viewed as valid therapeutic targets in human cancer diseases. Small molecules, from naturally occurring to synthetic, are exploited to specifically target G-quadruplexes and have proven to be a new class of anticancer agents. Notably, alkaloids are an important source of G-quadruplex ligands and have significant bioactivities in anticancer therapy. In this review, the authors provide a brief, up-to-date summary of heterocyclic alkaloids and their derivatives targeting G-quadruplexes.

G-quadruplexes in viruses: function and potential therapeutic applications

G-rich nucleic acids can form non-canonical G-quadruplex structures (G4s) in which four guanines fold in a planar arrangement through Hoogsteen hydrogen bonds. Although many biochemical and structural studies have focused on DNA sequences containing successive, adjacent guanines that spontaneously fold into G4s, evidence for their in vivo relevance has recently begun to accumulate. Complete sequencing of the human genome highlighted the presence of ∼300,000 sequences that can potentially form G4s. Likewise, the presence of putative G4-sequences has been reported in various viruses genomes [e.g., Human immunodeficiency virus (HIV-1), Epstein-Barr virus (EBV), papillomavirus (HPV)]. Many studies have focused on telomeric G4s and how their dynamics are regulated to enable telomere synthesis. Moreover, a role for G4s has been proposed in cellular and viral replication, recombination and gene expression control. In parallel, DNA aptamers that form G4s have been described as inhibitors and diagnostic tools to detect viruses [e.g., hepatitis A virus (HAV), EBV, cauliflower mosaic virus (CaMV), severe acute respiratory syndrome virus (SARS), simian virus 40 (SV40)]. Here, special emphasis will be given to the possible role of these structures in a virus life cycle as well as the use of G4-forming oligonucleotides as potential antiviral agents and innovative tools.

Abundance, arrangement, and function of sequence motifs in the chicken promoters

Background

Eukaryotic promoters are regions containing various sequence motifs necessary to control gene transcription. Much evidence has emerged showing that structural and/or contextual changes in regulatory elements can critically affect cis-regulatory activity. As sequence motifs can be key factors in maintaining complex promoter architectures, one effective approach to further understand the evolution of promoter regions in vertebrates is to compare the abundance and distribution patterns of sequence motifs in these regions between divergent species. When compared with mammals, the chicken (Gallus gallus) has a very different genome composition and sufficient genomic information to make it a good model for the exploration of promoter structure and evolution.

Results

More than 10% of chicken genes contained short tandem repeat (STR) in the region 2 kb upstream of promoters, but the total number of STRs observed in chicken is approximately half of that detected in human promoters. In terms of the STR motif frequencies, chicken promoter regions were more similar to other avian and mammalian promoters than these were to the entire chicken genome. Unlike other STRs, nearly half of the trinucleotide repeats found in promoters partly or entirely overlapped with CpG islands, indicating potential association with nucleosome positions. Moreover, the chicken promoters are abundant with sequence motifs such as poly-A, poly-G and G-quadruplexes, especially in the core region, that are otherwise rare in the genome. Most of sequence motifs showed strong functional enrichment for particular gene ontology (GO) categories, indicating roles in regulation of transcription and gene expression, as well as immune response and cognition.

Conclusions

Chicken promoter regions share some, but not all, of the structural features observed in mammalian promoters. The findings presented here provide empirical evidence suggesting that the frequencies and locations of STR motifs have been conserved through promoter evolution in a lineage-specific manner. Correlation analysis between GO categories and sequence motifs suggests motif-specific constraints acting on gene function.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-900) contains supplementary material, which is available to authorized users.

Impact of G-quadruplex structures and intronic polymorphisms rs17878362 and rs1642785 on basal and ionizing radiation-induced expression of alternative p53 transcripts

G-quadruplex (G4) structures in intron 3 of the p53 pre-mRNA modulate intron 2 splicing, altering the balance between the fully spliced p53 transcript (FSp53, encoding full-length p53) and an incompletely spliced transcript retaining intron 2 (p53I2 encoding the N-terminally truncated Δ40p53 isoform). The nucleotides forming G4s overlap the polymorphism rs17878362 (A1 wild-type allele, A2 16-base pair insertion) which is in linkage disequilibrium with rs1642785 in intron 2 (c.74+38 G>C). Biophysical and biochemical analyses show rs17878362 A2 alleles form similar G4 structures as A1 alleles although their position is shifted with respect to the intron 2 splice acceptor site. In addition basal FSp53 and p53I2 levels showed allele specific differences in both p53-null cells transfected with reporter constructs or lymphoblastoid cell lines. The highest FSp53 and p53I2 levels were associated with combined rs1642785-GG/rs17878362-A1A1 alleles, whereas the presence of rs1642785-C with either rs17878362 allele was associated with lower p53 pre-mRNA, total TP53, FSp53 and p53I2 levels, due to the lower stability of transcripts containing rs1642785-C. Treatment of lymphoblastoid cell with the G4 binding ligands 360A or PhenDC3 or with ionizing radiation increased FSp53 levels only in cells with rs17878362 A1 alleles, suggesting that under this G4 configuration full splicing is favoured. These results demonstrate the complex effects of intronic TP53 polymorphisms on G4 formation and identify a new role for rs1642785 on mRNA splicing and stability, and thus on the differential expression of isoform-specific transcripts of the TP53 gene.

DNA and RNA quadruplex-binding proteins

Four-stranded DNA structures were structurally characterized in vitro by NMR, X-ray and Circular Dichroism spectroscopy in detail. Among the different types of quadruplexes (i-Motifs, minor groove quadruplexes, G-quadruplexes, etc.), the best described are G-quadruplexes which are featured by Hoogsteen base-paring. Sequences with the potential to form quadruplexes are widely present in genome of all organisms. They are found often in repetitive sequences such as telomeric ones, and also in promoter regions and 5' non-coding sequences. Recently, many proteins with binding affinity to G-quadruplexes have been identified. One of the initially portrayed G-rich regions, the human telomeric sequence (TTAGGG)n, is recognized by many proteins which can modulate telomerase activity. Sequences with the potential to form G-quadruplexes are often located in promoter regions of various oncogenes. The NHE III1 region of the c-MYC promoter has been shown to interact with nucleolin protein as well as other G-quadruplex-binding proteins. A number of G-rich sequences are also present in promoter region of estrogen receptor alpha. In addition to DNA quadruplexes, RNA quadruplexes, which are critical in translational regulation, have also been predicted and observed. For example, the RNA quadruplex formation in telomere-repeat-containing RNA is involved in interaction with TRF2 (telomere repeat binding factor 2) and plays key role in telomere regulation. All these fundamental examples suggest the importance of quadruplex structures in cell processes and their understanding may provide better insight into aging and disease development.

DDX3X, the X homologue of AZFa gene DDX3Y, expresses a complex pattern of transcript variants only in the male germ line

DDX3X, the functional X homologue of the major AZFa gene, DDX3Y, belongs to the highly conserved PL10-subfamily of DEAD-box RNA helicase genes which are functionally conserved from yeast to man. They are mainly involved in cell cycle control and translation initiation control of gene transcripts with long 5'UTR extensions containing complex secondary structures. Interestingly, in humans both gene copies were found to be expressed at different phases of human spermatogenesis. Whereas DDX3Y transcripts are translated only in premeiotic male germ cells, the DDX3X protein is expressed only in postmeiotic spermatids. In this study, we found that the major class of DDX3X transcripts in human testis become activated first after meiosis and at a specific core promoter not active in somatic tissues and not present upstream of the DDX3Y homologue. Two alternative 5'UTR transcript lengths are subsequently produced by an additional testis-specific 5'UTR splicing event. Both transcripts are mainly processed for polyadenylation in their proximal 3'UTR. A minor transcript class starting at the same male germ line-specific core promoter produces primary transcripts with an extremely long 3'UTR (∼ 17 kb), which is subsequently spliced at distinct sites resulting in six short 3'UTR splice variants (I-VI). Comparative analyses of the DDX3X transcripts in mouse and primates revealed that this complex pattern of male germ line-specific transcript variants first evolved in primates. Our data thus suggest complex translational control mechanism(s) for the human DDX3X gene locus functioning only in the male germ line and resulting in expression of its protein only in the postmeiotic spermatids.

Lentivector integration sites in ependymal cells from a model of metachromatic leukodystrophy: non-B DNA as a new factor influencing integration

The blood–brain barrier controls the passage of molecules from the blood into the central nervous system (CNS) and is a major challenge for treatment of neurological diseases. Metachromatic leukodystrophy is a neurodegenerative lysosomal storage disease caused by loss of arylsulfatase A (ARSA) activity. Gene therapy via intraventricular injection of a lentiviral vector is a potential approach to rapidly and permanently deliver therapeutic levels of ARSA to the CNS. We present the distribution of integration sites of a lentiviral vector encoding human ARSA (LV-ARSA) in murine brain choroid plexus and ependymal cells, administered via a single intracranial injection into the CNS. LV-ARSA did not exhibit a strong preference for integration in or near actively transcribed genes, but exhibited a strong preference for integration in or near satellite DNA. We identified several genomic hotspots for LV-ARSA integration and identified a consensus target site sequence characterized by two G-quadruplex-forming motifs flanking the integration site. In addition, our analysis identified several other non-B DNA motifs as new factors that potentially influence lentivirus integration, including human immunodeficiency virus type-1 in human cells. Together, our data demonstrate a clinically favorable integration site profile in the murine brain and identify non-B DNA as a potential new host factor that influences lentiviral integration in murine and human cells.

Selection and characterization of a novel DNA aptamer for label-free fluorescence biosensing of ochratoxin A

Nucleic acid aptamers are emerging as useful molecular recognition tools for food safety monitoring. However, practical and technical challenges limit the number and diversity of available aptamer probes that can be incorporated into novel sensing schemes. This work describes the selection of novel DNA aptamers that bind to the important food contaminant ochratoxin A (OTA). Following 15 rounds of in vitro selection, sequences were analyzed for OTA binding. Two of the isolated aptamers demonstrated high affinity binding and selectivity to this mycotoxin compared to similar food adulterants. These sequences, as well as a truncated aptamer (minimal sequence required for binding), were incorporated into a SYBR® Green I fluorescence-based OTA biosensing scheme. This label-free detection platform is capable of rapid, selective, and sensitive OTA quantification with a limit of detection of 9 nM and linear quantification up to 100 nM.

Association of G-quadruplex forming sequences with human mtDNA deletion breakpoints

BACKGROUND

Mitochondrial DNA (mtDNA) deletions cause disease and accumulate during aging, yet our understanding of the molecular mechanisms underlying their formation remains rudimentary. Guanine-quadruplex (GQ) DNA structures are associated with nuclear DNA instability in cancer; recent evidence indicates they can also form in mitochondrial nucleic acids, suggesting that these non-B DNA structures could be associated with mtDNA deletions. Currently, the multiple types of GQ sequences and their association with human mtDNA stability are unknown.

RESULTS

Here, we show an association between human mtDNA deletion breakpoint locations (sites where DNA ends rejoin after deletion of a section) and sequences with G-quadruplex forming potential (QFP), and establish the ability of selected sequences to form GQ in vitro. QFP contain four runs of either two or three consecutive guanines (2G and 3G, respectively), and we identified four types of QFP for subsequent analysis: intrastrand 2G, intrastrand 3G, duplex derived interstrand (ddi) 2G, and ddi 3G QFP sequences. We analyzed the position of each motif set relative to either 5' or 3' unique mtDNA deletion breakpoints, and found that intrastrand QFP sequences, but not ddi QFP sequences, showed significant association with mtDNA deletion breakpoint locations. Moreover, a large proportion of these QFP sequences occur at smaller distances to breakpoints relative to distribution-matched controls. The positive association of 2G QFP sequences persisted when breakpoints were divided into clinical subgroups. We tested in vitro GQ formation of representative mtDNA sequences containing these 2G QFP sequences and detected robust GQ structures by UV-VIS and CD spectroscopy. Notably, the most frequent deletion breakpoints, including those of the "common deletion", are bounded by 2G QFP sequence motifs.

CONCLUSIONS

The potential for GQ to influence mitochondrial genome stability supports a high-priority investigation of these structures and their regulation in normal and pathological mitochondrial biology. These findings emphasize the potential importance of helicases that subsequently resolve GQ to maintain the stability of the mitochondrial genome.

Relevance of G-quadruplex structures to pharmacogenetics

G-quadruplexes are non-canonical secondary structures formed within nucleic acids that are involved in modulating cellular processes such as replication, gene regulation, recombination and epigenetics. Within genes, there is mounting evidence of G-quadruplex involvement in transcriptional and post transcriptional regulation. We report the presence of potential G-quadruplex motifs within relevant sites of some important pharmacogenes and discuss the possible implications of this on the function and expression of these genes. Appreciating the location and potential functions of these motifs may be of value when considering the impacts of some pharmacogenetic variants. G-quadruplexes are also the focus of drug development efforts in oncology and we highlight the broader pharmacological implications of treatment strategies that may target G-quadruplexes.

Selection of DNA aptamers against uropathogenic Escherichia coli NSM59 by quantitative PCR controlled Cell-SELEX

In order to better control nosocomial infections, and facilitate the most prudent and effective use of antibiotics, improved strategies for the rapid detection and identification of problematic bacterial pathogens are required. DNA aptamers have much potential in the development of diagnostic assays and biosensors to address this important healthcare need, but further development of aptamers targeting common pathogens, and the strategies used to obtain specific aptamers are required. Here we demonstrate the application of a quantitative PCR (qPCR) controlled Cell-SELEX process, coupled with downstream secondary-conformation-based aptamer profiling. We used this approach to identify and select DNA aptamers targeted against uropathogenic Escherichia coli, for which specific aptamers are currently lacking, despite the prevalence of these infections. The use of qPCR to monitor the Cell-SELEX process permitted a minimal number of SELEX cycles to be employed, as well as the cycle-by-cycle optimisation of standard PCR amplification of recovered aptamer pools at each round. Identification of useful aptamer candidates was also facilitated by profiling of secondary conformations and selection based on putative aptamer secondary structure. One aptamer selected this way (designated EcA5-27), displaying a guanine-quadruplex sequence motif, was shown to have high affinity and specificity for target cells, and the potential to discriminate between distinct strains of E. coli, highlighting the possibility for development of aptamers selectively recognising pathogenic strains. Overall, the identified aptamers hold much potential for the development of rapid diagnostic assays for nosocomial urinary tract infections caused by E. coli.

Impact of alternative DNA structures on DNA damage, DNA repair, and genetic instability

Repetitive genomic sequences can adopt a number of alternative DNA structures that differ from the canonical B-form duplex (i.e. non-B DNA). These non-B DNA-forming sequences have been shown to have many important biological functions related to DNA metabolic processes; for example, they may have regulatory roles in DNA transcription and replication. In addition to these regulatory functions, non-B DNA can stimulate genetic instability in the presence or absence of DNA damage, via replication-dependent and/or replication-independent pathways. This review focuses on the interactions of non-B DNA conformations with DNA repair proteins and how these interactions impact genetic instability.

Optimal antisense target reducing INS intron 1 retention is adjacent to a parallel G quadruplex

Splice-switching oligonucleotides (SSOs) have been widely used to inhibit exon usage but antisense strategies that promote removal of entire introns to increase splicing-mediated gene expression have not been developed. Here we show reduction of INS intron 1 retention by SSOs that bind transcripts derived from a human haplotype expressing low levels of proinsulin. This haplotype is tagged by a polypyrimidine tract variant rs689 that decreases the efficiency of intron 1 splicing and increases the relative abundance of mRNAs with extended 5' untranslated region (5' UTR), which curtails translation. Co-expression of haplotype-specific reporter constructs with SSOs bound to splicing regulatory motifs and decoy splice sites in primary transcripts revealed a motif that significantly reduced intron 1-containing mRNAs. Using an antisense microwalk at a single nucleotide resolution, the optimal target was mapped to a splicing silencer containing two pseudoacceptor sites sandwiched between predicted RNA guanine (G) quadruplex structures. Circular dichroism spectroscopy and nuclear magnetic resonance of synthetic G-rich oligoribonucleotide tracts derived from this region showed formation of a stable parallel 2-quartet G-quadruplex on the 3' side of the antisense retention target and an equilibrium between quadruplexes and stable hairpin-loop structures bound by optimal SSOs. This region interacts with heterogeneous nuclear ribonucleoproteins F and H that may interfere with conformational transitions involving the antisense target. The SSO-assisted promotion of weak intron removal from the 5' UTR through competing noncanonical and canonical RNA structures may facilitate development of novel strategies to enhance gene expression.

Analysis of FMRP mRNA target datasets reveals highly associated mRNAs mediated by G-quadruplex structures formed via clustered WGGA sequences

Fragile X syndrome, a common cause of intellectual disability and a well-known cause of autism spectrum disorder, is the result of loss or dysfunction of fragile X mental retardation protein (FMRP), a highly selective RNA-binding protein and translation regulator. A major research priority has been the identification of the mRNA targets of FMRP, particularly as recent studies suggest an excess of FMRP targets among genes implicated in idiopathic autism and schizophrenia. Several large-scale studies have attempted to identify mRNAs bound by FMRP through several methods, each generating a list of putative target genes, leading to distinct hypotheses by which FMRP recognizes its targets; namely, by RNA structure or sequence. However, no in depth analyses have been performed to identify the level of consensus among the studies. Here, we analyze four large FMRP target datasets to generate high-confidence consensus lists, and examine all datasets for sequence elements within the target RNAs to validate reported FMRP binding motifs (GACR, ACUK and WGGA). We found GACR to be highly enriched in FMRP datasets, while ACUK was not. The WGGA pattern was modestly enriched in several, but not all datasets. The previous association between FMRP and G-quadruplexes prompted the analysis of the distribution of WGGA in the target genes. Consistent with the requirements for G-quadruplex formation, we observed highly clustered WGGA motifs in FMRP targets compared with other genes, implicating both RNA structure and sequence in the recognition motif of FMRP. In addition, we generate a list of the top 40 FMRP targets associated with FXS-related phenotypes.

Divergent LIN28-mRNA associations result in translational suppression upon the initiation of differentiation

LIN28 function is fundamental to the activity and behavior of human embryonic stem cells (hESCs) and induced pluripotent stem cells. Its main roles in these cell types are the regulation of translational efficiency and let-7 miRNA maturation. However, LIN28-associated mRNA cargo shifting and resultant regulation of translational efficiency upon the initiation of differentiation remain unknown. An RNA-immunoprecipitation and microarray analysis protocol, eRIP, that has high specificity and sensitivity was developed to test endogenous LIN28-associated mRNA cargo shifting. A combined eRIP and polysome analysis of early stage differentiation of hESCs with two distinct differentiation cues revealed close similarities between the dynamics of LIN28 association and translational modulation of genes involved in the Wnt signaling, cell cycle, RNA metabolism and proteasomal pathways. Our data demonstrate that change in translational efficiency is a major contributor to early stages of differentiation of hESCs, in which LIN28 plays a central role. This implies that eRIP analysis of LIN28-associated RNA cargoes may be used for rapid functional quality control of pluripotent stem cells under manufacture for therapeutic applications.

Palindrome-mediated and replication-dependent pathogenic structural rearrangements within the NF1 gene

Palindromic sequences can form hairpin structures or cruciform extrusions, which render them susceptible to genomic rearrangements. A 197-bp long palindromic AT-rich repeat (PATRR17) is located within intron 40 of the neurofibromatosis type 1 (NF1) gene (17q11.2). Through comprehensive NF1 analysis, we identified six unrelated patients with a rearrangement involving intron 40 (five deletions and one reciprocal translocation t(14;17)(q32;q11.2)). We hypothesized that PATRR17 may be involved in these rearrangements thereby causing NF1. Breakpoint cloning revealed that PATRR17 was indeed involved in all of the rearrangements. As microhomology was present at all breakpoint junctions of the deletions identified, and PATRR17 partner breakpoints were located within 7.1 kb upstream of PATRR17, fork stalling and template switching/microhomology-mediated break-induced replication was the most likely rearrangement mechanism. For the reciprocal translocation case, a 51 bp insertion at the translocation breakpoints mapped to a short sequence within PATRR17, proximal to the breakpoint, suggesting a multiple stalling and rereplication process, in contrast to previous studies indicating a purely replication-independent mechanism for PATRR-mediated translocations. In conclusion, we show evidence that PATRR17 is a hotspot for pathogenic intragenic deletions within the NF1 gene and suggest a novel replication-dependent mechanism for PATRR-mediated translocation.

QGRS-Conserve: a computational method for discovering evolutionarily conserved G-quadruplex motifs

BACKGROUND

Nucleic acids containing guanine tracts can form quadruplex structures via non-Watson-Crick base pairing. Formation of G-quadruplexes is associated with the regulation of important biological functions such as transcription, genetic instability, DNA repair, DNA replication, epigenetic mechanisms, regulation of translation, and alternative splicing. G-quadruplexes play important roles in human diseases and are being considered as targets for a variety of therapies. Identification of functional G-quadruplexes and the study of their overall distribution in genomes and transcriptomes is an important pursuit. Traditional computational methods map sequence motifs capable of forming G-quadruplexes but have difficulty in distinguishing motifs that occur by chance from ones which fold into G-quadruplexes.

RESULTS

We present Quadruplex forming 'G'-rich sequences (QGRS)-Conserve, a computational method for calculating motif conservation across exomes and supports filtering to provide researchers with more precise methods of studying G-quadruplex distribution patterns. Our method quantitatively evaluates conservation between quadruplexes found in homologous nucleotide sequences based on several motif structural characteristics. QGRS-Conserve also efficiently manages overlapping G-quadruplex sequences such that the resulting datasets can be analyzed effectively.

CONCLUSIONS

We have applied QGRS-Conserve to identify a large number of G-quadruplex motifs in the human exome conserved across several mammalian and non-mammalian species. We have successfully identified multiple homologs of many previously published G-quadruplexes that play post-transcriptional regulatory roles in human genes. Preliminary large-scale analysis identified many homologous G-quadruplexes in the 5'- and 3'-untranslated regions of mammalian species. An expectedly smaller set of G-quadruplex motifs was found to be conserved across larger phylogenetic distances. QGRS-Conserve provides means to build datasets that can be filtered and categorized in a variety of biological dimensions for more targeted studies in order to better understand the roles that G-quadruplexes play.

Fragile X mental retardation protein regulates translation by binding directly to the ribosome

Fragile X syndrome (FXS) is the most common form of inherited mental retardation, and it is caused by loss of function of the fragile X mental retardation protein (FMRP). FMRP is an RNA-binding protein that is involved in the translational regulation of several neuronal mRNAs. However, the precise mechanism of translational inhibition by FMRP is unknown. Here, we show that FMRP inhibits translation by binding directly to the L5 protein on the 80S ribosome. Furthermore, cryoelectron microscopic reconstruction of the 80S ribosome⋅FMRP complex shows that FMRP binds within the intersubunit space of the ribosome such that it would preclude the binding of tRNA and translation elongation factors on the ribosome. These findings suggest that FMRP inhibits translation by blocking the essential components of the translational machinery from binding to the ribosome.

U3 region in the HIV-1 genome adopts a G-quadruplex structure in its RNA and DNA sequence

Genomic regions rich in G residues are prone to adopt G-quadruplex structure. Multiple Sp1-binding motifs arranged in tandem have been suggested to form this structure in promoters of cancer-related genes. Here, we demonstrate that the G-rich proviral DNA sequence of the HIV-1 U3 region, which serves as a promoter of viral transcription, adopts a G-quadruplex structure. The sequence contains three binding elements for transcription factor Sp1, which is involved in the regulation of HIV-1 latency, reactivation, and high-level virus expression. We show that the three Sp1 binding motifs can adopt different forms of G-quadruplex structure and that the Sp1 protein can recognize and bind to its site folded into a G-quadruplex. In addition, a c-kit2 specific antibody, designated hf2, binds to two different G-quadruplexes formed in Sp1 sites. Since U3 is encoded at both viral genomic ends, the G-rich sequence is also present in the RNA genome. We demonstrate that the RNA sequence of U3 forms dimers with characteristics known for intermolecular G-quadruplexes. Together with previous reports showing G-quadruplex dimers in the gag and cPPT regions, these results suggest that integrity of the two viral genomes is maintained through numerous intermolecular G-quadruplexes formed in different RNA genome locations. Reconstituted reverse transcription shows that the potassium-dependent structure formed in U3 RNA facilitates RT template switching, suggesting that the G-quadruplex contributes to recombination in U3.

A bacterial DNA quadruplex with exceptional K+ selectivity and unique structural polymorphism

The G-quadruplex forming sequence d[(G₄CT)₃G₄] shows complete and continuous quadruplex interconversion upon increasing K⁺-concentrations and pronounced K⁺ selectivity.

The RNA helicase RHAU (DHX36) suppresses expression of the transcription factor PITX1

RNA Helicase associated with AU-rich element (RHAU) (DHX36) is a DEAH (Aspartic acid, Glumatic Acid, Alanine, Histidine)-box RNA helicase that can bind and unwind G4-quadruplexes in DNA and RNA. To detect novel RNA targets of RHAU, we performed an RNA co-immunoprecipitation screen and identified the PITX1 messenger RNA (mRNA) as specifically and highly enriched. PITX1 is a homeobox transcription factor with roles in both development and cancer. Primary sequence analysis identified three probable quadruplexes within the 3′-untranslated region of the PITX1 mRNA. Each of these sequences, when isolated, forms stable quadruplex structures that interact with RHAU. We provide evidence that these quadruplexes exist in the endogenous mRNA; however, we discovered that RHAU is tethered to the mRNA via an alternative non–quadruplex-forming region. RHAU knockdown by small interfering RNA results in significant increases in PITX1 protein levels with only marginal changes in mRNA, suggesting a role for RHAU in translational regulation. Involvement of components of the microRNA machinery is supported by similar and non-additive increases in PITX1 protein expression on Dicer and combined RHAU/Dicer knockdown. We also demonstrate a requirement of argonaute-2, a key RNA-induced silencing complex component, to mediate RHAU-dependent changes in PITX1 protein levels. These results demonstrate a novel role for RHAU in microRNA-mediated translational regulation at a quadruplex-containing 3′-untranslated region.