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

Recombination events among virulence genes in malaria parasites are associated with G-quadruplex-forming DNA motifs

Background

Malaria parasites of the genus Plasmodium possess large hyper-variable families of antigen-encoding genes. These are often variantly-expressed and are major virulence factors for immune evasion and the maintenance of chronic infections. Recombination and diversification of these gene families occurs readily, and may be promoted by G-quadruplex (G4) DNA motifs within and close to the variant genes. G4s have been shown to cause replication fork stalling, DNA breakage and recombination in model systems, but these motifs remain largely unstudied in Plasmodium.

Results

We examined the nature and distribution of putative G4-forming sequences in multiple Plasmodium genomes, finding that their co-distribution with variant gene families is conserved across different Plasmodium species that have different types of variant gene families. In P. falciparum, where a large set of recombination events that occurred over time in cultured parasites has been mapped, we found a strong spatial association between these recombination events and putative G4-forming sequences. Finally, we searched Plasmodium genomes for the three classes of helicase that can unwind G4s: Plasmodium spp. have no identifiable homologue of the highly efficient G4 helicase PIF1, but they do encode two putative RecQ helicases and one homologue of the RAD3-family helicase FANCJ.

Conclusions

Our analyses, conducted at the whole-genome level in multiple species of Plasmodium, support the concept that G4s are likely to be involved in recombination and diversification of antigen-encoding gene families in this important protozoan pathogen.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3183-3) contains supplementary material, which is available to authorized users.

Development of a DNA aptamer that binds specifically to group A Streptococcus serotype M3

Group A streptococcus (GAS) is an important Gram-positive pathogen that causes various human diseases ranging from peripheral lesions to invasive infections. The M protein is one of the main virulence factors present on the cell surface and is associated with invasive GAS infections. Compared with other M types, serotype M3 has a predominant role in lethal infections and demonstrates epidemic behaviors, including streptococcal toxic shock syndrome, bacteremia, and necrotizing fasciitis. Traditional methods for M typing are time-consuming, tedious, contradictory, and generally restricted to reference laboratories. Therefore, development of a new M-typing technique is needed. Aptamers with the ability to detect their target with a high degree of accuracy and specificity can be ideal candidates for specific M-typing of Streptococcus pyogenes. In this study DNA aptamers with a high binding affinity towards S. pyogenes serotype M3 were selected through 12 iterative rounds of the Systematic Evolution of Ligands by EXponential (SELEX) enrichment procedure using live cells as a target. We monitored the progress of the SELEX procedure by flow cytometry analysis. Of several aptamer sequences analyzed, 12L18A showed the highest binding efficiency towards S. pyogenes type M3, with an apparent dissociation constant (K_d) of 7.47 ± 1.72 pmol/L being the lowest. Therefore the isolated aptamer can be used in any tool, such as a biosensor, for the detection of S. pyogenes and can be used in the development of a novel M-typing system.

Zika Virus Genomic RNA Possesses Conserved G-Quadruplexes Characteristic of the Flaviviridae Family

Zika virus has emerged as a global concern because neither a vaccine nor antiviral compounds targeting it exist. A structure for the positive-sense RNA genome has not been established, leading us to look for potential G-quadruplex sequences (PQS) in the genome. The analysis identified >60 PQSs in the Zika genome. To minimize the PQS population, conserved sequences in the Flaviviridae family were found by sequence alignment, identifying seven PQSs in the prM, E, NS1, NS3, and NS5 genes. Next, alignment of 78 Zika strain genomes identified a unique PQS near the end of the 3'-UTR. Structural studies on the G-quadruplex sequences found four of the conserved Zika virus sequences to adopt stable, parallel-stranded folds that bind a G-quadruplex-specific compound, and one that was studied caused polymerase stalling when folded to a G-quadruplex. Targeting these PQSs with G-quadruplex binding molecules validated in previous clinical trials may represent a new approach for inhibiting viral replication.

Laboratory evolution of artificially expanded DNA gives redesignable aptamers that target the toxic form of anthrax protective antigen

Reported here is a laboratory in vitro evolution (LIVE) experiment based on an artificially expanded genetic information system (AEGIS). This experiment delivers the first example of an AEGIS aptamer that binds to an isolated protein target, the first whose structural contact with its target has been outlined and the first to inhibit biologically important activities of its target, the protective antigen from Bacillus anthracis. We show how rational design based on secondary structure predictions can also direct the use of AEGIS to improve the stability and binding of the aptamer to its target. The final aptamer has a dissociation constant of ∼35 nM. These results illustrate the value of AEGIS-LIVE for those seeking to obtain receptors and ligands without the complexities of medicinal chemistry, and also challenge the biophysical community to develop new tools to analyze the spectroscopic signatures of new DNA folds that will emerge in synthetic genetic systems replacing standard DNA and RNA as platforms for LIVE.

Effects of metal ions and cosolutes on G-quadruplex topology

Topologies of G-quadruplexes depend on oligonucleotide sequences and on environmental factors, and the diversity of G-quadruplex topologies complicates investigation of functions of these nucleic acid structures. To investigate how metal ions and cosolutes regulate topologies of G-quadruplexes, we stabilized the antiparallel conformation by insertion of 2'-deoxyxanthosine and 8-oxo-2'-deoxyguanosine into selected positions of an oligonucleotide. Thermodynamic analyses of the oligonucleotide revealed that Na⁺ stabilized the antiparallel G-quadruplex, whereas K⁺ destabilized this topology. This result suggests that metal ions selectively stabilize G-quadruplex topologies with cavities between G-quartet planes of certain sizes. In the presence of KCl in 20wt% poly(ethylene glycol) with average molecular weight of 200, the antiparallel basket-type G-quadruplex conformation was not stabilized compared with the dilute condition. In the presence of NaCl, the cosolute did stabilize the G-quadruplex with respect to the dilute condition. The presented data show that metal ions and cosolutes regulate topologies of G-quadruplexes through mechanisms that depend on sizes of metal ion cavities and hydration states.

A G-Rich Motif in the lncRNA Braveheart Interacts with a Zinc-Finger Transcription Factor to Specify the Cardiovascular Lineage

Long non-coding RNAs (lncRNAs) are an emerging class of transcripts that can modulate gene expression; however, their mechanisms of action remain poorly understood. Here, we experimentally determine the secondary structure of Braveheart (Bvht) using chemical probing methods and show that this ∼590 nt transcript has a modular fold. Using CRISPR/Cas9-mediated editing of mouse embryonic stem cells, we find that deletion of 11 nt in a 5' asymmetric G-rich internal loop (AGIL) of Bvht (bvht^dAGIL) dramatically impairs cardiomyocyte differentiation. We demonstrate a specific interaction between AGIL and cellular nucleic acid binding protein (CNBP/ZNF9), a zinc-finger protein known to bind single-stranded G-rich sequences. We further show that CNBP deletion partially rescues the bvht^dAGIL mutant phenotype by restoring differentiation capacity. Together, our work shows that Bvht functions with CNBP through a well-defined RNA motif to regulate cardiovascular lineage commitment, opening the door for exploring broader roles of RNA structure in development and disease.

Palindrome analyser - A new web-based server for predicting and evaluating inverted repeats in nucleotide sequences

DNA cruciform structures play an important role in the regulation of natural processes including gene replication and expression, as well as nucleosome structure and recombination. They have also been implicated in the evolution and development of diseases such as cancer and neurodegenerative disorders. Cruciform structures are formed by inverted repeats, and their stability is enhanced by DNA supercoiling and protein binding. They have received broad attention because of their important roles in biology. Computational approaches to study inverted repeats have allowed detailed analysis of genomes. However, currently there are no easily accessible and user-friendly tools that can analyse inverted repeats, especially among long nucleotide sequences. We have developed a web-based server, Palindrome analyser, which is a user-friendly application for analysing inverted repeats in various DNA (or RNA) sequences including genome sequences and oligonucleotides. It allows users to search and retrieve desired gene/nucleotide sequence entries from the NCBI databases, and provides data on length, sequence, locations and energy required for cruciform formation. Palindrome analyser also features an interactive graphical data representation of the distribution of the inverted repeats, with options for sorting according to the length of inverted repeat, length of loop, and number of mismatches. Palindrome analyser can be accessed at http://bioinformatics.ibp.cz.

DNA G-quadruplexes show strong interaction with DNA methyltransferases in vitro

The DNA methyltransferase enzymes (DNMTs) catalyzing cytosine methylation do so at specific locations of the genome, although with some level of redundancy. The de novo methyltransferases DNMT3A and 3B play a vital role in methylating the genome of the developing embryo in regions devoid of methylation marks. The ability of DNMTs to colocalize at sites of DNA damage is suggestive that recognition of mispaired bases and unusual structures is inherent to the function of these proteins. We provide evidence for G-quadruplex formation within imprinted gene promoters, and report high-affinity binding of recombinant human DNMTs to such DNA G-quadruplexes in vitro. These observations suggest a potential interaction of G-quadruplexes with the DNA methylation machinery, which may be of epigenetic and biological significance.

RNA-Binding Protein FXR1 Regulates p21 and TERC RNA to Bypass p53-Mediated Cellular Senescence in OSCC

RNA-binding proteins (RBP) regulate numerous aspects of co- and post-transcriptional gene expression in cancer cells. Here, we demonstrate that RBP, fragile X-related protein 1 (FXR1), plays an essential role in cellular senescence by utilizing mRNA turnover pathway. We report that overexpressed FXR1 in head and neck squamous cell carcinoma targets (G-quadruplex (G4) RNA structure within) both mRNA encoding p21 (Cyclin-Dependent Kinase Inhibitor 1A (CDKN1A, Cip1) and the non-coding RNA Telomerase RNA Component (TERC), and regulates their turnover to avoid senescence. Silencing of FXR1 in cancer cells triggers the activation of Cyclin-Dependent Kinase Inhibitors, p53, increases DNA damage, and ultimately, cellular senescence. Overexpressed FXR1 binds and destabilizes p21 mRNA, subsequently reduces p21 protein expression in oral cancer cells. In addition, FXR1 also binds and stabilizes TERC RNA and suppresses the cellular senescence possibly through telomerase activity. Finally, we report that FXR1-regulated senescence is irreversible and FXR1-depleted cells fail to form colonies to re-enter cellular proliferation. Collectively, FXR1 displays a novel mechanism of controlling the expression of p21 through p53-dependent manner to bypass cellular senescence in oral cancer cells.

Multifaceted analyses of the interactions between human heart type fatty acid binding protein and its specific aptamers

BACKGROUND

Aptamer-protein interaction studies have been mainly confined to dissociation constant (K_d) determination. A combinatorial approach involving limited proteolysis mass spectroscopy, molecular docking and CD studies is reported here to elucidate the specific interactions involved.

METHODS

To generate aptamers specific for human FABP3, SELEX was performed incorporating counter SELEX cycles against control FABPs and GST tag, followed by their characterization by EMSA, CD and SVD analysis. Based on computationally obtained aptamer-protein complex models, the interacting aptamer, and protein residues were predicted and supported by limited proteolysis experiments.

RESULTS

Two aptamers N13 and N53 specific for human fatty acid binding protein (FABP3) were isolated with corresponding K_d of 0.0743±0.0142μM and 0.3337±0.1485μM for FABP3 interactions. Both aptamers possess stable B-DNA structures at salt concentration of 100mM and pH range (6-9). The N13 aptamer led interaction involved 3 salt bridges and 2 hydrogen bonds, whereas N53 had 2 salt bridges with 8 hydrogen and 7 hydrophobic interactions.

CONCLUSIONS

The aptamers generated are the first to be reported against human FABP3. The higher interaction footprint of N53 incited synergistic conformational changes in both N53 and FABP3 during interaction, leading to a decline in binding affinity in comparison to N13 which corroborated to the calculated K_d values.

GENERAL SIGNIFICANCE

This combinatorial method may be used to retrieve the possible specific binding modes and interaction patterns involved in large aptamer-protein complexes. Thus the method can be exploited to identify the optimum aptamer length for in-depth structure-function studies and its tailored applications.

Detection of cis- and trans-acting Factors in DNA Structure-Induced Genetic Instability Using In silico and Cellular Approaches

Sequences that can adopt alternative DNA structures (i.e., non-B DNA) are very abundant in mammalian genomes, and recent studies have revealed many important biological functions of non-B DNA structures in chromatin remodeling, DNA replication, transcription, and genetic instability. Here, we provide results from an in silico web-based search engine coupled with cell-based experiments to characterize the roles of non-B DNA conformations in genetic instability in eukaryotes. The purpose of this article is to illustrate strategies that can be used to identify and interrogate the biological roles of non-B DNA structures, particularly on genetic instability. We have included unpublished data using a short H-DNA-forming sequence from the human c-MYC promoter region as an example, and identified two different mechanisms of H-DNA-induced genetic instability in yeast and mammalian cells: a DNA replication-related model of mutagenesis; and a replication-independent cleavage model. Further, we identified candidate proteins involved in H-DNA-induced genetic instability by using a yeast genetic screen. A combination of in silico and cellular methods, as described here, should provide further insight into the contributions of non-B DNA structures in biological functions, genetic evolution, and disease development.

Multiple G-quartet structures in pre-edited mRNAs suggest evolutionary driving force for RNA editing in trypanosomes

Mitochondrial transcript maturation in African trypanosomes requires a U-nucleotide specific RNA editing reaction. In its most extreme form hundreds of U's are inserted into and deleted from primary transcripts to generate functional mRNAs. Unfortunately, both origin and biological role of the process have remained enigmatic. Here we report a so far unrecognized structural feature of pre-edited mRNAs. We demonstrate that the cryptic pre-mRNAs contain numerous clustered G-nt, which fold into G-quadruplex (GQ) structures. We identified 27 GQ's in the different pre-mRNAs and demonstrate a positive correlation between the steady state abundance of guide (g)RNAs and the sequence position of GQ-elements. We postulate that the driving force for selecting G-rich sequences lies in the formation of DNA/RNA hybrid G-quadruplex (HQ) structures between the pre-edited transcripts and the non-template strands of mitochondrial DNA. HQ's are transcription termination/replication initiation sites and thus guarantee an unperturbed replication of the mt-genome. This is of special importance in the insect-stage of the parasite. In the transcription-on state, the identified GQ's require editing as a GQ-resolving activity indicating a link between replication, transcription and RNA editing. We propose that the different processes have coevolved and suggest the parasite life-cycle and the single mitochondrion as evolutionary driving forces.

RAID3--An interleukin-6 receptor-binding aptamer with post-selective modification-resistant affinity

Aptamers are an emerging class of highly specific targeting ligands. They can be selected in vitro for a large variety of targets, ranging from small molecules to whole cells. Most aptamers selected are nucleic acid-based, allowing chemical synthesis and easy modification. Although their properties make them interesting drug candidates for a broad spectrum of applications and an interesting alternative to antibodies or fusion proteins, they are not yet broadly used. One major drawback of aptamers is their susceptibility to abundant serum nucleases, resulting in their fast degradation in biological fluids. Using modified nucleic acids has become a common strategy to overcome these disadvantages, greatly increasing their half-life under cell culture conditions or even in vivo. Whereas pre-selective modifications of the initial library for aptamer selection are relatively easy to obtain, post-selective modifications of already selected aptamers are still generally very labor-intensive and often compromise the aptamers ability to bind its target molecule. Here we report the selection, characterization and post-selective modification of a 34 nucleotide (nt) RNA aptamer for a non-dominant, novel target site (domain 3) of the interleukin-6 receptor (IL-6R). We performed structural analyses and investigated the affinity of the aptamer to the membrane-bound and soluble forms (sIL-6R) of the IL-6R. Further, we performed structural analyses of the aptamer in solution using small-angle X-ray scattering and determined its overall shape and oligomeric state. Post-selective exchange of all pyrimidines against their 2'-fluoro analogs increased the aptamers stability significantly without compromising its affinity for the target protein. The resulting modified aptamer could be shortened to its minimal binding motif without loss of affinity.

Genome-wide discovery of G-quadruplex forming sequences and their functional relevance in plants

DNA, in addition to the canonical B-form, can acquire a variety of alternate structures, such as G-quadruplexes. These structures have been implicated in several cellular processes in animals. In this study, we identified different types of G-quadruplex forming sequences (GQSes) in 15 sequenced plants and analyzed their distribution in various genomic features, including gene body, coding, intergenic and promoter regions. G2-type GQSes were most abundant in all the plant species analyzed. A strong association of G3-type GQSes with intergenic, promoter and intronic regions was found. However, G2-type GQSes were enriched in genic, CDS, exonic and untranslated regions. Further, we identified GQSes present in the conserved genes among monocots and dicots. The genes involved in development, cell growth and size, transmembrane transporter, and regulation of gene expression were found to be significantly enriched. In the promoter region, we detected strong co-occurrence of Telobox, ERF, MYB, RAV1B and E2F motifs with GQSes. Further, we validated the structure formation of several plant GQSes, demonstrated their effect on stalling in-vitro replication and revealed their interaction with plant nuclear proteins. Our data provide insights into the prevalence of GQSes in plants, establish their association with different genomic features and functional relevance.

IFI16 Preferentially Binds to DNA with Quadruplex Structure and Enhances DNA Quadruplex Formation

Interferon-inducible protein 16 (IFI16) is a member of the HIN-200 protein family, containing two HIN domains and one PYRIN domain. IFI16 acts as a sensor of viral and bacterial DNA and is important for innate immune responses. IFI16 binds DNA and binding has been described to be DNA length-dependent, but a preference for supercoiled DNA has also been demonstrated. Here we report a specific preference of IFI16 for binding to quadruplex DNA compared to other DNA structures. IFI16 binds to quadruplex DNA with significantly higher affinity than to the same sequence in double stranded DNA. By circular dichroism (CD) spectroscopy we also demonstrated the ability of IFI16 to stabilize quadruplex structures with quadruplex-forming oligonucleotides derived from human telomere (HTEL) sequences and the MYC promotor. A novel H/D exchange mass spectrometry approach was developed to assess protein interactions with quadruplex DNA. Quadruplex DNA changed the IFI16 deuteration profile in parts of the PYRIN domain (aa 0–80) and in structurally identical parts of both HIN domains (aa 271–302 and aa 586–617) compared to single stranded or double stranded DNAs, supporting the preferential affinity of IFI16 for structured DNA. Our results reveal the importance of quadruplex DNA structure in IFI16 binding and improve our understanding of how IFI16 senses DNA. IFI16 selectivity for quadruplex structure provides a mechanistic framework for IFI16 in immunity and cellular processes including DNA damage responses and cell proliferation.

Apta-nanosensor preparation and in vitro assay for rapid Diazinon detection using a computational molecular approach

Aptamers (ss-DNA or ss-RNA), also known as artificial antibodies, have been selected in vitro median to bind target molecules with high affinity and selectivity. Diazinon is one of the most widely used organophosphorus insecticides in developing and underdeveloped countries as insecticide and acaricide. Diazinon is readily absorbed from the gastrointestinal system and rapidly distributed throughout the body. Thus, the design of clinical and laboratory diagnostics using nanobiosensors is necessary. A computational approach allows us to screen or rank receptor structure and predict interaction outcomes with a deeper understanding, and it is much more cost effective than laboratory attempts. In this research, the best sequence (high affinity bind Diazinon-ssDNA) was ranked among 12 aptamers isolated from SELEX experimentation. Docking results, as the first virtual screening stage and static technique, selected frequent conformation of each aptamer. Then, the quantity and quality of aptamer-Diazinon interaction were simulated using molecular dynamics as a mobility technique. RMSD, RMSF, radius of gyration, and the number of hydrogen bonds formed between Diazinon-aptamer were monitored to assess the quantity and quality of interactions. G-quadruplex DNA aptamer (DF20) showed to be a reliable candidate for Diazinon biosensing. The apta-nanosensor designed using simulation results allowed with linearity detection in the range of .141-.65 nM and a LOD of 17.903 nM, and it was validated using a computational molecular approach.

Improved Aptamers for the Diagnosis and Potential Treatment of HER2-Positive Cancer

Aptamers provide a potential source of alternative targeting molecules for existing antibody diagnostics and therapeutics. In this work, we selected novel DNA aptamers targeting the HER2 receptor by an adherent whole-cell SELEX approach. Individual aptamers were identified by next generation sequencing and bioinformatics analysis. Two aptamers, HeA2_1 and HeA2_3, were shown to bind the HER2 protein with affinities in the nanomolar range. In addition, both aptamers were able to bind with high specificity to HER2-overexpressing cells and HER2-positive tumor tissue samples. Furthermore, we demonstrated that aptamer HeA2_3 is being internalized into cancer cells and has an inhibitory effect on cancer cell growth and viability. In the end, we selected novel DNA aptamers with great potential for the diagnosis and possible treatment of HER2-positive cancer.

QuadBase2: web server for multiplexed guanine quadruplex mining and visualization

DNA guanine quadruplexes or G4s are non-canonical DNA secondary structures which affect genomic processes like replication, transcription and recombination. G4s are computationally identified by specific nucleotide motifs which are also called putative G4 (PG4) motifs. Despite the general relevance of these structures, there is currently no tool available that can allow batch queries and genome-wide analysis of these motifs in a user-friendly interface. QuadBase2 (quadbase.igib.res.in) presents a completely reinvented web server version of previously published QuadBase database. QuadBase2 enables users to mine PG4 motifs in up to 178 eukaryotes through the EuQuad module. This module interfaces with Ensembl Compara database, to allow users mine PG4 motifs in the orthologues of genes of interest across eukaryotes. PG4 motifs can be mined across genes and their promoter sequences in 1719 prokaryotes through ProQuad module. This module includes a feature that allows genome-wide mining of PG4 motifs and their visualization as circular histograms. TetraplexFinder, the module for mining PG4 motifs in user-provided sequences is now capable of handling up to 20 MB of data. QuadBase2 is a comprehensive PG4 motif mining tool that further expands the configurations and algorithms for mining PG4 motifs in a user-friendly way.

In silico analysis of regulatory and structural motifs of the ovine HSP90AA1 gene

Gene promoters are essential regions of DNA where the transcriptional molecular machinery to produce RNA molecules is recruited. In this process, DNA epigenetic modifications can acquire a fundamental role in the regulation of gene expression. Recently, in a previous work of our group, functional features and DNA methylation involved in the ovine HSP90AA1 gene expression regulation have been observed. In this work, we report a combination of methylation analysis by bisulfite sequencing in several tissues and at different developmental stages together with in silico bioinformatic analysis of putative regulating factors in order to identify regulative mechanisms both at the promoter and gene body. Our results show a "hybrid structure" (TATA box + CpG island) of the ovine HSP90AA1 gene promoter both in somatic and non-differentiated germ tissues, revealing the ability of the HSP90AA1 gene to be regulated both in an inducible and constitutive fashion. In addition, in silico analysis showed that several putative alternative spliced regulatory motifs, exonic splicing enhancers (ESEs), and G-quadruplex secondary structures were somehow related to the DNA methylation pattern found. The results obtained here could help explain the differences in cell-type transcripts, tissue expression rate, and transcription silencing mechanisms found in this gene.

Constitutive translation of human α-synuclein is mediated by the 5'-untranslated region

Genetic and biochemical studies have established a central role for α-synuclein (SNCA) accumulation in the pathogenesis of Parkinson's disease. Uncovering and subsequently interfering with physiological mechanisms that control SNCA expression is one approach to limit disease progression. To this end, the long and GC-rich 5'-untranslated region (UTR) of SNCA, which is predicted to fold into stable hairpin and G-quadruplex RNA motifs, was investigated for its role in mRNA translation. Inclusion of SNCA 5'-UTR significantly induced expression of both SNCA and luciferase ORF constructs. This effect was not associated with a change in mRNA levels or differential nucleocytoplasmic shuttling. Further, the presence of the 5'-UTR enhanced SNCA synthesis when cap-dependent translation was attenuated with rapamycin treatment. Analysis using multiple methodologies revealed that the 5'-UTR harbours an internal ribosome entry site (IRES) element that spans most of its nucleotide sequence. Signals such as plasma-membrane depolarization, serum starvation and oxidative stress stimulated SNCA protein translation via its 5'-UTR as well as enhanced its IRES activity. Taken together, these data support the idea that the 5'-UTR is an important positive regulator of SNCA synthesis under diverse physiological and pathological conditions, explaining in part the abundance of SNCA in healthy neurons and its accumulation in degenerative cells.

Quantitative analysis and prediction of G-quadruplex forming sequences in double-stranded DNA

G-quadruplex (GQ) is a four-stranded DNA structure that can be formed in guanine-rich sequences. GQ structures have been proposed to regulate diverse biological processes including transcription, replication, translation and telomere maintenance. Recent studies have demonstrated the existence of GQ DNA in live mammalian cells and a significant number of potential GQ forming sequences in the human genome. We present a systematic and quantitative analysis of GQ folding propensity on a large set of 438 GQ forming sequences in double-stranded DNA by integrating fluorescence measurement, single-molecule imaging and computational modeling. We find that short minimum loop length and the thymine base are two main factors that lead to high GQ folding propensity. Linear and Gaussian process regression models further validate that the GQ folding potential can be predicted with high accuracy based on the loop length distribution and the nucleotide content of the loop sequences. Our study provides important new parameters that can inform the evaluation and classification of putative GQ sequences in the human genome.

DNA Instability Maintains the Repeat Length of the Yeast RNA Polymerase II C-terminal Domain

The C-terminal domain (CTD) of RNA polymerase II in eukaryotes is comprised of tandemly repeating units of a conserved seven-amino acid sequence. The number of repeats is, however, quite variable across different organisms. Furthermore, previous studies have identified evidence of rearrangements within the CTD coding region, suggesting that DNA instability may play a role in regulating or maintaining CTD repeat number. The work described here establishes a clear connection between DNA instability and CTD repeat number in Saccharomyces cerevisiae First, analysis of 36 diverse S. cerevisiae isolates revealed evidence of numerous past rearrangements within the DNA sequence that encodes the CTD. Interestingly, the total number of CTD repeats was relatively static (24-26 repeats in all strains), suggesting a balancing act between repeat expansion and contraction. In an effort to explore the genetic plasticity within this region, we measured the rates of repeat expansion and contraction using novel reporters and a doxycycline-regulated expression system for RPB1 In efforts to determine the mechanisms leading to CTD repeat variability, we identified the presence of DNA secondary structures, specifically G-quadruplex-like DNA, within the CTD coding region. Furthermore, we demonstrated that mutating PIF1, a G-quadruplex-specific helicase, results in increased CTD repeat length polymorphisms. We also determined that RAD52 is necessary for CTD repeat expansion but not contraction, identifying a role for recombination in repeat expansion. Results from these DNA rearrangements may help explain the CTD copy number variation seen across eukaryotes, as well as support a model of CTD expansion and contraction to maintain CTD integrity and overall length.

TDP-43 binds and transports G-quadruplex-containing mRNAs into neurites for local translation

Growth and differentiation of the neurites depends on long-distance transport of a specific set of mRNAs to restricted area and their local translation. Here, we found that a TAR DNA-binding protein of 43 kDa in size (TDP-43) plays an essential role in intracellular transport of mRNA. For identification of target RNAs recognized by TDP-43, we purified TDP-43 in soluble dimer form and subjected to in vitro systematic evolution of ligands by exponential enrichment (SELEX) screening. All the TDP-43-bound RNAs were found to contain G-quadruplex (G4). Using a double-fluorescent probe system, G4-containing RNAs were found to be transported, together with TDP-43, into the distal neurites. Two lines of evidence indicated that loss of function of TDP-43 results in the neurodegenerative disorder: (i) amyotrophic lateral sclerosis (ALS)-linked mutant TDP-43M337V lacks the activity of binding and transport of G4-containing mRNAs; and (ii) RNA containing G4-forming GGGGCC repeat expansion from the ALS-linked C9orf72 gene absorbs TDP-43, thereby reducing the intracellular pool of functional TDP-43. Taken together, we propose that TDP-43 within neurons plays an essential role of mRNA transport into distal neurites for local translation, and thus, dysfunctions of TDP-43 cause neural diseases such as ALS and frontotemporal lobar degeneration.

Transcriptional Regulation of Atp-Dependent Chromatin Remodeling Factors: Smarcal1 and Brg1 Mutually Co-Regulate Each Other

The ATP-dependent chromatin remodeling factors regulate gene expression. However, it is not known whether these factors regulate each other. Given the ability of these factors to regulate the accessibility of DNA to transcription factors, we postulate that one ATP-dependent chromatin remodeling factor should be able to regulate the transcription of another ATP-dependent chromatin remodeling factor. In this paper, we show that BRG1 and SMARCAL1, both members of the ATP-dependent chromatin remodeling protein family, regulate each other. BRG1 binds to the SMARCAL1 promoter, while SMARCAL1 binds to the brg1 promoter. During DNA damage, the occupancy of SMARCAL1 on the brg1 promoter increases coinciding with an increase in BRG1 occupancy on the SMARCAL1 promoter, leading to increased brg1 and SMARCAL1 transcripts respectively. This is the first report of two ATP-dependent chromatin remodeling factors regulating each other.

Downregulation of the WT1 gene expression via TMPyP4 stabilization of promoter G-quadruplexes in leukemia cells

The WT1 gene is an important oncogene, and its overexpression is considered as an effective target for anticancer therapy. Regulation of its gene transcription is one way for WT1-targeting drug design. Recently, in silico analysis of some oncogene promoters like WT1 showed some guanine-rich regions with the ability to form G-quadruplex structures. Ligands like 5,10,15,20-tetra(N-methyl-4-pyridyl)-porphine (TMPyP4) have predominant effect on G-quadruplex stabilization. The aim of this study was to understand the effect of TMPyP4 on WT1 gene transcription via stabilization of promoter G-quadruplexes. We examined the formation of new G-quadruplex motifs in WT1 promoter in the presence of TMPyP4. In order to understand the nature of its interaction with WT1 promoter quadruplexes, differential pulse voltammetry (DPV), circular dichroism (CD), polyacrylamide gel electrophoresis, electrophoretic mobility shift assay (EMSA), polymerase chain reaction (PCR) stop assays, and quantitative RT-PCR were performed. According to the results, the WT1 promoter can form stable intramolecular parallel G-quadruplexes. In addition, after 48 and 96 h of incubation, 100 μM TMPyP4 reduced the WT1 transcription to 9 and 0.4 %, respectively, compare to control. We report that ligand-mediated stabilization of G-quadruplexes within the WT1 promoter can silence WT1 expression. This study might offer the basis for the reasonable design and improvement of new porphyrin derivatives as effective anti-leukemia agents for cancer therapy.

G-quadruplexes as novel cis-elements controlling transcription during embryonic development

G-quadruplexes are dynamic structures folded in G-rich single-stranded DNA regions. These structures have been recognized as a potential nucleic acid based mechanism for regulating multiple cellular processes such as replication, transcription and genomic maintenance. So far, their transcriptional role in vivo during vertebrate embryonic development has not yet been addressed. Here, we performed an in silico search to find conserved putative G-quadruplex sequences (PQSs) within proximal promoter regions of human, mouse and zebrafish developmental genes. Among the PQSs able to fold in vitro as G-quadruplex, those present in nog3, col2a1 and fzd5 promoters were selected for further studies. In cellulo studies revealed that the selected G-quadruplexes affected the transcription of luciferase controlled by the SV40 nonrelated promoter. G-quadruplex disruption in vivo by microinjection in zebrafish embryos of either small ligands or DNA oligonucleotides complementary to the selected PQSs resulted in lower transcription of the targeted genes. Moreover, zebrafish embryos and larvae phenotypes caused by the presence of complementary oligonucleotides fully resembled those ones reported for nog3, col2a1 and fzd5 morphants. To our knowledge, this is the first work revealing in vivo the role of conserved G-quadruplexes in the embryonic development, one of the most regulated processes of the vertebrates biology.

Epigenetics of the myotonic dystrophy-associated DMPK gene neighborhood

Aim:

Identify epigenetic marks in the vicinity of DMPK (linked to myotonic dystrophy, DM1) that help explain tissue-specific differences in its expression.

Materials & methods:

At DMPK and its flanking genes (DMWD, SIX5, BHMG1 and RSPH6A), we analyzed many epigenetic and transcription profiles from myoblasts, myotubes, skeletal muscle, heart and 30 nonmuscle samples.

Results:

In the DMPK gene neighborhood, muscle-associated DNA hypermethylation and hypomethylation, enhancer chromatin, and CTCF binding were seen. Myogenic DMPK hypermethylation correlated with high expression and decreased alternative promoter usage. Testis/sperm hypomethylation of BHMG1 and RSPH6A was associated with testis-specific expression. G-quadruplex (G4) motifs and sperm-specific hypomethylation were found near the DM1-linked CTG repeats within DMPK.

Conclusion:

Tissue-specific epigenetic features in DMPK and neighboring genes help regulate its expression. G4 motifs in DMPK DNA and RNA might contribute to DM1 pathology.

Characterization of MNK1b DNA Aptamers That Inhibit Proliferation in MDA-MB231 Breast Cancer Cells

Elevated expression levels of eukaryotic initiation factor 4E (eIF4E) promote cancer development and progression. MAP kinase interacting kinases (MNKs) modulate the function of eIF4E through the phosphorylation that is necessary for oncogenic transformation. Therefore, pharmacologic MNK inhibitors may provide a nontoxic and effective anticancer strategy. MNK1b is a truncated isoform of MNK1a that is active in the absence of stimuli. Using in vitro selection, high-affinity DNA aptamers to MNK1b were selected from a library of ssDNA. Selection was monitored using the enzyme-linked oligonucleotide assay (ELONA), and the selected aptamer population was cloned and sequenced. Four groups of aptamers were identified, and the affinities of one representative for rMNK1b were determined using ELONA and quantitative polymerase chain reaction. Two aptamers, named apMNK2F and apMNK3R, had a lower Kd in the nmol/l range. The secondary structure of the selected aptamers was predicted using mFold, and the QGRS Mapper indicated the presence of potential G-quadruplex structures in both aptamers. The selected aptamers were highly specific against MNK1, showing higher affinity to MNK1b than to MNK1a. Interestingly, both aptamers were able to produce significant translation inhibition and prevent tumor cell proliferation and migration and colony formation in breast cancer cells. These results indicate that MNK1 aptamers have an attractive therapeutic potential.

RNA Helicase Associated with AU-rich Element (RHAU/DHX36) Interacts with the 3'-Tail of the Long Non-coding RNA BC200 (BCYRN1)

RNA helicase associated with AU-rich element (RHAU) is an ATP-dependent RNA helicase that demonstrates high affinity for quadruplex structures in DNA and RNA. To elucidate the significance of these quadruplex-RHAU interactions, we have performed RNA co-immunoprecipitation screens to identify novel RNAs bound to RHAU and characterize their function. In the course of this study, we have identified the non-coding RNA BC200 (BCYRN1) as specifically enriched upon RHAU immunoprecipitation. Although BC200 does not adopt a quadruplex structure and does not bind the quadruplex-interacting motif of RHAU, it has direct affinity for RHAU in vitro. Specifically designed BC200 truncations and RNase footprinting assays demonstrate that RHAU binds to an adenosine-rich region near the 3'-end of the RNA. RHAU truncations support binding that is dependent upon a region within the C terminus and is specific to RHAU isoform 1. Tests performed to assess whether BC200 interferes with RHAU helicase activity have demonstrated the ability of BC200 to act as an acceptor of unwound quadruplexes via a cytosine-rich region near the 3'-end of the RNA. Furthermore, an interaction between BC200 and the quadruplex-containing telomerase RNA was confirmed by pull-down assays of the endogenous RNAs. This leads to the possibility that RHAU may direct BC200 to bind and exert regulatory functions at quadruplex-containing RNA or DNA 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.

Recognition of an important G-quadruplex in the HIV-1 promoter with natural small molecules

Targeting a G-quadruplex with chemical small molecules is a useful strategy for gene therapy for disease. The guanine-rich sequence d(5′-TG1G2CCTG3G4G5CG6G7G8ACTG9G10G11-3′) in the HIV-1 promoter can form a G-quadruplex structure. In this study, circular dichroism was performed to study the conformation and thermal stability of the HIV-1 G-quadruplex before and after adding small molecules. A DMS footprinting assay was used to identify which guanosine can be integrated into the G-quadruplex structure. Electrospray ionization mass spectrometry was used to evaluate the binding affinities of the small molecules with the G-quadruplex. Our results showed that G1, G2, G3, G4, G7, G8, G9, and G10 of the above oligonucleotides formed a two G-tetrad antiparallel G-quadrulex, and nitidine chloride was found to have the highest binding affinity toward the HIV-1 G-quadruplex among the eight studied small molecules. The Tm value of the G-quadruplex was enhanced from 56.6 to 63.2 °C when fourfold nitidine chloride was added. This is potentially a novel approach for anti-HIV-1 drug development.

G-quadruplex DNA structures can interfere with uracil glycosylase activity in vitro

Genome sequences that contain tandem repeats of guanine can form stable four-stranded structures known as G-quadruplex, or G4 DNA. While the molecular mechanisms are not fully defined, such guanine-rich loci are prone to mutagenesis and recombination. Various repair pathways function to reduce the potential for genome instability by correcting base damage and replication errors; however, it is not yet fully defined how well these processes function at G4 DNA. One frequent form of base damage occurs from cytidine deamination, resulting in deoxyuracil and UG mismatches. In duplex and single-stranded DNA, uracil bases are recognised and excised by uracil glycosylases. Here, we tested the efficiency of uracil glycosylase activity in vitro on uracil bases located directly adjacent to guanine repeats and G4 DNA. We show that uracil excision by bacterial UDG and human hUNG2 is reduced at uracils positioned directly 5' or 3' of a guanine tetrad. Control reactions using oligonucleotides disrupted for G4 formation or reaction conditions that do not favour G4 formation resulted in full uracil excision activity. Based on these in vitro results, we suggest that folding of guanine-rich DNA into G4 DNA results in a DNA conformation that is resistant to uracil glycosylase-initiated repair and this has the potential to increase the risk of instability at guanine repeats in the genome.

hnRNP-Q1 represses nascent axon growth in cortical neurons by inhibiting Gap-43 mRNA translation

A novel posttranscriptional mechanism for regulating the neuronal protein GAP-43 is reported. The mRNA-binding protein hnRNP-Q1 represses Gap-43 mRNA translation by a mechanism involving a 5′ untranslated region G-quadruplex structure, which affects GAP-43 function, as demonstrated by a GAP-43–dependent increase in neurite length and number with hnRNP-Q1 knockdown.

Posttranscriptional regulation of gene expression by mRNA-binding proteins is critical for neuronal development and function. hnRNP-Q1 is an mRNA-binding protein that regulates mRNA processing events, including translational repression. hnRNP-Q1 is highly expressed in brain tissue, suggesting a function in regulating genes critical for neuronal development. In this study, we have identified Growth-associated protein 43 (Gap-43) mRNA as a novel target of hnRNP-Q1 and have demonstrated that hnRNP-Q1 represses Gap-43 mRNA translation and consequently GAP-43 function. GAP-43 is a neuronal protein that regulates actin dynamics in growth cones and facilitates axonal growth. Previous studies have identified factors that regulate Gap-43 mRNA stability and localization, but it remains unclear whether Gap-43 mRNA translation is also regulated. Our results reveal that hnRNP-Q1 knockdown increased nascent axon length, total neurite length, and neurite number in mouse embryonic cortical neurons and enhanced Neuro2a cell process extension; these phenotypes were rescued by GAP-43 knockdown. Additionally, we have identified a G-quadruplex structure in the 5′ untranslated region of Gap-43 mRNA that directly interacts with hnRNP-Q1 as a means to inhibit Gap-43 mRNA translation. Therefore hnRNP-Q1–mediated repression of Gap-43 mRNA translation provides an additional mechanism for regulating GAP-43 expression and function and may be critical for neuronal development.

SMARCAL1 Negatively Regulates C-Myc Transcription By Altering The Conformation Of The Promoter Region

SMARCAL1, a member of the SWI2/SNF2 protein family, stabilizes replication forks during DNA damage. In this manuscript, we provide the first evidence that SMARCAL1 is also a transcriptional co-regulator modulating the expression of c-Myc, a transcription factor that regulates 10-15% genes in the human genome. BRG1, SMARCAL1 and RNAPII were found localized onto the c-myc promoter. When HeLa cells were serum starved, the occupancy of SMARCAL1 on the c-myc promoter increased while that of BRG1 and RNAPII decreased correlating with repression of c-myc transcription. Using Active DNA-dependent ATPase A Domain (ADAAD), the bovine homolog of SMARCAL1, we show that the protein can hydrolyze ATP using a specific region upstream of the CT element of the c-myc promoter as a DNA effector. The energy, thereby, released is harnessed to alter the conformation of the promoter DNA. We propose that SMARCAL1 negatively regulates c-myc transcription by altering the conformation of its promoter region during differentiation.

The linear plastid chromosomes of maize: terminal sequences, structures, and implications for DNA replication

The structure of a chromosomal DNA molecule may influence the way in which it is replicated and inherited. For decades plastid DNA (ptDNA) was believed to be circular, with breakage invoked to explain linear forms found upon extraction from the cell. Recent evidence indicates that ptDNA in vivo consists of linear molecules with discrete termini, although these ends were not characterized. We report the sequences of two terminal regions, End1 and End2, for maize (Zea mays L.) ptDNA. We describe structural features of these terminal regions and similarities found in other plant ptDNAs. The terminal sequences are within inverted repeat regions (leading to four genomic isomers) and adjacent to origins of replication. Conceptually, stem-loop structures may be formed following melting of the double-stranded DNA ends. Exonuclease digestion indicates that the ends in maize are unobstructed, but tobacco (Nicotiana tabacum L.) ends may have a 5'-protein. If the terminal structure of ptDNA molecules influences the retention of ptDNA, the unprotected molecular ends in mature leaves of maize may be more susceptible to degradation in vivo than the protected ends in tobacco. The terminal sequences and cumulative GC skew profiles are nearly identical for maize, wheat (Triticum aestivum L.) and rice (Oryza sativa L.), with less similarity among other plants. The linear structure is now confirmed for maize ptDNA and inferred for other plants and suggests a virus-like recombination-dependent replication mechanism for ptDNA. Plastid transformation vectors containing the terminal sequences may increase the chances of success in generating transplastomic cereals.

PAM multiplicity marks genomic target sites as inhibitory to CRISPR-Cas9 editing

In CRISPR-Cas9 genome editing, the underlying principles for selecting guide RNA (gRNA) sequences that would ensure for efficient target site modification remain poorly understood. Here we show that target sites harbouring multiple protospacer adjacent motifs (PAMs) are refractory to Cas9-mediated repair in situ. Thus we refine which substrates should be avoided in gRNA design, implicating PAM density as a novel sequence-specific feature that inhibits in vivo Cas9-driven DNA modification.

Fragile X mental retardation protein interactions with a G quadruplex structure in the 3'-untranslated region of NR2B mRNA

Fragile X syndrome, the most common cause of inherited intellectual disability, is caused by a trinucleotide CGG expansion in the 5'-untranslated region of the FMR1 gene, which leads to the loss of expression of the fragile X mental retardation protein (FMRP). FMRP, an RNA-binding protein that regulates the translation of specific mRNAs, has been shown to bind a subset of its mRNA targets by recognizing G quadruplex structures. It has been suggested that FMRP controls the local protein synthesis of several protein components of the post synaptic density (PSD) in response to specific cellular needs. We have previously shown that the interactions between FMRP and mRNAs of the PSD scaffold proteins PSD-95 and Shank1 are mediated via stable G-quadruplex structures formed within the 3'-untranslated regions of these mRNAs. In this study we used biophysical methods to show that a comparable G quadruplex structure forms in the 3'-untranslated region of the glutamate receptor subunit NR2B mRNA encoding for a subunit of N-methyl-d-aspartate (NMDA) receptors that is recognized specifically by FMRP, suggesting a common theme for FMRP recognition of its dendritic mRNA targets.

Identification of RNA Oligonucleotides Binding to Several Proteins from Potential G-Quadruplex Forming Regions in Transcribed Pre-mRNA

G-quadruplexes (G4s) are noncanonical DNA/RNA structures formed by guanine-rich sequences. Recently, G4s have been found not only in aptamers but also in the genomic DNA and transcribed RNA. In this study, we identified new RNA oligonucleotides working as aptamers by focusing on G4-forming RNAs located within the pre-mRNA. We showed that the G4 in the 5′ UTR and first intron of VEGFA bound to the protein encoded in VEGFA gene, VEGF165, with high affinity. Moreover, G4-forming RNAs located within the PDGFA and the PDGFB introns bound to PDGF-AA and PDGF-BB, respectively, indicating that G4 in the pre-mRNA could be an aptamer. It had been reported that the putative G4-forming RNA sequences are located in some parts of most genes, thus our strategy for aptamer identification could be applicable to other proteins. It has been reported that some G4-forming RNAs in 5′ UTRs are involved in translation control; however, G4-forming excised intronic RNA function has not been revealed previously. Therefore, these findings could not only contribute to the identification of RNA aptamers but also provide new insights into the biological functioning of G4-forming RNAs located within intronic RNA sequences.

Intragenic G-quadruplex structure formed in the human CD133 and its biological and translational relevance

Cancer stem cells (CSCs) have been identified in several solid malignancies and are now emerging as a plausible target for drug discovery. Beside the questionable existence of CSCs specific markers, the expression of CD133 was reported to be responsible for conferring CSC aggressiveness. Here, we identified two G-rich sequences localized within the introns 3 and 7 of the CD133 gene able to form G-quadruplex (G4) structures, bound and stabilized by small molecules. We further showed that treatment of patient-derived colon CSCs with G4-interacting agents triggers alternative splicing that dramatically impairs the expression of CD133. Interestingly, this is strongly associated with a loss of CSC properties, including self-renewing, motility, tumor initiation and metastases dissemination. Notably, the effects of G4 stabilization on some of these CSC properties are uncoupled from DNA damage response and are fully recapitulated by the selective interference of the CD133 expression.In conclusion, we provided the first proof of the existence of G4 structures within the CD133 gene that can be pharmacologically targeted to impair CSC aggressiveness. This discloses a class of potential antitumoral agents capable of targeting the CSC subpopulation within the tumoral bulk.

Hidden Genetic Variation in LCA9-Associated Congenital Blindness Explained by 5'UTR Mutations and Copy-Number Variations of NMNAT1

Leber congenital amaurosis (LCA) is a severe autosomal‐recessive retinal dystrophy leading to congenital blindness. A recently identified LCA gene is NMNAT1, located in the LCA9 locus. Although most mutations in blindness genes are coding variations, there is accumulating evidence for hidden noncoding defects or structural variations (SVs). The starting point of this study was an LCA9‐associated consanguineous family in which no coding mutations were found in the LCA9 region. Exploring the untranslated regions of NMNAT1 revealed a novel homozygous 5′UTR variant, c.‐70A>T. Moreover, an adjacent 5′UTR variant, c.‐69C>T, was identified in a second consanguineous family displaying a similar phenotype. Both 5′UTR variants resulted in decreased NMNAT1 mRNA abundance in patients’ lymphocytes, and caused decreased luciferase activity in human retinal pigment epithelial RPE‐1 cells. Second, we unraveled pseudohomozygosity of a coding NMNAT1 mutation in two unrelated LCA patients by the identification of two distinct heterozygous partial NMNAT1 deletions. Molecular characterization of the breakpoint junctions revealed a complex Alu‐rich genomic architecture. Our study uncovered hidden genetic variation in NMNAT1‐associated LCA and emphasized a shift from coding to noncoding regulatory mutations and repeat‐mediated SVs in the molecular pathogenesis of heterogeneous recessive disorders such as hereditary blindness.

DNA aptamers against the Cry j 2 allergen of Japanese cedar pollen for biosensing applications

Sensing pollen allergens is required to prevent allergic disorders such as pollinosis. Aptamers, which bind to specific molecules, offer great potential as useful tools for detecting pollen allergens as measures against allergic disorders. Here, we report the identification of DNA aptamers binding to Cry j 2, one of the major allergens in Japanese cedar pollen, and the histochemical sensing of Cry j 2 in ruptured Japanese cedar pollen. DNA aptamers were selected by systematic evolution of ligands by exponential enrichment (SELEX) using nitrocellulose membranes. Through four rounds of SELEX, we identified aptamers binding to Cry j 2. The aptamers generated staining in ruptured Japanese cedar pollen on glass slides without extraction, similar to anti-Cry j 2 antibodies. The staining was compatible with starch localization, in which Cry j 2 is present. An aptamer, CJ2-06, which had high and specific binding ability to Cry j 2 (K(d)=24 nM), detected an amount of Cry j 2 equivalent to that in several tens of micrograms of pollen. Cry j 2 contained in house dust was detected in a spike test. The aptamers identified in this study can be powerful tools for allergen recognition in the practical biosensing of Cry j 2, leading to preventive measures against allergic disorders caused by Japanese cedar pollen.

Binding of Multiple Rap1 Proteins Stimulates Chromosome Breakage Induction during DNA Replication

Telomeres, the ends of linear eukaryotic chromosomes, have a specialized chromatin structure that provides a stable chromosomal terminus. In budding yeast Rap1 protein binds to telomeric TG repeat and negatively regulates telomere length. Here we show that binding of multiple Rap1 proteins stimulates DNA double-stranded break (DSB) induction at both telomeric and non-telomeric regions. Consistent with the role of DSB induction, Rap1 stimulates nearby recombination events in a dosage-dependent manner. Rap1 recruits Rif1 and Rif2 to telomeres, but neither Rif1 nor Rif2 is required for DSB induction. Rap1-mediated DSB induction involves replication fork progression but inactivation of checkpoint kinase Mec1 does not affect DSB induction. Rap1 tethering shortens artificially elongated telomeres in parallel with telomerase inhibition, and this telomere shortening does not require homologous recombination. These results suggest that Rap1 contributes to telomere homeostasis by promoting chromosome breakage.

Telomere length is maintained primarily through equilibrium between telomerase-mediated lengthening and the loss of telomeric sequence through the end-replication problem. In budding yeast Rap1 protein binds to telomeric TG repeat and negatively regulates telomerase recruitment in a dosage-dependent manner. In this paper we provide evidence suggesting an alternative Rap1-dependent telomere shortening mechanism in which binding of multiple Rap1 proteins mediates DNA break induction during DNA replication. This process does not involve recombination events; therefore, it is distinct from loop-mediated telomere trimming.

High-throughput sequencing of DNA G-quadruplex structures in the human genome

G-quadruplexes (G4s) are nucleic acid secondary structures that form within guanine-rich DNA or RNA sequences. G4 formation can affect chromatin architecture and gene regulation and has been associated with genomic instability, genetic diseases and cancer progression. Here we present a high-resolution sequencing-based method to detect G4s in the human genome. We identified 716,310 distinct G4 structures, 451,646 of which were not predicted by computational methods. These included previously uncharacterized noncanonical long loop and bulged structures. We observed a high G4 density in functional regions, such as 5' untranslated regions and splicing sites, as well as in genes previously not predicted to contain these structures (such as BRCA2). G4 formation was significantly associated with oncogenes, tumor suppressors and somatic copy number alterations related to cancer development. The G4s identified in this study may therefore represent promising targets for cancer intervention.

Phosphoproteomic Analysis of KSHV-Infected Cells Reveals Roles of ORF45-Activated RSK during Lytic Replication

Kaposi’s Sarcoma-Associated Herpesvirus (KSHV) is an oncogenic virus which has adapted unique mechanisms to modulate the cellular microenvironment of its human host. The pathogenesis of KSHV is intimately linked to its manipulation of cellular signaling pathways, including the extracellular signal-regulated kinase (ERK) mitogen-activated protein kinase (MAPK) pathway. We have previously shown that KSHV ORF45 contributes to the sustained activation of both ERK and p90 ribosomal S6 kinase (RSK, a major functional mediator of ERK/MAPK signaling) during KSHV lytic replication. ORF45-activated RSK is required for optimal KSHV lytic gene expression and progeny virion production, though the underlying mechanisms downstream of this activation are still unclear. We hypothesized that the activation of RSK by ORF45 causes differential phosphorylation of cellular and viral substrates, affecting biological processes essential for efficient KSHV lytic replication. Accordingly, we observed widespread and significant differences in protein phosphorylation upon induction of lytic replication. Mass-spectrometry-based phosphoproteomic screening identified putative substrates of ORF45-activated RSK in KSHV-infected cells. Bioinformatic analyses revealed that nuclear proteins, including several transcriptional regulators, were overrepresented among these candidates. We validated the ORF45/RSK-dependent phosphorylation of several putative substrates by employing KSHV BAC mutagenesis, kinase inhibitor treatments, and/or CRISPR-mediated knockout of RSK in KSHV-infected cells. Furthermore, we assessed the consequences of knocking out these substrates on ORF45/RSK-dependent regulation of gene expression and KSHV progeny virion production. Finally, we show data to support that ORF45 regulates the translational efficiency of a subset of viral/cellular genes with complex secondary structure in their 5’ UTR. Altogether, these data shed light on the mechanisms by which KSHV ORF45 manipulates components of the host cell machinery via modulation of RSK activity. Thus, this study has important implications for the pathobiology of KSHV and other diseases in which RSK activity is dysregulated.

Kaposi’s sarcoma-associated herpesvirus (KSHV) is a human tumor virus which hijacks the host signaling pathways in order to maintain persistent infection. We previously discovered that the KSHV protein ORF45 binds to and activates the cellular kinase RSK (p90 ribosomal S6 kinase), and that this activation is vital for optimal KSHV gene expression and virion production. Here, we performed a phosphoproteomic analysis of KSHV-infected cells to further characterize the specific substrates of ORF45-activated RSK. Bioinformatic analyses provided insights into the functional roles of these substrates. We verified the ORF45/RSK-dependent phosphorylation of a subset of these substrates by various means. Finally, we used genome editing to knock out RSK, as well as several cellular substrates identified by our screening, and characterized the consequent effect(s) on regulation of gene expression and virion production. Thus, this work further elucidates one of the key signaling nodes modulated by KSHV, and implicates ORF45-mediated activation of RSK in the regulation of viral and host gene expression during KSHV lytic replication.

Selection and Characterization of DNA Aptamers Targeting All Four Serotypes of Dengue Viruses

Dengue viruses (DENVs) are members of Flaviviridae family, which are associated with human disease. The envelope (E) protein plays an important role in viral infection. However, there is no effective antibody for clinical treatment due to antibody dependent enhancement of infection. In this study, using Systematic Evolution of Ligands by Exponential Enrichment (SELEX), we demonstrated the first aptamer (S15) that can bind to DENV-2 envelop protein domain III (ED3) with a high binding affinity. S15 was found to form a parallel quadruplex based on Quadfinder prediction, gel mobility assay and circular dichroism studies. Both the quadruplex structure and the sequence on 5’-end were necessary for the binding activity of S15. NMR titration experiments indicated that S15 bound to a highly conserved loop between β_A and β_B strands of ED3. Moreover, S15 can neutralize the infections by all four serotypes of DENVs. Our result provides a new opportunity in the development of DNA aptamers against DENVs in the future.

G4RNA: an RNA G-quadruplex database

G-quadruplexes (G4) are tetrahelical structures formed from planar arrangement of guanines in nucleic acids. A simple, regular motif was originally proposed to describe G4-forming sequences. More recently, however, formation of G4 was discovered to depend, at least in part, on the contextual backdrop of neighboring sequences. Prediction of G4 folding is thus becoming more challenging as G4 outlier structures, not described by the originally proposed motif, are increasingly reported. Recent observations thus call for a comprehensive tool, capable of consolidating the expanding information on tested G4s, in order to conduct systematic comparative analyses of G4-promoting sequences. The G4RNA Database we propose was designed to help meet the need for easily-retrievable data on known RNA G4s. A user-friendly, flexible query system allows for data retrieval on experimentally tested sequences, from many separate genes, to assess G4-folding potential. Query output sorts data according to sequence position, G4 likelihood, experimental outcomes and associated bibliographical references. G4RNA also provides an ideal foundation to collect and store additional sequence and experimental data, considering the growing interest G4s currently generate.

Database URL:scottgroup.med.usherbrooke.ca/G4RNA